var('x')
y=function('y',x)
print desolve(diff(y,x)-y*(1+3*(x^(-1)))-x-2,[y,x],ics=(1,exp(1)-1))

SageDays/Kreyszig-1.6-16 (last edited 2010-08-09 13:59:48 by SWATI AHUJA)