x = var('x')   
y = function('y',x)   
sol=desolve (diff(y,x)+3*y -sin(x) , [y,x], ics=(pi/2,0.3))
sol

SageDays/Kreyszig-1.6-19 (last edited 2010-08-11 10:50:45 by SWATI AHUJA)