t=var('t')
Y=function('Y',t)
DE = diff(diff(Y,t),t) + 2*diff(Y,t) + 5*Y + 13*sin(3*t) 
de=desolve (DE,[Y,t])
DE0=diff(diff(Y,t),t) + 2*diff(Y,t) + 5*Y
de0=desolve(DE0,[Y,t])
print 'Steady-state solution is'
show(de-de0)

SageDays/Kreyszig-2.11-2 (last edited 2010-08-11 14:03:33 by hardikgajera)