t,w=var('t,w')
Y=function('Y',t)
DE = diff(diff(Y,t),t) + Y-cos(w*t) 
de=desolve (DE,[Y,t])
print 'Transient solution is'
show(de)
print ',where w is not equal to 1'

SageDays/Kreyszig-2.11-9 (last edited 2010-08-11 11:56:12 by hardikgajera)