t=var('t)
y=function('y',t)
de=8*diff(y,t,2)-2*diff(y,t)-y
desolve(de,[y,t],ics=[0,-0.2,-0.325])

SageDays/Kreyszig-2.2-13 (last edited 2010-08-10 13:30:34 by aitians)