{VERSION 4 0 "IBM INTEL NT" "4.0" } {USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 1 0 0 0 0 1 }{CSTYLE "" -1 256 "" 1 18 0 0 0 0 0 1 0 0 0 0 0 0 0 0 } {CSTYLE "" -1 257 "" 1 14 0 0 0 0 0 0 0 0 0 0 0 0 0 0 }{PSTYLE "Normal " -1 0 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }1 1 0 0 0 0 1 0 1 0 2 2 0 1 }} {SECT 0 {EXCHG {PARA 0 "" 0 "" {TEXT 256 35 "General Motion Coupled Os cillators " }{TEXT 257 77 "of two masses linearly connected with three springs (Figure 11.3.1 in Fowles)" }}{PARA 0 "" 0 "" {TEXT -1 64 " \+ We can choose whatever starting positions/velocities we wish." }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "restart:with(plots):" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 88 "m:=1: k:=5: kp :=0.2: x1i := 0: x2i := 1: v1i := 0: v2i := 0: # CHANGE THESE VALUES ONLY" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 40 "w1 := sqrt(k/m): w2 := sqrt( (k+2*kp)/m):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 89 "a1:=0.5*(x1 i+x2i): a3:=0.5*(x1i-x2i): a2:=0.5*(1/w1)*(v1i+v2i): a4:=0.5*(1/w2)*(v 1i-v2i):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 58 "x1 := a1*cos(w1 *t)+a2*sin(w1*t)+a3*cos(w2*t)+a4*sin(w2*t):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 58 "x2 := a1*cos(w1*t)+a2*sin(w1*t)-a3*cos(w2*t)-a4*si n(w2*t):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 171 "plot([-3+x1,3+ x2],t=0..50,view=-5..5,thickness=[2,2],numpoints=1000,discont=true,col or=[red,blue],labels=[\"t\",\"X\"],title =\"X1=Red, X2=Blue\",resoluti on=50); #,view=-1.0..+1" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 5 " " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}{MARK "10" 0 } {VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }