{VERSION 4 0 "SUN SPARC SOLARIS" "4.0" } {USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 1 0 0 0 0 1 }{CSTYLE "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 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 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 } {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 Output" 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 10 "MATH 2280 " }}{PARA 258 "" 0 " " {TEXT 264 22 "PROJECT 2: EARTHQUAKES" }}{PARA 259 "" 0 "" {TEXT 270 10 "March 2001" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 240 "This is the Earthquake project on pages 333-334 of Edwar ds-Penney. You are mostly on your own for this project, but here is a small example of a spring system worked out on Maple, so that you can get an idea about useful commands to use. " }{TEXT 265 1 " " }} {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 204 "Let's st art with Example 1 on page 323 of Edwards-Penney. Initially it is an \+ unforced system with two masses and two springs, as you can see from t he description on page 323. We can write the system as " }{TEXT 257 7 "Mx''=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 v ector. 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 project" }}}{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 6 "x''=Ax" }{TEXT -1 176 ", and the eigenvectors of A determ ine fundamental modes, and the corresponding negative eigenvalues are \+ the (opposites) of the squares of the corresponding angular frequencie s:" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "eigenvects(A);" }} {PARA 11 "" 1 "" {XPPMATH 20 "6$7%!$+\"\"\"\"<#-%'vectorG6#7$F%!\"\"7% !#DF%<#-F(6#7$F%\"\"#" }}}{PARA 0 "" 0 "" {TEXT -1 286 "Therefore, the natural frequencies of this system are the 10 and 5, and the two fund amental modes correspond to the masses moving in opposite directions ( with equal amplitudes and angular frequency 10) and in parallel direc tions (with amplitude ratio of two and angular frequency 5). " }} {PARA 0 "" 0 "" {TEXT -1 164 " Now, let's consider the forced syst em with force vector equal to cos(wt)[0,50], i.e. the second mass is b eing forced periodically. In other words, the system " }{TEXT 262 11 "Mx''=Kx + F" }{TEXT -1 269 ", where F=cos(wt)[0,50] ; this is Example 3 on page 433. We follow the method described on that page to find \+ a particular solution to the forced oscillation problem, of the form g iven by equation (31). Here is the Maple version of the details summa rized in the text:" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 452 "F0:=e valm(inverse(M)&*vector([0,50]));\n #The F0 in the normalized equat ion (32), page 329\nIden:=array(1..2,1..2,identity);\n #the 2 by 2 \+ identity matrix\nAleft:=omega->evalm(A + omega^2*Iden);\n #the matr ix function multiplying\n #c on the left side of (32)\nc:=omega->ev alm(-inverse(Aleft(omega))&*F0);\n #the solution vector c(omega) to (32),\n #obtained by multiplying both sides of equation\n #(34) 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#>%&AleftGR6#%&omegaG6\"6$%)operatorG%&arrowGF(-%&evalmG 6#,&%\"AG\"\"\"*&)9$\"\"#F1%%IdenGF1F1F(F(F(" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"cGR6#%&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 329" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#-%'vectorG6#7$,$*&\"\"\"F),(\"%+DF)*& \"$D\"F))%&omegaG\"\"#F)!\"\"*$)F/\"\"%F)F)F1\"%]7,$*&,&!#vF)*$F.F)F)F )F*F1!#]" }}}{PARA 11 "" 1 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 636 "The vector c(w) above, times the oscillation cos(wt), is a par ticular solution to the forced oscillation problem we are considering. If we assume that our actual problem has a small amount of damping, \+ then we expect that this particular solution is very close to the stea dy periodic solution to the damped problem. See the discussion on pag e 330. We can study resonance phenomena for these slightly damped pro blems by plotting the maximum amplitude 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 maxi mum amplitude:" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 15 "norm(c(ome ga));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%$maxG6$,$*&\"\"\"F(-%$absG6 #,(\"%+DF(*&\"$D\"F()%&omegaG\"\"#F(!\"\"*$)F1\"\"%F(F(F3\"%]7,$-F*6#* &,&!#vF(*$F0F(F(F(F,F3\"#]" }}}{PARA 0 "" 0 "" {TEXT -1 229 "(If you w anted the actual amplitude of the solution when written in amplitude-p hase form you would want the square root of c1^2 + c2^2, but that func tion and the one above both blow up similarly near the two natural fre quencies.)" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 85 "plot(norm(c(om ega)),omega=0..15,amplitude=0..15,\n numpoints=200,color=`black` );" }}{PARA 13 "" 1 "" {GLPLOT2D 352 182 182 {PLOTDATA 2 "6&-%'CURVESG 6#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!)=IOKF?F ,7$$\"3x!He;G2=r#F,$\"3=*zjwU())o?F,7$$\"3?ze<&=Ezy#F,$\"31&\\uT<2c6#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['el# 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.t Y%F,7$$\"3v#oOt\")R6A%F,$\"3'*ybCx%[S%[F,7$$\"3mPw_0uN'H%F,$\"3&evhMB7 -I&F,7$$\"3E^-0gzhoVF,$\"3AiTaC&)>TeF,7$$\"3')Ge;LNtYWF,$\"3$*p5V%4Dae 'F,7$$\"3(3@U%QL=@XF,$\"3ezH$)[<>@vF,7$$\"3m=Qw_#Q&*f%F,$\"3#o'=*)Q3N# )))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,$\"3s 39[)\\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$$\"3F @T2rLj_\\F,$\"3%e>:BnyG4(F__l7$$\"3JrU&e4Bv&\\F,$\"34@$HA\"\\(H!zF__l7 $$\"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 183 "This is qualitatively the pi cture on page 330, figure 5.3.10. Notice the peaks at angular frequen cy 5 and 10, corresponding to approximate resonance near the two funda mental modes. " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 92 " We can get a plot of resonance as a function of per iod by recalling that 2*Pi/T=omega:" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 102 "res:=T->norm(c(2*Pi/T));\nplot(res(period),period=0. 1..3,amplitude=0..15,\nnumpoints=200,color=`black`);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$resGR6#%\"TG6\"6$%)operatorG%&arrowGF(-%%normG6#- %\"cG6#,$*&%#PiG\"\"\"9$!\"\"\"\"#F(F(F(" }}{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&ze D(=F*$\"3NRUe!>RDl%F-7$$\"3`lIh7w/;?F*$\"36'y[nQy_V&F-7$$\"3#yc8F&RWk@ 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$$\"3O8Fao% 3%\\FF*$\"3CWC()>p0l5F*7$$\"3UMoOt+!4!HF*$\"35AHf\"**R??\"F*7$$\"3HLmK 0.*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+>LFOv O9RF*$\"3Rv&z#GM!**[#F*7$$\"3I9HeE3-eSF*$\"3Vi'em$p.XFF*7$$\"3.f;LmP\" 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^*)ydDg ]#y%F*$\"39K\"ek$G*Qb%F*7$$\"3[V'Gds*\\F\\F*$\"3O*)y:eJ?&3&F*7$$\"3*49 GcS#Hx]F*$\"33G<2:A]PdF*7$$\"3f#['H>`$QA&F*$\"3k)e'pN2U;lF*7$$\"3*Q$oO BJil`F*$\"3MX/?(e\"3iuF*7$$\"39T\"GcIaI_&F*$\"3Gvn[@fLd))F*7$$\"3;AY# \\A8Xm&F*$\"3#oM1,C$\\j5!#<7$$\"3_w`2:3b:eF*$\"3ZF8F$zv'f8Fhu7$$\"3zmN rn?)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$$\"3qPx /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?#p O'\\e@g_Fhu7$$\"3NU([Z4I>Y'F*$\"3K0_1&=\\XH$Fhu7$$\"3i>T#[cpQ`'F*$\"34 w>=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+dHo F*$\"3td$*38**[k8Fhu7$$\"3&**)zfRx#z(pF*$\"3!GROH3V!y6Fhu7$$\"3CT\"Gc$ 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#4AW) 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*$\"3 r*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]+,A1W# *)*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$$\"3 uU&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.o9 QFhu7$$\"3c6BY;u=16Fhu$\"3g+dTMs`VVFhu7$$\"3M\\)pzT]/7\"Fhu$\"3+pC853_ F\\Fhu7$$\"3a!4='3q.N6Fhu$\"3wi()*fi(4mcFhu7$$\"3i9He&>r)\\6Fhu$\"3MkW v'RfQi'Fhu7$$\"3PE_/+n]j6Fhu$\"3SLCil!*3txFhu7$$\"3Zc7D1\"Fhu$\"3be[JW,J$=\"F\\x7$$\"3n$pQxl&*y?\"Fh u$\"3OC$p\\/\"f'f\"F\\x7$$\"3g;LmGOq97Fhu$\"3Fvijs!Rc(=F\\x7$$\"3vRze* f6:A\"Fhu$\"3x$*Q%QXaGE#F\\x7$$\"35>Qw%em$H7Fhu$\"3$y%pB1\"R)\\HF\\x7$ $\"3X)pR*p:AP7Fhu$\"338\">CK&p#>%F\\x7$$\"3_X7Fhu$\"3Y7``J 6F+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[7 Fhu$\"3([]$R)3#o[5F[z7$$\"3?lIOxki\\7Fhu$\"3e#Q[$*4pF=\"F[z7$$\"3fU5x4 &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 869 "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 here's one way to think about it. Maybe your instructor can help you more if it's still confusing. Anyhow, think of the grou nd as the zeroth story. In the rest frame it is shaking with oscillat ion Ecos(wt). And so its acceleration is its second time derivative, \+ namely -E*w^2*cos(wt). If you write down the inhomogeneous system of \+ EIGHT second order DE's for the accelerations of stories zero thru sev en, 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 sys tem " }{TEXT 266 1 "y" }{TEXT -1 78 "(t), then see what the shaking lo oks like to someone on the ground by letting " }}{PARA 0 "" 0 "" {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 other seven components will satisfy equation (2) on 334, exactly as the authors claim." }}{PARA 0 "" 0 " " {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 301 "(2) For large matrice s the eigenvect command won't work well unless you enter at least one \+ decimal number; if all entries are rational numbers (expressed without decimal points), Maple tries to find the eigenvalues and eigenvectors algebraically and exactly, instead of numerically, and often fails. \+ " }{TEXT 269 74 "Make sure at least one of your matrix entries has a d ecimal point in it. " }}}{MARK "30 3" 163 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }