x=var('x')
y=function('y',x)
de=diff(y,x,2)+2*diff(y,x)+5*y-1.25*exp(0.5*x)-40*cos(4*x)-55*sin(4*x)
desolve(de,[y,x],ics=[0,0.2,60.1])

SageDays/Kreyszig-2.9-4-S (last edited 2010-08-12 10:32:58 by aitians)