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

Kreyszig-1.3-20-U (last edited 2010-12-18 11:18:32 by sarvani.kotha)