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