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

SageDays/Kreyszig-2.9-3-S (last edited 2010-08-12 10:36:15 by aitians)