t=var('t')
y2=function('y2',t)
deq=diff(y2,t,2)-y2+8*exp(3*t)==0
desolve(deq,[y2,t])
y1=function('y1',t)
deq1=diff(y1,t,1)-desolve(deq,[y2,t])-exp(3*t)==0
desolve(deq1,[y1,t])

Kreyszig-3.6-3-U (last edited 2010-12-18 06:43:54 by kenguva.monica)