t=var('t')
Y=function('Y',t)
DE = 4*diff(diff(Y,t),t) +8*diff(Y,t)+ 3*Y-425*sin(2*t) 
de(t)=desolve (DE,[Y,t])
de1(t)=diff(de,t)
k1,k2=var('k1,k2')
solve([de(0)==-16,    
       de1(0)==-26],k1,k2)
de(t)=de(t).subs(k1 =0, k2= 0)
print 'The motion of the given mass-spring system is' 
print(de(t))
plot(de(t),[0,20])

SageDays/Kreyszig-2.11-15 (last edited 2010-08-11 14:16:56 by hardikgajera)