t=var('t')
y=function('y',t)
deq1=diff(y,t,2)-4*y-2==0
desolve(deq1,[y,t])
y2=function('y2',t)
deq2=diff(y2,t,1)-2*(desolve(deq1,[y,t]))-t==0
desolve(deq2,[y2,t])

Kreyszig-3.6-2-U (last edited 2010-12-17 12:26:40 by sarvani.kotha)