{VERSION 3 0 "IBM INTEL NT" "3.0" } {USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 1 0 0 0 0 }{CSTYLE "2D Math" -1 2 "Times" 0 1 0 0 0 0 0 0 2 0 0 0 0 0 0 }{CSTYLE "2D Comment" 2 18 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 } {CSTYLE "2D Output" 2 20 "" 0 1 0 0 255 1 0 0 0 0 0 0 0 0 0 }{CSTYLE " " -1 256 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 257 "" 0 14 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 258 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 259 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE " " -1 260 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 261 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 262 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 263 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 264 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 265 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 266 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 267 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 268 "" 0 14 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 269 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 270 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 271 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 272 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 273 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 274 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 275 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 276 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 277 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 278 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 279 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 280 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 281 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 282 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 283 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 284 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 285 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 286 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 287 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 288 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 289 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 290 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 291 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 292 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 293 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 294 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 295 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 296 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 297 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 298 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 299 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 300 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 301 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 302 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 303 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 304 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 305 "" 1 18 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 306 "" 1 18 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 307 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE " " -1 308 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 309 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 310 "" 0 1 0 0 0 0 0 1 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 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 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 }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 }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 }3 1 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "" 0 257 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }} {SECT 0 {EXCHG {PARA 256 "" 0 "" {TEXT -1 60 "Lab 1: Solving Different ial Equations, Mechanical Vibrations" }}{PARA 257 "" 0 "" {TEXT -1 42 "Click on the + symbol to open that section" }}{PARA 0 "" 0 "" {TEXT 273 78 "Evaluate all input (in red) by placing the cursor anywhere in \+ red and enter. " }{TEXT -1 1 " " }}}{SECT 1 {PARA 3 "" 0 "" {TEXT -1 33 "1. Solving Differential Equations" }{TEXT 260 20 ": The dsolve com mand" }}{PARA 0 "" 0 "" {TEXT -1 71 "Differential equations and initia l value problems are solved using the " }{TEXT 272 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 condition(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 sections wit h the " }{TEXT 278 7 "restart" }{TEXT -1 55 " command to clear all pre vious 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 solution to t he first order differential equation: " }{TEXT 287 8 "y' = 3 y" }} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 46 "sol1 := dsolve(diff(y(x),x) \+ = 3 * y(x), y(x));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%sol1G/-%\"yG6 #%\"xG*&%$_C1G\"\"\"-%$expG6#,$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 ": Conside r the initial value problem " }{TEXT 279 8 "y' = 3 y" }{TEXT -1 5 " , \+ " }{TEXT 280 12 "y(0) = 1/2. " }}{PARA 0 "" 0 "" {TEXT -1 1 "(" } {TEXT 281 1 "a" }{TEXT -1 21 ") Find the solution. " }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 58 "sol2 := dsolve(\{diff(y(x),x) = 3*y(x), y(0 ) = 1/2\}, y(x));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%sol2G/-%\"yG6# %\"xG,$-%$expG6#,$F)\"\"$#\"\"\"\"\"#" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 94 "Notice: Maple automatically solves for the arbitrary cons tant based on the initial condition. " }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 1 "(" }{TEXT 282 1 "b" }{TEXT -1 41 ") Plot the solution from x = 0 to x = 1. " }}{PARA 0 "" 0 "" {TEXT -1 96 "First we want to acc ess the right hand side (rhs) of sol2 and make it a function of x call ed y2." }}{PARA 0 "" 0 "" {TEXT -1 22 "This is done with the " }{TEXT 283 7 "unapply" }{TEXT -1 5 " and " }{TEXT 284 3 "rhs" }{TEXT -1 24 " \+ commands as follows. " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 27 " y2 := unapply(rhs(sol2),x);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#y2GR 6#%\"xG6\"6$%)operatorG%&arrowGF(,$-%$expG6#,$9$\"\"$#\"\"\"\"\"#F(F(F (" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 49 "You can plot the solution fr om x = 0 to x = 1 by " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 84 "pl ot(y2(x),x=0..1,labels=[\"x\",\"y2(x)\"],title=\"The solution to y' = \+ 3y, y(0) = 1/2\");" }}{PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "6&-%'CURVESG6$7S7$\"\"!$\"1+++++++]!#;7$$\"1nmm;arz@!#<$ \"1h&3XH%)yL&F+7$$\"1LL$e9ui2%F/$\"1gX!p&))R]cF+7$$\"1nmm\"z_\"4iF/$\" 1z`$G5lP-'F+7$$\"1mmmT&phN)F/$\"1RB>Jp_CkF+7$$\"1LLe*=)H\\5F+$\"1;4C8U &)\\oF+7$$\"1nm\"z/3uC\"F+$\"1%QT)pBIpsF+7$$\"1++DJ$RDX\"F+$\"1sDH'H-2 t(F+7$$\"1nm\"zR'ok;F+$\"1TK6T/rQ#)F+7$$\"1++D1J:w=F+$\"1ulmc%4$y()F+7 $$\"1LLL3En$4#F+$\"1R:)HB\\-P*F+7$$\"1nm;/RE&G#F+$\"1*Hj'fbgC**F+7$$\" 1+++D.&4]#F+$\"1%HMq#=!)e5!#:7$$\"1+++vB_BO\"Fco7$$\"1LLLLY.KNF+$\"1\\SlCLiU9Fco7$$\"1++D \"o7Tv$F+$\"1bio*y4?a\"Fco7$$\"1LLL$Q*o]RF+$\"1&pc#f;oN;Fco7$$\"1++D\" =lj;%F+$\"1*oVFk8]u\"Fco7$$\"1++vV&R#*[l8z>Fco7$$\"1LLeR\"3Gy%F+$\"1c&fjo\"\\*4#Fco7$$\"1nm;/T1 &*\\F+$\"1,(=\"=)GvB#Fco7$$\"1mm\"zRQb@&F+$\"1t\"=N=G0R#Fco7$$\"1***\\ (=>Y2aF+$\"1FI&fN2A`#Fco7$$\"1mm;zXu9cF+$\"1bM&**>rYp#Fco7$$\"1+++]y)) GeF+$\"1@)>rfmM(GFco7$$\"1****\\i_QQgF+$\"1i:#*ew&)fIFco7$$\"1***\\7y% 3TiF+$\"1k&)yJ5q^KFco7$$\"1****\\P![hY'F+$\"1D![d " 0 "" {MPLTEXT 1 0 16 "subs(x=1,y2(x) );" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,$-%$expG6#\"\"$#\"\"\"\"\"#" }} }{EXCHG {PARA 0 "" 0 "" {TEXT -1 53 "and get a numerical value to 6 si gnificant digits by " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "eva lf(%,6);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"'G/5!\"%" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 80 "The % symbol references Maple's last outp ut. Not necessarily the line above it." }}{PARA 0 "" 0 "" {TEXT -1 1 "(" }{TEXT 286 2 "d)" }{TEXT -1 74 " Find the value of x for which y2 \+ is 100. Ie. Solve y2(x) = 100 for x. " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "solve(y2(x) = 100, x);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,$-%#lnG6#\"$+##\"\"\"\"\"$" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "evalf(%,6);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"'5m " 0 "" {MPLTEXT 1 0 43 "ode3 := diff(y(x),x$2)+diff(y(x),x) +y(x)=0;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%ode3G/,(-%%diffG6$-%\"y G6#%\"xG-%\"$G6$F-\"\"#\"\"\"-F(6$F*F-F2F*F2\"\"!" }}}{EXCHG {PARA 0 " " 0 "" {TEXT -1 28 "and then solve it using the " }{TEXT 274 6 "dsolve " }{TEXT -1 8 " command" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 " dsolve(ode3, y(x));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"yG6#%\"xG, &*(%$_C1G\"\"\"-%$expG6#,$F'#!\"\"\"\"#F+-%$cosG6#,$*&-%%sqrtG6#\"\"$ \"\"\"F'F+#F+F2F+F+*(%$_C2GF+F,F<-%$sinGF5F+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 solutio n to the initial value problem: " }{TEXT 289 37 "y'' + y' + y = 0, y(0 ) = 1, y'(0) = 0" }{TEXT -1 1 "." }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 38 "dsolve(\{ode3,y(0)=1,D(y)(0)=0\}, y(x));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"yG6#%\"xG,&*&-%$expG6#,$F'#!\"\"\"\"#\"\"\"-%$cosG 6#,$*&-%%sqrtG6#\"\"$\"\"\"F'F1#F1F0F1F1*(F7F;F*F;-%$sinGF4F1#F1F:" }} }{EXCHG {PARA 0 "" 0 "" {TEXT -1 82 "Again, Maple solves for the arbi trary constants based on the initial conditions. " }}{PARA 0 "" 0 "" {TEXT -1 55 "Define y4 as the right hand side of the previous output" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "y4 := unapply(rhs(%),x); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#y4GR6#%\"xG6\"6$%)operatorG%&ar rowGF(,&*&-%$expG6#,$9$#!\"\"\"\"#\"\"\"-%$cosG6#,$*&-%%sqrtG6#\"\"$\" \"\"F2F6#F6F5F6F6*(F " 0 "" {MPLTEXT 1 0 100 "plot(y4(x),x=0..10,labels=[\"x\",\"y4\"],ti tle=\"The solution to y'' + y' + y = 0, y(0) = 1, y'(0) = 0\");" }} {PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "6&-%'CURVESG6$7Y7$ \"\"!$\"\"\"F(7$$\"1LL$3FWYs#!#<$\"1#ygq'=K'***!#;7$$\"1mmmT&)G\\aF.$ \"1mzdzAU&)**F17$$\"1++]7G$R<)F.$\"1$3o>K.v'**F17$$\"1LLL3x&)*3\"F1$\" 1xE>(ymF%**F17$$\"1++]ilyM;F1$\"1?z:lekt)*F17$$\"1nmm;arz@F1$\"1#3@A\" Gmz(*F17$$\"1++D\"y%*z7$F1$\"1Mm-v[bh&*F17$$\"1LL$e9ui2%F1$\"1x^lF3@\" G*F17$$\"1nmm\"z_\"4iF1$\"1w.tAgRk%)F17$$\"1ommT&phN)F1$\"1Vlg^B!>X(F1 7$$\"1LLe*=)H\\5!#:$\"1:^i>\"*eKjF17$$\"1nm\"z/3uC\"F\\o$\"1XO.[-Da_F1 7$$\"1++DJ$RDX\"F\\o$\"1E%oQ'[!f9%F17$$\"1nm\"zR'ok;F\\o$\"1JQvaN%40$F 17$$\"1++D1J:w=F\\o$\"1V2Psq4X?F17$$\"1MLL3En$4#F\\o$\"1Qm8=`(\\7\"F17 $$\"1nm;/RE&G#F\\o$\"1;[/#oBtB%F.7$$\"1+++D.&4]#F\\o$!1+pX2R*>O#F.7$$ \"1+++vB_c%*fc3Q\"F17$$\"1MLL347TLF\\o$!1\"y0(yQ2d:F17$$\"1MLLLY .KNF\\o$!1b6n*z`Ei\"F17$$\"1++D\"o7Tv$F\\o$!1X[K[x$yh\"F17$$\"1LLL$Q*o ]RF\\o$!1y*e\\+hVb\"F17$$\"1,+D\"=lj;%F\\o$!1=Nz_&\\cV\"F17$$\"1++vV&R Y2aF\\o$!1FsuENugSF.7$$\"1nm;zXu9cF\\o$!1 V?,(y%4VDF.7$$\"1+++]y))GeF\\o$!1$o&G1Hbq6F.7$$\"1++]i_QQgF\\o$!1@tbjT -=P!#>7$$\"1,+D\"y%3TiF\\o$\"1/1W&Q7*o&)!#=7$$\"1++]P![hY'F\\o$\"1$z[z /)e?;F.7$$\"1LLL$Qx$omF\\o$\"1j:5,Y^7@F.7$$\"1+++v.I%)oF\\o$\"1mP?6Et_ CF.7$$\"1mm\"zpe*zqF\\o$\"1y[SH:![h#F.7$$\"1,++D\\'QH(F\\o$\"1deYx=.cE F.7$$\"1LLe9S8&\\(F\\o$\"1f::cee(e#F.7$$\"1,+D1#=bq(F\\o$\"1Xl)e)GcGCF .7$$\"1LLL3s?6zF\\o$\"1:\"fikE(3AF.7$$\"1++DJXaE\")F\\o$\"1wyY\\<]K>F. 7$$\"1ommm*RRL)F\\o$\"1_xwmw\"4k\"F.7$$\"1om;a<.Y&)F\\o$\"1FrgL]jL8F.7 $$\"1NLe9tOc()F\\o$\"1?G'fquO.\"F.7$$\"1,++]Qk\\*)F\\o$\"1&eVk)GlS'f@pas(Fgv7$$\"1,+v=5s#y*F\\o$!1E " 0 "" {MPLTEXT 1 0 58 " ode3 :=5*diff(y(x),x,x)+10*sqrt(2)*diff(y(x),x)+10*y(x)=0;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%ode3G/,(-%%diffG6$-%\"yG6#%\"xG-%\"$G6$F- \"\"#\"\"&*&-%%sqrtG6#F1\"\"\"-F(6$F*F-\"\"\"\"#5F*F;\"\"!" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "dsolve(\{ode3\},y(x));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"yG6#%\"xG*&%$_C1G\"\"\"-%$expG6#,$*&-%%sqr tG6#\"\"#\"\"\"F'F*!\"\"F*" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 76 "HEY ! That's not right! There should be two linearly independent solutions . " }}{PARA 0 "" 0 "" {TEXT -1 50 "Check the roots of the characteris tic equation by " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 37 "solve(5 *r^2 + 10*sqrt(2)*r + 10=0,r);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6$,$*$ -%%sqrtG6#\"\"#\"\"\"!\"\"F#" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 31 "T here is a double root: r = " }{XPPEDIT 18 0 "-sqrt(2);" "6#,$-%%sqr tG6#\"\"#!\"\"" }{TEXT -1 52 " so a second (linearly independent) sol ution is x " }{XPPEDIT 18 0 "exp(-sqrt(2)*x);" "6#-%$expG6#,$*&-%%sqr tG6#\"\"#\"\"\"%\"xGF,!\"\"" }{TEXT -1 3 " . " }}{PARA 0 "" 0 "" {TEXT -1 42 "Maple did not find this second solution. " }}}{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 motion of a mass on \+ a spring" }}{PARA 256 "" 0 "" {TEXT -1 8 "m u'' + " }{XPPEDIT 18 0 "ga mma;" "6#%&gammaG" }{TEXT -1 16 " u' + k u = F(t)" }}{PARA 0 "" 0 "" {TEXT -1 74 "Here, u is the diplacement from equilibrium, m is the mas s of the object, " }{XPPEDIT 18 0 "gamma;" "6#%&gammaG" }{TEXT -1 136 " is the damping constant, k is the spring constant, and F(t) is some \+ external forcing term. All of the constants are positive: a larger " } {XPPEDIT 18 0 "gamma;" "6#%&gammaG" }{TEXT -1 66 " implies greater dam ping 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 291 6 "damped" }{TEXT -1 7 ", i f " }{XPPEDIT 18 0 "gamma" "6#%&gammaG" }{TEXT -1 19 " = 0 the motion is " }{TEXT 290 9 "undamped." }}{PARA 0 "" 0 "" {TEXT -1 1 " " }} {PARA 0 "" 0 "" {TEXT -1 90 "We can define the left side (lhs) of this ordinary differential equation (ode) as follows " }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 55 "lhsode := m*diff(u(t),t$2) + gam*diff(u(t),t) + k*u(t);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%'lhsodeG,(*&%\"mG\"\" \"-%%diffG6$-%\"uG6#%\"tG-%\"$G6$F/\"\"#F(F(*&%$gamGF(-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 22 "dsolve(lhsode=0,u(t));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"uG6#%\"tG,&*&%$_C1G\"\"\"-%$expG6#,$*&*&,& %$gamGF+*$-%%sqrtG6#,&*$)F3\"\"#\"\"\"F+*&%\"mGF+%\"kGF+!\"%F!\"\"#!\"\"F;F+F+*&%$_C2GF+-F-6#,$*&*&,&F3FCF4F+F+F'FFA#F+F;F +F+" }}}{PARA 0 "" 0 "" {TEXT -1 205 "This is a little messy but you c an see how the parameters m, gam, and k effect our solution. If the t erm under the radical sign is negative the solution will contain sine \+ and cosine terms of the form cos(" }{XPPEDIT 18 0 "omega;" "6#%&omegaG " }{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.%\"mGF. F2" }{TEXT -1 18 ". In this case, " }{XPPEDIT 18 0 "omega;" "6#%&ome gaG" }{TEXT -1 23 " (omega) is called the " }{TEXT 293 9 "frequency" } {TEXT -1 31 " of the sine and cosine terms. " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }{TEXT 292 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 := 2:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 7 "k := 5:" }}} {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 25 "dsolve(lhsode=F(t),u(t));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"uG6#%\"tG,&*(%$_C1G\"\"\"-%$expG6#,$F'!\"\"F+-%$si nG6#,$F'\"\"#F+F+*(%$_C2GF+F,\"\"\"-%$cosGF3F+F+" }}}{EXCHG {PARA 0 " " 0 "" {TEXT -1 12 "Notice, the " }{TEXT 256 9 "frequency" }{TEXT -1 24 " of the trig terms is " }{XPPEDIT 18 0 "sqrt(abs(2^2-4(1)(5)))/2 (1);" "6#*&-%%sqrtG6#-%$absG6#,&*$\"\"#\"\"#\"\"\"--\"\"%6#\"\"\"6#\" \"&!\"\"F.-\"\"#6#\"\"\"F6" }{TEXT -1 80 " = 2. However, the solutio n is not periodic because of the exponential term. " }}{PARA 0 "" 0 " " {TEXT -1 44 "In this case the solution is said to have a " }{TEXT 294 15 "quasi-frequency" }{TEXT -1 8 " of 2. " }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 191 "Suppose that we displace the mass 1 unit beyond eq uilibrium and let it go. The initial conditions 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 "uo := 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 solution to the initi al value problem by " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 62 "s ol := dsolve(\{lhsode = 0,u(0) = uo, D(u)(0) = uprimeo\},u(t));" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%$solG/-%\"uG6#%\"tG,&*&-%$expG6#,$F) !\"\"\"\"\"-%$sinG6#,$F)\"\"#F1#F1F6*&F,\"\"\"-%$cosGF4F1F1" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "u1 := unapply(rhs(sol),t);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#u1GR6#%\"tG6\"6$%)operatorG%&arro wGF(,&*&-%$expG6#,$9$!\"\"\"\"\"-%$sinG6#,$F2\"\"#F4#F4F9*&F.\"\"\"-%$ cosGF7F4F4F(F(F(" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 100 "plot(u 1(t),t=0..10,labels=[\"t\",\"u(t)\"],title=\"the solution to u'' + 2u' + 5u = 0, u(0)=1, u'(0)=0\");" }}{PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "6&-%'CURVESG6$7co7$\"\"!$\"\"\"F(7$$\"1LLL3x&)*3\"!#; $\"1`]R,s%[s*F.7$$\"1nmm;arz@F.$\"1VsJ![,s)*)F.7$$\"1++D\"y%*z7$F.$\"1 dHe6(e-2)F.7$$\"1LL$e9ui2%F.$\"1B>(z]RC)pF.7$$\"1n;H2Q\\4YF.$\"1>)=&39 xBjF.7$$\"1++voMrU^F.$\"1'[s?Kflk&F.7$$\"1L$3-8Lfn&F.$\"1?!H0-+='\\F.7 $$\"1nmm\"z_\"4iF.$\"1=Eg`&z%zUF.7$$\"1nm;zp!fu'F.$\"19obK(\\Tg$F.7$$ \"1nmmm6m#G(F.$\"1;s6F9R[HF.7$$\"1om;a`T>yF.$\"1$3))3re\">BF.7$$\"1omm T&phN)F.$\"1zgsFfLA*Fgo7$$\"1nm\"z/3uC\" F[p$!1K+wsa#oU\"F.7$$\"1LLe*ot*\\8F[p$!1Y$*oQdn*y\"F.7$$\"1++DJ$RDX\"F [p$!1VF)y)p\\+?F.7$$\"1LLe9x0z9F[p$!1ESlr1UK?F.7$$\"1nm\"z4wb]\"F[p$!1 qr+Sttb?F.7$$\"1++D\"[%4K:F[p$!1aENi<\"32#F.7$$\"1LLekGhe:F[p$!1m$fL3= !y?F.7$$\"1nm\"zCJ^e\"F[p$!1j3zd!Qx2#F.7$$\"1++DJ'\\;h\"F[p$!1KxY]$f.2 #F.7$$\"1LLe9!o\"Q;F[p$!1D?\"Q!RFc?F.7$$\"1nm\"zR'ok;F[p$!1r:2vj(e.#F. 7$$\"1LL3_(>/x\"F[p$!1`CPh`j**=F.7$$\"1++D1J:w=F[p$!1qka\\y4%p\"F.7$$ \"1MLL3En$4#F[p$!11S3t>)3:\"F.7$$\"1nm;/RE&G#F[p$!1yt!H*p\"[Z'Fgo7$$\" 1+++D.&4]#F[p$!1L#fCAP&)e\"Fgo7$$\"1+++vB_(R*G-\"*=Fgo7$$\" 1+++Dg(=#GF[p$\"1[5?az)))*HFgo7$$\"1+++v'Hi#HF[p$\"12n&oMoYC%Fgo7$$\"1MLL347TLF[p$\"1J9^L7G\\RFgo7$$\"1 MLLLY.KNF[p$\"19]S]Rd1JFgo7$$\"1++D\"o7Tv$F[p$\"1Z$>kw@b*=Fgo7$$\"1LLL $Q*o]RF[p$\"1;'Q]ak#)p)!#=7$$\"1,+D\"=lj;%F[p$!11e:bHH=E!#>7$$\"1++vV& Rx))Fex7$$\"1nm;/T1&*\\F[p$!10?)H;59`(Fex7$$\"1nm\"zRQb@& F[p$!1>qDO`B+_Fex7$$\"1++v=>Y2aF[p$!1lGG2^[5IFex7$$\"1nm;zXu9cF[p$!15k 06]!oF*F[y7$$\"1+++]y))GeF[p$\"1]-y3Pn'['F[y7$$\"1++]i_QQgF[p$\"1)QR. \"fFW:Fex7$$\"1,+D\"y%3TiF[p$\"1G84Np$*e=Fex7$$\"1++]P![hY'F[p$\"1*>Px 7w/t\"Fex7$$\"1LLL$Qx$omF[p$\"1T+Enl6a8Fex7$$\"1+++v.I%)oF[p$\"1=\"*e0 4?j%)F[y7$$\"1mm\"zpe*zqF[p$\"1X\")Hbk`J\"F[y7$$\" 1NLe9tOc()F[p$!1vt[O!>x+%!#?7$$\"1,++]Qk\\*)F[p$\"16pr')[/nAF[`l7$$\"1 NL$3dg6<*F[p$\"1#\\NF1;Vc'F[`l7$$\"1ommmxGp$*F[p$\"1$GF)Q+b0!)F[`l7$$ \"1++D\"oK0e*F[p$\"1+\"R;vu>j(F[`l7$$\"1,+v=5s#y*F[p$\"1>7<*\\J\"3hF[` l7$$\"#5F($\"1@sO2A2DRF[`l-%'COLOURG6&%$RGBG$Fgal!\"\"F(F(-%+AXESLABEL SG6$Q\"t6\"Q%u(t)Fdbl-%&TITLEG6#QTthe~solution~to~u''~+~2u'~+~5u~=~0,~ u(0)=1,~u'(0)=0Fdbl-%%VIEWG6$;F(Ffal%(DEFAULTG" 1 2 0 1 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 }}}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 38 "Voila! A damped oscillating solution. " }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }{TEXT 277 33 "Undamped Oscillator with Forcing:" } {TEXT -1 20 " m u'' + k u = F(t)" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 43 "restart: #clears all variable definitions. " }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 36 "lhsode := m*diff(u(t),t$2) + k*u(t);" }}{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 general solution of the homogeneous equation by" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "dsolve(lhsode=0,u(t));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-% \"uG6#%\"tG,&*&%$_C1G\"\"\"-%$sinG6#*&-%%sqrtG6#*&%\"kG\"\"\"%\"mG!\" \"F5F'F+F+F+*&%$_C2GF+-%$cosGF.F+F+" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 85 "In this case, the solution of the homogeneous equation is perio dic with frequency is " }{XPPEDIT 18 0 "sqrt(k/m);" "6#-%%sqrtG6#*&%\" kG\"\"\"%\"mG!\"\"" }{TEXT -1 1 "." }}{PARA 0 "" 0 "" {TEXT -1 93 "The frequency of the solution to the homogeneous solution in the undamped case is called the " }{TEXT 295 19 "natural frequency. " }}{PARA 0 " " 0 "" {TEXT -1 46 "Suppose k = 4 and m = 1. Let's see what we get" }} }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "m := 1: " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 7 "k := 4:" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 53 "The solution to the homogeneous equation is found by " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "dsolve(lhsode=0,u(t));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"uG6#%\"tG,&*&%$_C1G\"\"\"-%$sinG6#,$F'\"\" #F+F+*&%$_C2GF+-%$cosGF.F+F+" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 12 "Notice: The " }{TEXT 271 17 "natural fr equency" }{TEXT -1 4 " is " }{XPPEDIT 18 0 "sqrt(abs(0^2-4(1)(4)))/2(1 );" "6#*&-%%sqrtG6#-%$absG6#,&*$\"\"!\"\"#\"\"\"--\"\"%6#\"\"\"6#\"\"% !\"\"F.-\"\"#6#\"\"\"F6" }{TEXT -1 6 " = 2. " }}{PARA 0 "" 0 "" {TEXT -1 42 "Now assign the function F(t) = sin(3t) by " }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 20 "F := t -> sin(3*t): " }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 138 "Notice: The frequency of the forcing term is 3 and \+ the natural frequency is 2. Suppose we start the spring from the resti ng 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 59 "sol := dsolve(\{lhsode=F(t),u(0)=uo , D(u)(0)=uprimeo\},u(t));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$solG/ -%\"uG6#%\"tG,(*&,&-%$cosG6#,$F)\"\"&#!\"\"\"#?-F.F(#F3\"\"%\"\"\"-%$s inG6#,$F)\"\"#F8F8*&,&-F:F(F6-F:F/#F8F4F8-F.F;F8F8F9#\"\"$\"#5" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "u2 := unapply( rhs(sol), t); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#u2GR6#%\"tG6\"6$%)operatorG%&ar rowGF(,(*&,&-%$cosG6#,$9$\"\"&#!\"\"\"#?-F06#F3#F6\"\"%\"\"\"-%$sinG6# ,$F3\"\"#FF9F:-F>F1#F " 0 "" {MPLTEXT 1 0 20 "plot(u2(t),t=0..60);" }} {PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "6%-%'CURVESG6$7[dm 7$\"\"!F(7$$\"1+++DJdpK!#;$\"1lHUP!*GH;!#<7$$\"1+++]i9RlF,$\"1t,PZ?xZ5 F,7$$\"1++]7G$R<)F,$\"1btIh:p@(3)*F,$\"1\\wc;F/7$$\"1++voozw=FB$!1P\">CyFB$!1 #fK^xgL?\"F,7$$\"1+vVt(Q!>?FB$!1I:,&f9^*=F,7$$\"1]7yD(f,4#FB$!1vORZa') fDF,7$$\"1+]7y1Gh@FB$!1$3KGf\\N<$F,7$$\"1](o/j,CB#FB$!1r4L%y0Jr$F,7$$ \"1+D\"GeANI#FB$!1f^MuZUdTF,7$$\"1]i:NNkuBFB$!1mh:s(Q#)[%F,7$$\"1++]([ kdW#FB$!1!*Qc%)=\"4p%F,7$$\"1v$4@'fv&[#FB$!1D`%=t;Xu%F,7$$\"1](=nVZd_# FB$!1InR`-1`ZF,7$$\"1D\"G8\"*Qdc#FB$!1/uRR`&er%F,7$$\"1+v$fQIdg#FB$!1S KX&=GFj%F,7$$\"1]i:NLr&o#FB$!1&4W1191L%F,7$$\"1+]P%G'plFFB$!1O\"4Jz&pb QF,7$$\"1]PfL#zc%GFB$!1n9Q#=H\\A$F,7$$\"1+D\"G=ic#HFB$!1uA'eDQDY#F,7$$ \"1]7.K^k0IFB$!1wd:y:$*)f\"F,7$$\"1++D\"3Gc3$FB$!1Pgp7&))Hp'F/7$$\"1v$ feb>c7$FB$!1hl[z>C;>F/7$$\"1\\(o/.6c;$FB$\"1p_c%F,7$$\" 1++](=Wv/%FB$\"1(y>XXkbt$F,7$$\"1+]i:t0GTFB$\"1M<]\"QpI7$F,7$$\"1++vV/ d3UFB$\"1:&4yS1\"=CF,7$$\"1+](=d$3*G%FB$\"1>=*H5'Ra;F,7$$\"1,+++nfpVFB $\"1jEPzRdp')F/7$$\"1,]7G)4,X%FB$\"1]N?K3cJ!*!#=7$$\"1,+DcHiIXFB$!1,bz zY7KkF/7$$\"1,]P%3O6h%FB$!16=p2V108F,7$$\"1,+]7#\\;p%FB$!1\"pO/3F<(=F, 7$$\"1++voan_[FB$!1#)\\>\"4&>cEF,7$$\"1,++DDw#F,7$$\"1+]7yNAM`FB$!1S=-K#))*fAF,7$$\"1,vo/X[%\\ &FB$!1WF#*)3-)*f\"F,7$$\"1++DJauacFB$!1W(34Jl]^*F/7$$\"1,]P%Gn_(fFB$!1 ;toO%H=P\"F/7$$\"1++]P\"*y&H'FB$\"1(pJN1(*4+\"!#@7$$\"1+++DR&Hf'FB$\"1 x8yj'ecR\"F/7$$\"1++]7(=,*oFB$\"1\"3]oo#*4t)F/7$$\"1++D16qQqFB$\"1G.2R ))\\f9F,7$$\"1*******\\$G(=(FB$\"1]m$Qw#R$3#F,7$$\"1***\\P*e'eL(FB$\"1 s#H(\\zO9EF,7$$\"1****\\(G[W[(FB$\"1!)**HZ_i7HF,7$$\"1+]7./\"H_(FB$\"1 ?XiF,7$$ \"1+++vO*f%zFB$\"1uw=2^Kv8F,7$$\"1***\\i!z\"H-)FB$\"1*yO3PQ:a(F/7$$\"1 ++]P@%)*4)FB$\"1=MYpOKwjFaz7$$\"1)**\\(ojww\")FB$!1#4#*R&\\Z;nF/7$$\"1 *******f!p`#)FB$!1C9^`&3XU\"F,7$$\"1***\\7$[hI$)FB$!1yVIm92l@F,7$$\"1+ +]i!RvS)FB$!1y/9hhuiGF,7$$\"1++v$HjW[)FB$!18SFMz^([$F,7$$\"1,++DvQh&)F B$!1UwqRSP6SF,7$$\"1,+Dc$ oYF,7$$\"1,+]P9!R&*)FB$!1es4Tx!3S%F,7$$\"1++](ecM.*FB$!1>(y26Z+'RF,7$$ \"1****\\P<,8\"*FB$!1@LUIBFhLF,7$$\"1++]()oc#>*FB$!1KA)pl'*pi#F,7$$\"1 ,+]P?7s#*FB$!1,P\\cU-'y\"F,7$$\"1++](=x;N*FB$!1c\\h*)>T@()F/7$$\"1++]i ZX\"R*FB$!16#e/OvR*RF/7$$\"1****\\PBBJ%*FB$\"1))zgMD!\\u(Faz7$$\"1++]7 *45Z*FB$\"102N%Q)*\\`&F/7$$\"1++]([(y5&*FB$\"1]D/Sc'Q-\"F,7$$\"1,+]PEM !f*FB$\"1Du;QqPG>F,7$$\"1++](y(*)p'*FB$\"12M'yb>Uv#F,7$$\"1****\\PHX\\ (*FB$\"1<3o#*\\JoMF,7$$\"1++](33!H)*FB$\"1cZmJ()zUSF,7$$\"1,+]PKc3**FB $\"1(R6em!GcWF,7$$\"1++](Q=\"))**FB$\"1V_cJax%p%F,7$$\"1]7`Rox-5!#9$\" 1MX2#>4hu%F,7$$\"1+DJS=u15Fi[m$\"1&[[)*HhBv%F,7$$\"1]P4Toq55Fi[m$\"1*z kyjRSr%F,7$$\"1+](=%=n95Fi[m$\"17=t]+9KYF,7$$\"1+vVV=gA5Fi[m$\"1DL)y'H #RM%F,7$$\"1+++X=`I5Fi[m$\"1(>2B\"=l/RF,7$$\"1+DcY=YQ5Fi[m$\"1)=9!fnUP LF,7$$\"1+]7[=RY5Fi[m$\"1qs$)HkCqEF,7$$\"1+vo\\=Ka5Fi[m$\"1$>E0'*HX$>F ,7$$\"1++D^=Di5Fi[m$\"1k()Hj=]j6F,7$$\"1+D\"G&==q5Fi[m$\"1%y[J#yE0RF/7 $$\"1+]Pa=6y5Fi[m$!1\"=()*GBvX;uO5F,7$$ \"1++]d=(R4\"Fi[m$!1y1>K#\\vj\"F,7$$\"1+]ig=$)46Fi[m$!1$oTP*)yy^#F,7$$ \"1++vj=pD6Fi[m$!1<0ERxx8HF,7$$\"1]iSEe+U6Fi[m$!1#o&[h4wPGF,7$$\"1+D1* y>$e6Fi[m$!1LK*>3klQ#F,7$$\"1](=\"Fi[m$!1VW<)>^#f5F,7$$\"1+voRcdB7Fi[m$!1RZ>FK*=o\"F/7$$\"1+++lN?c7 Fi[m$!1h#=5ckI2%!#B7$$\"1]PfQ#)Qj7Fi[m$\"1Zsw7$$\"1+v=7Hdq7F i[m$\"1gMYPXAO8Faz7$$\"1]7y&edxF\"Fi[m$\"1'\\VDbubd%Faz7$$\"1+]PfA%\\G \"Fi[m$\"1Ndv'*\\(f2\"F/7$$\"1+Dc1;J*H\"Fi[m$\"1&z%y#z?dW$F/7$$\"1++v` 4o88Fi[m$\"1<]%R8FrY(F/7$$\"1+v$4I]!G8Fi[m$\"1].f]m](G\"F,7$$\"1+]7['> CM\"Fi[m$\"12`#*f?\\#*=F,7$$\"1+DJ&**)yc8Fi[m$\"1JGXsP(3X#F,7$$\"1++]U $e6P\"Fi[m$\"1*Qc\"yP!e$GF,7$$\"1D1*GY-_P\"Fi[m$\"1O*Q;^ug*GF,7$$\"1]7 G$eY#z8Fi[m$\"1[mI!e@3$HF,7$$\"1v=n.2H$Q\"Fi[m$\"1\"y$*fOe!QHF,7$$\"1+ D1C[L(Q\"Fi[m$\"1P4-U35;HF,7$$\"1]P%[1BaR\"Fi[m$\"1V>X'y&ezFF,7$$\"1+] i08^.9Fi[m$\"1(\\Py`v]^#F,7$$\"1]iSY&*f69Fi[m$\"1E\\0!=eH7#F,7$$\"1+v= (y(o>9Fi[m$\"10LDQDe5;F,7$$\"1](oz-wxU\"Fi[m$\"1Uu'))*fNA**F/7$$\"1++v oU'eV\"Fi[m$\"1'o,BXir)GF/7$$\"1]7`4D&RW\"Fi[m$!1)z[*RY![t%F/7$$\"1+DJ ]2/_9Fi[m$!1oxXl6Lj7F,7$$\"1]P4\"**G,Y\"Fi[m$!1JoQfeuY?F,7$$\"1+](=B<# o9Fi[m$!1nf)4,Z#)y#F,7$$\"1]ilsaIw9Fi[m$!1h^EBgx_MF,7$$\"1+vV8PR%[\"Fi [m$!1Pf#eMGv+%F,7$$\"1](=U&>[#\\\"Fi[m$!1_Lx=arBWF,7$$\"1+++&>q0]\"Fi[ m$!1Obs*Ra\"yYF,7$$\"1+v$4#4j/:Fi[m$!15%o6-(\\RZF,7$$\"1+](ok\"p3:Fi[m $!1$[F+mFXv%F,7$$\"1+D\"GP_F^\"Fi[m$!1eAxcRXAZF,7$$\"1++v)48o^\"Fi[m$! 1Xe_%pFIk%F,7$$\"1+]i]X$\\_\"Fi[m$!1efLG)))QM%F,7$$\"1++]-g0L:Fi[m$!1T +3T([h'QF,7$$\"1+]PauEC?F,7$$\"1+]7e.Ud:Fi[m$!1Jz')>=@u:F ,7$$\"1+D1%3\"[h:Fi[m$!1G01[\\N26F,7$$\"1+++5=al:Fi[m$!1r$z'p&zlG'F/7$ $\"1+v$f`-'p:Fi[m$!1!fx'pyKK9F/7$$\"1+](=EjOd\"Fi[m$\"1^=:$ofsV$F/7$$ \"1+D\"y)Rsx:Fi[m$\"1MfO)p[+F)F/7$$\"1++v8Zy\"e\"Fi[m$\"1\\dn%pN9I\"F, 7$$\"1+]ilh!**e\"Fi[m$\"19eodbo.AF,7$$\"1++]_FkB2%p$F,7$$\"1++D@0F9;Fi[m$\"1.`.o![8A%F,7$ $\"1+]7t>RA;Fi[m$\"1$>!pq5KuXF,7$$\"1+++DM^I;Fi[m$\"1Y$*4dMSTZF,7$$\"1 ++]sk;Y;Fi[m$\"1W#zZstg`%F,7$$\"1+++?&>=m\"Fi[m$\"1&>I1DZeq$F,7$$\"1++ vVgkp;Fi[m$\"1iz47QF1JF,7$$\"1++]nDZx;Fi[m$\"1^?*z7'3?CF,7$$\"1++D\"4* H&o\"Fi[m$\"1m[dACHy;F,7$$\"1+++:c7$p\"Fi[m$\"11a;)4Z18*F/7$$\"1++vQ@& 4q\"Fi[m$\"1u5nl+#=c\"F/7$$\"1++]i'y(3F,7$$\"1+](=n&f7=Fi[m$!10Gi/ZWE8F,7$$\"1]PfQ,\"o#=Fi[m$!1)H< 7_G$RyF/7$$\"1+DJ0Y-T=Fi[m$!1Aq!**>xFi[m$\"1**Gf\">!R**HFaz7$$\"1]P4m1;?>Fi[m$\"1lXJ& >=6,#F/7$$\"1+](=/jq$>Fi[m$\"1')e.ikF0fF/7$$\"1]ilFi[m$\"1y\"Qs] E+>\"F,7$$\"1+vV$zn3(>Fi[m$\"1Ci`wG#z*=F,7$$\"1](=#p,x()>Fi[m$\"19\"Gh W$3RDF,7$$\"1+++XDn/?Fi[m$\"1+g[/(f'3HF,7$$\"1]PMx@D3?Fi[m$\"1#)=-P7/M HF,7$$\"1+vo4=$=,#Fi[m$\"1'\\))Q$HiPHF,7$$\"1]7.U9T:?Fi[m$\"1r`Hb9E=HF ,7$$\"1+]Pu5**=?Fi[m$\"1^[.U>-vGF,7$$\"1+D1R.:E?Fi[m$\"1:tSe'4Vr#F,7$$ \"1++v.'4L.#Fi[m$\"1*Rd&G:v_CF,7$$\"1+vVo)o//#Fi[m$\"1$*Gi$He<4#F,7$$ \"1+]7L\"Gw/#Fi[m$\"1`!))yGQqj\"F,7$$\"1+D\"yR(ya?Fi[m$\"1Tl.7yb)4\"F, 7$$\"1++]im%>1#Fi[m$\"1mFN`%oI!\\F/7$$\"1+v=Ff5p?Fi[m$!1W0!Q:P7q\"F/7$ $\"1+](=>li2#Fi[m$!1^vf!o$Q@')F/7$$\"1+DccWU$3#Fi[m$!1:6y1Q)Gc\"F,7$$ \"1++D@Pe!4#Fi[m$!1b#y,)f6[AF,7$$\"1+v$f)Hu(4#Fi[m$!12G868:$*GF,7$$\"1 +]i]A!\\5#Fi[m$!1k$e**phQZ$F,7$$\"1+DJ::17@Fi[m$!1%Q5'Q/nnRF,7$$\"1+++ !y?#>@Fi[m$!1O3^!e0XN%F,7$$\"1voz1([v7#Fi[m$!1h`g?\"fzk%F,7$$\"1]PfLm( e8#Fi[m$!1l770r/bZF,7$$\"1D1RgX?W@Fi[m$!1EHdD^KlYF,7$$\"1+v=([KD:#Fi[m $!1xt#poWyP%F,7$$\"1vV)RTg3;#Fi[m$!1Ki^cDR,RF,7$$\"1]7yS$)=p@Fi[m$!1(H T$Hg>aKF,7$$\"1D\"yvE;v<#Fi[m$!1#>dvY0JY#F,7$$\"1+]P%>We=#Fi[m$!1TDZth Hi:F,7$$\"1QMxd\"3+>#Fi[m$!1]C`?i,$3\"F,7$$\"1v=<@@<%>#Fi[m$!13W;h)*f: fF/7$$\"17.d%3O$)>#Fi[m$!1.6u88mW$*Faz7$$\"1](oz/+D?#Fi[m$\"1PXefw=dSF /7$$\"1)=n8,km?#Fi[m$\"1o$*>iU<.!*F/7$$\"1Dcwuz#3@#Fi[m$\"1egiR'*y%Q\" F,7$$\"1jS;Q>*\\@#Fi[m$\"1%3i%Hkq`=F,7$$\"1+Dc,f:>AFi[m$\"1m))fu%R=I#F ,7$$\"1v$f$GQ[FAFi[m$\"1UqOTHC;JF,7$$\"1]i:b<\"eB#Fi[m$\"1^%48]LEz$F,7 $$\"1DJ&>oRTC#Fi[m$\"1R^%RoyEI%F,7$$\"1++v3wY_AFi[m$\"1lu_#4hki%F,7$$ \"1jS\">U`hD#Fi[m$\"1A!z!er>2ZF,7$$\"1D\"y]BR)fAFi[m$\"1t1UbO$*[ZF,7$$ \"1)=U#[]_jAFi[m$\"1SIlN`t^ZF,7$$\"1]iSh3@nAFi[m$\"1p&)4H*Qgr%F,7$$\"1 vVt([#euAFi[m$\"1#pF$[uqKXF,7$$\"1+D19T&>G#Fi[m$\"1X3fWsZ4UF,7$$\"1D1R SdK*G#Fi[m$\"1\\tElc%>w$F,7$$\"1](=nO(p'H#Fi[m$\"1%GI&RG,5KF,7$$\"1vo/ $**oSI#Fi[m$\"1pn5rN&pd#F,7$$\"1+]P>1W6BFi[m$\"1YAQq>O))=F,7$$\"1DJqXA \")=BFi[m$\"1I%*oJ&=5<\"F,7$$\"1]7.sQ=EBFi[m$\"1h8!yqBt^%F/7$$\"1v$f$) \\bNL#Fi[m$!1RQmu;)yV#F/7$$\"1+voCr#4M#Fi[m$!157$p[`)>*)F/7$$\"1Dc,^() H[BFi[m$!1K`C\"fvBZ\"F,7$$\"1]PMx.nbBFi[m$!1vt?&e1$o>F,7$$\"1v=n.?/jBF i[m$!1S_!3L;vO#F,7$$\"1+++IOTqBFi[m$!1\"HrbLsCm#F,7$$\"1v=Us9]yBFi[m$! 1j#)\\HM3jGF,7$$\"1]P%[J*e'Q#Fi[m$!1W6XIF+QHF,7$$\"1DcEdrn%R#Fi[m$!1Ck C/tW'*GF,7$$\"1+vo**\\w-CFi[m$!1f!)*f&\\+`FF,7$$\"1]7`%oS*=CFi[m$!10fr '=L%QAF,7$$\"1+]Ppj6NCFi[m$!1kP9^Dep:F,7$$\"1](=U0#H^CFi[m$!1[9'ebKv>* F/7$$\"1+D1RxYnCFi[m$!155*[#Hy#=%F/7$$\"1]i!RUVO[#Fi[m$!1:nwQE4G7F/7$$ \"1++v3\">)*\\#Fi[m$!1U*>'f:m.7Faz7$$\"1+v$f=A_^#Fi[m$\"1tx#y'*)[&p$Fi ]l7$$\"1+]7j_iIDFi[m$\"1G;Imj9hDFaz7$$\"1+DJS$Gga#Fi[m$\"1[+-yKjP;F/7$ $\"1++]<9VhDFi[m$\"1KC`.s\\)y%F/7$$\"1+vo%\\Mod#Fi[m$\"1,#>.S8Ay*F/7$$ \"1+](=dPAf#Fi[m$\"1[)fkplPg\"F,7$$\"1+D1\\1k2EFi[m$\"1MA5BeqRAF,7$$\" 1++DEP/BEFi[m$\"1(zWK*>VNFF,7$$\"1D\"GQZz5j#Fi[m$\"1Zdb4.j')GF,7$$\"1] iS@_6REFi[m$\"1@-C%\\)*)QHF,7$$\"1vV)*o4:ZEFi[m$\"1Y@TWWXxGF,7$$\"1+Dc ;n=bEFi[m$\"1g2B+Ek#p#F,7$$\"1D19kCAjEFi[m$\"1-xBCG$3Q#F,7$$\"1](=<@e7 n#Fi[m$\"1x6u(R**\\%>F,7$$\"1voHfRHzEFi[m$\"1l]DOr$\\R\"F,7$$\"1+](oqH to#Fi[m$\"1X.887#4Z(F/7$$\"1DJXaaO&p#Fi[m$\"13IXc+m+CFaz7$$\"1]7.-7S.F Fi[m$!1!fSdZ+gY(F/7$$\"1v$4'\\pV6FFi[m$!1!pZt/rI`\"F,7$$\"1+v=(ps%>FFi [m$!1)z)zB.J,BF,7$$\"1DcwW%3vs#Fi[m$!1#>q)**3Y;IF,7$$\"1]PM#>Wbt#Fi[m$ !1m@C)3uYk$F,7$$\"1v=#*R*zNu#Fi[m$!1\"ycm^t[:%F,7$$\"1++](o:;v#Fi[m$!1 QcYapT?XF,7$$\"1DcE7]**eFFi[m$!1oqf!\\%z5ZF,7$$\"1]7.PVPmFFi[m$!1\"yB \")G\"z]ZF,7$$\"1vozhOvtFFi[m$!1i)=bO%>NYF,7$$\"1+Dc')H8\"y#Fi[m$!1!f# =dkskVF,7$$\"1D\"G8J7&)y#Fi[m$!1ARY(pR$F,7$$\"1v$f3'4F.GFi[m$!1u%y[-b+s#F,7$$\"1+]i&G]1\"GFi[m$!1Tz$Q=fJ &>F,7$$\"1D1R5'H!=GFi[m$!1NMgH$yw6\"F,7$$\"1]i:N*3a#GFi[m$!1tfGR_EGCF/ 7$$\"1v=#*f#)yKGFi[m$\"1n$o=*f!fS'F/7$$\"1+vo%en,%GFi[m$\"10^;.3U,:F,7 $$\"1DJX4paZGFi[m$\"1?U7#eD%4BF,7$$\"1](=UBE\\&GFi[m$\"1'Gx))***\\OIF, 7$$\"1vV)*ebIiGFi[m$\"1(>1h0)zdOF,7$$\"1++v$)[opGFi[m$\"1Y\"R(HSr_TF,7 $$\"1v=U#[Ww(GFi[m$\"1qUu7'zq_%F,7$$\"1]P4\"3/c)GFi[m$\"1c/P=5Yv%F,7$$\"1DcwzOc$*GFi[m$\"1JGeFq#)RZF,7 $$\"1j:5zMa(*GFi[m$\"1yEO>$R0o%F,7$$\"1+vVyK_,HFi[m$\"1&HW0H\"*zd%F,7$ $\"1v$4r(G[4HFi[m$\"1^^#)HjO]UF,7$$\"1]7yvCW6)HFi[m$!1$=>45uhy\"F,7$$\"1D\"GQ') y!*)HFi[m$!1gK2'3P[D#F,7$$\"1++]i%Qq*HFi[m$!1WUW5k9/EF,7$$\"1Dc^[iI0IF i[m$!1n)pe%eFNGF,7$$\"1]7`MSd8IFi[m$!1z4^*e&)R$HF,7$$\"1i!Rv#zqE#F,7$$\"1+]i]2=jIFi[m$!1NbbT,O!e\"F,7$$\"1]ilAjrzIFi[m$!1hT/Xx()e\" *F/7$$\"1+vo%*=D'4$Fi[m$!1H^d[#G\"oSF/7$$\"1](=nY(y7JFi[m$!1I[#H_v=8\" F/7$$\"1++vQIKHJFi[m$!1.\"4'H&3b9*F[cm7$$\"1+D\"GI>**fSF/7$$\"1+D1fVH,KFi[m$\"1k8qxc)* y$)F/7$$\"1+]7B')o:KFi[m$\"1%*oO-f^)R\"F,7$$\"1+v=()G3IKFi[m$\"1i`x&*y .0?F,7$$\"1++D^rZWKFi[m$\"19#GUz68a#F,7$$\"1Cc,^-D_KFi[m$\"1)Gg3i\\kv# F,7$$\"1]7y]L-gKFi[m$\"1f/!)*3mW*GF,7$$\"1iSm+*4RE$Fi[m$\"1[#R-Z5#HHF, 7$$\"1voa]kznKFi[m$\"1yN#[Wo'QHF,7$$\"1(oH/+$orKFi[m$\"1?qCk*=8#HF,7$$ \"1+DJ]&pbF$Fi[m$\"1i;\")3f#f(GF,7$$\"1]P%)\\d6\"H$Fi[m$\"1m1LFi[m$\"1;w\"\\'4[l9F,7$$\"1C19\\]V9LFi[m$\"1b%G1a/M])F /7$$\"1]i!*[\"3AK$Fi[m$\"1`z:U'3Dh\"F/7$$\"1v=n[7)*HLFi[m$!1pswk)HSx&F /7$$\"1+vV[VvPLFi[m$!1f\"\\'z#ovL\"F,7$$\"1DJ?[u_XLFi[m$!1xD(['\\v)3#F ,7$$\"1](oza+LN$Fi[m$!1#3H5;H%*z#F,7$$\"1vVtZO2hLFi[m$!1]4^j**\\QMF,7$ $\"1++]Zn%)oLFi[m$!1qp[$[Co(RF,7$$\"1]i:?r(oP$Fi[m$!1)=*egH')*R%F,7$$ \"1,D\"G\\2\\Q$Fi[m$!1!GMO)y#[m%F,7$$\"1E19zE#*)Q$Fi[m$!1am7oogKZF,7$$ \"1](oa'y$HR$Fi[m$!1w')[H@GbZF,7$$\"1voz^I&pR$Fi[m$!15L*f(***>t%F,7$$ \"1+]7Q#o4S$Fi[m$!1h*fNxCCm%F,7$$\"1]7y5')**3MFi[m$!1\\w\")p9p&Q%F,7$$ \"1+vV$)*GqT$Fi[m$!1K_->'RH$RF,7$$\"1]P4c$f]U$Fi[m$!1Dm5.P??LF,7$$\"1+ +vG(*3LMFi[m$!1q([L`H5d#F,7$$\"1]iS,,7TMFi[m$!1Z;Cj\">ar\"F,7$$\"1+D1u /:\\MFi[m$!1T(Q\\e&R%)yF/7$$\"1E1Rgc;`MFi[m$!1Fgo()>E,JF/7$$\"1](=n%3= dMFi[m$\"1B#H/1MVr\"F/7$$\"1vo/Lg>hMFi[m$\"1&GWd)R)>^'F/7$$\"1+]P>7@lM Fi[m$\"1822$[\\T7\"F,7$$\"1]7.#fTKZ$Fi[m$\"1#eoW9,+.#F,7$$\"1+vok>F\"[ $Fi[m$\"1Q#zlpA:&GF,7$$\"1]PMPBI*[$Fi[m$\"1^YOz=FbNF,7$$\"1+++5FL(\\$F i[m$\"1[W4*)3S8TF,7$$\"1](oH'))=0NFi[m$\"1A$f0z:$)\\%F,7$$\"1+v$f,XI^$ Fi[m$\"1PfX9VL6ZF,7$$\"1v=U#4tp^$Fi[m$\"1wc]8dd^ZF,7$$\"1]i!*o6!4_$Fi[ m$\"1^SpuHkZZF,7$$\"1D1RX#H[_$Fi[m$\"10&\\LBR,q%F,7$$\"1+](=Kd(GNFi[m$ \"1&[iI`M,h%F,7$$\"1]P%[Z8m`$Fi[m$\"1#y0)Gt74VF,7$$\"1+D\"yipWa$Fi[m$ \"13.>syshQF,7$$\"1]7y!yDBb$Fi[m$\"1(QDGft3H$F,7$$\"1++vL>=gNFi[m$\"1/ 8;^*3Ti#F,7$$\"1](=n3Q!oNFi[m$\"1TGQ!\\h@*=F,7$$\"1+voRU*ed$Fi[m$\"1U3 `?TRF6F,7$$\"1]il#R]Pe$Fi[m$\"1#)ptJ[:AOF/7$$\"1+]iXlg\"f$Fi[m$!1N!Q@K 'oCPF/7$$\"1]Pf)pi%*f$Fi[m$!1x8Y#>$o[5F,7$$\"1,Dc^)=tg$Fi[m$!1#QKi6>Ek \"F,7$$\"1]7`/]<:OFi[m$!1J&y?9:c8#F,7$$\"1++]d6.BOFi[m$!1^9CB.%[^#F,7$ $\"1DJX*QK1j$Fi[m$!1.E9C6;nFF,7$$\"1^iS@OBQOFi[m$!1cut^&Gi!HF,7$$\"18G QPU.UOFi[m$!1(RHL\"[cMHF,7$$\"1v$fL&[$ek$Fi[m$!1omA[6sOHF,7$$\"1QfLpaj \\OFi[m$!1&G71f0S\"HF,7$$\"1+DJ&3OMl$Fi[m$!1,K?fD(z'GF,7$$\"1](=#\\&Q' oOFi[m$!1z#=^?x()[#F,7$$\"1+]785%Qo$Fi[m$!1gaNmm`/>F,7$$\"1]7.xM/*p$Fi [m$!1#Gh%*>*>l7F,7$$\"1+v$4%fC9PFi[m$!1;ik4'[J,(F/7$$\"1]P%[S[%HPFi[m$ !1gp')[ShtHF/7$$\"1++vo3lWPFi[m$!1:\\?iK**HxFaz7$$\"1(o/[!3()[PFi[m$!1 m%)p')[oCXFaz7$$\"1v$f3u!4`PFi[m$!1bpn9E-OBFaz7$$\"1iS\"pn5tv$Fi[m$!1c Q2,ED+**F[cm7$$\"1](oHhI:w$Fi[m$!1BG#ow\"fHHF[cm7$$\"1D\"y][q*pPFi[m$ \"1DnG#ozE/\"!#D7$$\"1+v=d.TyPFi[m$\"1?UQD1MbIF[cm7$$\"1]iS,,H&z$Fi[m$ \"14^bwRHOyFaz7$$\"1+]iX)p@\"QFi[m$\"1+F9Uu)QN$F/7$$\"1]P%)*e\\!HQFi[m $\"1&>g-;?==)F/7$$\"1+D1M$Hf%QFi[m$\"13oI:N8z9F,7$$\"1\\7Gy!4G'QFi[m$ \"1E=#F,7$$\"1++]A))ozQFi[m$\"129\\4miNFF,7$$\"1\\(oHUss)QFi[m$ \"1-GOA**p!)GF,7$$\"1*\\PM-c[*QFi[m$\"1!ejL:s%QHF,7$$\"1]i!RiRC!RFi[m$ \"1oUzwq='*GF,7$$\"1+]PCK-5RFi[m$\"1Mnhc`1XFF,7$$\"1]P%[#og A,?'35[#F,7$$\"1+DJD/>DRFi[m$\"1')yt)f$30@F,7$$\"1\\7yDSxKRFi[m$\"1aYm ?^vB;F,7$$\"1++DEwNSRFi[m$\"1n<(>L\\)[5F,7$$\"1](=nATz%RFi[m$\"1H)*p*R AC(RF/7$$\"1+v=F[_bRFi[m$!1G[uY@<(4$F/7$$\"1]ilF%3J'RFi[m$!1(\\@iaEq/ \"F,7$$\"1+]7G?pqRFi[m$!1U&QFzbqy\"F,7$$\"1]PfGcFyRFi[m$!1;Zy6Mt+DF,7$ $\"1+D1H#fe)RFi[m$!1&zN:NO)eJF,7$$\"1]7`HGW$*RFi[m$!18!G5fALt$F,7$$\"1 +++Ik-,SFi[m$!1jGke9k)>%F,7$$\"1+v=UN74SFi[m$!1*H(4MMS]XF,7$$\"1+]Pa1A ^TSFi[m$!1Iv_3%Fi[m$!1[4y>'>j[#F/ 7$$\"1](oz5Yg3%Fi[m$\"1iJ:@wypBF/7$$\"1+D1kY4!4%Fi[m$\"1#3QXQu1?(F/7$$ \"1^i:?K9%4%Fi[m$\"1?p\"\\x'\\&>\"F,7$$\"1++Dw<>)4%Fi[m$\"1w\"f%**HAe; F,7$$\"1+vV))))G1TFi[m$\"1yU\")H`AEDF,7$$\"1+]i+gQ9TFi[m$\"1)>JG8syG$F ,7$$\"1+D\"G6$[ATFi[m$\"1+%z5b2@\"RF,7$$\"1+++D-eITFi[m$\"14S$32)[uVF, 7$$\"1v$4@T9C(zu%4ZF,7$$\"1+vVt*G*fTFi[m$\"1^.1wJ'*=XF, 7$$\"1voaghEnTFi[m$\"1CK%>P_0>%F,7$$\"1]ilZLguTFi[m$\"1[31Pw&)RPF,7$$ \"1DcwM0%>=%Fi[m$\"1OnlBpn'=$F,7$$\"1+](=sx#*=%Fi[m$\"1UFW]u2aDF,7$$\" 1vV)*3\\h'>%Fi[m$\"1%R:*R0On=F,7$$\"1\\P4'4_R?%Fi[m$\"1O/A.I&H:\"F,7$$ \"1DJ?$G*G6UFi[m$\"1o.4_:wsVF/7$$\"1+DJqki=UFi[m$!1A#\\ykgNa#F/7$$\"1v =UdO'fA%Fi[m$!1xxg`Tg()*)F/7$$\"1]7`W3ILUFi[m$!1-)G$)*Ryv9F,7$$\"1D1kJ !Q1C%Fi[m$!1N'**H_T!p>F,7$$\"1++v=_(zC%Fi[m$!1?+do'fkO#F,7$$\"1DcE-n*f D%Fi[m$!1$p?crUEo#F,7$$\"1]7y&==SE%Fi[m$!1zX!f)eEsGF,7$$\"1i!Rv#*G!oUF i[m$!1C\"4#>8U?HF,7$$\"1voHp'R?F%Fi[m$!1cuS@9oQHF,7$$\"1)oa5T]gF%Fi[m$ !1r.>F.SGHF,7$$\"1+D\"G:h+G%Fi[m$!14*R5:b7*GF,7$$\"1]P%)>T5'H%Fi[m$!1` TeAv(o^#F,7$$\"1+](o3Z@J%Fi[m$!1:3=kz00>F,7$$\"1^i!R0!>GVFi[m$!1j)\\#z l`J7F,7$$\"1+v$4-LUM%Fi[m$!1vQ7['F/7$$\"1](oz)fFgVFi[m$!1`IA3>X'[# F/7$$\"1,++b*=jP%Fi[m$!1CwQC!es4&Faz7$$\"1v$feamQQ%Fi[m$!17e3oyzh9Faz7 $$\"1](=n899R%Fi[m$!17(eAj:#y:F[cm7$$\"1E\"yvsh*)R%Fi[m$\"1G&3AU[6'>!# A7$$\"1+vV=$4lS%Fi[m$\"1fvHzcGDGF[cm7$$\"1^i:+Xg@WFi[m$\"1Fr(oNs>;'Faz 7$$\"1+](=o*pOWFi[m$\"1'[K\\_.De#F/7$$\"1]Pfj[z^WFi[m$\"1\\\"*=2J$yM'F /7$$\"1+DJX+*oY%Fi[m$\"1Y(p2TMg<\"F,7$$\"1]7.F_)>[%Fi[m$\"1L#45UPr!=F, 7$$\"1++v3/3(\\%Fi[m$\"1&f)>@*QbS#F,7$$\"1+]7B#fG^%Fi[m$\"1WCTb.tPGF,7 $$\"1++]P!Q'GXFi[m$\"1$*p*\\5t'=HF,7$$\"1+](=&oTWXFi[m$\"19=(*zvbUDF,7 $$\"1++Dmc>gXFi[m$\"1IIm$GPeo\"F,7$$\"1+vVt]3oXFi[m$\"14BtVCm'4\"F,7$$ \"1+]i![ufd%Fi[m$\"1qsM_7-HUF/7$$\"1+D\"y)Q'Qe%Fi[m$!1fN_OA%o6$F/7$$\" 1+++&H`#>m/&fKF,7$$\"1++vB4JBY Fi[m$!1[4Ye-0QQF,7$$\"1^i!*)fP(QYFi[m$!1#p+FTi8f%F,7$$\"1+D1uU;aYFi[m$ !1F#4@w))[s%F,7$$\"1D1k6w(=m%Fi[m$!1y?;lPyPXF,7$$\"1](=#\\4fpYFi[m$!1% )z(y_zQ=%F,7$$\"1voz'G/tn%Fi[m$!1uXGMi!Rn$F,7$$\"1+]PCw,&o%Fi[m$!1*Gr \\FB`-$F,7$$\"1DJ&>'4t#p%Fi[m$!1q;E#F,7$$\"1]7`*HW/q%Fi[m$!1qR:XX K69F,7$$\"1v$4rjd\"3ZFi[m$!1dgIC)px1&F/7$$\"1+vou4(er%Fi[m$\"1QGovu!H< %F/7$$\"1EcE7VeBZFi[m$\"1ub>xOGD8F,7$$\"1^P%)\\wHJZFi[m$\"1>.&R?>k;OF,7$$ \"1](=UiuGw%Fi[m$\"1FcEcm(o`%F,7$$\"1+vVB\\-zZFi[m$\"10'zesSZt%F,7$$\" 1vo/t+5(y%Fi[m$\"1AV]oc=eXF,7$$\"1]ilA_<&z%Fi[m$\"1CmUGqE7UF,7$$\"1CcE s.D.[Fi[m$\"1`'>qF@rr$F,7$$\"1+](=_D8\"[Fi[m$\"1h3G273*4$F,7$$\"1vV[r1 S>[Fi[m$\"1X-bL`?*Q#F,7$$\"1]P4@eZF[Fi[m$\"1#))>q5&f@;F,7$$\"1DJqq4bN[ Fi[m$\"1&=***>\"[fJ)F/7$$\"1+DJ?hiV[Fi[m$\"1ze0o(GuR&Faz7$$\"1v=#*p7q^ [Fi[m$!1hOX*=e!)y'F/7$$\"1]7`>kxf[Fi[m$!1]33eE:Q8F,7$$\"1D19p:&y'[Fi[m $!1'eF/RT1!>F,7$$\"1++v=n#f([Fi[m$!1r$fbp!**[BF,7$$\"1D\"yv//P)[Fi[m$! 1y5omf+jEF,7$$\"1]iSw8[\"*[Fi[m$!1HJ$[.E!eGF,7$$\"18.#3/q`*[Fi[m$!1K\" Q+B&[6HF,7$$\"1vVB0(e#**[Fi[m$!1$=:umyl$HF,7$$\"1Q%['pt9.\\Fi[m$!1yfj> *fW$HF,7$$\"1+D1Mg.2\\Fi[m$!1(ebplvl!HF,7$$\"1](=!*)!#?7$$\"1^i:!>J,/&Fi[m$\"1F%4V1V!Q7Faz7$$\"1, v=KY3[]Fi[m$\"1M-NftdX[Faz7$$\"1](=U2Qg0&Fi[m$\"1Phor]K67F/7$$\"1++D;: *R1&Fi[m$\"1m!p,9&f$R#F/7$$\"1+DJ+%)*)z]Fi[m$\"1MMU\\bb\"G'F/7$$\"1+]P %G0e4&Fi[m$\"17+`OM&**>\"F,7$$\"1*\\P%o@r6^Fi[m$\"1!32#46_m=F,7$$\"1++ ]_!>w7&Fi[m$\"1n)yY5K?[#F,7$$\"1Cc,Om]N^Fi[m$\"1%4+C(GP:FF,7$$\"1]7`>U RV^Fi[m$\"1iP!*)HVM(GF,7$$\"1i!*G6!Qt9&Fi[m$\"1nO%=Py%=HF,7$$\"1vo/.=G ^^Fi[m$\"11\"Gxum\"QHF,7$$\"1)o/[fD_:&Fi[m$\"1s!QpT73$HF,7$$\"1+Dc'Qp \"f^Fi[m$\"1566%=5]*GF,7$$\"1]Pf`X%\\<&Fi[m$\"1d>\\1t!4X#F,7$$\"1+]i?( >2>&Fi[m$\"1#eS&H/#3`\"F,7$$\"1E19/tg)>&Fi[m$\"1%4Tm%e\\g\"*F/7$$\"1]i l()[\\1_Fi[m$\"1W9`$Q$pGAF/7$$\"1v=6T!H`PF,7$$\"1+]P\\L!=J&Fi[m$!1Eh<- D?kJF,7$$\"1DJX>70>`Fi[m$!1#p+*)H2!pCF,7$$\"1^7`*3*HE`Fi[m$!1!)=Z:YU!p \"F,7$$\"1w$4'fpaL`Fi[m$!1%f.R&R-X&)F/7$$\"1,voH[zS`Fi[m$\"1`:t_@#y/\" Faz7$$\"1Ecw*pU![`Fi[m$\"1y'y!3A-^()F/7$$\"1^P%)p0Hb`Fi[m$\"1(G.>^o*4< F,7$$\"1w=#*R%QDO&Fi[m$\"1'4k78Po[#F,7$$\"1,++5jyp`Fi[m$\"1'f]U!oszJF, 7$$\"1^7.(=$4y`Fi[m$\"1bKN4*e9%QF,7$$\"1,D1k+S'Q&Fi[m$\"19jak'4dL%F,7$ $\"1^P4Tpq%R&Fi[m$\"1\"=j&p%QMk%F,7$$\"1,]7=Q,.aFi[m$\"1L@wU$))[v%F,7$ $\"1^i:&p?8T&Fi[m$\"1j[[eQ$*pYF,7$$\"1,v=svi>aFi[m$\"1z(\\(HF!zR%F,7$$ \"1^(=#\\W$zU&Fi[m$\"1^W9'=()o&RF,7$$\"1,+DE8COaFi[m$\"1aFR&y!fsLF,7$$ \"1^7G.#[XW&Fi[m$\"1Mf\"3)GywEF,7$$\"1,DJ!3bGX&Fi[m$\"1m%*ph]X0>F,7$$ \"1_PMd>;haFi[m$\"10rT$=!z'4\"F,7$$\"1-]PM)o%paFi[m$\"1fR*=gT1*GF/7$$ \"1_iS6dxxaFi[m$!1#f68HqR\"[F/7$$\"1-vV)e#3'[&Fi[m$!1O^)=w(3#=\"F,7$$ \"1^(oaY*Q%\\&Fi[m$!1rJab`\"fy\"F,7$$\"1,+]Ujp-bFi[m$!1=W1oyUsAF,7$$\" 1,v=#)eb$[dFi[m$\"1'y.z7JV?#F,7$$\"1+D1MP[jdFi[m$\"1 Mv)*Ge$HF,7$$\"1+vo%)>\"Qz&Fi[m$\"1 jin%*RivFF,7$$\"1+++5h(*3eFi[m$\"1d!)\\')zF/7$$\"1^7GBk)o%eFi[m$!1CY J\"=kjJ*F/7$$\"1+v$f[oW&eFi[m$!1z+o_!zHn\"F,7$$\"1]Pf[00ieFi[m$!1%4vbZ UDR#F,7$$\"1++D6EjpeFi[m$!1x\"\\gud51$F,7$$\"1wo/t0yxeFi[m$!15J&3@z-p$ F,7$$\"1]P%[`Gf)eFi[m$!1q6l'=&z&>%F,7$$\"1D1k'\\wS*eFi[m$!1@-\"zc:.b%F ,7$$\"1+vVeWA-fFi[m$!1jFC#f=Ht%F,7$$\"1]7.#Q?&=fFi[m$!1Q[Q*R7s/.HF,7$$\"1 ](=#HA6^fFi[m$!1e!o;^9e2#F,7$$\"1)=<,@'=bfFi[m$!1g#ffFi[m$!1WJ,lq$)f6F,7$$\"1iS\">F/7$$\"1Q4rL@[rfFi[m$\"1z%>\"4YQXHF/7$$\"1v$4Y6cb(f Fi[m$\"1Z/B;e%**z(F/7$$\"17y]&4I'zfFi[m$\"1Y6D].0d7F,7$$\"1]iSwSq$)fFi [m$\"10wyPKe? " 0 "" {MPLTEXT 1 0 17 "# end of section " } }}}{SECT 0 {PARA 3 "" 0 "" {TEXT -1 10 "Assignment" }}{PARA 0 "" 0 "" {TEXT 306 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 b elow. " }}{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 303 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 304 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 305 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 308 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 309 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 310 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 297 51 "- You may work in groups of no more than 4 people. " }} {PARA 0 "" 0 "" {TEXT 296 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 307 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 298 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 299 51 "- I suggest you hand in a print ed Maple document. " }}{PARA 0 "" 0 "" {TEXT 300 84 "- Only hand in t he requested information. I do not want to see the Maple commands. " }{TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT 301 80 "- To enter the text edi tor in Maple, click on the \"T\" button from the tool bar. " }}{PARA 0 "" 0 "" {TEXT 302 72 "- To get a Maple prompt, click on the \"[>\" b utton from the tool bar. " }{TEXT -1 0 "" }}}}{MARK "3 11 0" 47 } {VIEWOPTS 1 1 0 1 1 1803 }