(*********************************************************************** Mathematica-Compatible Notebook This notebook can be used on any computer system with Mathematica 4.0, MathReader 4.0, or any compatible application. 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[ 53392, 1510]*) (*NotebookOutlinePosition[ 54277, 1538]*) (* CellTagsIndexPosition[ 54233, 1534]*) (*WindowFrame->Normal*) Notebook[{ Cell[CellGroupData[{ Cell["Business cycle tools", "Subtitle"], Cell[TextData[{ "Luci Ellis\nReserve Bank of Australia and Australian National University\n \ ", ButtonBox["elisha@dot.net.au", ButtonData:>{ URL[ "mailto:elisha@dot.net.au"], None}, Active->True, ButtonStyle->"Hyperlink", ButtonNote->"Send mail to Luci Ellis at home"] }], "Author"], Cell[CellGroupData[{ Cell["Finding peaks in data", "Section"], Cell[BoxData[ \(Clear[findPeaks]\)], "Input", CellLabel->"In[3]:="], Cell[CellGroupData[{ Cell[BoxData[ \(findPeaks::usage = "\"\)], "Input", CellLabel->"In[4]:="], Cell[BoxData[ \("findPeaks[data,consec] finds the points in data that are an increase \ from the previous point, and following which there are consec consecutive \ declines. This function can be used to identify recessions in macroeconomic \ data. The parameter consec is optional. If omitted, it takes the value 2. \ findPeaks[data,{consecups,consecdowns}] finds the points in data that \ complete a sequence of consecups consecutive rises and following which there \ are consecdowns consecutive falls."\)], "Output", CellLabel->"Out[4]="] }, Open ]], Cell[BoxData[ \(\(Options[findPeaks] = {DataInLevels \[Rule] True};\)\)], "Input", CellLabel->"In[5]:="], Cell[BoxData[ \(findPeaks[data : {__?NumericQ}, consec_: 2, opts___Rule] /; \ \((Head[consec] \[Equal] Integer && Length[data] > 4\ && \ consec > 0\ && Length[data] > consec + 2)\) := \[IndentingNewLine]Module[{difffirst, ddata, \ nested, \ negtest, posns}, \ \[IndentingNewLine]difffirst\ = \ \(DataInLevels /. \ {opts}\) /. Options[findPeaks]; \[IndentingNewLine]ddata = If[difffirst, ListCorrelate[{\(-1\), 1}, data], data]; \ (*\ differencing\ data\ if\ required\ *) \[IndentingNewLine]nested = Transpose[\ \(Drop[#, \(-consec\)] &\)\ /@ \ \((\(RotateRight[ ddata, #] &\)\ /@ \ Range[consec - 1, \(-1\), \(-1\)])\)]; \[IndentingNewLine]negtest = Join[{_?Positive}, Table[_?Negative, {consec}]]; \[IndentingNewLine]posns = Cases[Flatten[ Position[nested, negtest, {1}, Heads \[Rule] False]] - \((consec - If[difffirst, 2, 1])\), _? Positive]; \[IndentingNewLine]Transpose[{data\ \[LeftDoubleBracket]posns\[RightDoubleBracket], posns}]]\)], "Input", CellLabel->"In[6]:="], Cell[BoxData[ \(findPeaks[ data : {__?NumericQ}, {consecups_Integer?Positive, consecdowns_Integer?Positive}, opts___Rule] /; \ \((Length[data] > 4\ \ && Length[data] > consecups + consecdowns + 2)\) := \[IndentingNewLine]Module[{difffirst, ddata, \ nested, \ negtest, posns}, \ \[IndentingNewLine]difffirst\ = \ \(DataInLevels /. \ {opts}\) /. Options[findPeaks]; \[IndentingNewLine]ddata = If[difffirst, ListCorrelate[{\(-1\), 1}, data], data]; \ (*\ differencing\ data\ if\ required\ *) \[IndentingNewLine]nested = Transpose[\ \(Drop[#, \(-consecups\) - consecdowns] &\)\ /@ \ \((\(RotateRight[ ddata, #] &\)\ /@ \ Range[consecdowns - 1, \(-consecups\), \(-1\)])\)]; \ \[IndentingNewLine]negtest = Join[Table[_?Positive, {consecups}], Table[_?Negative, {consecdowns}]]; \[IndentingNewLine]posns = Cases[Flatten[ Position[nested, negtest, {1}, Heads \[Rule] False]] + consecups - 1 - \((\(+consecdowns\) - If[difffirst, 2, 1])\), _? Positive]; \[IndentingNewLine]Transpose[{data\ \[LeftDoubleBracket]posns\[RightDoubleBracket], posns}]]\)], "Input", CellLabel->"In[49]:="], Cell[CellGroupData[{ Cell["Error cases", "Subsubsection"], Cell[CellGroupData[{ Cell[BoxData[ \(findPeaks::negcon = "\"\)], "Input", CellLabel->"In[8]:="], Cell[BoxData[ \("The second parameter to this function must be positive."\)], "Output", CellLabel->"Out[8]="] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(findPeaks::conlen\ = \ "\"\)], "Input", CellLabel->"In[9]:="], Cell[BoxData[ \("The number of consecutive falls to be found `1` is larger than the \ length of the data series, `2`."\)], "Output", CellLabel->"Out[9]="] }, Open ]], Cell[BoxData[ \(findPeaks[data : {__?NumericQ}, consec_, opts___Rule] /; consec < 1 := Message[findPeaks::negcon]\)], "Input", CellLabel->"In[10]:="], Cell[BoxData[ \(findPeaks[data : {__?NumericQ}, consec_, opts___Rule] /; \((consec + 2)\) > Length[data] := Message[findPeaks::conlen, consec, Length[data]]\)], "Input", CellLabel->"In[11]:="] }, Closed]], Cell[CellGroupData[{ Cell["Demonstration", "Subsubsection"], Cell[BoxData[ \(\(testdata = NestList[#\ 0.95\ - 0.5\ + Random[] &, 0.3, 200];\)\)], "Input", CellLabel->"In[12]:="], Cell[BoxData[ \(Show[ ListPlot[Reverse /@ \ findPeaks[testdata, 3], Frame \[Rule] True, PlotStyle \[Rule] {Hue[0], AbsolutePointSize[6]}, DisplayFunction \[Rule] Identity], ListPlot[testdata, PlotJoined \[Rule] True, DisplayFunction \[Rule] Identity], DisplayFunction \[Rule] $DisplayFunction]\)], "Input", CellLabel->"In[57]:="] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Assymmetric smoothing", "Section"], Cell[TextData[{ "Reference: Julian Allwood and ", ButtonBox["David Shepherd", ButtonData:>{ URL[ "mailto:d.shepherd@ecomfac.unimelb.edu.au"], None}, Active->True, ButtonStyle->"Hyperlink", ButtonNote->"Send mail to David Shepherd"], " (1999), \"Alternative Detrending Procedures for Macroeconomic Time Series\ \", University of Melbourne Department of Economics Research Paper Number \ 698, June 1999" }], "Text"], Cell[BoxData[ \(Needs["\"]\)], "Input", CellLabel->"In[59]:="], Cell[BoxData[ \(assymetricDetrending[data : {__?NumericQ}, consec_Integer?Positive] /; \((\ Head[consec] \[Equal] Integer && Length[data] > 4\ && \ consec > 0\ && Length[data] > consec + 2)\) := \[IndentingNewLine]With[{points = findPeaks[data, consec]}, SplineFit[points, Cubic]]\)], "Input", CellLabel->"In[70]:="], Cell[BoxData[ \(assymetricDetrending[ data : {__?NumericQ}, {consecups_Integer?Positive, consecdowns_Integer?Positive}]\ /; \((\ Length[data] > 4\ && Length[data] > consecups\ + consecdowns + 2)\) := \[IndentingNewLine]With[{points = findPeaks[data, {consecups, consecdowns}]}, SplineFit[points, Cubic]]\)], "Input", CellLabel->"In[77]:="], Cell[CellGroupData[{ Cell["Demonstration", "Subsubsection"], Cell[CellGroupData[{ Cell[BoxData[ \(assDTtest = assymetricDetrending[testdata, 2]\)], "Input", CellLabel->"In[72]:="], Cell[BoxData[ InterpretationBox[\("SplineFunction["\[InvisibleSpace]Cubic\ \[InvisibleSpace]", "\[InvisibleSpace]{0.`, 21.`}\[InvisibleSpace]", <>]"\), SequenceForm[ "SplineFunction[", Cubic, ", ", {0.0, 21.0}, ", <>]"], Editable->False]], "Output", CellLabel->"Out[72]="] }, Open ]], Cell["\<\ Note that you different time series result in SplineFunctions with \ different ranges, so you may have to tweak the number in blue in the next \ cell.\ \>", "Text"], Cell[BoxData[ RowBox[{"Show", "[", RowBox[{\(ListPlot[Reverse /@ \ findPeaks[testdata, 2], PlotStyle \[Rule] {Hue[0], AbsolutePointSize[6]}, DisplayFunction \[Rule] Identity]\), ",", \(ListPlot[testdata, PlotJoined \[Rule] True, DisplayFunction \[Rule] Identity]\), ",", RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"Reverse", "/@", RowBox[{"Table", "[", RowBox[{\(assDTtest[x]\), ",", RowBox[{"{", RowBox[{"x", ",", "0", ",", StyleBox["21", FontColor->RGBColor[0, 0, 1]], ",", "0.2"}], "}"}]}], "]"}]}], ",", \(PlotJoined \[Rule] True\), ",", \(PlotStyle \[Rule] Hue[0.6]\), ",", \(DisplayFunction \[Rule] Identity\)}], "]"}], ",", \(DisplayFunction \[Rule] $DisplayFunction\)}], "]"}]], "Input", CellLabel->"In[73]:="], Cell["This suggests a further possible function:", "Text"], Cell[BoxData[ \(plotAssymetricDetrendedData[data : {__?NumericQ}, consec_: Integer] /; \((\ Head[consec] \[Equal] Integer && Length[data] > 4\ && \ consec > 0\ && Length[data] > consec + 2)\) := \[IndentingNewLine]With[{spliner = assymetricDetrending[data, consec], dots = Reverse /@ \ findPeaks[data, consec], T = Length[data]}, \[IndentingNewLine]With[{splinedots = Reverse /@ Table[spliner[x], Evaluate[ Join[{x}, spliner\[LeftDoubleBracket]2\[RightDoubleBracket], \ {spliner\[LeftDoubleBracket]2, 2\[RightDoubleBracket]/ T}]]]}, \[IndentingNewLine]Show[ ListPlot[dots, PlotStyle \[Rule] {Hue[0], AbsolutePointSize[6]}, Frame \[Rule] True, DisplayFunction \[Rule] Identity], ListPlot[data, PlotJoined \[Rule] True, Frame \[Rule] True, DisplayFunction \[Rule] Identity], ListPlot[splinedots, PlotJoined \[Rule] True, PlotStyle \[Rule] Hue[0.6], Frame \[Rule] True, DisplayFunction \[Rule] Identity], DisplayFunction \[Rule] $DisplayFunction]\ ]]\)], "Input", CellLabel->"In[83]:="], Cell[BoxData[ \(plotAssymetricDetrendedData[ data : {__?NumericQ}, {consecups_Integer?Positive, consecdowns_Integer?Positive}] /; \((\ Length[data] > 4\ && \ Length[data] > consecups + consecdowns + 2)\) := \[IndentingNewLine]With[{spliner = assymetricDetrending[data, {consecups, consecdowns}], dots = Reverse /@ \ findPeaks[data, {consecups, consecdowns}], T = Length[data]}, \[IndentingNewLine]With[{splinedots = Reverse /@ Table[spliner[x], Evaluate[ Join[{x}, spliner\[LeftDoubleBracket]2\[RightDoubleBracket], \ {spliner\[LeftDoubleBracket]2, 2\[RightDoubleBracket]/ T}]]]}, \[IndentingNewLine]Show[ ListPlot[dots, PlotStyle \[Rule] {Hue[0], AbsolutePointSize[6]}, Frame \[Rule] True, DisplayFunction \[Rule] Identity], ListPlot[data, PlotJoined \[Rule] True, Frame \[Rule] True, DisplayFunction \[Rule] Identity], ListPlot[splinedots, PlotJoined \[Rule] True, PlotStyle \[Rule] Hue[0.6], Frame \[Rule] True, DisplayFunction \[Rule] Identity], DisplayFunction \[Rule] $DisplayFunction]\ ]]\)], "Input", CellLabel->"In[84]:="], Cell[CellGroupData[{ Cell[BoxData[ \(plotAssymetricDetrendedData[testdata, {2, 2}] // Timing\)], "Input", CellLabel->"In[87]:="], 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.00473821 0.275249 0.143564 [ [.02381 -0.0125 -3 -9 ] [.02381 -0.0125 3 0 ] [.26072 -0.0125 -6 -9 ] [.26072 -0.0125 6 0 ] [.49763 -0.0125 -9 -9 ] [.49763 -0.0125 9 0 ] [.73454 -0.0125 -9 -9 ] [.73454 -0.0125 9 0 ] [.97145 -0.0125 -9 -9 ] [.97145 -0.0125 9 0 ] [ 0 0 -0.125 0 ] [-0.0125 .13168 -12 -4.5 ] [-0.0125 .13168 0 4.5 ] [-0.0125 .27525 -6 -4.5 ] [-0.0125 .27525 0 4.5 ] [-0.0125 .41881 -6 -4.5 ] [-0.0125 .41881 0 4.5 ] [-0.0125 .56238 -6 -4.5 ] [-0.0125 .56238 0 4.5 ] [ 0 0 -0.125 0 ] [ 0 .61803 .125 0 ] [ 1 0 .125 0 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .02381 0 m .02381 .00625 L s [(0)] .02381 -0.0125 0 1 Mshowa .26072 0 m .26072 .00625 L s [(50)] .26072 -0.0125 0 1 Mshowa .49763 0 m .49763 .00625 L s [(100)] .49763 -0.0125 0 1 Mshowa .73454 0 m .73454 .00625 L s [(150)] .73454 -0.0125 0 1 Mshowa .97145 0 m .97145 .00625 L s [(200)] .97145 -0.0125 0 1 Mshowa .125 Mabswid .07119 0 m .07119 .00375 L s .11857 0 m .11857 .00375 L s .16596 0 m .16596 .00375 L s .21334 0 m .21334 .00375 L s .3081 0 m .3081 .00375 L s .35548 0 m .35548 .00375 L s .40287 0 m .40287 .00375 L s .45025 0 m .45025 .00375 L s .54501 0 m .54501 .00375 L s .5924 0 m .5924 .00375 L s .63978 0 m .63978 .00375 L s .68716 0 m .68716 .00375 L s .78192 0 m .78192 .00375 L s .82931 0 m .82931 .00375 L s .87669 0 m .87669 .00375 L s .92407 0 m .92407 .00375 L s .25 Mabswid 0 0 m 1 0 L s 0 .13168 m .00625 .13168 L s [(-1)] -0.0125 .13168 1 0 Mshowa 0 .27525 m .00625 .27525 L s [(0)] -0.0125 .27525 1 0 Mshowa 0 .41881 m .00625 .41881 L s [(1)] -0.0125 .41881 1 0 Mshowa 0 .56238 m .00625 .56238 L s [(2)] -0.0125 .56238 1 0 Mshowa .125 Mabswid 0 .1604 m .00375 .1604 L s 0 .18911 m .00375 .18911 L s 0 .21782 m .00375 .21782 L s 0 .24654 m .00375 .24654 L s 0 .30396 m .00375 .30396 L s 0 .33267 m .00375 .33267 L s 0 .36139 m .00375 .36139 L s 0 .3901 m .00375 .3901 L s 0 .44753 m .00375 .44753 L s 0 .47624 m .00375 .47624 L s 0 .50495 m .00375 .50495 L s 0 .53366 m .00375 .53366 L s 0 .10297 m .00375 .10297 L s 0 .07426 m .00375 .07426 L s 0 .04555 m .00375 .04555 L s 0 .01683 m .00375 .01683 L s 0 .59109 m .00375 .59109 L s .25 Mabswid 0 0 m 0 .61803 L s .02381 .61178 m .02381 .61803 L s .26072 .61178 m .26072 .61803 L s .49763 .61178 m .49763 .61803 L s .73454 .61178 m .73454 .61803 L s .97145 .61178 m .97145 .61803 L s .125 Mabswid .07119 .61428 m .07119 .61803 L s .11857 .61428 m .11857 .61803 L s .16596 .61428 m .16596 .61803 L s .21334 .61428 m .21334 .61803 L s .3081 .61428 m .3081 .61803 L s .35548 .61428 m .35548 .61803 L s .40287 .61428 m .40287 .61803 L s .45025 .61428 m .45025 .61803 L s .54501 .61428 m .54501 .61803 L s .5924 .61428 m .5924 .61803 L s .63978 .61428 m .63978 .61803 L s .68716 .61428 m .68716 .61803 L s .78192 .61428 m .78192 .61803 L s .82931 .61428 m .82931 .61803 L s .87669 .61428 m .87669 .61803 L s .92407 .61428 m .92407 .61803 L s .25 Mabswid 0 .61803 m 1 .61803 L s .99375 .13168 m 1 .13168 L s .99375 .27525 m 1 .27525 L s .99375 .41881 m 1 .41881 L s .99375 .56238 m 1 .56238 L s .125 Mabswid .99625 .1604 m 1 .1604 L s .99625 .18911 m 1 .18911 L s .99625 .21782 m 1 .21782 L s .99625 .24654 m 1 .24654 L s .99625 .30396 m 1 .30396 L s .99625 .33267 m 1 .33267 L s .99625 .36139 m 1 .36139 L s .99625 .3901 m 1 .3901 L s .99625 .44753 m 1 .44753 L s .99625 .47624 m 1 .47624 L s .99625 .50495 m 1 .50495 L s .99625 .53366 m 1 .53366 L s .99625 .10297 m 1 .10297 L s .99625 .07426 m 1 .07426 L s .99625 .04555 m 1 .04555 L s .99625 .01683 m 1 .01683 L s .99625 .59109 m 1 .59109 L s .25 Mabswid 1 0 m 1 .61803 L s 0 .27525 m 1 .27525 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath 1 0 0 r 6 Mabswid .08067 .50099 Mdot .09962 .56601 Mdot .12805 .56029 Mdot .17069 .2477 Mdot .2086 .30956 Mdot .35075 .25303 Mdot .41234 .39396 Mdot .47394 .26072 Mdot .54975 .45443 Mdot .71085 .3779 Mdot .73454 .3699 Mdot .80088 .3298 Mdot .94776 .3263 Mdot 0 g .5 Mabswid .02855 .31832 m .03329 .31561 L .03802 .37828 L .04276 .36552 L .0475 .37101 L .05224 .39645 L .05698 .39801 L .06172 .438 L .06645 .39742 L .07119 .41498 L .07593 .47216 L .08067 .50099 L .08541 .49842 L .09014 .44654 L .09488 .50631 L .09962 .56601 L .10436 .54144 L .1091 .48724 L .11384 .48792 L .11857 .53302 L .12331 .54437 L .12805 .56029 L .13279 .48078 L .13753 .46493 L .14226 .40542 L .147 .35706 L .15174 .2948 L .15648 .22409 L .16122 .23839 L .16596 .23994 L .17069 .2477 L .17543 .20499 L .18017 .18091 L .18491 .23345 L .18965 .24914 L .19439 .23948 L .19912 .27673 L .20386 .29341 L .2086 .30956 L .21334 .2974 L .21808 .27001 L .22281 .25614 L .22755 .19109 L .23229 .2009 L .23703 .22588 L .24177 .25661 L .24651 .32859 L .25124 .2653 L .25598 .26531 L .26072 .26414 L Mistroke .26546 .32213 L .2702 .25713 L .27493 .31758 L .27967 .38587 L .28441 .36 L .28915 .40896 L .29389 .41764 L .29863 .36132 L .30336 .39486 L .3081 .38127 L .31284 .42765 L .31758 .42585 L .32232 .36732 L .32706 .29378 L .33179 .29903 L .33653 .24602 L .34127 .1907 L .34601 .19217 L .35075 .25303 L .35548 .24583 L .36022 .19125 L .36496 .18153 L .3697 .17163 L .37444 .24196 L .37918 .20191 L .38391 .18254 L .38865 .18482 L .39339 .25586 L .39813 .30725 L .40287 .30119 L .4076 .35395 L .41234 .39396 L .41708 .32882 L .42182 .25751 L .42656 .18899 L .4313 .15964 L .43603 .229 L .44077 .15983 L .44551 .1618 L .45025 .11385 L .45499 .1705 L .45973 .1788 L .46446 .1949 L .4692 .22538 L .47394 .26072 L .47868 .20665 L .48342 .16416 L .48815 .20281 L .49289 .20921 L .49763 .25253 L Mistroke .50237 .28189 L .50711 .31992 L .51185 .34183 L .51658 .31119 L .52132 .35534 L .52606 .41754 L .5308 .42199 L .53554 .38419 L .54027 .42231 L .54501 .44305 L .54975 .45443 L .55449 .4147 L .55923 .38332 L .56397 .38785 L .5687 .42519 L .57344 .45564 L .57818 .48272 L .58292 .5276 L .58766 .45332 L .5924 .46538 L .59713 .46611 L .60187 .40694 L .60661 .40769 L .61135 .45381 L .61609 .49867 L .62082 .47129 L .62556 .51647 L .6303 .51267 L .63504 .43687 L .63978 .41816 L .64452 .38263 L .64925 .34422 L .65399 .30504 L .65873 .33662 L .66347 .35718 L .66821 .42259 L .67294 .47569 L .67768 .41702 L .68242 .46237 L .68716 .4128 L .6919 .3585 L .69664 .32222 L .70137 .36219 L .70611 .36844 L .71085 .3779 L .71559 .3185 L .72033 .27954 L .72507 .30897 L .7298 .33939 L .73454 .3699 L Mistroke .73928 .34518 L .74402 .29127 L .74876 .31474 L .75349 .28533 L .75823 .22144 L .76297 .27854 L .76771 .26659 L .77245 .2783 L .77719 .26291 L .78192 .32022 L .78666 .3175 L .7914 .30354 L .79614 .30622 L .80088 .3298 L .80561 .30225 L .81035 .24938 L .81509 .28422 L .81983 .22653 L .82457 .2618 L .82931 .19344 L .83404 .19609 L .83878 .18496 L .84352 .11841 L .84826 .09528 L .853 .05105 L .85774 .02579 L .86247 .01472 L .86721 .09599 L .87195 .04333 L .87669 .03583 L .88143 .10782 L .88616 .10138 L .8909 .13547 L .89564 .13845 L .90038 .13399 L .90512 .20988 L .90986 .16268 L .91459 .21095 L .91933 .1559 L .92407 .21044 L .92881 .20143 L .93355 .23322 L .93828 .21601 L .94302 .2648 L .94776 .3263 L .9525 .29515 L .95724 .27413 L .96198 .27919 L .96671 .32142 L .97145 .25521 L Mistroke .97619 .18795 L Mfstroke 0 .4 1 r .08067 .50099 m .08175 .5042 L .08284 .50742 L .08393 .51067 L .08502 .51397 L .08611 .51732 L .08721 .52075 L .08832 .52426 L .08943 .52788 L .09055 .53161 L .09168 .53547 L .09282 .53947 L .09398 .54363 L .09514 .54797 L .09631 .55249 L .0975 .55721 L .09871 .56216 L .09993 .56733 L .10117 .57269 L .10243 .57809 L .10373 .58336 L .10507 .58835 L .10645 .59289 L .10789 .59682 L .10939 .59998 L .11095 .6022 L .1126 .60332 L .11432 .60318 L .11614 .60161 L .11805 .59846 L .12007 .59356 L .1222 .58675 L .12444 .57787 L .12682 .56674 L .12932 .55324 L .13195 .53742 L .13467 .51963 L .13748 .5002 L .14034 .47945 L .14323 .45773 L .14613 .43536 L .14901 .41268 L .15185 .39003 L .15463 .36774 L .15732 .34615 L .1599 .32558 L .16234 .30637 L .16463 .28886 L .16674 .27338 L .16864 .26026 L Mistroke .17031 .24984 L .17174 .24241 L .17296 .23786 L .17403 .2359 L .17502 .23622 L .17598 .23851 L .17699 .24248 L .17809 .24781 L .17936 .25421 L .18086 .26137 L .18266 .26898 L .1848 .27675 L .18736 .28437 L .1904 .29153 L .19398 .29794 L .19817 .30328 L .20302 .30726 L .2086 .30956 L .21495 .30999 L .222 .30869 L .22967 .30592 L .23787 .30194 L .2465 .29699 L .25549 .29133 L .26474 .2852 L .27417 .27887 L .28369 .27257 L .29321 .26658 L .30264 .26113 L .31189 .25647 L .32088 .25287 L .32952 .25057 L .33772 .24983 L .34539 .25089 L .35245 .25401 L .35885 .25927 L .36463 .26639 L .36985 .27509 L .37456 .28506 L .37882 .29599 L .38268 .30759 L .38619 .31956 L .3894 .33159 L .39237 .34339 L .39516 .35465 L .3978 .36507 L .40037 .37436 L .40291 .3822 L .40547 .3883 L .40812 .39236 L Mistroke .41089 .39407 L .41385 .39316 L .417 .38963 L .42034 .38384 L .42383 .37613 L .42745 .36684 L .43118 .35634 L .435 .34497 L .43889 .33308 L .44282 .32102 L .44677 .30914 L .45072 .29779 L .45464 .28732 L .45852 .27809 L .46232 .27044 L .46604 .26472 L .46964 .26128 L .4731 .26047 L .47641 .2626 L .47959 .26755 L .48271 .27498 L .48581 .28458 L .48895 .296 L .49219 .30891 L .49557 .32298 L .49916 .33788 L .50301 .35327 L .50718 .36883 L .51172 .38422 L .51667 .39911 L .52211 .41317 L .52808 .42606 L .53465 .43746 L .54185 .44702 L .54975 .45443 L .55838 .45943 L .56766 .46216 L .57747 .46283 L .58772 .46165 L .5983 .45884 L .60911 .45462 L .62003 .4492 L .63096 .4428 L .64179 .43563 L .65243 .42791 L .66275 .41986 L .67267 .41169 L .68206 .40362 L .69083 .39586 L .69887 .38863 L Mistroke .70606 .38215 L .71232 .37662 L .71759 .37218 L .72194 .36873 L .72546 .36619 L .72823 .36444 L .73034 .3634 L .73187 .36296 L .73291 .36303 L .73355 .3635 L .73386 .36428 L .73394 .36527 L .73386 .36636 L .73372 .36747 L .7336 .36848 L .73358 .36931 L .73375 .36985 L .73419 .37 L .73499 .36967 L .73619 .36884 L .73778 .36755 L .73976 .36585 L .74213 .3638 L .74488 .36144 L .74801 .35884 L .75151 .35603 L .75539 .35308 L .75962 .35003 L .76422 .34694 L .76918 .34385 L .77448 .34082 L .78014 .3379 L .78614 .33514 L .79248 .3326 L .79916 .33032 L .80616 .32836 L .81348 .3267 L .8211 .32534 L .82899 .32426 L .83714 .32344 L .84553 .32284 L .85413 .32247 L .86293 .32229 L .87191 .32229 L .88105 .32245 L .89032 .32274 L .89971 .32315 L .90921 .32366 L .91878 .32425 L .9284 .3249 L Mistroke .93807 .32559 L .94776 .3263 L Mfstroke % End of Graphics MathPictureEnd \ \>"], "Graphics", CellLabel->"From In[87]:=", ImageSize->{288, 177.938}, ImageMargins->{{35, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHg=_N`03001_Nfmk025_N`03 001_Nfmk009_N`004Fmk00<006mkKg/0hVmk0P008fmk00<006mkKg/00Vmk000AKg/00`00Kg]_N`3Q Kg/00`00Kg/0000SKg/00`00Kg]_N`02Kg/0015_N`03001_Nfmk0>5_N`03001_N`0002=_N`03001_ Nfmk009_N`004Fmk00<006mkKg/0hFmk00<006mk00000Vmk00<006mkKg/07Vmk00<006mkKg/00Vmk 000AKg/00`00Kg]_N`3QKg/00`00Kg/00002Kg/00`00Kg]_N`0NKg/00`00Kg]_N`02Kg/0015_N`03 001_Nfmk0>1_N`05001_Nfmk001_N`02000PKg/00`00Kg]_N`02Kg/0015_N`800>5_N`05001_Nfmk 001_N`02000OKg/20004Kg/0015_N`03001_Nfmk0>1_N`05001_Nfmk001_N`02000PKg/00`00Kg]_ N`02Kg/0015_N`03001_Nfmk0>1_N`05001_Nfmk001_N`02000PKg/00`00Kg]_N`02Kg/0015_N`03 001_Nfmk0>1_N`04001_NfmkKg/200000fmk001_N`0NKg/00`00Kg]_N`02Kg/0015_N`03001_Nfmk 0>1_N`04001_NfmkKg/200000fmk001_N`0NKg/00`00Kg]_N`02Kg/0015_N`03001_Nfmk0>1_N`04 001_NfmkKg/200000fmk001_N`0NKg/00`00Kg]_N`02Kg/0015_N`03001_Nfmk0>1_N`04001_Nfmk Kg/200000fmk001_N`0NKg/00`00Kg]_N`02Kg/0015_N`03001_Nfmk0>1_N`05001_NfmkKg/00002 Kg/00`00Kg]_N`0MKg/00`00Kg]_N`02Kg/0015_N`800>1_N`03001_Nfmk009_N`04001_Nfmk000N Kg/20004Kg/0015_N`03001_Nfmk0=m_N`03001_Nfmk009_N`04001_Nfmk000OKg/00`00Kg]_N`02 Kg/0015_N`03001_Nfmk0=m_N`03001_Nfmk009_N`04001_Nfmk000OKg/00`00Kg]_N`02Kg/0015_ N`03001_Nfmk0=m_N`03001_Nfmk009_N`04001_Nfmk000OKg/00`00Kg]_N`02Kg/0015_N`03001_ Nfmk0=m_N`03001_Nfmk009_N`04001_Nfmk000OKg/00`00Kg]_N`02Kg/0015_N`03001_Nfmk0=m_ N`03001_Nfmk009_N`05001_NfmkKg/0000NKg/00`00Kg]_N`02Kg/0015_N`03001_Nfmk0=m_N`03 001_Nfmk00I_N`03001_Nfmk01a_N`03001_Nfmk009_N`004Fmk0P00gfmk00<006mkKg/01fmk0P00 76mk0P0016mk000AKg/00`00Kg]_N`3NKg/00`00Kg]_N`07Kg/2000MKg/00`00Kg]_N`02Kg/0015_ N`03001_Nfmk0=i_N`03001_Nfmk00M_N`8001e_N`03001_Nfmk009_N`004Fmk00<006mkKg/0MFmk 00<006mkKg/0IFmk00<006mkKg/02Fmk00<006mkKg/06fmk00<006mkKg/00Vmk000AKg/00`00Kg]_ N`1eKg/00`00Kg]_N`1UKg/00`00Kg]_N`09Kg/00`00Kg]_N`0KKg/00`00Kg]_N`02Kg/000Q_N`<0 00I_N`03001_Nfmk07E_N`03001_Nfmk06E_N`03001_Nfmk00Y_N`03001_Nfmk01Y_N`03001_Nfmk 009_N`002Fmk00<006mkKg/01Fmk00<006mkKg/0MFmk00<006mkKg/0IFmk00<006mkKg/02Vmk00<0 06mkKg/06Vmk00<006mkKg/00Vmk00000fmk000000030003Kg/00`00Kg]_N`05Kg/00`00Kg]_N`1d Kg/2001WKg/00`00Kg]_N`0:Kg/00`00Kg]_N`0JKg/00`00Kg]_N`02Kg/000U_N`03001_Nfmk00E_ N`8007E_N`8006M_N`03001_Nfmk00Y_N`03001_Nfmk01Q_N`<000A_N`002Fmk00<006mkKg/01Fmk 00<006mkKg/0M6mk0P00Ifmk00<006mkKg/02Vmk00@006mk000001U_N`03001_Nfmk009_N`0026mk 0P001fmk00<006mkKg/0M6mk0P00Ifmk00<006mkKg/02fmk00<006mk00006Fmk00<006mkKg/00Vmk 0009Kg/00`00Kg]_N`05Kg/00`00Kg]_N`1dKg/00`00Kg/0001VKg/00`00Kg]_N`0=Kg/00`00Kg]_ N`0GKg/00`00Kg]_N`02Kg/0015_N`03001_Nfmk07A_N`03001_N`0006I_N`03001_Nfmk00e_N`03 001_Nfmk01M_N`03001_Nfmk009_N`004Fmk00<006mkKg/0M6mk00<006mk0000IVmk00<006mkKg/0 3Fmk00<006mkKg/05fmk00<006mkKg/00Vmk000AKg/00`00Kg]_N`1cKg/01000Kg]_N`00IFmk00<0 06mkKg/03Vmk00<006mkKg/00Vmk00<006mkKg/04Vmk00<006mkKg/00Vmk000AKg/00`00Kg]_N`1` Kg/01000Kg/000000Vmk00<006mkKg/0Hfmk00<006mkKg/03Vmk00<006mkKg/00Vmk00<006mkKg/0 4Vmk00<006mkKg/00Vmk000AKg/2001aKg/01000Kg/000000Vmk00<006mkKg/0Hfmk00<006mkKg/0 3Vmk00<006mk00000Vmk00<006mkKg/04Fmk0P0016mk000AKg/00`00Kg]_N`1`Kg/00`00Kg/00003 Kg/00`00Kg]_N`05Kg/00`00Kg]_N`1KKg/00`00Kg]_N`0>Kg/00`00Kg/00002Kg/00`00Kg]_N`0B Kg/00`00Kg]_N`02Kg/0015_N`03001_Nfmk06m_N`800003Kg/006mk009_N`03001_Nfmk00E_N`03 001_Nfmk05]_N`03001_Nfmk00i_N`03001_N`00009_N`03001_Nfmk019_N`03001_Nfmk009_N`00 4Fmk00<006mkKg/0Gfmk00<006mkKg/03Fmk0P0000=_N`00Kg/00Vmk00<006mkKg/01Fmk00<006mk Kg/0Ffmk00<006mkKg/03Vmk00@006mk001_N`8001A_N`03001_Nfmk009_N`004Fmk00<006mkKg/0 Gfmk00<006mkKg/03Fmk0P0000=_N`00Kg/00fmk00<006mkKg/016mk00<006mkKg/0Ffmk00<006mk Kg/03fmk0P0000=_N`00000056mk00<006mkKg/00Vmk000AKg/00`00Kg]_N`0]Kg/00`00Kg]_N`0^ Kg/2000?Kg/200000fmk001_N`03Kg/00`00Kg]_N`03Kg/2001MKg/00`00Kg]_N`0?Kg/200000fmk 0000000DKg/00`00Kg]_N`02Kg/0015_N`03001_Nfmk02e_N`03001_Nfmk02i_N`8000=_N`8000U_ N`05001_N`00Kg/00004Kg/00`00Kg]_N`03Kg/2001MKg/00`00Kg]_N`0?Kg/2000016mk001_N`00 4fmk00<006mkKg/00Vmk000AKg/2000^Kg/00`00Kg]_N`0]Kg/00`00Kg/00003Kg/20009Kg/01@00 Kg/006mk00001Fmk00<006mkKg/00Vmk00<006mk0000G6mk00<006mkKg/03fmk00H006mk000006mk 000=Kg/00`00Kg]_N`02Kg/20004Kg/0015_N`03001_Nfmk02e_N`03001_Nfmk00U_N`03001_Nfmk 01a_N`8000=_N`04001_Nfmk0002Kg/20009Kg/01@00Kg/006mk00001Fmk00<006mkKg/00Vmk00<0 06mk0000FVmk0P004Vmk00H006mk000006mk000=Kg/00`00Kg]_N`03Kg/00`00Kg]_N`02Kg/0015_ N`03001_Nfmk02a_N`8000]_N`8001e_N`8000=_N`08001_Nfmk001_N`00Kg/000U_N`04001_Nfmk 0007Kg/01`00Kg]_Nfmk001_N`0005Y_N`80019_N`06001_N`00001_N`003Fmk00<006mkKg/00fmk 00<006mkKg/00Vmk000AKg/00`00Kg]_N`0/Kg/2000;Kg/00`00Kg/0000LKg/20003Kg/02000Kg]_ N`00Kg/006mk0009Kg/01000Kg]_N`001fmk00@006mkKg/0009_N`03001_Nfmk05Q_N`03001_Nfmk 015_N`03001_N`00009_N`03001_Nfmk00]_N`03001_Nfmk00=_N`03001_Nfmk009_N`004Fmk00<0 06mkKg/0;6mk0P002fmk00<006mk000076mk0P000fmk00P006mkKg/006mk001_N`002Fmk00@006mk Kg/000M_N`04001_Nfmk0002Kg/00`00Kg]_N`1HKg/00`00Kg]_N`0AKg/00`00Kg/00002Kg/2000< Kg/00`00Kg]_N`03Kg/00`00Kg]_N`02Kg/0015_N`03001_Nfmk02a_N`8000]_N`03001_N`0001]_ N`03001_N`00009_N`07001_NfmkKg/006mk00000Vmk00<006mkKg/01Vmk00@006mkKg/000M_N`04 001_Nfmk0002Kg/00`00Kg]_N`1HKg/00`00Kg]_N`0AKg/00`00Kg/00002Kg/2000Kg/00`00Kg]_N`07Kg/00`6OKg]_N`04Kg/00`00Kg]_N`04Kg/0 0`00Kg]_N`07Kg/00`6OKg]_N`03Kg/00`6OKg]_N`03Kg/00`00Kg]_N`0gKg/00`00Kg]_N`05Kg/0 1000Kg]_N`000Vmk0P001fmk00<006mk00008Fmk00@006mkKg]_N`800003Kg/006mk00I_N`03001_ Nfmk009_N`0026mk00@006mkKg/000E_N`03001_Nfmk02E_N`05001_NfmkKg/1W`04Kg/01@6OKg]_ Nfmk000016mk00<006mkKg/00fmk00@006mkKg/0009_N`030001W`0000m_N`03001_Nfmk00Q_N`03 0Im_Nfmk00=_N`03001_Nfmk00A_N`03001_Nfmk00M_N`030Im_Nfmk00A_N`030Im_Nfmk009_N`03 001_Nfmk03M_N`03001_Nfmk00E_N`03001_Nfmk00A_N`03001_Nfmk00E_N`03001_N`00025_N`08 001_Nfmk001_N`00Kg/000M_N`03001_Nfmk009_N`0026mk00@006mkKg/000E_N`03001_Nfmk02E_ N`05001_NfmkKg/1W`05Kg/0106OKg]_N`000fmk00<006mkKg/016mk00<006mk00000Vmk0P6O00=_ N`00Kg/03Vmk00<006mkKg/02Fmk00<1WfmkKg/00fmk00<006mkKg/00fmk00<006mkKg/01Vmk00<1 WfmkKg/01Fmk00<1WfmkKg/00Vmk00<006mkKg/0=fmk00<006mkKg/01Fmk00<006mkKg/016mk00<0 06mkKg/01Fmk00<006mk00008Fmk00P006mkKg/006mk001_N`001fmk00<006mkKg/00Vmk0009Kg/2 0006Kg/00`00Kg]_N`0UKg/01@00Kg]_Nfmk0Il01Fmk00D1WfmkKg]_N`00009_N`03001_Nfmk00A_ N`03001_N`000081W`03Kg/006mk008000i_N`03001_Nfmk00U_N`030Im_Nfmk00=_N`03001_Nfmk 00=_N`03001_Nfmk00I_N`030Im_Nfmk00E_N`030Im_Nfmk009_N`03001_Nfmk03M_N`03001_Nfmk 00=_N`03001_N`0000I_N`03001_Nfmk00A_N`03001_Nfmk029_N`08001_Nfmk001_N`00Kg/000M_ N`03001_Nfmk009_N`004Fmk00<006mkKg/09Fmk00D006mkKg]_N`6O00I_N`040Im_Nfmk0002Kg/0 0`00Kg]_N`05Kg/00`000Im_N`02Kg/200000fmk001_N`06mk0P001Fmk00<006mk 00001fmk00<006mkKg/00fmk00<006mkKg/08Vmk00<006mk00000Vmk00<006mk00001fmk00<006mk Kg/00Vmk000AKg/00`00Kg]_N`0UKg/01@00Kg]_Nfmk0Il01Vmk00@1WfmkKg/0009_N`03001_Nfmk 00A_N`81W`04001_NfmkKg/200000fmk001_N`06mk00<006mk000016mk00<006mk 00001fmk00<006mkKg/00fmk00<006mkKg/08fmk0P000Vmk00<006mk00001fmk00<006mkKg/00Vmk 000AKg/00`00Kg]_N`0UKg/01000Kg]_N`6O26mk00<1Wfmk00000W`000<006mkKg/00Vmk0P6O00=_ N`0000000fmk0P0000=_N`00Kg/036mk0P0036mk00<1WfmkKg/00Vmk0P000fmk00<006mkKg/01Vmk 00<1WfmkKg/01fmk00D1WfmkKg]_N`0003Q_N`03001_N`0000A_N`03001_N`0000M_N`03001_Nfmk 00=_N`03001_Nfmk02=_N`8000=_N`8000M_N`03001_Nfmk009_N`004Fmk0P009Vmk00@006mkKg/1 W`U_N`040Iml00000002O0000fmk0Il1W`03Kg/20003Kg/200000fmk001_N`0Kg/0106OKg]_N`00 0fmk00<006mkKg/016mk00<1WfmkKg/02fmk00D1WfmkKg/006mk0080029_N`03001_N`0000U_N`03 001_Nfmk009_N`05001_NfmkKg/00002Kg/00`00Kg]_N`0;Kg/2O0000`000Il1W`040IlNKg/30Il3 O00:Kg/00`00Kg]_N`02Kg/0015_N`8000M_N`03001_Nfmk01]_N`04001_Nfmk0IlNKg/00`00Kg]_ N`0;Kg/00`00Kg]_N`0>Kg/01@6OKg]_Nfmk00000Vmk00<006mkKg/00fmk00<1WfmkKg/036mk00<1 WfmkKg/00P0000=_N`00Kg/08Fmk00<006mk00002Fmk00<006mkKg/00Vmk00D006mkKg]_N`00009_ N`03001_Nfmk00Y_N`@1W`=l025_N`Il00U_N`8000A_N`004Fmk00<006mkKg/01Vmk00<006mkKg/0 6fmk00@006mkKg/1Wai_N`03001_Nfmk00Y_N`03001_Nfmk00m_N`050Im_NfmkKg/00002Kg/00`00 Kg]_N`03Kg/00`6OKg]_N`0=Kg/01P6OKg]_N`00Kg/0029_N`03001_N`0000U_N`8000=_N`05001_ NfmkKg/00002Kg/00`00Kg]_N`09Kg/00`6OKg]l0005O00RKg/4O00;Kg/00`00Kg]_N`02Kg/0015_ N`03001_Nfmk00I_N`03001_Nfmk01]_N`04001_Nfmk0IlNKg/00`00Kg]_N`0:Kg/00`00Kg]_N`0? Kg/01@6OKg]_Nfmk00000Vmk00<006mkKg/00fmk00<1WfmkKg/03Fmk00H1WfmkKg/006mk000RKg/0 0`00Kg/00009Kg/20003Kg/01@00Kg]_Nfmk00000Vmk00<006mkKg/01fmk0P6O0fmk17`0Kg/00`6OKg]_N`02Kg/00`00Kg]_N`0PKg/01000Kg]_N`001fmk00<006mk00000fmk00D0 06mkKg]_N`00009_N`03001_Nfmk00=_N`030Im_Nfmk03a_N`03001_Nfmk009_N`004Fmk00<006mk Kg/01fmk00<006mkKg/06Vmk00@006mkKg/1Wai_N`03001_Nfmk00Y_N`03001_Nfmk011_N`040Im_ Nfmk0002Kg/00`00Kg]_N`02Kg/00`6OKg]_N`0?Kg/00`6OKg]_N`02Kg/00`00Kg]_N`0OKg/01000 Kg]_N`001fmk00<006mk00000fmk00D006mkKg]_N`00009_N`03001_Nfmk009_N`030Im_Nfmk03e_ N`03001_Nfmk009_N`004Fmk00<006mkKg/01fmk00<006mkKg/06Vmk00@006mkKg/1Wam_N`03001_ Nfmk00U_N`03001_Nfmk011_N`040Im_Nfmk0002Kg/01@00Kg]_Nfmk0Il04Vmk00<1WfmkKg/00Vmk 00<006mkKg/07fmk00@006mkKg/000M_N`06001_Nfmk001_N`0016mk00<006mk000016mk00<1Wfmk Kg/0?Vmk00<006mkKg/00Vmk000AKg/00`00Kg]_N`07Kg/00`00Kg]_N`0JKg/01000Kg]_N`6O7fmk 00<006mkKg/02Fmk00<006mkKg/046mk00@1WfmkKg/0009_N`05001_NfmkKg/1W`0BKg/00`6OKg]_ N`02Kg/00`00Kg]_N`0NKg/00`00Kg]_N`02Kg/00`00Kg]_N`04Kg/01P00Kg]_N`00Kg/000E_N`80 009_N`81Wd5_N`03001_Nfmk009_N`004Fmk0P0026mk00<006mkKg/06Vmk00<006mk0Il086mk00<0 06mk00000Vmk00<006mkKg/016mk00<006mkKg/04Fmk00<1Wfmk00000Vmk00D006mkKg]_N`6O019_ N`030Im_Nfmk009_N`03001_Nfmk01i_N`03001_Nfmk009_N`03001_Nfmk00A_N`06001_Nfmk001_ N`0016mk00El0000001l006O049_N`8000A_N`004Fmk00<006mkKg/01fmk0P006fmk00<006mk0Il0 86mk00<006mk00000Vmk00<006mkKg/016mk00<006mkKg/04Fmk00D1Wfmk001_N`0000=_N`030Im_ Nfmk019_N`050Im_NfmkKg/0000PKg/00`00Kg]_N`02Kg/00`00Kg]_N`04Kg/01P00Kg]_N`00Kg/0 00=_N`<1W`03O001Wg`004=_N`03001_Nfmk009_N`004Fmk00<006mkKg/01fmk0P006fmk00<006mk 0Il086mk00<006mk00000Vmk00<006mkKg/016mk00<006mkKg/04Fmk00<1WfmkKg/00P000fmk00<1 WfmkKg/04Vmk00D1WfmkKg]_N`00021_N`03001_Nfmk009_N`03001_Nfmk00A_N`04001_NfmkKg/2 0002O0000`6OO01l00020Il2O013Kg/00`00Kg]_N`02Kg/0015_N`03001_Nfmk00M_N`<001Y_N`03 001_N`6O021_N`<0009_N`8000E_N`03001_Nfmk015_N`030Im_Nfmk008000=_N`030Im_Nfmk019_ N`050Im_NfmkKg/0000PKg/00`00Kg]_N`02Kg/00`00Kg]_N`04Kg/01000Kg]_Ng`00P0000=l006O O0001W`0@fmk00<006mkKg/00Vmk000AKg/00`00Kg]_N`07Kg/00`00Kg/0000JKg/00`00Kg/1W`0P Kg/30002Kg/20005Kg/00`00Kg]_N`0AKg/00`6OKg]_N`020002Kg/00`6OKg]_N`0CKg/01@6OKg]_ Nfmk00007fmk00<006mkKg/00fmk00<006mkKg/016mk00<006mkKg/00W`000<0006OO0001g`0@fmk 00<006mkKg/00Vmk000AKg/00`00Kg]_N`07Kg/01000Kg]_N`006Fmk00<006mk0Il086mk0`000Vmk 0P0000=_N`00Kg/00Vmk00<006mkKg/04Fmk00<1WfmkKg/00P000Vmk00<1WfmkKg/056mk00@1Wfmk Kg/001m_N`03001_Nfmk00=_N`03001_Nfmk00A_N`03001_Nfmk009l00050Iml07`0O01_N`04O014 Kg/00`00Kg]_N`02Kg/0015_N`03001_Nfmk00Y_N`03001_Nfmk01I_N`03001_N`6O025_N`03001_ N`00009_N`800003Kg/006mk009_N`03001_Nfmk019_N`060Iml0000001l006O5fmk00@1WfmkKg/0 009_N`03001_Nfmk00A_N`03001_Nfmk01=_N`03001_Nfmk00=_N`03001_Nfmk00A_N`03001_Nfmk 0081W`Al04U_N`03001_Nfmk009_N`004Fmk00<006mkKg/02Vmk00<006mkKg/05Vmk00<006mk0Il0 8Fmk00D006mkKg/006mk00800005Kg/006mkKg/0000EKg/01W`00Iml0000O001WaM_N`040Im_Nfmk 0002Kg/00`00Kg]_N`04Kg/2000DKg/00`00Kg]_N`03Kg/00`00Kg]_N`03Kg/01@00Kg]_N`6OKg/0 17`0BVmk00<006mkKg/00Vmk000AKg/2000I_N`03001_Nfmk009_N`00 4Fmk00<006mkKg/05Vmk00@1W`00Kg/000=_N`04001_N`000IoVKg/00`00Kg]_N`02Kg/0015_N`03 001_Nfmk01I_N`040Il006mk0004Kg/200000`6OKg]_N`3TKg/00`00Kg]_N`02Kg/0015_N`03001_ Nfmk01I_N`040Il006mk0004Kg/00`000Im_N`3VKg/00`00Kg]_N`02Kg/000M_N`D000E_N`03001_ Nfmk01I_N`030Il0000000A_N`04O000006OO03VKg/00`00Kg]_N`02Kg/000Q_N`04001_Nfmk0005 Kg/00`00Kg]_N`0FKg/017`00Il007`00Vmk0g`000<1Wg`0O000iFmk00<006mkKg/00Vmk0009Kg/0 0`00Kg]_N`05Kg/00`00Kg]_N`0EKg/2O0001@6O001l07`0Kg/00g`000<1Wg`0O000iFmk00<006mk Kg/00Vmk000:Kg/00`00Kg]_N`04Kg/2000FKg/2O0001@6OO01l07`0Kg/00g`000<1Wg`0O000hfmk 0`0016mk0007Kg/01@00Kg]_Nfmk00001Fmk00<006mkKg/05Fmk0W`000D1Wg`0O01l06mk00=l0003 0Iml07`00>E_N`03001_Nfmk009_N`001fmk00D006mkKg]_N`0000E_N`03001_Nfmk01E_N`=l0003 0Iml07`0009_N`9l00030Iml06mk0>E_N`03001_Nfmk009_N`0026mk0`001Vmk00<006mkKg/05Vmk 0W`000<1Wg`0Kg/00fmk00<1WfmkKg/0iVmk00<006mkKg/00Vmk000AKg/00`00Kg]_N`0HKg/00`6O Kg]_N`03Kg/00`6OKg]_N`3VKg/00`00Kg]_N`02Kg/0015_N`03001_Nfmk01Q_N`030Im_Nfmk00=_ N`030Im_Nfmk0>I_N`03001_Nfmk009_N`004Fmk00<006mkKg/066mk00<1WfmkKg/00Vmk00<1Wfmk Kg/0ifmk00<006mkKg/00Vmk000AKg/2000JKg/01@6OKg]_Nfmk0Il0j6mk0P0016mk000AKg/00`00 Kg]_N`0IKg/01@6OKg]_Nfmk0Il0jFmk00<006mkKg/00Vmk000AKg/00`00Kg]_N`0JKg/00`6OKg/1 W`3ZKg/00`00Kg]_N`02Kg/0015_N`03001_Nfmk01Y_N`030Im_N`6O0>Y_N`03001_Nfmk009_N`00 4Fmk00<006mkKg/06fmk00<1WfmkKg/0jFmk00<006mkKg/00Vmk000AKg/00`00Kg]_N`3oKg/8Kg/0 0`00Kg]_N`02Kg/0015_N`03001_Nfmk00=_N`03001_Nfmk00Y_N`03001_Nfmk00U_N`03001_Nfmk 00Y_N`03001_Nfmk00Y_N`03001_Nfmk00U_N`03001_Nfmk00Y_N`03001_Nfmk00U_N`03001_Nfmk 00Y_N`03001_Nfmk00Y_N`03001_Nfmk00U_N`03001_Nfmk00Y_N`03001_Nfmk00Y_N`03001_Nfmk 00U_N`03001_Nfmk00Y_N`03001_Nfmk00U_N`03001_Nfmk00Y_N`03001_Nfmk00Y_N`03001_Nfmk 00U_N`03001_Nfmk00Y_N`03001_Nfmk00Y_N`03001_Nfmk00A_N`03001_Nfmk009_N`004Fmko`00 300016mk0000\ \>"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {-18.7307, -2.23972, \ 0.792407, 0.0261527}}], Cell[BoxData[ RowBox[{"{", RowBox[{\(0.43333333333339397`\ Second\), ",", TagBox[\(\[SkeletonIndicator] Graphics \[SkeletonIndicator]\), False, Editable->False]}], "}"}]], "Output", CellLabel->"Out[87]="] }, Open ]] }, Open ]] }, Open ]] }, Open ]] }, FrontEndVersion->"4.0 for Macintosh", ScreenRectangle->{{0, 1152}, {0, 850}}, WindowSize->{944, 749}, WindowMargins->{{Automatic, 68}, {Automatic, 14}}, StyleDefinitions -> "ArticleModern.nb", MacintoshSystemPageSetup->"\<\ 01h0001804P000000aP2@_ogooP3@09:8085N`?P0000005X0FP000003gP;BP5/ 038;BTLH05000@4100000BL?00400@000000000000000000000000000040I000 00000000000@X@804@000000001C>30`\>" ] (*********************************************************************** 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[1739, 51, 40, 0, 34, "Subtitle"], Cell[1782, 53, 309, 9, 76, "Author"], Cell[CellGroupData[{ Cell[2116, 66, 40, 0, 63, "Section"], Cell[2159, 68, 73, 2, 42, "Input"], Cell[CellGroupData[{ Cell[2257, 74, 568, 8, 138, "Input"], Cell[2828, 84, 546, 8, 101, "Output"] }, Open ]], Cell[3389, 95, 110, 2, 42, "Input"], Cell[3502, 99, 1260, 23, 186, "Input"], Cell[4765, 124, 1383, 26, 202, "Input"], Cell[CellGroupData[{ Cell[6173, 154, 36, 0, 42, "Subsubsection"], Cell[CellGroupData[{ Cell[6234, 158, 140, 3, 27, "Input"], Cell[6377, 163, 115, 2, 70, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[6529, 170, 187, 3, 43, "Input"], Cell[6719, 175, 160, 3, 70, "Output"] }, Open ]], Cell[6894, 181, 162, 3, 27, "Input"], Cell[7059, 186, 216, 4, 43, "Input"] }, Closed]], Cell[CellGroupData[{ Cell[7312, 195, 38, 0, 26, "Subsubsection"], Cell[7353, 197, 133, 3, 42, "Input"], Cell[7489, 202, 382, 8, 74, "Input"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[7920, 216, 40, 0, 63, "Section"], Cell[7963, 218, 442, 11, 40, "Text"], Cell[8408, 231, 95, 2, 42, "Input"], Cell[8506, 235, 385, 7, 74, "Input"], Cell[8894, 244, 449, 10, 90, "Input"], Cell[CellGroupData[{ Cell[9368, 258, 38, 0, 42, "Subsubsection"], Cell[CellGroupData[{ Cell[9431, 262, 103, 2, 42, "Input"], Cell[9537, 266, 290, 5, 41, "Output"] }, Open ]], Cell[9842, 274, 174, 4, 24, "Text"], Cell[10019, 280, 975, 20, 90, "Input"], Cell[10997, 302, 58, 0, 24, "Text"], Cell[11058, 304, 1321, 25, 170, "Input"], Cell[12382, 331, 1383, 26, 186, "Input"], Cell[CellGroupData[{ Cell[13790, 361, 113, 2, 42, "Input"], Cell[13906, 365, 39183, 1131, 194, 10874, 776, "GraphicsData", "PostScript", \ "Graphics"], Cell[53092, 1498, 248, 6, 41, "Output"] }, Open ]] }, Open ]] }, Open ]] }, Open ]] } ] *) (*********************************************************************** End of Mathematica Notebook file. ***********************************************************************)