{VERSION 6 0 "IBM INTEL NT" "6.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 Comment" 2 18 "" 0 1 0 0 0 0 0 0 0 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 } {CSTYLE "" -1 256 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 257 "" 0 14 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 258 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 259 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 260 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 } {CSTYLE "" -1 261 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 262 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 263 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 264 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 265 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 } {CSTYLE "" -1 266 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 267 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 268 "" 0 14 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 269 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 270 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 } {CSTYLE "" -1 271 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 272 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 273 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 274 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 275 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 } {CSTYLE "" -1 276 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 277 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 278 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 279 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 280 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 } {CSTYLE "" -1 281 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 282 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 283 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 284 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 285 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 } {CSTYLE "" -1 286 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 287 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 288 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 289 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 290 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 } {CSTYLE "" -1 291 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 292 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 293 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 294 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 295 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 } {CSTYLE "" -1 296 "" 1 18 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 297 "" 1 18 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 298 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 299 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 300 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 } {CSTYLE "" -1 301 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 302 "" 0 1 0 0 0 0 0 1 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 }{PSTYLE "Heading 1" -1 3 1 {CSTYLE "" -1 -1 "Times " 1 18 0 0 0 1 2 1 2 2 2 2 1 1 1 1 }1 1 0 0 8 4 1 0 1 0 2 2 0 1 } {PSTYLE "Maple Output" -1 11 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }3 3 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Maple Plot " -1 13 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 } 3 1 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Normal" -1 256 1 {CSTYLE "" -1 -1 "Times" 1 18 0 0 0 1 2 1 2 2 2 2 1 1 1 1 }3 1 0 0 0 0 1 0 1 0 2 2 0 1 }} {SECT 0 {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{SECT 1 {PARA 3 "" 0 "" {TEXT -1 33 "1. Solving Differential Equations" }{TEXT 260 20 ": The dsolve command" }}{PARA 0 "" 0 "" {TEXT -1 71 "Differential equ ations and initial value problems are solved using the " }{TEXT 271 6 "dsolve" }{TEXT -1 10 " command. " }}{PARA 0 "" 0 "" {TEXT -1 7 "usage : " }}{PARA 0 "" 0 "" {TEXT 257 7 "dsolve(" }{TEXT -1 0 "" }{TEXT 258 1 " " }{TEXT -1 25 "the differential equation" }{TEXT 266 1 " " } {TEXT -1 24 ", the dependent variable" }{TEXT 267 4 ") " }{TEXT -1 2 "or" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }{TEXT 261 10 "dsolve( \{ " } {TEXT -1 0 "" }{TEXT 262 45 "a differential equation, initial conditio n(s)" }{TEXT -1 0 "" }{TEXT 263 1 " " }{TEXT -1 0 "" }{TEXT 264 4 "\}, " }{TEXT -1 22 "the dependent variable" }{TEXT 259 2 " )" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 31 "I start most sect ions with the " }{TEXT 275 7 "restart" }{TEXT -1 55 " command to clear all previous variable definitions by " }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "restart;" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 " " 0 "" {TEXT 265 9 "Example 1" }{TEXT -1 71 ": Find the general soluti on to the first order differential equation: " }{TEXT 280 8 "y' = 3 y " }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 47 "sol1 := dsolve(diff(y(x) ,x) = 3 * y(x), y(x));\n" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%sol1G/- %\"yG6#%\"xG*&%$_C1G\"\"\"-%$expG6#,$*&\"\"$F,F)F,F,F," }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 12 "Notice: The " }{XPPEDIT 18 0 "_C1;" "6#%$ _C1G" }{TEXT -1 34 " indicates an arbitrary constant. " }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT 256 9 "Example 2" }{TEXT -1 39 ": Consider the initial value problem " }{TEXT 276 8 "y' = 3 y " }{TEXT -1 5 " , " }{TEXT 277 12 "y(0) = 1/2. " }}{PARA 0 "" 0 "" {TEXT -1 1 "(" }{TEXT 278 1 "a" }{TEXT -1 21 ") Find the solution. " } }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 59 "sol2 := dsolve(\{diff(y(x), x) = 3*y(x), y(0) = 1/2\}, y(x));\n" }}{PARA 11 "" 1 "" {XPPMATH 20 "6 #>%%sol2G/-%\"yG6#%\"xG,$*&#\"\"\"\"\"#F--%$expG6#,$*&\"\"$F-F)F-F-F-F -" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 94 "Notice: Maple automatically \+ solves for the arbitrary constant based on the initial condition. " }} }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 28 "y2 := unapply(rhs(sol2),x);\n" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#y2Gf*6#%\"xG6\"6$%)operatorG%&arrowGF(,$*&#\"\"\"\" \"#F/-%$expG6#,$*&\"\"$F/9$F/F/F/F/F(F(F(" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 49 "You can plot the solution from x = 0 to x = 1 by " }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 85 "plot(y2(x),x=0..1,labels=[\" x\",\"y2(x)\"],title=\"The solution to y' = 3y, y(0) = 1/2\");\n" }} {PARA 13 "" 1 "" {GLPLOT2D 485 485 485 {PLOTDATA 2 "6&-%'CURVESG6$7S7$ $\"\"!F)$\"3++++++++]!#=7$$\"3emmm;arz@!#>$\"3Rh&3XH%)yL&F,7$$\"3[LL$e 9ui2%F0$\"3WgX!p&))R]cF,7$$\"3nmmm\"z_\"4iF0$\"39\"QNG5lP-'F,7$$\"3[mm mT&phN)F0$\"3(eL#>Jp_CkF,7$$\"3CLLe*=)H\\5F,$\"3O74C8U&)\\oF,7$$\"3gmm \"z/3uC\"F,$\"3Wz8%)pBIpsF,7$$\"3%)***\\7LRDX\"F,$\"3noDH'H-2t(F,7$$\" 3]mm\"zR'ok;F,$\"39OK6T/rQ#)F,7$$\"3w***\\i5`h(=F,$\"3Ywlmc%4$y()F,7$$ \"3WLLL3En$4#F,$\"3AT:)HB\\-P*F,7$$\"3qmm;/RE&G#F,$\"3$)*Hj'fbgC**F,7$ $\"3\")*****\\K]4]#F,$\"3g$HMq#=!)e5!#<7$$\"3$******\\PAvr#F,$\"3w(3;F &y()H6Fdo7$$\"3)******\\nHi#HF,$\"3(z)zTdO)G?\"Fdo7$$\"3jmm\"z*ev:JF,$ \"3/!=V!)yeKF\"Fdo7$$\"3?LLL347TLF,$\"3mf=()\\*>BO\"Fdo7$$\"3,LLLLY.KN F,$\"3'*[SlCLiU9Fdo7$$\"3w***\\7o7Tv$F,$\"3Ybio*y4?a\"Fdo7$$\"3'GLLLQ* o]RF,$\"3=&pc#f;oN;Fdo7$$\"3A++D\"=lj;%F,$\"3v*oVFk8]u\"Fdo7$$\"31++vV &R#*[l8z>Fdo7$$\"3GLLe R\"3Gy%F,$\"3Sb&fjo\"\\*4#Fdo7$$\"3cmm;/T1&*\\F,$\"3k*p=\"=)GvB#Fdo7$$ \"3&em;zRQb@&F,$\"3It\"=N=G0R#Fdo7$$\"3\\***\\(=>Y2aF,$\"3rEI&fN2A`#Fd o7$$\"39mm;zXu9cF,$\"3;bM&**>rYp#Fdo7$$\"3l******\\y))GeF,$\"3A@)>rfmM (GFdo7$$\"3'*)***\\i_QQgF,$\"3Oh:#*ew&)fIFdo7$$\"3@***\\7y%3TiF,$\"33j &)yJ5q^KFdo7$$\"35****\\P![hY'F,$\"3+D![d " 0 "" {MPLTEXT 1 0 17 "subs(x=1,y2(x));\n " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,$*&#\"\"\"\"\"#F&-%$expG6#\"\"$F& F&" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 53 "and get a numerical value t o 6 significant digits by " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "evalf(%,6);\n" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"'G/5!\"%" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 23 "solve(y2(x) = 100, x);\n" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,$*&#\"\"\"\"\"$F&-%#lnG6#\"$+#F&F&" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "evalf(%,6);\n" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"'5m " 0 "" {MPLTEXT 1 0 44 "ode3 := diff(y(x),x$2)+diff( y(x),x)+y(x)=0;\n" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%ode3G/,(-%%dif fG6$-%\"yG6#%\"xG-%\"$G6$F-\"\"#\"\"\"-F(6$F*F-F2F*F2\"\"!" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 28 "and then solve it using the " }{TEXT 272 6 "dsolve" }{TEXT -1 8 " command" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "dsolve(ode3, y(x));\n" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"yG6#%\"xG,&*(%$_C1G\"\"\"-%$expG6#,$*&\"\"#!\"\"F'F+F2F+-%$ sinG6#,$*(F1F2\"\"$#F+F1F'F+F+F+F+*(%$_C2GF+F,F+-%$cosGF5F+F+" }}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 9 "Again, " }{XPPEDIT 18 0 "_C1;" "6 #%$_C1G" }{TEXT -1 5 " and " }{XPPEDIT 18 0 "_C2;" "6#%$_C2G" }{TEXT -1 31 " indicate arbitrary constants. " }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT 270 9 "Example 4" }{TEXT -1 59 ": Find and \+ plot the solution to the initial value problem: " }{TEXT 282 37 "y'' + y' + y = 0, y(0) = 1, y'(0) = 0" }{TEXT -1 1 "." }}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 39 "dsolve(\{ode3,y(0)=1,D(y)(0)=0\}, y(x));\n" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"yG6#%\"xG,&*&#\"\"\"\"\"$F+*(F,#F +\"\"#-%$expG6#,$*&F/!\"\"F'F+F5F+-%$sinG6#,$*(F/F5F,F.F'F+F+F+F+F+*&F 0F+-%$cosGF8F+F+" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "y4 := u napply(rhs(%),x);\n" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#y4Gf*6#%\"xG 6\"6$%)operatorG%&arrowGF(,&*&#\"\"\"\"\"$F/*(F0#F/\"\"#-%$expG6#,$*&# F/F3F/9$F/!\"\"F/-%$sinG6#,$*&F2F/*&F0F2F:F/F/F/F/F/F/*&F4F/-%$cosGF>F /F/F(F(F(" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 101 "plot(y4(x),x= 0..10,labels=[\"x\",\"y4\"],title=\"The solution to y'' + y' + y = 0, \+ y(0) = 1, y'(0) = 0\");\n" }}{PARA 13 "" 1 "" {GLPLOT2D 485 485 485 {PLOTDATA 2 "6&-%'CURVESG6$7Y7$$\"\"!F)$\"\"\"F)7$$\"3ALL$3FWYs#!#>$\" 3I#ygq'=K'***!#=7$$\"3WmmmT&)G\\aF/$\"3olzdzAU&)**F27$$\"3m****\\7G$R< )F/$\"33$3o>K.v'**F27$$\"3GLLL3x&)*3\"F2$\"3ExE>(ymF%**F27$$\"3$***** \\ilyM;F2$\"3a>z:lekt)*F27$$\"3emmm;arz@F2$\"3e\"3@A\"Gmz(*F27$$\"3v** *\\7y%*z7$F2$\"3QMm-v[bh&*F27$$\"3[LL$e9ui2%F2$\"3Cv^lF3@\"G*F27$$\"3n mmm\"z_\"4iF2$\"3]!QIF-'Rk%)F27$$\"39ommT&phN)F2$\"3CXlg^B!>X(F27$$\"3 KLLe*=)H\\5!#<$\"3&\\6D'>\"*eKjF27$$\"3smm\"z/3uC\"F]o$\"3/XO.[-Da_F27 $$\"3!****\\7LRDX\"F]o$\"3%fUoQ'[!f9%F27$$\"3%om;zR'ok;F]o$\"31JQvaN%4 0$F27$$\"33++D1J:w=F]o$\"3YV2Psq4X?F27$$\"3oLLL3En$4#F]o$\"3:Qm8=`(\\7 \"F27$$\"3#pmmT!RE&G#F]o$\"3p5[/#oBtB%F/7$$\"3D+++D.&4]#F]o$!3k**oX2R* >O#F/7$$\"3;+++vB_(3pBL;Q6F 27$$\"3&om;z*ev:JF]o$!3A(>c%*fc3Q\"F27$$\"3_LLL347TLF]o$!3I!y0(yQ2d:F2 7$$\"3nLLLLY.KNF]o$!3qb6n*z`Ei\"F27$$\"33++D\"o7Tv$F]o$!3hW[K[x$yh\"F2 7$$\"3?LLL$Q*o]RF]o$!3;y*e\\+hVb\"F27$$\"3m++D\"=lj;%F]o$!3p=Nz_&\\cV \"F27$$\"3S++vV&RY2aF]o$!3SEsuENugSF/ 7$$\"3Znm;zXu9cF]o$!3EU?,(y%4VDF/7$$\"34+++]y))GeF]o$!3O$o&G1Hbq6F/7$$ \"3H++]i_QQgF]o$!35'GdN;C!=P!#@7$$\"3b++D\"y%3TiF]o$\"3-,1W&Q7*o&)!#?7 $$\"3+++]P![hY'F]o$\"30#z[z/)e?;F/7$$\"3iKLL$Qx$omF]o$\"3Uk:5,Y^7@F/7$ $\"3Y+++v.I%)oF]o$\"3InP?6Et_CF/7$$\"3?mm\"zpe*zqF]o$\"3by[SH:![h#F/7$ $\"3;,++D\\'QH(F]o$\"3KeeYx=.cEF/7$$\"3%HL$e9S8&\\(F]o$\"3Ef::cee(e#F/ 7$$\"3s++D1#=bq(F]o$\"3wWl)e)GcGCF/7$$\"3\"HLL$3s?6zF]o$\"3T:\"fikE(3A F/7$$\"3a***\\7`Wl7)F]o$\"3(f(yY\\<]K>F/7$$\"3enmmm*RRL)F]o$\"3>_xwmw \"4k\"F/7$$\"3%zmmTvJga)F]o$\"3qFrgL]jL8F/7$$\"3]MLe9tOc()F]o$\"3s?G'f quO.\"F/7$$\"31,++]Qk\\*)F]o$\"3M#fVk)Gl " 0 "" {MPLTEXT 1 0 59 "ode3 :=5*diff(y(x ),x,x)+10*sqrt(2)*diff(y(x),x)+10*y(x)=0;\n" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%ode3G/,(*&\"\"&\"\"\"-%%diffG6$-%\"yG6#%\"xG-%\"$G6$ F0\"\"#F)F)*(\"#5F)F4#F)F4-F+6$F-F0F)F)*&F6F)F-F)F)\"\"!" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 38 "dsolve(\{ode3,y(0)=0,D(y)(0)=1\},y( x));\n" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"yG6#%\"xG*&-%$expG6#,$* &\"\"##\"\"\"F.F'F0!\"\"F0F'F0" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "y5 := unapply(rhs(%),x);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>% #y5Gf*6#%\"xG6\"6$%)operatorG%&arrowGF(*&-%$expG6#,$*&\"\"##\"\"\"F29$ F4!\"\"F4F5F4F(F(F(" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "plot (y5(x),x=0..10);" }}{PARA 13 "" 1 "" {GLPLOT2D 485 485 485 {PLOTDATA 2 "6%-%'CURVESG6$7go7$$\"\"!F)F(7$$\"3gmmTN@Ki8!#>$\"3j$*[!3mEjL\"F-7$ $\"3ALL$3FWYs#F-$\"3qo(y+Xa;i#F-7$$\"3%)***\\iSmp3%F-$\"3%\\!*4&)3Vu&Q F-7$$\"3WmmmT&)G\\aF-$\"3OD)3\"e^6X]F-7$$\"3m****\\7G$R<)F-$\"3sE$3-_> ;G(F-7$$\"3GLLL3x&)*3\"!#=$\"3:I$fsdQ=M*F-7$$\"3YmmTN@Ki8FG$\"3?]UZx.f B6FG7$$\"3$*****\\ilyM;FG$\"3E:m(3DVtH\"FG7$$\"3RLLe*)4D2>FG$\"3Y[<=^_ Nc9FG7$$\"3emmm;arz@FG$\"3#z;XH^#\\,;FG7$$\"3;L$e*)4bQl#FG$\"3)[^Ygb(R B=FG7$$\"3v***\\7y%*z7$FG$\"3kyFG$\"3q&e]%Q2s(e#FG7$$\"3_om\"zW#z(3)FG$\"3'\\>$o; R&od#FG7$$\"39ommT&phN)FG$\"3=%*yHVs?jDFG7$$\"3A,+v=ddC%*FG$\"3%pV0$>j _&[#FG7$$\"3KLLe*=)H\\5!#<$\"3Uu9Kc_AzBFG7$$\"3-++v=JN[6Ffs$\"3mVi\"RE mME#FG7$$\"3smm\"z/3uC\"Ffs$\"3V!o8p_/x\"Ffs$ \"3OPEs)>2yW\"FG7$$\"33++D1J:w=Ffs$\"3'o-tgy#=@8FG7$$\"3+n;HdG\"\\)>Ff s$\"3qO*Q$z$)\\)>\"FG7$$\"3oLLL3En$4#Ffs$\"3)y6HHFUR3\"FG7$$\"3#pmmT!R E&G#Ffs$\"3nwwAF_&Fbz7$$\"3emm;/T1&*\\Ffs$\"3i[Tk))Q:sUF bz7$$\"3=nm\"zRQb@&Ffs$\"35P:+x^#eE$Fbz7$$\"3:++v=>Y2aFfs$\"3KnvvgZ8\" e#Fbz7$$\"3Znm;zXu9cFfs$\"33NsvBF6**>Fbz7$$\"34+++]y))GeFfs$\"3eU1xSg4 L:Fbz7$$\"3H++]i_QQgFfs$\"3Q^+/Pv&4=\"Fbz7$$\"3b++D\"y%3TiFfs$\"3Ag@$p ShQ;*!#@7$$\"3+++]P![hY'Ffs$\"3v'*z\"*e661pF`]l7$$\"3iKLL$Qx$omFfs$\"3 ;AeON>f]`F`]l7$$\"3Y+++v.I%)oFfs$\"3KG`\\3YGqSF`]l7$$\"3?mm\"zpe*zqFfs $\"3%4NhY)R8uJF`]l7$$\"3;,++D\\'QH(Ffs$\"3xJY80+V;CF`]l7$$\"3%HL$e9S8& \\(Ffs$\"3A'oxI'4,o=F`]l7$$\"3s++D1#=bq(Ffs$\"3jTH1Y;AE9F`]l7$$\"3\"HL L$3s?6zFfs$\"3MUA[XUq%4\"F`]l7$$\"3a***\\7`Wl7)Ffs$\"3GC,!f4MGH)!#A7$$ \"3enmmm*RRL)Ffs$\"33)G=ZB8EM'F^`l7$$\"3%zmmTvJga)Ffs$\"3:1()3zid=[F^` l7$$\"3]MLe9tOc()Ffs$\"3c'RSUs]om$F^`l7$$\"31,++]Qk\\*)Ffs$\"3i#G.J*4Y ^GF^`l7$$\"3![LL3dg6<*Ffs$\"3d%Hj=hfh8#F^`l7$$\"3%ymmmw(Gp$*Ffs$\"3/]( G:WK!\\;F^`l7$$\"3C++D\"oK0e*Ffs$\"3U0FnS]u]7F^`l7$$\"35,+v=5s#y*Ffs$ \"3%*y'y`m\"G&f*!#B7$$\"#5F)$\"3%o8np_TN@(Fgbl-%'COLOURG6&%$RGBG$Fjbl! \"\"F(F(-%+AXESLABELSG6$Q\"x6\"Q!Fgcl-%%VIEWG6$;F(Fibl%(DEFAULTG" 1 2 0 1 10 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1" }}}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 38 "solve(5*r^2 + 10*sqrt(2)*r + 10=0,r);\n" }}{PARA 11 "" 1 "" {XPPMATH 20 "6$,$*$\"\"##\"\"\"F%!\"\" F#" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 31 "At most one relative extrem um. " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "# end of section" } }}}{SECT 1 {PARA 3 "" 0 "" {TEXT -1 24 "2. Mechanical Vibrations" }} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 37 "restart; # clears previous \+ variables" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 86 "Here we consider the differential equation describing the motio n of a mass on a spring" }}{PARA 256 "" 0 "" {TEXT -1 8 "m u'' + " } {XPPEDIT 18 0 "gamma;" "6#%&gammaG" }{TEXT -1 16 " u' + k u = F(t)" }} {PARA 0 "" 0 "" {TEXT -1 74 "Here, u is the diplacement from equilibri um, m is the mass of the object, " }{XPPEDIT 18 0 "gamma;" "6#%&gammaG " }{TEXT -1 136 " is the damping constant, k is the spring constant, a nd F(t) is some external forcing term. All of the constants are positi ve: a larger " }{XPPEDIT 18 0 "gamma;" "6#%&gammaG" }{TEXT -1 66 " imp lies greater damping and a larger k implies a stiffer spring. " }} {PARA 0 "" 0 "" {TEXT -1 4 "If " }{XPPEDIT 18 0 "gamma" "6#%&gammaG" }{TEXT -1 31 " > 0, the motion is considered " }{TEXT 284 6 "damped" } {TEXT -1 7 ", if " }{XPPEDIT 18 0 "gamma" "6#%&gammaG" }{TEXT -1 19 " = 0 the motion is " }{TEXT 283 9 "undamped." }}{PARA 0 "" 0 "" {TEXT -1 1 " " }}{PARA 0 "" 0 "" {TEXT -1 90 "We can define the left s ide (lhs) of this ordinary differential equation (ode) as follows " }} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 56 "lhsode := m*diff(u(t),t$2) + gam*diff(u(t),t) + k*u(t);\n" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%'lh sodeG,(*&%\"mG\"\"\"-%%diffG6$-%\"uG6#%\"tG-%\"$G6$F/\"\"#F(F(*&%$gamG F(-F*6$F,F/F(F(*&%\"kGF(F,F(F(" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 72 "We can get a look at the general solution of the homogeneous equation by" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "dsolve(lhsode=0,u(t) );\n" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"uG6#%\"tG,&*&%$_C1G\"\"\" -%$expG6#,$**\"\"#!\"\"%\"mGF2,&%$gamGF2*$,&*$)F5F1F+F+*(\"\"%F+%\"kGF +F3F+F2#F+F1F+F+F'F+F+F+F+*&%$_C2GF+-F-6#,$**F1F2,&F5F+F6F+F+F3F2F'F+F 2F+F+" }}}{PARA 0 "" 0 "" {TEXT -1 205 "This is a little messy but you can see how the parameters m, gam, and k effect our solution. If the term under the radical sign is negative the solution will contain sin e and cosine terms of the form cos(" }{XPPEDIT 18 0 "omega;" "6#%&omeg aG" }{TEXT -1 12 " t) and sin(" }{XPPEDIT 18 0 "omega;" "6#%&omegaG" } {TEXT -1 10 " t) where " }{XPPEDIT 18 0 "omega;" "6#%&omegaG" }{TEXT -1 3 " = " }{XPPEDIT 18 0 "sqrt(abs(gam^2-4*mk))/(2*m)" "6#*&-%%sqrtG6 #-%$absG6#,&*$%$gamG\"\"#\"\"\"*&\"\"%F.%#mkGF.!\"\"F.*&F-F.%\"mGF.F2 " }{TEXT -1 18 ". In this case, " }{XPPEDIT 18 0 "omega;" "6#%&omega G" }{TEXT -1 23 " (omega) is called the " }{TEXT 286 9 "frequency" } {TEXT -1 31 " of the sine and cosine terms. " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }{TEXT 285 30 "Damped Oscillat or, no Forcing:" }}{EXCHG {PARA 0 "" 0 "" {TEXT -1 55 "Suppose m=gamma =k=1 and F(t) = 0. We assign these by " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 7 "m := 1:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 9 "g am := 1:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 7 "k := 1:" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "F := t -> 0:" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 73 "Let's take a look at the general solution of this homogeneous equation. " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "dsolve(lhsode=F(t),u(t));\n" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"uG6#%\"tG,&*(%$_C1G\"\"\"-%$expG6#,$*&\"\"#!\"\"F' F+F2F+-%$sinG6#,$*(F1F2\"\"$#F+F1F'F+F+F+F+*(%$_C2GF+F,F+-%$cosGF5F+F+ " }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 191 "Suppose that we displace th e mass 1 unit beyond equilibrium and let it go. The initial condition s corresponding to this statement are u(0)=1 and u'(0) = 0. These can \+ be assigned values by " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 9 "u o := 1: " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "uprimeo := 0:" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 72 "Now we can find and plot the so lution to the initial value problem by " }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 63 "sol := dsolve(\{lhsode = 0,u(0) = uo, D(u)(0) = upr imeo\},u(t));\n" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$solG/-%\"uG6#%\" tG,&*&#\"\"\"\"\"$F-*(F.#F-\"\"#-%$expG6#,$*&F1!\"\"F)F-F7F--%$sinG6#, $*(F1F7F.F0F)F-F-F-F-F-*&F2F--%$cosGF:F-F-" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 27 "u1 := unapply(rhs(sol),t);\n" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#u1Gf*6#%\"tG6\"6$%)operatorG%&arrowGF(,&*&#\"\"\"\" \"$F/*(F0#F/\"\"#-%$expG6#,$*&#F/F3F/9$F/!\"\"F/-%$sinG6#,$*&F2F/*&F0F 2F:F/F/F/F/F/F/*&F4F/-%$cosGF>F/F/F(F(F(" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 101 "plot(u1(t),t=0..10,labels=[\"t\",\"u(t)\"],title=\"t he solution to u'' + 2u' + 5u = 0, u(0)=1, u'(0)=0\");\n" }}{PARA 13 " " 1 "" {GLPLOT2D 485 485 485 {PLOTDATA 2 "6&-%'CURVESG6$7Y7$$\"\"!F)$ \"\"\"F)7$$\"3ALL$3FWYs#!#>$\"3I#ygq'=K'***!#=7$$\"3WmmmT&)G\\aF/$\"3o lzdzAU&)**F27$$\"3m****\\7G$R<)F/$\"33$3o>K.v'**F27$$\"3GLLL3x&)*3\"F2 $\"3ExE>(ymF%**F27$$\"3$*****\\ilyM;F2$\"3a>z:lekt)*F27$$\"3emmm;arz@F 2$\"3e\"3@A\"Gmz(*F27$$\"3v***\\7y%*z7$F2$\"3QMm-v[bh&*F27$$\"3[LL$e9u i2%F2$\"3Cv^lF3@\"G*F27$$\"3nmmm\"z_\"4iF2$\"3]!QIF-'Rk%)F27$$\"39ommT &phN)F2$\"3CXlg^B!>X(F27$$\"3KLLe*=)H\\5!#<$\"3&\\6D'>\"*eKjF27$$\"3sm m\"z/3uC\"F]o$\"3/XO.[-Da_F27$$\"3!****\\7LRDX\"F]o$\"3%fUoQ'[!f9%F27$ $\"3%om;zR'ok;F]o$\"31JQvaN%40$F27$$\"33++D1J:w=F]o$\"3YV2Psq4X?F27$$ \"3oLLL3En$4#F]o$\"3:Qm8=`(\\7\"F27$$\"3#pmmT!RE&G#F]o$\"3p5[/#oBtB%F/ 7$$\"3D+++D.&4]#F]o$!3k**oX2R*>O#F/7$$\"3;+++vB_(3pBL;Q6F27$$\"3&om;z*ev:JF]o$!3A(>c%*fc3Q\"F27$$ \"3_LLL347TLF]o$!3I!y0(yQ2d:F27$$\"3nLLLLY.KNF]o$!3qb6n*z`Ei\"F27$$\"3 3++D\"o7Tv$F]o$!3hW[K[x$yh\"F27$$\"3?LLL$Q*o]RF]o$!3;y*e\\+hVb\"F27$$ \"3m++D\"=lj;%F]o$!3p=Nz_&\\cV\"F27$$\"3S++vV&RY2aF]o$!3SEsuENugSF/7$$\"3Znm;zXu9cF]o$!3EU?,(y%4VDF/7$$\"34 +++]y))GeF]o$!3O$o&G1Hbq6F/7$$\"3H++]i_QQgF]o$!35'GdN;C!=P!#@7$$\"3b++ D\"y%3TiF]o$\"3-,1W&Q7*o&)!#?7$$\"3+++]P![hY'F]o$\"30#z[z/)e?;F/7$$\"3 iKLL$Qx$omF]o$\"3Uk:5,Y^7@F/7$$\"3Y+++v.I%)oF]o$\"3InP?6Et_CF/7$$\"3?m m\"zpe*zqF]o$\"3by[SH:![h#F/7$$\"3;,++D\\'QH(F]o$\"3KeeYx=.cEF/7$$\"3% HL$e9S8&\\(F]o$\"3Ef::cee(e#F/7$$\"3s++D1#=bq(F]o$\"3wWl)e)GcGCF/7$$\" 3\"HLL$3s?6zF]o$\"3T:\"fikE(3AF/7$$\"3a***\\7`Wl7)F]o$\"3(f(yY\\<]K>F/ 7$$\"3enmmm*RRL)F]o$\"3>_xwmw\"4k\"F/7$$\"3%zmmTvJga)F]o$\"3qFrgL]jL8F /7$$\"3]MLe9tOc()F]o$\"3s?G'fquO.\"F/7$$\"31,++]Qk\\*)F]o$\"3M#fVk)Gl< xF^w7$$\"3![LL3dg6<*F]o$\"3W.%*er(*fj\\F^w7$$\"3%ymmmw(Gp$*F]o$\"3-RE< 0A=yFF^w7$$\"3C++D\"oK0e*F]o$\"3\")=S'f@pas(Fhv7$$\"35,+v=5s#y*F]o$!3_ I " 0 "" {MPLTEXT 1 0 43 "restart: #clears all variable defin itions. " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 37 "lhsode := m*dif f(u(t),t$2) + k*u(t);\n" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%'lhsodeG, &*&%\"mG\"\"\"-%%diffG6$-%\"uG6#%\"tG-%\"$G6$F/\"\"#F(F(*&%\"kGF(F,F(F (" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 72 "We can get a look at the gen eral solution of the homogeneous equation by" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "dsolve(lhsode=0,u(t));\n" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"uG6#%\"tG,&*&%$_C1G\"\"\"-%$sinG6#*(%\"mG#!\"\"\" \"#%\"kG#F+F3F'F+F+F+*&%$_C2GF+-%$cosGF.F+F+" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 7 "m := 1:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 7 "k := 4:" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 53 "The solution to t he homogeneous equation is found by " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "dsolve(lhsode=0,u(t));\n" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"uG6#%\"tG,&*&%$_C1G\"\"\"-%$sinG6#,$*&\"\"#F+F'F+F +F+F+*&%$_C2GF+-%$cosGF.F+F+" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "F := t -> sin(3*t): " }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 138 "No tice: The frequency of the forcing term is 3 and the natural frequency is 2. Suppose we start the spring from the resting position by " }} }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "uo := 0:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "uprimeo := 0:" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 22 "and find the solution " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 60 "sol := dsolve(\{lhsode=F(t),u(0)=uo, D(u)(0)=uprimeo \},u(t));\n" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$solG/-%\"uG6#%\"tG,& *&#\"\"$\"#5\"\"\"-%$sinG6#,$*&\"\"#F/F)F/F/F/F/*&#F/\"\"&F/-F16#,$*&F -F/F)F/F/F/!\"\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "u2 := u napply( rhs(sol), t);\n" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#u2Gf*6#% \"tG6\"6$%)operatorG%&arrowGF(,&*&#\"\"$\"#5\"\"\"-%$sinG6#,$*&\"\"#F1 9$F1F1F1F1*&#F1\"\"&F1-F36#,$*&F/F1F8F1F1F1!\"\"F(F(F(" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "plot(u2(t),t=0..60);\n" }}{PARA 13 "" 1 "" {GLPLOT2D 485 485 485 {PLOTDATA 2 "6%-%'CURVESG6$7[dm7$$\"\"!F )F(7$$\"3')*****\\7t&pK!#=$\"3RiHUP!*GH;!#>7$$\"3s******\\i9RlF-$\"3Qs ,PZ?xZ5F-7$$\"3#*****\\7G$R<)F-$\"3'\\N28c\"p@(3)*F -$\"3[\\wc;F07$$\"3 1++voozw=FC$!3]J\">CyFC$!3S#fK^xgL?\"F-7$$\"3$* *\\PMxQ!>?FC$!3+I:,&f9^*=F-7$$\"33]7yD(f,4#FC$!3]vORZa')fDF-7$$\"3C+]7 y1Gh@FC$!3#G3KGf\\N<$F-7$$\"3S](o/j,CB#FC$!3hq4L%y0Jr$F-7$$\"36+D\"GeA NI#FC$!3yf^MuZUdTF-7$$\"3#)\\i:NNkuBFC$!3ilh:s(Q#)[%F-7$$\"3)*****\\([ kdW#FC$!3u*)Qc%)=\"4p%F-7$$\"3+v$4@'fv&[#FC$!3uD`%=t;Xu%F-7$$\"3-](=nV Zd_#FC$!3KInR`-1`ZF-7$$\"3/D\"G8\"*Qdc#FC$!3)RS(RR`&er%F-7$$\"31+v$fQI dg#FC$!3%4C`a=GFj%F-7$$\"35]i:NLr&o#FC$!3;&4W1191L%F-7$$\"39+]P%G'plFF C$!3CO\"4Jz&pbQF-7$$\"3=]PfL#zc%GFC$!3'oY\"Q#=H\\A$F-7$$\"3B+D\"G=ic#H FC$!3muA'eDQDY#F-7$$\"3G]7.K^k0IFC$!3Cxd:y:$*)f\"F-7$$\"3))***\\73Gc3$ FC$!3oVgp7&))Hp'F07$$\"3!\\Pfeb>c7$FC$!3yhl[z>C;>F07$$\"3[\\(o/.6c;$FC $\"3!pTbp_c%F-7$$\"3y* ***\\(=Wv/%FC$\"3O'y>XXkbt$F-7$$\"3++]i:t0GTFC$\"3(Ht,:QpI7$F-7$$\"3A+ +vV/d3UFC$\"3v:&4yS1\"=CF-7$$\"3W+](=d$3*G%FC$\"3g==*H5'Ra;F-7$$\"3m++ ++nfpVFC$\"3igEPzRdp')F07$$\"3)3+D\"G)4,X%FC$\"3OAN?K3cJ!*!#?7$$\"36,+ DcHiIXFC$!3j+bzzY7KkF07$$\"3L,]P%3O6h%FC$!3q6=p2V108F-7$$\"3m++]7#\\;p %FC$!3W\"pO/3F<(=F-7$$\"3A++voan_[FC$!3w#)\\>\"4&>cEF-7$$\"3m+++DADw#F-7$$\"3.+]7yNAM`FC$!3I P=-K#))*fAF-7$$\"3g+vo/X[%\\&FC$!37UF#*)3-)*f\"F-7$$\"3G++DJauacFC$!3( zt34Jl]^*F07$$\"3a+]P%Gn_(fFC$!3]8toO%H=P\"F07$$\"3!*****\\P\"*y&H'FC$ \"3oW))zjq*4+\"!#B7$$\"3w*****\\#R&Hf'FC$\"3gl8yj'ecR\"F07$$\"3i****\\ 7(=,*oFC$\"3n$3]oo#*4t)F07$$\"3c***\\i5,(QqFC$\"3SG.2R))\\f9F-7$$\"3[* ******\\$G(=(FC$\"3y\\m$Qw#R$3#F-7$$\"3U***\\P*e'eL(FC$\"3*)p#H(\\zO9E F-7$$\"3O****\\(G[W[(FC$\"3_x**HZ_i7HF-7$$\"3f**\\7./\"H_(FC$\"3]=XiF-7$$\"3a*****\\n$*f%zFC$\"3wtw=2^Kv8F-7$$\"37***\\i!z\"H-) FC$\"371o$3PQ:a(F07$$\"3e****\\P@%)*4)FC$\"3#>Dj%pOKwjFbz7$$\"3G)**\\( ojww\")FC$!31x?*R&\\Z;nF07$$\"3y)******f!p`#)FC$!3OC9^`&3XU\"F-7$$\"3C ***\\7$[hI$)FC$!3OzVIm92l@F-7$$\"3s****\\i!RvS)FC$!33y/9hhuiGF-7$$\"3? ++v$HjW[)FC$!3!4,uU$z^([$F-7$$\"3m+++DvQh&)FC$!3&>k2(RSP6SF-7$$\"39,+D c$oYF- 7$$\"3y++]P9!R&*)FC$!3Zes4Tx!3S%F-7$$\"3\"*****\\(ecM.*FC$!3%)=(y26Z+' RF-7$$\"3/****\\P<,8\"*FC$!3_?LUIBFhLF-7$$\"3%*****\\()oc#>*FC$!3$GB#) pl'*pi#F-7$$\"3'3++v.A@F*FC$!3g-P\\cU-'y\"F-7$$\"3)*****\\(=x;N*FC$!3w #)\\h*)>T@()F07$$\"3W++]iZX\"R*FC$!3Y$>e/OvR*RF07$$\"37****\\PBBJ%*FC$ \"3y()zgMD!\\u(Fbz7$$\"3e****\\7*45Z*FC$\"32z1N%Q)*\\`&F07$$\"3.++]([( y5&*FC$\"3K[D/Sc'Q-\"F-7$$\"3#4++vjU.f*FC$\"3yCu;QqPG>F-7$$\"31++](y(* )p'*FC$\"3A3M'yb>Uv#F-7$$\"3?****\\PHX\\(*FC$\"3+;3o#*\\JoMF-7$$\"35++ ](33!H)*FC$\"3zcZmJ()zUSF-7$$\"3,,+]PKc3**FC$\"3#zR6em!GcWF-7$$\"39++] (Q=\"))**FC$\"3YV_cJax%p%F-7$$\"3#*\\7`Rox-5!#;$\"3_LX2#>4hu%F-7$$\"3+ +DJS=u15Fj[m$\"3W&[[)*HhBv%F-7$$\"34]P4Toq55Fj[m$\"3j)zkyjRSr%F-7$$\"3 <+](=%=n95Fj[m$\"3M8=t]+9KYF-7$$\"3;+vVV=gA5Fj[m$\"3WAL)y'H#RM%F-7$$\" 3:+++X=`I5Fj[m$\"30'>2B\"=l/RF-7$$\"39+DcY=YQ5Fj[m$\"3&4>9!fnUPLF-7$$ \"37+]7[=RY5Fj[m$\"3cos$)HkCqEF-7$$\"37+vo\\=Ka5Fj[m$\"3O!>E0'*HX$>F-7 $$\"35++D^=Di5Fj[m$\"3Sk()Hj=]j6F-7$$\"35+D\"G&==q5Fj[m$\"30.)[J#yE0RF 07$$\"33+]Pa=6y5Fj[m$!3*Q@()*GBvX;uO 5F-7$$\"31++]d=(R4\"Fj[m$!3Wy1>K#\\vj\"F-7$$\"3.+]ig=$)46Fj[m$!3;&oTP* )yy^#F-7$$\"3,++vj=pD6Fj[m$!3U<0ERxx8HF-7$$\"33]iSEe+U6Fj[m$!3!po&[h4w PGF-7$$\"3)**\\i!*y>$e6Fj[m$!3;NK*>3klQ#F-7$$\"3')\\(=\"Fj[m$!3SSW<)>^#f5F-7$$\"32+voRcdB7Fj[m $!3Y3Z>FK*=o\"F07$$\"3/+++lN?c7Fj[m$!3Amik(ekI2%!#D7$$\"31]PfQ#)Qj7Fj[ m$\"3\"H/ix,')R`\"!#@7$$\"3))*\\(=7Hdq7Fj[m$\"3+gHYPXAO8Fbz7$$\"3!*\\7 y&edxF\"Fj[m$\"3G!pVDbubd%Fbz7$$\"3#***\\PfA%\\G\"Fj[m$\"338dv'*\\(f2 \"F07$$\"3%**\\ilg6$*H\"Fj[m$\"397[y#z?dW$F07$$\"3'****\\P&4o88Fj[m$\" 3g9]%R8FrY(F07$$\"3***\\P4I]!G8Fj[m$\"3i_.f]m](G\"F-7$$\"3,+]7['>CM\"F j[m$\"3m.`#*f?\\#*=F-7$$\"3/+DJ&**)yc8Fj[m$\"3SDGXsP(3X#F-7$$\"31++]U$ e6P\"Fj[m$\"3S*Qc\"yP!e$GF-7$$\"3-D1*GY-_P\"Fj[m$\"3-Q*Q;^ug*GF-7$$\"3 )*\\7G$eY#z8Fj[m$\"3T`mI!e@3$HF-7$$\"37v=n.2H$Q\"Fj[m$\"3X\"y$*fOe!QHF -7$$\"33+D1C[L(Q\"Fj[m$\"3LQ4-U35;HF-7$$\"3=]P%[1BaR\"Fj[m$\"3aV>X'y&e zFF-7$$\"36+]i08^.9Fj[m$\"3$*)\\Py`v]^#F-7$$\"3.]iSY&*f69Fj[m$\"3fE\\0 !=eH7#F-7$$\"38+v=(y(o>9Fj[m$\"3A1LDQDe5;F-7$$\"3B](oz-wxU\"Fj[m$\"3(G Un))*fNA**F07$$\"3:++voU'eV\"Fj[m$\"3Ck;I_C;()GF07$$\"33]7`4D&RW\"Fj[m $!3]a([*RY![t%F07$$\"3=+DJ]2/_9Fj[m$!3OtxXl6Lj7F-7$$\"3G]P4\"**G,Y\"Fj [m$!3cMoQfeuY?F-7$$\"3?+](=B<#o9Fj[m$!3Qhf)4,Z#)y#F-7$$\"37]ilsaIw9Fj[ m$!3ag^EBgx_MF-7$$\"3/+vV8PR%[\"Fj[m$!3UTf#eMGv+%F-7$$\"3'*\\(=U&>[#\\ \"Fj[m$!3i\\Lx=arBWF-7$$\"31+++&>q0]\"Fj[m$!3sLbs*Ra\"yYF-7$$\"3***\\P 4#4j/:Fj[m$!3w6%o6-(\\RZF-7$$\"35+](ok\"p3:Fj[m$!3=#[F+mFXv%F-7$$\"3/+ D\"GP_F^\"Fj[m$!3**eAxcRXAZF-7$$\"3(****\\()48o^\"Fj[m$!3HXe_%pFIk%F-7 $$\"3,+]i]X$\\_\"Fj[m$!3!*efLG)))QM%F-7$$\"30++]-g0L:Fj[m$!3!)Q+3T([h' QF-7$$\"3\"***\\PauEC?F-7$$\"3++]7e.Ud:Fj[m$! 3eNz')>=@u:F-7$$\"3%**\\iS3\"[h:Fj[m$!3VA01[\\N26F-7$$\"3')******4=al: Fj[m$!3!pSz'p&zlG'F07$$\"3)**\\Pf`-'p:Fj[m$!3#**ex'pyKK9F07$$\"3!***\\ (=EjOd\"Fj[m$\"3%[#=:$ofsV$F07$$\"3%)*\\7y)Rsx:Fj[m$\"33\")eO)p[+F)F07 $$\"3%****\\Pr%y\"e\"Fj[m$\"3Padn%pN9I\"F-7$$\"3!)**\\ilh!**e\"Fj[m$\" 3U9eodbo.AF-7$$\"3%)****\\RA;Fj[m$\"3o#>!pq5KuXF-7$$\"3$)*****\\U80j\"Fj[m$\"3jW$*4dMSTZF -7$$\"3')****\\sk;Y;Fj[m$\"3`Y#zZstg`%F-7$$\"3))******>&>=m\"Fj[m$\"3P +-j]s%eq$F-7$$\"33++vVgkp;Fj[m$\"3Gez47QF1JF-7$$\"3#*****\\nDZx;Fj[m$ \"3v]?*z7'3?CF-7$$\"3w***\\74*H&o\"Fj[m$\"3#)o[dACHy;F-7$$\"3%******\\ hDJp\"Fj[m$\"3Tj`;)4Z18*F07$$\"38++vQ@&4q\"Fj[m$\"3!o2rc1?=c\"F07$$\"3 (*****\\i'y(3F-7$$\"3/+](=n&f7=Fj[m$!3+2Gi/ZWE8F-7$$\"3')\\PfQ,\"o#=Fj[m$!3 H&H<7_G$RyF07$$\"3.+DJ0Y-T=Fj[m$!3hKq!**>xFj[m$\"3/w\"Qs]E+>\"F-7$$\"3/+vV$zn3(>Fj[m$\"3Q>i`wG#z*=F-7$$\"3 /](=#p,x()>Fj[m$\"3@9\"GhW$3RDF-7$$\"3.+++XDn/?Fj[m$\"3k/g[/(f'3HF-7$$ \"3))\\PMx@D3?Fj[m$\"3Dy=-P7/MHF-7$$\"33+vo4=$=,#Fj[m$\"3F'\\))Q$HiPHF -7$$\"3$*\\7.U9T:?Fj[m$\"3^u`Hb9E=HF-7$$\"3z**\\Pu5**=?Fj[m$\"3G`[.U>- vGF-7$$\"3%)*\\i!R.:E?Fj[m$\"3v;tSe'4Vr#F-7$$\"3!****\\Pg4L.#Fj[m$\"3# *)Rd&G:v_CF-7$$\"3&**\\P%o)o//#Fj[m$\"3g#*Gi$He<4#F-7$$\"3-+]7L\"Gw/#F j[m$\"3]`!))yGQqj\"F-7$$\"31+D\"yR(ya?Fj[m$\"3)*Ql.7yb)4\"F-7$$\"37++] im%>1#Fj[m$\"3vlFN`%oI!\\F07$$\"3=+v=Ff5p?Fj[m$!3y-0!Q:P7q\"F07$$\"3C+ ](=>li2#Fj[m$!3Vxvf!o$Q@')F07$$\"3%**\\ilXCM3#Fj[m$!3)[6\"y1Q)Gc\"F-7$ $\"3)****\\7s$e!4#Fj[m$!3m\\#y,)f6[AF-7$$\"3/+v$f)Hu(4#Fj[m$!336G868:$ *GF-7$$\"3u**\\i]A!\\5#Fj[m$!3Ok$e**phQZ$F-7$$\"3z*\\7`^h?6#Fj[m$!3/!Q 5'Q/nnRF-7$$\"3&)******z2A>@Fj[m$!3#z$3^!e0XN%F-7$$\"3ruoz1([v7#Fj[m$! 3wi`g?\"fzk%F-7$$\"3$*\\PfLm(e8#Fj[m$!3%RE@^5Z]v%F-7$$\"3:D1RgX?W@Fj[m $!3cDHdD^KlYF-7$$\"3-+v=([KD:#Fj[m$!3#oPFpoWyP%F-7$$\"3)[P%)RTg3;#Fj[m $!3=Ki^cDR,RF-7$$\"3u\\7yS$)=p@Fj[m$!3oH8MHg>aKF-7$$\"3gC\"yvE;v<#Fj[m $!3ynsbna5jCF-7$$\"3#)**\\P%>We=#Fj[m$!3%zdsM<'Hi:F-7$$\"3gPMxd\"3+>#F j[m$!31SC`?i,$3\"F-7$$\"3/v=<@@<%>#Fj[m$!3#HIk6')*f:fF07$$\"3Y7.d%3O$) >#Fj[m$!3!)R+u88mW$*Fbz7$$\"3!*\\(oz/+D?#Fj[m$\"3arUefw=dSF07$$\"3L(=n 8,km?#Fj[m$\"3Bo')>iU<.!*F07$$\"3wCcwuz#3@#Fj[m$\"37!*fiR'*y%Q\"F-7$$ \"3>iS;Q>*\\@#Fj[m$\"39=?YHkq`=F-7$$\"3i*\\i:!f:>AFj[m$\"36/))fu%R=I#F -7$$\"3\\u$f$GQ[FAFj[m$\"3w')pOTHC;JF-7$$\"3r\\i:b<\"eB#Fj[m$\"3aA%48] LEz$F-7$$\"3#\\7`>oRTC#Fj[m$\"3>R^%RoyEI%F-7$$\"3y***\\(3wY_AFj[m$\"3= lu_#4hki%F-7$$\"3AiS\">U`hD#Fj[m$\"3G2ZF-7$$\"3kC\"y]BR)fAFj[m$ \"39l1UbO$*[ZF-7$$\"3U(=U#[]_jAFj[m$\"3iNIlN`t^ZF-7$$\"3%)\\iSh3@nAFj[ m$\"3qs&)4H*Qgr%F-7$$\"3/vVt([#euAFj[m$\"3Y#pF$[uqKXF-7$$\"3))*\\iS6a> G#Fj[m$\"3[`3fWsZ4UF-7$$\"3uC1RSdK*G#Fj[m$\"3U\\tElc%>w$F-7$$\"3%*\\(= nO(p'H#Fj[m$\"3g*GI&RG,5KF-7$$\"39vo/$**oSI#Fj[m$\"3;en5rN&pd#F-7$$\"3 )***\\P>1W6BFj[m$\"3YXAQq>O))=F-7$$\"3$[7.dC7)=BFj[m$\"3wR%*oJ&=5<\"F- 7$$\"3/]7.sQ=EBFj[m$\"3/<8!yqBt^%F07$$\"3Cv$f$)\\bNL#Fj[m$!3wQQmu;)yV# F07$$\"34+voCr#4M#Fj[m$!3=j7$p[`)>*)F07$$\"3%\\i:5v)H[BFj[m$!31K`C\"fv BZ\"F-7$$\"3y\\PMx.nbBFj[m$!3+vt?&e1$o>F-7$$\"3ku=n.?/jBFj[m$!3AS_!3L; vO#F-7$$\"3%)******HOTqBFj[m$!39$HrbLsCm#F-7$$\"3wu=Us9]yBFj[m$!3cc#) \\HM3jGF-7$$\"3/]P%[J*e'Q#Fj[m$!3pM6XIF+QHF-7$$\"3'\\ils:xYR#Fj[m$!3#* =kC/tW'*GF-7$$\"3))*\\(o**\\w-CFj[m$!31f!)*f&\\+`FF-7$$\"3t\\7`%oS*=CF j[m$!3u:fr'=L%QAF-7$$\"3%***\\Ppj6NCFj[m$!3%3wV6b#ep:F-7$$\"39](=U0#H^ CFj[m$!3Y[9'ebKv>*F07$$\"3++D1RxYnCFj[m$!3AH5*[#Hy#=%F07$$\"3%)\\i!RUV O[#Fj[m$!3g9nwQE4G7F07$$\"3/++v3\">)*\\#Fj[m$!3*)G%>'f:m.7Fbz7$$\"30+v $f=A_^#Fj[m$\"3U)fe$o*)[&p$Fj]l7$$\"30+]7j_iIDFj[m$\"3Gs?Imj9hDFbz7$$ \"31+DJS$Gga#Fj[m$\"33t+-yKjP;F07$$\"31++]<9VhDFj[m$\"37\\C`.s\\)y%F07 $$\"31+vo%\\Mod#Fj[m$\"35f\">.S8Ay*F07$$\"31+](=dPAf#Fj[m$\"3YQ)fkplPg \"F-7$$\"31+D1\\1k2EFj[m$\"35AA5BeqRAF-7$$\"32++DEP/BEFj[m$\"3f%yWK*>V NFF-7$$\"3ED\"GQZz5j#Fj[m$\"3iOdb4.j')GF-7$$\"35]iS@_6REFj[m$\"3JE-C% \\)*)QHF-7$$\"3%\\P%)*o4:ZEFj[m$\"3g`@TWWXxGF-7$$\"3y*\\ilr'=bEFj[m$\" 3if2B+Ek#p#F-7$$\"3iC19kCAjEFj[m$\"3)**pPU#G$3Q#F-7$$\"3#)\\(=<@e7n#Fj [m$\"3`x6u(R**\\%>F-7$$\"3-voHfRHzEFj[m$\"3/j]DOr$\\R\"F-7$$\"3')**\\( oqHto#Fj[m$\"3S#HIJ@@4Z(F07$$\"3qCJXaaO&p#Fj[m$\"3%y+`k0g1S#Fbz7$$\"3! *\\7.-7S.FFj[m$!3#\\[SdZ+gY(F07$$\"34v$4'\\pV6FFj[m$!3'*zwMZ52L:F-7$$ \"3$**\\(=(ps%>FFj[m$!3d0))zB.J,BF-7$$\"3yCcwW%3vs#Fj[m$!3\\--()**3Y;I F-7$$\"3(*\\PM#>Wbt#Fj[m$!3)zNYF-7$$\"3=+Dc')H8\"y#Fj[m$!3CsD=dkskVF-7$$\"3ID\"G8J7&)y#Fj[m$!3=, RY(pR$F-7$$\"3>v$f3'4F.GFj[m$!3 UR%y[-b+s#F-7$$\"3K+]i&G]1\"GFj[m$!3;,z$Q=fJ&>F-7$$\"35D1R5'H!=GFj[m$! 3qCMgH$yw6\"F-7$$\"3B]i:N*3a#GFj[m$!3zwaGR_EGCF07$$\"3Nv=#*f#)yKGFj[m$ \"3xP(o=*f!fS'F07$$\"38+vo%en,%GFj[m$\"3C&4lJ!3U,:F-7$$\"3EDJX4paZGFj[ m$\"3-pU7#eD%4BF-7$$\"3Q](=UBE\\&GFj[m$\"3.?t())***\\OIF-7$$\"3;vV)*eb IiGFj[m$\"3M-i5c!)zdOF-7$$\"3G++v$)[opGFj[m$\"3g\\\"R(HSr_TF-7$$\"3)\\ (=U#[Ww(GFj[m$\"3HqUu7'zq_%F-7$$\"3/]P4\"3/c)GFj[m$\"3sd/P=5Yv%F-7$$\"35DcwzOc$*GFj[m$\"3DNGeFq#)RZF -7$$\"3ii:5zMa(*GFj[m$\"3K%oi$>$R0o%F-7$$\"3:+vVyK_,HFj[m$\"3W&HW0H\"* zd%F-7$$\"3@v$4r(G[4HFj[m$\"3iU^#)HjO]UF-7$$\"3F]7yvCW6)HFj [m$!3%R=>45uhy\"F-7$$\"33D\"GQ')y!*)HFj[m$!3;fK2'3P[D#F-7$$\"37++]i%Qq *HFj[m$!3WZUW5k9/EF-7$$\"3+Dc^[iI0IFj[m$!3.k)pe%eFNGF-7$$\"3))\\7`MSd8 IFj[m$!3Wt4^*e&)R$HF-7$$\"3]i!Rv#zqE#F-7$$\"3')** \\i]2=jIFj[m$!3MMbbT,O!e\"F-7$$\"3)*\\ilAjrzIFj[m$!3Q%>W]ux)e\"*F07$$ \"33+vo%*=D'4$Fj[m$!3e.^d[#G\"oSF07$$\"3%)\\(=nY(y7JFj[m$!36!yCH_v=8\" F07$$\"3'****\\(QIKHJFj[m$!3A2TgH&3b9*F\\cm7$$\"3C+D\"GI>**fSF07$$\"3&**\\i!fVH ,KFj[m$\"3Ey8qxc)*y$)F07$$\"3))**\\7B')o:KFj[m$\"3s*)oO-f^)R\"F-7$$\"3 \")*\\(=()G3IKFj[m$\"37o`x&*y.0?F-7$$\"3u***\\7:xWC$Fj[m$\"3_1#GUz68a# F-7$$\"3QCc,^-D_KFj[m$\"3W)Gg3i\\kv#F-7$$\"3s\\7y]L-gKFj[m$\"3E\\/!)*3 mW*GF-7$$\"3/iSm+*4RE$Fj[m$\"3UU#R-Z5#HHF-7$$\"33voa]kznKFj[m$\"3O$eB[ Wo'QHF-7$$\"3Q(oH/+$orKFj[m$\"3vFqCk*=8#HF-7$$\"3q*\\7.bpbF$Fj[m$\"3yh ;\")3f#f(GF-7$$\"3o\\P%)\\d6\"H$Fj[m$\"3!*>K;(z'\\+CF-7$$\"3n**\\P\\>m 1LFj[m$\"3#[h<\\'4[l9F-7$$\"3IC19\\]V9LFj[m$\"3-b%G1a/M])F07$$\"3l\\i! *[\"3AK$Fj[m$\"3fA!e@k3Dh\"F07$$\"3+v=n[7)*HLFj[m$!39;swk)HSx&F07$$\"3 j*\\P%[VvPLFj[m$!3Zf\"\\'z#ovL\"F-7$$\"3)\\7.#[u_XLFj[m$!3smD(['\\v)3# F-7$$\"3K](oza+LN$Fj[m$!3+!4H5;H%*z#F-7$$\"3'\\PMxkt5O$Fj[m$!3x\\4^j** \\QMF-7$$\"3K++]Zn%)oLFj[m$!3[hp[$[Co(RF-7$$\"35]i:?r(oP$Fj[m$!3c#>*eg H')*R%F-7$$\"3i+D\"G\\2\\Q$Fj[m$!39xUj$)y#[m%F-7$$\"3'eiS\"zE#*)Q$Fj[m $!3@am7oogKZF-7$$\"3S](oa'y$HR$Fj[m$!3Mz')[H@GbZF-7$$\"3%\\(oz^I&pR$Fj [m$!3-4L*f(***>t%F-7$$\"3?+]7Q#o4S$Fj[m$!35h*fNxCCm%F-7$$\"3**\\7y5')* *3MFj[m$!3C]w\")p9p&Q%F-7$$\"3z*\\PM)*GqT$Fj[m$!3:H_->'RH$RF-7$$\"3I]P 4c$f]U$Fj[m$!3/Im5.P??LF-7$$\"34++vG(*3LMFj[m$!31i([L`H5d#F-7$$\"3*)\\ iS,,7TMFj[m$!37Z;Cj\">ar\"F-7$$\"3S+D1u/:\\MFj[m$!3]T(Q\\e&R%)yF07$$\" 3mD1Rgc;`MFj[m$!3z(4'o()>E,JF07$$\"3?](=n%3=dMFj[m$\"3d;\"H/1MVr\"F07$ $\"3uuo/Lg>hMFj[m$\"3<&GWd)R)>^'F07$$\"3++]P>7@lMFj[m$\"3Z122$[\\T7\"F -7$$\"3y\\7.#fTKZ$Fj[m$\"3q\"foW9,+.#F-7$$\"3e*\\(ok>F\"[$Fj[m$\"3'yBz lpA:&GF-7$$\"33]PMPBI*[$Fj[m$\"3%4lk$z=FbNF-7$$\"3))******4FL(\\$Fj[m$ \"37XW4*)3S8TF-7$$\"3n\\(oH'))=0NFj[m$\"3%=Kf0z:$)\\%F-7$$\"3=+v$f,XI^ $Fj[m$\"3lOfX9VL6ZF-7$$\"3Wv=U#4tp^$Fj[m$\"3Kwc]8dd^ZF-7$$\"3)*\\i!*o6 !4_$Fj[m$\"3QZSpuHkZZF-7$$\"3_C1RX#H[_$Fj[m$\"3X3&\\LBR,q%F-7$$\"3y** \\(=Kd(GNFj[m$\"3W&[iI`M,h%F-7$$\"3H]P%[Z8m`$Fj[m$\"3(=y0)Gt74VF-7$$\" 33+D\"yipWa$Fj[m$\"3)pJ!>syshQF-7$$\"3))\\7y!yDBb$Fj[m$\"3;!QDGft3H$F- 7$$\"3Q++vL>=gNFj[m$\"3P68;^*3Ti#F-7$$\"3=](=n3Q!oNFj[m$\"3qIGQ!\\h@*= F-7$$\"3)**\\(oRU*ed$Fj[m$\"3HU3`?TRF6F-7$$\"3[]il#R]Pe$Fj[m$\"3s\")pt J[:AOF07$$\"3G+]iXlg\"f$Fj[m$!38rz8AjoCPF07$$\"33]Pf)pi%*f$Fj[m$!3su8Y #>$o[5F-7$$\"3e+Dc^)=tg$Fj[m$!3)QQKi6>Ek\"F-7$$\"3Q]7`/]<:OFj[m$!3$*H& y?9:c8#F-7$$\"3=++]d6.BOFj[m$!3C^9CB.%[^#F-7$$\"3MDJX*QK1j$Fj[m$!3]2E9 C6;nFF-7$$\"3^]iS@OBQOFj[m$!3'3YPH7 1f0S\"HF-7$$\"39+DJ&3OMl$Fj[m$!3E6K?fD(z'GF-7$$\"3w\\(=#\\&Q'oOFj[m$!3 #)y#=^?x()[#F-7$$\"33+]785%Qo$Fj[m$!3ElaNmm`/>F-7$$\"3S]7.xM/*p$Fj[m$! 3c\"Gh%*>*>l7F-7$$\"3-+v$4%fC9PFj[m$!3N0ik4'[J,(F07$$\"3l\\P%[S[%HPFj[ m$!3%\\$p')[ShtHF07$$\"3)****\\(o3lWPFj[m$!3D()[?iK**HxFbz7$$\"3S(o/[! 3()[PFj[m$!3%zV)p')[oCXFbz7$$\"3#[Pf3u!4`PFj[m$!3Fbpn9E-OBFbz7$$\"3CiS \"pn5tv$Fj[m$!3jbQ2,ED+**F\\cm7$$\"3o\\(oHhI:w$Fj[m$!3_hH#ow\"fHHF\\cm 7$$\"3_C\"y][q*pPFj[m$\"3([s'G#ozE/\"!#F7$$\"33+v=d.TyPFj[m$\"3I\")QPD 1MbIF\\cm7$$\"3[]iS,,H&z$Fj[m$\"3L3cbwRHOyFbz7$$\"3<+]iX)p@\"QFj[m$\"3 ]>F9Uu)QN$F07$$\"3')\\P%)*e\\!HQFj[m$\"37y,Eg,#==)F07$$\"3c*\\iSLHf%QF j[m$\"3g.oI:N8z9F-7$$\"3C\\7Gy!4G'QFj[m$\"3O4cwX(>E=#F-7$$\"3m****\\A) )ozQFj[m$\"38'R\"\\4miNFF-7$$\"3?\\(oHUss)QFj[m$\"39-GOA**p!)GF-7$$\"3 Y*\\PM-c[*QFj[m$\"3.'ejL:s%QHF-7$$\"3s\\i!RiRC!RFj[m$\"3mnUzwq='*GF-7$ $\"3)***\\PCK-5RFj[m$\"3YInhc`1XFF-7$$\"3D]P%[#ogA,?'35[#F -7$$\"3!)*\\7`U!>DRFj[m$\"3f')yt)f$30@F-7$$\"3N\\7yDSxKRFj[m$\"3$elk17 bPi\"F-7$$\"3i***\\iid.%RFj[m$\"3un<(>L\\)[5F-7$$\"3))\\(=nATz%RFj[m$ \"3Q!z*p*RAC(RF07$$\"3W*\\(=F[_bRFj[m$!3'*RSuY@<(4$F07$$\"3)*[ilF%3J'R Fj[m$!3o^8AYl-Z5F-7$$\"3C**\\7G?pqRFj[m$!3Uw%QFzbqy\"F-7$$\"3]\\PfGcFy RFj[m$!3v:Zy6Mt+DF-7$$\"3y*\\i!H#fe)RFj[m$!3o)yN:NO)eJF-7$$\"3-]7`HGW$ *RFj[m$!3B8!G5fALt$F-7$$\"3e******Hk-,SFj[m$!33rGke9k)>%F-7$$\"3#)*\\( =UN74SFj[m$!3?.t4MMS]XF-7$$\"3O**\\Pa1A^TSFj[m$!3Arv_CHV(pdJHF-7$$\"3S****\\FL!e1%Fj[m$!3%*eV@tI z8@F-7$$\"3=\\Pf$)=&)pSFj[m$!3=,h@O(R\"p;F-7$$\"3k*\\(oR/!R2%Fj[m$!3\\ IUh\"[zn?\"F-7$$\"37]7y&**[z2%Fj[m$!30\\Midj'fJ(F07$$\"3!***\\(=b(*>3% Fj[m$!39[4y>'>j[#F07$$\"3m\\(oz5Yg3%Fj[m$\"3_iJ:@wypBF07$$\"3U*\\iSm%4 !4%Fj[m$\"3#[?PXQu1?(F07$$\"3=\\i:?K9%4%Fj[m$\"3#>w;\\x'\\&>\"F-7$$\"3 m***\\ix\">)4%Fj[m$\"31;\"f%**HAe;F-7$$\"3\"**\\P%))))G1TFj[m$\"3Q$G9) H`AEDF-7$$\"3W**\\i+gQ9TFj[m$\"3KQ6$G8syG$F-7$$\"3p*\\7G6$[ATFj[m$\"3F Y$z5b2@\"RF-7$$\"3$******\\A!eITFj[m$\"3[4S$32)[uVF-7$$\"3uu$4@TP_0>%F-7$$\"3A]ilZLguTFj[m$\"3%)[31Pw&)RPF-7$$\" 3-DcwM0%>=%Fj[m$\"3]DnlBpn'=$F-7$$\"3%)**\\(=sx#*=%Fj[m$\"3E[FW]u2aDF- 7$$\"3kuV)*3\\h'>%Fj[m$\"3=/a\"*R0On=F-7$$\"3X\\P4'4_R?%Fj[m$\"3[N/A.I &H:\"F-7$$\"3(\\7.KG*G6UFj[m$\"3?3.4_:wsVF07$$\"3y*\\7.ZE'=UFj[m$!3#z: \\ykgNa#F07$$\"3fu=UdO'fA%Fj[m$!3[#y2O:/w)*)F07$$\"36]7`W3ILUFj[m$!3'4 !)G$)*Ryv9F-7$$\"3#\\iS;.Q1C%Fj[m$!3qM'**H_T!p>F-7$$\"3u***\\(=_(zC%Fj [m$!3#)>+do'fkO#F-7$$\"3CDcE-n*fD%Fj[m$!3M$p?crUEo#F-7$$\"3/]7y&==SE%F j[m$!3w&e/f)eEsGF-7$$\"33i!Rv#*G!oUFj[m$!3$)Q\"4#>8U?HF-7$$\"3#[(oHp'R ?F%Fj[m$!3rRuS@9oQHF-7$$\"3e(oa5T]gF%Fj[m$!3B!R!>F.SGHF-7$$\"3M+D\"G:h +G%Fj[m$!3#)*))R5:b7*GF-7$$\"3#*\\P%)>T5'H%Fj[m$!3?_TeAv(o^#F-7$$\"3@+ ](o3Z@J%Fj[m$!32C3=kz00>F-7$$\"3^]i!R0!>GVFj[m$!3Pt)\\#zl`J7F-7$$\"35+ v$4-LUM%Fj[m$!3-yQ7['F07$$\"3S](oz)fFgVFj[m$!3y0IA3>X'[#F07$$\"3q+ ++b*=jP%Fj[m$!3I<$*QC!es4&Fbz7$$\"3Av$feamQQ%Fj[m$!3O$)Q3oyzh9Fbz7$$\" 3U](=n899R%Fj[m$!3es$[Aj:#y:F\\cm7$$\"3kD\"yvsh*)R%Fj[m$\"3%z_3AU[6'>! #C7$$\"3;+vV=$4lS%Fj[m$\"34!\\(HzcGDGF\\cm7$$\"3g]i:+Xg@WFj[m$\"3eO(yo Ns>;'Fbz7$$\"3L+](=o*pOWFj[m$\"3!z_K\\_.De#F07$$\"31]Pfj[z^WFj[m$\"3= \\\"*=2J$yM'F07$$\"3]+DJX+*oY%Fj[m$\"3ob(p2TMg<\"F-7$$\"3A]7.F_)>[%Fj[ m$\"3#RA45UPr!=F-7$$\"3%****\\(3/3(\\%Fj[m$\"3.&f)>@*QbS#F-7$$\"3j**\\ 7B#fG^%Fj[m$\"3XWCTb.tPGF-7$$\"3.++]P!Q'GXFj[m$\"31,q*\\5t'=HF-7$$\"3U +](=&oTWXFj[m$\"3)p!=(*zvbUDF-7$$\"35++Dmc>gXFj[m$\"35LIm$GPeo\"F-7$$ \"3%**\\PM2&3oXFj[m$\"3A,BtVCm'4\"F-7$$\"3!)**\\i![ufd%Fj[m$\"3B&RZBD@ !HUF07$$\"3k*\\7y)Q'Qe%Fj[m$!3]m/&fKF-7$$\"3E++vB4JBYFj[m$!3sZ4Ye- 0QQF-7$$\"3`]i!*)fP(QYFj[m$!3?02q7CO\"f%F-7$$\"35+D1uU;aYFj[m$!3/D#4@w ))[s%F-7$$\"3AD1k6w(=m%Fj[m$!3Qy?;lPyPXF-7$$\"3k\\(=#\\4fpYFj[m$!3!R)z (y_zQ=%F-7$$\"3yuoz'G/tn%Fj[m$!3aoXGMi!Rn$F-7$$\"3#***\\PCw,&o%Fj[m$!3 F.8(\\FB`-$F-7$$\"3/DJ&>'4t#p%Fj[m$!3Kz;E#F-7$$\"3=]7`*HW/q%Fj[m$ !3mpR:XXK69F-7$$\"3hu$4rjd\"3ZFj[m$!3_SpIC)px1&F07$$\"3u*\\(ou4(er%Fj[ m$\"3sf?ovu!H<%F07$$\"3)[ilAJ%eBZFj[m$\"3-pa>xOGD8F-7$$\"3I\\P%)\\wHJZ Fj[m$\"3[\">]R?G]%>>k;OF-7$$\"3k\\(=UiuGw%Fj[m$\"36-cEcm(o`%F-7$$\" 3q*\\PM#\\-zZFj[m$\"3\\'fzesSZt%F-7$$\"33vo/t+5(y%Fj[m$\"31AV]oc=eXF-7 $$\"3w\\ilA_<&z%Fj[m$\"3;;mUGqE7UF-7$$\"3WCcEs.D.[Fj[m$\"3@t'>qF@rr$F- 7$$\"3#)**\\(=_D8\"[Fj[m$\"3+h3G273*4$F-7$$\"3@vV[r1S>[Fj[m$\"3OB-bL`? *Q#F-7$$\"3*)\\P4@eZF[Fj[m$\"38\"*)>q5&f@;F-7$$\"3dCJqq4bN[Fj[m$\"3*\\ =***>\"[fJ)F07$$\"3&**\\7.7EO%[Fj[m$\"3Qxt0o(GuR&Fbz7$$\"3Mv=#*p7q^[Fj [m$!3YjOX*=e!)y'F07$$\"3-]7`>kxf[Fj[m$!3l`33eE:Q8F-7$$\"3qC19p:&y'[Fj[ m$!3^'eF/RT1!>F-7$$\"33++v=n#f([Fj[m$!3Yr$fbp!**[BF-7$$\"35D\"yv//P)[F j[m$!39q5omf+jEF-7$$\"35]iSw8[\"*[Fj[m$!3yGJ$[.E!eGF-7$$\"3g7.#3/q`*[F j[m$!3(*R\"Q+B&[6HF-7$$\"35vVB0(e#**[Fj[m$!3Y*>:umyl$HF-7$$\"3gP%['pt9 .\\Fj[m$!3')ofj>*fW$HF-7$$\"36+D1Mg.2\\Fj[m$!3t'ebplvl!HF-7$$\"38](=

