{VERSION 5 0 "SUN SPARC SOLARIS" "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 } {CSTYLE "" -1 256 "" 1 14 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 257 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 258 "" 0 1 0 0 0 0 0 1 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 1 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 "" 1 14 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 266 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 } {CSTYLE "" -1 267 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 268 "" 0 1 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 0 }{CSTYLE "" -1 271 "fixed" 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 "fixed" 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 0 }{PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Maple O utput" 0 11 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 }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 1 }3 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "" 0 256 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 }3 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "" 11 257 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 }1 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "" 0 258 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 }3 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "" 0 259 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 }} {SECT 0 {PARA 256 "" 0 "" {TEXT 256 12 "MATH 2280-1 " }}{PARA 258 "" 0 "" {TEXT 264 27 "MAPLE PROJECT : EARTHQUAKES" }}{PARA 259 "" 0 "" {TEXT 270 13 "March 6, 2006" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 318 "This handout might help for the Earthquake pro ject in section 5.3 of the text, on pages 326-328. You are mostly on \+ your own for this project, but here is a small example of a spring sys tem worked out on Maple, so that you can get an idea about useful comm ands to use. This is the example we did by hand on Friday...." }} {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 195 "This is \+ example 1 on page 317 of Edwards-Penney. Initially it is an unforced \+ system with two masses and two springs, as you can see from the descri ption on page 317. We can write the system as " }{TEXT 257 2 "Mx" } {TEXT 271 2 "''" }{TEXT 272 3 "=Kx" }{TEXT -1 8 ", where " }{TEXT 258 1 "M" }{TEXT -1 25 " is the ``mass matrix'', " }{TEXT 259 1 "K" } {TEXT -1 31 " is the ``spring matrix'', and " }{TEXT 260 1 "x" }{TEXT -1 69 " is the displacement vector. Following the book's notation, we enter" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 58 "with(linalg):with(plots):with(DEtools): #tools for pr oject" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 83 "M:=matrix([[2,0],[ 0,1]]);\nK:=matrix([[-150,50],[50,-50]]);\nA:=evalm(inverse(M)&*K);" } }{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"MG-%'matrixG6#7$7$\"\"#\"\"!7$F+ \"\"\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"KG-%'matrixG6#7$7$!$]\" \"#]7$F+!#]" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"AG-%'matrixG6#7$7$! #v\"#D7$\"#]!#]" }}}{PARA 257 "" 1 "" {TEXT -1 39 "Then the system can also be written as " }{TEXT 261 1 "x" }{TEXT 273 2 "''" }{TEXT 274 3 "=Ax" }{TEXT -1 199 "; as we discussed on Friday, the eigenvectors of \+ A determine fundamental modes, and the corresponding negative eigenval ues are the (opposites) of the squares of the corresponding angular fr equencies:" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "eigenvectors(A );" }}{PARA 11 "" 1 "" {XPPMATH 20 "6$7%!$+\"\"\"\"<#-%'vectorG6#7$!\" \"F%7%!#DF%<#-F(6#7$F%\"\"#" }}}{PARA 0 "" 0 "" {TEXT -1 286 "Therefor e, the natural frequencies of this system are the 10 and 5, and the tw o fundamental modes correspond to the masses moving in opposite direct ions (with equal amplitudes and angular frequency 10) and in parallel directions (with amplitude ratio of two and angular frequency 5). " }}{PARA 0 "" 0 "" {TEXT -1 164 " Now, let's consider the forced sy stem with force vector equal to cos(wt)[0,50], i.e. the second mass is being forced periodically. In other words, the system " }{TEXT 262 11 "Mx''=Kx + F" }{TEXT -1 307 ", where F=cos(wt)[0,50] ; this is Exam ple 3 on page 323, and we also discussed this on Friday. We follow t he method described on that page to find a particular solution to the \+ forced oscillation problem, of the form given by equation (31). Here \+ is the Maple version of the details summarized in the text:" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 452 "F0:=evalm(inverse(M)&*vector([0,50 ]));\n #The F0 in the normalized equation (32), page 323\nIden:=arr ay(1..2,1..2,identity);\n #the 2 by 2 identity matrix\nAleft:=omega ->evalm(A + omega^2*Iden);\n #the matrix function multiplying\n \+ #c on the left side of (32)\nc:=omega->evalm(-inverse(Aleft(omega))&*F 0);\n #the solution vector c(omega) to (32),\n #obtained by mult iplying both sides of equation\n #(32) on the left, by the inverse \+ to Aleft" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#F0G-%'vectorG6#7$\"\"! \"#]" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%IdenG-%&arrayG6&%)identityG ;\"\"\"\"\"#F)7\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%&AleftGf*6#%&om egaG6\"6$%)operatorG%&arrowGF(-%&evalmG6#,&%\"AG\"\"\"*&)9$\"\"#F1%%Id enGF1F1F(F(F(" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"cGf*6#%&omegaG6\" 6$%)operatorG%&arrowGF(-%&evalmG6#,$-%#&*G6$-%(inverseG6#-%&AleftG6#9$ %#F0G!\"\"F(F(F(" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 37 "c(omega ); #see equation (35) page 323" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'v ectorG6#7$,$*&\"%]7\"\"\",(\"%+DF**&\"$D\"F*)%&omegaG\"\"#F*!\"\"*$)F0 \"\"%F*F*F2F*,$*(\"#]F*,&\"#vF2*$F/F*F*F*F+F2F2" }}}{PARA 11 "" 1 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 636 "The vector c(w) above, t imes the oscillation cos(wt), is a particular solution to the forced o scillation problem we are considering. If we assume that our actual p roblem has a small amount of damping, then we expect that this particu lar solution is very close to the steady periodic solution to the damp ed problem. See the discussion on page 324. We can study resonance p henomena for these slightly damped problems by plotting the maximum am plitude of the steady state solutions to the undamped problems. That \+ would be the maximum absolute value of c1 and c2 above. Use the Maple command ``norm'' to measure this maximum amplitude:" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 15 "norm(c(omega));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%$maxG6$,$*&\"%]7\"\"\"-%$absG6#,(\"%+DF)*&\"$D\"F))%& omegaG\"\"#F)!\"\"*$)F2\"\"%F)F)F4F),$*&\"#]F)-F+6#*&,&\"#vF4*$F1F)F)F )F-F4F)F)" }}}{PARA 0 "" 0 "" {TEXT -1 100 "Another way to measure the size of c(omega) is to take its Euclidean magnitude, which is the com mand" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "norm(c(omega),2);" } }{PARA 11 "" 1 "" {XPPMATH 20 "6#,$*&\"#]\"\"\",&*&\"$D'F&-%$absG6#,( \"%+DF&*&\"$D\"F&)%&omegaG\"\"#F&!\"\"*$)F2\"\"%F&F&!\"#F&*$)-F+6#*&,& \"#vF4*$F1F&F&F&F-F4F3F&F&#F&F3F&" }}}{PARA 0 "" 0 "" {TEXT -1 436 "(Y ou will use the first command in the Earthquake project, which perhaps makes the most sense since it will be measuring the maximum amplitude that any floor oscillates.) The following picture illustrates that t he maximum amplitude of the particular solution blows up when omega is near the two natural angular frequences. Thus, in the slightly dampe d problem, one would experience practical resonance in the steady peri odic solution." }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 85 "plot(norm( c(omega)),omega=0..15,amplitude=0..15,\n numpoints=200,color=`bl ack`);" }}{PARA 13 "" 1 "" {GLPLOT2D 352 182 182 {PLOTDATA 2 "6&-%'CUR VESG6#7[\\l7$$\"\"!F)$\"3++++++++:!#<7$$\"3$=W)oPvS')y!#>$\"3Qu&>!e@M+ :F,7$$\"3X>RycG$[Z\"!#=$\"3=JRmSt>,:F,7$$\"3A!3;KkFlC#F6$\"3nb,%[F\"y- :F,7$$\"3G!4=OsPL-$F6$\"3QZ%QGMX]]\"F,7$$\"3%>`18wbkz$F6$\"3T,n(G6sz] \"F,7$$\"3J07C[^B8XF6$\"3T<]2WIH6:F,7$$\"3+Z'Hfo=aD&F6$\"3IPzd9lN::F,7 $$\"3!yb6B'f)H-'F6$\"3JTG4u!R-_\"F,7$$\"3L;NqS;4)y'F6$\"3'3*p\"H,3e_\" F,7$$\"3Tf=Puo4vvF6$\"3un7&fh%GK:F,7$$\"398Gc7:Ho#)F6$\"3[:W'4\\L'Q:F, 7$$\"3k6BY#\\k'[!*F6$\"39\"f-L(H_Y:F,7$$\"3ApRzeRy?#zF\"F,$\"3'>c:_,_ff \"F,7$$\"3vrV([#>Fe8F,$\"3C_tTIrO4;F,7$$\"32W)oP:&RH9F,$\"3-\\sH9P9A;F ,7$$\"3?d9H3'Gu]\"F,$\"3-t$3=q$=P;F,7$$\"3Z2:I5\\t\"e\"F,$\"3\"H*e-0\" \\Dl\"F,7$$\"35Ig?TdEf;F,$\"3$*[.02bsp;F,7$$\"377C[YKYIRjok&>F,$\"3/L>O3&*oZ)\\\">F,7$$\"3!HkGd9-3\\#F,$\"3LClFOn7^ >F,7$$\"3(zmLn4$fhDF,$\"3Z[1!\\!z?')>F,7$$\"3'HhAXI')*QEF,$\"3Ex!)=IOK F?F,7$$\"3x!He;G2=r#F,$\"3=*zjwU())o?F,7$$\"3?ze<&=Ezy#F,$\"31&\\uT<2c 6#F,7$$\"3XlJjEjMiGF,$\"3*pu6y/W[;#F,7$$\"3u&>RyJd-%HF,$\"32bTz*\\l0A# F,7$$\"3jd:Ji[H:IF,$\"3nuMRFxryAF,7$$\"3)eBZ%*eJ?4$F,$\"3.'*)3pUsKM#F, 7$$\"3a$pQx*G8oJF,$\"3%R`vrd))HT#F,7$$\"317C['4i!QKF,$\"3RCHhT3n#[#F,7 $$\"3kV([(\\)3#=LF,$\"3#y[7\"o\"H+d#F,7$$\"3bKlIhI*)*Q$F,$\"3ET))[l['e l#F,7$$\"3%R!3;#QBjY$F,$\"3]')GM*Gqlv#F,7$$\"3FPv]^pZRNF,$\"3@)QEfHoJ' GF,7$$\"3eBZ%*)e/2i$F,$\"3/t_&H8r^*HF,7$$\"3i7E_ao)3TXI7 $F,7$$\"3\\9Hem](3x$F,$\"3<<#*GJlH'G$F,7$$\"3M&4>QYMO%QF,$\"3f'\\nMmOF,7$$\"3US#['HA*=*RF,$\"3Q%H%fx) )[vQF,7$$\"3(Gb5@#)*4qSF,$\"3+9j:]+)=:%F,7$$\"3+(Qxa\\Wc9%F,$\"3mQT-*z .tY%F,7$$\"3v#oOt\")R6A%F,$\"3'*ybCx%[S%[F,7$$\"3mPw_0uN'H%F,$\"3&evhM B7-I&F,7$$\"3E^-0gzhoVF,$\"3AiTaC&)>TeF,7$$\"3')Ge;LNtYWF,$\"3$*p5V%4D ae'F,7$$\"3(3@U%QL=@XF,$\"3ezH$)[<>@vF,7$$\"3m=Qw_#Q&*f%F,$\"3#o'=*)Q3 N#)))F,7$$\"31kHf=$o/n%F,$\"33$4nNAAt1\"!#;7$$\"3N*)zf>$f)[ZF,$\"3ie03 $4*)HQ\"F__l7$$\"3)3@U%QU%R#[F,$\"3QwFaY:**[>F__l7$$\"3'fLn%=`Rh[F,$\" 3s39[)\\Z0Y#F__l7$$\"3%>Y#\\)RY))*[F,$\"34'y4a*=#4N$F__l7$$\"3xkIh(H0% =\\F,$\"3UbnS'3D39%F__l7$$\"3[oOt'>kz$\\F,$\"33tmoTO#)GaF__l7$$\"3T>Q^ @R&G%\\F,$\"32*[Zab%e))eF__l7$$\"3MqRHYOuZ\\F,$\"3kJg&[M(QMkF__l7$$\"3 F@T2rLj_\\F,$\"3%e>:BnyG4(F__l7$$\"3JrU&e4Bv&\\F,$\"34@$HA\"\\(H!zF__l 7$$\"3N@Wj?GTi\\F,$\"3c%oR#**H%Q#*)F__l7$$\"3GsXTXDIn\\F,$\"3eKE'GA/]- \"!#:7$$\"3uZY!ySZ(p\\F,$\"3[^UQD2R26F\\cl7$$\"3@BZ>qA>s\\F,$\"3Q$\\91 8kU?\"F\\cl7$$\"3n)z%eKrju\\F,$\"3o$R<#[Q\")>8F\\cl7$$\"39u[(\\*>3x\\F ,$\"3$QAI0(p,g9F\\cl7$$\"3D2!*[-QLz\\F,$\"3@C@$f?'\\=;F\\cl7$$\"3PSJ+5 ce\")\\F,$\"3+o0!yzMd\"=F\\cl7$$\"3\\ts^%*z\"F^gl7$$\"3drW9&o.T+&F,$\"3%[Q^F,Ds6)F\\cl7$$ \"3p/'eE\\bj+&F,$\"3Z:o,D=DR_F\\cl7$$\"3\"yts,I2'3]F,$\"3`;N#QREr'QF\\ cl7$$\"3$4(oo2\"f3,&F,$\"3Ex\"e/wkS1$F\\cl7$$\"30/5?:468]F,$\"3_y+@$yb o`#F\\cl7$$\"3;P^rAFO:]F,$\"3_zECn#*>k@F\\cl7$$\"3Gq#H-`9w,&F,$\"3&oRn z9Ao)=F\\cl7$$\"3S.MuPj')>]F,$\"3%\\U)fo_Ks;F\\cl7$$\"3_OvDX\"=@-&F,$ \"3i%z!elM],:F\\cl7$$\"3kp;x_*pV-&F,$\"379\"z!z)\\AO\"F\\cl7$$\"3w-eGg )=vdB7uF__l 7$$\"3%[8Fa$)R\"\\]F,$\"3?\\+w#*H#ys'F__l7$$\"3N>T2rWOe]F,$\"3M*p;Wm*p bcF__l7$$\"3u/6s1\"*en]F,$\"3l#>iW&fAw[F__l7$$\"38!4oBu8o2&F,$\"3&G!>P #3oRG%F__l7$$\"3ju],y$Qg3&F,$\"3r[c@](3(=QF__l7$$\"3`W!4$\\w[/^F,$\"30 i#3^aTY8$F__l7$$\"3U9Ig?p$H7&F,$\"3f(f(fD>*el#F__l7$$\"35**)zfHxO;&F,$ \"3^.w;wZ+\")>F__l7$$\"3z$yc8ny_F, $\"3[B=Y#))zE9\"F__l7$$\"3lV([(*pRON&F,$\"3s$e9Xi8-())F,7$$\"3?v],`#o. V&F,$\"3F/aFXeN*=(F,7$$\"3Rf=PCn*3]&F,$\"37!=F,7$$\"3m*******z`_S'F,$\"3wJ9n^7<(z\"F,7$$\"3::JiucK&[' F,$\"35GF*Q<+jm\"F,7$$\"3'GrU&3,halF,$\"3q=g6bwYj:F,7$$\"3BV([(*z&=MmF ,$\"3/hIA5Zbb9F,7$$\"3%>d9H8\"**3nF,$\"3z-]hb%)pi8F,7$$\"3(GqS\"GR2$y' F,$\"3PE%35[RyF\"F,7$$\"3')=Rycc@doF,$\"3kCG8&>'3*>\"F,7$$\"3!*F67$$\"3E\\+,- %H>p(F,$\"3-es<*[>(f*)F67$$\"3,d:J7'[0w(F,$\"3kwk#4^#p@*)F67$$\"3AtZ&4 4=x$yF,$\"3Q!=lt0Mh*))F67$$\"3&Q(\\**)4*R:zF,$\"3-eCL6z.*)))F67$$\"3-> Qw-4r#*zF,$\"33zG=OH+,*)F67$$\"3H&Gd9%))Qk!)F,$\"3GTeQ'*\\nH*)F67$$\"3 uIhA&>2'Q\")F,$\"3[NE3bW/y*)F67$$\"3/?V'G#RO:#)F,$\"3@_G'QUj#\\!*F67$$ \"3!p^.2\\u=H)F,$\"3qc%GNc!HV\"*F67$$\"3v].29]dq$)F,$\"3`dBNn*QhE*F67$ $\"3-Y%*)yZ%*)R%)F,$\"3KN>/-&e%)R*F67$$\"3K&RzexJz^)F,$\"3wx%R@F*[x&*F 67$$\"3%=c7D]*G'f)F,$\"3Gz\"=?WYLz*F67$$\"3S'Gd9p,=n)F,$\"3)3:\"[%G^S+ \"F,7$$\"3+!3;K*RPS()F,$\"3%Gd(400LI5F,7$$\"395BY#48>#))F,$\"3q:=M]k8n 5F,7$$\"3iB^-0u)4*))F,$\"3>L]mT.*Q5\"F,7$$\"3pOv]^sLr*)F,$\"3J&[yOC@W: \"F,7$$\"3*)4?S![gC/*F,$\"3,%ePeX]w?\"F,7$$\"3zAY#\\$R\\?\"*F,$\"3eWxO r_`x7F,7$$\"3urY$pB+[>*F,$\"33w`x.yDe8F,7$$\"3.&>Ry1JBF*F,$\"304b]d.yh 9F,7$$\"3Fxb6t&GNM*F,$\"3n%y/-()4+e\"F,7$$\"3]D]+,ZK?%*F,$\"3![MI3Ql>u \"F,7$$\"3)ohBZH%4+&*F,$\"3'\\0whqO^'>F,7$$\"3GkFbgR`p&*F,$\"3Z3P0$45) GAF,7$$\"3p\\+,-2`W'*F,$\"3Y>owPixJEF,7$$\"3\"y&=Pu'4?s*F,$\"3h\")*fj4 3(zKF,7$$\"33He;Lx!yz*F,$\"3!efIoUL$)R%F,7$$\"3]hD^_i9r)*F,$\"3[\"z\"= `VHRnF,7$$\"3s4AW87'=\"**F,$\"3PVJt!e$=C(*F,7$$\"3$z&=Puhd_**F,$\"3`@i H!o')Qy\"F__l7$$\"33U%)=DAsh**F,$\"3)))f*f*)Rp.AF__l7$$\"3YC]+w#o3(**F ,$\"3))=$[RhQr)GF__l7$$\"3k:LT,8Wv**F,$\"3?L!*ye6x>MF__l7$$\"3$og@oK9+ )**F,$\"3__K6mb;'>%F__l7$$\"3U_d_R3I#)**F,$\"3=;#3^6>[t%F__l7$$\"3-)*) HAN(e%)**F,$\"3IJg***F,$\"3Y<4V c$*R'4#F\\cl7$$\"3erZXGkI)***F,$\"3wSVLjv?B\\F\\cl7$$\"3s\"*e6%Hf++\"F __l$\"3(e`u2.EeS\"F^gl7$$\"3G1jQXzG+5F__l$\"3]_ID]<4)*GF\\cl7$$\"3%3sc mf;0+\"F__l$\"3Jl\"R%yN9<;F\\cl7$$\"3dNr#zCX2+\"F__l$\"3C+P&\\A:A7\"F \\cl7$$\"3J]v>**Q(4+\"F__l$\"3CfGC\"*Q$of)F__l7$$\"3(['zY]D?,5F__l$\"3 sj$yGAa)ppF__l7$$\"3Vz$QsW*)=g$G1\"F__l$\"3)Q^S![rvk;F,7$$\"3 -S!3;p402\"F__l$\"3R[MEL/l8:F,7$$\"3y&=PC#)>\"y5F__l$\"3q*pIq$p*>R\"F, 7$$\"3jd:JUF6&3\"F__l$\"3/.,-O%G&)H\"F,7$$\"3%3>QwTFJ4\"F__l$\"3a>-N_# zx?\"F,7$$\"3xpRzQeH+6F__l$\"3U8krx_=Q6F,7$$\"3t'Rz3()Qz5\"F__l$\"3.oH ![%[]t5F,7$$\"3QqS\"yAa_6\"F__l$\"34'QQK]e!>5F,7$$\"3y)yd:*pPB6F__l$\" 3Z9h2'R&zb'*F67$$\"3uxb632TI6F__l$\"3G383M5*>C*F67$$\"3#zg@$RSRQ6F__l$ \"3&pdA@Lw'=))F67$$\"3(eAX!z*pc9\"F__l$\"3MF/FlKRp%)F67$$\"3&f@VOrGO: \"F__l$\"3')*)y)[kb87)F67$$\"3cS\"Gcv&\\g6F__l$\"3ds#ed)H\"e%yF67$$\"3 nrV([^;$o6F__l$\"3w(R\"G@:IcvF67$$\"3\"F__l$\"3!3Goh;4 q$oF67$$\"3oT$o'G$o\")>\"F__l$\"3]:v0SUdPmF67$$\"3a**)zf))zf?\"F__l$\" 3^E/3a**fNkF67$$\"3cPv]m[U87F__l$\"3'Rq,6ag\\D'F67$$\"3_)pRzNg7A\"F__l $\"3$*y,wFQ6wgF67$$\"3)GhAXO`$G7F__l$\"3Dl$o*))=EBfF67$$\"3S:JikC>O7F_ _l$\"3'H;CMIIMw&F67$$\"3lPv]c4qV7F__l$\"3ZR=Hpn_=cF67$$\"3%Gc7D<\">^7F __l$\"3ONzVs3^F67$$\"3n%**) z*H[QE]qmH\"F__l$\"3'*H(y1,0\\x%F67$$\"33C['zNVVI\"F__l$\"3w8t8ZDFt YF67$$\"3]-050iR68F__l$\"3;$[6whINe%F67$$\"3!Rze@vKM] \"\\%F67$$\"3W-05!>vjK\"F__l$\"3mYb\"H:fOS%F67$$\"3\"*>S!3,cVL\"F__l$ \"3AWc_#*fH8VF67$$\"3dLnMR()RT8F__l$\"3=;`Ff\\YOUF67$$\"3*=S!3'eC&\\8F __l$\"3*o$GK'[25:%F67$$\"3&pS\"GOv&oN\"F__l$\"3)yT)o2qlwSF67$$\"3>b5@2 J6k8F__l$\"3E*3?lOKb+%F67$$\"3'3P\"F__l$\"3'zXzG8:;$RF67$$\"3w&= Pu!)[(z8F__l$\"3p>kI'yH*fQF67$$\"3U4>Q^i(oQ\"F__l$\"3iPN'QS&z'z$F67$$ \"39c7D!GXWR\"F__l$\"3L&Q$fN%==t$F67$$\"3W;Lm7>$=S\"F__l$\"3)4N-HVi.n$ F67$$\"3*ziD,5R)49F__l$\"3%3*QnTu\"eg$F67$$\"3n(eQE6DZU\"F__l$\"3!zx4_d(>\"\\$F67$$\"3%Qxafk0AV\"F__l$\" 3C+2k5Z2OMF67$$\"3%pQxa#RhR9F__l$\"3c2WuV)>IQ$F67$$\"3O3o)o\\h;!G$F67$$\"3N(\\**[m0E Y\"F__l$\"3?YxFZiRFKF67$$\"3rg@Vm^&*p9F__l$\"3Ky#3P#)H.=$F67$$\"32He;o +@x9F__l$\"3lSg\"zHb]8$F67$$\"3)Gf=(Qa'\\[\"F__l$\"3#>g1Gd3z3$F67$$\"3 QPv]ERq#\\\"F__l$\"3#z+niR,@/$F67$$\"#:F)$\"3))**************HF6-%+AXE SLABELSG6$Q&omega6\"Q*amplitudeFf`p-%'COLOURG6&%$RGBGF)F)F)-%%VIEWG6$; F(F^`pF_ap" 1 2 0 1 10 0 2 6 1 4 2 1.000000 45.000000 45.000000 0 0 "C urve 1" }}}}{PARA 0 "" 0 "" {TEXT -1 208 "This is qualitatively the pi cture on page 323, figure 5.3.10, although they plotted the Euclidean \+ magnitude of c(omega) rather than the maximum amplitude. Notice how w e get Maple to label the axes as desired" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 92 " We can get a plot of resonan ce as a function of period by recalling that 2*Pi/T=omega:" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 86 "plot(norm(c(2*Pi/period)),period=0. 1..3,amplitude=0..15,\nnumpoints=200,color=`black`);" }}{PARA 13 "" 1 "" {GLPLOT2D 305 183 183 {PLOTDATA 2 "6&-%'CURVESG6#7i[l7$$\"3/+++++++ 5!#=$\"3qDJ1$QnGG\"!#>7$$\"3FV'GdaqC:\"F*$\"3_K:?ZR?6p8#F-7$$\"3IU%)ox'GVV\"F*$\"3>vh+ByFwEF-7$$\"3](\\**)f >^%e\"F*$\"3ZIu#\\cXgG$F-7$$\"3mHf=Z9)Rt\"F*$\"3=!Q76XP>'RF-7$$\"3;LmK &zeD(=F*$\"3NRUe!>RDl%F-7$$\"3`lIh7w/;?F*$\"36'y[nQy_V&F-7$$\"3#yc8F&R Wk@F*$\"3O\\w$\\\\J1K'F-7$$\"3nY$pQQkBJ#F*$\"3%H%)pH\")eRG(F-7$$\"3RE_ /H(=XY#F*$\"3q8#QC(=,k$)F-7$$\"3Uxa4fj`)f#F*$\"3/o\"RJej[R*F-7$$\"3O8F ao%3%\\FF*$\"3CWC()>p0l5F*7$$\"3UMoOt+!4!HF*$\"35AHf\"**R??\"F*7$$\"3H LmK0.*o/$F*$\"38'G**R!GBX8F*7$$\"3I,.1AMYzJF*$\"3Xl#*fW9f&[\"F*7$$\"3E :Ig!o0rL$F*$\"3)GO?+;'\\m;F*7$$\"33*yd:N\\1Z$F*$\"3D]QS`z$G$=F*7$$\"35 _/4)Q#*fi$F*$\"3+>LFOvO9RF*$\"3Rv&z#GM!**[#F*7$$\"3I9HeE3-eSF*$\"3Vi'em$p.XFF*7$$\"3.f;Lm P\"z?%F*$\"3uP4yP]IRIF*7$$\"3-kE`wAcXVF*$\"3*oh)px2DRLF*7$$\"3]hAX!zMS \\%F*$\"33&*Q\"))yO5q$F*7$$\"3Xu[(\\+c#[YF*$\"3<22$3D\"QGTF*7$$\"3^*)y dDg]#y%F*$\"39K\"ek$G*Qb%F*7$$\"3[V'Gds*\\F\\F*$\"3O*)y:eJ?&3&F*7$$\"3 *49GcS#Hx]F*$\"33G<2:A]PdF*7$$\"3f#['H>`$QA&F*$\"3k)e'pN2U;lF*7$$\"3*Q $oOBJil`F*$\"3MX/?(e\"3iuF*7$$\"39T\"GcIaI_&F*$\"3Gvn[@fLd))F*7$$\"3;A Y#\\A8Xm&F*$\"3#oM1,C$\\j5!#<7$$\"3_w`2:3b:eF*$\"3ZF8F$zv'f8Fhu7$$\"3z mNrn?)R)eF*$\"3k`20HXQi:Fhu7$$\"3=eBFhu7$$\"3B>Pu)=S?5'F*$\"3t&4!zgo\\xJFhu7$$\"3iu[ (**zLC<'F*$\"3U1'eDXOq-&Fhu7$$\"3!*Gg?6u#GC'F*$\"33c!H-'4fG8!#;7$$\"3q Px/7^-_iF*$\"3SuN+,ar6]F'F*$\"3vFdl\"4vPW'F\\x7$$\"3%H'Gd9#='ziF*$\"3g7Rn$\\Y5Z\"! #:7$$\"3*yr$*\\1`jF*$\"3%=4['*4vD$yFhu7$$\"3@mLnC1***Q'F*$\"3? #pO'\\e@g_Fhu7$$\"3NU([Z4I>Y'F*$\"3K0_1&=\\XH$Fhu7$$\"3i>T#[cpQ`'F*$\" 34w>=RE]hCFhu7$$\"3r[**)H_$=4mF*$\"3>&=(o?-7()>Fhu7$$\"3zxd:\"[(\\%o'F *$\"3#G?\")yz\"3%p\"Fhu7$$\"3pGd9uP.dnF*$\"3@sA,wnl-:Fhu7$$\"3\\yc8n+d HoF*$\"3td$*38**[k8Fhu7$$\"3&**)zfRx#z(pF*$\"3!GROH3V!y6Fhu7$$\"3CT\"G c$p0DrF*$\"3(f4e]RHa1\"Fhu7$$\"3\\jE`'Q`-E(F*$\"3D)R2[7Y#z**F*7$$\"3?/ 4=wP?:uF*$\"3LaJ!Gzy6Z*F*7$$\"3^'Hf=D$z`vF*$\"3S'**e$zL*3=*F*7$$\"3#4A W)Q&e:q(F*$\"3M\">='3-B%**)F*7$$\"33tX\"Hx))H%yF*$\"3u.pqd^v1*)F*7$$\" 3&*)zf>()G++)F*$\"3e#y],u$4$*))F*7$$\"3'zrV(eS,O\")F*$\"3/;F))f4!4%*)F *7$$\"3-,.1i%e.H)F*$\"3#\\yG\"=\"Fhu7$$\"3C+,-%)4H1$*F*$ \"3r*zRg$4:88Fhu7$$\"3d=QwAZ*fW*F*$\"3%eZu=SBJV\"Fhu7$$\"39.17k\"=qf*F *$\"3o.[(=&*3Vd\"Fhu7$$\"3c2;KaW&4u*F*$\"3!H4[)\\,)=s\"Fhu7$$\"3E]+,A1 W#*)*F*$\"3=!z`6BxL*=Fhu7$$\"3_18E3s&H+\"Fhu$\"3wtmP.\\^l?Fhu7$$\"3qW* )y,G6=5Fhu$\"3P&f9^%f,yAFhu7$$\"3!4;K%>#HE.\"Fhu$\"3q\\-^93$y]#Fhu7$$ \"3uU&3PI5r/\"Fhu$\"3Wn![265#oFFhu7$$\"3<4=O_eBi5Fhu$\"3'fdDuS553$Fhu7 $$\"3DY#\\[.nh2\"Fhu$\"3o!)3\\n[,:MFhu7$$\"3+$e;8ZM/4\"Fhu$\"3E\"*[Cq. o9QFhu7$$\"3c6BY;u=16Fhu$\"3g+dTMs`VVFhu7$$\"3M\\)pzT]/7\"Fhu$\"3+pC85 3_F\\Fhu7$$\"3a!4='3q.N6Fhu$\"3wi()*fi(4mcFhu7$$\"3i9He&>r)\\6Fhu$\"3M kWv'RfQi'Fhu7$$\"3PE_/+n]j6Fhu$\"3SLCil!*3txFhu7$$\"3Zc7D1\"Fhu$\"3be[JW,J$=\"F\\x7$$\"3n$pQxl&*y?\" Fhu$\"3OC$p\\/\"f'f\"F\\x7$$\"3g;LmGOq97Fhu$\"3Fvijs!Rc(=F\\x7$$\"3vRz e*f6:A\"Fhu$\"3x$*Q%QXaGE#F\\x7$$\"35>Qw%em$H7Fhu$\"3$y%pB1\"R)\\HF\\x 7$$\"3X)pR*p:AP7Fhu$\"338\">CK&p#>%F\\x7$$\"3_X7Fhu$\"3Y7` `J6F+sF\\x7$$\"38Y#*4=A3Y7Fhu$\"3-Q>.#H#)\\\"yF\\x7$$\"3!4?:HGopC\"Fhu $\"3c4E=b,OU&)F\\x7$$\"3nb6tZV&yC\"Fhu$\"3WKd(zm/lT*F\\x7$$\"3V5ra7/u[ 7Fhu$\"3([]$R)3#o[5F[z7$$\"3?lIOxki\\7Fhu$\"3e#Q[$*4pF=\"F[z7$$\"3fU5x 4&p+D\"Fhu$\"3#\\=J!z@a70D\"Fhu$\"3c]'>g@acN\"F[z7$$ \"3N(*peub&4D\"Fhu$\"3PIz;0(=BY\"F[z7$$\"3'\\(\\*pg)R^7Fhu$\"3%39H&*4D qe\"F[z7$$\"3+3`Skp$=D\"Fhu$\"3;(Q2:UsIt\"F[z7$$\"3$3k:=KvAD\"Fhu$\"3( R-:*QaZ3>F[z7$$\"3ltfAzOr_7Fhu$\"3[ScsSQ2B@F[z7$$\"3[1jjO?:`7Fhu$\"3Q& >eV)zl\"R#F[z7$$\"3IRm/%R!f`7Fhu$\"3cd#o,VIvt#F[z7$$\"3MspX^(GSD\"Fhu$ \"3#z6yW#)e'*>$F[z7$$\"3R0t')3rYa7Fhu$\"36\\\"p1h)\\[QF[z7$$\"3@QwFma! \\D\"Fhu$\"3'*eps!3\\e#[F[z7$$\"3/rzoBQMb7Fhu$\"3C`kabtwlkF[z7$$\"33/$ )4\"=#yb7Fhu$\"3wD>V\"4*[(y*F[z7$$\"38P'3&Q0Ac7Fhu$\"3a\"3#*R$445?!#97 $$\"3&*p*=f*)emD\"Fhu$\"3'yp&42*zp$Q!#87$$\"3y-$HLD(4d7Fhu$\"3*HEH>%fn @=F]bm7$$\"3#ejR2hNvD\"Fhu$\"3-G]rhLnN$*F[z7$$\"3()o*\\\"oR(zD\"Fhu$\" 3ql/5yNhyiF[z7$$\"3p,.cDBTe7Fhu$\"3i1[f\"*QGJZF[z7$$\"3^M1(Ho])e7Fhu$ \"3EF%4F:\"y'z$F[z7$$\"3cn4QS!*Gf7Fhu$\"3xkEfvY@rJF[z7$$\"3h+8z(RF(f7F hu$\"3;*G^Tn=Js#F[z7$$\"3VL;?bd;g7Fhu$\"3_#)4??NN'Q#F[z7$$\"3Dm>h7Tgg7 Fhu$\"3#*pv'[=8S7#F[z7$$\"33*HA+ZU5E\"Fhu$\"3^E7.%)3)Q\">F[z7$$\"3!>jK u#3[h7Fhu$\"3]Tt\"oC#yT>E\"Fhu$\"3a'y#)>;[#)f\"F[z7 $$\"3*zH`AadBE\"Fhu$\"3Mz*o,(>rw9F[z7$$\"3\"3jj'**ezi7Fhu$\"3Nj-*)[lZs 8F[z7$$\"3kjR2dUBj7Fhu$\"3)[X(3'>$4#G\"F[z7$$\"3Y'H%[9Enj7Fhu$\"3w#o( \\4F(H?\"F[z7$$\"3Li\\IH$\\XE\"Fhu$\"3/s]JyC.r5F[z7$$\"3?Gc7WgUl7Fhu$ \"3:2-H$pZTl*F\\x7$$\"3&y30UzojE\"Fhu$\"3InvI()F\\x7$$\"3tZXGW:Jn7 Fhu$\"3p&yv%R.ZqzF\\x7$$\"3g2SO%Ha#o7Fhu$\"35z@Gd-eLtF\\x7$$\"3EnMWWq> p7Fhu$\"3;4'ye;+Bz'F\\x7$$\"3,(Q-Ya#3r7Fhu$\"30*3KB(Gp@fF\\x7$$\"3a18w W!oHF\"Fhu$\"395)z*3\">@D&F\\x7$$\"3#e9z]/RnF\"Fhu$\"3.#\\)HYoz*G%F\\x 7$$\"3([)pRX+^!G\"Fhu$\"3]*Hz25'\\JOF\\x7$$\"3a**)z\\b#3)G\"Fhu$\"3j!* HnE=W'y#F\\x7$$\"3?9Gck]l&H\"Fhu$\"343=i(H2%pAF\\x7$$\"308E_*yM%48Fhu$ \"3:L7T7!p!4&oZ6F\\x7$$\"3xNrUI'HQN\"Fhu$\"3=2kTq+SW)*Fhu7$$\"3![#\\ )4iCsO\"Fhu$\"31^c$R\\c3!))Fhu7$$\"3e!4=YD4EQ\"Fhu$\"3/F#\\&ehHwyFhu7$ $\"3'3:I!>;2(R\"Fhu$\"39j-h$QM?>(Fhu7$$\"3TpQx#H%R69Fhu$\"3K_Uyh[cSmFh u7$$\"3!yb6.OGdU\"Fhu$\"3WiPmaBK#='Fhu7$$\"3m@V'=?[.W\"Fhu$\"3-B%p`bL( )y&Fhu7$$\"3&)\\**)*))\\%eX\"Fhu$\"3O'eo.QUZV&Fhu7$$\"3DKkGDS0q9Fhu$\" 3V\"e'R^0Yb^Fhu7$$\"3.\">Q'o,3%[\"Fhu$\"3NS#p=:,T\"\\Fhu7$$\"3+,-/&)Q2 *\\\"Fhu$\"3!>E-(eR3(o%Fhu7$$\"3)3;K9'\\.9:Fhu$\"3QQCh@*Gq[%Fhu7$$\"3) )[(\\fTju_\"Fhu$\"3,@3,[qPEVFhu7$$\"3anMp)QqLa\"Fhu$\"3hPE&>JWc:%Fhu7$ $\"3)f>Reg!oc:Fhu$\"31(zwg^vn-%Fhu7$$\"3lX\"H)*R[Ad\"Fhu$\"3=CIOsC\\5.j\"Fhu$\" 3K!>D)[y$3[$Fhu7$$\"3D05?TuDX;Fhu$\"39#)z,&)y[&R$Fhu7$$\"3Iv],w^6f;Fhu $\"3DDz$\\_$3ALFhu7$$\"3v=PudSYt;Fhu$\"3zKh2%=#G^KFhu7$$\"3\")oPv\"p.$ )o\"Fhu$\"3'erqgrtI=$Fhu7$$\"3tY$p[t&4.JFhu7$$\"3%[ $pQp6J=RKHrJQI\\>w\"Fhu$\"3%[V45#3l/HFh u7$$\"3[v],F$[lx\"Fhu$\"3N^;'zty*eGFhu7$$\"3AU%)oQc!)*y\"Fhu$\"3c;&=s: 8(>GFhu7$$\"3i8Fak)pb!=Fhu$\"3=7'>FVTbx#Fhu7$$\"3%3>Q;BC*==Fhu$\"3MFxj %\\,,u#Fhu7$$\"3_d9HN&eW$=Fhu$\"39$*)Hmy$)4q#Fhu7$$\"3Jb5@E!4#[=Fhu$\" 3;srf_B6oEFhu7$$\"3:Fa3%\\&Hj=Fhu$\"3w_ql/O!Qj#Fhu7$$\"3_.29z8mx=Fhu$ \"3VrPkTEs-EFhu7$$\"3)yd:Jn]E*=Fhu$\"3p`M/!oH=d#Fhu7$$\"3hyc8CbT1>Fhu$ \"3(4UW=h]Za#Fhu7$$\"3NQw_vFE@>Fhu$\"3*3,Y.#f$o^#Fhu7$$\"3a**)zp*[oO>F hu$\"34)y!3!op\"*[#Fhu7$$\"3E,-/**)4,&>Fhu$\"3FKt$)\\\">hY#Fhu7$$\"3Lw _0p#4Y'>Fhu$\"3,UXENvBUCFhu7$$\"3(fAXq`)ez>Fhu$\"3O%4e8,0'=CFhu7$$\"3e g?TGGC%*>Fhu$\"3gzu3XgW'R#Fhu7$$\"3P&4>)3;U3?Fhu$\"3$RD8NndeP#Fhu7$$\" 35E_/FZ;C?Fhu$\"3bpz(=s@RN#Fhu7$$\"33u[(*=1JQ?Fhu$\"3Cjv=VN*\\L#Fhu7$$ \"3j\\**)zP9M0#Fhu$\"3;RE'HeabJ#Fhu7$$\"3p(e<&G15n?Fhu$\"3#HK(G74f)H#F hu7$$\"3d$yc`Jj?3#Fhu$\"3!yXM-[:2G#Fhu7$$\"3K:IgP?9'4#Fhu$\"3lAY+\"p-X E#Fhu7$$\"3%*[(\\*e$e36#Fhu$\"3IZSn%*=:[AFhu7$$\"3%R#['HDY_7#Fhu$\"3CF 4BncsKAFhu7$$\"3U!*zfW!4.9#Fhu$\"354#**\\%)Qr@#Fhu7$$\"3szf>.j\"[:#Fhu $\"3Ub8!HWWE?#Fhu7$$\"3'3@U/2_'p@Fhu$\"3T>bxe7K)=#Fhu7$$\"3@E_/!*\\O%= #Fhu$\"3VHi_M:fu@Fhu7$$\"3Wyc8NY)y>#Fhu$\"3m\"[ZQ)HPi@Fhu7$$\"3H-05u'z L@#Fhu$\"3(H%R0!)e\")[@Fhu7$$\"3'>Mo;iQsA#Fhu$\"3!R'=SIk2P@Fhu7$$\"3'R $oO]^,UAFhu$\"3ap/p5K%\\7#Fhu7$$\"3ipQxt\"ehD#Fhu$\"3sH:o,[o8@Fhu7$$\" 35#RyO=i=F#Fhu$\"3w'z%)Q5s:5#Fhu7$$\"3*QycBqgaG#Fhu$\"3G\\#y5S*R\"4#Fh u7$$\"3`U%)oU^*3I#Fhu$\"34k!)=NK>!3#Fhu7$$\"3)oV([f>'\\J#Fhu$\"3Bf0Ta= Gq?Fhu7$$\"3W%)oPY)[.L#Fhu$\"3p.#Fhu7$$\"3\"3)43\"y*>Fhu7$$\"3&Gd9>uafW#Fhu$\"3-Fwp$f9**)>Fhu7$$\"3AFhu7$$\"3]=Pu/l\"[Z#Fhu$\"37dA_TPou>Fhu7$$\"3oc8F)4s **[#Fhu$\"3:Ic)=E#*p'>Fhu7$$\"36tX\"f^)[/DFhu$\"3!)[y*z_8)f>Fhu7$$\"3& \\&4>+'p*=DFhu$\"3$GmYjMHG&>Fhu7$$\"3;@U%)[^4MDFhu$\"3Qlh;FlrX>Fhu7$$ \"3Be;LJj-[DFhu$\"35P6fsVKR>Fhu7$$\"3)\\**)znPHiDFhu$\"3(y,IFhIH$>Fhu7 $$\"3JBZ%HrY!yDFhu$\"33HE.2X/E>Fhu7$$\"3ahAX9(4Bf#Fhu$\"3!>?/OKh*>>Fhu 7$$\"3_-050j*og#Fhu$\"3y16zDT)Q\">Fhu7$$\"3gE`1#\\IFhu7$$\"3OQw_'*fONEFhu$\"3CYv=:WU->Fhu7$$\"3ooOt-5!*\\EFhu$\"3%p!4lZ) pn*=Fhu7$$\"3PQw_n`KkEFhu$\"3QzJ+&*QG\"*=Fhu7$$\"3k06Aa\\vzEFhu$\"3.5z *3n\\b)=Fhu7$$\"3%>Nqg*3P$p#Fhu$\"3m0M)G(4g!)=Fhu7$$\"3m5@Um334FFhu$\" 3hp8&z)p,v=Fhu7$$\"3%pQxM!zDBFFhu$\"3*)[2YG+4q=Fhu7$$\"3SS!31M&GPFFhu$ \"3S._*3f;`'=Fhu7$$\"33(Rz=MpBv#Fhu$\"3]#oNf%HHg=Fhu7$$\"3TE_/hV^nFFhu $\"3#H(eTK&f`&=Fhu7$$\"3DD]+'3%H\"y#Fhu$\"35)[n?Ck4&=Fhu7$$\"3k&4>=aFf z#Fhu$\"3dB5K:7RY=Fhu7$$\"327C[k$3-\"GFhu$\"3ZmIv^'>?%=Fhu7$$\"3wZ&4p# *)oDGFhu$\"3`#G7vozt$=Fhu7$$\"3yOtY%HFhu$\"3I,>LD\\`0=Fhu7$$\"3B.17l%Rf&HFhu$\"3<&*3N \"= " 0 "" {MPLTEXT 1 0 27 "E*(omega)^2*cos(omega*t)*b;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#**%\"EG\"\"\")%&omegaG\"\"#F%-%$cosG6#*&F'F%%\"tGF %F%%\"bGF%" }}}{PARA 0 "" 0 "" {TEXT -1 852 "where b is the transpose \+ of [1,1,1,1,1,1,1]. They ask, ``how can the earthquake be forcing al l seven stories, it seems like it's just shaking the bottom one.'' We ll, the students are correct, but so is Edwards-Penney. The authors t alk about an ``opposite inertial force'' being the reason for this for cing term and there's a detailed discussion of this on page 4 of Frid ay's notes. here's a brief summary. Tink of the ground as the zeroth story. In the rest frame it is shaking with oscillation Ecos(wt). A nd so its acceleration is its second time derivative, namely -E*w^2*co s(wt). If you write down the inhomogeneous system of EIGHT second ord er DE's for the accelerations of stories zero thru seven, the forcing \+ (well, accelerating) term is -E*w^2*cos(wt)*[1,0,0,0,0,0,0,0], as you \+ would expect. Call the solution 8-vector to this system " }{TEXT 266 1 "y" }{TEXT -1 79 "(t), then see what the shaking looks like to someo ne on the ground by letting " }{TEXT 267 1 "x" }{TEXT -1 4 "(t)=" } {TEXT 268 1 "y" }{TEXT -1 194 "(t)-E*cos(wt)*[1,1,1,1,1,1,1,1]. Then \+ the zeroth story component of x(t) will be identically zero, and the o ther seven components will satisfy equation (2) on 327, exactly as the authors claim." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT 275 20 "Very important note:" }}{PARA 0 "" 0 "" {TEXT -1 301 "(2 ) For large matrices the eigenvect command won't work well unless you enter at least one decimal number; if all entries are rational number s (expressed without decimal points), Maple tries to find the eigenval ues and eigenvectors algebraically and exactly, instead of numerically , and often fails. " }{TEXT 269 74 "Make sure at least one of your mat rix entries has a decimal point in it. " }}}{MARK "32 0" 0 } {VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }