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

Kreyszig-1.6-20-U (last edited 2010-12-18 12:04:19 by 10)