(************** Content-type: application/mathematica ************** CreatedBy='Mathematica 5.2' Mathematica-Compatible Notebook This notebook can be used with any Mathematica-compatible application, such as Mathematica, MathReader or Publicon. The data for the notebook starts with the line containing 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[ 18785, 469]*) (*NotebookOutlinePosition[ 19464, 492]*) (* CellTagsIndexPosition[ 19420, 488]*) (*WindowFrame->Normal*) Notebook[{ Cell[CellGroupData[{ Cell[TextData[StyleBox["Sullivan Vortex Stability", "Subtitle"]], "Title"], Cell["\<\ Method 1: Fast, Using Known Inner Products on P.78 of my thesis \ appendix\ \>", "Subsubtitle"], Cell[BoxData[ \(EqA28\ := \ 2 j\ KroneckerDelta[j, k]\)], "Input", CellLabel->"In[23]:="], Cell[BoxData[ \(EqA29\ \ := \ If[j \[LessEqual] k, \ 1, \ 0]\)], "Input", CellLabel->"In[24]:="], Cell[BoxData[ \(EqA30\ := \ \(-2\^\(-\((j + k + 1)\)\)\)\ \(\((j + k)\)!\)\/\(\(j!\)\ \ \(k!\)\)\)], "Input", CellLabel->"In[25]:="], Cell[BoxData[ \(a[l_, m_] := \ \(a[l, m]\ = \ If[l \[Equal] 1\ || \ m \[Equal] 1, 2\^m - 1, a[l - 1, m] + a[l, m - 1]]\)\)], "Input", CellLabel->"In[26]:="], Cell[BoxData[ \(EqA31\ := \ \(-2\^\(-\((j + k + 1)\)\)\)\ a[j + 1, k + 1]\)], "Input", CellLabel->"In[27]:="], Cell[BoxData[ \(A\ := \ \(-EqA28\)\ + 6\ \((EqA29\ + \ EqA30\ + EqA31)\)\ // Simplify\)], "Input", CellLabel->"In[28]:="], Cell["\<\ Ndim is set to 10 here to show the matrix output in a reasonable \ size, but you can easily go to N>100\ \>", "Note"], Cell[BoxData[ \(Ndim := 10\)], "Input", CellLabel->"In[29]:="], Cell[BoxData[ \(\(\(Ajk\)\(\ \)\(:=\)\(\ \)\(Table[ A, \ {j, 0, Ndim}, \ {k, 0, Ndim}]\)\(\ \)\)\)], "Input", CellLabel->"In[30]:="], Cell[CellGroupData[{ Cell[BoxData[ \(Ajk\ // \ MatrixForm\)], "Input", CellLabel->"In[31]:="], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0"}, {\(-3\), \(-\(1\/2\)\), \(3\/4\), \(3\/8\), \(3\/16\), \(3\/32\), \ \(3\/64\), \(3\/128\), \(3\/256\), \(3\/512\), \(3\/1024\)}, {\(-\(3\/2\)\), \(-3\), \(-\(17\/8\)\), \(9\/8\), \(21\/32\), \(3\ \/8\), \(27\/128\), \(15\/128\), \(33\/512\), \(9\/256\), \(39\/2048\)}, {\(-\(3\/4\)\), \(-\(15\/8\)\), \(-3\), \(-\(63\/16\)\), \(87\/64\ \), \(111\/128\), \(69\/128\), \(21\/64\), \(201\/1024\), \(237\/2048\), \(69\ \/1024\)}, {\(-\(3\/8\)\), \(-\(9\/8\)\), \(-\(33\/16\)\), \(-3\), \ \(-\(745\/128\)\), \(195\/128\), \(33\/32\), \(87\/128\), \(897\/2048\), \ \(567\/2048\), \(705\/4096\)}, {\(-\(3\/16\)\), \(-\(21\/32\)\), \(-\(87\/64\)\), \ \(-\(279\/128\)\), \(-3\), \(-\(1981\/256\)\), \(843\/512\), \(1191\/1024\), \ \(3279\/4096\), \(4413\/8192\), \(5823\/16384\)}, {\(-\(3\/32\)\), \(-\(3\/8\)\), \(-\(111\/128\)\), \ \(-\(195\/128\)\), \(-\(579\/256\)\), \(-3\), \(-\(9909\/1024\)\), \ \(1785\/1024\), \(10419\/8192\), \(927\/1024\), \(20655\/32768\)}, {\(-\(3\/64\)\), \(-\(27\/128\)\), \(-\(69\/128\)\), \ \(-\(33\/32\)\), \(-\(843\/512\)\), \(-\(2379\/1024\)\), \(-3\), \ \(-\(23815\/2048\)\), \(29847\/16384\), \(44679\/32768\), \(32667\/32768\)}, {\(-\(3\/128\)\), \(-\(15\/128\)\), \(-\(21\/64\)\), \(-\(87\/128\ \)\), \(-\(1191\/1024\)\), \(-\(1785\/1024\)\), \(-\(4857\/2048\)\), \(-3\), \ \(-\(445289\/32768\)\), \(61839\/32768\), \(47253\/32768\)}, {\(-\(3\/256\)\), \(-\(33\/512\)\), \(-\(201\/1024\)\), \ \(-\(897\/2048\)\), \(-\(3279\/4096\)\), \(-\(10419\/8192\)\), \ \(-\(29847\/16384\)\), \(-\(78999\/32768\)\), \(-3\), \ \(-\(1019505\/65536\)\), \(254649\/131072\)}, {\(-\(3\/512\)\), \(-\(9\/256\)\), \(-\(237\/2048\)\), \ \(-\(567\/2048\)\), \(-\(4413\/8192\)\), \(-\(927\/1024\)\), \ \(-\(44679\/32768\)\), \(-\(61839\/32768\)\), \(-\(160143\/65536\)\), \(-3\), \ \(-\(4595015\/262144\)\)} }, RowSpacings->1, ColumnSpacings->1, ColumnAlignments->{Left}], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]], "Output", CellLabel->"Out[31]//MatrixForm="] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Sort[Eigenvalues[\(-Ajk\)\ // N]]\ // MatrixForm\)], "Input", CellLabel->"In[32]:="], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", TagBox[GridBox[{ {"0.`"}, {"1.643043501284098`"}, {"3.148566939666461`"}, {"4.77927899071293`"}, {"6.538698571012099`"}, {"8.357368430515805`"}, {"10.207854472048982`"}, {"12.089628613028907`"}, {\(\(\(13.658136745802251`\)\(\[InvisibleSpace]\)\) - 4.124581429298101`\ \[ImaginaryI]\)}, {\(\(\(13.658136745802251`\)\(\[InvisibleSpace]\)\) + 4.124581429298101`\ \[ImaginaryI]\)}, {"14.019701266249344`"} }, RowSpacings->1, ColumnAlignments->{Left}], Column], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]], "Output", CellLabel->"Out[32]//MatrixForm="] }, Open ]], Cell["\<\ Here's a one line version what I had to use a messy FORTRAN program \ to do back in 92 or 93... for N=100.\ \>", "Note"], Cell[CellGroupData[{ Cell[BoxData[ \(Sort[ Eigenvalues[ Table[\(-A\), \ {j, 0, 300}, \ {k, 0, 300}]\ // N]]\)], "Input", CellLabel->"In[33]:="], Cell[BoxData[ \({0.`, 1.6428866828103132`, 3.1483919458683918`, 4.7835014548081825`, 6.5572364810547645`, 8.402252988209401`, 10.286107860990413`, 12.19445626281926`, 14.119707610779471`, 16.057260379775773`, 18.00409767031128`, 19.95814167094741`, 21.91790884063662`, 23.88230875172764`, 25.850520122740193`, 27.821911318166784`, 29.79598768287191`, 31.77235566582875`, 33.75069777229529`, 35.730754692006805`, 37.712312297562114`, 39.695192020829175`, 41.67924361803611`, 43.664339653820605`, 45.65037124194308`, 47.63724471836508`, 49.62487901438844`, 51.61320356350321`, 53.602156616254355`, 55.59168387601221`, 57.5817373813115`, 59.572274588477214`, 61.56325761004214`, 63.554652579922625`, 65.54642912238356`, 67.53855990224949`, 69.53102024581402`, 71.52378781659186`, 73.51684233864707`, 75.51016535874035`, 77.50374004121025`, 79.49755099006843`, 81.49158409506278`, 83.48582639663134`, 85.48026596818521`, 87.47489181287683`, 89.46969377240976`, 91.46466244637958`, 93.45978912188299`, 95.45506570952544`, 97.45048468734026`, 99.44603905125089`, 101.44172226962125`, 103.43752824274745`, 105.4334512680934`, 107.42948600622856`, 109.42562745264543`, 111.42187091068901`, 113.41821196823734`, 115.41464647501103`, 117.41117052385178`, 119.4077804325701`, 121.40447272668138`, 123.40124412594132`, 125.39809152951888`, 127.39501200429007`, 129.39200277292818`, 131.38906120340107`, 133.3861848000966`, 135.38337119353807`, 137.38061813360244`, 139.37792348136728`, 141.37528520172972`, 143.37270135791886`, 145.37017010474437`, 147.36768968322323`, 149.36525841608852`, 151.3628747019799`, 153.3605370117797`, 155.35824388452164`, 157.35599392291678`, 159.35378579049834`, 161.35161820786433`, 163.34948994984904`, 165.3473998423027`, 167.3453467595236`, 169.34332962205383`, 171.34134739351583`, 173.33939907936863`, 175.33748372368495`, 177.3356004084247`, 179.33374825046818`, 181.33192640079506`, 183.33013404182762`, 185.3283703872241`, 187.32663467919517`, 189.32492618780722`, 191.32324420942126`, 193.32158806580054`, 195.31995710260605`, 197.31835068860306`, 199.3167682138928`, 201.31520909065767`, 203.31367275016962`, 205.31215864299432`, 207.31066623858868`, 209.30919502354996`, 211.30774450128501`, 213.30631419151854`, 215.30490362940483`, 217.30351236470983`, 219.30213996190952`, 221.30078599831896`, 223.2994500655733`, 225.29813176642236`, 227.2968307169433`, 229.29554654434938`, 231.29427888677967`, 233.29302739373534`, 235.29179172440914`, 237.29057154851267`, 239.28936654495024`, 241.2881764019343`, 243.2870008168769`, 245.28583949514797`, 247.28469215102595`, 249.28355850644107`, 251.28243829072002`, 253.28133124133834`, 255.28023710214825`, 257.2791556244444`, 259.27808656586734`, 261.2770296905941`, 263.27598476913096`, 265.2749515780484`, 267.2739298990907`, 269.27291952061813`, 271.27192023568114`, 273.27093184270933`, 275.2699541457395`, 277.2689869526582`, 279.26803007744763`, 281.2670833377322`, 283.26614655556403`, 285.2652195581006`, 287.264302176303`, 289.2633942446746`, 291.2624956024733`, 293.26160609242953`, 295.2607255605317`, 297.2598538573677`, 299.2589908360486`, 301.25813635372504`, 303.257290270056`, 305.2564524490047`, 307.25562275677635`, 309.25480106237956`, 311.2539872387954`, 313.25318116114715`, 315.2523827064345`, 317.2515917567696`, 319.25080819432765`, 321.25003190482636`, 323.2492627772307`, 325.2485007012429`, 327.24774557046936`, 329.24699727944767`, 331.2462557263519`, 333.2455208099715`, 335.2447924321388`, 337.2440704968176`, 339.2433549089593`, 341.24264557706886`, 343.24194240934196`, 345.24124531827425`, 347.24055421603`, 349.23986901768814`, 351.2391896402562`, 353.23851600079126`, 355.23784802010147`, 357.2371856193812`, 359.2365287212927`, 361.23587725030893`, 363.2352311335972`, 365.23459029670664`, 367.233954670726`, 369.233324184578`, 371.2326987717423`, 373.2320783636926`, 375.2314628969833`, 377.2308523073202`, 379.2302465318684`, 381.2296455102519`, 383.22904918326856`, 385.2284574920373`, 387.22787038081924`, 389.227287792815`, 391.2267096757681`, 393.22613597477766`, 395.22556663981743`, 397.22500161848814`, 399.224440861995`, 401.22388431832314`, 403.2233319408777`, 405.2227836764064`, 407.2222394757833`, 409.2216992851348`, 411.2211630487356`, 413.2206307081674`, 415.2201021967601`, 417.2195774460355`, 419.2190563723314`, 421.21853888644154`, 423.21802488308106`, 425.2175142401982`, 427.2170068144778`, 429.2165024409717`, 431.2160009240689`, 433.2155020338603`, 435.21500550453015`, 437.21451102216065`, 439.21401822474235`, 441.2135266923577`, 443.21303594197315`, 445.2125454216131`, 447.2120545053406`, 449.21156248444055`, 451.21106857105934`, 453.2105718884576`, 455.2100714712831`, 457.20956628062385`, 459.2090551912505`, 461.208537016744`, 463.20801052348065`, 465.2074744525293`, 467.20692754746455`, 469.20636860519255`, 471.2057965173672`, 473.20521034174664`, 475.20460937737`, 477.203993260637`, 479.20336207866995`, 481.20271649668297`, 483.2020579240494`, 485.2013886703346`, 487.20071217095875`, 489.2000331969483`, 491.19935811458777`, 493.1986951773926`, 495.19805482048764`, 497.19745001572056`, 499.1968966371812`, 501.19641385453536`, 503.19602456365214`, 505.1957558477895`, 507.1956394533068`, 509.1957123218499`, 511.1960171341633`, 513.1966029200136`, 515.1975257065673`, 517.198849240363`, 519.2006457951813`, 521.2029970930486`, 523.2059953664399`, 525.2097445898081`, 527.214361930002`, 529.2199794402201`, 531.2267460010237`, 533.2348295399744`, 535.2444194142217`, 537.2557288285867`, 539.2689969292798`, 541.2844899900517`, 543.3025006514257`, 545.323343585317`, 547.3473451703927`, 549.3748238015073`, 551.4060568167615`, 553.441230246822`, 555.4803702704006`, 557.523261760811`, 559.5693704949622`, 561.6177981422181`, 563.6673041691276`, \(\(564.4358923978023`\)\(\[InvisibleSpace]\)\) - 19.4551255414929`\ \[ImaginaryI], \(\(564.4358923978023`\)\(\ \[InvisibleSpace]\)\) + 19.4551255414929`\ \[ImaginaryI], 565.7164150590254`, 567.7636068590984`, 569.8075099492614`, 571.8470724414565`, 573.8816409893841`, 575.9109569351588`, 577.9350945670029`, 579.9543752428024`, 581.9692825642661`, 583.9803911924794`, 585.9883124025613`, 587.9936546621971`, 589.9969961024348`, 591.9988660139066`, 593.9997332421949`}\)], "Output", CellLabel->"Out[33]="] }, Open ]], Cell[CellGroupData[{ Cell["Method 2: Slow, Using Integration", "Subsubtitle"], Cell[BoxData[ \(z := \ r\^2\/2\)], "Input", CellLabel->"In[34]:="], Cell[CellGroupData[{ Cell[BoxData[ \(\[Omega][k_] = \ LaguerreL[k, 0, z]\ \[ExponentialE]\^\(-z\)\)], "Input",\ CellLabel->"In[35]:="], Cell[BoxData[ \(\[ExponentialE]\^\(-\(r\^2\/2\)\)\ LaguerreL[k, r\^2\/2]\)], "Output", CellLabel->"Out[35]="] }, Open ]], Cell[BoxData[ \(InnerProductTest[j_, k_]\ := \ \[Integral]\_0\%\[Infinity] \[Omega][j] \[Omega][k] r\ \(\[ExponentialE]\^z\) \[DifferentialD]r\)], "Input", CellLabel->"In[36]:="], Cell[CellGroupData[{ Cell[BoxData[ \(InnerProductTest[1, 1]\)], "Input", CellLabel->"In[37]:="], Cell[BoxData[ \(1\)], "Output", CellLabel->"Out[37]="] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(InnerProductTest[1, 2]\)], "Input", CellLabel->"In[38]:="], Cell[BoxData[ \(0\)], "Output", CellLabel->"Out[38]="] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(L\[Degree]\[Omega][k_] = 1\/r\ \[PartialD]\_r\ \((r\ \[PartialD]\_r \[Omega][k])\) + \ 1\/r\ \[PartialD]\_r\((r\^2\ \[Omega][k])\)\ // \ Simplify\)], "Input", CellLabel->"In[39]:="], Cell[BoxData[ \(\[ExponentialE]\^\(-\(r\^2\/2\)\)\ \((r\^2\ LaguerreL[\(-2\) + k, 2, r\^2\/2] + \((\(-2\) + r\^2)\)\ LaguerreL[\(-1\) + k, 1, r\^2\/2])\)\)], "Output", CellLabel->"Out[39]="] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(L\[Omega][k_]\ = \ L\[Degree]\[Omega][k]\ - \ 6\/r\ \((1 - \[ExponentialE]\^\(-z\))\)\ \[PartialD]\_r \[Omega][ k]\ - \ 6\ \[ExponentialE]\^\(-z\)\ \[Omega][k]\ // \ Simplify\)], "Input", CellLabel->"In[40]:="], Cell[BoxData[ \(\[ExponentialE]\^\(-r\^2\)\ \((6\ \((\(-2\) + \ \[ExponentialE]\^\(r\^2\/2\))\)\ LaguerreL[k, r\^2\/2] + \[ExponentialE]\^\(r\^2\/2\)\ r\^2\ LaguerreL[\(-2\) \ + k, 2, r\^2\/2] + \((\(-6\) + \[ExponentialE]\^\(r\^2\/2\)\ \((4 + r\^2)\))\)\ LaguerreL[\(-1\) + k, 1, r\^2\/2])\)\)], "Output", CellLabel->"Out[40]="] }, Open ]], Cell[BoxData[ \(\[Omega]L\[Omega][j_, k_]\ := \(\[Omega]L\[Omega][j, k]\ = \ \[Integral]\_0\%\[Infinity] \[Omega][j] L\[Omega][k] r\ \(\[ExponentialE]\^z\) \[DifferentialD]r\)\)], "Input", CellLabel->"In[41]:="], Cell["\<\ In Mathematica the integration above takes O(N^2) time, so we have \ to keep Ndim low. \ \>", "Note"], Cell[BoxData[ \(Ndim\ := \ 6\)], "Input", CellLabel->"In[42]:="], Cell[BoxData[ \(Ajk\ := \ Table[\[Omega]L\[Omega][j, k], \ {j, 0, Ndim}, \ {k, 0, Ndim}]\)], "Input", CellLabel->"In[43]:="], Cell[CellGroupData[{ Cell[BoxData[ \(Ajk\ // MatrixForm\)], "Input", CellLabel->"In[44]:="], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"0", "0", "0", "0", "0", "0", "0"}, {\(-3\), \(-\(1\/2\)\), \(3\/4\), \(3\/8\), \(3\/16\), \(3\/32\), \ \(3\/64\)}, {\(-\(3\/2\)\), \(-3\), \(-\(17\/8\)\), \(9\/8\), \(21\/32\), \(3\ \/8\), \(27\/128\)}, {\(-\(3\/4\)\), \(-\(15\/8\)\), \(-3\), \(-\(63\/16\)\), \(87\/64\ \), \(111\/128\), \(69\/128\)}, {\(-\(3\/8\)\), \(-\(9\/8\)\), \(-\(33\/16\)\), \(-3\), \ \(-\(745\/128\)\), \(195\/128\), \(33\/32\)}, {\(-\(3\/16\)\), \(-\(21\/32\)\), \(-\(87\/64\)\), \ \(-\(279\/128\)\), \(-3\), \(-\(1981\/256\)\), \(843\/512\)}, {\(-\(3\/32\)\), \(-\(3\/8\)\), \(-\(111\/128\)\), \ \(-\(195\/128\)\), \(-\(579\/256\)\), \(-3\), \(-\(9909\/1024\)\)} }, RowSpacings->1, ColumnSpacings->1, ColumnAlignments->{Left}], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]], "Output", CellLabel->"Out[44]//MatrixForm="] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Sort[Eigenvalues[\(-Ajk\)\ // N\ ]]\ // MatrixForm\)], "Input", CellLabel->"In[45]:="], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", TagBox[GridBox[{ {"0.`"}, {"1.6425300919629242`"}, {"3.116636738328387`"}, {"4.627141431831342`"}, {"6.195673550188847`"}, {\(\(\(7.107934875094251`\)\(\[InvisibleSpace]\)\) - 3.1764422221995474`\ \[ImaginaryI]\)}, {\(\(\(7.107934875094251`\)\(\[InvisibleSpace]\)\) + 3.1764422221995474`\ \[ImaginaryI]\)} }, RowSpacings->1, ColumnAlignments->{Left}], Column], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]], "Output", CellLabel->"Out[45]//MatrixForm="] }, Open ]] }, Open ]] }, Open ]] }, FrontEndVersion->"5.2 for Macintosh", ScreenRectangle->{{4, 1920}, {0, 1178}}, WindowSize->{786, 1085}, WindowMargins->{{14, Automatic}, {Automatic, 21}}, StyleDefinitions -> "ArticleModern.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[1776, 53, 74, 0, 62, "Title"], Cell[1853, 55, 105, 3, 28, "Subsubtitle"], Cell[1961, 60, 97, 2, 42, "Input"], Cell[2061, 64, 103, 2, 42, "Input"], Cell[2167, 68, 139, 3, 59, "Input"], Cell[2309, 73, 184, 4, 42, "Input"], Cell[2496, 79, 116, 2, 44, "Input"], Cell[2615, 83, 141, 3, 42, "Input"], Cell[2759, 88, 127, 3, 22, "Note"], Cell[2889, 93, 68, 2, 42, "Input"], Cell[2960, 97, 146, 3, 42, "Input"], Cell[CellGroupData[{ Cell[3131, 104, 79, 2, 42, "Input"], Cell[3213, 108, 2315, 40, 302, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[5565, 153, 108, 2, 42, "Input"], Cell[5676, 157, 898, 23, 212, "Output"] }, Open ]], Cell[6589, 183, 130, 3, 22, "Note"], Cell[CellGroupData[{ Cell[6744, 190, 143, 4, 42, "Input"], Cell[6890, 196, 7009, 103, 506, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[13936, 304, 56, 0, 28, "Subsubtitle"], Cell[13995, 306, 72, 2, 60, "Input"], Cell[CellGroupData[{ Cell[14092, 312, 120, 3, 42, "Input"], Cell[14215, 317, 115, 2, 61, "Output"] }, Open ]], Cell[14345, 322, 205, 4, 57, "Input"], Cell[CellGroupData[{ Cell[14575, 330, 80, 2, 42, "Input"], Cell[14658, 334, 60, 2, 42, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[14755, 341, 80, 2, 42, "Input"], Cell[14838, 345, 60, 2, 42, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[14935, 352, 231, 5, 57, "Input"], Cell[15169, 359, 226, 4, 61, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[15432, 368, 279, 6, 57, "Input"], Cell[15714, 376, 388, 7, 110, "Output"] }, Open ]], Cell[16117, 386, 252, 5, 57, "Input"], Cell[16372, 393, 112, 3, 22, "Note"], Cell[16487, 398, 71, 2, 42, "Input"], Cell[16561, 402, 149, 4, 42, "Input"], Cell[CellGroupData[{ Cell[16735, 410, 77, 2, 42, "Input"], Cell[16815, 414, 1035, 22, 202, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[17887, 441, 110, 2, 42, "Input"], Cell[18000, 445, 745, 19, 148, "Output"] }, Open ]] }, Open ]] }, Open ]] } ] *) (******************************************************************* End of Mathematica Notebook file. *******************************************************************)