x=var('x')
y=function('y',x)
deq1=(2*x*y)*diff(y,x)+(x-1)*y^2-x^2*e^x==0
r=desolve(deq1,[y,x],y)
show(r.simplify_full())

Kreyszig-1.6-38-U (last edited 2010-12-18 06:26:06 by 10)