(*********************************************************************** 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[ 64388, 1774]*) (*NotebookOutlinePosition[ 65285, 1803]*) (* CellTagsIndexPosition[ 65241, 1799]*) (*WindowFrame->Normal*) Notebook[{ Cell[CellGroupData[{ Cell[TextData[ "Shooting schemes\nfor the shallow water equations\non a \[Beta]-plane"], "Title", TextAlignment->Center], Cell[CellGroupData[{ Cell["Functions for solving the ODE", "Section", InitializationCell->True], Cell[BoxData[ \(soln[\[Omega]_, k_, LR_, Lf_, L_] := NDSolve[{\n\t\t\t \(\(v'\)'\)[y] + \((\(-\((1 + y/Lf)\)^2\) - k^2 + \[Omega]^2 - \((k/\[Omega])\)*\((LR/Lf)\))\) v[y] == 0, \n\t\t v[0] == 0, \n\t\t\t\(v'\)[0] == 1}, \n\t\tv, {y, 0, L}]\)], "Input",\ InitializationCell->True], Cell[BoxData[ \(vAtL[\[Omega]_, k_, LR_, Lf_, L_] := v[L] /. \(soln[\[Omega], k, LR, Lf, L]\)[\([1]\)]\)], "Input", InitializationCell->True] }, Closed]], Cell[CellGroupData[{ Cell["Some parameters we don't change much", "Section", InitializationCell->True], Cell[BoxData[ \(\({k, LR, Lf, L} = { .1, 1, 10, 1}; \)\)], "Input", InitializationCell->True] }, Closed]], Cell[CellGroupData[{ Cell["Try it!", "Section"], Cell[CellGroupData[{ Cell[BoxData[ \(trysoln = soln[\(- .001\), k, LR, Lf, L]\)], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{"v", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.`, 1.`}}, "<>"]\), False, Editable->False]}], "}"}], "}"}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(ParametricPlot[{v[y] /. trysoln[\([1]\)], y}, {y, 0, 1}, AspectRatio -> 2, Ticks -> {{}, Automatic}]; \)\)], "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 2.84708 0.047619 1.90476 [ [.01131 .42857 -18 -4.5 ] [.01131 .42857 0 4.5 ] [.01131 .80952 -18 -4.5 ] [.01131 .80952 0 4.5 ] [.01131 1.19048 -18 -4.5 ] [.01131 1.19048 0 4.5 ] [.01131 1.57143 -18 -4.5 ] [.01131 1.57143 0 4.5 ] [.01131 1.95238 -6 -4.5 ] [.01131 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 0 .04762 m 1 .04762 L s .02381 .42857 m .03006 .42857 L s [(0.2)] .01131 .42857 1 0 Mshowa .02381 .80952 m .03006 .80952 L s [(0.4)] .01131 .80952 1 0 Mshowa .02381 1.19048 m .03006 1.19048 L s [(0.6)] .01131 1.19048 1 0 Mshowa .02381 1.57143 m .03006 1.57143 L s [(0.8)] .01131 1.57143 1 0 Mshowa .02381 1.95238 m .03006 1.95238 L s [(1)] .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 .52381 m .02756 .52381 L s .02381 .61905 m .02756 .61905 L s .02381 .71429 m .02756 .71429 L s .02381 .90476 m .02756 .90476 L s .02381 1 m .02756 1 L s .02381 1.09524 m .02756 1.09524 L s .02381 1.28571 m .02756 1.28571 L s .02381 1.38095 m .02756 1.38095 L s .02381 1.47619 m .02756 1.47619 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 0 0 m 1 0 L 1 2 L 0 2 L closepath clip newpath .5 Mabswid .02381 .04762 m .13902 .12489 L .26268 .20916 L .37503 .2883 L .47798 .36442 L .58041 .44543 L .67052 .52342 L .75566 .60631 L .82601 .68617 L .88185 .763 L .90542 .80189 L .92759 .84473 L .94503 .88583 L .95765 .92343 L .96728 .96225 L .97111 .98366 L .97262 .99433 L .97376 1.00401 L .97464 1.01322 L .97537 1.02318 L .97568 1.02866 L .97589 1.0337 L .97605 1.03853 L .9761 1.04095 L .97613 1.04223 L .97615 1.04362 L .97618 1.04608 L .97619 1.04842 L .97619 1.05069 L .97618 1.0528 L .97616 1.05507 L .97612 1.05757 L .97606 1.06018 L .976 1.06264 L .97584 1.0671 L .97565 1.07118 L .97509 1.08039 L .97431 1.08987 L .97325 1.1 L .97072 1.11824 L .96672 1.13993 L .96215 1.15964 L .95068 1.19801 L .9345 1.23961 L .8915 1.32144 L .8367 1.40025 L .77241 1.47602 L .69249 1.5567 L .60557 1.63435 L .50386 1.71691 L Mistroke .39836 1.79643 L .29145 1.87292 L .17639 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]HGAYHf4PAg9QL6QYHgGoo00<007ooOol0?Woo000COol00`00Oomoo`0jOol00`00Oomo o`0mOol001=oo`03001oogoo03]oo`03001oogoo03aoo`004goo00<007ooOol0?7oo00<007ooOol0 >goo000COol00`00Oomoo`0mOol2000kOol0009oo`8000Aoo`8000=oo`@0009oo`03001oogoo03mo o`03001oogoo03Qoo`0000Eoo`00Oomoo`0000=oo`8000Aoo`03001oogoo009oo`03001oogoo041o o`03001oogoo03Moo`0000Eoo`00Oomoo`0000Yoo`04001oogooOol20012Ool00`00Oomoo`0fOol0 0005Ool007ooOol0000;Ool01000Oomoo`00A7oo00<007ooOol0=Goo00001Goo001oogoo000027oo 00@007ooOol0009oo`03001oogoo04=oo`8003Eoo`000Woo0P002Woo0P000goo00<007ooOol0AGoo 00<007ooOol0Ool001=oo`03001oogoo06Qoo`03001oogoo00moo`004goo 00<007ooOol0Igoo00<007ooOol047oo000COol2001WOol00`00Oomoo`0AOol001=oo`03001oogoo 06Ioo`03001oogoo015oo`004goo00<007ooOol0IGoo00<007ooOol04Woo000COol00`00Oomoo`1T Ool00`00Oomoo`0COol001=oo`03001oogoo06Aoo`03001oogoo01=oo`004goo00<007ooOol0Hgoo 00<007ooOol057oo000COol00`00Oomoo`1ROol00`00Oomoo`0EOol001=oo`03001oogoo065oo`03 001oogoo01Ioo`004goo00<007ooOol0H7oo00<007ooOol05goo000COol00`00Oomoo`1OOol00`00 Oomoo`0HOol001=oo`03001oogoo05ioo`03001oogoo01Uoo`004goo00<007ooOol0GGoo00<007oo Ool06Woo000COol2001MOol00`00Oomoo`0KOol001=oo`03001oogoo05]oo`03001oogoo01aoo`00 4goo00<007ooOol0FWoo00<007ooOol07Goo000COol00`00Oomoo`1IOol00`00Oomoo`0NOol001=o o`03001oogoo05Qoo`03001oogoo01moo`004goo00<007ooOol0Egoo00<007ooOol087oo000COol0 0`00Oomoo`1FOol00`00Oomoo`0QOol001=oo`03001oogoo05Eoo`03001oogoo029oo`004goo00<0 07ooOol0E7oo00<007ooOol08goo000COol00`00Oomoo`1COol00`00Oomoo`0TOol0009oo`8000Ao o`8000Aoo`8000=oo`03001oogoo055oo`8002Moo`0000Eoo`00Oomoo`0000=oo`8000=oo`04001o ogoo0002Ool00`00Oomoo`1@Ool00`00Oomoo`0WOol00005Ool007ooOol00008Ool01000Oomoo`00 0Woo0P00D7oo00<007ooOol0:7oo00001Goo001oogoo00002Goo0P000goo00<007ooOol0CWoo00<0 07ooOol0:Goo00001Goo001oogoo000027oo00@007ooOol0009oo`03001oogoo04eoo`03001oogoo 02Yoo`000Woo0P002Woo0P000goo00<007ooOol0C7oo00<007ooOol0:goo000COol00`00Oomoo`1; Ool00`00Oomoo`0/Ool001=oo`03001oogoo04Yoo`03001oogoo02eoo`004goo00<007ooOol0BGoo 00<007ooOol0;Woo000COol00`00Oomoo`18Ool00`00Oomoo`0_Ool001=oo`03001oogoo04Ioo`80 039oo`004goo00<007ooOol0AGoo00<007ooOol0Goo000COol00`00Oomoo`0mOol00`00Oomoo`0jOol001=oo`03001oogoo03ao o`03001oogoo03]oo`004goo00<007ooOol0>goo00<007ooOol0?7oo000COol00`00Oomoo`0iOol2 000oOol001=oo`03001oogoo03Qoo`03001oogoo03moo`004goo00<007ooOol0=Woo0P00@Woo000C Ool00`00Oomoo`0eOol00`00Oomoo`12Ool001=oo`03001oogoo03Aoo`03001oogoo04=oo`004goo 0P00Ool00`00Oomoo`02Ool00`00Oomoo`1jOol000ioo`03001oogoo009oo`03001oogoo 07Yoo`003Goo0P0017oo00<007ooOol0NWoo000COol00`00Oomoo`1jOol001=oo`03001oogoo07Yo o`004goo00<007ooOol0NWoo002@Ool0091oo`00T7oo002@Ool0091oo`00T7oo002@Ool0091oo`00 T7oo002@Ool0091oo`00T7oo002@Ool0091oo`00T7oo002@Ool0091oo`00\ \>"], ImageRangeCache->{{{0, 143}, {287, 0}} -> {-0.0543986, -0.0958895, 0.00277814, 0.00415254}}] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["Try shooting for v=0 at y=L:", "Section"], Cell["We see there are some high frequency solutions:", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(\(Plot[vAtL[\[Omega], k, LR, Lf, L], {\[Omega], \(-8\), \(- .001\)}]; \)\)], "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.97619 0.119048 0.10624 0.421305 [ [.02381 .09374 -6 -9 ] [.02381 .09374 6 0 ] [.2619 .09374 -6 -9 ] [.2619 .09374 6 0 ] [.5 .09374 -6 -9 ] [.5 .09374 6 0 ] [.7381 .09374 -6 -9 ] [.7381 .09374 6 0 ] [.96369 .02198 -24 -4.5 ] [.96369 .02198 0 4.5 ] [.96369 .1905 -18 -4.5 ] [.96369 .1905 0 4.5 ] [.96369 .27476 -18 -4.5 ] [.96369 .27476 0 4.5 ] [.96369 .35902 -18 -4.5 ] [.96369 .35902 0 4.5 ] [.96369 .44328 -18 -4.5 ] [.96369 .44328 0 4.5 ] [.96369 .52754 -6 -4.5 ] [.96369 .52754 0 4.5 ] [.96369 .61181 -18 -4.5 ] [.96369 .61181 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 .02381 .10624 m .02381 .11249 L s [(-8)] .02381 .09374 0 1 Mshowa .2619 .10624 m .2619 .11249 L s [(-6)] .2619 .09374 0 1 Mshowa .5 .10624 m .5 .11249 L s [(-4)] .5 .09374 0 1 Mshowa .7381 .10624 m .7381 .11249 L s [(-2)] .7381 .09374 0 1 Mshowa .125 Mabswid .08333 .10624 m .08333 .10999 L s .14286 .10624 m .14286 .10999 L s .20238 .10624 m .20238 .10999 L s .32143 .10624 m .32143 .10999 L s .38095 .10624 m .38095 .10999 L s .44048 .10624 m .44048 .10999 L s .55952 .10624 m .55952 .10999 L s .61905 .10624 m .61905 .10999 L s .67857 .10624 m .67857 .10999 L s .79762 .10624 m .79762 .10999 L s .85714 .10624 m .85714 .10999 L s .91667 .10624 m .91667 .10999 L s .25 Mabswid 0 .10624 m 1 .10624 L s .97619 .02198 m .98244 .02198 L s [(-0.2)] .96369 .02198 1 0 Mshowa .97619 .1905 m .98244 .1905 L s [(0.2)] .96369 .1905 1 0 Mshowa .97619 .27476 m .98244 .27476 L s [(0.4)] .96369 .27476 1 0 Mshowa .97619 .35902 m .98244 .35902 L s [(0.6)] .96369 .35902 1 0 Mshowa .97619 .44328 m .98244 .44328 L s [(0.8)] .96369 .44328 1 0 Mshowa .97619 .52754 m .98244 .52754 L s [(1)] .96369 .52754 1 0 Mshowa .97619 .61181 m .98244 .61181 L s [(1.2)] .96369 .61181 1 0 Mshowa .125 Mabswid .97619 .04304 m .97994 .04304 L s .97619 .06411 m .97994 .06411 L s .97619 .08517 m .97994 .08517 L s .97619 .12731 m .97994 .12731 L s .97619 .14837 m .97994 .14837 L s .97619 .16944 m .97994 .16944 L s .97619 .21157 m .97994 .21157 L s .97619 .23263 m .97994 .23263 L s .97619 .2537 m .97994 .2537 L s .97619 .29583 m .97994 .29583 L s .97619 .31689 m .97994 .31689 L s .97619 .33796 m .97994 .33796 L s .97619 .38009 m .97994 .38009 L s .97619 .40115 m .97994 .40115 L s .97619 .42222 m .97994 .42222 L s .97619 .46435 m .97994 .46435 L s .97619 .48541 m .97994 .48541 L s .97619 .50648 m .97994 .50648 L s .97619 .54861 m .97994 .54861 L s .97619 .56967 m .97994 .56967 L s .97619 .59074 m .97994 .59074 L s .97619 .00091 m .97994 .00091 L s .25 Mabswid .97619 0 m .97619 .61803 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath .5 Mabswid .02381 .15921 m .02846 .1596 L .03278 .15988 L .03753 .16011 L .04019 .16021 L .04261 .16027 L .04385 .16029 L .04516 .16031 L .04627 .16032 L .0475 .16032 L .04882 .16032 L .04956 .16032 L .05024 .16032 L .05156 .1603 L .05281 .16028 L .05396 .16026 L .05519 .16022 L .05741 .16015 L .06005 .16004 L .06244 .15992 L .06787 .15955 L .07286 .1591 L .08413 .15773 L .09407 .1561 L .10457 .15395 L .12421 .14875 L .14254 .14257 L .18144 .12567 L .2228 .10333 L .26264 .07957 L .30096 .05696 L .34174 .03592 L .36199 .02754 L .37191 .02409 L .381 .02135 L .38967 .01914 L .39916 .0172 L .40378 .01645 L .4087 .01579 L .41143 .0155 L .41391 .01527 L .41641 .01507 L .41875 .01493 L .42105 .01482 L .42228 .01478 L .42359 .01475 L .42425 .01474 L .42497 .01473 L .42628 .01472 L .42755 .01472 L Mistroke .42873 .01473 L .42993 .01475 L .43107 .01478 L .43208 .01481 L .43318 .01486 L .43566 .01499 L .43793 .01515 L .44276 .01561 L .44805 .0163 L .45357 .01725 L .45872 .01834 L .4692 .02117 L .47911 .02464 L .49768 .03321 L .51763 .0455 L .53606 .0597 L .57537 .09883 L .61317 .14691 L .65342 .20738 L .69215 .27189 L .73334 .34359 L .77301 .41172 L .81117 .47255 L .85177 .52833 L .89087 .56987 L .91225 .58646 L .92205 .59243 L .93241 .59754 L .93702 .59937 L .94197 .60101 L .9448 .60179 L .94739 .60238 L .94864 .60262 L .95 .60284 L .95127 .60302 L .95246 .60315 L .95352 .60324 L .95467 .6033 L .95594 .60332 L .95711 .60329 L .95826 .60321 L .95951 .60306 L .96083 .60282 L .96205 .6025 L .9632 .60209 L .9643 .60158 L .96528 .60101 L .96634 .60022 L .9675 .59909 L .96872 .59746 L Mistroke .96977 .59554 L .97093 .59249 L .97222 .58692 L .97341 .57732 L .97406 .56755 L .97477 .547 L .97607 .1289 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]HGAYHf4PAg9QL6QYHggoo00<007ooOol0Ngoo00<007ooOol027oo001HOol2000nOol00`00Oomo o`1jOol00`00Oomoo`08Ool005Ioo`80045oo`03001oogoo07Uoo`03001oogoo00Qoo`00E7oo0P00 A7oo00<007ooOol0N7oo00<007ooOol027oo001BOol20017Ool00`00Oomoo`1gOol20009Ool000]o o`80045oo`800003Ool007oo03ioo`<000Qoo`03001oogoo03Aoo`@003ioo`03001oogoo00Qoo`00 2Woo00@007ooOol003moo`04001oo`000011Ool00`00Oomoo`08Ool00`00Oomoo`0dOol00`00Oomo o`0nOol00`00Oomoo`08Ool000Aoo`D00005Ool007ooOol0000iOol500001Goo000007oo0000>Goo 1@0000=oo`0000000P002Woo00<007ooOol0;7oo1@000goo00<007ooOol0?Goo00<007ooOol027oo 000;Ool2000nOol50010Ool00`00Ool0000goo0P0017oo0`00@7oo0P003Woo00<007oo Ool0<7oo0P00?goo0P002Goo0016Ool2001JOol00`00Oomoo`1`Ool00`00Oomoo`08Ool004Eoo`03 001oogoo05]oo`03001oogoo06moo`03001oogoo00Qoo`00@goo0P00Ggoo00<007ooOol0KWoo00<0 07ooOol027oo0012Ool00`00Oomoo`1OOol00`00Oomoo`1^Ool00`00Oomoo`08Ool0041oo`8006=o o`03001oogoo06eoo`03001oogoo00Qoo`000Wooo`0070000goo0009Ool00`00Oomoo`0>Ool00`00 Oomoo`0=Ool00`00Oomoo`0>Ool3000>Ool00`00Oomoo`0>Ool00`00Oomoo`0=Ool00`00Oomoo`0> Ool00`00Oomoo`0>Ool00`00Oomoo`0>Ool00`00Oomoo`04Ool00`00Oomoo`07Ool00`00Oomoo`0= Ool00`00Oomoo`0>Ool00`00Oomoo`0>Ool00`00Oomoo`0>Ool00`00Oomoo`0=Ool00`00Oomoo`0> Ool00`00Oomoo`08Ool003Yoo`8006aoo`03001oogoo06Yoo`03001oogoo00Qoo`00>7oo0P00KWoo 00<007ooOol0JWoo00<007ooOol027oo000fOol2001aOol00`00Oomoo`1YOol00`00Oomoo`08Ool0 03Aoo`8007Aoo`03001oogoo06Qoo`8000Uoo`00Ool00`00Oomo o`14Ool00`00Oomoo`08Ool00goo00<007ooOol027oo003HOol00`00Oomoo`0jOol00`00Oomoo`08Ool00=Qoo`03 001oogoo03Yoo`03001oogoo00Qoo`00fGoo00<007ooOol0>Goo00<007ooOol027oo003JOol00`00 Oomoo`0hOol00`00Oomoo`08Ool00=Yoo`03001oogoo03Qoo`8000Uoo`00fgoo00<007ooOol0=goo 00<007ooOol027oo003KOol00`00Oomoo`0gOol00`00Oomoo`08Ool00=aoo`03001oogoo03Ioo`03 001oogoo00Qoo`00g7oo00<007ooOol0=Woo00<007ooOol027oo003MOol00`00Oomoo`0eOol00`00 Oomoo`08Ool00=ioo`03001oogoo03Aoo`8000Uoo`00gWoo00<007ooOol0=7oo00<007ooOol027oo 003OOol00`00Oomoo`0cOol00`00Oomoo`08Ool00>1oo`03001oogoo01moo`8000Aoo`8000Aoo`80 00Eoo`03001oogoo00Qoo`00h7oo00<007ooOol07Woo00@007ooOol000=oo`8000=oo`04001oogoo 0004Ool00`00Oomoo`08Ool00>5oo`03001oogoo01eoo`04001oogoo0008Ool01000Oomoo`0017oo 0`0027oo003ROol00`00Oomoo`0LOol01000Oomoo`002Goo0P001Goo00<007ooOol027oo003ROol0 0`00Oomoo`0LOol01000Oomoo`0027oo00@007ooOol000Aoo`03001oogoo00Qoo`00hgoo00<007oo Ool077oo0P002Woo0P001Goo00<007ooOol027oo003SOol00`00Oomoo`0_Ool00`00Oomoo`08Ool0 0>Aoo`03001oogoo02ioo`03001oogoo00Qoo`00iGoo00<007ooOol0;Goo0P002Goo003UOol00`00 Oomoo`0]Ool00`00Oomoo`08Ool00>Ioo`03001oogoo02aoo`03001oogoo00Qoo`00igoo00<007oo Ool0:goo00<007ooOol027oo003WOol00`00Oomoo`0[Ool00`00Oomoo`08Ool00>Qoo`03001oogoo 02Yoo`03001oogoo00Qoo`00jGoo00<007ooOol0:Goo0P002Goo003ZOol00`00Oomoo`0XOol00`00 Oomoo`08Ool00>Yoo`03001oogoo02Qoo`03001oogoo00Qoo`00jgoo00<007ooOol09goo00<007oo Ool027oo003/Ool00`00Oomoo`0VOol00`00Oomoo`08Ool00>eoo`03001oogoo02Eoo`03001oogoo 00Qoo`00kGoo00<007ooOol09Goo0P002Goo003^Ool00`00Oomoo`0TOol00`00Oomoo`08Ool00>mo o`03001oogoo02=oo`03001oogoo00Qoo`00kgoo00<007ooOol08goo00<007ooOol027oo003`Ool0 0`00Oomoo`0JOol30005Ool00`00Oomoo`08Ool00?5oo`03001oogoo01Yoo`03001oogoo00Aoo`03 001oogoo00Qoo`00lWoo00<007ooOol06Goo00<007ooOol017oo0`0027oo003bOol00`00Oomoo`0I Ool00`00Oomoo`04Ool00`00Oomoo`08Ool00?=oo`03001oogoo01Qoo`03001oogoo00Aoo`03001o ogoo00Qoo`00m7oo00<007ooOol05Woo0P001Woo00<007ooOol027oo003eOol00`00Oomoo`0MOol0 0`00Oomoo`08Ool00?Ioo`03001oogoo01aoo`03001oogoo00Qoo`00mgoo00<007ooOol06goo0P00 2Goo003hOol00`00Oomoo`0JOol00`00Oomoo`08Ool00?Uoo`03001oogoo01Uoo`03001oogoo00Qo o`00nWoo00<007ooOol067oo00<007ooOol027oo003kOol00`00Oomoo`0GOol00`00Oomoo`08Ool0 0?aoo`03001oogoo01Ioo`03001oogoo00Qoo`00oGoo0P005Woo0P002Goo003oOol00`00Oomoo`0C Ool00`00Oomoo`08Ool00?moo`5oo`03001oogoo015oo`8000Yoo`00ogoo0Woo00<007ooOol047oo 0P002Woo003oOol3Ool00`00Oomoo`0?Ool2000:Ool00?moo`Aoo`8000moo`8000Yoo`00ogoo1Woo 00<007ooOol037oo0`002Goo003oOol7Ool3000:Ool3000:Ool00?moo`Yoo`<000Moo`03001oo`00 00Yoo`00ogoo3Goo0`000Woo0P000Woo00<007ooOol027oo003oOol2Ool30004Ool20003Ool50003 Ool00`00Oomoo`08Ool00?moo`=oo`03001oogoo00=oo`8000Aoo`03001oogoo00Aoo`03001oogoo 00Qoo`00ogoo0goo00<007ooOol02Woo00<007ooOol00goo0`0027oo003oOol3Ool00`00Oomoo`0; Ool00`00Oomoo`02Ool00`00Oomoo`08Ool00?moo`=oo`03001oogoo00Qoo`04001oogoo0004Ool0 0`00Oomoo`08Ool00?moo`9oo`8000]oo`80011oo`00ogoo8Goo0000\ \>"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {-8.27678, -0.252176, 0.0298035, 0.00842156}}] }, Open ]], Cell["\<\ ...but also there are some low frequency, quasi-geostrophic \ solutions:\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(\(Plot[ vAtL[\[Omega], k, LR, Lf, L], {\[Omega], \(- .001\), \(- .0001\)}]; \)\)], "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.97619 952.381 0.384698 1.70309 [ [.02381 .3722 -18 -9 ] [.02381 .3722 18 0 ] [.21429 .3722 -21 -9 ] [.21429 .3722 21 0 ] [.40476 .3722 -21 -9 ] [.40476 .3722 21 0 ] [.59524 .3722 -21 -9 ] [.59524 .3722 21 0 ] [.78571 .3722 -21 -9 ] [.78571 .3722 21 0 ] [.96369 .04408 -24 -4.5 ] [.96369 .04408 0 4.5 ] [.96369 .12923 -30 -4.5 ] [.96369 .12923 0 4.5 ] [.96369 .21439 -24 -4.5 ] [.96369 .21439 0 4.5 ] [.96369 .29954 -30 -4.5 ] [.96369 .29954 0 4.5 ] [.96369 .46985 -24 -4.5 ] [.96369 .46985 0 4.5 ] [.96369 .55501 -18 -4.5 ] [.96369 .55501 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 .02381 .3847 m .02381 .39095 L s [(-0.001)] .02381 .3722 0 1 Mshowa .21429 .3847 m .21429 .39095 L s [(-0.0008)] .21429 .3722 0 1 Mshowa .40476 .3847 m .40476 .39095 L s [(-0.0006)] .40476 .3722 0 1 Mshowa .59524 .3847 m .59524 .39095 L s [(-0.0004)] .59524 .3722 0 1 Mshowa .78571 .3847 m .78571 .39095 L s [(-0.0002)] .78571 .3722 0 1 Mshowa .125 Mabswid .07143 .3847 m .07143 .38845 L s .11905 .3847 m .11905 .38845 L s .16667 .3847 m .16667 .38845 L s .2619 .3847 m .2619 .38845 L s .30952 .3847 m .30952 .38845 L s .35714 .3847 m .35714 .38845 L s .45238 .3847 m .45238 .38845 L s .5 .3847 m .5 .38845 L s .54762 .3847 m .54762 .38845 L s .64286 .3847 m .64286 .38845 L s .69048 .3847 m .69048 .38845 L s .7381 .3847 m .7381 .38845 L s .83333 .3847 m .83333 .38845 L s .88095 .3847 m .88095 .38845 L s .92857 .3847 m .92857 .38845 L s .25 Mabswid 0 .3847 m 1 .3847 L s .97619 .04408 m .98244 .04408 L s [(-0.2)] .96369 .04408 1 0 Mshowa .97619 .12923 m .98244 .12923 L s [(-0.15)] .96369 .12923 1 0 Mshowa .97619 .21439 m .98244 .21439 L s [(-0.1)] .96369 .21439 1 0 Mshowa .97619 .29954 m .98244 .29954 L s [(-0.05)] .96369 .29954 1 0 Mshowa .97619 .46985 m .98244 .46985 L s [(0.05)] .96369 .46985 1 0 Mshowa .97619 .55501 m .98244 .55501 L s [(0.1)] .96369 .55501 1 0 Mshowa .125 Mabswid .97619 .06111 m .97994 .06111 L s .97619 .07814 m .97994 .07814 L s .97619 .09517 m .97994 .09517 L s .97619 .1122 m .97994 .1122 L s .97619 .14627 m .97994 .14627 L s .97619 .1633 m .97994 .1633 L s .97619 .18033 m .97994 .18033 L s .97619 .19736 m .97994 .19736 L s .97619 .23142 m .97994 .23142 L s .97619 .24845 m .97994 .24845 L s .97619 .26548 m .97994 .26548 L s .97619 .28251 m .97994 .28251 L s .97619 .31657 m .97994 .31657 L s .97619 .33361 m .97994 .33361 L s .97619 .35064 m .97994 .35064 L s .97619 .36767 m .97994 .36767 L s .97619 .40173 m .97994 .40173 L s .97619 .41876 m .97994 .41876 L s .97619 .43579 m .97994 .43579 L s .97619 .45282 m .97994 .45282 L s .97619 .48688 m .97994 .48688 L s .97619 .50392 m .97994 .50392 L s .97619 .52095 m .97994 .52095 L s .97619 .53798 m .97994 .53798 L s .97619 .02705 m .97994 .02705 L s .97619 .01002 m .97994 .01002 L s .97619 .57204 m .97994 .57204 L s .97619 .58907 m .97994 .58907 L s .97619 .6061 m .97994 .6061 L s .25 Mabswid .97619 0 m .97619 .61803 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath .5 Mabswid .02381 .47597 m .05858 .43919 L .0965 .39832 L .13212 .35933 L .16637 .32145 L .20283 .28093 L .23792 .24204 L .27522 .20125 L .31116 .16303 L .34573 .12792 L .38251 .0933 L .41793 .06378 L .45198 .04032 L .47091 .02998 L .48824 .02257 L .49679 .01973 L .50131 .01847 L .50613 .01731 L .51075 .01639 L .51494 .01572 L .51697 .01546 L .51916 .01522 L .52124 .01503 L .52313 .0149 L .52415 .01484 L .52527 .01479 L .52644 .01475 L .52756 .01473 L .5286 .01472 L .5297 .01472 L .53064 .01473 L .53169 .01475 L .53228 .01477 L .5329 .01479 L .53402 .01484 L .53522 .01491 L .53653 .015 L .53888 .01522 L .54108 .01548 L .54588 .01623 L .55024 .01716 L .56003 .02008 L .56902 .02386 L .5773 .02833 L .59602 .04224 L .6129 .05974 L .63148 .08502 L .66779 .15494 L .70273 .25008 L .73988 .37842 L Mistroke .75845 .44813 L .77566 .5105 L .79233 .56224 L .79713 .57433 L .80165 .58415 L .80601 .59198 L .81008 .59761 L .81211 .59975 L .81429 .60151 L .81552 .60226 L .8167 .60279 L .81775 .60312 L .81889 .60331 L .82001 .60332 L .82064 .60325 L .82124 .60313 L .82186 .60295 L .82243 .60274 L .82372 .60208 L .82476 .60137 L .8259 .60041 L .82819 .59784 L .83038 .59459 L .83236 .59093 L .83681 .58022 L .84094 .56703 L .84482 .55163 L .8536 .5062 L .86287 .44319 L .87949 .3107 L .88095 .29969 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]HGAYHf4PAg9QL6QYHgOol2000I Ool00`00Oomoo`1/Ool20006Ool008aoo`8001aoo`8006aoo`03001oogoo00Eoo`00RWoo0P0087oo 00<007ooOol0JGoo00<007ooOol01Goo0028Ool2000SOol00`00Oomoo`1DOol20004Ool20003Ool4 0005Ool00`00Oomoo`05Ool008Ioo`8002Ioo`03001oogoo059oo`04001oogoo0003Ool20004Ool0 0`00Oomoo`05Ool00`00Oomoo`05Ool008Eoo`03001oogoo02Moo`8004aoo`D00005Ool007ooOol0 000:Ool00`00Oomoo`04Ool30005Ool008=oo`8002aoo`03001oogoo04moo`04001oogoo000;Ool0 0`00Oomoo`03Ool00`00Oomoo`05Ool0089oo`03001oogoo02eoo`03001oogoo04ioo`04001oogoo 0008Ool01000Oomoo`001Goo00<007ooOol01Goo0020Ool2000aOol00`00Oomoo`1>Ool2000:Ool2 0006Ool00`00Oomoo`05Ool007moo`03001oogoo039oo`03001oogoo065oo`03001oogoo00Eoo`00 OGoo0P00=Goo00<007ooOol0HGoo0P001Woo001lOol00`00Oomoo`0fOol00`00Oomoo`1POol00`00 Oomoo`05Ool007]oo`03001oogoo03Qoo`03001oogoo05moo`03001oogoo00Eoo`00NWoo00<007oo Ool0>Woo00<007ooOol0GWoo00<007ooOol01Goo001iOol00`00Oomoo`0kOol00`00Oomoo`1NOol2 0006Ool007Qoo`03001oogoo03eoo`03001oogoo05eoo`03001oogoo00Eoo`00Mgoo00<007ooOol0 ?goo00<007ooOol0G7oo00<007ooOol01Goo001fOol00`00Oomoo`10Ool00`00Oomoo`1LOol00`00 Oomoo`05Ool007Eoo`03001oogoo049oo`03001oogoo05]oo`03001oogoo00Eoo`00Lgoo0P00AGoo 00<007ooOol0Fgoo0P001Woo001bOol00`00Oomoo`16Ool00`00Oomoo`1JOol00`00Oomoo`05Ool0 075oo`03001oogoo04Moo`03001oogoo05Yoo`03001oogoo00Eoo`00L7oo00<007ooOol0BGoo00<0 07ooOol0FGoo00<007ooOol01Goo001_Ool00`00Oomoo`1:Ool00`00Oomoo`1IOol00`00Oomoo`05 Ool006ioo`03001oogoo04aoo`03001oogoo05Qoo`8000Ioo`00KGoo00<007ooOol0CGoo00<007oo Ool0F7oo00<007ooOol01Goo001/Ool00`00Oomoo`1?Ool00`00Oomoo`0mOol20004Ool20003Ool3 0003Ool30006Ool00`00Oomoo`05Ool006]oo`03001oogoo051oo`03001oogoo03aoo`04001oogoo 0003Ool20004Ool00`00Oomoo`05Ool00`00Oomoo`03Ool00`00Oomoo`05Ool006Yoo`03001oogoo 059oo`03001oogoo03Eoo`D00005Ool007ooOol00009Ool00`00Oomoo`05Ool00`00Oomoo`03Ool3 0005Ool006Yoo`03001oogoo059oo`03001oogoo03]oo`04001oogoo0009Ool00`00Oomoo`02Ool3 0006Ool00`00Oomoo`05Ool006Uoo`03001oogoo05Aoo`03001oogoo03Yoo`04001oogoo0009Ool0 0`00Oomoo`02Ool00`00Oomoo`06Ool00`00Oomoo`05Ool006Qoo`03001oogoo05Eoo`03001oogoo 03]oo`8000Uoo`8000Aoo`@000Eoo`03001oogoo00Eoo`00Igoo00<007ooOol0Egoo00<007ooOol0 E7oo00<007ooOol01Goo001VOol00`00Oomoo`1HOol00`00Oomoo`1DOol20006Ool006Eoo`03001o ogoo05Yoo`03001oogoo05=oo`03001oogoo00Eoo`00I7oo00<007ooOol0Fgoo00<007ooOol0Dgoo 00<007ooOol01Goo001SOol00`00Oomoo`1MOol00`00Oomoo`1BOol00`00Oomoo`05Ool0069oo`03 001oogoo05ioo`03001oogoo059oo`03001oogoo00Eoo`00HGoo00<007ooOol0H7oo00<007ooOol0 DGoo0P001Woo001POol00`00Oomoo`1QOol00`00Oomoo`1AOol00`00Oomoo`05Ool005moo`03001o ogoo069oo`03001oogoo055oo`03001oogoo00Eoo`00GWoo00<007ooOol0I7oo00<007ooOol0D7oo 00<007ooOol01Goo001MOol00`00Oomoo`1UOol00`00Oomoo`1@Ool00`00Oomoo`05Ool005aoo`03 001oogoo06Moo`03001oogoo04moo`8000Ioo`00Fgoo00<007ooOol0J7oo00<007ooOol0Cgoo00<0 07ooOol01Goo001JOol00`00Oomoo`1YOol00`00Oomoo`1?Ool00`00Oomoo`05Ool005Uoo`03001o ogoo06]oo`03001oogoo04ioo`03001oogoo00Eoo`00F7oo00<007ooOol0K7oo00<007ooOol0CWoo 0P001Woo001GOol00`00Oomoo`1]Ool00`00Oomoo`1>Ool00`00Oomoo`05Ool005Moo`03001oogoo 06ioo`03001oogoo04eoo`03001oogoo00Eoo`00EWoo00<007ooOol0Kgoo00<007ooOol0>Goo0P00 17oo0P000goo0`001Woo00<007ooOol01Goo001EOol00`00Oomoo`1aOol00`00Oomoo`0gOol01000 Oomoo`000goo0P0017oo00<007ooOol01Goo00<007ooOol01Goo001DOol00`00Oomoo`1bOol00`00 Oomoo`0aOol500001Goo001oogoo00002Goo00<007ooOol01Goo0`001Goo001COol00`00Oomoo`1c Ool00`00Oomoo`0gOol01000Oomoo`002Goo00<007ooOol01Goo00<007ooOol01Goo001BOol00`00 Oomoo`1eOol00`00Oomoo`0fOol01000Oomoo`002Goo00<007ooOol01Goo00<007ooOol01Goo001B Ool00`00Oomoo`1eOol00`00Oomoo`0gOol20009Ool20007Ool00`00Oomoo`05Ool0055oo`03001o ogoo07Moo`03001oogoo04Yoo`03001oogoo00Eoo`00D7oo00<007ooOol0N7oo00<007ooOol0BWoo 0P001Woo001?Ool00`00Oomoo`1iOol00`00Oomoo`1:Ool00`00Oomoo`05Ool004ioo`03001oogoo 07]oo`03001oogoo04Uoo`03001oogoo00Eoo`00CGoo00<007ooOol0O7oo00<007ooOol0BGoo00<0 07ooOol01Goo001goo00<007ooOol0U7oo00<007ooOol0:7oo00<007oo Ool067oo00<007ooOol01Goo000jOol00`00Oomoo`2EOol00`00Oomoo`0XOol00`00Oomoo`0HOol2 0006Ool003Uoo`03001oogoo09Ioo`03001oogoo02Qoo`03001oogoo01Qoo`03001oogoo00Eoo`00 >Goo00<007ooOol0Ugoo00<007ooOol09goo00<007ooOol067oo00<007ooOol01Goo000hOol00`00 Oomoo`2HOol00`00Oomoo`0WOol00`00Oomoo`0HOol00`00Oomoo`05Ool003Moo`03001oogoo09Uo o`03001oogoo02Ioo`03001oogoo01Uoo`03001oogoo00Eoo`00=Woo00<007ooOol0VWoo00<007oo Ool09Woo00<007ooOol06Goo0P001Woo000eOol00`00Oomoo`2LOol00`00Oomoo`0UOol00`00Oomo o`0IOol00`00Oomoo`05Ool003Aoo`03001oogoo09eoo`03001oogoo02Eoo`03001oogoo01Uoo`03 001oogoo00Eoo`0027oo0P0017oo0P0017oo0P0017oo0P000goo0`004Goo00<007ooOol00goo0P00 17oo0P0017oo0P0017oo0P0017oo0P0017oo0P0057oo0P0017oo0P0017oo0P0017oo0P0017oo0P00 17oo0P005Goo0P0017oo0P0017oo0P0017oo0P0017oo0P0017oo0`004Goo00@007oo000000Aoo`80 00Aoo`8000Aoo`8000Aoo`8000=oo`@000Eoo`03001oogoo01Uoo`03001oogoo00Eoo`001goo00@0 07ooOol000=oo`8000=oo`04001oogoo0002Ool01000Oomoo`000goo00<007ooOol03goo00<007oo Ool00goo00@007ooOol000=oo`8000=oo`04001oogoo0002Ool01000Oomoo`000Woo00@007ooOol0 009oo`04001oogoo000BOol01000Oomoo`000goo0P000goo00@007ooOol0009oo`04001oogoo0002 Ool01000Oomoo`000Woo00@007ooOol001=oo`04001oogoo0003Ool20003Ool01000Oomoo`000Woo 00@007ooOol0009oo`04001oogoo0004Ool00`00Oomoo`0AOol01000Oomoo`000goo0P000goo00@0 07ooOol0009oo`04001oogoo0002Ool01000Oomoo`000goo00<007ooOol01Goo00<007ooOol06Goo 0P001Woo00000goo0000000300001Goo001oogoo000027oo00@007ooOol0009oo`04001oogoo0003 Ool00`00Oomoo`0>Ool600001Goo001oogoo000027oo00@007ooOol0009oo`04001oogoo0002Ool0 1000Oomoo`000Woo00@007ooOol000aoo`D00005Ool007ooOol00008Ool01000Oomoo`000Woo00@0 07ooOol0009oo`04001oogoo0002Ool01000Oomoo`003Goo1@0000Eoo`00Oomoo`0000Qoo`04001o ogoo0002Ool01000Oomoo`000Woo00@007ooOol0009oo`@000aoo`D00005Ool007ooOol00008Ool0 1000Oomoo`000Woo00@007ooOol0009oo`04001oogoo0004Ool00`00Oomoo`04Ool00`00Oomoo`0I Ool00`00Oomoo`05Ool000Moo`04001oogoo0008Ool01000Oomoo`000Woo00@007ooOol000=oo`03 001oogoo00eoo`03001oogoo00Eoo`04001oogoo0008Ool01000Oomoo`000Woo00@007ooOol0009o o`04001oogoo0003Ool2000COol01000Oomoo`0027oo00@007ooOol0009oo`04001oogoo0002Ool0 1000Oomoo`000Woo0`0057oo00@007ooOol000Qoo`04001oogoo0002Ool01000Oomoo`000Woo00@0 07ooOol0009oo`03001oo`0001=oo`04001oogoo0008Ool01000Oomoo`000Woo00@007ooOol0009o o`04001oogoo0005Ool00`00Oomoo`03Ool00`00Oomoo`0IOol00`00Oomoo`05Ool000Moo`04001o ogoo0008Ool01000Oomoo`000Woo00@007ooOol000=oo`03001oogoo00aoo`03001oogoo00Ioo`04 001oogoo0008Ool01000Oomoo`000Woo00@007ooOol0009oo`04001oogoo0002Ool01000Oomoo`00 4Woo00@007ooOol000Qoo`04001oogoo0002Ool01000Oomoo`000Woo00@007ooOol0009oo`03001o ogoo01Aoo`04001oogoo0008Ool01000Oomoo`000Woo00@007ooOol0009oo`04001oogoo0002Ool0 0`00Ool0000COol01000Oomoo`0027oo00@007ooOol0009oo`04001oogoo0002Ool01000Oomoo`00 0Woo00@007ooOol000Eoo`03001oogoo01Uoo`03001oogoo00Eoo`0027oo0P002Woo0P0017oo0P00 0goo0P003Woo00<007ooOol01goo0P002Woo0P0017oo0P0017oo0P0017oo0P0057oo0P002Woo0P00 17oo0P0017oo0P0017oo0`0057oo0P002Woo0P0017oo0P0017oo0P0017oo0P0057oo0P002Woo0P00 17oo0P0017oo0P0017oo0P001Goo00<007ooOol06Woo00<007ooOol01Goo000^Ool00`00Oomoo`2U Ool00`00Oomoo`0ROol00`00Oomoo`0JOol20006Ool002eoo`03001oogoo0:Ioo`03001oogoo029o o`03001oogoo01Yoo`03001oogoo00Eoo`00;7oo00<007ooOol0Z7oo00<007ooOol08Goo00<007oo Ool06Woo00<007ooOol01Goo000[Ool00`00Oomoo`2YOol00`00Oomoo`0QOol00`00Oomoo`0JOol0 0`00Oomoo`05Ool002Yoo`03001oogoo0:Yoo`03001oogoo025oo`03001oogoo01Yoo`03001oogoo 00Eoo`002gooo`005P00000BOol00`00Oomoo`0:Ool00`00Oomoo`07Ool01000Oomoo`0037oo00<0 07ooOol02Woo00<007ooOol02Woo00<007ooOol02Woo00<007ooOol02Woo00<007ooOol02Woo00<0 07ooOol02goo00<007ooOol02Woo00<007ooOol02Woo00<007ooOol02Woo00<007ooOol02Woo00<0 07ooOol02Woo00<007ooOol02Woo00<007oo00002Woo00<007ooOol02goo00<007ooOol01goo00@0 07ooOol000aoo`03001oogoo00Yoo`03001oogoo00Eoo`004Woo00<007ooOol04goo00<007ooOol0 6goo00<007ooOol01oo`03001oogoo01=oo`03001oogoo01moo`8000Ioo`00hGoo00<007ooOol0 4Woo00<007ooOol07goo00<007ooOol01Goo003QOol00`00Oomoo`0BOol00`00Oomoo`0OOol00`00 Oomoo`05Ool00>5oo`03001oogoo019oo`03001oogoo01moo`03001oogoo00Eoo`00hGoo00<007oo Ool04Goo00<007ooOol087oo0P001Woo003ROol00`00Oomoo`0@Ool00`00Oomoo`0POol00`00Oomo o`05Ool00>9oo`03001oogoo011oo`03001oogoo021oo`03001oogoo00Eoo`00hWoo00<007ooOol0 47oo00<007ooOol087oo00<007ooOol01Goo003SOol00`00Oomoo`0?Ool00`00Oomoo`0POol00`00 Oomoo`05Ool00>=oo`03001oogoo00moo`03001oogoo021oo`8000Ioo`00hgoo00<007ooOol03goo 00<007ooOol087oo00<007ooOol01Goo003SOol00`00Oomoo`0>Ool00`00Oomoo`0QOol00`00Oomo o`05Ool00>Aoo`03001oogoo00eoo`03001oogoo00eoo`8000Aoo`8000=oo`<000Ioo`03001oogoo 00Eoo`00i7oo00<007ooOol03Goo00<007ooOol037oo00@007ooOol000=oo`8000Aoo`03001oogoo 00Eoo`03001oogoo00Eoo`00i7oo00<007ooOol03Goo00<007ooOol037oo00@007ooOol000Uoo`03 001oogoo00Eoo`<000Eoo`00iGoo00<007ooOol037oo00<007ooOol037oo00@007ooOol000Uoo`03 001oogoo00Eoo`03001oogoo00Eoo`00iGoo00<007ooOol02goo00<007ooOol03Goo00@007ooOol0 00Uoo`03001oogoo00Eoo`03001oogoo00Eoo`00iGoo00<007ooOol02goo00<007ooOol03Woo0P00 2Goo0P001goo00<007ooOol01Goo003VOol00`00Oomoo`0:Ool00`00Oomoo`0ROol00`00Oomoo`05 Ool00>Moo`03001oogoo00Uoo`03001oogoo029oo`8000Ioo`00igoo00<007ooOol027oo00<007oo Ool08goo00<007ooOol01Goo003XOol00`00Oomoo`07Ool00`00Oomoo`0SOol00`00Oomoo`05Ool0 0>Qoo`03001oogoo00Moo`03001oogoo02=oo`03001oogoo00Eoo`00j7oo00<007ooOol01Woo00<0 07ooOol097oo0P001Woo003YOol00`00Oomoo`04Ool00`00Oomoo`0UOol00`00Oomoo`05Ool00>Uo o`03001oogoo00Aoo`03001oogoo02Eoo`03001oogoo00Eoo`00jWoo0P000goo00<007ooOol09Woo 00<007ooOol01Goo003[Ool4000YOol00`00Oomoo`05Ool00?mooaUoo`8000Ioo`00ogoo6Goo00<0 07ooOol01Goo003oOolIOol00`00Oomoo`05Ool00?mooaUoo`03001oogoo00Eoo`00ogoo8Goo003o OolQOol00?moob5oo`00\ \>"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {-0.00106859, -0.232951, 3.81045*^-6, 0.00213083}}] }, Open ]], Cell[TextData["Let's find a root for\[Omega] exactly:"], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(aroot = FindRoot[vAtL[\[Omega], k, LR, Lf, L], {\[Omega], \(- .001\), \(- .0008\)}]\)], "Input"], Cell[BoxData[ \({\[Omega] \[Rule] \(-0.000910550669933957301`\)}\)], "Output"] }, Open ]], Cell["Now find v[y] that satisfies the boundary conditions:", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(hitsoln = soln[\[Omega] /. aroot[\([1]\)], k, LR, Lf, L]\)], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{"v", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.`, 1.`}}, "<>"]\), False, Editable->False]}], "}"}], "}"}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(ParametricPlot[{v[y] /. hitsoln[\([1]\)], y}, {y, 0, 1}, AspectRatio -> 2, Ticks -> {{}, Automatic}]; \)\)], "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 2.99209 0.047619 1.90476 [ [.01131 .42857 -18 -4.5 ] [.01131 .42857 0 4.5 ] [.01131 .80952 -18 -4.5 ] [.01131 .80952 0 4.5 ] [.01131 1.19048 -18 -4.5 ] [.01131 1.19048 0 4.5 ] [.01131 1.57143 -18 -4.5 ] [.01131 1.57143 0 4.5 ] [.01131 1.95238 -6 -4.5 ] [.01131 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 0 .04762 m 1 .04762 L s .02381 .42857 m .03006 .42857 L s [(0.2)] .01131 .42857 1 0 Mshowa .02381 .80952 m .03006 .80952 L s [(0.4)] .01131 .80952 1 0 Mshowa .02381 1.19048 m .03006 1.19048 L s [(0.6)] .01131 1.19048 1 0 Mshowa .02381 1.57143 m .03006 1.57143 L s [(0.8)] .01131 1.57143 1 0 Mshowa .02381 1.95238 m .03006 1.95238 L s [(1)] .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 .52381 m .02756 .52381 L s .02381 .61905 m .02756 .61905 L s .02381 .71429 m .02756 .71429 L s .02381 .90476 m .02756 .90476 L s .02381 1 m .02756 1 L s .02381 1.09524 m .02756 1.09524 L s .02381 1.28571 m .02756 1.28571 L s .02381 1.38095 m .02756 1.38095 L s .02381 1.47619 m .02756 1.47619 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 0 0 m 1 0 L 1 2 L 0 2 L closepath clip newpath .5 Mabswid .02381 .04762 m .14486 .12489 L .27457 .20916 L .39203 .2883 L .49911 .36442 L .60483 .44543 L .69683 .52342 L .78239 .60631 L .85145 .68617 L .90436 .763 L .92581 .80189 L .94514 .84473 L .95922 .88537 L .96417 .90357 L .96861 .92343 L .97059 .9342 L .97216 .94419 L .97352 .95454 L .97451 .96393 L .97526 .97316 L .97553 .97739 L .97577 .98186 L .97596 .98664 L .97603 .98877 L .97609 .99104 L .97614 .99349 L .97617 .99578 L .97618 .99717 L .97619 .9985 L .97619 1.00098 L .97618 1.00325 L .97615 1.00571 L .97611 1.00803 L .97606 1.01015 L .9759 1.01509 L .97566 1.02028 L .97536 1.02536 L .97502 1.03008 L .97403 1.0408 L .97174 1.05861 L .96839 1.07764 L .96374 1.09812 L .95749 1.12033 L .94296 1.16063 L .90103 1.24247 L .84557 1.32127 L .77916 1.39705 L .69552 1.47773 L .60379 1.55538 L .49589 1.63793 L Mistroke .38366 1.71745 L .26985 1.79394 L .1445 1.87534 L .02381 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]HGAYHf4PAg9QL6QYHgOol00`00Oomoo`1iOol001=oo`03001oogoo00moo`8007Uoo`004goo00<007ooOol04Goo00<0 07ooOol0MWoo000COol00`00Oomoo`0BOol2001fOol001=oo`8001Eoo`03001oogoo07=oo`004goo 00<007ooOol05Goo0P00Lgoo000COol00`00Oomoo`0GOol00`00Oomoo`1`Ool001=oo`03001oogoo 01Qoo`80071oo`004goo00<007ooOol06Woo00<007ooOol0KGoo000COol00`00Oomoo`0KOol2001] Ool001=oo`03001oogoo01eoo`03001oogoo06Yoo`004goo00<007ooOol07Woo0P00JWoo000COol0 0`00Oomoo`0POol00`00Oomoo`1WOol001=oo`03001oogoo025oo`8006Moo`004goo00<007ooOol0 8goo00<007ooOol0I7oo000COol00`00Oomoo`0TOol2001TOol001=oo`03001oogoo02Ioo`03001o ogoo065oo`004goo0P00:7oo0P00HGoo000COol00`00Oomoo`0YOol2001OOol001=oo`03001oogoo 02]oo`03001oogoo05aoo`004goo00<007ooOol0;7oo0P00G7oo000COol00`00Oomoo`0^Ool00`00 Oomoo`1IOol001=oo`03001oogoo02moo`8005Uoo`004goo00<007ooOol0Goo0P00Cgoo000COol00`00Oomoo`0kOol00`00Oomo o`1Woo00001Goo001oogoo00002goo00@007ooOol0051oo`03 001oogoo03Uoo`0000Eoo`00Oomoo`0000Qoo`04001oogoo0002Ool00`00Oomoo`1?Ool00`00Oomo o`0hOol0009oo`8000Yoo`8000=oo`03001oogoo051oo`8003Qoo`004goo00<007ooOol0DWoo00<0 07ooOol0=Goo000COol00`00Oomoo`1COol00`00Oomoo`0dOol001=oo`03001oogoo05Aoo`03001o ogoo03=oo`004goo00<007ooOol0EGoo00<007ooOol0Ool00`00Oomoo`0iOol00005Ool0 07ooOol00008Ool01000Oomoo`000Woo0P00CWoo00<007ooOol0>Woo00001Goo001oogoo00002Goo 0P000goo00<007ooOol0C7oo00<007ooOol0>goo00001Goo001oogoo000027oo00@007ooOol0009o o`03001oogoo04Yoo`8003ioo`000Woo0P002Woo0P000goo00<007ooOol0BGoo00<007ooOol0?Woo 000COol00`00Oomoo`18Ool00`00Oomoo`0oOol001=oo`03001oogoo04Moo`03001oogoo041oo`00 4goo00<007ooOol0AGoo0P00@goo000COol00`00Oomoo`14Ool00`00Oomoo`13Ool001=oo`03001o ogoo04=oo`03001oogoo04Aoo`004goo00<007ooOol0@Woo00<007ooOol0AGoo000COol00`00Oomo o`10Ool20018Ool001=oo`03001oogoo03moo`03001oogoo04Qoo`004goo00<007ooOol0?Goo0P00 Bgoo000COol00`00Oomoo`0lOol00`00Oomoo`1;Ool001=oo`8003]oo`8004ioo`004goo00<007oo Ool0>Goo00<007ooOol0CWoo000COol00`00Oomoo`0hOol00`00Oomoo`1?Ool001=oo`03001oogoo 03Ioo`80059oo`004goo00<007ooOol0=Goo00<007ooOol0DWoo000COol00`00Oomoo`0cOol2001E Ool001=oo`03001oogoo039oo`03001oogoo05Eoo`004goo00<007ooOol0<7oo0P00F7oo000COol0 0`00Oomoo`0_Ool00`00Oomoo`1HOol001=oo`03001oogoo02eoo`8005]oo`004goo00<007ooOol0 ;7oo00<007ooOol0Fgoo000COol00`00Oomoo`0ZOol2001NOol001=oo`03001oogoo02Qoo`80061o o`004goo00<007ooOol09goo00<007ooOol0H7oo000COol2000VOol2001SOol001=oo`03001oogoo 02Aoo`03001oogoo06=oo`004goo00<007ooOol08Woo0P00IWoo000COol00`00Oomoo`0QOol00`00 Oomoo`1VOol001=oo`03001oogoo01moo`8006Uoo`004goo00<007ooOol07Woo00<007ooOol0JGoo 000COol00`00Oomoo`0LOol2001/Ool001=oo`03001oogoo01]oo`03001oogoo06aoo`004goo00<0 07ooOol06Woo00<007ooOol0KGoo000COol00`00Oomoo`0HOol2001`Ool001=oo`03001oogoo01Mo o`03001oogoo071oo`004goo00<007ooOol05Goo0P00Lgoo000COol00`00Oomoo`0DOol00`00Oomo o`1cOol001=oo`8001Aoo`03001oogoo07Aoo`004goo00<007ooOol04Goo0P00Mgoo000COol00`00 Oomoo`0@Ool00`00Oomoo`1gOol001=oo`03001oogoo00ioo`8007Yoo`004goo00<007ooOol03Goo 00<007ooOol0NWoo000COol00`00Oomoo`0;Ool2001mOol001=oo`03001oogoo00Uoo`8007moo`00 4goo00<007ooOol027oo00<007ooOol0Ogoo000COol00`00Oomoo`06Ool20022Ool001=oo`03001o ogoo00Aoo`8008Aoo`004goo00<007ooOol00goo00<007ooOol0Q7oo000COol01000Oomoogoo0P00 Qgoo000=Ool30003Ool01000Ool00000RGoo000>Ool00`00Oomoo`02Ool2002;Ool000ioo`03001o ogoo009oo`8008]oo`003Woo00<007ooOol00Woo00<007ooOol0RWoo000>Ool00`00Oomoo`02Ool0 0`00Oomoo`2:Ool000eoo`8000Aoo`03001oogoo08Yoo`004goo00<007ooOol0RWoo000COol00`00 Oomoo`2:Ool001=oo`03001oogoo08Yoo`004goo00<007ooOol0RWoo002POol00:1oo`00X7oo002P Ool00:1oo`00X7oo002POol00:1oo`00X7oo002POol00:1oo`00X7oo002POol00:1oo`00X7oo002P Ool00001\ \>"], ImageRangeCache->{{{0, 159.062}, {319.125, 0}} -> {-0.046344, -0.0871419, 0.00234249, 0.0036797}}] }, Open ]], Cell["\<\ Clearly the fun is just beginning, but the rest is up to you!\ \>", "Text"] }, Closed]] }, Open ]] }, FrontEndVersion->"Macintosh 3.0", ScreenRectangle->{{0, 1152}, {0, 850}}, AutoGeneratedPackage->None, WindowSize->{519, 675}, WindowMargins->{{79, Automatic}, {Automatic, 8}}, 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, 125, 3, 107, "Title"], Cell[CellGroupData[{ Cell[1881, 58, 76, 1, 63, "Section", InitializationCell->True], Cell[1960, 61, 349, 8, 137, "Input", InitializationCell->True], Cell[2312, 71, 154, 3, 65, "Input", InitializationCell->True] }, Closed]], Cell[CellGroupData[{ Cell[2503, 79, 83, 1, 37, "Section", InitializationCell->True], Cell[2589, 82, 99, 2, 47, "Input", InitializationCell->True] }, Closed]], Cell[CellGroupData[{ Cell[2725, 89, 26, 0, 37, "Section"], Cell[CellGroupData[{ Cell[2776, 93, 73, 1, 47, "Input"], Cell[2852, 96, 224, 6, 52, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[3113, 107, 146, 2, 83, "Input"], Cell[3262, 111, 13183, 312, 304, 2676, 178, "GraphicsData", "PostScript", "Graphics"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[16494, 429, 47, 0, 37, "Section"], Cell[16544, 431, 63, 0, 53, "Text"], Cell[CellGroupData[{ Cell[16632, 435, 110, 2, 47, "Input"], Cell[16745, 439, 13625, 441, 194, 4819, 328, "GraphicsData", "PostScript", "Graphics"] }, Open ]], Cell[30385, 883, 97, 3, 70, "Text"], Cell[CellGroupData[{ Cell[30507, 890, 122, 3, 47, "Input"], Cell[30632, 895, 18398, 504, 194, 4763, 331, "GraphicsData", "PostScript", "Graphics"] }, Open ]], Cell[49045, 1402, 64, 0, 53, "Text"], Cell[CellGroupData[{ Cell[49134, 1406, 134, 3, 65, "Input"], Cell[49271, 1411, 82, 1, 52, "Output"] }, Open ]], Cell[49368, 1415, 69, 0, 53, "Text"], Cell[CellGroupData[{ Cell[49462, 1419, 89, 1, 47, "Input"], Cell[49554, 1422, 224, 6, 52, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[49815, 1433, 146, 2, 83, "Input"], Cell[49964, 1437, 14293, 327, 337, 2690, 179, "GraphicsData", "PostScript", "Graphics"] }, Open ]], Cell[64272, 1767, 88, 3, 53, "Text"] }, Closed]] }, Open ]] } ] *) (*********************************************************************** End of Mathematica Notebook file. ***********************************************************************)