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