t=var('t')
y1=function('y1',t)
y2=function('y2',t)
deq1=diff(y1,t)-4*y1+8*y2-2*cosh(t)
deq2=diff(y2,t)-2*y1+6*y2-cosh(t)-2*sinh(t)
sol=desolve_system([deq1,deq2],[y1,y2],ics=[0,0,0])
sol

SageDays/Kreyszig-3.6-10 (last edited 2010-08-12 13:32:58 by PraveenKumar)