{VERSION 5 0 "IBM INTEL NT" "5.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 "2D Math" -1 2 "Times" 0 1 0 0 0 0 0 0 2 0 0 0 0 0 0 1 }{CSTYLE "2D Output" 2 20 "" 0 1 0 0 255 1 0 0 0 0 0 0 0 0 0 1 } {PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Text Output" -1 2 1 {CSTYLE "" -1 -1 "Courier" 1 10 0 0 255 1 0 0 0 0 0 1 3 0 3 0 }1 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Heading 1" 0 3 1 {CSTYLE "" -1 -1 "" 1 18 0 0 0 0 0 1 0 0 0 0 0 0 0 0 }1 0 0 0 8 4 0 0 0 0 0 0 -1 0 } {PSTYLE "Warning" 2 7 1 {CSTYLE "" -1 -1 "" 0 1 0 0 255 1 0 0 0 0 0 0 1 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Maple Output" 0 11 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 }3 3 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Maple Plot" 0 13 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 }3 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Title" 0 18 1 {CSTYLE "" -1 -1 "" 1 18 0 0 0 0 0 1 1 0 0 0 0 0 0 0 } 3 0 0 -1 12 12 0 0 0 0 0 0 19 0 }{PSTYLE "" 18 256 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 2 0 0 0 0 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 } {PSTYLE "" 18 257 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 2 0 0 0 0 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }} {SECT 0 {PARA 256 "" 0 "" {TEXT -1 12 "Mass-Spring " }}{PARA 257 "" 0 "" {TEXT -1 6 "System" }}{EXCHG }{SECT 1 {PARA 3 "" 0 "" {TEXT -1 26 " Initialize Maple Worksheet" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 36 "restart:\nwith(plots):\nwith(DEtools):" }}{PARA 7 "" 1 "" {TEXT -1 50 "Warning, the name changecoords has been redefined\n" }}}}{SECT 1 {PARA 3 "" 0 "" {TEXT -1 26 "Set up physics (constants)" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 52 "m := 1;\nk := 1;\nb := 1;\nA := 1; \nB := 1;\nF := t -> 0;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"mG\"\" \"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"kG\"\"\"" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#>%\"bG\"\"\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"A G\"\"\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"BG\"\"\"" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#>%\"FG\"\"!" }}}}{SECT 1 {PARA 3 "" 0 "" {TEXT -1 36 "Set up o.d.e. and initial conditions" }}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 53 "ode := m*diff(x(t),t$2)+b*diff(x(t),t)+k*x(t) = F(t );" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$odeG/,(-%%diffG6$-%\"xG6#%\"t G-%\"$G6$F-\"\"#\"\"\"-F(6$F*F-F2F*F2\"\"!" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 30 "ics := x(0) = A, D(x)(0) = B ;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$icsG6$/-%\"xG6#\"\"!\"\"\"/--%\"DG6#F(F)F+" }}}} {SECT 1 {PARA 3 "" 0 "" {TEXT -1 27 "Solve o.d.e. (using dsolve)" }} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "dsolve( \{ode,ics\}, x(t) ); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"xG6#%\"tG,&*(\"\"$#\"\"\"\"\" #-%$expG6#,$*&F-!\"\"F'F,F3F,-%$sinG6#,$*(F-F3F*F+F'F,F,F,F,*&F.F,-%$c osGF6F,F," }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "sol := unapply (rhs(%),t);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$solGf*6#%\"tG6\"6$%) operatorG%&arrowGF(,&*(\"\"$#\"\"\"\"\"#-%$expG6#,$*&#F0F1F09$F0!\"\"F 0-%$sinG6#,$*&F/F0*&F.F/F8F0F0F0F0F0*&F2F0-%$cosGF " 0 "" {MPLTEXT 1 0 21 "plot(sol(t),t=0..10);" }}{PARA 13 " " 1 "" {GLPLOT2D 594 379 379 {PLOTDATA 2 "6%-%'CURVESG6$7^o7$$\"\"!F)$ \"\"\"F)7$$\"3GLLL3x&)*3\"!#=$\"3=G2&3DHt4\"!#<7$$\"3emmm;arz@F/$\"3)3 *Q!p)>Fs6F27$$\"3v***\\7y%*z7$F/$\"3.Q([**>2/A\"F27$$\"3[LL$e9ui2%F/$ \"3+[I.#RDPD\"F27$$\"3\"pm\"H2Q\\4YF/$\"35xXV'pYjE\"F27$$\"3z***\\(oMr U^F/$\"313Xves![F\"F27$$\"3Cm\"z%*HB$4aF/$\"3%>T7DTPvF\"F27$$\"3oK$3-8 Lfn&F/$\"3-#H2qq*Hz7F27$$\"3B+v$4'HaUfF/$\"3cGaMp$=,G\"F27$$\"3nmmm\"z _\"4iF/$\"3W_()3^w,!G\"F27$$\"3ImmT&))HvZ'F/$\"3y^iDxA,z7F27$$\"3/nm;z p!fu'F/$\"3[#QG#3Y7x7F27$$\"3ynm\"H2%G9qF/$\"3A\\JK\\$zVF\"F27$$\"3Unm mm6m#G(F/$\"3SV33h6!3F\"F27$$\"3ynm;a`T>yF/$\"3Q0&*HwOdB'>s6F27$$\"3smm\"z/3uC\"F2$\"31hl*eNC82\"F27$$\"3!** **\\7LRDX\"F2$\"3ca)zK_v,Y*F/7$$\"3[LLekGhe:F2$\"3KGg%)3cjd()F/7$$\"3% om;zR'ok;F2$\"3d0\"R/s'RK!)F/7$$\"3OLL3_(>/x\"F2$\"35gM+!4boH(F/7$$\"3 3++D1J:w=F2$\"3C:TTLbxdlF/7$$\"3+n;HdG\"\\)>F2$\"3j@UK!QoC!eF/7$$\"3oL LL3En$4#F2$\"3?Yw-Hz-g]F/7$$\"3_++Dc#o%*=#F2$\"3d`c6%fBBU%F/7$$\"3#pmm T!RE&G#F2$\"3!*=,pWP?/QF/7$$\"3D+++D.&4]#F2$\"3#H%HOvF2-DF/7$$\"3;+++v B_7$$\"3&om;z *ev:JF2$!3Ew@YsadzLFcu7$$\"3_LLL347TLF2$!3a.\\`FAgB5F/7$$\"3nLLLLY.KNF 2$!3eBuWNNTf9F/7$$\"33++D\"o7Tv$F2$!37+ic&[n5\"=F/7$$\"3?LLL$Q*o]RF2$! 3Ko6*HD1n*>F/7$$\"3m++D\"=lj;%F2$!3?&*R!eA7D3#F/7$$\"3S++vV&RF/7$$\"3]LLeR\"3Gy%F2$!3=kr YdB-G=F/7$$\"3emm;/T1&*\\F2$!3WsLvg@KI;F/7$$\"3=nm\"zRQb@&F2$!3GM()zli R&R\"F/7$$\"3:++v=>Y2aF2$!3!4O\\Y)H()y6F/7$$\"3Znm;zXu9cF2$!3K!G+C)\\' [V*Fcu7$$\"34+++]y))GeF2$!3!oE/&\\Om$3(Fcu7$$\"3H++]i_QQgF2$!3k2Z[2xsR \\Fcu7$$\"3b++D\"y%3TiF2$!3CIsO_z9lIFcu7$$\"3+++]P![hY'F2$!33)[WcoPXD \"Fcu7$$\"3iKLL$Qx$omF2$\"3G3\"e&>6h(4\"!#?7$$\"3Y+++v.I%)oF2$\"3]4xhW %=jG\"Fcu7$$\"3?mm\"zpe*zqF2$\"3BA[28XJ3@Fcu7$$\"3;,++D\\'QH(F2$\"39o- 'f\"H#ov#Fcu7$$\"3%HL$e9S8&\\(F2$\"3WD()p\\`<\\JFcu7$$\"3s++D1#=bq(F2$ \"3E&\\>f!H:gLFcu7$$\"3\"HLL$3s?6zF2$\"3cRO+-#4$)R$Fcu7$$\"3a***\\7`Wl 7)F2$\"3YDS**>sR#H$Fcu7$$\"3enmmm*RRL)F2$\"3k,)zJ#y1!3$Fcu7$$\"3%zmmTv Jga)F2$\"3yGd))y9(4y#Fcu7$$\"3]MLe9tOc()F2$\"3T2&=3,E)HCFcu7$$\"31,++] Qk\\*)F2$\"3K*ospA3-3#Fcu7$$\"3![LL3dg6<*F2$\"3IuChf!)Hp;Fcu7$$\"3%ymm mw(Gp$*F2$\"3S_h'eTT%38Fcu7$$\"3C++D\"oK0e*F2$\"3d`^=;%)3U%*F^[l7$$\"3 5,+v=5s#y*F2$\"3O3**)*fn-]iF^[l7$$\"#5F)$\"3`qLtwQO:KF^[l-%'COLOURG6&% $RGBG$F\\`l!\"\"F(F(-%+AXESLABELSG6$Q\"t6\"Q!Fi`l-%%VIEWG6$;F(F[`l%(DE FAULTG" 1 2 0 1 10 0 2 9 1 4 2 1.000000 42.000000 47.000000 0 0 "Curve 1" }}}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}{MARK "7" 0 } {VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }