{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 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 0 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 }{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 Outpu t" 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 1 }3 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }} {SECT 0 {PARA 256 "" 0 "" {TEXT 256 11 "MATH 2250-3" }}{PARA 258 "" 0 "" {TEXT 264 23 "SPRINGS AND EARTHQUAKES" }}{PARA 259 "" 0 "" {TEXT 269 17 "November 17, 2003" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 " " 0 "" {TEXT -1 353 "Your final Maple project for Math 2250 this semes ter is the Earthquake project on pages 437-438 of Edwards-Penney. The project is at our home page http://www.math.utah.edu/~korevaar/2250f all03/2250maple.html. In these notes we will work through the book ex amples from section 7.4, using illustrative Maple commands, as a warmu p for your project work." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 " " 0 "" {TEXT -1 245 "Let's start with example 1 on page 427 of Edwards -Penney, which we have also been working by hand. Initially it is an \+ unforced system with two masses and two springs, as you can see from t he description on page 427. 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" }}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 46 "with(linalg):with(plots):with(DEtools): #tools " }}}{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 a lso be written as " }{TEXT 261 6 "x''=Ax" }{TEXT -1 176 ", and the eig envectors of A determine fundamental modes, and the corresponding nega tive eigenvalues are the (opposites) of the squares of the correspondi ng angular frequencies:" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "e igenvects(A);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6$7%!$+\"\"\"\"<#-%'vec torG6#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 fundamental modes correspond to the masses moving in op posite directions (with equal amplitudes and angular frequency 10) an d in parallel directions (with amplitude ratio of two and angular freq uency 5). " }}{PARA 0 "" 0 "" {TEXT -1 164 " Now, let's consider \+ the forced system with force vector equal to cos(wt)[0,50], i.e. the s econd mass is being forced periodically. In other words, the system \+ " }{TEXT 262 11 "Mx''=Kx + F" }{TEXT -1 322 ", where F=cos(wt)[0,50] d iscussed on page 433. We follow the method described on that page to \+ find a particular solution to the forced oscillation problem, of the f orm given by equation (31). The details of this computation are expla ined in example 3 of the text and in our handwritten notes and here is a Maple version:" }}{PAGEBK }{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 325 "F0:=evalm(inverse(M)&*vector([0,50 ]));\n #The F0 in the normalized equation (30), page 433\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 on the left side \+ of (32)\nc:=omega->evalm(-inverse(Aleft(omega))&*F0);\n #the vecto r c(omega) in (32)" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#F0G-%'vectorG 6#7$\"\"!\"#]" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%IdenG-%&arrayG6&%) identityG;\"\"\"\"\"#F)7\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%&Aleft Gf*6#%&omegaG6\"6$%)operatorG%&arrowGF(-%&evalmG6#,&%\"AG\"\"\"*&)9$\" \"#F1%%IdenGF1F1F(F(F(" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"cGf*6#%& omegaG6\"6$%)operatorG%&arrowGF(-%&evalmG6#,$-%#&*G6$-%(inverseG6#-%&A leftG6#9$%#F0G!\"\"F(F(F(" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 37 "c(omega); #see equation (35) page 433" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'vectorG6#7$,$*&\"%]7\"\"\",(\"%+DF**&\"$D\"F*)%&omeg aG\"\"#F*!\"\"*$)F0\"\"%F*F*F2F*,$*(\"#]F*,&\"#vF2*$F/F*F*F*F+F2F2" }} }{PARA 11 "" 1 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 612 "The ve ctor c(w) above, times the oscillation cos(wt), is a particular soluti on 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 steady state solut ion to the damped problem. See the discussion on page 434. We can st udy resonance phenomena for these slightly damped problems by plotting the maximum amplitude of the steady state solutions to the undamped p roblems, much like you did in the Tacoma Narrows project. Use ``norm, infinity'' to measure this maximum amplitude:" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "norm(c(omega),infinity);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%$maxG6$,$*&\"#]\"\"\"-%$absG6#*&,&\"#v!\"\"*$)%&omega G\"\"#F)F)F),(\"%+DF)*&\"$D\"F)F2F)F0*$)F3\"\"%F)F)F0F)F),$*&\"%]7F)-F +6#F5F0F)" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 77 "plot(norm(c(om ega)),omega=0..15,y=0..15,\n numpoints=200,color=`black`);" }} {PARA 13 "" 1 "" {GLPLOT2D 282 153 153 {PLOTDATA 2 "6&-%'CURVESG6#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$$\"3J 07C[^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$$\"3v rV([#>Fe8F,$\"3C_tTIrO4;F,7$$\"32W)oP:&RH9F,$\"3-\\sH9P9A;F,7$$\"3?d9H 3'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$$\"3 x!He;G2=r#F,$\"3=*zjwU())o?F,7$$\"3?ze<&=Ezy#F,$\"31&\\uT<2c6#F,7$$\"3 XlJjEjMiGF,$\"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$pQ x*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$$\"3e BZ%*)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&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,$\"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$$\"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$\\918kU?\"F \\cl7$$\"3n)z%eKrju\\F,$\"3o$R<#[Q\")>8F\\cl7$$\"39u[(\\*>3x\\F,$\"3$Q AI0(p,g9F\\cl7$$\"3D2!*[-QLz\\F,$\"3@C@$f?'\\=;F\\cl7$$\"3PSJ+5ce\")\\ 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+@$ybo`#F\\cl 7$$\"3;P^rAFO:]F,$\"3_zECn#*>k@F\\cl7$$\"3Gq#H-`9w,&F,$\"3&oRnz9Ao)=F \\cl7$$\"3S.MuPj')>]F,$\"3%\\U)fo_Ks;F\\cl7$$\"3_OvDX\"=@-&F,$\"3i%z!e lM],:F\\cl7$$\"3kp;x_*pV-&F,$\"379\"z!z)\\AO\"F\\cl7$$\"3w-eGg)=vdB7uF__l7$$\"3% [8Fa$)R\"\\]F,$\"3?\\+w#*H#ys'F__l7$$\"3N>T2rWOe]F,$\"3M*p;Wm*pbcF__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,$\"30i#3^aTY 8$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,$\"3 5GF*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,$\"3 PE%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&44=x$y F,$\"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$$\"3uIhA& >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&*F67$$ \"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]k8n5F,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,$\"3eWxOr_`x 7F,7$$\"3urY$pB+[>*F,$\"33w`x.yDe8F,7$$\"3.&>Ry1JBF*F,$\"304b]d.yh9F,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\")*fj43(zK F,7$$\"33He;Lx!yz*F,$\"3!efIoUL$)R%F,7$$\"3]hD^_i9r)*F,$\"3[\"z\"=`VHR nF,7$$\"3s4AW87'=\"**F,$\"3PVJt!e$=C(*F,7$$\"3$z&=Puhd_**F,$\"3`@iH!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<4Vc$*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%3scmf;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$\"3sj$yG Aa)ppF__l7$$\"3Vz$QsW*)=g$G1\"F__l$\"3)Q^S![rvk;F,7$$\"3-S!3;p4 02\"F__l$\"3R[MEL/l8:F,7$$\"3y&=PC#)>\"y5F__l$\"3q*pIq$p*>R\"F,7$$\"3j d: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![%[]t5 F,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$$\"3nrV([^; $o6F__l$\"3w(R\"G@:IcvF67$$\"3\"F__l$\"3!3Goh;4q$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$\"3 ONzVs3^F67$$\"3n%**)z*H[QE]qmH\"F__l$\"3'*H(y1,0\\x%F67$$\"33C['zNVVI\"F__l$\"3w8t8ZDFtYF67$$ \"3]-050iR68F__l$\"3;$[6whINe%F67$$\"3!Rze@vKM]\"\\%F6 7$$\"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@2J6k8F__l $\"3E*3?lOKb+%F67$$\"3'3P\"F__l$\"3'zXzG8:;$RF67$$\"3w&=Pu!)[(z8 F__l$\"3p>kI'yH*fQF67$$\"3U4>Q^i(oQ\"F__l$\"3iPN'QS&z'z$F67$$\"39c7D!G XWR\"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+2k5Z2O MF67$$\"3%pQxa#RhR9F__l$\"3c2WuV)>IQ$F67$$\"3O3o)o\\h;!G$F67$$\"3N(\\**[m0EY\"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$$\"3QPv]ERq# \\\"F__l$\"3#z+niR,@/$F67$$\"#:F)$\"3))**************HF6-%+AXESLABELSG 6$Q&omega6\"Q\"yFf`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 "Curve 1" }}}} {PARA 0 "" 0 "" {TEXT -1 288 "This is the picture on page 434. It mea sures the maximum amount that one of the masses will be displaced from equilibrium, for the particular solution we derived. Notice the peak s are at angular frequency 5 and 10, reflecting the resonance which wi ll occur at the natural frequencies. " }}{PAGEBK }{PARA 0 "" 0 "" {TEXT -1 91 " We can get a plot of resonance as a function of peri od by recalling that 2*Pi/T=omega:" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 83 "res:=T->norm(c(2*Pi/T));\nplot(res(T),T=0.1..3,y=0..1 5,numpoints=200,color=`black`);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$ resGf*6#%\"TG6\"6$%)operatorG%&arrowGF(-%%normG6#-%\"cG6#,$*(\"\"#\"\" \"%#PiGF59$!\"\"F5F(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&zeD(=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()>p0l5 F*7$$\"3UMoOt+!4!HF*$\"35AHf\"**R??\"F*7$$\"3HLmK0.*o/$F*$\"38'G**R!GB X8F*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;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*49GcS#Hx]F*$\"33G<2:A]P dF*7$$\"3f#['H>`$QA&F*$\"3k)e'pN2U;lF*7$$\"3*Q$oOBJil`F*$\"3MX/?(e\"3i uF*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$$\"3zmNrn?)R)eF*$\"3k`20HXQi: Fhu7$$\"3=eBFh u7$$\"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?#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\"Fhu 7$$\"3pGd9uP.dnF*$\"3@sA,wnl-:Fhu7$$\"3\\yc8n+dHoF*$\"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*$\"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]+,A1W#*)*F*$\"3=!z`6BxL*=Fh u7$$\"3_18E3s&H+\"Fhu$\"3wtmP.\\^l?Fhu7$$\"3qW*)y,G6=5Fhu$\"3P&f9^%f,y AFhu7$$\"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$\"3 o!)3\\n[,:MFhu7$$\"3+$e;8ZM/4\"Fhu$\"3E\"*[Cq.o9QFhu7$$\"3c6BY;u=16Fhu $\"3g+dTMs`VVFhu7$$\"3M\\)pzT]/7\"Fhu$\"3+pC853_F\\Fhu7$$\"3a!4='3q.N6 Fhu$\"3wi()*fi(4mcFhu7$$\"3i9He&>r)\\6Fhu$\"3MkWv'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$$\"3vRze*f6:A\"Fhu$\"3x$*Q%QXaG E#F\\x7$$\"35>Qw%em$H7Fhu$\"3$y%pB1\"R)\\HF\\x7$$\"3X)pR*p:AP7Fhu$\"33 8\">CK&p#>%F\\x7$$\"3_X7Fhu$\"3Y7``J6F+sF\\x7$$\"38Y#*4=A3 Y7Fhu$\"3-Q>.#H#)\\\"yF\\x7$$\"3!4?:HGopC\"Fhu$\"3c4E=b,OU&)F\\x7$$\"3 nb6tZV&yC\"Fhu$\"3WKd(zm/lT*F\\x7$$\"3V5ra7/u[7Fhu$\"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$\"3PI z;0(=BY\"F[z7$$\"3'\\(\\*pg)R^7Fhu$\"3%39H&*4Dqe\"F[z7$$\"3+3`Skp$=D\" Fhu$\"3;(Q2:UsIt\"F[z7$$\"3$3k:=KvAD\"Fhu$\"3(R-:*QaZ3>F[z7$$\"3ltfAzO r_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$$\"3 R0t')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\"F hu$\"3-G]rhLnN$*F[z7$$\"3()o*\\\"oR(zD\"Fhu$\"3ql/5yNhyiF[z7$$\"3p,.cD BTe7Fhu$\"3i1[f\"*QGJZF[z7$$\"3^M1(Ho])e7Fhu$\"3EF%4F:\"y'z$F[z7$$\"3c n4QS!*Gf7Fhu$\"3xkEfvY@rJF[z7$$\"3h+8z(RF(f7Fhu$\"3;*G^Tn=Js#F[z7$$\"3 VL;?bd;g7Fhu$\"3_#)4??NN'Q#F[z7$$\"3Dm>h7Tgg7Fhu$\"3#*pv'[=8S7#F[z7$$ \"33*HA+ZU5E\"Fhu$\"3^E7.%)3)Q\">F[z7$$\"3!>jKu#3[h7Fhu$\"3]Tt\"oC#yT< F[z7$$\"3&\\'H%[=>>E\"Fhu$\"3a'y#)>;[#)f\"F[z7$$\"3*zH`AadBE\"Fhu$\"3M z*o,(>rw9F[z7$$\"3\"3jj'**ezi7Fhu$\"3Nj-*)[lZs8F[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&y 30UzojE\"Fhu$\"3InvI()F\\x7$$\"3tZXGW:Jn7Fhu$\"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$$\"3a18wW!oHF\"Fhu$\"395)z*3\">@ D&F\\x7$$\"3#e9z]/RnF\"Fhu$\"3.#\\)HYoz*G%F\\x7$$\"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$\"39 j-h$QM?>(Fhu7$$\"3TpQx#H%R69Fhu$\"3K_Uyh[cSmFhu7$$\"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%F hu7$$\"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(z wg^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$$\"3 v=PudSYt;Fhu$\"3zKh2%=#G^KFhu7$$\"3\")oPv\"p.$)o\"Fhu$\"3'erqgrtI=$Fhu 7$$\"3tY$p[t&4.JFhu7$$\"3%[$pQp6J=RKHrJQI\\>w\"Fhu$\"3%[V45#3l/HFhu7$$\"3[v],F$[lx\"Fhu$\" 3N^;'zty*eGFhu7$$\"3AU%)oQc!)*y\"Fhu$\"3c;&=s:8(>GFhu7$$\"3i8Fak)pb!=F hu$\"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>Fhu$\"34)y!3!op\"*[#Fhu7$$ \"3E,-/**)4,&>Fhu$\"3FKt$)\\\">hY#Fhu7$$\"3Lw_0p#4Y'>Fhu$\"3,UXENvBUCF hu7$$\"3(fAXq`)ez>Fhu$\"3O%4e8,0'=CFhu7$$\"3eg?TGGC%*>Fhu$\"3gzu3XgW'R #Fhu7$$\"3P&4>)3;U3?Fhu$\"3$RD8NndeP#Fhu7$$\"35E_/FZ;C?Fhu$\"3bpz(=s@R N#Fhu7$$\"33u[(*=1JQ?Fhu$\"3Cjv=VN*\\L#Fhu7$$\"3j\\**)zP9M0#Fhu$\"3;RE 'HeabJ#Fhu7$$\"3p(e<&G15n?Fhu$\"3#HK(G74f)H#Fhu7$$\"3d$yc`Jj?3#Fhu$\"3 !yXM-[:2G#Fhu7$$\"3K:IgP?9'4#Fhu$\"3lAY+\"p-XE#Fhu7$$\"3%*[(\\*e$e36#F hu$\"3IZSn%*=:[AFhu7$$\"3%R#['HDY_7#Fhu$\"3CF4BncsKAFhu7$$\"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'zL@#Fhu$\"3(H%R0!)e\")[@Fh u7$$\"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%)Q5 s:5#Fhu7$$\"3*QycBqgaG#Fhu$\"3G\\#y5S*R\"4#Fhu7$$\"3`U%)oU^*3I#Fhu$\"3 4k!)=NK>!3#Fhu7$$\"3)oV([f>'\\J#Fhu$\"3Bf0Ta=Gq?Fhu7$$\"3W%)oPY)[.L#Fh u$\"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'>Fh u7$$\"36tX\"f^)[/DFhu$\"3!)[y*z_8)f>Fhu7$$\"3&\\&4>+'p*=DFhu$\"3$GmYjM HG&>Fhu7$$\"3;@U%)[^4MDFhu$\"3Qlh;FlrX>Fhu7$$\"3Be;LJj-[DFhu$\"35P6fsV KR>Fhu7$$\"3)\\**)znPHiDFhu$\"3(y,IFhIH$>Fhu7$$\"3JBZ%HrY!yDFhu$\"33HE .2X/E>Fhu7$$\"3ahAX9(4Bf#Fhu$\"3!>?/OKh*>>Fhu7$$\"3_-050j*og#Fhu$\"3y1 6zDT)Q\">Fhu7$$\"3gE`1#\\IFhu7$$\"3OQw_'*fONEFhu$ \"3CYv=:WU->Fhu7$$\"3ooOt-5!*\\EFhu$\"3%p!4lZ)pn*=Fhu7$$\"3PQw_n`KkEFh u$\"3QzJ+&*QG\"*=Fhu7$$\"3k06Aa\\vzEFhu$\"3.5z*3n\\b)=Fhu7$$\"3%>Nqg*3 P$p#Fhu$\"3m0M)G(4g!)=Fhu7$$\"3m5@Um334FFhu$\"3hp8&z)p,v=Fhu7$$\"3%pQx M!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$$\"3 DD]+'3%H\"y#Fhu$\"35)[n?Ck4&=Fhu7$$\"3k&4>=aFfz#Fhu$\"3dB5K:7RY=Fhu7$$ \"327C[k$3-\"GFhu$\"3ZmIv^'>?%=Fhu7$$\"3wZ&4p#*)oDGFhu$\"3`#G7vozt$=Fh u7$$\"3yOtY%HFhu $\"3I,>LD\\`0=Fhu7$$\"3B.17l%Rf&HFhu$\"3<&*3N\"= " 0 "" {MPLTEXT 1 0 27 "E*(omega)^2*cos(om ega*t)*b;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#**%\"EG\"\"\")%&omegaG\" \"#F%-%$cosG6#*&F'F%%\"tGF%F%%\"bGF%" }}}{PARA 0 "" 0 "" {TEXT -1 795 "where b is the transpose of [1,1,1,1,1,1,1]. They ask, ``how can th e earthquake be forcing all seven stories, it seems like it's just sha king the bottom one.'' Well, the students are correct, but so is Edwa rds-Penney. The authors talk about an ``opposite inertial force'' bei ng the reason for this forcing term and here's one way to think about it: Think of the ground as the zeroth story. In the rest frame it is shaking with oscillation Ecos(wt). And so its acceleration is its se cond time derivative, namely -E*w^2*cos(wt). If you write down the in homogeneous system of EIGHT second order 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 solutio n 8-vector to this system " }{TEXT 265 1 "y" }{TEXT -1 78 "(t), then s ee what the shaking looks like to someone on the ground by letting " } }{PARA 0 "" 0 "" {TEXT 266 1 "x" }{TEXT -1 4 "(t)=" }{TEXT 267 1 "y" } {TEXT -1 213 "(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 compon ents will satisfy equation (2) on the bottom of page 303, exactly as t he authors claim." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 379 "(2) For large matrices the eigenvect command won't work well unless Maple knows you want a decimal approximation: If all ent ries are rational numbers (expressed without decimal points), Maple tr ies to find the eigenvalues and eigenvectors algebraically and exactly , instead of numerically, and often fails. The way around this is to \+ either ask for evalf(eigenvects(A)), or to " }{TEXT 268 74 "Make sure \+ at least one of your matrix entries has a decimal point in it. " }}} {MARK "13" 0 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }