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

Kreyszig-1.6-16-U (last edited 2010-12-17 10:41:24 by mlec)