Differences between revisions 1 and 2
⇤ ← Revision 1 as of 2010-08-11 13:32:32 →
Size: 607
Editor: hardikgajera
Comment:
← Revision 2 as of 2010-08-11 14:16:56 → ⇥
Size: 614
Editor: hardikgajera
Comment:
Deletions are marked like this. Additions are marked like this.
Line 17: Line 17:
print 'The motion of the given mass-system is' print 'The motion of the given mass-spring system is'

  • Book

    Advanced Engineering Mathematics

    Author

    Erwin Kreyszig

    Edition

    8th Edition 

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])

  • Solution by:

    • Hardik Gajera, student, IISER Pune
    • Lokesh Pimpale, student, IISER Pune

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