x=var('x')
y=function('y',x)
deq=diff(y,x,1)-(y+4*x)^2==0
desolve(deq,[y,x],contrib_ode=True)