!*) !#A7$$\"3!3Dc,>J,/&Fj[m$\"3,F%4V1V!Q7Fbz7$$\"3k+v=KY3[]Fj[m$\"3'o_[$ft dX[Fbz7$$\"3[](=U2Qg0&Fj[m$\"3zOhor]K67F07$$\"3I++D;:*R1&Fj[m$\"3)z:p, 9&f$R#F07$$\"3)**\\7.S)*)z]Fj[m$\"33PMU\\bb\"G'F07$$\"3O+]P%G0e4&Fj[m$ \"3#z.IlV`**>\"F-7$$\"3u+vVo@r6^Fj[m$\"3/Mr?46_m=F-7$$\"3T++]_!>w7&Fj[ m$\"3Em)yY5K?[#F-7$$\"3wDc,Om]N^Fj[m$\"3[A,SsGP:FF-7$$\"3S]7`>URV^Fj[m $\"3;mP!*)HVM(GF-7$$\"3Qi!*G6!Qt9&Fj[m$\"36aO%=Py%=HF-7$$\"3/vo/.=G^^F j[m$\"3DC\"Gxum\"QHF-7$$\"3s(o/[fD_:&Fj[m$\"3#G2QpT73$HF-7$$\"3S+Dc'Qp \"f^Fj[m$\"3W'46T=5]*GF-7$$\"3S]Pf`X%\\<&Fj[m$\"3#)f>\\1t!4X#F-7$$\"3Q +]i?(>2>&Fj[m$\"3t%eS&H/#3`\"F-7$$\"3uD19/tg)>&Fj[m$\"3-\"4Tm%e\\g\"*F 07$$\"3P]il()[\\1_Fj[m$\"3#f^JNQ$pGAF07$$\"3,v=6T!Hn5,V* >`PF-7$$\"3]+]P\\L!=J&Fj[m$!3'\\8w@]-U;$F-7$$\"3;DJX>70>`Fj[m$!3[42!*) H2!pCF-7$$\"3`]7`*3*HE`Fj[m$!3ud=Z:YU!p\"F-7$$\"3\"fP4'fpaL`Fj[m$!3QNP !R&R-X&)F07$$\"3d+voH[zS`Fj[m$\"3!>UHF:Ay/\"Fbz7$$\"3CDcw*pU![`Fj[m$\" 3uVz23A-^()F07$$\"3g]P%)p0Hb`Fj[m$\"3b\"=.>^o*4aFj[m$\"39k(\\(HF!zR%F-7$$ \"3]](=#\\W$zU&Fj[m$\"3=)\\Wh=()o&RF-7$$\"3f++DE8COaFj[m$\"3'Q#GR&y!fs LF-7$$\"3n]7G.#[XW&Fj[m$\"3]lf\"3)GywEF-7$$\"3/+DJ!3bGX&Fj[m$\"3A&f*ph ]X0>F-7$$\"37]PMd>;haFj[m$\"3MZsT$=!z'4\"F-7$$\"3?+]PM)o%paFj[m$\"3I%* \\*=gT1*GF07$$\"3G]iS6dxxaFj[m$!3lW.J\"HqR\"[F07$$\"3O+vV)e#3'[&Fj[m$! 3Y@])=w(3#=\"F-7$$\"3W](oaY*Q%\\&Fj[m$!3?@Jab`\"fy\"F-7$$\"3`++]Ujp-bF j[m$!3k@W1oyUsAF-7$$\"31+v=#)eb6`M/zd*F07$$\"36+DJ?Jr1cFj[m$!37Z3.a`_(y%F07$$\"3l******fE d@cFj[m$!37cn(=p,fr\"F07$$\"3A+v=(RfKl&Fj[m$!3'*4=q0'=h2#Fj]l7$$\"33+] PMh%\\o&Fj[m$\"3=DVsG%fCG\"F07$$\"3O](oH]*y+dFj[m$\"3Gh%Q-[1<@%F07$$\" 3$**\\i:(Gj;dFj[m$\"3,>p9i\"\\'=\"*F07$$\"3A]i:SiZKdFj[m$\"3YMlQC:NY:F -7$$\"3]++v3'>$[dFj[m$\"3b&y.z7JV?#F-7$$\"3[+D1MP[jdFj[m$\"3h9E([D!p1F F-7$$\"3Y+]PfykydFj[m$\"3wuT`()*Ge$HF-7$$\"3U+vo%)>\"Qz&Fj[m$\"3=win%* RivFF-7$$\"3S+++5h(*3eFj[m$\"3))e!)\\')zF07$$\"3q]7GBk)o%eFj[ m$!3CAYJ\"=kjJ*F07$$\"3M+v$f[oW&eFj[m$!3s*3!o_!zHn\"F-7$$\"3)*\\Pf[00i eFj[m$!37:^dvCa#R#F-7$$\"3K++D6EjpeFj[m$!3Cx\"\\gud51$F-7$$\"3mvo/t0yx eFj[m$!3/HJ&3@z-p$F-7$$\"3G]P%[`Gf)eFj[m$!37h6l'=&z&>%F-7$$\"3!\\iSm\\ wS*eFj[m$!3GI-\"zc:.b%F-7$$\"3B+vVeWA-fFj[m$!3paFC#f=Ht%F-7$$\"3>]7.#Q ?&=fFj[m$!3/O[Qg#ffFj[m$!3+WJ,lq$)f6F-7$$\"3UiS \">!\\C1I1oF07$$\"33+D\"G:3u'fFj[m$!3a=%oEK;4%>F07$$\"3 wP4rL@[rfFj[m$\"3;s$>\"4YQXHF07$$\"3Uv$4Y6cb(fFj[m$\"3o&eIi\"e%**z(F07 $$\"3Q7y]&4I'zfFj[m$\"3mK6D].0d7F-7$$\"3/]iSwSq$)fFj[m$\"3d9wyPKe? " 0 "" {MPLTEXT 1 0 17 "# end of sectio n " }}}}{SECT 1 {PARA 3 "" 0 "" {TEXT -1 10 "Assignment" }}{PARA 0 "" 0 "" {TEXT 297 3 "(1)" }{TEXT -1 83 " Consider the Undamped Oscillator with Forcing defined by the initial value problem" }}{PARA 256 "" 0 " " {TEXT -1 0 "" }{TEXT 257 43 "u'' + 4u = F(t), u(0) = 0 and u'(0) = 0" }}{PARA 0 "" 0 "" {TEXT -1 41 "with forcing terms F(t) described below. " }}{PARA 0 "" 0 "" {TEXT -1 1 " " }{TEXT 256 4 "(a) " }{TEXT -1 77 "The forcing term is periodic with frequency equal to the natura l frequency. (" }{TEXT 294 9 "Resonance" }{TEXT -1 1 ")" }}{PARA 0 "" 0 "" {TEXT 258 0 "" }{TEXT 259 5 " (b) " }{TEXT -1 80 "The forcing ter m is periodic with frequency \"close to\" the natural frequency. (" } {TEXT 295 5 "Beats" }{TEXT -1 1 ")" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }} {PARA 0 "" 0 "" {TEXT -1 15 "For each case: " }}{PARA 0 "" 0 "" {TEXT -1 51 "(i) Define a forcing term matching the description." }}{PARA 0 "" 0 "" {TEXT -1 25 "(ii) Find the solution. " }}{PARA 0 "" 0 "" {TEXT -1 82 "(iii) Plot the solution over a long enough time period th at a pattern is apparent." }}{PARA 0 "" 0 "" {TEXT -1 84 "(iv) Describ e the behavior of the solution in terms of oscillations and amplitude s." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }} {PARA 0 "" 0 "" {TEXT 296 3 "(2)" }{TEXT -1 81 " Consider the Damped O scillator with Forcing defined by the initial value problem" }}{PARA 256 "" 0 "" {TEXT -1 0 "" }{TEXT 256 49 "u'' + 2u' + 5u = F(t), u(0 ) = 0 and u'(0) = 0" }}{PARA 0 "" 0 "" {TEXT -1 41 "with forcing term s F(t) described below. " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 " " 0 "" {TEXT 299 3 "(a)" }{TEXT -1 99 " Find a nonzero forcing term wh ich results in a solution that tends to zero as t goes to infinity. " }}{PARA 0 "" 0 "" {TEXT 300 3 "(b)" }{TEXT -1 94 " Find a forcing term which results in a solution that becomes unbounded as t goes to infin ity." }}{PARA 0 "" 0 "" {TEXT 301 3 "(c)" }{TEXT -1 133 " Find a forci ng term which results in an oscillatory solution that does not tend to zero or become unbounded at t goes to infinity. " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 15 "For each case: " }}{PARA 0 "" 0 "" {TEXT -1 87 "(i) Define a forcing term which results in a so lution with the requested properties. " }}{PARA 0 "" 0 "" {TEXT -1 80 "(ii) Plot the solution over a long enough period that the behavior is apparent. " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT 288 51 "- You may work in groups of no more than 2 people. " }} {PARA 0 "" 0 "" {TEXT 287 53 "- Five of the points will be based on pr esentation. " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }{TEXT 256 34 "- You sh ould hand in three pages: " }}{PARA 0 "" 0 "" {TEXT 257 28 " one pag e for problem 1(a)" }}{PARA 0 "" 0 "" {TEXT 258 3 " " }{TEXT 298 25 "one page for problem 1(b)" }}{PARA 0 "" 0 "" {TEXT 259 98 " one pag e for problem 2 (you should be able to fit all three functions and gra phs on one page). " }}{PARA 0 "" 0 "" {TEXT 289 162 "- You may write u p the lab in any text editor you want but make sure you can cut and pa ste Maple graphics into the document and describe mathematical equatio ns. " }}{PARA 0 "" 0 "" {TEXT 290 51 "- I suggest you hand in a print ed Maple document. " }}{PARA 0 "" 0 "" {TEXT 291 84 "- Only hand in t he requested information. I do not want to see the Maple commands. " }{TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT 292 80 "- To enter the text edi tor in Maple, click on the \"T\" button from the tool bar. " }}{PARA 0 "" 0 "" {TEXT 293 72 "- To get a Maple prompt, click on the \"[>\" b utton from the tool bar. " }{TEXT -1 0 "" }}}}{MARK "3" 0 } {VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }