(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 9.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 157, 7] NotebookDataLength[ 42999, 1045] NotebookOptionsPosition[ 41214, 977] NotebookOutlinePosition[ 41557, 992] CellTagsIndexPosition[ 41514, 989] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"f", "[", RowBox[{"n_", ",", "p_"}], "]"}], ":=", RowBox[{ RowBox[{"n", RowBox[{"(", RowBox[{"1", "-", FractionBox["n", "6"]}], ")"}]}], "-", FractionBox[ RowBox[{"n", "*", "p"}], RowBox[{"(", RowBox[{"4", "+", RowBox[{"1.2", "n"}], "+", "p"}], ")"}]]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"g", "[", RowBox[{"n_", ",", "p_"}], "]"}], ":=", RowBox[{ FractionBox[ RowBox[{"n", "*", "p"}], RowBox[{"(", RowBox[{"4", "+", RowBox[{"1.2", "n"}], "+", "p"}], ")"}]], "-", "p"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Solve", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"n", RowBox[{"(", RowBox[{"1", "-", FractionBox["n", "6"]}], ")"}]}], "-", FractionBox[ RowBox[{"n", "*", "p"}], RowBox[{"(", RowBox[{"4", "+", RowBox[{"1.2", "n"}], "+", "p"}], ")"}]]}], "\[Equal]", "0"}], "&&", RowBox[{ RowBox[{ FractionBox[ RowBox[{"n", "*", "p"}], RowBox[{"(", RowBox[{"4", "+", RowBox[{"1.2", "n"}], "+", "p"}], ")"}]], "-", "p"}], "\[Equal]", "0"}]}], ",", RowBox[{"{", RowBox[{"n", ",", "p"}], "}"}], ",", "Reals"}], "]"}], " "}]}], "Input", CellChangeTimes->{{3.6621414745188074`*^9, 3.6621415351147385`*^9}, { 3.6621632908146877`*^9, 3.662163503752867*^9}, {3.662163577785102*^9, 3.6621636655481215`*^9}, {3.6621746376582565`*^9, 3.662174639002925*^9}, { 3.6621766448111167`*^9, 3.6621766662498894`*^9}, 3.6622597909397426`*^9, { 3.6623078340201416`*^9, 3.662307986570867*^9}, {3.6623466129361596`*^9, 3.6623466536210513`*^9}, {3.66234674857813*^9, 3.6623468019241004`*^9}, { 3.6623473061761103`*^9, 3.6623473316248884`*^9}, {3.662347369712449*^9, 3.6623474500133495`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"n", "\[Rule]", RowBox[{"-", "2.479473661428269`"}]}], ",", RowBox[{"p", "\[Rule]", RowBox[{"-", "3.5041052677143463`"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"n", "\[Rule]", "0"}], ",", RowBox[{"p", "\[Rule]", "0"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"n", "\[Rule]", "6.`"}], ",", RowBox[{"p", "\[Rule]", "0"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"n", "\[Rule]", "9.67947366142827`"}], ",", RowBox[{"p", "\[Rule]", RowBox[{"-", "5.935894732285654`"}]}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{ 3.662307784210293*^9, 3.6623079912581353`*^9, 3.6623466557930174`*^9, 3.6623468049867287`*^9, {3.6623474388566256`*^9, 3.6623474514352713`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"jacobianMatrix", "[", RowBox[{"n_", ",", "p_"}], "]"}], "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"f", "[", RowBox[{"n", ",", "p"}], "]"}], ",", "n"}], "]"}], ",", RowBox[{"D", "[", RowBox[{ RowBox[{"f", "[", RowBox[{"n", ",", "p"}], "]"}], ",", "p"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"g", "[", RowBox[{"n", ",", "p"}], "]"}], ",", "n"}], "]"}], ",", RowBox[{"D", "[", RowBox[{ RowBox[{"g", "[", RowBox[{"n", ",", "p"}], "]"}], ",", "p"}], "]"}]}], "}"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"MatrixForm", "[", RowBox[{"jacobianMatrix", "[", RowBox[{"n", ",", "p"}], "]"}], "]"}]}], "Input", CellChangeTimes->{{3.662163732707963*^9, 3.6621638021829367`*^9}, { 3.662163971916645*^9, 3.66216397899605*^9}, {3.662164182409684*^9, 3.66216420244083*^9}}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{"1", "-", FractionBox["n", "3"], "+", FractionBox[ RowBox[{"1.2`", " ", "n", " ", "p"}], SuperscriptBox[ RowBox[{"(", RowBox[{"4", "+", RowBox[{"1.2`", " ", "n"}], "+", "p"}], ")"}], "2"]], "-", FractionBox["p", RowBox[{"4", "+", RowBox[{"1.2`", " ", "n"}], "+", "p"}]]}], RowBox[{ FractionBox[ RowBox[{"n", " ", "p"}], SuperscriptBox[ RowBox[{"(", RowBox[{"4", "+", RowBox[{"1.2`", " ", "n"}], "+", "p"}], ")"}], "2"]], "-", FractionBox["n", RowBox[{"4", "+", RowBox[{"1.2`", " ", "n"}], "+", "p"}]]}]}, { RowBox[{ RowBox[{"-", FractionBox[ RowBox[{"1.2`", " ", "n", " ", "p"}], SuperscriptBox[ RowBox[{"(", RowBox[{"4", "+", RowBox[{"1.2`", " ", "n"}], "+", "p"}], ")"}], "2"]]}], "+", FractionBox["p", RowBox[{"4", "+", RowBox[{"1.2`", " ", "n"}], "+", "p"}]]}], RowBox[{ RowBox[{"-", "1"}], "-", FractionBox[ RowBox[{"n", " ", "p"}], SuperscriptBox[ RowBox[{"(", RowBox[{"4", "+", RowBox[{"1.2`", " ", "n"}], "+", "p"}], ")"}], "2"]], "+", FractionBox["n", RowBox[{"4", "+", RowBox[{"1.2`", " ", "n"}], "+", "p"}]]}]} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output", CellChangeTimes->{3.6623077842592955`*^9, 3.6623080044488897`*^9, 3.6623466610588617`*^9, 3.6623468080181055`*^9, 3.6623474548885403`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ss1", "=", RowBox[{"jacobianMatrix", "[", RowBox[{"0", ",", "0"}], "]"}]}]], "Input", CellChangeTimes->{{3.6623077576697745`*^9, 3.662307758129801*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"1.`", ",", "0.`"}], "}"}], ",", RowBox[{"{", RowBox[{"0.`", ",", RowBox[{"-", "1.`"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{ 3.6623077842782965`*^9, 3.66230802421402*^9, {3.662308708231144*^9, 3.6623087112483163`*^9}, 3.662346663980857*^9, 3.6623468118307643`*^9, 3.6623474589512115`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Eigenvalues", "[", "ss1", "]"}]], "Input", CellChangeTimes->{{3.6621746301969433`*^9, 3.6621746534039087`*^9}, { 3.6623077607139487`*^9, 3.6623077631270866`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"-", "1.`"}], ",", "1.`"}], "}"}]], "Output", CellChangeTimes->{3.6623077842972975`*^9, 3.6623087688489356`*^9, 3.6623466654184256`*^9, 3.662346813080804*^9, 3.662347461232566*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ss2", "=", RowBox[{"jacobianMatrix", "[", RowBox[{"2.479473661428269`", ",", RowBox[{"-", "3.5041052677143463`"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.6623089054271755`*^9, 3.6623089058483763`*^9}, { 3.6623089521180573`*^9, 3.6623089557372637`*^9}, {3.6623467090764847`*^9, 3.6623467103265257`*^9}, 3.662346817752885*^9, {3.6623477143194027`*^9, 3.662347721134907*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0.31771751526248626`", ",", RowBox[{"-", "1.4353293929785926`"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "0.14420873573857595`"}], ",", "0.4353293929785925`"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{ 3.66230895941887*^9, {3.6623466874974623`*^9, 3.6623467171236954`*^9}, 3.662346819159191*^9, 3.662347722654809*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Eigenvalues", "[", "ss2", "]"}]], "Input", CellChangeTimes->{{3.6623089671096835`*^9, 3.6623089675308847`*^9}, { 3.66234671268601*^9, 3.662346712920395*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"0.8352664928413875`", ",", RowBox[{"-", "0.08221958460030876`"}]}], "}"}]], "Output", CellChangeTimes->{3.662347746294318*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ss3", "=", RowBox[{"jacobianMatrix", "[", RowBox[{"6", ",", "0"}], "]"}]}]], "Input", CellChangeTimes->{{3.6623477645258694`*^9, 3.662347782620373*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "1.`"}], ",", RowBox[{"-", "0.5357142857142857`"}]}], "}"}], ",", RowBox[{"{", RowBox[{"0.`", ",", RowBox[{"-", "0.4642857142857143`"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{{3.662347778948344*^9, 3.6623477839954343`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Eigenvalues", "[", "ss3", "]"}]], "Input", CellChangeTimes->{{3.6623477958240476`*^9, 3.6623477960740576`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"-", "1.`"}], ",", RowBox[{"-", "0.4642857142857143`"}]}], "}"}]], "Output", CellChangeTimes->{3.6623477969959707`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ss4", "=", RowBox[{"jacobianMatrix", "[", RowBox[{"9.67947366142827`", ",", "5.935894732285654`"}], "]"}]}]], "Input", CellChangeTimes->{{3.662347818137476*^9, 3.6623478339662604`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "2.3534748433422843`"}], ",", RowBox[{"-", "0.325430698319373`"}]}], "}"}], ",", RowBox[{"{", RowBox[{"0.12698362286619402`", ",", RowBox[{"-", "0.6745693016806271`"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.6623478366069927`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Eigenvalues", "[", "ss4", "]"}]], "Input", CellChangeTimes->{{3.6623478045900373`*^9, 3.662347804840047*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"-", "2.32848912856166`"}], ",", RowBox[{"-", "0.6995550164612516`"}]}], "}"}]], "Output", CellChangeTimes->{3.662347838638327*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"u1", "=", RowBox[{"NDSolve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"f", "[", RowBox[{ RowBox[{"n", "[", "t", "]"}], ",", RowBox[{"p", "[", "t", "]"}]}], "]"}], "\[Equal]", RowBox[{ RowBox[{"n", "'"}], "[", "t", "]"}]}], "&&", RowBox[{ RowBox[{"g", "[", RowBox[{ RowBox[{"n", "[", "t", "]"}], ",", RowBox[{"p", "[", "t", "]"}]}], "]"}], "\[Equal]", RowBox[{ RowBox[{"p", "'"}], "[", "t", "]"}]}], "&&", RowBox[{ RowBox[{"n", "[", "0", "]"}], "\[Equal]", "1"}], "&&", RowBox[{ RowBox[{"p", "[", "0", "]"}], "\[Equal]", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"n", "[", "t", "]"}], ",", RowBox[{"p", "[", "t", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "10"}], "}"}]}], "]"}]}], ";"}], "\n", RowBox[{"plot1", "=", RowBox[{"ParametricPlot", "[", RowBox[{ RowBox[{"Evaluate", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"n", "[", "t", "]"}], ",", RowBox[{"p", "[", "t", "]"}]}], "}"}], "/.", "u1"}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "10"}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", "All"}], ",", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{"Blue", ",", "Thick"}], "}"}]}], ",", RowBox[{"AxesOrigin", "\[Rule]", RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}]}], ",", RowBox[{"AxesLabel", "\[Rule]", RowBox[{"{", RowBox[{"\"\\"", ",", "\"\\""}], "}"}]}]}], "]"}]}]}], "Input", CellChangeTimes->{{3.662259841723117*^9, 3.662259862583359*^9}, { 3.6622598957722406`*^9, 3.6622599666189413`*^9}, 3.662260120609727*^9, { 3.662309007186154*^9, 3.662309007654155*^9}}], Cell[BoxData[ GraphicsBox[{{}, {}, {RGBColor[0, 0, 1], Thickness[Large], LineBox[CompressedData[" 1:eJw1z3k01PsbB/AZWy5iNPENlaWQFgmpG3p/aC+XSglJsuSmlDK0SJSSNZVE dLlKXbeUrRSibClRXUuRnWiMkWUMhiZ37vmd33POc97ndd7n+ePRdD2200OM RqMlifa//N8M4f+5pUw/xVd6CK1bGbF8ziCYvIXf+hlDKFy4Vn6qYRBdGg8v czSHUK3uUrYyYxCrLrw20Fg/BO2L3ZW/OwzCMdB/29HwITg3BewpyPuORAl2 WjI1DB23QKu8swP49Vz6tQcaw7hevNYx9tAAjtam7S7UG8bgocD2c3YDELMY 8ew3HQbr1Emt8wYD8HWuabjkMox3kxzrFT1cpNtyG9IeDKOcs9LTZAcXH1Qq ZsdajqC8stEg1KgfycHPfpP7bQQxQ8st7mv0wyY0fX/knhH0eS9+2jCzH0eN 9zGTjoxAzmhjrCubA9/QoWZa/Aj2BZfmGSZzIDB0mu7ljqBs8L1CqhwHf3KP cBf+wcNd/paOfaNs6HS9OrY3nYcE92fBO7rZeKa953t8Lg9JPQ8Gd9aycVVm SadWFQ8/Us9KB2WxMaChlxM8zoPenoLuM95s7F5hnuRkO4qmxh9kvO8bpBqf 2CQq8FH8vS4nh9uLTYYaXq9V+TD8HD9q3daLlg0+NePafFRdeSs59qEXGRYT x7zM+Dg8OoPlntsLv0Wmy2IO8RFp33+l5IzI8Xo1n8r5MAvUMvGT7YWG+94V LkFjSB823HvZuAfdjjc7VkaOITO5d4G3bg+oaLs2RvwYRrjvs51Ue9DcF1/x 6fEYbCJmCe1pPZg29Fh/p20MoZuDfCqqv0L89nfOKozj4fYwWU/Pr5BcUrbF SWICNM1vFxLudsM7o/bpbcYE3M4787vju+Hxkm/cNXcCpRfOf10V1Y0zyoFe Z1dOAIoe8wWsbngOFaUNHJzAl8mWTR0bu1Ed5futtGoCdwfcH5pyu/DmbmeF 8U0B6OmnlgrMuuD1teyi8l0BZvcU3Nhk2IWeHe8thZkC2BmVZfyh2wWjNqP8 prcCdM69VOcyqwuPvy0+/V4oAOuTp40luxMH124OvX9wEubRlrlH4joR91PL /RezKciz3B+a8DsQQne3ebx1Cm2Bys+TOB14fieeb+8whWhLd550Rwfsl6w+ /NJ/Cmz17i9TVR1Q3SD9+mP2FOy12rMsUzvQThnHiOv9gLiE/Q1rmw7oWmbP ua0mhDDcdbmLaztqG0oejy4WQuexapuaTTt25eS+2rlGCDfXDvSatsPi3lJr dQch7pn39N9XagdLfuZbwU0hmIoy1z69bQOjepdsDvMnRvO9Dn01boOu0px5 8xWn0XvAbJWHUiuk0g/HKGhMg+3MzJaRbEW4IDBIcvk0fq8eflg42gJmROb7 n1bT2ForRif1LbA2m+rUCZ/G/gAd/X9iW6DcPrBdm04jwo9TJTOVWrD8Xq/t twAa+WkvZput3gxNHPC+F0Ij/N2DD54ymlEk5pPkEUkj64YPHC6nN6OQuXI5 N5FGmu9FSdK/fkFQXpXMrAIa+TuA/YKf/gUu1VfMX07QyOYgc93Xxl+g/tSz 98ZJOtEOenG2yrYJjq9dxs2D6eSVsopB7/omnJ4qWNsXRidhn56NyJk0QZ+W 6rkxkU4mOlKDA+Y0oa3gZrdqEZ3wUk8UNrc24rPBupjF4mKkoNEx5k+vRuwQ 2yuovypGmJ/yZkxFfkZ/uNdITbY4MSrOZilwGxAdar+/skCcXKlMqepqb4BW b6ZdaZk4ieII8ovrGjBHZdGXFw3i5HC66+arhQ0ocVTrLxOIk3wDyulWZAMq s/n9lKUEcWmOUzdf1oA0g4SJ8joJohsbd/WWfz0cEph7bk5KkutXJxrNNOog udrP1lxcijAXxKzOVq6DmnuRXI+sFLFv1v9rycw6vJh/86HJPCmSJShLMZys xbCsYFcXpMhu8V8vhdXXInPquL7DJSmSqJhmtz2iFoUttwpzGTOIhO3S4hWT /8DxyH7J64ukyYGpLoutnI9gX826fOOyDFEqz22NL67Bpuf5PiZP5YjLzxUn HM5UIaFCLr+1Qp70yVtVWedWIsx4k6y8L4MU5FxTWq1UgZQPnMMf9iuSbkaa 5hNeKdzf9WXO2zaL3HbaPtdm4BVq1+RpqqoxSXyJBiN8XzG8fVpMjN4wiZTn g1q/gEJon7WpzPaeTZaOr2nNyn0OZ6vSIIfR2WT+I55BiFsevl35ULTjhBKZ U9/4ldJ5gmVxN5glzUrE+68XWt7Lc6Apa78udbkyabWzjawKz8Ln0zbWz88q k66NqYvcWh4B5W6UR74yCaJfYNToPoRDpNmBkn5lsomxTdXhVTq0rWJ9/lKh SFbQ3wq73e8jJfhUcYg5RZrmP260OpiG0MsfGT/2UWRugUnT0PE7SLPnOR87 RZGXYS0s17w/scTGco18NEXGTdUMhIHJECrX7i9MoYiKl1uFSVISrtBvbNr4 mCIJbAutNL1bCF7n3Le+gCK9A0ciKiZv4v7F38uWVFDkYIrD2di6G3ipozSR WUORgXdOtx1psSh+e8jySANFfrt/ssPG+BoG9/3x7nQzRbaWUhoZTjGIE+wc O9dBEbm9Ygd97kbhY+tYiHIPRWzMMzOLeOG4lMLwfMKmyHSJG8s34zK0KOE1 Ca7ov8XbRlYdvQST2gG/0u8U+ef+lV1Xj4XA9fqG54NDFFk4mphpUXweXMWZ MrwRiuir2e2VYwaj5qgBL3qUIhrhKiXmOwLx5HLyvj4+Rayl9aniN2eA1E+N y8YpUrll96WTm09B9vh+ruEERYQa1XLqev4I8Y1ctErkWh3Wgtka/ijNUnQ1 FTlvgyJXjvJH/UXjmnUiq9rKR/8i6Q+20S8Xd4msYmXVrd/ph885SX5+Ips9 uR2xKNEPXaZOfk9F9jUDqZH3Q5TNYZd8kV/Twx6slPLDIS35LUUiO3/2j08V sqCZocSsEJk/ksAL4bJg2eIYUS9y3N/60SeqWKApGc7giZwlljxjLJSFsKT1 1WP/3Vv2D8edY+HjGcmYSZFnXNPd/Ks/C9KchJliAopwEl8aRniw8IZaOcgQ +dX8R4271rOgMMxOnf1fH5szqWDGQtsx+51zRD4Svy3wvREL/hLHM9RFnh0e kbx7AQvHvefZLRCZEe+wZr4aC3u0ooQ6Is/ylV/BmcVCSXXxncUi0x+Zn8+X YWGz8rMN+iKX1QfoRYmx8C8EjKa6 "]]}}, Axes->True, AxesLabel->{ FormBox["\"N\"", TraditionalForm], FormBox["\"P\"", TraditionalForm]}, AxesOrigin->{0, 0}, Method->{}, PlotRange->{All, All}, PlotRangeClipping->True, PlotRangePadding->{Automatic, Automatic}]], "Output", CellChangeTimes->{3.662307784372302*^9, 3.662308978591304*^9, 3.6623090087149568`*^9, 3.6623466919195204`*^9, 3.6623468226749635`*^9, 3.662347841716572*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"u2", "=", RowBox[{"NDSolve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"f", "[", RowBox[{ RowBox[{"n", "[", "t", "]"}], ",", RowBox[{"p", "[", "t", "]"}]}], "]"}], "\[Equal]", RowBox[{ RowBox[{"n", "'"}], "[", "t", "]"}]}], "&&", RowBox[{ RowBox[{"g", "[", RowBox[{ RowBox[{"n", "[", "t", "]"}], ",", RowBox[{"p", "[", "t", "]"}]}], "]"}], "\[Equal]", RowBox[{ RowBox[{"p", "'"}], "[", "t", "]"}]}], "&&", RowBox[{ RowBox[{"n", "[", "0", "]"}], "\[Equal]", "3"}], "&&", RowBox[{ RowBox[{"p", "[", "0", "]"}], "\[Equal]", "4"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"n", "[", "t", "]"}], ",", RowBox[{"p", "[", "t", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "10"}], "}"}]}], "]"}]}], ";"}], "\n", RowBox[{"plot2", "=", RowBox[{"ParametricPlot", "[", RowBox[{ RowBox[{"Evaluate", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"n", "[", "t", "]"}], ",", RowBox[{"p", "[", "t", "]"}]}], "}"}], "/.", "u2"}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "10"}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", "All"}], ",", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{"Red", ",", "Dashed"}], "}"}]}], ",", RowBox[{"AxesOrigin", "\[Rule]", RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}]}], ",", RowBox[{"AxesLabel", "\[Rule]", RowBox[{"{", RowBox[{"\"\\"", ",", "\"\\""}], "}"}]}]}], "]"}]}]}], "Input", CellChangeTimes->{{3.6622599740255003`*^9, 3.662260002417307*^9}, 3.6622601176408434`*^9}], Cell[BoxData[ GraphicsBox[{{}, {}, {RGBColor[1, 0, 0], Dashing[{Small, Small}], LineBox[CompressedData[" 1:eJwl1Is31HkfB3CGGb/fjDEM+aljpchD6SG5bcr3E61SbIsS7cq1JKIL2TLR RSkVxTaD3JNb6WJcxprC0Nio3JZBQm6FpkVsyQ7P7zn7Pud93uf1D7xX+YW5 HqTIyclFkP3//hsM/l1V6Gu/GiulYMB95DF5UMqC+v34bB8VA5XOh8HdvSzg SLxT6xgYqPq3bRkuZ0HlxYUPOQQGZXvczBNDWGC75Dr61ASDjrARJes+FeBQ oqq9fDGYXO07d6OBCcyOxaNxARjwF/9Si+Qz4dK1otnSQAx0XAIqgnOZcC/f MIoVioE0ihMeeY4JWx/l7e0/g8EZFNJJbGGCd2pFbPtvGISbGkvCBcpwftGh T6UJA1u6he2hSgZInc4d+f0VBpRrawMO5DNg+5c3i4FtGIimZUN+txmw/9EV ixfdGJSetKq/eZIBCiM21XffY6BVrb4QbcqAAo7hOgEVh9tOLqL0B3RwHofb 7vY4yIo+ONwpweHy2PFS2I5DXnfKDCMTh3zXscfGu3Co2Sx5EZeAwxuPD050 Nxwcki50ZYfhQHM4PfjGD4dI6X+D9pnhMH1aPURyAQfviE15r6sw+LomiPOf BhwIw8+Zxa1KkGd0Vpf6Bw7FszuWfOqVIKHte4eRZhweL99toluhBM6hBeyC DhwOhO4eFtxRgpLr7GTbYRy6ciT9roFKEL1R2bGMQocdPhM5enJKcD1y1Jxv R4eVLMozZWsa5DxMOlUvpsP0uH75/RpFePOHPHukiQ5RMyKVzDJFCC7y7Ke2 0GHVF/wEr0gRdLpvUZ0ldKgpMdyalawIvffEx8bf00F5oPYE5bAiiPS1Wzxx BuS2bXvEZivCyIWejWrODNAwTYixDlaAtVt3EJkSBlALo58lm1Lg0xSr07eP AT6vLXTcDShQcVC+0OAdAx7bdwSu1KbAdNrDOP4EA7z2l9wVK1HgB21nWq+M AYYzbmsDB+TBUXi8zkNfGX7L0JJMJsrDgqf5BvUTyqD4vnL5wJwcaO9vYq5T ZUL2oHyEsHAJqdIaozo0mOBy/NfGhJQlRPvaEs9ZzgQ972la8JUlZFziGNe6 mglzy+XNth5eQq/cn3THWDCBcdGK42+4hCRn5ju+7WcCr3tBdrp4Ed3Kuf/O M58JbNvazJpyGRrJW808Y6sCwpLoZJv3C8jk3LUw9nUWrMg22JkV/BW566/Y qXFEFb417mmhxswhYaLcD5ryahB546W6l95nJGuN7jmZrQY16w1HIianUMJ5 g4AgEzbY3NxN2735E0qfMDfoLGOD7vncW1dXTSKvH6OKFLaow/wI7rPS+ANa /LhY68NXh41Cwmjo6ShiZFnWvzLSAEPe5+ehR4YR8uXclyVrwHET7Zn2Q++Q uONo7fpPGrDOnn/xxqkBJL6puqbMYRmYXtyitG9zH8q9VNiocmsZnOjc+i3o eg/iCzYMWfQsg8/lH8tCZrrQ2NV7zmZampAwantw4/0/UewbC+F1J03Qu9N0 qE7UjhwWm/wWojUhQDoSMDTWigpnTzwpL9SEbitl77cPXyMFj76P2q80oX5D m1bmRDNyE0hzl09pgiR2AjlbvEDxm7U8/1EmIAulu/lLxcgs3bOvyJAApgEt Ei9tQPPBcSEpQMBqgj80KBChCcN8+QN7CMgwLu4OUKlDU2Zi1V8CCZCmu74I ljxDagXtoaxIAmR6VL85dSHKO2kUI44lIGExodYiowo5ck9vl0siIC+8heFh Won6zi6Of0kn4MuKyHm6Xjnat6sqJr6AABOmg6xCh48+nY3luD4mQOUzi/00 6Aky8MtapiQgIGattc/exkdoVekQ/2wNAQLL723GxQ8Q1bE7qOI5AcpOx3b9 Y1OMPiUtzQibCQhIWakqrS5AlQlR2gptBPx+uHXM++I9JGeenEk+HdQx7OhB OXfR21tnL1X2ELBvW1R2gTgHubg8WP/TWwIK08aebzLPRgrph6uiBwkgfv4i +5iSgUTLsNC0YQKc69X0l9rTkLeo4W/LMQI4vREF0l9T0NWVd3E70mvQMbfx 0BRkI4onnEmHUdO+Gw9IQUmpMfr+pN1osdS5n1JQao+9QiLpt9SPM1uMUlBI hkfXKGm7S9RnSb089HKyfyjxPQF9VOvAzcBDtMpTP6eR5hXeuCmw5KFDO31b 8kjfyZ/KtV7PQ3t5VkVVpEv/PsnZtoKHvgYZrBkm/V1lbGf8LBe9nplHFh8I uGwdIWgt4iL9X8ShiHSPm05pfDYXXWrAUh1JJ3cYRe3gcZGuz9cBL9JH7WTh LbFc9NKv1uoy6eqS4l66DxddCThjd5M0XRi1etidi+rkhDvSSP/I7zGuceYi /67+7Q9Jj2W2Xbhgw0WFpxpAQPqAuWVzkBkX6cVsshCRZs8JW/YYcZGdnv2a l6QTNzQn2OtyUTh3ULWLtO8mMc2S4KL/Ad/oPw0= "]]}}, Axes->True, AxesLabel->{ FormBox["\"N\"", TraditionalForm], FormBox["\"P\"", TraditionalForm]}, AxesOrigin->{0, 0}, Method->{}, PlotRange->{All, All}, PlotRangeClipping->True, PlotRangePadding->{Automatic, Automatic}]], "Output", CellChangeTimes->{3.6623077844093037`*^9, 3.6623089822885103`*^9, 3.6623466933727055`*^9, 3.6623468241906567`*^9, 3.662347843669788*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"u3", "=", RowBox[{"NDSolve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"f", "[", RowBox[{ RowBox[{"n", "[", "t", "]"}], ",", RowBox[{"p", "[", "t", "]"}]}], "]"}], "\[Equal]", RowBox[{ RowBox[{"n", "'"}], "[", "t", "]"}]}], "&&", RowBox[{ RowBox[{"g", "[", RowBox[{ RowBox[{"n", "[", "t", "]"}], ",", RowBox[{"p", "[", "t", "]"}]}], "]"}], "\[Equal]", RowBox[{ RowBox[{"p", "'"}], "[", "t", "]"}]}], "&&", RowBox[{ RowBox[{"n", "[", "0", "]"}], "\[Equal]", "5"}], "&&", RowBox[{ RowBox[{"p", "[", "0", "]"}], "\[Equal]", "6"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"n", "[", "t", "]"}], ",", RowBox[{"p", "[", "t", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "10"}], "}"}]}], "]"}]}], ";"}], "\n", RowBox[{"plot3", "=", RowBox[{"ParametricPlot", "[", RowBox[{ RowBox[{"Evaluate", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"n", "[", "t", "]"}], ",", RowBox[{"p", "[", "t", "]"}]}], "}"}], "/.", "u3"}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "10"}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", "All"}], ",", RowBox[{"PlotStyle", "\[Rule]", "Green"}], ",", RowBox[{"AxesOrigin", "\[Rule]", RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}]}], ",", RowBox[{"AxesLabel", "\[Rule]", RowBox[{"{", RowBox[{"\"\\"", ",", "\"\\""}], "}"}]}]}], "]"}]}]}], "Input", CellChangeTimes->{{3.6622600160741262`*^9, 3.662260043356512*^9}, { 3.6622600965462265`*^9, 3.662260111749984*^9}}], Cell[BoxData[ GraphicsBox[{{}, {}, {RGBColor[0, 1, 0], LineBox[CompressedData[" 1:eJwl13c81d8fB/Br3GFkfC6uFlIRDWSXnLesrFBGdkZCSEYLFWlIRsnKulxF IiqUWUZSaKlsFVIUQrL5ne/jd//5PJ6P87jnfj73nPN+vT8bXI4fOMJOIpFo bCTSf9f/f4Tg/1dR4GufLqxapMM5J38n/ikGrK9Wywido0O9eLeO1BADrtBf 6/JO0WH349ZirxYGfLtcpmE/RIfX7q7pJikMUP78+U1dKx1MfANGTikxIEHy 4froTDpM5eZ/BD8R6DZMP7tnPx0qIkJZ1YtCcDZE8UaFER0OSggvZf4RgpSM kdjdBnR4xZqlXBsUgmb5QmUDHTyf4RrRoBYhkP27/tHNXXRo8aHl3k0TAi2O IxNd0nRYeftVKUBTCNITy9dJkenwQLK+KO4KHcL/BVSIctBBsSXpuVMIHTL3 mzfxsdHhvVX2sPIJOqSta5AiLxHw+waP/4wtHRxYJTFrpgkoVBDoebMDf89G TPD7dwLWRUuINLcTsFfTJrT/FQGT7ffmO3cQoMxmkHaoiYABr9XXdm4mQH5H tfP7RgKUmHatcWsJMPrT/vl1PQHbrsU6O9II8N9tt6W3moDvEbsO7+kXhLci lfahjwigMRsfUZMEYSH8n0paOgEhgdsmSjkFoTR3W7lLGp4/RtG5ZkEAlDv3 3d6aSkBRm9TqlkkBSPSr9WtIJoD1sujgzFcBqBrKUiFuEaBvM51yvRqPRxUe XRdFgOO1t6yiUwJgw0G52nuWgL48AcGLX/ghpLPuxZ8z+Hn2BKlsf8sPoU/k +CnYbj5adV9r+OFop1qNyikCepL1/rlm8INo54B9cQABJwTOyXc58MO7yci1 HD4EnL0QO3Ghjw98DPbqqXgTcP6Y0YeaN3wgIhvtfOwYAe//cj4jP+MDa48M 125P/L9Njt99kMkHi+NcXh/cCVihiK/xduYDsR+9kmudCcgTKOyMG1oFL+Ye pXofJmB7ldS0ffsqiBOU/f3MiQBv3duNCk2rwEmubpW3IwH/SGNB0/dWAeye lu6yw/NH1R4f9FkFxdk6brzWBFjaecjbz/ACK82ZedaKAO6JCbmHP3mh/+OB tBFLvD6c1Jt8Xbygu0x689aCgA7FWM/+Kl7IkCr98+AAARHFah2j4bzwxKHh b81+Ahi9MzqxArxQYkeK1sTu27p0UJmNFyR/aa9+bkIA4Tb5dXCCB97VVFBe GRMQELXK7nAbD9jzm1AmDQmQNtu8qSmRB07nNvKFYS/T4wPLr/DAixPPlwSx qw8XCZacxk66dlXdgIBTCkmrX9jyAI3cJpmpT4DBTFqppTgPHObOk9+D/cVg 9nmcAA9cU73H6NUjwM7R4OdHNh7YyHnZSxI7wtPP//QgNwjt5Ziu0yGAz8zo 9dZ8bsh1jk8/jm0418pdlMoNBRkLgmLYXunRKrujuQF5vrS+oE3ArMh9VV8/ bug2WNpqu5eAxoB6TU1VbqiXVqUKY3tGBirTZLhBRbI9670Wvp/kVcNda7jB y2Izvxn2tts8freXueA9x3KGMxAQ31FnH9bIBR/WehvKYq87jm7HPeWCSvHZ vL8I79/hZve8fC7480H3TAy2S+iWxIkYLrivUbn7syYB7O3BhcPWXHCPa51l AXa3Bke9hiEXjNw4IByBvcGvTj1Zgwskb6hFqGFHTWkreWzgApckmrwwtkq+ RGEPHY+b94RM7SHAT0jlhA2FC+60fJkqwd52JvTFsV80kCKrrU3GPvvgpyFb Hw1y+wZ7Q7AneFfIme9ocPnW2V2u2P8UHCZ06mmwySVazgh7vCKTOlVKgy/h DyuVsHV1XbXv5dFAdGdwrzi2va9O5tFUGoTJmPXMaxDgQXzIWrxAAyHqSPUv 7Ofx4x2PA2nw3Gd6oRc7OsxK0N+DBr4FhgnvsY3S7+ip2tOAk6i72ogdlBcc yGFGA+9RxU9V2JfdTiV+1qZBVPTpMyXYCzdV84pVaWCk5RFYiL1/2vbOja00 iLfubsjFnje+EHlGHI+L53uysIPt9Mw86DSw/VHgzMQ27jeadqDS4N5y5f0M 7OFsuUDbBSpQbUv0/rNaSFqTwzgVvLTCFDKxO3q9xo4OUKFllP9YFvYR2dBv p9upEJxoMpmDvaOuND2umQpN/tKN97DDyTNiRc+oQKq49r0Iu9hio/PHx1QI YLmZPsGmCaw/TMqjgp97Nvtz7PakRoZiGhX0diksv8L2MWaP8I6jgo7VMnzC Fo15xCyIoEL5xL9X37C3VJd7Tp2mwl9tWvI49i/1+S7woYJ2hHj+8n+/52b1 55YzFWaGZVb48Po0phffGbekgkOBUKIEdojur39mhlRoZG/1UcTeUPXzy1NN KpzfBJf0sdldouylFanQsc+u2x77RFKpR5o0FSRyeY4HYFdUGZBE1+H7Gftq mIO95pB6rSSZCvf0vsbVYL+Jsel/OEeB+lIZgS7sFBWfc/pjFHhcZlomhPfr A9tZkcvtFFDSJ3crYvMay/HItVCAeu+OkgX2pvokz77nFADnW+HJ2Kdi9SzN 8inQ0DEcWIWttX+0i8ikQObp7cnfsMNjHCu74ymQHBDtI4fP10V359ALoRSQ lhJS68R23bUjVdicAmzFmokUfD6fvPa9zq1HgbtDAeuVsZ+HL5STd1MgLo2X lYDNt6eigm8zBYLmbY+74fN+lkn89Jojg16F5vxt7P1K82Y3xshw6L7B3Q/Y wolF7dUDZDjpOuCgi+tHaqKOnMwbMqguRMqp4XrjYe6Qto9FBkJiOzkUu6yu 9kBlEhnkT/NursdeNfCbrnydDGK3zmoexPVKyUSOtSeIDDvV/pSE6xIgGRPt k61PhrgSs/L32C6aZpPGGmScO8NLG3A9tM2RsFuWJ4PnFnntRuzbX4RfBq0h g1qgYMb6fXg93Y5cLxvlhNYlhYiT2Eo/7EuT+jkh+4F71jvsYwtjGefbOUGH PzjkKq7PIfk7sl1rOaHlumUCxYgAod87RoMTOMEiZMV8N86HzzqvSoc1OYEv xkoiF9uKZfxJSYkTOn4F04VMCRD7luJ7WYYTPgoxXSewRd/0LmoIccK7noNR 5eYEFKiI8c4Mc4C5zZn5czivsomw9amJHNB8kR6ziD20RrLF6DoH7P3dpBSM 880rhVLMHs4BzHLrJ8E4/zbcOEm96MMBvowXAZE2BKTdnL/YqsMBk569Xz47 4P5B9OPWy3/Z4YePWrQLzlNW5O7kxhF2iHx8R3Ece1enaR3vN3aoqK47yoPz OCpNHfJb2UFw4JC5qQsBh3y3npbJZQfzlIfGnEf/O3+84cp27NDbEmy47gQB S8I5aV4v2aDq8ibNF9iOxtYJdjVswKfRynvcnwAJa7EO81I2uM8eqvoC9xfB Hsvb97PYQGZgJC/4JM77CVb5pfNswMubKy8QQoD4QT/dUTU2OBTTmvP+Cs53 df+Ey0Uk0CHRN1+5ivNObIR0P5cEl149MtoTSYDCjH1+WyYJiqsMI+9fw/kS +qtAIY4Eh5m9j2Kj8XqVPN4icYIEjTtTpU/G4/oeKVfZoEiCQ4z2gbJMfJ7v 81GUT68gefvQ/vNMApJt91266LmCxjO9/Q2ycD1jj//WabuCPj29o92XjfdL E0ssY88KstBaSeG7S8BNviLpCI4V5KSsdi2ygICnwUYOdXHLaOLL2dWD5Xj/ STVmDxUvodSq0MB3FQRkjS+qrGQtocElcdfqStw/nRFxkIhfQvJ15rPJuH9M 9dhlEBK0hEp62katagmIfKfFdXvXErr8MOPwKO5H7d7KeXg0LqKBhkR0voOA KbScKTu4gO7u+7qXf46AsbnPWyVV51GvCYcw1zwB/Q6HTQS3zqMEnUbgWCDg mVtzE0V8HlHuLmXPLuLnfedSzUmdR0cTRsN/kOgQaeW4w7h9Dulr5al/otGh W/raO/KZOSQmzc89KkqHnTU84jy1syhJU+zFpDodnlHGvJNdZpD1QHrdHO7f BWYXG1QOzaAkm02zJA06eOf8c+wymUFLnhlcApp00NDoMFZQn0EWQ48r5PfS YfBVYvV6gRkUU/B6NtgQv294CTmM1fxDe726VlTscN++I3/QRuIfMuUPvLUl lA6NzA3mTyf+ooinzj+1mnDfby2Wuy5jEm3UzhYc5hOCTo/TPJtK/yCnEqP9 jbZCcO/Up4GG+lFkMa3h6ssSAqi7YaqWNoLc55+0vcXvIS1WQofkST9RodXa J5e2CENAbN2rwJDvSDRWvOOgmzD0rm7mprgPIDv1xaAtWcKwS3FMf6fLN2Tm H3yuvEMY/GUM1VvRFzQ15B7XRhOByDZVo4+mPch19QX2GTURWO/jE+sc1Iku VERs13IRgdEn/XEkoh3Rp7ZOBFwRgU6/B0m+0x9RGI1j+mCBCIStOPIwv39A 7RqPlu83i0CUDY+ooux7VCyRLnF1RASM2GajWJNvkOXFDuIXmQErTX5DA/+a kYJTxU1iPQO2qZzwdQp4hWRrs1KaFRlwfmJKuGWxEdmICltM6zFA2i1Ud/Zn AxrYGJ64bM2AvJ3RGyP561G1VP2PRx4MOJBt801esxYZXNqu2HCSATJv1UPZ vJ6hnG4vj0sRDGiUSpYlu1chqVNefqJxDDCh5/N7d5ajnh29ex1vM2AL8mEL iH+Csvpk671zGHCmsjqt3qoUKWocmB8oYECiaawm4f8YddtnZMuUMiDNoY29 M+EhyiFKck9WMeDtYa1unb4iVGdSdzS9jgFzzbdOEEQhWhtAsWM1MUCflbHI aZ6PltIi2spaGeCYqGH+cCEXlZ2LMt71gQGxV6oCJ77fQbXJwtp5nxnQb7F3 c1w7Cz1ry+Y80sWA8pSx1xars9H3K2HfA7GlWqbifwpmI7eKiGcR2CMphd7n ubORahDZKQdbcv+gftlCFhrt1jMewFbUlwnb15eFkoMP9Dt1M6A+ymRsiZWF dM5UCVj0MCCYc71qjkIWOkr+stYVe48R980k2SzEbnlRwh9biRr2KXpjFgp9 aS8Sh41eKgheFM5CD5+wSpqxM+qL1ePmmGhDGbVGq5cBo1G5SKmOiSgRTBcz bAH3yhuWlUw0qDZMcsK+Zf336ZkSJvL+miAbgq3huRj2+i4TtQwW7CzD/p4d 1BMdxURpecJ3GrBr4+tKX0Qw0eYjuoJt2HSfcFvSOSbSNye3j2On1IaxnfNj IlPdVPllbOKhBVutJxNNVnBE8PYxIPpE7UuKKxMdvaX1fg22bPWNg6b2TDT0 2FpUBpuyP5Z125KJ/gfK2/E3 "]]}}, Axes->True, AxesLabel->{ FormBox["\"N\"", TraditionalForm], FormBox["\"P\"", TraditionalForm]}, AxesOrigin->{0, 0}, Method->{}, PlotRange->{All, All}, PlotRangeClipping->True, PlotRangePadding->{Automatic, Automatic}]], "Output", CellChangeTimes->{3.662307784441306*^9, 3.662308985252516*^9, 3.6623466950758915`*^9, 3.662346826237613*^9, 3.6623478461073895`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Show", "[", RowBox[{"plot1", ",", "plot2", ",", "plot3"}], "]"}]], "Input", CellChangeTimes->{{3.66226004852859*^9, 3.662260062888573*^9}}], Cell[BoxData[ GraphicsBox[{{{}, {}, {RGBColor[0, 0, 1], Thickness[Large], LineBox[CompressedData[" 1:eJw1z3k01PsbB/AZWy5iNPENlaWQFgmpG3p/aC+XSglJsuSmlDK0SJSSNZVE dLlKXbeUrRSibClRXUuRnWiMkWUMhiZ37vmd33POc97ndd7n+ePRdD2200OM RqMlifa//N8M4f+5pUw/xVd6CK1bGbF8ziCYvIXf+hlDKFy4Vn6qYRBdGg8v czSHUK3uUrYyYxCrLrw20Fg/BO2L3ZW/OwzCMdB/29HwITg3BewpyPuORAl2 WjI1DB23QKu8swP49Vz6tQcaw7hevNYx9tAAjtam7S7UG8bgocD2c3YDELMY 8ew3HQbr1Emt8wYD8HWuabjkMox3kxzrFT1cpNtyG9IeDKOcs9LTZAcXH1Qq ZsdajqC8stEg1KgfycHPfpP7bQQxQ8st7mv0wyY0fX/knhH0eS9+2jCzH0eN 9zGTjoxAzmhjrCubA9/QoWZa/Aj2BZfmGSZzIDB0mu7ljqBs8L1CqhwHf3KP cBf+wcNd/paOfaNs6HS9OrY3nYcE92fBO7rZeKa953t8Lg9JPQ8Gd9aycVVm SadWFQ8/Us9KB2WxMaChlxM8zoPenoLuM95s7F5hnuRkO4qmxh9kvO8bpBqf 2CQq8FH8vS4nh9uLTYYaXq9V+TD8HD9q3daLlg0+NePafFRdeSs59qEXGRYT x7zM+Dg8OoPlntsLv0Wmy2IO8RFp33+l5IzI8Xo1n8r5MAvUMvGT7YWG+94V LkFjSB823HvZuAfdjjc7VkaOITO5d4G3bg+oaLs2RvwYRrjvs51Ue9DcF1/x 6fEYbCJmCe1pPZg29Fh/p20MoZuDfCqqv0L89nfOKozj4fYwWU/Pr5BcUrbF SWICNM1vFxLudsM7o/bpbcYE3M4787vju+Hxkm/cNXcCpRfOf10V1Y0zyoFe Z1dOAIoe8wWsbngOFaUNHJzAl8mWTR0bu1Ed5futtGoCdwfcH5pyu/DmbmeF 8U0B6OmnlgrMuuD1teyi8l0BZvcU3Nhk2IWeHe8thZkC2BmVZfyh2wWjNqP8 prcCdM69VOcyqwuPvy0+/V4oAOuTp40luxMH124OvX9wEubRlrlH4joR91PL /RezKciz3B+a8DsQQne3ebx1Cm2Bys+TOB14fieeb+8whWhLd550Rwfsl6w+ /NJ/Cmz17i9TVR1Q3SD9+mP2FOy12rMsUzvQThnHiOv9gLiE/Q1rmw7oWmbP ua0mhDDcdbmLaztqG0oejy4WQuexapuaTTt25eS+2rlGCDfXDvSatsPi3lJr dQch7pn39N9XagdLfuZbwU0hmIoy1z69bQOjepdsDvMnRvO9Dn01boOu0px5 8xWn0XvAbJWHUiuk0g/HKGhMg+3MzJaRbEW4IDBIcvk0fq8eflg42gJmROb7 n1bT2ForRif1LbA2m+rUCZ/G/gAd/X9iW6DcPrBdm04jwo9TJTOVWrD8Xq/t twAa+WkvZput3gxNHPC+F0Ij/N2DD54ymlEk5pPkEUkj64YPHC6nN6OQuXI5 N5FGmu9FSdK/fkFQXpXMrAIa+TuA/YKf/gUu1VfMX07QyOYgc93Xxl+g/tSz 98ZJOtEOenG2yrYJjq9dxs2D6eSVsopB7/omnJ4qWNsXRidhn56NyJk0QZ+W 6rkxkU4mOlKDA+Y0oa3gZrdqEZ3wUk8UNrc24rPBupjF4mKkoNEx5k+vRuwQ 2yuovypGmJ/yZkxFfkZ/uNdITbY4MSrOZilwGxAdar+/skCcXKlMqepqb4BW b6ZdaZk4ieII8ovrGjBHZdGXFw3i5HC66+arhQ0ocVTrLxOIk3wDyulWZAMq s/n9lKUEcWmOUzdf1oA0g4SJ8joJohsbd/WWfz0cEph7bk5KkutXJxrNNOog udrP1lxcijAXxKzOVq6DmnuRXI+sFLFv1v9rycw6vJh/86HJPCmSJShLMZys xbCsYFcXpMhu8V8vhdXXInPquL7DJSmSqJhmtz2iFoUttwpzGTOIhO3S4hWT /8DxyH7J64ukyYGpLoutnI9gX826fOOyDFEqz22NL67Bpuf5PiZP5YjLzxUn HM5UIaFCLr+1Qp70yVtVWedWIsx4k6y8L4MU5FxTWq1UgZQPnMMf9iuSbkaa 5hNeKdzf9WXO2zaL3HbaPtdm4BVq1+RpqqoxSXyJBiN8XzG8fVpMjN4wiZTn g1q/gEJon7WpzPaeTZaOr2nNyn0OZ6vSIIfR2WT+I55BiFsevl35ULTjhBKZ U9/4ldJ5gmVxN5glzUrE+68XWt7Lc6Apa78udbkyabWzjawKz8Ln0zbWz88q k66NqYvcWh4B5W6UR74yCaJfYNToPoRDpNmBkn5lsomxTdXhVTq0rWJ9/lKh SFbQ3wq73e8jJfhUcYg5RZrmP260OpiG0MsfGT/2UWRugUnT0PE7SLPnOR87 RZGXYS0s17w/scTGco18NEXGTdUMhIHJECrX7i9MoYiKl1uFSVISrtBvbNr4 mCIJbAutNL1bCF7n3Le+gCK9A0ciKiZv4v7F38uWVFDkYIrD2di6G3ipozSR WUORgXdOtx1psSh+e8jySANFfrt/ssPG+BoG9/3x7nQzRbaWUhoZTjGIE+wc O9dBEbm9Ygd97kbhY+tYiHIPRWzMMzOLeOG4lMLwfMKmyHSJG8s34zK0KOE1 Ca7ov8XbRlYdvQST2gG/0u8U+ef+lV1Xj4XA9fqG54NDFFk4mphpUXweXMWZ MrwRiuir2e2VYwaj5qgBL3qUIhrhKiXmOwLx5HLyvj4+Rayl9aniN2eA1E+N y8YpUrll96WTm09B9vh+ruEERYQa1XLqev4I8Y1ctErkWh3Wgtka/ijNUnQ1 FTlvgyJXjvJH/UXjmnUiq9rKR/8i6Q+20S8Xd4msYmXVrd/ph885SX5+Ips9 uR2xKNEPXaZOfk9F9jUDqZH3Q5TNYZd8kV/Twx6slPLDIS35LUUiO3/2j08V sqCZocSsEJk/ksAL4bJg2eIYUS9y3N/60SeqWKApGc7giZwlljxjLJSFsKT1 1WP/3Vv2D8edY+HjGcmYSZFnXNPd/Ks/C9KchJliAopwEl8aRniw8IZaOcgQ +dX8R4271rOgMMxOnf1fH5szqWDGQtsx+51zRD4Svy3wvREL/hLHM9RFnh0e kbx7AQvHvefZLRCZEe+wZr4aC3u0ooQ6Is/ylV/BmcVCSXXxncUi0x+Zn8+X YWGz8rMN+iKX1QfoRYmx8C8EjKa6 "]]}}, {{}, {}, {RGBColor[1, 0, 0], Dashing[{Small, Small}], LineBox[CompressedData[" 1:eJwl1Is31HkfB3CGGb/fjDEM+aljpchD6SG5bcr3E61SbIsS7cq1JKIL2TLR RSkVxTaD3JNb6WJcxprC0Nio3JZBQm6FpkVsyQ7P7zn7Pud93uf1D7xX+YW5 HqTIyclFkP3//hsM/l1V6Gu/GiulYMB95DF5UMqC+v34bB8VA5XOh8HdvSzg SLxT6xgYqPq3bRkuZ0HlxYUPOQQGZXvczBNDWGC75Dr61ASDjrARJes+FeBQ oqq9fDGYXO07d6OBCcyOxaNxARjwF/9Si+Qz4dK1otnSQAx0XAIqgnOZcC/f MIoVioE0ihMeeY4JWx/l7e0/g8EZFNJJbGGCd2pFbPtvGISbGkvCBcpwftGh T6UJA1u6he2hSgZInc4d+f0VBpRrawMO5DNg+5c3i4FtGIimZUN+txmw/9EV ixfdGJSetKq/eZIBCiM21XffY6BVrb4QbcqAAo7hOgEVh9tOLqL0B3RwHofb 7vY4yIo+ONwpweHy2PFS2I5DXnfKDCMTh3zXscfGu3Co2Sx5EZeAwxuPD050 Nxwcki50ZYfhQHM4PfjGD4dI6X+D9pnhMH1aPURyAQfviE15r6sw+LomiPOf BhwIw8+Zxa1KkGd0Vpf6Bw7FszuWfOqVIKHte4eRZhweL99toluhBM6hBeyC DhwOhO4eFtxRgpLr7GTbYRy6ciT9roFKEL1R2bGMQocdPhM5enJKcD1y1Jxv R4eVLMozZWsa5DxMOlUvpsP0uH75/RpFePOHPHukiQ5RMyKVzDJFCC7y7Ke2 0GHVF/wEr0gRdLpvUZ0ldKgpMdyalawIvffEx8bf00F5oPYE5bAiiPS1Wzxx BuS2bXvEZivCyIWejWrODNAwTYixDlaAtVt3EJkSBlALo58lm1Lg0xSr07eP AT6vLXTcDShQcVC+0OAdAx7bdwSu1KbAdNrDOP4EA7z2l9wVK1HgB21nWq+M AYYzbmsDB+TBUXi8zkNfGX7L0JJMJsrDgqf5BvUTyqD4vnL5wJwcaO9vYq5T ZUL2oHyEsHAJqdIaozo0mOBy/NfGhJQlRPvaEs9ZzgQ972la8JUlZFziGNe6 mglzy+XNth5eQq/cn3THWDCBcdGK42+4hCRn5ju+7WcCr3tBdrp4Ed3Kuf/O M58JbNvazJpyGRrJW808Y6sCwpLoZJv3C8jk3LUw9nUWrMg22JkV/BW566/Y qXFEFb417mmhxswhYaLcD5ryahB546W6l95nJGuN7jmZrQY16w1HIianUMJ5 g4AgEzbY3NxN2735E0qfMDfoLGOD7vncW1dXTSKvH6OKFLaow/wI7rPS+ANa /LhY68NXh41Cwmjo6ShiZFnWvzLSAEPe5+ehR4YR8uXclyVrwHET7Zn2Q++Q uONo7fpPGrDOnn/xxqkBJL6puqbMYRmYXtyitG9zH8q9VNiocmsZnOjc+i3o eg/iCzYMWfQsg8/lH8tCZrrQ2NV7zmZampAwantw4/0/UewbC+F1J03Qu9N0 qE7UjhwWm/wWojUhQDoSMDTWigpnTzwpL9SEbitl77cPXyMFj76P2q80oX5D m1bmRDNyE0hzl09pgiR2AjlbvEDxm7U8/1EmIAulu/lLxcgs3bOvyJAApgEt Ei9tQPPBcSEpQMBqgj80KBChCcN8+QN7CMgwLu4OUKlDU2Zi1V8CCZCmu74I ljxDagXtoaxIAmR6VL85dSHKO2kUI44lIGExodYiowo5ck9vl0siIC+8heFh Won6zi6Of0kn4MuKyHm6Xjnat6sqJr6AABOmg6xCh48+nY3luD4mQOUzi/00 6Aky8MtapiQgIGattc/exkdoVekQ/2wNAQLL723GxQ8Q1bE7qOI5AcpOx3b9 Y1OMPiUtzQibCQhIWakqrS5AlQlR2gptBPx+uHXM++I9JGeenEk+HdQx7OhB OXfR21tnL1X2ELBvW1R2gTgHubg8WP/TWwIK08aebzLPRgrph6uiBwkgfv4i +5iSgUTLsNC0YQKc69X0l9rTkLeo4W/LMQI4vREF0l9T0NWVd3E70mvQMbfx 0BRkI4onnEmHUdO+Gw9IQUmpMfr+pN1osdS5n1JQao+9QiLpt9SPM1uMUlBI hkfXKGm7S9RnSb089HKyfyjxPQF9VOvAzcBDtMpTP6eR5hXeuCmw5KFDO31b 8kjfyZ/KtV7PQ3t5VkVVpEv/PsnZtoKHvgYZrBkm/V1lbGf8LBe9nplHFh8I uGwdIWgt4iL9X8ShiHSPm05pfDYXXWrAUh1JJ3cYRe3gcZGuz9cBL9JH7WTh LbFc9NKv1uoy6eqS4l66DxddCThjd5M0XRi1etidi+rkhDvSSP/I7zGuceYi /67+7Q9Jj2W2Xbhgw0WFpxpAQPqAuWVzkBkX6cVsshCRZs8JW/YYcZGdnv2a l6QTNzQn2OtyUTh3ULWLtO8mMc2S4KL/Ad/oPw0= "]]}}, {{}, {}, {RGBColor[0, 1, 0], LineBox[CompressedData[" 1:eJwl13c81d8fB/Br3GFkfC6uFlIRDWSXnLesrFBGdkZCSEYLFWlIRsnKulxF IiqUWUZSaKlsFVIUQrL5ne/jd//5PJ6P87jnfj73nPN+vT8bXI4fOMJOIpFo bCTSf9f/f4Tg/1dR4GufLqxapMM5J38n/ikGrK9Wywido0O9eLeO1BADrtBf 6/JO0WH349ZirxYGfLtcpmE/RIfX7q7pJikMUP78+U1dKx1MfANGTikxIEHy 4froTDpM5eZ/BD8R6DZMP7tnPx0qIkJZ1YtCcDZE8UaFER0OSggvZf4RgpSM kdjdBnR4xZqlXBsUgmb5QmUDHTyf4RrRoBYhkP27/tHNXXRo8aHl3k0TAi2O IxNd0nRYeftVKUBTCNITy9dJkenwQLK+KO4KHcL/BVSIctBBsSXpuVMIHTL3 mzfxsdHhvVX2sPIJOqSta5AiLxHw+waP/4wtHRxYJTFrpgkoVBDoebMDf89G TPD7dwLWRUuINLcTsFfTJrT/FQGT7ffmO3cQoMxmkHaoiYABr9XXdm4mQH5H tfP7RgKUmHatcWsJMPrT/vl1PQHbrsU6O9II8N9tt6W3moDvEbsO7+kXhLci lfahjwigMRsfUZMEYSH8n0paOgEhgdsmSjkFoTR3W7lLGp4/RtG5ZkEAlDv3 3d6aSkBRm9TqlkkBSPSr9WtIJoD1sujgzFcBqBrKUiFuEaBvM51yvRqPRxUe XRdFgOO1t6yiUwJgw0G52nuWgL48AcGLX/ghpLPuxZ8z+Hn2BKlsf8sPoU/k +CnYbj5adV9r+OFop1qNyikCepL1/rlm8INo54B9cQABJwTOyXc58MO7yci1 HD4EnL0QO3Ghjw98DPbqqXgTcP6Y0YeaN3wgIhvtfOwYAe//cj4jP+MDa48M 125P/L9Njt99kMkHi+NcXh/cCVihiK/xduYDsR+9kmudCcgTKOyMG1oFL+Ye pXofJmB7ldS0ffsqiBOU/f3MiQBv3duNCk2rwEmubpW3IwH/SGNB0/dWAeye lu6yw/NH1R4f9FkFxdk6brzWBFjaecjbz/ACK82ZedaKAO6JCbmHP3mh/+OB tBFLvD6c1Jt8Xbygu0x689aCgA7FWM/+Kl7IkCr98+AAARHFah2j4bzwxKHh b81+Ahi9MzqxArxQYkeK1sTu27p0UJmNFyR/aa9+bkIA4Tb5dXCCB97VVFBe GRMQELXK7nAbD9jzm1AmDQmQNtu8qSmRB07nNvKFYS/T4wPLr/DAixPPlwSx qw8XCZacxk66dlXdgIBTCkmrX9jyAI3cJpmpT4DBTFqppTgPHObOk9+D/cVg 9nmcAA9cU73H6NUjwM7R4OdHNh7YyHnZSxI7wtPP//QgNwjt5Ziu0yGAz8zo 9dZ8bsh1jk8/jm0418pdlMoNBRkLgmLYXunRKrujuQF5vrS+oE3ArMh9VV8/ bug2WNpqu5eAxoB6TU1VbqiXVqUKY3tGBirTZLhBRbI9670Wvp/kVcNda7jB y2Izvxn2tts8freXueA9x3KGMxAQ31FnH9bIBR/WehvKYq87jm7HPeWCSvHZ vL8I79/hZve8fC7480H3TAy2S+iWxIkYLrivUbn7syYB7O3BhcPWXHCPa51l AXa3Bke9hiEXjNw4IByBvcGvTj1Zgwskb6hFqGFHTWkreWzgApckmrwwtkq+ RGEPHY+b94RM7SHAT0jlhA2FC+60fJkqwd52JvTFsV80kCKrrU3GPvvgpyFb Hw1y+wZ7Q7AneFfIme9ocPnW2V2u2P8UHCZ06mmwySVazgh7vCKTOlVKgy/h DyuVsHV1XbXv5dFAdGdwrzi2va9O5tFUGoTJmPXMaxDgQXzIWrxAAyHqSPUv 7Ofx4x2PA2nw3Gd6oRc7OsxK0N+DBr4FhgnvsY3S7+ip2tOAk6i72ogdlBcc yGFGA+9RxU9V2JfdTiV+1qZBVPTpMyXYCzdV84pVaWCk5RFYiL1/2vbOja00 iLfubsjFnje+EHlGHI+L53uysIPt9Mw86DSw/VHgzMQ27jeadqDS4N5y5f0M 7OFsuUDbBSpQbUv0/rNaSFqTwzgVvLTCFDKxO3q9xo4OUKFllP9YFvYR2dBv p9upEJxoMpmDvaOuND2umQpN/tKN97DDyTNiRc+oQKq49r0Iu9hio/PHx1QI YLmZPsGmCaw/TMqjgp97Nvtz7PakRoZiGhX0diksv8L2MWaP8I6jgo7VMnzC Fo15xCyIoEL5xL9X37C3VJd7Tp2mwl9tWvI49i/1+S7woYJ2hHj+8n+/52b1 55YzFWaGZVb48Po0phffGbekgkOBUKIEdojur39mhlRoZG/1UcTeUPXzy1NN KpzfBJf0sdldouylFanQsc+u2x77RFKpR5o0FSRyeY4HYFdUGZBE1+H7Gftq mIO95pB6rSSZCvf0vsbVYL+Jsel/OEeB+lIZgS7sFBWfc/pjFHhcZlomhPfr A9tZkcvtFFDSJ3crYvMay/HItVCAeu+OkgX2pvokz77nFADnW+HJ2Kdi9SzN 8inQ0DEcWIWttX+0i8ikQObp7cnfsMNjHCu74ymQHBDtI4fP10V359ALoRSQ lhJS68R23bUjVdicAmzFmokUfD6fvPa9zq1HgbtDAeuVsZ+HL5STd1MgLo2X lYDNt6eigm8zBYLmbY+74fN+lkn89Jojg16F5vxt7P1K82Y3xshw6L7B3Q/Y wolF7dUDZDjpOuCgi+tHaqKOnMwbMqguRMqp4XrjYe6Qto9FBkJiOzkUu6yu 9kBlEhnkT/NursdeNfCbrnydDGK3zmoexPVKyUSOtSeIDDvV/pSE6xIgGRPt k61PhrgSs/L32C6aZpPGGmScO8NLG3A9tM2RsFuWJ4PnFnntRuzbX4RfBq0h g1qgYMb6fXg93Y5cLxvlhNYlhYiT2Eo/7EuT+jkh+4F71jvsYwtjGefbOUGH PzjkKq7PIfk7sl1rOaHlumUCxYgAod87RoMTOMEiZMV8N86HzzqvSoc1OYEv xkoiF9uKZfxJSYkTOn4F04VMCRD7luJ7WYYTPgoxXSewRd/0LmoIccK7noNR 5eYEFKiI8c4Mc4C5zZn5czivsomw9amJHNB8kR6ziD20RrLF6DoH7P3dpBSM 880rhVLMHs4BzHLrJ8E4/zbcOEm96MMBvowXAZE2BKTdnL/YqsMBk569Xz47 4P5B9OPWy3/Z4YePWrQLzlNW5O7kxhF2iHx8R3Ece1enaR3vN3aoqK47yoPz OCpNHfJb2UFw4JC5qQsBh3y3npbJZQfzlIfGnEf/O3+84cp27NDbEmy47gQB S8I5aV4v2aDq8ibNF9iOxtYJdjVswKfRynvcnwAJa7EO81I2uM8eqvoC9xfB Hsvb97PYQGZgJC/4JM77CVb5pfNswMubKy8QQoD4QT/dUTU2OBTTmvP+Cs53 df+Ey0Uk0CHRN1+5ivNObIR0P5cEl149MtoTSYDCjH1+WyYJiqsMI+9fw/kS +qtAIY4Eh5m9j2Kj8XqVPN4icYIEjTtTpU/G4/oeKVfZoEiCQ4z2gbJMfJ7v 81GUT68gefvQ/vNMApJt91266LmCxjO9/Q2ycD1jj//WabuCPj29o92XjfdL E0ssY88KstBaSeG7S8BNviLpCI4V5KSsdi2ygICnwUYOdXHLaOLL2dWD5Xj/ STVmDxUvodSq0MB3FQRkjS+qrGQtocElcdfqStw/nRFxkIhfQvJ15rPJuH9M 9dhlEBK0hEp62katagmIfKfFdXvXErr8MOPwKO5H7d7KeXg0LqKBhkR0voOA KbScKTu4gO7u+7qXf46AsbnPWyVV51GvCYcw1zwB/Q6HTQS3zqMEnUbgWCDg mVtzE0V8HlHuLmXPLuLnfedSzUmdR0cTRsN/kOgQaeW4w7h9Dulr5al/otGh W/raO/KZOSQmzc89KkqHnTU84jy1syhJU+zFpDodnlHGvJNdZpD1QHrdHO7f BWYXG1QOzaAkm02zJA06eOf8c+wymUFLnhlcApp00NDoMFZQn0EWQ48r5PfS YfBVYvV6gRkUU/B6NtgQv294CTmM1fxDe726VlTscN++I3/QRuIfMuUPvLUl lA6NzA3mTyf+ooinzj+1mnDfby2Wuy5jEm3UzhYc5hOCTo/TPJtK/yCnEqP9 jbZCcO/Up4GG+lFkMa3h6ssSAqi7YaqWNoLc55+0vcXvIS1WQofkST9RodXa J5e2CENAbN2rwJDvSDRWvOOgmzD0rm7mprgPIDv1xaAtWcKwS3FMf6fLN2Tm H3yuvEMY/GUM1VvRFzQ15B7XRhOByDZVo4+mPch19QX2GTURWO/jE+sc1Iku VERs13IRgdEn/XEkoh3Rp7ZOBFwRgU6/B0m+0x9RGI1j+mCBCIStOPIwv39A 7RqPlu83i0CUDY+ooux7VCyRLnF1RASM2GajWJNvkOXFDuIXmQErTX5DA/+a kYJTxU1iPQO2qZzwdQp4hWRrs1KaFRlwfmJKuGWxEdmICltM6zFA2i1Ud/Zn AxrYGJ64bM2AvJ3RGyP561G1VP2PRx4MOJBt801esxYZXNqu2HCSATJv1UPZ vJ6hnG4vj0sRDGiUSpYlu1chqVNefqJxDDCh5/N7d5ajnh29ex1vM2AL8mEL iH+Csvpk671zGHCmsjqt3qoUKWocmB8oYECiaawm4f8YddtnZMuUMiDNoY29 M+EhyiFKck9WMeDtYa1unb4iVGdSdzS9jgFzzbdOEEQhWhtAsWM1MUCflbHI aZ6PltIi2spaGeCYqGH+cCEXlZ2LMt71gQGxV6oCJ77fQbXJwtp5nxnQb7F3 c1w7Cz1ry+Y80sWA8pSx1xars9H3K2HfA7GlWqbifwpmI7eKiGcR2CMphd7n ubORahDZKQdbcv+gftlCFhrt1jMewFbUlwnb15eFkoMP9Dt1M6A+ymRsiZWF dM5UCVj0MCCYc71qjkIWOkr+stYVe48R980k2SzEbnlRwh9biRr2KXpjFgp9 aS8Sh41eKgheFM5CD5+wSpqxM+qL1ePmmGhDGbVGq5cBo1G5SKmOiSgRTBcz bAH3yhuWlUw0qDZMcsK+Zf336ZkSJvL+miAbgq3huRj2+i4TtQwW7CzD/p4d 1BMdxURpecJ3GrBr4+tKX0Qw0eYjuoJt2HSfcFvSOSbSNye3j2On1IaxnfNj IlPdVPllbOKhBVutJxNNVnBE8PYxIPpE7UuKKxMdvaX1fg22bPWNg6b2TDT0 2FpUBpuyP5Z125KJ/gfK2/E3 "]]}}}, Axes->True, AxesLabel->{ FormBox["\"N\"", TraditionalForm], FormBox["\"P\"", TraditionalForm]}, AxesOrigin->{0, 0}, Method->{}, PlotRange->{All, All}, PlotRangeClipping->True, PlotRangePadding->{Automatic, Automatic}]], "Output", CellChangeTimes->{ 3.662307784461307*^9, {3.662308988996522*^9, 3.6623090148145676`*^9}, 3.662346696966625*^9, 3.662346828034563*^9, 3.662347849201268*^9}] }, Open ]] }, WindowSize->{1546, 847}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (January 25, 2013)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[579, 22, 1971, 59, 118, "Input"], Cell[2553, 83, 852, 24, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[3442, 112, 1103, 33, 52, "Input"], Cell[4548, 147, 2126, 62, 82, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[6711, 214, 184, 4, 31, "Input"], Cell[6898, 220, 394, 11, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[7329, 236, 189, 3, 31, "Input"], Cell[7521, 241, 235, 5, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[7793, 251, 423, 8, 31, "Input"], Cell[8219, 261, 434, 12, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[8690, 278, 184, 3, 31, "Input"], Cell[8877, 283, 171, 4, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[9085, 292, 184, 4, 31, "Input"], Cell[9272, 298, 342, 10, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[9651, 313, 136, 2, 31, "Input"], Cell[9790, 317, 175, 5, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[10002, 327, 217, 4, 31, "Input"], Cell[10222, 333, 350, 10, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[10609, 348, 134, 2, 31, "Input"], Cell[10746, 352, 187, 5, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[10970, 362, 1922, 54, 52, "Input"], Cell[12895, 418, 3795, 69, 109, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[16727, 492, 1826, 53, 52, "Input"], Cell[18556, 547, 2998, 55, 275, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[21591, 607, 1793, 51, 52, "Input"], Cell[23387, 660, 5820, 102, 386, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[29244, 767, 167, 3, 31, "Input"], Cell[29414, 772, 11784, 202, 416, "Output"] }, Open ]] } ] *) (* End of internal cache information *)