t=var('t')
y1=function('y1',t)
y2=function('y2',t)
deq1=diff(y1,t)-3*y1-y2+3*sin(3*t)
deq2=diff(y2,t)-7*y1+3*y2-9*cos(3*t)+16*sin(3*t)
sol=desolve_system([deq1,deq2],[y1,y2])
sol

SageDays/Kreyszig-3.6-5 (last edited 2010-08-12 12:49:52 by PraveenKumar)