(*********************************************************************** Mathematica-Compatible Notebook This notebook can be used on any computer system with Mathematica 3.0, MathReader 3.0, or any compatible application. The data for the notebook starts with the line of stars above. To get the notebook into a Mathematica-compatible application, do one of the following: * Save the data starting with the line of stars above into a file with a name ending in .nb, then open the file inside the application; * Copy the data starting with the line of stars above to the clipboard, then use the Paste menu command inside the application. Data for notebooks contains only printable 7-bit ASCII and can be sent directly in email or through ftp in text mode. Newlines can be CR, LF or CRLF (Unix, Macintosh or MS-DOS style). NOTE: If you modify the data for this notebook not in a Mathematica- compatible application, you must delete the line below containing the word CacheID, otherwise Mathematica-compatible applications may try to use invalid cache data. For more information on notebooks and Mathematica-compatible applications, contact Wolfram Research: web: http://www.wolfram.com email: info@wolfram.com phone: +1-217-398-0700 (U.S.) Notebook reader applications are available free of charge from Wolfram Research. ***********************************************************************) (*CacheID: 232*) (*NotebookFileLineBreakTest NotebookFileLineBreakTest*) (*NotebookOptionsPosition[ 119488, 3682]*) (*NotebookOutlinePosition[ 120385, 3711]*) (* CellTagsIndexPosition[ 120341, 3707]*) (*WindowFrame->Normal*) Notebook[{ Cell[CellGroupData[{ Cell["\<\ Internal Gravity Waves in an Atmosphere\ \>", "Title", TextAlignment->Center], Cell[TextData[{ "This notebook was used to find the trapped wave solution shown in Figure 1 \ of \"Anelastic Equations for the Atmosophere\". The equation we solve is:\n\n\ ", Cell[BoxData[ \(TraditionalForm\`d\^2\/dz\^2\)]], Cell[BoxData[ \(TraditionalForm\`\(w\&^\)\)]], " -", Cell[BoxData[ \(TraditionalForm\`1\/H\)]], Cell[BoxData[ \(TraditionalForm\`d\/dz\)]], Cell[BoxData[ \(TraditionalForm\`\(w\&^\)\)]], " -", Cell[BoxData[ \(TraditionalForm\`\(k\^2\)(1 - N\^2\/\[Omega]\^2)\)]], Cell[BoxData[ \(TraditionalForm\`\(w\&^\)\)]], "=0\n\nThis notebook could benefit from more comments and explanations. \ Basically, we are solving an ODE with boundary conditions of ", Cell[BoxData[ \(TraditionalForm\`\(w\&^\)\)]], "=0 at z=0 and z=lid , but solving it as an initial value problem. This \ means we keep \"guessing\" or improving the eigenvalue \[Omega] until the \ shooting scheme hits the target of ", Cell[BoxData[ \(TraditionalForm\`\(w\&^\)\)]], "=0 at z=lid." }], "Text"], Cell[CellGroupData[{ Cell["Functions for solving the ODE", "Section", InitializationCell->True], Cell[BoxData[ \(soln[\[Omega]_, k_] := NDSolve[{\n\t\t\t \(\(w'\)'\)[z] - \(w'\)[z]/H - k^2*\((1 - BuoyFreq^2/\[Omega]^2)\)*w[z] == 0, \n\t\t\t w[0] == 0, \n\t\t\t\(w'\)[0] == 1}, \n\t\tw, {z, 0, lid}]\)], "Input", InitializationCell->True], Cell[BoxData[ \(wAtlid[\[Omega]_, k_] := w[lid] /. \(soln[\[Omega], k]\)[\([1]\)]\)], "Input", InitializationCell->True] }, Closed]], Cell[CellGroupData[{ Cell["Some parameters we don't change much", "Section", InitializationCell->True], Cell[BoxData[ \(\(H = 10. ; \)\)], "Input", InitializationCell->True], Cell[BoxData[ \(\(lid = 20; \)\)], "Input", InitializationCell->True] }, Closed]], Cell[CellGroupData[{ Cell["Try it with constant BuoyFreq", "Section"], Cell[BoxData[ \(\(BuoyFreq = .01; \)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(trysoln = soln[ .006, 1]\)], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{"w", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.`, 20.`}}, "<>"]\), False, Editable->False]}], "}"}], "}"}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(ParametricPlot[{w[z] /. trysoln[\([1]\)], z}, {z, 0, lid}, AspectRatio -> 2]; \)\)], "Input"], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: 2 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Scaling calculations 0.472988 0.24704 0.047619 0.0952381 [ [.22595 .03512 -6 -9 ] [.22595 .03512 6 0 ] [.72003 .03512 -3 -9 ] [.72003 .03512 3 0 ] [.96707 .03512 -3 -9 ] [.96707 .03512 3 0 ] [.46049 .52381 -6 -4.5 ] [.46049 .52381 0 4.5 ] [.46049 1 -12 -4.5 ] [.46049 1 0 4.5 ] [.46049 1.47619 -12 -4.5 ] [.46049 1.47619 0 4.5 ] [.46049 1.95238 -12 -4.5 ] [.46049 1.95238 0 4.5 ] [ 0 0 0 0 ] [ 1 2 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .22595 .04762 m .22595 .05387 L s [(-1)] .22595 .03512 0 1 Mshowa .72003 .04762 m .72003 .05387 L s [(1)] .72003 .03512 0 1 Mshowa .96707 .04762 m .96707 .05387 L s [(2)] .96707 .03512 0 1 Mshowa .125 Mabswid .27536 .04762 m .27536 .05137 L s .32476 .04762 m .32476 .05137 L s .37417 .04762 m .37417 .05137 L s .42358 .04762 m .42358 .05137 L s .5224 .04762 m .5224 .05137 L s .5718 .04762 m .5718 .05137 L s .62121 .04762 m .62121 .05137 L s .67062 .04762 m .67062 .05137 L s .76944 .04762 m .76944 .05137 L s .81884 .04762 m .81884 .05137 L s .86825 .04762 m .86825 .05137 L s .91766 .04762 m .91766 .05137 L s .17654 .04762 m .17654 .05137 L s .12713 .04762 m .12713 .05137 L s .07772 .04762 m .07772 .05137 L s .02832 .04762 m .02832 .05137 L s .25 Mabswid 0 .04762 m 1 .04762 L s .47299 .52381 m .47924 .52381 L s [(5)] .46049 .52381 1 0 Mshowa .47299 1 m .47924 1 L s [(10)] .46049 1 1 0 Mshowa .47299 1.47619 m .47924 1.47619 L s [(15)] .46049 1.47619 1 0 Mshowa .47299 1.95238 m .47924 1.95238 L s [(20)] .46049 1.95238 1 0 Mshowa .125 Mabswid .47299 .14286 m .47674 .14286 L s .47299 .2381 m .47674 .2381 L s .47299 .33333 m .47674 .33333 L s .47299 .42857 m .47674 .42857 L s .47299 .61905 m .47674 .61905 L s .47299 .71429 m .47674 .71429 L s .47299 .80952 m .47674 .80952 L s .47299 .90476 m .47674 .90476 L s .47299 1.09524 m .47674 1.09524 L s .47299 1.19048 m .47674 1.19048 L s .47299 1.28571 m .47674 1.28571 L s .47299 1.38095 m .47674 1.38095 L s .47299 1.57143 m .47674 1.57143 L s .47299 1.66667 m .47674 1.66667 L s .47299 1.7619 m .47674 1.7619 L s .47299 1.85714 m .47674 1.85714 L s .25 Mabswid .47299 0 m .47299 2 L s 0 0 m 1 0 L 1 2 L 0 2 L closepath clip newpath .5 Mabswid .47299 .04762 m .56798 .08523 L .60926 .10439 L .64338 .12489 L .65049 .13049 L .65624 .13579 L .66091 .14095 L .66436 .14571 L .66687 .15023 L .66785 .1525 L .66869 .15501 L .6693 .15752 L .66954 .15893 L .66969 .16022 L .66979 .16273 L .66967 .16511 L .66949 .16654 L .66925 .16788 L .66861 .17039 L .66765 .1731 L .6663 .17603 L .66302 .18131 L .65907 .18617 L .64665 .19735 L .62875 .20916 L .59083 .22812 L .54016 .24844 L .43521 .2851 L .32916 .32473 L .3058 .33558 L .28495 .34711 L .27083 .35696 L .26475 .36225 L .25955 .36783 L .25601 .37275 L .25446 .37555 L .25331 .37814 L .2525 .38054 L .25198 .38272 L .25164 .38523 L .25156 .38761 L .25176 .39019 L .25197 .39149 L .25229 .39293 L .25307 .39546 L .25405 .39784 L .25716 .40318 L .26186 .40893 L .27384 .41927 L .28922 .42909 L Mistroke .32765 .44752 L .44195 .48721 L .56164 .52401 L .624 .5452 L .67185 .56513 L .6884 .57384 L .70253 .58303 L .70847 .58784 L .71379 .59311 L .71752 .5979 L .71882 .60003 L .72001 .60237 L .721 .6049 L .72141 .60629 L .72171 .60759 L .72204 .60991 L .72211 .61241 L .72202 .61374 L .72182 .61519 L .72151 .6167 L .72111 .61811 L .72017 .62063 L .71881 .62333 L .71507 .62868 L .71028 .63372 L .69808 .64315 L .68211 .65244 L .66413 .66099 L .61324 .68029 L .47374 .72109 L .33672 .7592 L .27911 .77799 L .23582 .79581 L .21776 .80564 L .2101 .81083 L .20343 .81634 L .19904 .82097 L .19717 .82342 L .19551 .82607 L .19433 .82849 L .19352 .83072 L .193 .83284 L .19272 .83508 L .19273 .83754 L .19288 .83893 L .19312 .84021 L .19378 .84257 L .19482 .84506 L .19604 .84729 L .19743 .84938 L Mistroke .20148 .85414 L .20723 .85927 L .21371 .86396 L .23248 .87451 L .27926 .89334 L .33715 .91123 L .64068 .98766 L .70734 1.00794 L .73394 1.01789 L .75417 1.02703 L .77201 1.03761 L .77886 1.04317 L .78165 1.04601 L .78411 1.04904 L .78577 1.05162 L .78692 1.05396 L .78774 1.05632 L .78821 1.05855 L .78836 1.06115 L .78813 1.06355 L .78747 1.06619 L .78697 1.06752 L .78629 1.06899 L .78322 1.07369 L .78088 1.07636 L .77831 1.07884 L .76354 1.0892 L .74568 1.09815 L .72169 1.10771 L .64988 1.12958 L .5686 1.14944 L .37664 1.19107 L .29617 1.20935 L .2195 1.22944 L .18988 1.23876 L .16315 1.24877 L .14325 1.25821 L .12994 1.26684 L .12706 1.26928 L .12443 1.27188 L .12252 1.27413 L .12084 1.27656 L .12001 1.27805 L .11938 1.27942 L .1189 1.28073 L .11852 1.28211 L .11818 1.28451 L .11826 1.28714 L .1188 1.28967 L .11967 1.29196 L Mistroke .12031 1.29323 L .12113 1.29461 L .12296 1.29712 L .12773 1.30195 L .13356 1.30641 L .15005 1.31572 L .17412 1.32574 L .23289 1.34385 L .42146 1.38509 L .62634 1.42449 L .72307 1.44531 L .76851 1.45674 L .80489 1.46751 L .83316 1.47788 L .84513 1.48334 L .8555 1.4891 L .86235 1.49396 L .86526 1.49653 L .86783 1.49929 L .86967 1.50183 L .87092 1.50415 L .87176 1.50647 L .87216 1.50868 L .87219 1.51087 L .87178 1.51329 L .87086 1.51581 L .86955 1.51819 L .86793 1.52039 L .8661 1.52242 L .86082 1.52698 L .85381 1.53162 L .84426 1.53669 L .81921 1.54691 L .78557 1.55747 L .74715 1.56744 L .66036 1.58611 L .42266 1.628 L .21737 1.66462 L .12389 1.68543 L .08782 1.69561 L .06136 1.70491 L .05052 1.70961 L .04074 1.71473 L .0335 1.71957 L .02855 1.72406 L .02642 1.72676 L .02498 1.72929 L .02417 1.73157 L .02381 1.73403 L .02383 1.73531 L Mistroke .02402 1.73672 L .02441 1.73819 L .02494 1.73956 L .0263 1.74202 L .02836 1.74464 L .03123 1.74746 L .03451 1.75007 L .04411 1.75596 L .05504 1.76115 L .06914 1.76667 L .1408 1.78672 L .23323 1.8055 L .35324 1.82573 L .59424 1.8623 L .83079 1.90245 L .87526 1.9122 L .91491 1.92262 L .93368 1.92861 L .94802 1.93404 L .95904 1.93914 L .96823 1.94466 L .97079 1.94662 L .97309 1.94868 L .97619 1.95238 L Mfstroke % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{144, 288}, ImageMargins->{{43, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHg7oo0013 Ool00`00Oomoo`0BOol2000fOol004=oo`03001oogoo01Aoo`03001oogoo03=oo`00@goo00<007oo Ool05Goo00<007ooOol0goo0013Ool2000; Ool3000mOol004=oo`03001oogoo00Qoo`80041oo`00@goo00<007ooOol01Woo0P00@Woo0013Ool0 0`00Oomoo`03Ool30014Ool004=oo`03001oogoo00<004Moo`00@goo0`00BWoo0010Ool4001Woo0`001Woo00<007ooOol0BWoo000gOol30009Ool00`00 Oomoo`1:Ool003Aoo`<000aoo`03001oogoo04Yoo`00Goo10001Woo 00<007ooOol0BWoo000fOol3000:Ool00`00Oomoo`1:Ool0039oo`@000eoo`03001oogoo04Yoo`00 ;goo0`004Goo00<007ooOol0BWoo000/Ool3000DOol00`00Oomoo`1:Ool002Uoo`<001Moo`03001o ogoo04Yoo`009goo0P006Woo00<007ooOol0BWoo000TOol3000LOol00`00Oomoo`1:Ool0025oo`<0 01moo`03001oogoo04Yoo`0087oo00<007ooOol087oo00<007ooOol0BWoo000NOol2000SOol2001; Ool001eoo`03001oogoo02=oo`03001oogoo04Yoo`0077oo00<007ooOol097oo00<007ooOol0BWoo 000LOol00`00Oomoo`0TOol00`00Oomoo`1:Ool001aoo`03001oogoo02Aoo`03001oogoo04Yoo`00 77oo00<007ooOol097oo00<007ooOol0BWoo000LOol2000UOol00`00Oomoo`1:Ool001ioo`03001o ogoo029oo`03001oogoo04Yoo`007Woo0`008Woo00<007ooOol0BWoo000QOol2000POol00`00Oomo o`1:Ool002=oo`<001eoo`03001oogoo04Yoo`009Woo0P006goo00<007ooOol0BWoo000XOol3000H Ool00`00Oomoo`1:Ool002]oo`@001Aoo`03001oogoo04Yoo`00;goo100047oo0P00Bgoo000cOol4 000goo100017oo00<007ooOol0 BWoo000oOol5001Ool004=oo`03001oogoo03aoo`@000Yoo`00@goo 00<007ooOol0@7oo0P0027oo000fOol40003Ool20004Ool00`00Oomoo`12Ool00`00Oomoo`05Ool0 03Moo`03001oogoo009oo`04001oogoo0003Ool00`00Oomoo`13Ool20005Ool003Qoo`05001oogoo Ool00002Ool01000Oomoogoo0P00AGoo00<007ooOol00goo000iOol01000Oomoo`000Woo00D007oo Oomoo`0004aoo`00=Woo00@007ooOol0009oo`04001oogoo0003Ool00`00Oomoo`1:Ool003Moo`80 00Aoo`8000Aoo`03001oogoo04Yoo`00@goo00<007ooOol0BWoo0013Ool00`00Oomoo`1:Ool004=o o`03001oogoo04Yoo`00\ \>"], ImageRangeCache->{{{0, 143}, {287, 0}} -> {-1.93632, -0.799267, 0.0286107, 0.0742138}}] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["Try shooting for w=0 at the top", "Section"], Cell["\<\ The roots in the following plot \"hit the target\" of w=0 at \ z=lid:\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(\(Plot[wAtlid[\[Omega], 1], {\[Omega], .006, .007}]; \)\)], "Input"], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Scaling calculations -5.69048 952.381 0.287956 0.118317 [ [.21429 .27546 -18 -9 ] [.21429 .27546 18 0 ] [.40476 .27546 -18 -9 ] [.40476 .27546 18 0 ] [.59524 .27546 -18 -9 ] [.59524 .27546 18 0 ] [.78571 .27546 -18 -9 ] [.78571 .27546 18 0 ] [.97619 .27546 -15 -9 ] [.97619 .27546 15 0 ] [.01131 .05132 -12 -4.5 ] [.01131 .05132 0 4.5 ] [.01131 .16964 -12 -4.5 ] [.01131 .16964 0 4.5 ] [.01131 .40627 -6 -4.5 ] [.01131 .40627 0 4.5 ] [.01131 .52459 -6 -4.5 ] [.01131 .52459 0 4.5 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .21429 .28796 m .21429 .29421 L s [(0.0062)] .21429 .27546 0 1 Mshowa .40476 .28796 m .40476 .29421 L s [(0.0064)] .40476 .27546 0 1 Mshowa .59524 .28796 m .59524 .29421 L s [(0.0066)] .59524 .27546 0 1 Mshowa .78571 .28796 m .78571 .29421 L s [(0.0068)] .78571 .27546 0 1 Mshowa .97619 .28796 m .97619 .29421 L s [(0.007)] .97619 .27546 0 1 Mshowa .125 Mabswid .07143 .28796 m .07143 .29171 L s .11905 .28796 m .11905 .29171 L s .16667 .28796 m .16667 .29171 L s .2619 .28796 m .2619 .29171 L s .30952 .28796 m .30952 .29171 L s .35714 .28796 m .35714 .29171 L s .45238 .28796 m .45238 .29171 L s .5 .28796 m .5 .29171 L s .54762 .28796 m .54762 .29171 L s .64286 .28796 m .64286 .29171 L s .69048 .28796 m .69048 .29171 L s .7381 .28796 m .7381 .29171 L s .83333 .28796 m .83333 .29171 L s .88095 .28796 m .88095 .29171 L s .92857 .28796 m .92857 .29171 L s .25 Mabswid 0 .28796 m 1 .28796 L s .02381 .05132 m .03006 .05132 L s [(-2)] .01131 .05132 1 0 Mshowa .02381 .16964 m .03006 .16964 L s [(-1)] .01131 .16964 1 0 Mshowa .02381 .40627 m .03006 .40627 L s [(1)] .01131 .40627 1 0 Mshowa .02381 .52459 m .03006 .52459 L s [(2)] .01131 .52459 1 0 Mshowa .125 Mabswid .02381 .07498 m .02756 .07498 L s .02381 .09865 m .02756 .09865 L s .02381 .12231 m .02756 .12231 L s .02381 .14598 m .02756 .14598 L s .02381 .1933 m .02756 .1933 L s .02381 .21697 m .02756 .21697 L s .02381 .24063 m .02756 .24063 L s .02381 .26429 m .02756 .26429 L s .02381 .31162 m .02756 .31162 L s .02381 .33528 m .02756 .33528 L s .02381 .35895 m .02756 .35895 L s .02381 .38261 m .02756 .38261 L s .02381 .42994 m .02756 .42994 L s .02381 .4536 m .02756 .4536 L s .02381 .47726 m .02756 .47726 L s .02381 .50093 m .02756 .50093 L s .02381 .02766 m .02756 .02766 L s .02381 .00399 m .02756 .00399 L s .02381 .54825 m .02756 .54825 L s .02381 .57192 m .02756 .57192 L s .02381 .59558 m .02756 .59558 L s .25 Mabswid .02381 0 m .02381 .61803 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath .5 Mabswid .02381 .52896 m .02499 .52891 L .02605 .52885 L .02846 .52867 L .03068 .52844 L .03279 .52816 L .0375 .52732 L .04262 .5261 L .05293 .52262 L .06244 .51824 L .0829 .50514 L .10458 .48603 L .14241 .44129 L .1827 .38104 L .22147 .31522 L .2627 .24225 L .30241 .17444 L .3406 .1162 L .36192 .08814 L .38125 .06608 L .40163 .04673 L .42038 .03275 L .42955 .02732 L .43962 .02246 L .44975 .01875 L .45528 .01722 L .46044 .01612 L .46287 .01572 L .46545 .01536 L .46679 .0152 L .46828 .01506 L .46964 .01495 L .47089 .01486 L .47209 .0148 L .47335 .01475 L .47462 .01472 L .47533 .01472 L .476 .01472 L .47721 .01473 L .47832 .01475 L .47959 .0148 L .48079 .01486 L .48304 .01502 L .48543 .01525 L .4897 .01584 L .49429 .01671 L .49935 .01796 L .50951 .02137 L .5188 .02552 L .53969 .03841 L Mistroke .55871 .05425 L .57945 .07567 L .62015 .12881 L .65933 .191 L .70096 .26461 L .74108 .33825 L .77969 .40729 L .82074 .47427 L .86028 .52859 L .89831 .56836 L .91785 .58331 L .92788 .58942 L .93878 .59482 L .94877 .59863 L .95316 .59995 L .958 .60116 L .96256 .60205 L .96683 .60267 L .96918 .60293 L .97134 .60311 L .97247 .60318 L .9737 .60324 L .97499 .60329 L .97619 .60332 L Mfstroke % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{35, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgOol000moo`03001oogoo06aoo`<0015oo`8008ao o`003goo0P00Jgoo0P005Woo0P00RWoo000?Ool00`00Oomoo`1YOol00`00Oomoo`0HOol20028Ool0 00moo`03001oogoo06Moo`8001eoo`8008Ioo`003goo00<007ooOol0IGoo0P008Goo0P00Q7oo000? Ool00`00Oomoo`1TOol00`00Oomoo`0SOol00`00Oomoo`21Ool000Moo`@000Aoo`03001oogoo06=o o`03001oogoo02Eoo`03001oogoo081oo`0027oo00<007ooOol017oo00<007ooOol0HWoo00<007oo Ool09goo00<007ooOol0Ogoo00000goo000000030003Ool00`00Oomoo`03Ool3001QOol00`00Oomo o`0YOol00`00Oomoo`1nOol000Yoo`03001oogoo009oo`03001oogoo061oo`03001oogoo02]oo`03 001oogoo07eoo`001goo00@007ooOol000Aoo`03001oogoo05moo`03001oogoo02eoo`03001oogoo 07aoo`0027oo0P001Goo00<007ooOol0GWoo00<007ooOol0;goo00<007ooOol0Ngoo000?Ool00`00 Oomoo`1MOol00`00Oomoo`0aOol00`00Oomoo`1jOol000moo`03001oogoo05eoo`03001oogoo039o o`03001oogoo07Uoo`003goo0P00GGoo00<007ooOol0=7oo00<007ooOol0N7oo000?Ool00`00Oomo o`1KOol00`00Oomoo`0eOol00`00Oomoo`1hOol000moo`03001oogoo05Yoo`03001oogoo03Moo`03 001oogoo07Moo`003goo00<007ooOol0FGoo00<007ooOol0>Goo00<007ooOol0MWoo000?Ool00`00 Oomoo`1HOol00`00Oomoo`0kOol00`00Oomoo`1eOol000moo`03001oogoo05Qoo`03001oogoo03ao o`03001oogoo07Aoo`003goo0P00F7oo00<007ooOol0?Goo00<007ooOol0M7oo000?Ool00`00Oomo o`1FOol00`00Oomoo`0oOol00`00Oomoo`1cOol000moo`03001oogoo05Eoo`03001oogoo045oo`03 001oogoo079oo`003goo00<007ooOol0E7oo00<007ooOol0@goo00<007ooOol0LGoo000?Ool00`00 Oomoo`1COol00`00Oomoo`14Ool00`00Oomoo`1aOol000moo`03001oogoo059oo`03001oogoo04Io o`03001oogoo071oo`003goo00<007ooOol0DWoo00<007ooOol0Agoo00<007ooOol0Kgoo000?Ool2 001BOol00`00Oomoo`19Ool00`00Oomoo`1^Ool000moo`03001oogoo051oo`03001oogoo04]oo`03 001oogoo06eoo`003goo00<007ooOol0D7oo00<007ooOol0Bgoo00<007ooOol0KGoo000?Ool00`00 Oomoo`1?Ool00`00Oomoo`1=Ool00`00Oomoo`1/Ool000moo`03001oogoo04moo`03001oogoo04eo o`03001oogoo06aoo`003goo00<007ooOol0CWoo00<007ooOol0Cgoo00<007ooOol0Jgoo000?Ool2 001>Ool00`00Oomoo`1AOol00`00Oomoo`1ZOol000moo`03001oogoo04eoo`03001oogoo055oo`03 001oogoo06Yoo`003goo00<007ooOol0C7oo00<007ooOol0Dgoo00<007ooOol0JGoo000?Ool00`00 Oomoo`1;Ool00`00Oomoo`1DOol00`00Oomoo`1YOol000Moo`<000Eoo`03001oogoo04]oo`03001o ogoo05Eoo`03001oogoo06Qoo`0027oo00<007ooOol017oo00<007ooOol0BWoo00<007ooOol0Egoo 00<007ooOol0Igoo00000goo000000030002Ool00`00Oomoo`04Ool3001:Ool00`00Oomoo`1GOol0 0`00Oomoo`1WOol000Qoo`03001oogoo00Aoo`03001oogoo04Uoo`03001oogoo05Uoo`03001oogoo 06Ioo`0027oo00<007ooOol017oo00<007ooOol0B7oo00<007ooOol0FWoo00<007ooOol0IWoo0007 Ool20006Ool00`00Oomoo`18Ool00`00Oomoo`1KOol00`00Oomoo`1UOol000moo`03001oogoo04Mo o`03001oogoo05eoo`03001oogoo06Aoo`003goo00<007ooOol0AWoo00<007ooOol0GWoo00<007oo Ool0I7oo000?Ool00`00Oomoo`16Ool00`00Oomoo`1OOol00`00Oomoo`1SOol000moo`8004Ioo`03 001oogoo061oo`03001oogoo06=oo`003goo00<007ooOol0AGoo00<007ooOol0HGoo00<007ooOol0 HWoo000?Ool00`00Oomoo`14Ool00`00Oomoo`1ROol00`00Oomoo`1ROol000moo`03001oogoo04=o o`03001oogoo06Aoo`03001oogoo065oo`003goo00<007ooOol0@goo00<007ooOol0I7oo00<007oo Ool0HGoo000?Ool00`00Oomoo`12Ool00`00Oomoo`1VOol00`00Oomoo`1POol000moo`8004=oo`03 001oogoo06Moo`03001oogoo05moo`003goo00<007ooOol0@Goo00<007ooOol0J7oo00<007ooOol0 Ggoo000?Ool00`00Oomoo`10Ool00`00Oomoo`1ZOol00`00Oomoo`1NOol000moo`03001oogoo041o o`03001oogoo06Yoo`03001oogoo05ioo`003goo00<007ooOol0?goo00<007ooOol0K7oo00<007oo Ool0GGoo000?Ool00`00Oomoo`0oOol00`00Oomoo`1/Ool00`00Oomoo`1MOol000moo`03001oogoo 03ioo`03001oogoo06ioo`03001oogoo05aoo`003goo0P00?Woo00<007ooOol0Kgoo00<007ooOol0 G7oo000?Ool00`00Oomoo`0mOol00`00Oomoo`1`Ool00`00Oomoo`1KOol000moo`03001oogoo021o o`8000Aoo`8000Aoo`8000Aoo`8000Aoo`80009oo`D001=oo`8000Aoo`8000Aoo`8000Aoo`8000Ao o`8000Aoo`<0019oo`8000Aoo`8000Aoo`8000Aoo`8000Aoo`8000Aoo`8000Yoo`03001oogoo00Io o`8000Aoo`8000Aoo`8000Aoo`8000Aoo`8000Aoo`8001Moo`8000Aoo`8000Aoo`8000Aoo`8000Ao o`03001oogoo009oo`003goo00<007ooOol07goo00@007ooOol000=oo`8000=oo`04001oogoo0002 Ool01000Oomoo`000Woo00P007ooOol007oo001oo`0057oo00@007ooOol000=oo`8000=oo`04001o ogoo0002Ool01000Oomoo`000Woo00@007ooOol000Aoo`03001oogoo011oo`04001oogoo0003Ool2 0003Ool01000Oomoo`000Woo00@007ooOol0009oo`04001oogoo0002Ool01000Oomoo`002Goo00<0 07ooOol01Goo00@007ooOol000=oo`8000=oo`04001oogoo0002Ool01000Oomoo`000Woo00@007oo Ool0009oo`04001oogoo000EOol01000Oomoo`000goo0P000goo00@007ooOol0009oo`04001oogoo 0003Ool00`00Oomoo`02Ool000moo`03001oogoo01moo`04001oogoo0008Ool01000Oomoo`000Woo 00@007ooOol0009oo`03001oogoo008000=oo`03001oogoo015oo`04001oogoo0008Ool01000Oomo o`000Woo00@007ooOol0009oo`04001oogoo0002Ool4000AOol01000Oomoo`0027oo00@007ooOol0 009oo`04001oogoo0002Ool01000Oomoo`000Woo00@007ooOol000Yoo`03001oogoo00Aoo`04001o ogoo0008Ool01000Oomoo`000Woo00@007ooOol0009oo`04001oogoo0002Ool01000Oomoo`005Goo 00@007ooOol000Qoo`04001oogoo0002Ool01000Oomoo`0017oo00<007ooOol00Goo000?Ool00`00 Oomoo`0OOol01000Oomoo`0027oo00@007ooOol0009oo`04001oogoo0002Ool300000goo001oo`03 Ool00`00Oomoo`0@Ool01000Oomoo`0027oo00@007ooOol0009oo`04001oogoo0002Ool30003Ool0 0`00Ool0000BOol01000Oomoo`0027oo00@007ooOol0009oo`04001oogoo0002Ool30003Ool3000; Ool00`00Oomoo`04Ool01000Oomoo`0027oo00@007ooOol0009oo`04001oogoo0002Ool30004Ool2 000FOol01000Oomoo`0027oo00@007ooOol0009oo`04001oogoo0004Ool00`00Oomoo`01Ool000mo o`80021oo`04001oogoo0008Ool01000Oomoo`000Woo00@007ooOol0009oo`04001oogoo0002Ool0 1000Oomoo`004Woo00@007ooOol000Qoo`04001oogoo0002Ool01000Oomoo`000Woo00<007ooOol0 0goo00<007oo00004Woo00@007ooOol000Qoo`04001oogoo0002Ool01000Oomoo`000Woo00<007oo Ool00goo00<007ooOol037oo00<007ooOol00goo00@007ooOol000Qoo`04001oogoo0002Ool01000 Oomoo`000Woo00<007ooOol00goo00@007ooOol001Eoo`04001oogoo0008Ool01000Oomoo`000Woo 00@007ooOol000Eoo`40005oo`5oo`003goo00<007ooOol087oo0P002Woo0P0017oo0P0017oo0`00 0goo0P0057oo0P002Woo0P0017oo0P0017oo0`000goo0P004goo0P002Woo0P0017oo0P0017oo0`00 0goo0`002goo00<007ooOol017oo0P002Woo0P0017oo0P0017oo0`000goo0P005goo0P002Woo0P00 17oo0P000goo10000Woo000?Ool00`00Oomoo`0iOol00`00Oomoo`1hOol00`00Oomoo`1GOol000mo o`03001oogoo03Qoo`03001oogoo07Uoo`03001oogoo05Moo`003goo00<007ooOol0>7oo00<007oo Ool0NWoo00<007ooOol0EWoo000?Ool00`00Oomoo`0gOol00`00Oomoo`1kOol00`00Oomoo`1FOol0 00Qoool0010000Uoo`003goo00<007ooOol02Woo00<007ooOol02Woo00<007ooOol02Goo00<007oo Ool02Woo00<007ooOol00goo00<007ooOol017oo00<007ooOol02Woo00<007ooOol02Woo00<007oo Ool02Woo00<007ooOol02Goo00<007ooOol02Woo00<007ooOol02Woo00<007ooOol02Woo00<007oo Ool02Woo00<007ooOol02Woo00<007ooOol00goo00<007ooOol00goo00<007ooOol02Woo00<007oo Ool02Woo00<007ooOol02Woo00<007ooOol02Woo00<007ooOol02Woo00<007ooOol03Goo000?Ool0 0`00Oomoo`0`Ool00`00Oomoo`03Ool00`00Oomoo`0[Ool00`00Oomoo`0`Ool00`00Oomoo`0MOol0 0`00Oomoo`0@Ool00`00Oomoo`0aOol00`00Oomoo`0=Ool000moo`03001oogoo03Eoo`03001oogoo 081oo`03001oogoo05=oo`003goo00<007ooOol0=7oo00<007ooOol0PGoo00<007ooOol0Dgoo000? Ool00`00Oomoo`0dOol00`00Oomoo`22Ool00`00Oomoo`1BOol000moo`03001oogoo03=oo`03001o ogoo08=oo`03001oogoo059oo`003goo0P00=7oo00<007ooOol0Q7oo00<007ooOol0DGoo000?Ool0 0`00Oomoo`0bOol00`00Oomoo`25Ool00`00Oomoo`1AOol000moo`03001oogoo039oo`03001oogoo 08Ioo`03001oogoo051oo`003goo00<007ooOol0Woo000? Ool00`00Oomoo`0IOol00`00Oomoo`2eOol00`00Oomoo`0jOol000moo`03001oogoo01Qoo`03001o ogoo0;Moo`03001oogoo03Uoo`003goo00<007ooOol067oo00<007ooOol0]goo00<007ooOol0>Goo 000?Ool00`00Oomoo`0GOol00`00Oomoo`2iOol00`00Oomoo`0hOol000moo`8001Moo`03001oogoo 0;]oo`03001oogoo03Moo`003goo00<007ooOol05Goo00<007ooOol0_7oo00<007ooOol0=goo000? Ool00`00Oomoo`0DOol00`00Oomoo`2nOol00`00Oomoo`0fOol000moo`03001oogoo019oo`800<9o o`03001oogoo03Eoo`003goo00<007ooOol04Goo00<007ooOol0`goo00<007ooOol0=7oo000?Ool0 0`00Oomoo`0@Ool00`00Oomoo`34Ool00`00Oomoo`0dOol000moo`03001oogoo00moo`03001oogoo 01oo`03001oogoo02]oo`003goo00<007oo Ool0hGoo00<007ooOol0:Woo000?Ool00`00Oomoo`3ROol00`00Oomoo`0YOol000moo`800>Aoo`03 001oogoo02Qoo`003goo00<007ooOol0i7oo00<007ooOol09goo000?Ool00`00Oomoo`3UOol00`00 Oomoo`0VOol000moo`03001oogoo0>Ioo`03001oogoo02Eoo`003goo00<007ooOol0igoo00<007oo Ool097oo000?Ool00`00Oomoo`3XOol00`00Oomoo`0SOol000moo`03001oogoo0>Uoo`8002=oo`00 3goo0P00k7oo00<007ooOol087oo000?Ool00`00Oomoo`3/Ool00`00Oomoo`0OOol000moo`03001o ogoo0>eoo`03001oogoo01ioo`003goo00<007ooOol0kWoo00<007ooOol07Goo000?Ool00`00Oomo o`3_Ool00`00Oomoo`0LOol000moo`03001oogoo0?1oo`<001]oo`003goo0P00m7oo0`0067oo000? Ool00`00Oomoo`3fOol3000EOol000moo`03001oogoo0?Uoo`H000moo`003goo00<007ooOol0ogoo 3goo000?Ool00`00Oomoo`3oOol?Ool000moo`03001oogoo0?moo`moo`003goo00<007ooOol0ogoo 3goo003oOolQOol00?moob5oo`00ogoo8Goo003oOolQOol00?moob5oo`00\ \>"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {0.00594011, -2.59691, 3.89667*^-6, 0.0313659}}] }, Open ]], Cell[TextData["Let's find one eigenvalue \[Omega] exactly:"], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(aroot = FindRoot[wAtlid[\[Omega], 1], {\[Omega], .006, .0064}]\)], "Input"], Cell[BoxData[ \({\[Omega] \[Rule] 0.00622375355877907399`}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(hitsoln = soln[\[Omega] /. aroot[\([1]\)], 1]\)], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{"w", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.`, 20.`}}, "<>"]\), False, Editable->False]}], "}"}], "}"}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(ParametricPlot[{w[z] /. hitsoln[\([1]\)], z}, {z, 0, lid}, AspectRatio -> 2]; \)\)], "Input"], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: 2 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Scaling calculations 0.529717 0.248775 0.047619 0.0952381 [ [.03217 .03512 -6 -9 ] [.03217 .03512 6 0 ] [.28094 .03512 -6 -9 ] [.28094 .03512 6 0 ] [.77849 .03512 -3 -9 ] [.77849 .03512 3 0 ] [.51722 .52381 -6 -4.5 ] [.51722 .52381 0 4.5 ] [.51722 1 -12 -4.5 ] [.51722 1 0 4.5 ] [.51722 1.47619 -12 -4.5 ] [.51722 1.47619 0 4.5 ] [.51722 1.95238 -12 -4.5 ] [.51722 1.95238 0 4.5 ] [ 0 0 0 0 ] [ 1 2 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .03217 .04762 m .03217 .05387 L s [(-2)] .03217 .03512 0 1 Mshowa .28094 .04762 m .28094 .05387 L s [(-1)] .28094 .03512 0 1 Mshowa .77849 .04762 m .77849 .05387 L s [(1)] .77849 .03512 0 1 Mshowa .125 Mabswid .08192 .04762 m .08192 .05137 L s .13168 .04762 m .13168 .05137 L s .18143 .04762 m .18143 .05137 L s .23119 .04762 m .23119 .05137 L s .3307 .04762 m .3307 .05137 L s .38045 .04762 m .38045 .05137 L s .43021 .04762 m .43021 .05137 L s .47996 .04762 m .47996 .05137 L s .57947 .04762 m .57947 .05137 L s .62923 .04762 m .62923 .05137 L s .67898 .04762 m .67898 .05137 L s .72874 .04762 m .72874 .05137 L s .82825 .04762 m .82825 .05137 L s .878 .04762 m .878 .05137 L s .92776 .04762 m .92776 .05137 L s .97751 .04762 m .97751 .05137 L s .25 Mabswid 0 .04762 m 1 .04762 L s .52972 .52381 m .53597 .52381 L s [(5)] .51722 .52381 1 0 Mshowa .52972 1 m .53597 1 L s [(10)] .51722 1 1 0 Mshowa .52972 1.47619 m .53597 1.47619 L s [(15)] .51722 1.47619 1 0 Mshowa .52972 1.95238 m .53597 1.95238 L s [(20)] .51722 1.95238 1 0 Mshowa .125 Mabswid .52972 .14286 m .53347 .14286 L s .52972 .2381 m .53347 .2381 L s .52972 .33333 m .53347 .33333 L s .52972 .42857 m .53347 .42857 L s .52972 .61905 m .53347 .61905 L s .52972 .71429 m .53347 .71429 L s .52972 .80952 m .53347 .80952 L s .52972 .90476 m .53347 .90476 L s .52972 1.09524 m .53347 1.09524 L s .52972 1.19048 m .53347 1.19048 L s .52972 1.28571 m .53347 1.28571 L s .52972 1.38095 m .53347 1.38095 L s .52972 1.57143 m .53347 1.57143 L s .52972 1.66667 m .53347 1.66667 L s .52972 1.7619 m .53347 1.7619 L s .52972 1.85714 m .53347 1.85714 L s .25 Mabswid .52972 0 m .52972 2 L s 0 0 m 1 0 L 1 2 L 0 2 L closepath clip newpath .5 Mabswid .52972 .04762 m .62587 .08523 L .66861 .10439 L .70534 .12489 L .7201 .13579 L .73023 .14571 L .73393 .15049 L .73556 .15301 L .73705 .15568 L .73813 .15801 L .73911 .16058 L .73978 .16292 L .74024 .16511 L .74044 .16652 L .74056 .16782 L .74062 .16926 L .7406 .17079 L .7405 .17225 L .74033 .17362 L .73984 .17619 L .73911 .17873 L .73823 .18106 L .73551 .18633 L .73176 .19162 L .72648 .19743 L .71221 .20916 L .69558 .21964 L .67671 .22957 L .63425 .24818 L .51698 .29001 L .41239 .32645 L .36172 .34719 L .3409 .35735 L .3246 .36667 L .31155 .37576 L .30074 .3857 L .29663 .39088 L .29376 .39569 L .29257 .39836 L .29175 .40078 L .2914 .40212 L .2911 .40358 L .29079 .40622 L .29076 .40883 L .29099 .41127 L .29124 .41267 L .29153 .41395 L .29245 .41683 L .29365 .41959 L .29532 .4226 L Mistroke .29933 .42808 L .30414 .43314 L .31039 .43853 L .32451 .44823 L .36593 .46886 L .4253 .49114 L .55557 .5315 L .68074 .57067 L .73807 .59289 L .7599 .60354 L .77659 .61361 L .78856 .62319 L .79355 .6286 L .79547 .63118 L .79699 .6336 L .79813 .63579 L .79899 .6378 L .79972 .64003 L .80024 .64243 L .80048 .64462 L .80049 .647 L .80021 .64949 L .79969 .6518 L .79881 .65437 L .79751 .65714 L .79423 .66217 L .7899 .66707 L .78492 .67158 L .77016 .68184 L .74956 .69269 L .7025 .71161 L .639 .73192 L .50433 .76857 L .36063 .80814 L .29484 .83051 L .26982 .84119 L .25077 .85123 L .23791 .86005 L .23236 .86501 L .22836 .86956 L .22677 .87184 L .22536 .87427 L .22469 .87567 L .22418 .87694 L .22341 .87943 L .22296 .88199 L .22287 .88344 L .22288 .88477 L .22314 .88722 L .22378 .88982 L Mistroke .22475 .89236 L .22594 .89468 L .22967 .89991 L .23443 .90475 L .24009 .90935 L .25669 .91963 L .28076 .93085 L .33975 .95151 L .40672 .97023 L .5823 1.01218 L .66776 1.0325 L .74905 1.05395 L .80616 1.07226 L .83231 1.08278 L .8513 1.09236 L .85887 1.0971 L .86569 1.10225 L .87072 1.10712 L .87415 1.11162 L .87551 1.11408 L .87646 1.11638 L .87713 1.11891 L .87734 1.12033 L .87743 1.12163 L .87731 1.12411 L .87707 1.12551 L .87675 1.12681 L .87578 1.12945 L .87505 1.13092 L .87427 1.13226 L .87061 1.13705 L .86601 1.14146 L .85136 1.15142 L .83367 1.1602 L .81031 1.1695 L .74945 1.18839 L .66756 1.20868 L .49441 1.24531 L .31048 1.28484 L .22653 1.30719 L .19469 1.31787 L .17051 1.32791 L .15287 1.33763 L .14559 1.34311 L .14281 1.34571 L .14063 1.34814 L .13898 1.35037 L .13777 1.35239 L .13675 1.35464 L .13604 1.35705 L Mistroke .13573 1.35937 L .13581 1.36188 L .13628 1.36426 L .13707 1.36645 L .13828 1.36881 L .14001 1.37134 L .14426 1.37592 L .15118 1.38137 L .15927 1.38636 L .18246 1.39718 L .21137 1.40747 L .28132 1.42663 L .37495 1.44706 L .57205 1.48383 L .78206 1.52376 L .87684 1.54626 L .91239 1.55697 L .93904 1.56697 L .95708 1.57593 L .96466 1.58095 L .96994 1.58553 L .97223 1.58814 L .97389 1.59054 L .9752 1.59319 L .97565 1.59452 L .97599 1.596 L .97619 1.59846 L .97609 1.59984 L .97587 1.60111 L .97514 1.60345 L .9739 1.60592 L .9721 1.60849 L .97001 1.61084 L .96378 1.6161 L .95609 1.62097 L .94718 1.62557 L .92146 1.63589 L .88503 1.64708 L .80518 1.66599 L .69572 1.68681 L .43808 1.72889 L .30712 1.7505 L .19669 1.77103 L .11622 1.78945 L .08062 1.79988 L .05486 1.80954 L .04468 1.81443 L .03589 1.81971 L .03264 1.82213 L .02976 1.82469 L Mistroke .02771 1.8269 L .026 1.82927 L .0248 1.83158 L .02411 1.83368 L .02381 1.83579 L .02389 1.83779 L .02442 1.84012 L .02534 1.84227 L .02683 1.84463 L .02894 1.84713 L .03444 1.85187 L .04256 1.85696 L .06294 1.86617 L .09644 1.87713 L .13547 1.88717 L .30013 1.91833 L .52972 1.95238 L Mfstroke % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{144, 288}, ImageMargins->{{35, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgWoo0`003goo00<007ooOol0@Goo000mOol3000Ool3000` Ool004aoo`03001oogoo015oo`<002eoo`00C7oo00<007ooOol057oo0`00:Woo001goo 001Woo00<007ooOol017oo001Goo00<007ooOol01Goo 001Goo00<007ooOol0@Goo000; Ool4000mOol00`00Oomoo`11Ool000Yoo`03001oogoo03moo`03001oogoo045oo`001goo0`00@Woo 00<007ooOol0@Goo0006Ool00`00Oomoo`13Ool00`00Oomoo`11Ool000Eoo`03001oogoo04Aoo`03 001oogoo045oo`0017oo00<007ooOol0AGoo00<007ooOol0@Goo0004Ool20016Ool00`00Oomoo`11 Ool000Eoo`03001oogoo04Aoo`03001oogoo045oo`001Woo0`00@goo00<007ooOol0@Goo0009Ool0 0`00Oomoo`10Ool00`00Oomoo`11Ool000Yoo`@003ioo`03001oogoo045oo`003Woo0`00>goo00<0 07ooOol0@Goo000AOol3000hOol00`00Oomoo`11Ool001Aoo`@003Aoo`03001oogoo045oo`0067oo 1P00;Woo00<007ooOol0@Goo000NOol5000YOol00`00Oomoo`11Ool002=oo`H002=oo`03001oogoo 045oo`00:Goo1P007Goo00<007ooOol0@Goo000_Ool7000FOol00`00Oomoo`11Ool003Ioo`L000mo o`03001oogoo045oo`00?Goo1P000goo0P0017oo00<007ooOol0@Goo0010Ool00`00Oomoo`070002 Ool00`00Oomoo`11Ool0045oo`05001oogooOol00002Ool00`00Ool000030012Ool0049oo`04001o ogoo0002Ool01@00Oomoogoo0000@goo000oOol01000Oomoo`000Woo00@007ooOol000=oo`03001o ogoo045oo`00@7oo0P0017oo0P0017oo00<007ooOol0@Goo001"], ImageRangeCache->{{{0, 143}, {287, 0}} -> {-2.17143, -0.799267, 0.0284112, 0.0742138}}] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["Trapped Waves", "Section"], Cell[BoxData[ \(BuoyFreq := 0.01*Sech[8*\((z - lid/2)\)/lid]^2\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(\(buoyfreqplot = ParametricPlot[{BuoyFreq, z}, {z, 0, lid}, AspectRatio -> 2, Ticks -> {{ .005, .01}, Automatic}, AxesLabel -> {"\", "\"}]; \)\)], "Input"], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: 2 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Scaling calculations 0.0238095 95.2387 0.047619 0.0952381 [ [.5 .03512 -15 -9 ] [.5 .03512 15 0 ] [.9762 .03512 -12 -9 ] [.9762 .03512 12 0 ] [1.025 .04762 0 -6 ] [1.025 .04762 10 6 ] [.01131 .52381 -6 -4.5 ] [.01131 .52381 0 4.5 ] [.01131 1 -12 -4.5 ] [.01131 1 0 4.5 ] [.01131 1.47619 -12 -4.5 ] [.01131 1.47619 0 4.5 ] [.01131 1.95238 -12 -4.5 ] [.01131 1.95238 0 4.5 ] [.02381 2.025 -5 0 ] [.02381 2.025 5 12 ] [ 0 0 0 0 ] [ 1 2 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .5 .04762 m .5 .05387 L s [(0.005)] .5 .03512 0 1 Mshowa .9762 .04762 m .9762 .05387 L s [(0.01)] .9762 .03512 0 1 Mshowa 0 .04762 m 1 .04762 L s gsave 1.025 .04762 -61 -10 Mabsadd m 1 1 Mabs scale currentpoint translate 0 20 translate 1 -1 scale gsave 0.000000 0.000000 0.000000 setrgbcolor 1.000000 setlinewidth gsave newpath 61.000000 16.000000 moveto 558.000000 16.000000 lineto 558.000000 4.000000 lineto 61.000000 4.000000 lineto 61.000000 16.000000 lineto closepath clip newpath 63.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor 0.000000 0.000000 rmoveto 63.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor (N) show 69.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor 0.000000 0.000000 rmoveto 1.000000 setlinewidth grestore grestore %%DocumentNeededResources: font Courier %%DocumentSuppliedResources: %%DocumentNeededFonts: Courier %%DocumentSuppliedFonts: %%DocumentFonts: font Courier grestore .02381 .52381 m .03006 .52381 L s [(5)] .01131 .52381 1 0 Mshowa .02381 1 m .03006 1 L s [(10)] .01131 1 1 0 Mshowa .02381 1.47619 m .03006 1.47619 L s [(15)] .01131 1.47619 1 0 Mshowa .02381 1.95238 m .03006 1.95238 L s [(20)] .01131 1.95238 1 0 Mshowa .125 Mabswid .02381 .14286 m .02756 .14286 L s .02381 .2381 m .02756 .2381 L s .02381 .33333 m .02756 .33333 L s .02381 .42857 m .02756 .42857 L s .02381 .61905 m .02756 .61905 L s .02381 .71429 m .02756 .71429 L s .02381 .80952 m .02756 .80952 L s .02381 .90476 m .02756 .90476 L s .02381 1.09524 m .02756 1.09524 L s .02381 1.19048 m .02756 1.19048 L s .02381 1.28571 m .02756 1.28571 L s .02381 1.38095 m .02756 1.38095 L s .02381 1.57143 m .02756 1.57143 L s .02381 1.66667 m .02756 1.66667 L s .02381 1.7619 m .02756 1.7619 L s .02381 1.85714 m .02756 1.85714 L s .25 Mabswid .02381 0 m .02381 2 L s gsave .02381 2.025 -66 -4 Mabsadd m 1 1 Mabs scale currentpoint translate 0 20 translate 1 -1 scale gsave 0.000000 0.000000 0.000000 setrgbcolor 1.000000 setlinewidth gsave newpath 61.000000 16.000000 moveto 558.000000 16.000000 lineto 558.000000 4.000000 lineto 61.000000 4.000000 lineto 61.000000 16.000000 lineto closepath clip newpath 63.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor 0.000000 0.000000 rmoveto 63.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor (z) show 69.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor 0.000000 0.000000 rmoveto 1.000000 setlinewidth grestore grestore %%DocumentNeededResources: font Courier %%DocumentSuppliedResources: %%DocumentNeededFonts: Courier %%DocumentSuppliedFonts: %%DocumentFonts: font Courier grestore 0 0 m 1 0 L 1 2 L 0 2 L closepath clip newpath .5 Mabswid .02509 .04762 m .02556 .08523 L .02625 .12489 L .02723 .16511 L .02795 .18791 L .02876 .20916 L .03078 .24999 L .03332 .28715 L .0369 .32544 L .04232 .36701 L .04954 .40668 L .06023 .44874 L .07401 .48781 L .09142 .52441 L .11744 .5649 L .1493 .60195 L .19211 .63992 L .25326 .68137 L .32762 .72072 L .42572 .76266 L .66676 .84582 L .7918 .88756 L .84157 .90586 L .88972 .92596 L .92549 .94405 L .94153 .95401 L .9541 .96345 L .9598 .96858 L .96507 .97417 L .96869 .97881 L .97187 .98393 L .97407 .98875 L .97495 .99139 L .97557 .99387 L .97595 .99618 L .97615 .99829 L .97619 1.00061 L .97604 1.00306 L .97568 1.00553 L .97539 1.00692 L .97507 1.0082 L .97434 1.01052 L .97334 1.01306 L .97108 1.01748 L .968 1.02215 L .95847 1.03268 L .94643 1.04254 L .91772 1.06025 L .87795 1.07928 L .64319 1.16183 L Mistroke .41551 1.24135 L .25459 1.31784 L .19575 1.35726 L .14814 1.39924 L .11232 1.44214 L .08774 1.48251 L .06985 1.52274 L .05778 1.55971 L .04868 1.59744 L .04145 1.63879 L .03654 1.67789 L .03279 1.71974 L .03138 1.74014 L .0301 1.7623 L .0283 1.80257 L .02703 1.8422 L .02607 1.88424 L .02549 1.91941 L .02509 1.95238 L Mfstroke % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{144, 288}, ImageMargins->{{43, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHg7oo00@007ooOol000=oo`8000=oo`04001oogoo0002Ool01000Oomoo`001Goo 00<007ooOol07Woo00@007ooOol000=oo`8000=oo`04001oogoo0003Ool00`00Oomoo`06Ool003Qo o`04001oogoo0008Ool01000Oomoo`000Woo00@007ooOol000Eoo`03001oogoo01ioo`04001oogoo 0008Ool01000Oomoo`000goo00<007ooOol01Woo000=Ool00`00Oomoo`0XOol01000Oomoo`0027oo 00@007ooOol0009oo`04001oogoo0002Ool3000QOol01000Oomoo`0027oo00@007ooOol000=oo`03 001oogoo00Ioo`003Goo00<007ooOol0:7oo00@007ooOol000Qoo`04001oogoo0002Ool01000Oomo o`000Woo00<007ooOol08Goo00@007ooOol000Qoo`04001oogoo0003Ool00`00Oomoo`06Ool000eo o`03001oogoo02Uoo`8000Yoo`8000Aoo`8000=oo`@0025oo`8000Yoo`8000=oo`8000Qoo`003Goo 00<007ooOol0Mgoo0P000Woo00<007ooOol00Woo000=Ool00`00Oomoo`1hOol01000Ool0000017oo 000:Oomi0005Ool01000Ool0000017oo000=Ool00`00Oomoo`0fOol00`00Oomoo`0fOol00`00Oomo o`06Ool200000goo001oo`03Ool000eoo`03001oogoo07Qoo`800003Ool007oo00=oo`003Goo00<0 07ooOol0Mgoo0P000Woo0P000goo000=Ool00`00Oomoo`20Ool000eoo`03001oogoo081oo`003Goo 00<007ooOol0P7oo000=Ool00`00Oomoo`20Ool000eoo`03001oogoo081oo`003Goo00<007ooOol0 P7oo000=Ool00`00Oomoo`20Ool000eoo`03001oogoo081oo`003Goo00<007ooOol0P7oo000=Ool0 0`00Oomoo`20Ool000eoo`03001oogoo081oo`003Goo00<007ooOol0P7oo000=Ool00`00Oomoo`20 Ool000eoo`03001oogoo081oo`003Goo00<007ooOol0P7oo000=Ool20021Ool000eoo`80085oo`00 3Goo0P00PGoo000=Ool20021Ool000eoo`80085oo`003Goo0P00PGoo000=Ool20021Ool000eoo`80 085oo`003Goo0P00PGoo000=Ool20021Ool000eoo`80085oo`003Goo0P00PGoo000=Ool20021Ool0 00eoo`03001oo`00081oo`003Goo00<007oo0000P7oo000=Ool00`00Ool00020Ool000eoo`03001o o`00081oo`003Goo00<007oo0000P7oo000=Ool00`00Ool00020Ool000eoo`03001oo`00081oo`00 3Goo00<007oo0000P7oo000=Ool00`00Ool00020Ool000eoo`04001oogoo001oOol000eoo`04001o ogoo001oOol000eoo`04001oogoo001oOol000eoo`04001oogoo001oOol000eoo`04001oogoo001o Ool000eoo`05001oogooOol0001nOol000eoo`05001oogooOol0001nOol000eoo`05001oogooOol0 001nOol000eoo`05001oogooOol0001nOol000eoo`03001oogoo009oo`03001oogoo07]oo`003Goo 00<007ooOol00Woo00<007ooOol0Ngoo000=Ool00`00Oomoo`02Ool00`00Oomoo`1kOol000eoo`03 001oogoo00=oo`03001oogoo07Yoo`003Goo00<007ooOol00goo00<007ooOol0NWoo0007Ool30003 Ool00`00Oomoo`04Ool00`00Oomoo`1iOol000Yoo`04001oogoo0006Ool00`00Oomoo`1iOol000Yo o`03001oogoo008000Ioo`03001oogoo07Qoo`001goo0`000goo00<007ooOol01Goo00<007ooOol0 N7oo0007Ool00`00Oomoo`03Ool00`00Oomoo`06Ool00`00Oomoo`1gOol000Moo`@0009oo`03001o ogoo00Moo`03001oogoo07Ioo`003Goo00<007ooOol01goo00<007ooOol0MWoo000=Ool00`00Oomo o`08Ool00`00Oomoo`1eOol000eoo`03001oogoo00Uoo`03001oogoo07Aoo`003Goo00<007ooOol0 2Woo00<007ooOol0Lgoo000=Ool00`00Oomoo`0;Ool00`00Oomoo`1bOol000eoo`03001oogoo00ao o`03001oogoo075oo`003Goo00<007ooOol03Goo00<007ooOol0L7oo000=Ool00`00Oomoo`0>Ool0 0`00Oomoo`1_Ool000eoo`03001oogoo00moo`03001oogoo06ioo`003Goo00<007ooOol047oo00<0 07ooOol0KGoo000=Ool00`00Oomoo`0AOol2001]Ool000eoo`03001oogoo01=oo`03001oogoo06Yo o`003Goo00<007ooOol057oo0P00JWoo000=Ool00`00Oomoo`0FOol00`00Oomoo`1WOol000eoo`03 001oogoo01Moo`03001oogoo06Ioo`003Goo00<007ooOol067oo0P00IWoo000=Ool00`00Oomoo`0J Ool2001TOol000eoo`03001oogoo01aoo`80069oo`003Goo00<007ooOol07Woo00<007ooOol0Ggoo 000=Ool00`00Oomoo`0OOol2001OOol000eoo`03001oogoo025oo`8005eoo`003Goo00<007ooOol0 8goo0`00FWoo000=Ool00`00Oomoo`0VOol2001HOol000eoo`03001oogoo02Qoo`8005Ioo`003Goo 00<007ooOol0:Woo0`00Dgoo000=Ool00`00Oomoo`0]Ool2001AOol000eoo`03001oogoo02moo`<0 04ioo`003Goo00<007ooOol0goo0`00@Woo000=Ool00`00Oomoo`0nOol3000oOol0 00eoo`03001oogoo045oo`<003aoo`003Goo00<007ooOol0A7oo0`00>Goo000=Ool00`00Oomoo`17 Ool3000fOol000eoo`03001oogoo04Yoo`<003=oo`003Goo00<007ooOol0CGoo1000;goo000=Ool0 0`00Oomoo`1AOol3000/Ool000eoo`03001oogoo05Aoo`@002Qoo`003Goo00<007ooOol0F7oo0`00 9Goo000=Ool00`00Oomoo`1KOol2000SOol000eoo`03001oogoo05eoo`80025oo`003Goo00<007oo Ool0Ggoo0P007goo000=Ool00`00Oomoo`1QOol3000LOol000eoo`03001oogoo06Aoo`8001Yoo`00 3Goo00<007ooOol0IWoo0`005goo000=Ool00`00Oomoo`1YOol2000EOol000eoo`03001oogoo06]o o`03001oogoo019oo`003Goo00<007ooOol0Jgoo00<007ooOol04Woo000=Ool00`00Oomoo`1/Ool2 000BOol000eoo`03001oogoo06eoo`80015oo`0000Aoo`00000000Aoo`8000=oo`03001oogoo06io o`03001oogoo00moo`000Woo00<007ooOol00Woo00@007ooOol0009oo`03001oogoo06ioo`80011o o`000Woo00<007ooOol00Woo00@007ooOol0009oo`80071oo`03001oogoo00ioo`000Woo00<007oo Ool00Woo00@007ooOol0009oo`03001oogoo06ioo`80011oo`000Woo00<007ooOol00Woo00@007oo Ool0009oo`03001oogoo06ioo`03001oogoo00moo`0000=oo`0000001Goo0P000goo00<007ooOol0 KWoo00<007ooOol03goo000=Ool00`00Oomoo`1/Ool2000BOol000eoo`03001oogoo06Yoo`8001Ao o`003Goo00<007ooOol0JGoo00<007ooOol057oo000=Ool00`00Oomoo`1WOol2000GOol000eoo`03 001oogoo06Aoo`<001Uoo`003Goo00<007ooOol0HWoo0P0077oo000=Ool00`00Oomoo`1OOol3000N Ool000eoo`03001oogoo05aoo`<0025oo`003Goo00<007ooOol0FGoo0`0097oo000=Ool00`00Oomo o`1FOol3000WOol000eoo`03001oogoo05Aoo`8002Yoo`003Goo00<007ooOol0DGoo0`00;7oo000= Ool00`00Oomoo`1>Ool3000_Ool000eoo`03001oogoo04]oo`<0039oo`003Goo00<007ooOol0B7oo 0`00=Goo000=Ool00`00Oomoo`16Ool2000hOol000eoo`03001oogoo04=oo`<003Yoo`003Goo00<0 07ooOol0@7oo0`00?Goo000=Ool00`00Oomoo`0nOol20010Ool000eoo`03001oogoo03]oo`<0049o o`003Goo00<007ooOol0>7oo0`00AGoo000=Ool00`00Oomoo`0eOol30018Ool000eoo`03001oogoo 03=oo`8004]oo`003Goo00<007ooOol0<7oo0`00CGoo000=Ool00`00Oomoo`0]Ool3001@Ool000eo o`03001oogoo02]oo`8005=oo`003Goo00<007ooOol0:Goo0P00EGoo000=Ool00`00Oomoo`0WOol2 001GOol000eoo`03001oogoo02Eoo`8005Uoo`003Goo00<007ooOol08Woo0`00Fgoo000=Ool00`00 Oomoo`0POol2001NOol000eoo`03001oogoo01ioo`80061oo`003Goo00<007ooOol077oo0P00HWoo 000=Ool00`00Oomoo`0JOol2001TOol000eoo`03001oogoo01Qoo`8006Ioo`003Goo00<007ooOol0 5goo00<007ooOol0IWoo000=Ool00`00Oomoo`0EOol2001YOol000eoo`03001oogoo01Aoo`03001o ogoo06Uoo`003Goo00<007ooOol04goo00<007ooOol0JWoo000=Ool00`00Oomoo`0AOol2001]Ool0 00eoo`03001oogoo011oo`03001oogoo06eoo`003Goo00<007ooOol03goo00<007ooOol0KWoo000= Ool00`00Oomoo`0>Ool00`00Oomoo`1_Ool000eoo`03001oogoo00eoo`03001oogoo071oo`003Goo 00<007ooOol037oo00<007ooOol0LGoo000=Ool00`00Oomoo`0"], ImageRangeCache->{{{0, 143}, {287, 0}} -> {-0.00118688, -2.1755, 8.79652*^-5, 0.0879657}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(Display["\", buoyfreqplot, "\"]; \)\)], "Input"], Cell[OutputFormData["\<\ Graphics[{{Line[{{0.00001340951576393732, 8.333333333333332*^-7}, {0.00001838776398338612, 0.3949634152604035}, {0.00002564690798088263, 0.8113398314583159}, {0.00003593604480395857, 1.233631650861948}, {0.00004350737136797001, 1.473093217150288}, {0.00005198582875441859, 1.696175997187473}, {0.0000731792675188476, 2.124933078177035}, {0.00009984778460256594, 2.515029157291579}, {0.000137470985590519, 2.917109255604063}, {0.0001943555230595457, 3.353554648854}, {0.0002701811903409296, 3.770132622869616}, {0.0003824635231465732, 4.211752471874737}, {0.0005270680571806353, 4.621981284491925}, {0.0007099431319227855, 5.006289293020456}, {0.0009831328470784622, 5.431480382090784}, {0.001317686099548589, 5.820498445624491}, {0.001767148743787187, 6.219189937319843}, {0.002409220395264368, 6.654379929686842}, {0.003189959860422914, 7.067514482012424}, {0.004220075550273834, 7.507933745207509}, {0.006750898424020856, 8.381159892186492}, {0.008063866345255751, 8.819375903436599}, {0.00858640169022054, 9.011580542663294}, {0.009092024229026941, 9.222555538227966}, {0.009467549543182187, 9.412549201579945}, {0.009636041033599331, 9.517141293485588}, {0.009767998680493848, 9.616223787601795}, {0.009827880408022743, 9.670112673203112}, {0.009883180101164279, 9.728732237303455}, {0.009921201132781138, 9.777492507109071}, {0.00995461010343175, 9.831314408904146}, {0.00997771292981757, 9.881889243978784}, {0.009986936447993087, 9.909601782005712}, {0.009993370997484284, 9.935618647965331}, {0.009997420615810028, 9.959845393898302}, {0.009999484974882694, 9.982058379368929}, {0.009999933532357484, 10.00644534410668}, {0.009998349400378591, 10.03212068333193}, {0.009994611892058342, 10.05804117174831}, {0.009991557439924379, 10.0726606622486}, {0.009988137230508327, 10.08614000130273}, {0.009980515609475355, 10.11042457217502}, {0.009969969481355908, 10.13713765790263}, {0.009946324482815921, 10.18348770704075}, {0.009913962606418531, 10.23255935727303}, {0.009813896814687527, 10.3431887234329}, {0.009687420742950227, 10.44669143714493}, {0.009386010278503093, 10.63263725519093}, {0.008968452246633949, 10.83239781343575}, {0.006503490432451081, 11.69920609989365}, {0.004112857557955896, 12.53418388541404}, {0.002423173205102296, 13.33733116999691}, {0.001805395511676562, 13.75127743710995}, {0.001305480516108419, 14.19198128697562}, {0.0009293449222580296, 14.64247381638821}, {0.0006712563563819785, 15.06630373541264}, {0.0004834690786093283, 15.48877936760137}, {0.0003567312073293161, 15.87696518197464}, {0.0002611219451160392, 16.27307672710099}, {0.0001852062864603217, 16.70729895999495}, {0.0001337090479605842, 17.11781169076203}, {0.00009426037610035994, 17.55730506947358}, {0.00007947546626036896, 17.77150086315877}, {0.00006602087426641118, 18.00419314727982}, {0.00004711953985327765, 18.42698351041053}, {0.00003380118617294877, 18.84305828972868}, {0.00002375618312916946, 19.28450378186829}, {0.00001768520096086855, 19.65377712511259}, {0.00001340951576393733, 19.99999916666666}}]}}, {PlotRange -> Automatic, AspectRatio -> 2, DisplayFunction :> $DisplayFunction, ColorOutput -> Automatic, Axes -> Automatic, AxesOrigin -> Automatic, PlotLabel -> None, AxesLabel -> {\"N\", \"z\"}, Ticks -> {{0.005, 0.01}, Automatic}, GridLines -> None, Prolog -> {}, Epilog -> {}, AxesStyle -> Automatic, Background -> Automatic, DefaultColor -> Automatic, DefaultFont :> $DefaultFont, RotateLabel -> True, Frame -> False, FrameStyle -> Automatic, FrameTicks -> Automatic, FrameLabel -> None, PlotRegion -> Automatic, ImageSize -> Automatic, TextStyle :> $TextStyle, FormatType :> $FormatType}]\ \>", "\<\ -Graphics-\ \>"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(wAtlid[ .001, 1]\)], "Input"], Cell[BoxData[ \(\(-23379.4113751917986`\)\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(Plot[wAtlid[\[Omega], 1], {\[Omega], .002, .008}]; \)\)], "Input"], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Scaling calculations -0.293651 158.73 0.456541 2.77143e-08 [ [.18254 .44404 -15 -9 ] [.18254 .44404 15 0 ] [.34127 .44404 -15 -9 ] [.34127 .44404 15 0 ] [.5 .44404 -15 -9 ] [.5 .44404 15 0 ] [.65873 .44404 -15 -9 ] [.65873 .44404 15 0 ] [.81746 .44404 -15 -9 ] [.81746 .44404 15 0 ] [.97619 .44404 -15 -9 ] [.97619 .44404 15 0 ] [.01131 .04083 -46.5 -6.25 ] [.01131 .04083 0 6.25 ] [.01131 .1794 -34.5 -6.25 ] [.01131 .1794 0 6.25 ] [.01131 .31797 -34.5 -6.25 ] [.01131 .31797 0 6.25 ] [.01131 .59511 -28.5 -6.25 ] [.01131 .59511 0 6.25 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .18254 .45654 m .18254 .46279 L s [(0.003)] .18254 .44404 0 1 Mshowa .34127 .45654 m .34127 .46279 L s [(0.004)] .34127 .44404 0 1 Mshowa .5 .45654 m .5 .46279 L s [(0.005)] .5 .44404 0 1 Mshowa .65873 .45654 m .65873 .46279 L s [(0.006)] .65873 .44404 0 1 Mshowa .81746 .45654 m .81746 .46279 L s [(0.007)] .81746 .44404 0 1 Mshowa .97619 .45654 m .97619 .46279 L s [(0.008)] .97619 .44404 0 1 Mshowa .125 Mabswid .05556 .45654 m .05556 .46029 L s .0873 .45654 m .0873 .46029 L s .11905 .45654 m .11905 .46029 L s .15079 .45654 m .15079 .46029 L s .21429 .45654 m .21429 .46029 L s .24603 .45654 m .24603 .46029 L s .27778 .45654 m .27778 .46029 L s .30952 .45654 m .30952 .46029 L s .37302 .45654 m .37302 .46029 L s .40476 .45654 m .40476 .46029 L s .43651 .45654 m .43651 .46029 L s .46825 .45654 m .46825 .46029 L s .53175 .45654 m .53175 .46029 L s .56349 .45654 m .56349 .46029 L s .59524 .45654 m .59524 .46029 L s .62698 .45654 m .62698 .46029 L s .69048 .45654 m .69048 .46029 L s .72222 .45654 m .72222 .46029 L s .75397 .45654 m .75397 .46029 L s .78571 .45654 m .78571 .46029 L s .84921 .45654 m .84921 .46029 L s .88095 .45654 m .88095 .46029 L s .9127 .45654 m .9127 .46029 L s .94444 .45654 m .94444 .46029 L s .25 Mabswid 0 .45654 m 1 .45654 L s .02381 .04083 m .03006 .04083 L s gsave .01131 .04083 -46.5 -4.25 Mabsadd m 1 1 Mabs scale /Courier findfont 10 scalefont setfont (-1.5) show /Math1Mono findfont 10 scalefont setfont (\\264) show /Courier findfont 10 scalefont setfont (10) show 0 5.25 rmoveto /Courier findfont 7.5 scalefont setfont (7) show grestore .02381 .1794 m .03006 .1794 L s gsave .01131 .1794 -34.5 -4.25 Mabsadd m 1 1 Mabs scale /Courier findfont 10 scalefont setfont (-1) show /Math1Mono findfont 10 scalefont setfont (\\264) show /Courier findfont 10 scalefont setfont (10) show 0 5.25 rmoveto /Courier findfont 7.5 scalefont setfont (7) show grestore .02381 .31797 m .03006 .31797 L s gsave .01131 .31797 -34.5 -4.25 Mabsadd m 1 1 Mabs scale /Courier findfont 10 scalefont setfont (-5) show /Math1Mono findfont 10 scalefont setfont (\\264) show /Courier findfont 10 scalefont setfont (10) show 0 5.25 rmoveto /Courier findfont 7.5 scalefont setfont (6) show grestore .02381 .59511 m .03006 .59511 L s gsave .01131 .59511 -28.5 -4.25 Mabsadd m 1 1 Mabs scale /Courier findfont 10 scalefont setfont (5) show /Math1Mono findfont 10 scalefont setfont (\\264) show /Courier findfont 10 scalefont setfont (10) show 0 5.25 rmoveto /Courier findfont 7.5 scalefont setfont (6) show grestore .125 Mabswid .02381 .06854 m .02756 .06854 L s .02381 .09626 m .02756 .09626 L s .02381 .12397 m .02756 .12397 L s .02381 .15168 m .02756 .15168 L s .02381 .20711 m .02756 .20711 L s .02381 .23483 m .02756 .23483 L s .02381 .26254 m .02756 .26254 L s .02381 .29026 m .02756 .29026 L s .02381 .34568 m .02756 .34568 L s .02381 .3734 m .02756 .3734 L s .02381 .40111 m .02756 .40111 L s .02381 .42883 m .02756 .42883 L s .02381 .48426 m .02756 .48426 L s .02381 .51197 m .02756 .51197 L s .02381 .53968 m .02756 .53968 L s .02381 .5674 m .02756 .5674 L s .02381 .01311 m .02756 .01311 L s .25 Mabswid .02381 0 m .02381 .61803 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath .5 Mabswid .02381 .4525 m .03279 .46179 L .03793 .46672 L .04262 .47018 L .04384 .47087 L .04517 .47151 L .04631 .47196 L .04758 .47235 L .04827 .47252 L .04892 .47264 L .0496 .47273 L .05034 .4728 L .05102 .47282 L .05167 .47281 L .05288 .47271 L .05403 .47251 L .05525 .47219 L .05641 .47179 L .05745 .47135 L .06002 .46996 L .06244 .46827 L .07312 .45769 L .08439 .44476 L .08903 .44016 L .09401 .43622 L .09687 .43453 L .09948 .4334 L .10073 .433 L .1014 .43283 L .1021 .43267 L .10339 .43247 L .10458 .43238 L .10574 .43237 L .10695 .43245 L .10799 .43259 L .10913 .43283 L .11036 .43316 L .11167 .43363 L .11403 .43472 L .11868 .43776 L .1229 .44144 L .14313 .46622 L .15245 .47794 L .16232 .48787 L .16747 .49155 L .1704 .4931 L .1718 .4937 L .17314 .49418 L .17434 .49453 L .17565 .49483 L Mistroke .17684 .49503 L .17795 .49515 L .17926 .49522 L .18045 .4952 L .18111 .49515 L .18182 .49508 L .18311 .49489 L .18436 .49462 L .18549 .49431 L .18804 .49338 L .19041 .49223 L .19264 .49091 L .19765 .48712 L .20309 .48184 L .22417 .45328 L .24359 .42347 L .2539 .4096 L .25931 .40339 L .265 .39783 L .27006 .39382 L .27274 .39209 L .2756 .39054 L .27823 .38939 L .28062 .38858 L .28178 .38828 L .28303 .38801 L .28421 .38781 L .28529 .38768 L .28639 .3876 L .28757 .38756 L .28869 .38758 L .2897 .38765 L .29089 .38777 L .29219 .38798 L .29341 .38823 L .29454 .38852 L .29679 .38925 L .29928 .39029 L .30432 .39314 L .30923 .39684 L .3144 .40167 L .32364 .41253 L .34449 .44534 L .38257 .51838 L .40379 .55561 L .4231 .58188 L .43272 .5914 L .43779 .59535 L .44316 .59869 L .44785 .60088 L Mistroke .45052 .60182 L .45301 .60249 L .45419 .60274 L .4553 .60294 L .45632 .60308 L .4574 .6032 L .45857 .60328 L .45965 .60332 L .46094 .60332 L .46211 .60327 L .46331 .60317 L .46444 .60304 L .46567 .60285 L .46699 .60259 L .46969 .6019 L .47214 .60106 L .47688 .5989 L .48137 .5962 L .49135 .58797 L .5022 .57566 L .52199 .54488 L .543 .50205 L .58077 .40638 L .62099 .29269 L .6597 .18774 L .67936 .1405 L .70086 .09579 L .72163 .06084 L .7405 .03711 L .7506 .02777 L .76013 .02122 L .76536 .01856 L .77019 .01673 L .77271 .016 L .77548 .01539 L .77688 .01516 L .77838 .01496 L .7791 .01488 L .77979 .01483 L .78108 .01475 L .78234 .01472 L .78304 .01472 L .78369 .01473 L .78488 .01477 L .78613 .01487 L .78726 .01498 L .78832 .01512 L .79071 .01554 L .79329 .01616 L .79565 .01688 L Mistroke .80094 .01902 L .80548 .02144 L .8104 .02469 L .81935 .03221 L .82878 .04244 L .83918 .05643 L .8602 .0934 L .8796 .13767 L .90048 .19594 L .94169 .34189 L .97619 .49353 L Mfstroke % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{35, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgaoo`007oo00<007ooOol07Woo000aOol00`00Oomoo`2@Ool00`00Oomoo`0hOol00`00Oomoo`0NOol0 035oo`03001oogoo08moo`03001oogoo03Uoo`03001oogoo01ioo`00Goo00<007ooOol07Woo000aOol00`00Oomoo`2?Ool00`00Oomoo`0jOol00`00Oomoo`0M Ool0035oo`03001oogoo08ioo`03001oogoo03]oo`03001oogoo01eoo`00goo00<007ooOol07Goo000aOol00`00Oomoo`2=Ool00`00Oomoo`0mOol00`00Oomoo`0LOol0 035oo`03001oogoo08eoo`03001oogoo03eoo`03001oogoo01aoo`00 Ool00`00Oomoo`0EOol0035oo`03001oogoo089oo`03001oogoo04moo`03001oogoo01Eoo`00Woo00<007ooOol0EGoo00<007ooOol04Woo000aOol00`00Oomoo`0EOol20004Ool20004Ool20004 Ool20003Ool30003Ool00`00Oomoo`05Ool20004Ool2000017oo001oogoo0P0017oo0P0017oo0`00 2Woo0P0017oo0P0017oo0P0017oo0P000goo0`001Goo00<007ooOol00goo0P0017oo0P0017oo0P00 17oo0P0017oo0P002goo0P0017oo0P0017oo0P0017oo0P0017oo00<007ooOol02goo0P0017oo0`00 0goo0P0017oo0P0017oo0P000goo000aOol00`00Oomoo`0DOol01000Oomoo`000goo0P000goo00@0 07ooOol0009oo`04001oogoo0005Ool00`00Ool00007Ool01000Oomoo`000goo0P0000Aoo`00Ool0 009oo`04001oogoo0002Ool00`00Oomoo`02Ool00`00Oomoo`08Ool01000Oomoo`000goo0P000goo 00@007ooOol0009oo`04001oogoo0005Ool00`00Oomoo`02Ool00`00Oomoo`02Ool01000Oomoo`00 0goo0P000goo00@007ooOol0009oo`04001oogoo0002Ool01000Oomoo`002Goo00@007ooOol000=o o`8000=oo`04001oogoo0002Ool01000Oomoo`000goo00<007ooOol02Woo00@007ooOol000=oo`<0 009oo`04001oogoo0002Ool01000Oomoo`000Woo00@007ooOol0009oo`00Ool0035oo`80 01ioo`03001oogoo00Eoo`8002Ioo`03001oogoo02Qoo`03001oogoo061oo`03001oogoo00ioo`00 aoo`00;7oo00<007ooOol00Woo00<007ooOol0k7oo000]Ool20002Ool0 0`00Oomoo`3/Ool00?moob5oo`00ogoo8Goo003oOolQOol00?moob5oo`00ogoo8Goo003oOolQOol0 0?moob5oo`00ogoo8Goo003oOolQOol00?moob5oo`00ogoo8Goo003oOolQOol00?moob5oo`00ogoo 8Goo003oOolQOol00?moob5oo`00\ \>"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {0.000667724, -1.89131*^7, 2.6957*^-5, 154393}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(aroot = FindRoot[wAtlid[\[Omega], 1] == 0, {\[Omega], .003, .0035}]\)], "Input"], Cell[BoxData[ \(FindRoot::"frmp" \( : \ \) "Machine precision is insufficient to achieve the accuracy \ \!\(1.00000000000000066`*^-6\)."\)], "Message"], Cell[BoxData[ \({\[Omega] \[Rule] 0.00324882306888461958`}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(hit = \[Omega] /. aroot\)], "Input"], Cell[BoxData[ \(0.00324882306888461958`\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(hitsoln = soln[hit, 1]\)], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{"w", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.`, 20.`}}, "<>"]\), False, Editable->False]}], "}"}], "}"}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(wmax = Max[Table[w[z] /. hitsoln[\([1]\)], {z, 0, lid, .001*lid}]]\)], "Input"], Cell[BoxData[ \(999.820124732531922`\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(wplot = ParametricPlot[{\((1/wmax)\) w[z] /. hitsoln[\([1]\)], z}, {z, 0, lid}, AspectRatio -> 2, AxesLabel -> {"\", "\"}]; \)\)], "Input"], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: 2 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Scaling calculations 0.550609 0.425625 0.047619 0.0952381 [ [.12498 .03512 -6 -9 ] [.12498 .03512 6 0 ] [.3378 .03512 -12 -9 ] [.3378 .03512 12 0 ] [.76342 .03512 -9 -9 ] [.76342 .03512 9 0 ] [.97623 .03512 -3 -9 ] [.97623 .03512 3 0 ] [1.025 .04762 0 -6 ] [1.025 .04762 10.125 6 ] [.53811 .52381 -6 -4.5 ] [.53811 .52381 0 4.5 ] [.53811 1 -12 -4.5 ] [.53811 1 0 4.5 ] [.53811 1.47619 -12 -4.5 ] [.53811 1.47619 0 4.5 ] [.53811 1.95238 -12 -4.5 ] [.53811 1.95238 0 4.5 ] [.55061 2.025 -5 0 ] [.55061 2.025 5 12 ] [ 0 0 0 0 ] [ 1 2 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .12498 .04762 m .12498 .05387 L s [(-1)] .12498 .03512 0 1 Mshowa .3378 .04762 m .3378 .05387 L s [(-0.5)] .3378 .03512 0 1 Mshowa .76342 .04762 m .76342 .05387 L s [(0.5)] .76342 .03512 0 1 Mshowa .97623 .04762 m .97623 .05387 L s [(1)] .97623 .03512 0 1 Mshowa .125 Mabswid .16755 .04762 m .16755 .05137 L s .21011 .04762 m .21011 .05137 L s .25267 .04762 m .25267 .05137 L s .29523 .04762 m .29523 .05137 L s .38036 .04762 m .38036 .05137 L s .42292 .04762 m .42292 .05137 L s .46548 .04762 m .46548 .05137 L s .50805 .04762 m .50805 .05137 L s .59317 .04762 m .59317 .05137 L s .63573 .04762 m .63573 .05137 L s .6783 .04762 m .6783 .05137 L s .72086 .04762 m .72086 .05137 L s .80598 .04762 m .80598 .05137 L s .84855 .04762 m .84855 .05137 L s .89111 .04762 m .89111 .05137 L s .93367 .04762 m .93367 .05137 L s .08242 .04762 m .08242 .05137 L s .03986 .04762 m .03986 .05137 L s .25 Mabswid 0 .04762 m 1 .04762 L s gsave 1.025 .04762 -61 -10 Mabsadd m 1 1 Mabs scale currentpoint translate 0 20 translate 1 -1 scale gsave 0.000000 0.000000 0.000000 setrgbcolor 1.000000 setlinewidth gsave newpath 61.000000 16.000000 moveto 558.000000 16.000000 lineto 558.000000 4.000000 lineto 61.000000 4.000000 lineto 61.000000 16.000000 lineto closepath clip newpath 63.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor 0.000000 0.000000 rmoveto 63.125000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor (w) show 69.125000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor 0.000000 0.000000 rmoveto 1.000000 setlinewidth grestore grestore %%DocumentNeededResources: font Courier %%DocumentSuppliedResources: %%DocumentNeededFonts: Courier %%DocumentSuppliedFonts: %%DocumentFonts: font Courier grestore .55061 .52381 m .55686 .52381 L s [(5)] .53811 .52381 1 0 Mshowa .55061 1 m .55686 1 L s [(10)] .53811 1 1 0 Mshowa .55061 1.47619 m .55686 1.47619 L s [(15)] .53811 1.47619 1 0 Mshowa .55061 1.95238 m .55686 1.95238 L s [(20)] .53811 1.95238 1 0 Mshowa .125 Mabswid .55061 .14286 m .55436 .14286 L s .55061 .2381 m .55436 .2381 L s .55061 .33333 m .55436 .33333 L s .55061 .42857 m .55436 .42857 L s .55061 .61905 m .55436 .61905 L s .55061 .71429 m .55436 .71429 L s .55061 .80952 m .55436 .80952 L s .55061 .90476 m .55436 .90476 L s .55061 1.09524 m .55436 1.09524 L s .55061 1.19048 m .55436 1.19048 L s .55061 1.28571 m .55436 1.28571 L s .55061 1.38095 m .55436 1.38095 L s .55061 1.57143 m .55436 1.57143 L s .55061 1.66667 m .55436 1.66667 L s .55061 1.7619 m .55436 1.7619 L s .55061 1.85714 m .55436 1.85714 L s .25 Mabswid .55061 0 m .55061 2 L s gsave .55061 2.025 -66 -4 Mabsadd m 1 1 Mabs scale currentpoint translate 0 20 translate 1 -1 scale gsave 0.000000 0.000000 0.000000 setrgbcolor 1.000000 setlinewidth gsave newpath 61.000000 16.000000 moveto 558.000000 16.000000 lineto 558.000000 4.000000 lineto 61.000000 4.000000 lineto 61.000000 16.000000 lineto closepath clip newpath 63.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor 0.000000 0.000000 rmoveto 63.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor (z) show 69.000000 13.000000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000000 0.000000 0.000000 setrgbcolor 0.000000 0.000000 rmoveto 1.000000 setlinewidth grestore grestore %%DocumentNeededResources: font Courier %%DocumentSuppliedResources: %%DocumentNeededFonts: Courier %%DocumentSuppliedFonts: %%DocumentFonts: font Courier grestore 0 0 m 1 0 L 1 2 L 0 2 L closepath clip newpath .5 Mabswid .55061 .04762 m .55079 .08523 L .55101 .12489 L .55132 .16511 L .55156 .18791 L .55183 .20916 L .55217 .23015 L .55257 .24999 L .55358 .28715 L .5543 .30654 L .55524 .32709 L .55758 .36397 L .55919 .3828 L .56134 .40304 L .56667 .43964 L .57067 .45985 L .57617 .48191 L .58281 .50299 L .59028 .52209 L .6114 .56144 L .62725 .58304 L .64571 .60338 L .69458 .64354 L .7569 .68048 L .84424 .72133 L .88999 .74127 L .92892 .75946 L .94565 .76848 L .96069 .77828 L .96616 .78277 L .97082 .78754 L .97374 .79162 L .97493 .79394 L .9757 .79611 L .97615 .79841 L .97619 .80092 L .9757 .80354 L .97472 .80598 L .97334 .80827 L .97167 .81035 L .96644 .81502 L .96256 .81764 L .95837 .82005 L .94657 .82551 L .9165 .8354 L .87731 .84475 L .8233 .85479 L .69699 .87293 L .5167 .89435 L .35428 .91421 L Mistroke .2899 .92402 L .26362 .92905 L .24132 .93447 L .23168 .93758 L .22822 .93895 L .22498 .94044 L .2226 .94176 L .22057 .94318 L .21902 .94467 L .21812 .94606 L .21775 .94742 L .21794 .94889 L .21867 .95027 L .21983 .95156 L .2233 .95398 L .2288 .95656 L .23587 .95904 L .24384 .96131 L .26694 .96643 L .32549 .97568 L .41151 .98596 L .5174 .99676 L .70925 1.0159 L .7975 1.02641 L .83259 1.03157 L .85912 1.03632 L .87867 1.04082 L .88637 1.04309 L .89312 1.04558 L .89808 1.04808 L .90005 1.04948 L .90135 1.05076 L .90248 1.05324 L .90189 1.0556 L .89998 1.0577 L .89656 1.05994 L .89252 1.06186 L .88687 1.06397 L .87003 1.06867 L .84913 1.07308 L .79039 1.08246 L .70829 1.09271 L .5132 1.1134 L .33562 1.13242 L .19682 1.15039 L .13774 1.16023 L .08826 1.17091 L .07121 1.17562 L .05623 1.18061 L .04497 1.18529 L Mistroke .0369 1.18957 L .03009 1.19453 L .02743 1.19726 L .0256 1.19986 L .02449 1.20225 L .02393 1.20447 L .02381 1.20691 L .02423 1.2095 L .02472 1.21099 L .02532 1.21235 L .02719 1.21543 L .0294 1.21815 L .03226 1.22102 L .04674 1.23145 L .06532 1.24138 L .08524 1.25043 L .13621 1.27076 L .23242 1.30744 L .32651 1.34886 L .39292 1.38717 L .44004 1.42436 L .45868 1.44338 L .47538 1.46386 L .48903 1.4842 L .49939 1.50284 L .51497 1.53941 L .5213 1.55903 L .52617 1.57728 L .53085 1.59859 L .53444 1.61869 L .53733 1.63835 L .53988 1.65964 L .54322 1.69682 L .54451 1.71577 L .54568 1.73661 L .54667 1.75861 L .54742 1.7787 L .54852 1.81763 L .54935 1.85942 L .54992 1.8975 L .55044 1.93898 L .5506 1.95238 L Mfstroke % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{160.062, 320.125}, ImageMargins->{{35, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgOol0051oo`8004ioo`00D7oo0P00CWoo001@Ool2 001>Ool0051oo`8004ioo`00D7oo0P00CWoo001@Ool2001>Ool0051oo`8004ioo`00D7oo0P00CWoo 001@Ool00`00Ool0001=Ool0051oo`03001oo`0004eoo`00D7oo00<007oo0000CGoo001@Ool00`00 Ool0001=Ool0051oo`03001oo`0004eoo`00D7oo00<007oo0000CGoo001@Ool00`00Ool0001=Ool0 051oo`04001oogoo001Woo001@Ool00`00Oomoo`0COol00`00Oomoo`0g Ool0051oo`03001oogoo01Aoo`8003Moo`00D7oo00<007ooOol05Woo00<007ooOol0=7oo001@Ool0 0`00Oomoo`0GOol2000dOol0051oo`03001oogoo01Uoo`03001oogoo035oo`00D7oo00<007ooOol0 6Woo0P00Goo0P004Woo001@Ool00`00Oomoo`0jOol00`00Oomoo`0@Ool0051oo`03001oogoo 03Yoo`03001oogoo011oo`00D7oo00<007ooOol0>Woo00<007ooOol047oo001@Ool00`00Oomoo`0i Ool2000BOol0051oo`03001oogoo03Moo`8001Aoo`00D7oo00<007ooOol0=7oo0`005Woo001@Ool0 0`00Oomoo`0aOol3000IOol0051oo`03001oogoo02ioo`<001aoo`00D7oo00<007ooOol0:7oo1P00 7goo001@Ool00`00Oomoo`0QOol7000UOol0051oo`03001oogoo01]oo`H002aoo`00D7oo00<007oo Ool05Goo1P00Ool7000hOol0051oo`03001oogoo00Eoo`T003moo`00 Cgoo2@00B7oo0017Ool800000goo001oo`1>Ool003moo`P000Uoo`03001oogoo04eoo`00=goo2000 4Goo00<007ooOol0CGoo000_Ool8000IOol00`00Oomoo`1=Ool002Qoo`L0025oo`03001oogoo04eo o`009Goo0`00:7oo00<007ooOol0CGoo000POol5000[Ool00`00Oomoo`1=Ool001moo`8002moo`03 001oogoo04eoo`007goo0P00;goo00<007ooOol0CGoo000QOol5000ZOol00`00Oomoo`1=Ool002Io o`<002Moo`03001oogoo04eoo`00:Goo1@008Woo00<007ooOol0CGoo000^Ool8000JOol00`00Oomo o`1=Ool003Ioo`X000=oo`<000Aoo`8000Aoo`03001oogoo04eoo`00@7oo200000Eoo`00Oomoo`00 00=oo`03001oogoo04eoo`00A7oo00@007ooOomoo`/004eoo`00A7oo00<007ooOol00Woo00@007oo Ool000=oo`03001oogoo00h003moo`00A7oo00<007ooOol00Woo00@007ooOol000=oo`03001oogoo 00ioo`X003Eoo`00@goo0P001Goo0P0017oo00<007ooOol067oo1P00;goo001@Ool00`00Oomoo`0N Ool6000YOol0051oo`03001oogoo02Aoo`P0025oo`00D7oo00<007ooOol0;7oo0P007goo001@Ool0 0`00Oomoo`0^Ool2000MOol0051oo`03001oogoo02moo`8001aoo`00D7oo00<007ooOol0;7oo1000 7Goo001@Ool00`00Oomoo`0ZOol2000QOol0051oo`03001oogoo02Aoo`H002=oo`00D7oo00<007oo Ool077oo2000:Goo001@Ool00`00Oomoo`0FOol6000aOol0051oo`03001oogoo00ioo`P003Moo`00 D7oo00<007ooOol01Goo2@00?goo001?Ool90018Ool004Aoo`/00003Ool007oo04ioo`00=goo3@00 37oo00<007ooOol0CGoo000]Ool:000IOol00`00Oomoo`1=Ool002Ioo`L002=oo`03001oogoo04eo o`007goo1`00:Woo00<007ooOol0CGoo000HOol7000aOol00`00Oomoo`1=Ool0019oo`H003Qoo`03 001oogoo04eoo`003Woo1000?Woo00<007ooOol0CGoo000;Ool30012Ool00`00Oomoo`1=Ool000Mo o`@004Eoo`03001oogoo04eoo`001Woo00<007ooOol0Agoo00<007ooOol0CGoo0004Ool2001:Ool0 0`00Oomoo`1=Ool000=oo`03001oogoo04Yoo`03001oogoo04eoo`000goo00<007ooOol0BWoo00<0 07ooOol0CGoo0003Ool00`00Oomoo`1:Ool00`00Oomoo`1=Ool000Aoo`03001oogoo04Uoo`03001o ogoo04eoo`001Goo00<007ooOol0B7oo00<007ooOol0CGoo0006Ool30017Ool00`00Oomoo`1=Ool0 00Uoo`<004Aoo`03001oogoo04eoo`0037oo0P00@Woo00<007ooOol0CGoo000>Ool3000oOol00`00 Oomoo`1=Ool0015oo`8003eoo`03001oogoo04eoo`004goo0P00>goo00<007ooOol0CGoo000EOol3 000hOol00`00Oomoo`1=Ool001Qoo`8003Ioo`03001oogoo04eoo`006Woo0P00=7oo00<007ooOol0 CGoo000LOol3000aOol00`00Oomoo`1=Ool001moo`8002moo`03001oogoo04eoo`008Goo0P00;Goo 00<007ooOol0CGoo000SOol3000ZOol00`00Oomoo`1=Ool002Ioo`8002Qoo`03001oogoo04eoo`00 :7oo0P009Woo00<007ooOol0CGoo000ZOol3000SOol00`00Oomoo`1=Ool002eoo`80025oo`03001o ogoo04eoo`00;goo0P007goo00<007ooOol0CGoo000aOol2000MOol00`00Oomoo`1=Ool003=oo`80 01]oo`03001oogoo04eoo`00=Goo0P006Goo00<007ooOol0CGoo000gOol2000GOol00`00Oomoo`1= Ool003Uoo`03001oogoo01Aoo`03001oogoo04eoo`00>Woo00<007ooOol04goo00<007ooOol0CGoo 000kOol00`00Oomoo`0BOol00`00Oomoo`1=Ool003aoo`03001oogoo015oo`03001oogoo04eoo`00 ?Goo00<007ooOol047oo00<007ooOol0CGoo000nOol00`00Oomoo`0?Ool00`00Oomoo`1=Ool003mo o`8000moo`03001oogoo04eoo`00@Goo00<007ooOol037oo00<007ooOol0CGoo0012Ool00`00Oomo o`0;Ool00`00Oomoo`1=Ool004=oo`03001oogoo00Yoo`03001oogoo04eoo`00A7oo00<007ooOol0 2Goo00<007ooOol0CGoo0013Ool30003Ool30004Ool00`00Oomoo`1=Ool004Aoo`8000Ioo`05001o ogooOol0001?Ool004Aoo`03001oo`0000Eoo`05001oogooOol0001?Ool004Aoo`05001oogoo001o o`030004Ool00`00Oomoo`1=Ool004Aoo`06001oogoo001oo`001Woo00<007ooOol0CGoo0013Ool2 0003Ool50003Ool00`00Oomoo`1=Ool004Qoo`03001oogoo00Eoo`03001oogoo04eoo`00B7oo00<0 07ooOol01Goo00<007ooOol0CGoo0019Ool00`00Oomoo`04Ool00`00Oomoo`1=Ool004Uoo`03001o ogoo00Aoo`03001oogoo04eoo`00BGoo00<007ooOol017oo00<007ooOol0CGoo001:Ool00`00Oomo o`03Ool00`00Oomoo`1=Ool004Yoo`03001oogoo00=oo`03001oogoo04eoo`00Bgoo00<007ooOol0 0Woo00<007ooOol0CGoo001;Ool00`00Oomoo`02Ool00`00Oomoo`1=Ool004]oo`03001oogoo009o o`03001oogoo04eoo`00C7oo00D007ooOomoo`0004moo`00C7oo00D007ooOomoo`0004moo`00C7oo 00D007ooOomoo`0004moo`00CGoo00@007ooOol004moo`00CGoo00@007ooOol004moo`00CGoo00@0 07ooOol004moo`00CGoo00@007ooOol004moo`00CGoo00@007ooOol004moo`00CGoo00@007ooOol0 04moo`00CWoo00<007oo0000Cgoo001>Ool00`00Ool0001?Ool004ioo`03001oo`0004moo`00CWoo 00<007oo0000Cgoo001>Ool00`00Ool0001?Ool004ioo`03001oo`0004moo`00CWoo00<007oo0000 Cgoo001>Ool00`00Ool0001?Ool004ioo`03001oo`0004moo`00CWoo00<007oo0000Cgoo001>Ool0 0`00Ool0001?Ool004moo`8004moo`00Cgoo0P00Cgoo001?Ool2001?Ool004moo`8004moo`00Cgoo 0P00Cgoo001?Ool2001?Ool004moo`8004moo`00Cgoo0P00Cgoo001?Ool2001?Ool004moo`8004mo o`00Cgoo0P00Cgoo001?Ool2001?Ool004moo`8004moo`00Cgoo0P00Cgoo001?Ool2001?Ool004mo o`8004moo`00Cgoo0P00Cgoo001?Ool2001?Ool004moo`8004moo`00Cgoo0P00Cgoo001?Ool2001? Ool004moo`8004moo`00Cgoo0P00Cgoo001?Ool2001?Ool004moo`8004moo`00Cgoo0P00Cgoo001? Ool2001?Ool004moo`8004moo`00Cgoo0P00Cgoo001?Ool2001?Ool004moo`8004moo`00Cgoo0P00 Cgoo001?Ool2001?Ool004=oo`@000=oo`8000=oo`8004moo`00A7oo00<007ooOol00Woo00@007oo Ool0009oo`8004moo`00AGoo00D007ooOomoo`00009oo`05001oogooOol0001?Ool004Ioo`04001o ogoo0002Ool01@00Oomoogoo0000Cgoo0013Ool01000Oomoo`000Woo00@007ooOol000=oo`03001o ogoo04eoo`00A7oo0P0017oo0P0017oo00<007ooOol0CGoo001@Ool00`00Oomoo`1=Ool0051oo`03 001oogoo04eoo`00D7oo00<007ooOol0CGoo001@Ool00`00Oomoo`1=Ool00:1oo`00X7oo002POol0 0:1oo`00X7oo002POol00:1oo`00CWoo1000CWoo001?Ool00`00Oomoo`1>Ool004moo`8004moo`00 D7oo00<007ooOol0CGoo001>Ool4001>Ool00:1oo`00X7oo002POol00:1oo`00X7oo002POol00:1o o`00X7oo0000\ \>"], ImageRangeCache->{{{0, 159.062}, {319.125, 0}} -> {-1.29366, -1.10637, 0.0161695, 0.0722626}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(Display["\", wplot, "\"]; \)\)], "Input"], Cell[OutputFormData["\<\ Graphics[{{Line[{{8.334835016741514*^-10, 8.333333333333332*^-7}, {0.0004134977501991956, 0.3949634152604035}, {0.0009411449951950292, 0.8113398314583159}, {0.001672440551194225, 1.233631650861948}, {0.002226851029490147, 1.473093217150288}, {0.002871396498822547, 1.696175997187473}, {0.003664786450166736, 1.916581853576525}, {0.004596318840403102, 2.124933078177035}, {0.006979960336810846, 2.515029157291579}, {0.008664269523032089, 2.718633784923599}, {0.01088671187533187, 2.934477757769428}, {0.01637702322420161, 3.321724147916491}, {0.02016289630686067, 3.519357107081514}, {0.02521173981406187, 3.73195296053368}, {0.03773765833556398, 4.116260969062211}, {0.04713338030091423, 4.328457981888415}, {0.06005356421997989, 4.560080723571348}, {0.07564952812703202, 4.781349632136697}, {0.0932053627071011, 4.98197295406207}, {0.1428262013493113, 5.395107506387652}, {0.1800723591580786, 5.621887618121903}, {0.2234383103344827, 5.835526769582737}, {0.3382643478857791, 6.257164806817227}, {0.4846797409323519, 6.645091914686702}, {0.6898931431768789, 7.073951100393891}, {0.797371841088871, 7.2833788314849}, {0.88883749450492, 7.474329391248982}, {0.928135668141628, 7.568989026667469}, {0.9634698149527299, 7.671958433145478}, {0.976340319698375, 7.719111076337184}, {0.9872887462965479, 7.76920483168689}, {0.9941459144138963, 7.811971188215083}, {0.9969302318850424, 7.836318973919053}, {0.9987557779030927, 7.859176969691879}, {0.9998124372051343, 7.883272750205294}, {0.9998984717736009, 7.909630552626772}, {0.9987362610184389, 7.937148035700091}, {0.996451159776906, 7.962749621396223}, {0.9931999256104887, 7.986805309506316}, {0.9892821907295285, 8.008663058202401}, {0.9769892847067477, 8.057713171756148}, {0.9678852904516045, 8.085197715308453}, {0.9580191207754451, 8.110538781830728}, {0.9303017058727599, 8.167803822124302}, {0.8596607865000937, 8.271736366873755}, {0.7675704957830466, 8.369876673332537}, {0.6406806825573776, 8.475308258783279}, {0.3439289878564343, 8.665766217456551}, {-0.07966055305965115, 8.890672850286593}, {-0.4612775015790412, 9.099199990415019}, {-0.612525321406037, 9.202243744533017}, {-0.6742894609168032, 9.254978538521354}, {-0.7266674465240229, 9.311948270826066}, {-0.7493210336247063, 9.344629983131701}, {-0.7574531879602661, 9.358955410730397}, {-0.7650725629365469, 9.37465911650621}, {-0.770653899571625, 9.388524231683309}, {-0.7754190234482286, 9.40336817161185}, {-0.779073255108026, 9.419043415945452}, {-0.7811895049626348, 9.433658391546856}, {-0.7820403041035404, 9.447865440226848}, {-0.7816053014096446, 9.463299312905155}, {-0.7798862480375237, 9.477874615831856}, {-0.7771628749973563, 9.491341895661275}, {-0.7690158370234509, 9.516827374677932}, {-0.7560885952338868, 9.543866864659417}, {-0.7394835932981024, 9.569917763956603}, {-0.7207451189058073, 9.59377181442475}, {-0.6664769527997432, 9.64749881835487}, {-0.5289071978040241, 9.744661712008253}, {-0.3268093712007861, 9.85255956965075}, {-0.07801517081697126, 9.966008276872916}, {0.3727237572651019, 10.16698558325906}, {0.5800753098439503, 10.27735340444549}, {0.6625185994897168, 10.33152731907228}, {0.7248504534406868, 10.38141051616153}, {0.7707864046096601, 10.42863792020772}, {0.7888738216867072, 10.4524327612051}, {0.8047332700393157, 10.47857265444744}, {0.8163707706474051, 10.50480642435676}, {0.8210089274180994, 10.51952300243286}, {0.8240524857241544, 10.53295211895978}, {0.8267237863253839, 10.55903372440318}, {0.8253226051479286, 10.58376154063997}, {0.8208523373304277, 10.60588511086966}, {0.8128110363494316, 10.6294077581851}, {0.80330792344471, 10.6494895009833}, {0.7900418829119604, 10.67167393315626}, {0.7504795817275679, 10.72101861053452}, {0.7013722040590649, 10.76738170843477}, {0.5633671485374494, 10.86578108474501}, {0.3704753532685561, 10.97350555604479}, {-0.08789599172526864, 11.19069565670389}, {-0.5051221282789473, 11.39042696081777}, {-0.8312144985163339, 11.57913945645481}, {-0.9700238816735735, 11.68246354663047}, {-1.086276784537598, 11.79456218265954}, {-1.126336284327286, 11.84403067284447}, {-1.16153473552512, 11.89644650918065}, {-1.188002798519181, 11.94555426192358}, {-1.206945937741074, 11.99046156839616}, {-1.222961124429312, 12.04255641882053}, {-1.229196262135353, 12.07127819521466}, {-1.233496482193146, 12.09856720776512}, {-1.236111546506153, 12.1236767761238}, {-1.237433404407047, 12.14697662898429}, {-1.237709858797188, 12.17254637937265}, {-1.23671145741426, 12.19979112902364}, {-1.235565132271463, 12.21537686562873}, {-1.234154732732037, 12.22970927036655}, {-1.22976748475605, 12.2620563572804}, {-1.224585203719995, 12.29057269862364}, {-1.217860356494308, 12.32068408131197}, {-1.183845248244461, 12.43024186739148}, {-1.140187907371456, 12.53448297716591}, {-1.09337114531561, 12.62947936092251}, {-0.9736154821836585, 12.84295597221756}, {-0.7475787686331773, 13.22806759291131}, {-0.5265111934933564, 13.66302182566885}, {-0.3704809953721853, 14.06523798134776}, {-0.2597694430185631, 14.45583077890885}, {-0.2159829270000973, 14.65547186823866}, {-0.1767491833782854, 14.87057786884671}, {-0.1446850287587797, 15.08410481777934}, {-0.120345409272845, 15.27978457269191}, {-0.08372236632048631, 15.66375958686647}, {-0.06887013463682958, 15.86983433295298}, {-0.05741930142655284, 16.06146473530071}, {-0.04642474283614166, 16.28518554809921}, {-0.03797915130739642, 16.49628596780286}, {-0.03120324704395413, 16.70263322229755}, {-0.02521064076509822, 16.92620920997048}, {-0.01734962828988923, 17.31665417926005}, {-0.01432686026056343, 17.51561316522503}, {-0.01159244578186552, 17.73435832739694}, {-0.00924775698533331, 17.96538219334312}, {-0.007499008345086533, 18.17636705363388}, {-0.004915881569010654, 18.58507256213444}, {-0.0029599429307068, 19.02392175852851}, {-0.001624387964575414, 19.42371732004379}, {-0.0004066997330600284, 19.85931829394395}, {-0.0000297303929012573, 19.99999916666666}}]}}, {PlotRange -> Automatic, AspectRatio -> 2, DisplayFunction :> $DisplayFunction, ColorOutput -> Automatic, Axes -> Automatic, AxesOrigin -> Automatic, PlotLabel -> None, AxesLabel -> {\"w\", \"z\"}, Ticks -> Automatic, GridLines -> None, Prolog -> {}, Epilog -> {}, AxesStyle -> Automatic, Background -> Automatic, DefaultColor -> Automatic, DefaultFont :> $DefaultFont, RotateLabel -> True, Frame -> False, FrameStyle -> Automatic, FrameTicks -> Automatic, FrameLabel -> None, PlotRegion -> Automatic, ImageSize -> Automatic, TextStyle :> $TextStyle, FormatType :> $FormatType}]\ \>", "\<\ -Graphics-\ \>"], "Output"] }, Open ]] }, Closed]] }, Open ]] }, FrontEndVersion->"Macintosh 3.0", ScreenRectangle->{{0, 1152}, {0, 850}}, AutoGeneratedPackage->None, WindowSize->{519, 675}, WindowMargins->{{Automatic, 98}, {Automatic, 7}}, StyleDefinitions -> "ok.nb", MacintoshSystemPageSetup->"\<\ 00<0001804P000000]P2:?oQon82n@960dL5:0?l0080001804P000000]P2:001 0000I00000400`<300000BL?00400@00000000000000060801T1T000000@0000 00000000004000000000000000000000\>" ] (*********************************************************************** Cached data follows. If you edit this Notebook file directly, not using Mathematica, you must remove the line containing CacheID at the top of the file. The cache data will then be recreated when you save this file from within Mathematica. ***********************************************************************) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[1731, 51, 90, 4, 78, "Title"], Cell[1824, 57, 1069, 30, 230, "Text"], Cell[CellGroupData[{ Cell[2918, 91, 76, 1, 63, "Section", InitializationCell->True], Cell[2997, 94, 291, 7, 137, "Input", InitializationCell->True], Cell[3291, 103, 129, 3, 47, "Input", InitializationCell->True] }, Closed]], Cell[CellGroupData[{ Cell[3457, 111, 83, 1, 37, "Section", InitializationCell->True], Cell[3543, 114, 75, 2, 47, "Input", InitializationCell->True], Cell[3621, 118, 75, 2, 47, "Input", InitializationCell->True] }, Closed]], Cell[CellGroupData[{ Cell[3733, 125, 48, 0, 37, "Section"], Cell[3784, 127, 54, 1, 47, "Input"], Cell[CellGroupData[{ Cell[3863, 132, 57, 1, 47, "Input"], Cell[3923, 135, 225, 6, 52, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[4185, 146, 122, 2, 65, "Input"], Cell[4310, 150, 16430, 547, 304, 6292, 417, "GraphicsData", "PostScript", "Graphics"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[20789, 703, 50, 0, 37, "Section"], Cell[20842, 705, 94, 3, 53, "Text"], Cell[CellGroupData[{ Cell[20961, 712, 90, 1, 47, "Input"], Cell[21054, 715, 15569, 433, 194, 4165, 288, "GraphicsData", "PostScript", "Graphics"] }, Open ]], Cell[36638, 1151, 69, 0, 53, "Text"], Cell[CellGroupData[{ Cell[36732, 1155, 100, 2, 47, "Input"], Cell[36835, 1159, 76, 1, 52, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[36948, 1165, 78, 1, 47, "Input"], Cell[37029, 1168, 225, 6, 52, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[37291, 1179, 122, 2, 65, "Input"], Cell[37416, 1183, 16488, 542, 304, 6168, 410, "GraphicsData", "PostScript", "Graphics"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[53953, 1731, 32, 0, 37, "Section"], Cell[53988, 1733, 79, 1, 47, "Input"], Cell[CellGroupData[{ Cell[54092, 1738, 209, 4, 101, "Input"], Cell[54304, 1744, 14263, 414, 304, 5491, 301, "GraphicsData", "PostScript", "Graphics"] }, Open ]], Cell[CellGroupData[{ Cell[68604, 2163, 98, 2, 65, "Input"], Cell[68705, 2167, 4225, 87, 50, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[72967, 2259, 49, 1, 47, "Input"], Cell[73019, 2262, 59, 1, 52, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[73115, 2268, 90, 1, 47, "Input"], Cell[73208, 2271, 20073, 614, 194, 6838, 446, "GraphicsData", "PostScript", "Graphics"] }, Open ]], Cell[CellGroupData[{ Cell[93318, 2890, 112, 3, 65, "Input"], Cell[93433, 2895, 162, 3, 45, "Message"], Cell[93598, 2900, 76, 1, 52, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[93711, 2906, 56, 1, 47, "Input"], Cell[93770, 2909, 57, 1, 52, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[93864, 2915, 55, 1, 47, "Input"], Cell[93922, 2918, 225, 6, 52, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[94184, 2929, 103, 2, 65, "Input"], Cell[94290, 2933, 54, 1, 52, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[94381, 2939, 185, 3, 101, "Input"], Cell[94569, 2944, 17140, 566, 337, 7562, 443, "GraphicsData", "PostScript", "Graphics"] }, Open ]], Cell[CellGroupData[{ Cell[111746, 3515, 81, 1, 47, "Input"], Cell[111830, 3518, 7618, 159, 50, "Output"] }, Open ]] }, Closed]] }, Open ]] } ] *) (*********************************************************************** End of Mathematica Notebook file. ***********************************************************************)