(*********************************************************************** 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[ 22911, 816]*) (*NotebookOutlinePosition[ 23966, 852]*) (* CellTagsIndexPosition[ 23922, 848]*) (*WindowFrame->Normal*) Notebook[{ Cell[CellGroupData[{Cell[TextData["CHAPTER 2"], "Section", Evaluatable->False, AspectRatioFixed->True], Cell[CellGroupData[{Cell[TextData["1"], "Subsection", Evaluatable->False, AspectRatioFixed->True], Cell[CellGroupData[{Cell[TextData["1a"], "Subsubsection", Evaluatable->False, AspectRatioFixed->True], Cell[TextData[{ StyleBox["D[f,x]", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[" differentiates ", Evaluatable->False, AspectRatioFixed->True], StyleBox["f", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[" with respect to ", Evaluatable->False, AspectRatioFixed->True], StyleBox["x", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[":", Evaluatable->False, AspectRatioFixed->True] }], "Text", Evaluatable->False, AspectRatioFixed->True], Cell[TextData["n[t_]:=a Exp[r t]"], "Input", AspectRatioFixed->True], Cell[TextData["D[n[t],t]"], "Input", AspectRatioFixed->True], Cell[TextData[{ StyleBox["Verify that this derivative is equal to ", Evaluatable->False, AspectRatioFixed->True], StyleBox["rn[t]", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[" and that ", Evaluatable->False, AspectRatioFixed->True], StyleBox["n[t]", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[" has the correct initial value:", Evaluatable->False, AspectRatioFixed->True] }], "Text", Evaluatable->False, AspectRatioFixed->True], Cell[TextData["r n[t]"], "Input", AspectRatioFixed->True], Cell[TextData["n[0]"], "Input", AspectRatioFixed->True]}, Open]], Cell[CellGroupData[{Cell[TextData["1b"], "Subsubsection", Evaluatable->False, AspectRatioFixed->True], Cell[TextData[{ StyleBox["DSolve[]", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[ " solves a differential equation; you must first clear the existing \ definition of ", Evaluatable->False, AspectRatioFixed->True], StyleBox["n", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[":", Evaluatable->False, AspectRatioFixed->True] }], "Text", Evaluatable->False, AspectRatioFixed->True], Cell[TextData["Clear[n]"], "Input", AspectRatioFixed->True], Cell[TextData["DSolve[{n'[t]==r n[t],n[0]==a},n[t],t]"], "Input", AspectRatioFixed->True]}, Open]]}, Open]], Cell[CellGroupData[{Cell[TextData["2"], "Subsection", Evaluatable->False, AspectRatioFixed->True], Cell[TextData[ "I first calculate and plot the logarithms of the numbers."], "Text", Evaluatable->False, AspectRatioFixed->True], Cell[TextData["logdata=N[Log[{8,30,81,282,705,1325}]];"], "Input", AspectRatioFixed->True], Cell[TextData[ "logplot=ListPlot[logdata,PlotRange->{{0,6},{1.5,7.5}},\n\ AxesOrigin->{0,1.5}]"], "Input", AspectRatioFixed->True], Cell[TextData["I now calculate and plot the best-fitting line:"], "Text", Evaluatable->False, AspectRatioFixed->True], Cell[TextData["regline=Fit[logdata,{1,t},t]"], "Input", AspectRatioFixed->True], Cell[TextData[ "Show[logplot,Plot[regline,{t,1,6},\n\t\t\t\t\tDisplayFunction->Identity]]"], "Input", AspectRatioFixed->True], Cell[TextData[{ StyleBox["Notes: ", Evaluatable->False, AspectRatioFixed->True], StyleBox["Log[", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox["list", Evaluatable->False, AspectRatioFixed->True, FontSlant->"Italic"], StyleBox["]", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[" will return a list of the logarithms of the elements of ", Evaluatable->False, AspectRatioFixed->True], StyleBox["list", Evaluatable->False, AspectRatioFixed->True, FontSlant->"Italic"], StyleBox[", but in this case ", Evaluatable->False, AspectRatioFixed->True], StyleBox["N", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[ " is needed to convert integers to real numbers. Note the use of ", Evaluatable->False, AspectRatioFixed->True], StyleBox["Plot", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[" and ", Evaluatable->False, AspectRatioFixed->True], StyleBox["ListPlot", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[ " for plotting continuous curves and lists of data respectively, and of ", Evaluatable->False, AspectRatioFixed->True], StyleBox["Show", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[" for combining two or more graphics in a single plot. ", Evaluatable->False, AspectRatioFixed->True], StyleBox["Fit", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[" calculates regressions by least squares. In this case ", Evaluatable->False, AspectRatioFixed->True], StyleBox["{1,t}", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[" fits a linear regression, with an intercept (called for by ", Evaluatable->False, AspectRatioFixed->True], StyleBox["1", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[") and a slope (called for by ", Evaluatable->False, AspectRatioFixed->True], StyleBox["t", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[").", Evaluatable->False, AspectRatioFixed->True] }], "Text", Evaluatable->False, AspectRatioFixed->True]}, Open]], Cell[CellGroupData[{Cell[TextData["3"], "Subsection", Evaluatable->False, AspectRatioFixed->True], Cell[TextData[{ StyleBox["Solve the logistic equation as in Exercise 1b (giving the name ", Evaluatable->False, AspectRatioFixed->True], StyleBox["logisticsoln", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[ " to the rule for the solution and calling the initial value a and the \ parameters r and k) by writing ", Evaluatable->False, AspectRatioFixed->True], StyleBox["logisticsoln=DSolve[..]", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[". ", Evaluatable->False, AspectRatioFixed->True] }], "Text", Evaluatable->False, AspectRatioFixed->True], Cell[TextData[{ StyleBox["Note that ", Evaluatable->False, AspectRatioFixed->True], StyleBox["logisticsoln ", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox["is a replacement rule ", Evaluatable->False, AspectRatioFixed->True], StyleBox["n[t]\[Dash]>xxx", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox["; if we write ", Evaluatable->False, AspectRatioFixed->True], StyleBox["n[t]/.logisticsoln", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[", where ", Evaluatable->False, AspectRatioFixed->True], StyleBox["/. ", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[ "(pronounced \"slash dot\" and meaning \"evaluated at\") is the replacement \ operator, then ", Evaluatable->False, AspectRatioFixed->True], StyleBox["n[t]", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[" is temporarily replaced by ", Evaluatable->False, AspectRatioFixed->True], StyleBox["xxx", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"] }], "Text", Evaluatable->False, AspectRatioFixed->True], Cell[TextData["n[t]/.logisticsoln"], "Input", AspectRatioFixed->True], Cell[TextData[ "The solution with values of r, k and a temporarily assigned the values 0.1, \ 1000 and 1500 can be plotted thus"], "Text", Evaluatable->False, AspectRatioFixed->True], Cell[TextData[ "Plot[n[t]/.logisticsoln/.{r->0.1,k->1000,a->1500},\n\t\ {t,0,30},PlotRange->All]"], "Input", AspectRatioFixed->True]}, Open]], Cell[CellGroupData[{Cell[TextData["4"], "Subsection", Evaluatable->False, AspectRatioFixed->True], Cell[TextData[ "Choose parameter values giving the same carrying capacity and initial rate \ of increase for the two models and plot the curves. You can superimpose the \ two curves making one solid and the other dotted as in the following example: \ "], "Text", Evaluatable->False, AspectRatioFixed->True], Cell[TextData[ "Plot[{Sin[x],Cos[x]},{x,0,4Pi},\n\t\ PlotStyle->{{Thickness[0.004]},{Dashing[{.05,.05}]}}]\n\t\t"], "Input", AspectRatioFixed->True]}, Open]], Cell[CellGroupData[{Cell[TextData["5"], "Subsection", Evaluatable->False, AspectRatioFixed->True], Cell[CellGroupData[{Cell[TextData["5a"], "Subsubsection", Evaluatable->False, AspectRatioFixed->True], Cell[TextData[ "gausedata={2,10,10,11,21,56,104,137,165,194,217,\n\t\t\t\ 199,201,182,192,179,190,206,209,196,\n\t\t\t195,234,210,210,180};"], "Input", AspectRatioFixed->True], Cell[TextData[ "The days are not in exact sequence so that to plot the data against time it \ is necessary to compute a matrix {{0,2},{2,10},{3,10},...,{25,180}}:"], "Text",\ Evaluatable->False, AspectRatioFixed->True], Cell[TextData[ "pairsdata=Transpose[{Range[25],gausedata}];\npairsdata[[1,1]]=0;\n\ gauseplot=ListPlot[pairsdata]"], "Input", AspectRatioFixed->True]}, Open]], Cell[CellGroupData[{Cell[TextData["5b"], "Subsubsection", Evaluatable->False, AspectRatioFixed->True], Cell[TextData[{ StyleBox[ "We first show that (K\[Dash]n)/n is proportional to exp(\[Dash]rt). (I \ assume that ", Evaluatable->False, AspectRatioFixed->True], StyleBox["logisticsoln", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[ " is defined from Ex 3 as the replacement rule for the logistic curve. If \ not, you must define it by entering ", Evaluatable->False, AspectRatioFixed->True], StyleBox["logisticsoln=n[t]\[Dash]>xxx", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox["):", Evaluatable->False, AspectRatioFixed->True] }], "Text", Evaluatable->False, AspectRatioFixed->True], Cell[TextData["(k-n[t])/n[t]/.logisticsoln"], "Input", AspectRatioFixed->True], Cell[TextData["Simplify[%]"], "Input", AspectRatioFixed->True], Cell[TextData["We now estimate K and r:"], "Text", Evaluatable->False, AspectRatioFixed->True], Cell[TextData["k=Sum[gausedata[[i]],{i,10,25}]/16 //N"], "Input", AspectRatioFixed->True], Cell[TextData[ "data9=Take[gausedata,9]\n\ pairs9=Transpose[{Range[9],Log[(k\[Dash]data9)/data9]}]//N"], "Input", AspectRatioFixed->True], Cell[TextData["pairs9[[1,1]]=0;"], "Input", AspectRatioFixed->True], Cell[TextData["Fit[pairs9,{1,t},t]"], "Input", AspectRatioFixed->True], Cell[TextData[ "Assign r the value from the negative slope of this line, then plot the \ logistic curve and superimpose it on the data. (Since r and k have been \ assigned values, they should not be included in the temporary replacement \ rule in plotting the logistic curve.)"], "Text", Evaluatable->False, AspectRatioFixed->True]}, Open]]}, Open]], Cell[CellGroupData[{Cell[TextData["6"], "Subsection", Evaluatable->False, AspectRatioFixed->True], Cell[TextData["Adapt Exercise 5"], "Text", Evaluatable->False, AspectRatioFixed->True]}, Open]], Cell[CellGroupData[{Cell[TextData["7"], "Subsection", Evaluatable->False, AspectRatioFixed->True], Cell[TextData[{ StyleBox["Clear the existing values of ", Evaluatable->False, AspectRatioFixed->True], StyleBox["r", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[" and ", Evaluatable->False, AspectRatioFixed->True], StyleBox["k", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[ ", define and differentiate the function, simplify it since it looks rather \ complicated, and then evaluate it at the equilibrium values. ", Evaluatable->False, AspectRatioFixed->True], StyleBox["%", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[ " denotes the output from the previous calculation; this notation should be \ used sparingly.", Evaluatable->False, AspectRatioFixed->True] }], "Text", Evaluatable->False, AspectRatioFixed->True], Cell[TextData["Clear[r,k]\nf[n_]:=r (1-n/k) n"], "Input", AspectRatioFixed->True], Cell[TextData["D[f[n],n]"], "Input", AspectRatioFixed->True], Cell[TextData["Simplify[%]"], "Input", AspectRatioFixed->True], Cell[TextData["%/.{{n->0},{n->k}}"], "Input", AspectRatioFixed->True]}, Open]], Cell[CellGroupData[{Cell[TextData["8"], "Subsection", Evaluatable->False, AspectRatioFixed->True], Cell[TextData["Adapt Exercise 7."], "Text", Evaluatable->False, AspectRatioFixed->True]}, Open]], Cell[CellGroupData[{Cell[TextData["9"], "Subsection", Evaluatable->False, AspectRatioFixed->True], Cell[TextData["F[n_]:=n Exp[r (1-n/k)]"], "Input", AspectRatioFixed->True], Cell[TextData[ "Evaluate this function at 0 and k, and then evaluate its derivatives at \ these points."], "Text", Evaluatable->False, AspectRatioFixed->True]}, Open]], Cell[CellGroupData[{Cell[TextData["10"], "Subsection", Evaluatable->False, AspectRatioFixed->True], Cell[TextData[{ StyleBox["The operator ", Evaluatable->False, AspectRatioFixed->True], StyleBox["NestList[]", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[" provides a powerful method of iteration. For example", Evaluatable->False, AspectRatioFixed->True] }], "Text", Evaluatable->False, AspectRatioFixed->True], Cell[TextData["r=2.5;\nk=1000;\nNestList[F,800,10]"], "Input", AspectRatioFixed->True], Cell[TextData[{ StyleBox[ "returns a list of the values obtained by iterating F(n) 10 times, starting \ with the initial value 800. (Thus there are 11 values altogether.) The \ parameters in the function must be assigned values before ", Evaluatable->False, AspectRatioFixed->True], StyleBox["NestList[]", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[" is used. This list of values can be passed for plotting to ", Evaluatable->False, AspectRatioFixed->True], StyleBox["ListPlot", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[":", Evaluatable->False, AspectRatioFixed->True] }], "Text", Evaluatable->False, AspectRatioFixed->True], Cell[TextData["ListPlot[NestList[F,800,10],PlotJoined->True]"], "Input", AspectRatioFixed->True], Cell[TextData[{ StyleBox["The graphics option ", Evaluatable->False, AspectRatioFixed->True], StyleBox["PlotJoined\[Dash]>True", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[ " joins the points by lines. Some experimentation with further graphics \ options such as ", Evaluatable->False, AspectRatioFixed->True], StyleBox["AxesOrigin", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[" and ", Evaluatable->False, AspectRatioFixed->True], StyleBox["PlotRange", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[" may be needed to make the graph look pretty.", Evaluatable->False, AspectRatioFixed->True] }], "Text", Evaluatable->False, AspectRatioFixed->True], Cell[TextData[ " Use this method to investigate the behavior of the iterates of \ F(n) with different values of r and k and with different initial values."], "Text", Evaluatable->False, AspectRatioFixed->True]}, Open]], Cell[CellGroupData[{Cell[TextData["11"], "Subsection", Evaluatable->False, AspectRatioFixed->True], Cell[TextData[{ StyleBox["First enter the data, ", Evaluatable->False, AspectRatioFixed->True], StyleBox["data={. , ... , .}", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[". Then calculate", Evaluatable->False, AspectRatioFixed->True] }], "Text", Evaluatable->False, AspectRatioFixed->True], Cell[TextData[ "logdata=N[Log[data]];\nvar1=Rest[logdata];\nvar2=Drop[logdata,-1];\n\ pairsdata=Transpose[{var2,var1}];"], "Input", AspectRatioFixed->True], Cell[TextData[{ StyleBox["var1", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[" is the list of log numbers with its first element dropped, ", Evaluatable->False, AspectRatioFixed->True], StyleBox["var2", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[" is the list with its last element dropped, and ", Evaluatable->False, AspectRatioFixed->True], StyleBox["pairsdata ", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox["is the list of pairs of values {{x", Evaluatable->False, AspectRatioFixed->True], StyleBox["1", Evaluatable->False, AspectRatioFixed->True, FontVariations->{"CompatibilityType"->"Subscript"}], StyleBox[",x", Evaluatable->False, AspectRatioFixed->True], StyleBox["2", Evaluatable->False, AspectRatioFixed->True, FontVariations->{"CompatibilityType"->"Subscript"}], StyleBox["},{x", Evaluatable->False, AspectRatioFixed->True], StyleBox["2", Evaluatable->False, AspectRatioFixed->True, FontVariations->{"CompatibilityType"->"Subscript"}], StyleBox[",x", Evaluatable->False, AspectRatioFixed->True], StyleBox["3", Evaluatable->False, AspectRatioFixed->True, FontVariations->{"CompatibilityType"->"Subscript"}], StyleBox["},...,{x", Evaluatable->False, AspectRatioFixed->True], StyleBox["n\[Dash]1", Evaluatable->False, AspectRatioFixed->True, FontVariations->{"CompatibilityType"->"Subscript"}], StyleBox[",x", Evaluatable->False, AspectRatioFixed->True], StyleBox["n", Evaluatable->False, AspectRatioFixed->True, FontVariations->{"CompatibilityType"->"Subscript"}], StyleBox["}}, where x", Evaluatable->False, AspectRatioFixed->True], StyleBox["t ", Evaluatable->False, AspectRatioFixed->True, FontVariations->{"CompatibilityType"->"Subscript"}], StyleBox["= log n", Evaluatable->False, AspectRatioFixed->True], StyleBox["t", Evaluatable->False, AspectRatioFixed->True, FontVariations->{"CompatibilityType"->"Subscript"}], StyleBox[". Thus ", Evaluatable->False, AspectRatioFixed->True], StyleBox["ListPlot[pairsdata]", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[" will give the required plot, and ", Evaluatable->False, AspectRatioFixed->True], StyleBox["Fit[pairsdata,{1,x},x]", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[" will give the ordinary least squares fit.", Evaluatable->False, AspectRatioFixed->True] }], "Text", Evaluatable->False, AspectRatioFixed->True], Cell[TextData[{ StyleBox[" The package ", Evaluatable->False, AspectRatioFixed->True], StyleBox["MonteCarlo", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[" contains a procedure called ", Evaluatable->False, AspectRatioFixed->True], StyleBox["Pollard", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[ " which will implement the test of Pollard et al. (1987) for density \ dependence. You must first read in the package by using the command ", Evaluatable->False, AspectRatioFixed->True], StyleBox["<False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[ ". (This may not work if the package is not in the same folder as the ", Evaluatable->False, AspectRatioFixed->True], StyleBox["Mathematica", Evaluatable->False, AspectRatioFixed->True, FontSlant->"Italic"], StyleBox[" program.) The command ", Evaluatable->False, AspectRatioFixed->True], StyleBox["Pollard[data,99] ", Evaluatable->False, AspectRatioFixed->True, FontFamily->"Courier", FontWeight->"Bold"], StyleBox[ "will then calculate 99 Monte Carlo simulations and return the rank of the \ observed slope among the 100 values (99 simulated + 1 observed); ", Evaluatable->False, AspectRatioFixed->True], StyleBox["?Pollard", Evaluatable->False, AspectRatioFixed->True, FontWeight->"Bold"], StyleBox[" gives information about the procedure.", Evaluatable->False, AspectRatioFixed->True] }], "Text", Evaluatable->False, AspectRatioFixed->True], Cell[TextData["<True], Cell[TextData["?Pollard"], "Input", AspectRatioFixed->True], Cell[TextData["Pollard[data,99]"], "Input", AspectRatioFixed->True]}, Open]]}, Open]] }, FrontEndVersion->"Macintosh 3.0", ScreenRectangle->{{0, 640}, {0, 460}}, WindowToolbars->{}, CellGrouping->Manual, WindowSize->{520, 365}, WindowMargins->{{44, Automatic}, {30, Automatic}}, PrivateNotebookOptions->{"ColorPalette"->{RGBColor, -1}}, ShowCellLabel->True, ShowCellTags->False, RenderingOptions->{"ObjectDithering"->True, "RasterDithering"->False}, MacintoshSystemPageSetup->"\<\ AVU/IFiQKFD000000V:^/09R]g0000000OVaH097bCP0AP1Y06`0I@1^0642HZj` 0V:gT0000001nK500TO9>000000000000000009R[[0000000000000000000000 00000000000000000000000000000000\>" ] (*********************************************************************** 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, 86, 2, 70, "Section", Evaluatable->False], Cell[CellGroupData[{ Cell[1840, 55, 81, 2, 70, "Subsection", Evaluatable->False], Cell[CellGroupData[{ Cell[1944, 59, 85, 2, 70, "Subsubsection", Evaluatable->False], Cell[2032, 63, 596, 24, 70, "Text", Evaluatable->False], Cell[2631, 89, 70, 1, 70, "Input"], Cell[2704, 92, 62, 1, 70, "Input"], Cell[2769, 95, 552, 20, 70, "Text", Evaluatable->False], Cell[3324, 117, 59, 1, 70, "Input"], Cell[3386, 120, 57, 1, 70, "Input"] }, Closed]], Cell[CellGroupData[{ Cell[3475, 123, 85, 2, 70, "Subsubsection", Evaluatable->False], Cell[3563, 127, 490, 19, 70, "Text", Evaluatable->False], Cell[4056, 148, 61, 1, 70, "Input"], Cell[4120, 151, 91, 1, 70, "Input"] }, Closed]] }, Closed]], Cell[CellGroupData[{ Cell[4252, 154, 81, 2, 70, "Subsection", Evaluatable->False], Cell[4336, 158, 132, 3, 70, "Text", Evaluatable->False], Cell[4471, 163, 92, 1, 70, "Input"], Cell[4566, 166, 132, 3, 70, "Input"], Cell[4701, 171, 121, 2, 70, "Text", Evaluatable->False], Cell[4825, 175, 81, 1, 70, "Input"], Cell[4909, 178, 130, 3, 70, "Input"], Cell[5042, 183, 2397, 86, 70, "Text", Evaluatable->False] }, Closed]], Cell[CellGroupData[{ Cell[7471, 271, 81, 2, 70, "Subsection", Evaluatable->False], Cell[7555, 275, 668, 22, 70, "Text", Evaluatable->False], Cell[8226, 299, 1267, 47, 70, "Text", Evaluatable->False], Cell[9496, 348, 71, 1, 70, "Input"], Cell[9570, 351, 186, 4, 70, "Text", Evaluatable->False], Cell[9759, 357, 134, 3, 70, "Input"] }, Closed]], Cell[CellGroupData[{ Cell[9925, 362, 81, 2, 70, "Subsection", Evaluatable->False], Cell[10009, 366, 310, 6, 70, "Text", Evaluatable->False], Cell[10322, 374, 150, 3, 70, "Input"] }, Closed]], Cell[CellGroupData[{ Cell[10504, 379, 81, 2, 70, "Subsection", Evaluatable->False], Cell[CellGroupData[{ Cell[10608, 383, 85, 2, 70, "Subsubsection", Evaluatable->False], Cell[10696, 387, 177, 3, 70, "Input"], Cell[10876, 392, 223, 5, 70, "Text", Evaluatable->False], Cell[11102, 399, 151, 3, 70, "Input"] }, Closed]], Cell[CellGroupData[{ Cell[11285, 404, 85, 2, 70, "Subsubsection", Evaluatable->False], Cell[11373, 408, 706, 24, 70, "Text", Evaluatable->False], Cell[12082, 434, 80, 1, 70, "Input"], Cell[12165, 437, 64, 1, 70, "Input"], Cell[12232, 440, 98, 2, 70, "Text", Evaluatable->False], Cell[12333, 444, 93, 1, 70, "Input"], Cell[12429, 447, 139, 3, 70, "Input"], Cell[12571, 452, 69, 1, 70, "Input"], Cell[12643, 455, 72, 1, 70, "Input"], Cell[12718, 458, 335, 6, 70, "Text", Evaluatable->False] }, Closed]] }, Closed]], Cell[CellGroupData[{ Cell[13094, 466, 81, 2, 70, "Subsection", Evaluatable->False], Cell[13178, 470, 90, 2, 70, "Text", Evaluatable->False] }, Closed]], Cell[CellGroupData[{ Cell[13300, 474, 81, 2, 70, "Subsection", Evaluatable->False], Cell[13384, 478, 895, 31, 70, "Text", Evaluatable->False], Cell[14282, 511, 83, 1, 70, "Input"], Cell[14368, 514, 62, 1, 70, "Input"], Cell[14433, 517, 64, 1, 70, "Input"], Cell[14500, 520, 71, 1, 70, "Input"] }, Closed]], Cell[CellGroupData[{ Cell[14603, 523, 81, 2, 70, "Subsection", Evaluatable->False], Cell[14687, 527, 91, 2, 70, "Text", Evaluatable->False] }, Closed]], Cell[CellGroupData[{ Cell[14810, 531, 81, 2, 70, "Subsection", Evaluatable->False], Cell[14894, 535, 76, 1, 70, "Input"], Cell[14973, 538, 162, 4, 70, "Text", Evaluatable->False] }, Closed]], Cell[CellGroupData[{ Cell[15167, 544, 82, 2, 70, "Subsection", Evaluatable->False], Cell[15252, 548, 378, 13, 70, "Text", Evaluatable->False], Cell[15633, 563, 88, 1, 70, "Input"], Cell[15724, 566, 753, 23, 70, "Text", Evaluatable->False], Cell[16480, 591, 98, 1, 70, "Input"], Cell[16581, 594, 825, 29, 70, "Text", Evaluatable->False], Cell[17409, 625, 225, 5, 70, "Text", Evaluatable->False] }, Closed]], Cell[CellGroupData[{ Cell[17666, 632, 82, 2, 70, "Subsection", Evaluatable->False], Cell[17751, 636, 358, 13, 70, "Text", Evaluatable->False], Cell[18112, 651, 157, 3, 70, "Input"], Cell[18272, 656, 2747, 94, 70, "Text", Evaluatable->False], Cell[21022, 752, 1662, 54, 70, "Text", Evaluatable->False], Cell[22687, 808, 67, 1, 70, "Input"], Cell[22757, 811, 61, 1, 70, "Input"], Cell[22821, 814, 69, 1, 70, "Input"] }, Closed]] }, Open ]] } ] *) (*********************************************************************** End of Mathematica Notebook file. ***********************************************************************)