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

Kreyszig-1.6-12 (last edited 2010-12-17 09:10:01 by punchagan)