(*********************************************************************** 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[ 33821, 1137]*) (*NotebookOutlinePosition[ 34525, 1163]*) (* CellTagsIndexPosition[ 34481, 1159]*) (*WindowFrame->Normal*) Notebook[{ Cell[CellGroupData[{ Cell["Linear Water Waves", "Title", TextAlignment->Center], Cell[TextData[{ "We solve for very small ampitude waves in a layer of \"water\" between \ z=-H and z=0. All fields are assumed to be proportional to ", Cell[BoxData[ \(TraditionalForm\`e\^\(i(kx - \[Omega]t)\)\)]], ". The velocity potential is \[Phi] and the surface displacement is \ \[Eta]. We know that a separable solution for \[Phi] is \[Phi]=f(z)", Cell[BoxData[ \(TraditionalForm\`e\^\(i(kx - \[Omega]t)\)\)]], "where f(z)=A", Cell[BoxData[ \(TraditionalForm\`e\^kz\)]], "+B", Cell[BoxData[ \(TraditionalForm\`e\^\(-kz\)\)]], ". We assume \[Eta]=C", Cell[BoxData[ \(TraditionalForm\`e\^\(i(kx - \[Omega]t)\)\)]], ". A, B and C are complex constants. We find \[Omega], A and B. See \ Kundu, page 190." }], "Text"], Cell[BoxData[ \(\[Phi] := f*phase\)], "Input"], Cell[BoxData[ \(\[Eta] := C*phase\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(phase = Exp[I \((k*x - \[Omega]*t)\)]\)], "Input"], Cell[BoxData[ \(E\^\(I\ \((k\ x - t\ \[Omega])\)\)\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(f = A*Exp[k*z] + B*Exp[\(-k\)*z]\)], "Input"], Cell[BoxData[ \(B\ E\^\(\(-k\)\ z\) + A\ E\^\(k\ z\)\)], "Output"] }, Open ]], Cell[TextData[{ "The bottom boundary condition is w=0, or ", Cell[BoxData[ \(TraditionalForm\`\[PartialD]\[Phi]\/\[PartialD]z\)]], "=0, at z=-H." }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(bottombc = \((D[\[Phi], z] == 0)\) /. z -> \(-H\)\)], "Input"], Cell[BoxData[ \(E\^\(I\ \((k\ x - t\ \[Omega])\)\)\ \((A\ E\^\(\(-H\)\ k\)\ k - B\ E\^\(H\ k\)\ k)\) == 0\)], "Output"] }, Open ]], Cell[TextData[{ "The top dynamic boundary condition is ", Cell[BoxData[ \(TraditionalForm\`\[PartialD]\[Phi]\/\[PartialD]t\)]], "=-g\[Eta], at z=0." }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(dynbc = \((D[\[Phi], t] == \(-g\)*\[Eta])\) /. z -> 0\)], "Input"], Cell[BoxData[ \(\(-I\)\ \((A + B)\)\ E\^\(I\ \((k\ x - t\ \[Omega])\)\)\ \[Omega] == \(-C\)\ E\^\(I\ \((k\ x - t\ \[Omega])\)\)\ g\)], "Output"] }, Open ]], Cell[TextData[{ "The top kinematic bottom boundary ", Cell[BoxData[ \(TraditionalForm\`\[PartialD]\[Eta]\/\[PartialD]t\)]], "=", Cell[BoxData[ \(TraditionalForm\`\[PartialD]\[Phi]\/\[PartialD]z\)]], ", at z=0." }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(kinbc = \((D[\[Phi], z] == D[\[Eta], t])\) /. z -> 0\)], "Input"], Cell[BoxData[ \(E\^\(I\ \((k\ x - t\ \[Omega])\)\)\ \((A\ k - B\ k)\) == \(-I\)\ C\ E\^\(I\ \((k\ x - t\ \[Omega])\)\)\ \[Omega]\)], "Output"] }, Open ]], Cell[TextData[ "The three boundary condition equationshave a trivial solution A=B=C=0. We \ seek the nontrivial solution that occurs when the eiqenvalue \[Omega] assumes \ certain values. We first eliminate two of the unknowns A and B and reduce \ the three equations to one. "], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(oneeq = Simplify[Eliminate[{bottombc, dynbc, kinbc}, {A, B}]]\)], "Input"], Cell[BoxData[ \(\(-C\)\ E\^\(\(-H\)\ k + 2\ I\ \((k\ x - t\ \[Omega])\)\)\ \((\(-\((\(-1\) + E\^\(2\ H\ k\))\)\)\ g\ k + \((1 + E\^\(2\ H\ k\))\)\ \[Omega]\^2)\) == 0\)], "Output"] }, Open ]], Cell[TextData[ "We have no use for the C=0 solution in the above equation, so we go on to \ solve for \[Omega]:"], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(dis1 = Solve[oneeq, \[Omega]]\)], "Input"], Cell[BoxData[ \({{\[Omega] \[Rule] \(-\(\(\@\(\(-1\) + E\^\(2\ H\ k\)\)\ \@g\ \@k\)\/\@\(1 + E\^\(2\ H\ k\)\)\)\)}, { \[Omega] \[Rule] \(\@\(\(-1\) + E\^\(2\ H\ k\)\)\ \@g\ \@k\)\/\@\(1 + E\^\(2\ H\ k\)\)}}\)], "Output"] }, Open ]], Cell["\<\ I will next try converting the above rules to a more familiar \ dispersion relation:\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(dis2 = \[Omega]^2 == \((\[Omega]^2 /. dis1[\([1, 1]\)])\)\)], "Input"], Cell[BoxData[ \(\[Omega]\^2 == \(\((\(-1\) + E\^\(2\ H\ k\))\)\ g\ k\)\/\(1 + E\^\(2\ H\ k\)\)\)], "Output"] }, Open ]], Cell["\<\ That still does not look too familiar. Let's beat on it some more:\ \ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(dis3 = FullSimplify[ExpToTrig[dis2]]\)], "Input"], Cell[BoxData[ \(\[Omega]\^2 == g\ k\ Tanh[H\ k]\)], "Output"] }, Open ]], Cell["\<\ What do the waves look like? Let's eliminate C between two \ equations and find A and B in terms of C:\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(AB = Flatten[Solve[{bottombc, kinbc}, {A, B}]]\)], "Input"], Cell[BoxData[ \({A \[Rule] \(-\(\(I\ C\ E\^\(2\ H\ k\)\ \[Omega]\)\/\(\((\(-1\) + E\^\(2\ H\ k\))\)\ k\)\)\), B \[Rule] \(-\(\(I\ C\ \[Omega]\)\/\(\((\(-1\) + E\^\(2\ H\ k\))\)\ k\)\)\)} \)], "Output"] }, Open ]], Cell["Now find f(z):", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(f = FullSimplify[ExpToTrig[Factor[f /. AB]]]\)], "Input"], Cell[BoxData[ \(\(-\(\(I\ C\ \[Omega]\ Cosh[k\ \((H + z)\)]\ Csch[H\ k]\)\/k\)\)\)], "Output"] }, Open ]], Cell[TextData[ "Recall \[Phi] was defined with a delayed assignment, \[Phi]:=f*phase, so we \ can use \[Phi] and the above f will be used:"], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(w = D[\[Phi], z]\)], "Input"], Cell[BoxData[ \(\(-I\)\ C\ E\^\(I\ \((k\ x - t\ \[Omega])\)\)\ \[Omega]\ Csch[H\ k]\ Sinh[k\ \((H + z)\)]\)], "Output"] }, Open ]], Cell[BoxData[ \(TextForm \`We\ have\ use\ for \[PartialD]\[Eta]\/\[PartialD]t\ \(also : \)\)], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(detadt = D[\[Eta], t]\)], "Input"], Cell[BoxData[ \(\(-I\)\ C\ E\^\(I\ \((k\ x - t\ \[Omega])\)\)\ \[Omega]\)], "Output"] }, Open ]], Cell[BoxData[ FormBox[ RowBox[{"We", " ", "normalize", " ", "w", " ", "by", RowBox[{ FormBox[\(\[PartialD]\[Eta]\/\[PartialD]t\), "TextForm"], ":"}]}], TextForm]], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(wnorm = w/detadt\)], "Input"], Cell[BoxData[ \(Csch[H\ k]\ Sinh[k\ \((H + z)\)]\)], "Output"] }, Open ]], Cell["\<\ We also normalize the height, expressing height by a dimensionless \ Z:\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(wnondim = wnorm /. z -> Z*H\)], "Input"], Cell[BoxData[ \(Csch[H\ k]\ Sinh[k\ \((H + H\ Z)\)]\)], "Output"] }, Open ]], Cell["\<\ We next make a plot of w(z) for a moderate wavelength k=7/H and a \ long wavelength k=1/H.\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(\(ParametricPlot[{{wnondim /. k -> 7/H, Z}, {wnondim /. k -> 1/H, Z}}, \n\t{Z, \(-1\), 0}, PlotStyle -> {{RGBColor[1. , 0. , 0. ]}, {RGBColor[0. , 0. , 1. ]}}]; \)\)], "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.0238095 0.952381 0.603319 0.588604 [ [.21429 .59082 -9 -9 ] [.21429 .59082 9 0 ] [.40476 .59082 -9 -9 ] [.40476 .59082 9 0 ] [.59524 .59082 -9 -9 ] [.59524 .59082 9 0 ] [.78571 .59082 -9 -9 ] [.78571 .59082 9 0 ] [.97619 .59082 -3 -9 ] [.97619 .59082 3 0 ] [.01131 .01472 -12 -4.5 ] [.01131 .01472 0 4.5 ] [.01131 .13244 -24 -4.5 ] [.01131 .13244 0 4.5 ] [.01131 .25016 -24 -4.5 ] [.01131 .25016 0 4.5 ] [.01131 .36788 -24 -4.5 ] [.01131 .36788 0 4.5 ] [.01131 .4856 -24 -4.5 ] [.01131 .4856 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 .60332 m .21429 .60957 L s [(0.2)] .21429 .59082 0 1 Mshowa .40476 .60332 m .40476 .60957 L s [(0.4)] .40476 .59082 0 1 Mshowa .59524 .60332 m .59524 .60957 L s [(0.6)] .59524 .59082 0 1 Mshowa .78571 .60332 m .78571 .60957 L s [(0.8)] .78571 .59082 0 1 Mshowa .97619 .60332 m .97619 .60957 L s [(1)] .97619 .59082 0 1 Mshowa .125 Mabswid .07143 .60332 m .07143 .60707 L s .11905 .60332 m .11905 .60707 L s .16667 .60332 m .16667 .60707 L s .2619 .60332 m .2619 .60707 L s .30952 .60332 m .30952 .60707 L s .35714 .60332 m .35714 .60707 L s .45238 .60332 m .45238 .60707 L s .5 .60332 m .5 .60707 L s .54762 .60332 m .54762 .60707 L s .64286 .60332 m .64286 .60707 L s .69048 .60332 m .69048 .60707 L s .7381 .60332 m .7381 .60707 L s .83333 .60332 m .83333 .60707 L s .88095 .60332 m .88095 .60707 L s .92857 .60332 m .92857 .60707 L s .25 Mabswid 0 .60332 m 1 .60332 L s .02381 .01472 m .03006 .01472 L s [(-1)] .01131 .01472 1 0 Mshowa .02381 .13244 m .03006 .13244 L s [(-0.8)] .01131 .13244 1 0 Mshowa .02381 .25016 m .03006 .25016 L s [(-0.6)] .01131 .25016 1 0 Mshowa .02381 .36788 m .03006 .36788 L s [(-0.4)] .01131 .36788 1 0 Mshowa .02381 .4856 m .03006 .4856 L s [(-0.2)] .01131 .4856 1 0 Mshowa .125 Mabswid .02381 .04415 m .02756 .04415 L s .02381 .07358 m .02756 .07358 L s .02381 .10301 m .02756 .10301 L s .02381 .16187 m .02756 .16187 L s .02381 .1913 m .02756 .1913 L s .02381 .22073 m .02756 .22073 L s .02381 .27959 m .02756 .27959 L s .02381 .30902 m .02756 .30902 L s .02381 .33845 m .02756 .33845 L s .02381 .39731 m .02756 .39731 L s .02381 .42674 m .02756 .42674 L s .02381 .45617 m .02756 .45617 L s .02381 .51503 m .02756 .51503 L s .02381 .54446 m .02756 .54446 L s .02381 .57389 m .02756 .57389 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 1 0 0 r .5 Mabswid .02381 .01472 m .02431 .03859 L .0249 .06463 L .02555 .08909 L .02632 .11261 L .02736 .13765 L .02865 .16175 L .03047 .18736 L .0328 .21204 L .03578 .23578 L .04002 .26104 L .04548 .28536 L .05245 .30874 L .06233 .33364 L .07505 .3576 L .09319 .38308 L .11671 .40762 L .14682 .43123 L .18965 .45635 L .24492 .48053 L .32397 .50623 L .42675 .53099 L .55874 .55482 L .74689 .58016 L .97619 .60332 L s 0 0 1 r .02381 .01472 m .05669 .03859 L .09262 .06463 L .12648 .08909 L .15922 .11261 L .1943 .13765 L .22835 .16175 L .26493 .18736 L .3006 .21204 L .33538 .23578 L .37293 .26104 L .4097 .28536 L .44567 .30874 L .48471 .33364 L .52306 .3576 L .56474 .38308 L .60585 .40762 L .64634 .43123 L .69053 .45635 L .73422 .48053 L .78197 .50623 L .82934 .53099 L .87626 .55482 L .92771 .58016 L .97619 .60332 L s % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{35, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgGoOI`7eI>ooooiVIVH:2P]7AdTi>C_ooohVHZWB1T<[AfDI=E_oooe1FIDA9EJf`]b@W;_oooiBD U:f][B4Q8Kfm_MkNgV=SHoXC@3_j=55Anooo003oool0od92@^7QhLS8b:n_[iJFUWemOFATI4];Bc8b <[no_gmoOdaL1000K0@03?`4101T10TWW0@006`400cl10@0K0@03B@410>@1000K0@03?`4101`1 0TWT0@006`400cl00@0N0@03B@410>41000K0@03?`4101l10TWQ0@006`400cl10@0Q0@03B@410=h1 000K0@03?`41028100=90@40g@4001/100042BLP1 000K0@03004o03X100=90@40a@4001/100<00Cl0>`42BLD1000K0@03004o03d100=90@40`P4001/1 0P000cl10@0l0@03B@410<41000K0@03004o03l10TW10@006`400`01?`110@03B@410;h1000K0@03 004o048100=90@40_@4001/100<00Cl0@`42BKd1000K0@03004o04D100=90@40^P4001/100<00Cl0 AP400dT10@2i0@006`400`01?`170@99^@4001/10P000cl10@170@99]`4001/100@00@4oBP400dT1 0@2d0@006`4010010Cm;0@99]04001/100@00@4oC@42BK81000K0@040041?dl100=90@40[`4001/1 00@00@4oD042BJl100080@8010400`010@030@801@4010010CmB0@03B@410:`100070@04004100P1 00@00@40104010010CmC0@03B@410:/100070@04004100P100@00@4010420081000@006`400`010@060@03?`4106T10TV>0@006`400`010@060@03?`4106/1 00=90@40R`4001/100<00@401`400cl10@1[0@99R`4001/100<00@401`400cl10@1]0@03B@4108P1 000K0@802@400cl10@1]0@99R04001/100<00@4020400cl10@1_0@99QP4001/100<00@402@400cl1 0@1`0@99Q04001/100<00@402@400cl10@1b0@99PP4001/100<00@402P400cl10@1c0@99P04001/1 00<00@402P400cl10@1e0@03B@4107d1000K0@03004100/10004001/10P1E0@DoL@42BCH1000K0@03004105T1 13m_0@99=04001/100<00@40G@45?f`10TTb0@006`400`010@1R0@@oJP42BC01000K0@03004106H1 1CmW0@99;P4001/100<00@40J`45?f@10TT/0@006`400`010@1`0@DoH@42BBX1000K0@03004107D1 1CmN0@99:04001/10P1k0@DoF`42BBH1000K0@03004107l11CmH0@99904001/100<00@40Q045?eD1 0TTR0@006`400`010@290@LoD042BB01000K0@030041090123m:0@997P4001/100<00@40V048?d@1 0TTL0@006`400`010@2P0@To?@42BAX1000K0@03004102P10P040@0300410081100T0@8010400`01 0@040@03004102<10P040@03004100<10P0=0@Po3`4200@100<00@400`4201X10TT>0@@01P4001/1 0P0X0@04004100P100<00@40904010010@0:0@030041028100@00@40204010010@0D0@Po1P401001 0@080@04004101/10TT>0@03004100D1000K0@03004102L100@00@402@400`010@0S0@04004100L1 1@0S0@04004100P100@00@407049?`03004100H100@00@407@42B@`100<00@401@4001/100<00@40 9`4010010@0:0@030041028100@00@401`4010010@0T0@04004100P10`0S0@03004100To00<10000 8042B@X100<00@401@4001/100<00@409`4010010@0;0@030041024100@00@4020400`01000T0@04 004100P100<00@408`4010010@080@Po7@42B@P100<00@401@4001/100<00@409`4010010@080@04 004102<100@00@402@4202@100@00@4020400`010@0S0@04004100P100@00@401049?aH10TT40@<0 1`4001/100<00@40:04200X10P0U0@802`400`010@0S0@802P4302<10P0:0@803P49?`l10TT40@03 004100D1000K0@0300410>D123l90@992P4001/100<00@40k@49?`810TT80@00507o004013l00dT0 000500006`400`010@090@03004100X100<00@402P400`010@090@03004100X100<00@402P400`01 0@090@03004100X100<00@402P400`010@090@03004100X100<00@402P400`010@090@03004100X1 00<00@402P400`010@090@03004100X100<00@402P400`010@090@03004100X100<00@401@4001/1 00<00@40o`430@006`400`010@3o0@<1000K0@0300410?l10`400?l18@400?l18@400?l18@400?l1 8@400?l18@400001\ \>"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {-0.107864, -1.0669, 0.00394728, 0.00638683}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(u = D[\[Phi], x]\)], "Input"], Cell[BoxData[ \(C\ E\^\(I\ \((k\ x - t\ \[Omega])\)\)\ \[Omega]\ Cosh[k\ \((H + z)\)]\ Csch[H\ k]\)], "Output"] }, Open ]], Cell["\<\ We next make a plot of u(z) for a moderate wavelength k=7/H and a \ long wavelength k=1/H .But first we make a dimensionless profile:\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(unondim = u/\((\[Eta]*\[Omega])\) /. z -> Z*H\)], "Input"], Cell[BoxData[ \(General::"spell1" \( : \ \) "Possible spelling error: new symbol name \"\!\(unondim\)\" is similar \ to existing symbol \"\!\(wnondim\)\"."\)], "Message"], Cell[BoxData[ \(Cosh[k\ \((H + H\ Z)\)]\ Csch[H\ k]\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(ParametricPlot[{{unondim /. k -> 7/H, Z}, {unondim /. k -> 1/H, Z}}, \n\t{Z, \(-1\), 0}, PlotStyle -> {{RGBColor[1. , 0. , 0. ]}, {RGBColor[0. , 0. , 1. ]}}]; \)\)], "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.0238095 0.725328 0.603319 0.588604 [ [.16888 .59082 -9 -9 ] [.16888 .59082 9 0 ] [.31394 .59082 -9 -9 ] [.31394 .59082 9 0 ] [.45901 .59082 -9 -9 ] [.45901 .59082 9 0 ] [.60407 .59082 -9 -9 ] [.60407 .59082 9 0 ] [.74914 .59082 -3 -9 ] [.74914 .59082 3 0 ] [.8942 .59082 -9 -9 ] [.8942 .59082 9 0 ] [.01131 .01472 -12 -4.5 ] [.01131 .01472 0 4.5 ] [.01131 .13244 -24 -4.5 ] [.01131 .13244 0 4.5 ] [.01131 .25016 -24 -4.5 ] [.01131 .25016 0 4.5 ] [.01131 .36788 -24 -4.5 ] [.01131 .36788 0 4.5 ] [.01131 .4856 -24 -4.5 ] [.01131 .4856 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 .16888 .60332 m .16888 .60957 L s [(0.2)] .16888 .59082 0 1 Mshowa .31394 .60332 m .31394 .60957 L s [(0.4)] .31394 .59082 0 1 Mshowa .45901 .60332 m .45901 .60957 L s [(0.6)] .45901 .59082 0 1 Mshowa .60407 .60332 m .60407 .60957 L s [(0.8)] .60407 .59082 0 1 Mshowa .74914 .60332 m .74914 .60957 L s [(1)] .74914 .59082 0 1 Mshowa .8942 .60332 m .8942 .60957 L s [(1.2)] .8942 .59082 0 1 Mshowa .125 Mabswid .06008 .60332 m .06008 .60707 L s .09634 .60332 m .09634 .60707 L s .13261 .60332 m .13261 .60707 L s .20514 .60332 m .20514 .60707 L s .24141 .60332 m .24141 .60707 L s .27767 .60332 m .27767 .60707 L s .35021 .60332 m .35021 .60707 L s .38647 .60332 m .38647 .60707 L s .42274 .60332 m .42274 .60707 L s .49527 .60332 m .49527 .60707 L s .53154 .60332 m .53154 .60707 L s .56781 .60332 m .56781 .60707 L s .64034 .60332 m .64034 .60707 L s .6766 .60332 m .6766 .60707 L s .71287 .60332 m .71287 .60707 L s .7854 .60332 m .7854 .60707 L s .82167 .60332 m .82167 .60707 L s .85794 .60332 m .85794 .60707 L s .93047 .60332 m .93047 .60707 L s .96674 .60332 m .96674 .60707 L s .25 Mabswid 0 .60332 m 1 .60332 L s .02381 .01472 m .03006 .01472 L s [(-1)] .01131 .01472 1 0 Mshowa .02381 .13244 m .03006 .13244 L s [(-0.8)] .01131 .13244 1 0 Mshowa .02381 .25016 m .03006 .25016 L s [(-0.6)] .01131 .25016 1 0 Mshowa .02381 .36788 m .03006 .36788 L s [(-0.4)] .01131 .36788 1 0 Mshowa .02381 .4856 m .03006 .4856 L s [(-0.2)] .01131 .4856 1 0 Mshowa .125 Mabswid .02381 .04415 m .02756 .04415 L s .02381 .07358 m .02756 .07358 L s .02381 .10301 m .02756 .10301 L s .02381 .16187 m .02756 .16187 L s .02381 .1913 m .02756 .1913 L s .02381 .22073 m .02756 .22073 L s .02381 .27959 m .02756 .27959 L s .02381 .30902 m .02756 .30902 L s .02381 .33845 m .02756 .33845 L s .02381 .39731 m .02756 .39731 L s .02381 .42674 m .02756 .42674 L s .02381 .45617 m .02756 .45617 L s .02381 .51503 m .02756 .51503 L s .02381 .54446 m .02756 .54446 L s .02381 .57389 m .02756 .57389 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 1 0 0 r .5 Mabswid .02513 .01472 m .02513 .01544 L .02513 .0161 L .02513 .01687 L .02513 .01759 L .02513 .01887 L .02514 .02026 L .02514 .0218 L .02514 .02342 L .02515 .02634 L .02515 .02935 L .02516 .03217 L .02519 .03859 L .02522 .04514 L .02526 .05114 L .02537 .06463 L .0255 .07637 L .02569 .08918 L .02592 .10235 L .02617 .1143 L .02684 .13848 L .02779 .16324 L .02903 .18707 L .03082 .21241 L .03314 .23681 L .03611 .26028 L .04035 .28526 L .04581 .30931 L .05361 .33487 L .06374 .3595 L .07673 .38319 L .09523 .40839 L .11912 .43266 L .14959 .45599 L .19283 .48084 L .24842 .50475 L .32773 .53018 L .43048 .55466 L .56194 .57822 L .74884 .60328 L .74914 .60332 L s 0 0 1 r .641 .01472 m .641 .01544 L .64101 .0161 L .64101 .01687 L .64101 .01759 L .64102 .01887 L .64103 .02026 L .64105 .0218 L .64107 .02342 L .64112 .02634 L .64119 .02935 L .64128 .03217 L .64151 .03859 L .64183 .04514 L .64219 .05114 L .64323 .06463 L .64439 .07637 L .64595 .08918 L .64986 .1143 L .6547 .13848 L .66101 .16418 L .66824 .18894 L .67627 .21277 L .68599 .23811 L .69651 .26251 L .70772 .28598 L .72084 .31096 L .73466 .33501 L .75065 .36057 L .76735 .3852 L .78465 .40889 L .8044 .43409 L .82478 .45836 L .84565 .48169 L .86929 .50654 L .89347 .53045 L .92076 .55587 L .94862 .58036 L .97619 .60332 L s % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{35, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgGoOI`7eI>ooooiVIVH:2P]7AdTi>C_ooohVHZWB1T<[AfDI=E_oooe1FIDA9EJf`]b@W;_oooiBD U:f][B4Q8Kfm_MkNgV=SHoXC@3_j=55Anooo003oool0od92@^7QhLS8b:n_[iJFUWemOFATI4];Bc8b <[no_gmoOda0@006`4010010Cn`0@03B@4104h1000K0@0500410Cl0/0400dT10@1=0@00 6`42008100`4001/100<00@403@400cl1 0@2d0@03B@4103/1000K0@803`400cl10@2d0@03B@4103X1000K0@03004100l100@4001/100<00@4040400cl10@2d0@03B@4103P1 000K0@0300410141000@006`400`010@0K0@04004100T1 00<00@405`4010010@070@D05`4010010@080@04004101411Sl20@03004100H100@00@407`400`01 0@0M0@03004100P100<00@401P400dT10@0=0@006`400`010@0K0@04004100X100<00@405P401001 0@070@04004101P100@00@40204301L100<0?cl01Cl50@8080400`010@0M0@03004100T100<00@40 1P400dT10@0<0@006`400`010@0K0@04004100/100<00@405@4010010@080@03004001P100@00@40 20400`010@0G0@04004100@11cl00`010@0M0@03004101d100<00@402P400`010@060@03B@4100/1 000K0@03004101/100@00@40204010010@0G0@04004100T10P0H0@04004100P100<00@405`401001 0@080@03004100Lo5`4301d10`090@04004100T100=90@402P4001/100<00@40704200X10P0I0@80 2`400`010@0G0@802P4301L10P0:0@801`48?a4100<00@407@400`010@080@802`42B@X1000K0@03 00410:d11cm40@03B@4100L1000K0@0300410;@11cln0@03B@4100H1000D0LD013lk0003B@0000D0 000K0@03004100H100<00@401`400`010@070@03004100H100<00@401`400`010@070@03004100H1 00<00@401`400`010@070@03004100H100<00@401`400`010@070@03004100H100<00@401`400`01 0@070@03004100H100<00@401`400`010@060@03004100L100<00@401`400`010@060@03004100L1 00<00@401`400`010@060@03004100L100<00@401`400`010@070@006`400`010@3o0@<1000K0@03 00410?l10`4001/100<00@40o`430@00o`4Q0@00o`4Q0@00o`4Q0@00o`4Q0@00o`4Q0@00\ \>"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {-0.14163, -1.0669, 0.00518292, 0.00638683}}] }, Open ]], Cell[TextData[{ "What have we learned that is new? Hopefully a lot about ", StyleBox["Mathematica", FontSlant->"Italic"], " ! But also in the plots we are seeing how the assumptions of shallow \ water theory are becoming apparent in the long-wavelength wave. Feel free to \ play with the above plot. Try a very long wavelength, just for fun." }], "Text"] }, Open ]] }, FrontEndVersion->"X 3.0", ScreenRectangle->{{0, 1280}, {0, 1024}}, WindowToolbars->{}, WindowSize->{783, 659}, WindowMargins->{{Automatic, -76}, {Automatic, 60}}, Magnification->1, StyleDefinitions -> "ok.nb" ] (*********************************************************************** 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, 60, 1, 47, "Title"], Cell[1794, 54, 779, 20, 104, "Text"], Cell[2576, 76, 50, 1, 48, "Input"], Cell[2629, 79, 50, 1, 48, "Input"], Cell[CellGroupData[{ Cell[2704, 84, 70, 1, 48, "Input"], Cell[2777, 87, 68, 1, 56, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[2882, 93, 65, 1, 48, "Input"], Cell[2950, 96, 70, 1, 56, "Output"] }, Open ]], Cell[3035, 100, 169, 5, 56, "Text"], Cell[CellGroupData[{ Cell[3229, 109, 82, 1, 48, "Input"], Cell[3314, 112, 132, 2, 56, "Output"] }, Open ]], Cell[3461, 117, 172, 5, 56, "Text"], Cell[CellGroupData[{ Cell[3658, 126, 86, 1, 48, "Input"], Cell[3747, 129, 155, 2, 56, "Output"] }, Open ]], Cell[3917, 134, 245, 8, 56, "Text"], Cell[CellGroupData[{ Cell[4187, 146, 85, 1, 48, "Input"], Cell[4275, 149, 153, 2, 56, "Output"] }, Open ]], Cell[4443, 154, 289, 4, 87, "Text"], Cell[CellGroupData[{ Cell[4757, 162, 97, 2, 48, "Input"], Cell[4857, 166, 206, 3, 56, "Output"] }, Open ]], Cell[5078, 172, 123, 2, 53, "Text"], Cell[CellGroupData[{ Cell[5226, 178, 62, 1, 48, "Input"], Cell[5291, 181, 288, 6, 79, "Output"] }, Open ]], Cell[5594, 190, 108, 3, 53, "Text"], Cell[CellGroupData[{ Cell[5727, 197, 90, 1, 48, "Input"], Cell[5820, 200, 122, 3, 73, "Output"] }, Open ]], Cell[5957, 206, 93, 3, 53, "Text"], Cell[CellGroupData[{ Cell[6075, 213, 69, 1, 48, "Input"], Cell[6147, 216, 65, 1, 56, "Output"] }, Open ]], Cell[6227, 220, 127, 3, 53, "Text"], Cell[CellGroupData[{ Cell[6379, 227, 79, 1, 48, "Input"], Cell[6461, 230, 253, 6, 76, "Output"] }, Open ]], Cell[6729, 239, 30, 0, 53, "Text"], Cell[CellGroupData[{ Cell[6784, 243, 77, 1, 48, "Input"], Cell[6864, 246, 101, 2, 70, "Output"] }, Open ]], Cell[6980, 251, 149, 2, 53, "Text"], Cell[CellGroupData[{ Cell[7154, 257, 49, 1, 48, "Input"], Cell[7206, 260, 130, 2, 56, "Output"] }, Open ]], Cell[7351, 265, 113, 3, 70, "Text"], Cell[CellGroupData[{ Cell[7489, 272, 54, 1, 48, "Input"], Cell[7546, 275, 89, 1, 56, "Output"] }, Open ]], Cell[7650, 279, 212, 5, 70, "Text"], Cell[CellGroupData[{ Cell[7887, 288, 49, 1, 48, "Input"], Cell[7939, 291, 66, 1, 53, "Output"] }, Open ]], Cell[8020, 295, 95, 3, 53, "Text"], Cell[CellGroupData[{ Cell[8140, 302, 60, 1, 48, "Input"], Cell[8203, 305, 69, 1, 53, "Output"] }, Open ]], Cell[8287, 309, 114, 3, 53, "Text"], Cell[CellGroupData[{ Cell[8426, 316, 215, 4, 67, "Input"], Cell[8644, 322, 11161, 351, 194, 3691, 254, "GraphicsData", "PostScript", "Graphics"] }, Open ]], Cell[CellGroupData[{ Cell[19842, 678, 49, 1, 48, "Input"], Cell[19894, 681, 122, 2, 56, "Output"] }, Open ]], Cell[20031, 686, 157, 4, 70, "Text"], Cell[CellGroupData[{ Cell[20213, 694, 78, 1, 48, "Input"], Cell[20294, 697, 181, 3, 47, "Message"], Cell[20478, 702, 69, 1, 53, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[20584, 708, 215, 4, 67, "Input"], Cell[20802, 714, 12620, 410, 194, 4447, 305, "GraphicsData", "PostScript", "Graphics"] }, Open ]], Cell[33437, 1127, 368, 7, 87, "Text"] }, Open ]] } ] *) (*********************************************************************** End of Mathematica Notebook file. ***********************************************************************)