#Kreyszig Pg.108 Problem set 2.9 Question 20
var('x')
y=function('y',x)
print desolve(diff(y,x,2)+2*diff(y,x)+5*y==1.25*exp(0.5*x)+40*cos(4*x)-55*sin(4*x),[y,x],ics=[0,0,1])

SageDays/Kreyszig-2.9-20 (last edited 2010-08-11 11:52:57 by S Micheal Ammal)