t=var('t')
Y=function('Y',t)
DE = 2*diff(diff(Y,t),t) + 200*Y-220*4*cos(4*t) 
d(t)=desolve (DE,[Y,t])
d1(t)=diff(d(t),t)
k1,k2=var('k1,k2')
solve([d(0)==0,
       d1(0)==0],k1,k2)
d=d(t).subs(k1 = 0, k2 = (-110/21))
d

SageDays/Kreyszig-2.12-15 (last edited 2010-08-12 09:20:57 by hardikgajera)