t=var('t')
Y=function('Y',t)
DE = diff(diff(Y,t),t) + 25*Y - 24*sin(t) 
de(t)=desolve (DE,[Y,t])
print de(t)
de1(t)=diff(de(t),t)
k1,k2=var('k1,k2')
solve([de(0)==1,    
       de1(0)==1],k1,k2)
de(t)=de(t).subs(k1 =0, k2= 1)
print 'The motion of the given mass-spring system is' 
print(de(t))
plot(de(t),[0,20])

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