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

SageDays/Kreyszig-2.2-11 (last edited 2010-08-10 13:36:01 by prabhat)