(* 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[ 44323, 1021] NotebookOptionsPosition[ 42986, 969] NotebookOutlinePosition[ 43329, 984] CellTagsIndexPosition[ 43286, 981] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"f", "[", RowBox[{"n_", ",", "p_"}], "]"}], ":=", RowBox[{ RowBox[{"1.7", "n", RowBox[{"(", RowBox[{"20", "-", "n"}], ")"}]}], "-", FractionBox[ RowBox[{"2", "n", "*", "p"}], RowBox[{"(", RowBox[{"5", "+", RowBox[{"2.2", SuperscriptBox["n", "2"]}]}], ")"}]]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"g", "[", RowBox[{"n_", ",", "p_"}], "]"}], ":=", RowBox[{ FractionBox[ RowBox[{"2", "n", "*", "p"}], RowBox[{"(", RowBox[{"5", "+", RowBox[{"2.2", SuperscriptBox["n", "2"]}]}], ")"}]], "-", RowBox[{"0.8", "p"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Solve", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"1.7", "n", RowBox[{"(", RowBox[{"20", "-", "n"}], ")"}]}], "-", FractionBox[ RowBox[{"2", "n", "*", "p"}], RowBox[{"(", RowBox[{"5", "+", RowBox[{"2.2", SuperscriptBox["n", "2"]}]}], ")"}]]}], "\[Equal]", "0"}], "&&", RowBox[{ RowBox[{ FractionBox[ RowBox[{"2", "n", "*", "p"}], RowBox[{"(", RowBox[{"5", "+", RowBox[{"2.2", SuperscriptBox["n", "2"]}]}], ")"}]], "-", RowBox[{"0.8", "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.6623479198448353`*^9, 3.662347978284768*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"n", "\[Rule]", "0"}], ",", RowBox[{"p", "\[Rule]", "0"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"n", "\[Rule]", "20.`"}], ",", RowBox[{"p", "\[Rule]", "0.`"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.662307784210293*^9, 3.6623079912581353`*^9, 3.6623466557930174`*^9, 3.6623479798785834`*^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[{ RowBox[{"1.7`", " ", RowBox[{"(", RowBox[{"20", "-", "n"}], ")"}]}], "-", RowBox[{"1.7`", " ", "n"}], "+", FractionBox[ RowBox[{"8.8`", " ", SuperscriptBox["n", "2"], " ", "p"}], SuperscriptBox[ RowBox[{"(", RowBox[{"5", "+", RowBox[{"2.2`", " ", SuperscriptBox["n", "2"]}]}], ")"}], "2"]], "-", FractionBox[ RowBox[{"2", " ", "p"}], RowBox[{"5", "+", RowBox[{"2.2`", " ", SuperscriptBox["n", "2"]}]}]]}], RowBox[{"-", FractionBox[ RowBox[{"2", " ", "n"}], RowBox[{"5", "+", RowBox[{"2.2`", " ", SuperscriptBox["n", "2"]}]}]]}]}, { RowBox[{ RowBox[{"-", FractionBox[ RowBox[{"8.8`", " ", SuperscriptBox["n", "2"], " ", "p"}], SuperscriptBox[ RowBox[{"(", RowBox[{"5", "+", RowBox[{"2.2`", " ", SuperscriptBox["n", "2"]}]}], ")"}], "2"]]}], "+", FractionBox[ RowBox[{"2", " ", "p"}], RowBox[{"5", "+", RowBox[{"2.2`", " ", SuperscriptBox["n", "2"]}]}]]}], RowBox[{ RowBox[{"-", "0.8`"}], "+", FractionBox[ RowBox[{"2", " ", "n"}], RowBox[{"5", "+", RowBox[{"2.2`", " ", SuperscriptBox["n", "2"]}]}]]}]} }, 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.6623479834412217`*^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[{"34.`", ",", "0.`"}], "}"}], ",", RowBox[{"{", RowBox[{"0.`", ",", RowBox[{"-", "0.8`"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{ 3.6623077842782965`*^9, 3.66230802421402*^9, {3.662308708231144*^9, 3.6623087112483163`*^9}, 3.662346663980857*^9, 3.6623479851288013`*^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[{"34.`", ",", RowBox[{"-", "0.8`"}]}], "}"}]], "Output", CellChangeTimes->{3.6623077842972975`*^9, 3.6623087688489356`*^9, 3.6623466654184256`*^9, 3.662347986910127*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ss2", "=", RowBox[{"jacobianMatrix", "[", RowBox[{"20", ",", "0"}], "]"}]}]], "Input", CellChangeTimes->{{3.6623089054271755`*^9, 3.6623089058483763`*^9}, { 3.6623089521180573`*^9, 3.6623089557372637`*^9}, {3.6623467090764847`*^9, 3.6623467103265257`*^9}, 3.6623479908946695`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "34.`"}], ",", RowBox[{"-", "0.04519774011299434`"}]}], "}"}], ",", RowBox[{"{", RowBox[{"0.`", ",", RowBox[{"-", "0.7548022598870057`"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{ 3.66230895941887*^9, {3.6623466874974623`*^9, 3.6623467171236954`*^9}, 3.662347992113467*^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[{ RowBox[{"-", "34.`"}], ",", RowBox[{"-", "0.7548022598870057`"}]}], "}"}]], "Output", CellChangeTimes->{ 3.6623089686384864`*^9, {3.662346689200658*^9, 3.6623467182174892`*^9}, 3.662347993988551*^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:eJw11Hc8Vu/7AHCrKCFKKYq0kJUjLVzH6GNmy3iMT8p6VAqFSqRkpcJHQoOQ WUJGVsrOikhSVjzL83iMrFC/c8739Tv/3K/36x7Xue/rel27XLwtXLk4ODgU ODk48PF/3xT8//i5S0nraeY0WApcK+GeZcPnD1WGEn6z4F15aI0OhQ2uovmO oxfmwE/EXPpVPxsCZy9uupm0ANt+N44h7WzoldILz1tYgh2n1RM/1LJBuF07 WClmBXQWrrfrv2FD+hHNzGE1TtQnItVgJZkNKbHqLrrDXGjHy6puwXtssGF/ WVrN50HnrxrIjIWwIUAx0s25jA9NTzj9ttqVDcpcrTtFF9ej7p9IBpU2bGDk 5FnvchFAfVsfV543wNYfVjUusBdGhVciZhcOsOFPkJVsfKcIarD9cE+vBBvW fvJRzrXbjMY37VDyE2AD/0md64Z5W1FvVf17PJOTYPhAcX/GHkm0JPJUwXz+ JPDcGn4w0LIfvTttwf187yTM3FTe8W+zDOrv2Sh1ZfMkvK9p4n5eK4t+OCN1 VId7EjLuZIZA8gF0IFNhuWOYBZMvRZvGAuVRQeq1+AedLLiALAmMmiqgFVNr ec1rWNBwlDy0+4ciCgff0LtTWDCaoudAiVdCVdaf0HsYxQJPcbudy5LKaK62 0RmHQBZUrHfcXXhRGbWaDVje68GCsry5Y3JVymhG6YG3s6dYcEk65k0Y10F0 ONn2UOMJFnTGfZZg6hxE/Vzkl56osuCjRqZ4avlB1MpSCXEWYcEZiWU/B+ZB tE+TFmjAyQKTDX/m/oiroAIb9uQem2LCQMBr4QY9FXRn5l0RZIgJD52aIpje KmhTzcdFlQ4mSFZ9vVP0nwp6Tls75Vg1E3YNJglalqqgTgefHDTMZ0KScpYK X68Kqie9wDydwoQZs6D+ddMqqLHZ7qWbUUyYlTb2jVmPoH27+IJyA5lwU4w/ tGEXgl5WbU4e8GDCFTfuhhE1BI159fqiqC0T5CsNlIUMEVSyLmqbrR4TXnrV H/cnYfPqJvnpakwIyXpsoeWFoEWW3zQW9zLB3a1pb2wggr75eH/IWpQJXZa1 9yLuIOiGy6ppVTxMUDmc6XIoHkGP9rpbNc9MwAzCik14iqD+WlODg8MT0JHw K64lG0G7nWONOTonwCNDMJZdiKAuDZey5KsngCnq4iJdgaAfh/0WzuRNgJYk Z93V91i8WUmdF0kT4KIn1iLYjKD83Yb3Z8MnAEUqHsy0I6h0fuKQ0ZUJ6FJu Uz3yGUG1b19WKzg7ARecNefZfQjaHAsJOywnIPGr2mOZ7wgadCFw5aHWBDy6 rabJHkLQY70lXuLKExC65vwP059Y/LXskdydE0C+F8dlTUHQ1BwrxxMCE6DX K/svHx1BOWUyBxjLDDigHu3jM4GghdardikMBhjtP78uhYWgyk67ek/1M8CM o4ASwUbQgXcnDXY0M8B3UstVexpBi29FF02WMuB6oTqjbQZBw8bJG1ozGdBI okQo/EJQEy5ri6L/GKCQHf7MYw5BrYMdgjNuMUBEhK4cOo+g2RKe99J8GDC6 NrQ1aAFB27QUruacZoAayTTRaRFBxWPlNarMGPC8pyJg/xKCPhxq6+oHBvBl zTO/Yf7+4R8VLiUGJJpIGV77jaA/kQ+2h3Yy4NX591Hrl7H8j6LgK8AA/q/x d6Mwt1eUD1Wu0GGTPT91BXM2jVtViEmHC6+iNp1ZQVDf2alD5wfowOv6wKEG c8x+jaHej3RwulTCEFpFUPu0JFmDCjqgHpIfbTHD/EvBphw6nJeRHEvCnETe EmGaRIdKw8vCPZgzbcj3RyPoMErOXM/3B0EP/D28MySADi+bI4LUMB9mrleU 8aDDmKbnqjNm0n2P8m82dBBs19schtnNr6c0UY8Oe8xpR15gDlvukHY+TIdl 0dW+Osyz4n1zyvvp8HVzo/kgZq86sqzAVjo8LE9fnsP8Q4e3bHYtHZJ67u3m /4ugAQJyj8bmafDWrWrDTsz3dpo3DlFoIFJ+YV4Rs+jv38fGv9BA5o+cqCZm Z/ls9q9GGgjFTCYbYW7v6f8hWEaDcoWuYhvM48c3/VXJokFkWmSYC+blom/G pxNpcHLgruo5zLwSZbVJ4TSoDl4z5If52SnE6rs/DTj9KlOv4679ziHrQQM/ 10vxtzBrbNeuDbalgTOi2haJuersQvSwPg28fp51vo/53VSso9FRGgR/rnP6 D/Nfcq5irSwNckvj+x9h/uUyuKi5nQZx958PPsZMN2wraVpPA1Ki1NVUzD2r W53tl6nw+umJsnTMlskek/MTVJCTiUp9gVmp0s3xyXcqlB5R2JyD+czkCzHz DiqEUMuM8zA/+IfPVaCWChwXb5q9xHxF3VOju5AKR9UC9xRgjt326vGzdCoY W+i1v8ZsyKgOupxAhXXad0yLMAfTgwcsw6nguGv1VTFmrpKh18cCqVBd9mvy DWaBrPZFWS8q8L8uFSnFnGFz7NUuRyr86ym3owyz2Bextl2mVHjsnbOxHPPQ GZeTclpU+B4Swcbtd1hY4ThCBYMby5VvMYcES7lZ7aXCEqvMvwLz84zY2ctb qTCicbeuEvNaNGeha5kCOT4xp6owR9U7nN0wSQEhx8hB3PNXHgibDlMgvemF QzVm9pD576RuChh40rpwiz57vY1VT4FvT8ShBvOen3Vk/TIKJHmtz8T9zSCD lpdDAROSqtM7zCWlyjbR9yhA6sgtwi1ZsE2T9yYFDvFactVilrU21Iz2pYDE MGcC7hGNrf55thRInGF9xe3NXvNCz4gCx63zt73H/Gdy7NuEBgXqk+vicJO2 H9Ux3k2B3XUlrbj9rbac59tCgVqNWK4P+PyJisg2PgpU6jWRcdNVjSPOTY6D /Zsh2TrMrs0jovwV40DuvHC1HvOLM042NRbj8OxtYRruUJ8zAQsnxsGuNKIZ 9znQ0zpydBza57eKNmA+P2ZZ/klyHJRP7M3FnR2rWuXAGgPZ+xeaGzGbmynt /RE5Bt87Cxm4pZ3eO7UFjUGCZdeGJswfXna++nBpDILrc81w83S036y3G4MP RylfcF++lBIsITcGLvKOg82Y2zpWtLa3/YS+KyPfP+L1/M7imbLoT8hbffu0 AzN/NvufyYoROETLovVi3mEqPRR9fBDk4gsrBvB8uTgW7RH4Br4v26RHMLfQ A/1TDnyBTqOqrHHMyu5GYMbshpm++hN0PJ8rPSNrLD/BUgZnOhN/nzaOacqn Vqheem7JxuspMy2mbbYJ+BtKkqYxJ1AVOtIL6qF1OLtsFq/n7T9sv2e8BzfD D+Q5vB+Ff4/N86qB6beklXnMU7nBQbaUt5C3m3xyEbOuQWJoxo5SyORzJS3h /cX4XTCypRhibGJP/cb7QcJV7hrd1/CtXkdyGfPetuD+BN58YIal9OM+H+qx f/lrFniLlV9ewdx3ePgk2S4Ddv3eN4P7murt7a1308BTi9ttFfNHhWm5O+NP QCy3pR43kkp2el2VBII0buE/mHevWkmUHXkIi11i9rgjTbhl+mfiQE2yPh73 NgfmkUyrB9DWIPwR98sMdf+etmgo5IleINy12HipMhxawtJ2/MWsmI62TgXe hrAb3pq4SU6OlQ0+N+GTlJAtbntZZ3mWwA3Iat/gjVvm4WJ831wgpBYshuHW HOFZe3HvFeCNcE7A3Znb0qUh5guhx90ycZ/1/+p9XfsiBHBlFOF+xokafF7x AjmfCcKFRQoicwNeQBXVKMZdcb7IT6PSCyhz/YRvD5qiWle94M2XtSW4HY/G TAr9JoMAmlWG+wRX6JoXHGS4LbpC2MxbdEBi1BNyIk3LcRdzCYWW1XlCMmmK 8KMD9Ks24Z5wqki8AreRu7mR5kZPYG7NJOzxbV9ZaK8HeFzcU4mbtVa9xizF A7KvbanCXR88+fK2jAe4+z4m7Pe3sFx6wh3M3shU444I0S12eOUOPNk+hH98 aShAHrjDVZlawlmbx7Ye9HEHEVGPGtz+Jx+6PlVzB2p6HeGF1kCawDZ3GOPZ +47YfyRJOH3ZDf6jLxE+F657R+O9G2j0uNXiviK1bK6d4QYNHJ8JyztQBc6F u0HFvxrvce9fpSZXkN2g7v5zwkMFfAkKJm4gvpLbiDvUZGn56KOz4K/f0EHc /5jF/aSa09CL6vTiHr+0RmNFxQl2J/gR1imWvGCX4gih9iKEUxM5n75e6wjP nz3pwS1s37r3crUD/KM79xk3j6ym+4yfA5AP3iWcWJxaPGfmAF+dNhMud82/ vE/BAaKqurtx+47HdZqKOwBJKYZwUnSr9+11DlC3Zr6LuP8lj7NCFBKMp1wh XHCgXarpAwnODJkQpv5qCch7RgKZWDfCpVu6TBTvkUCZI41wUZ7ckarrJCD7 cBHn1x8wU7T2IkGTwn3CJYKVckt2JEju1CD+b+6N4L5sfRKslf5N2LZbV8zl MLY/OIa4v8cN84b5zSTgf+T4hai3PTI+XdwkqLZW6CPysXFoqXTGHtIGqr/i /jvOCnr+yR4eT/j34+7wCjuX/c4e4oL2fMOtvvBp+9tX9uCrU0t4X8Hfmz1P 7CGQqTeAe7FBI+X3XXsYzRD4jvtx/6NmNy97eCFFIszToHNu1cAeHsUvEesP RWzQyuG1B8VVI8Kxj4MurjTYwernXYSjR52FTSLtwLa6nvCM60KS4kk7CK/a TZzHfbF40VLJDjqOhhKuUlLOfixsB085+gir7S+y4ftlC1vKd/4g6uvcIP3h F1uwsE4eJOr9l9DcMNsGoHOScFe3emhYmQ1US2gPEf0k1VmpINgGOEMohJ8M kqXuCtvA3wD5Edz9NQdv6B0/Bff9bxFu1h5VvMt1CnhL+glP0tYkn/5oDZW7 Q0Zxn4z9T/AtyRoM5jJ+Eu+X/z535rYVXLgjMYZbazq5oUbCCuKPRRK+98WC llRqCR1/Fwlvt+i1um5hCXLdAeO4uco07tlPWcC1qT+E9yTcFA9+aAEBn+Mo uF1uhsibalrA8ZNUwm4DT6Z7RSzAoUuXStRrY/Oveqo5LAjkE47a5CgkXm0O +wt30HD3xJdu74szh3WeyYStxA5OLXiYQ/lPaTpuMR8Zz+tgDndqvQjrpw69 +ClnDuadZYSVfBadzmwxB3a/HQN3Gn0lN49lBvqcrwlfvAVzvv1mYJ+9boJ4 377cfVYNZqC104GwoFHoYeNCM/iknkz4vG7+PocnZvDTqobwgfctp+oum4H5 6CnCInJrhCNNzEB9ve//3CJl9MDEBGyEyEzcup9CJz7JGoOoiygLd+OKC2/c RiPYHH6NcF/eiIQBlyGEbPyf09PEwySDDeCn4mnCoZ+zRaL/6AO3diNhkSvt 71Tf6cPKdo9J3EIOqtM3QvSBQ3ic8G3HcH4PbX2woO9jE/0pOu76JI8+eO95 Nk3EH373r/obXWjs2zxD9Cuo2xmlpwvVZrmEXcneZpt+6ICvLfcs7krZlu6c bB0IkXcnzGvXzW7x0wGT4i+EnRQ2Ui9p6YB9s8ov3JJky5RUQR1IsdMk5ofK hY9H62vDpcZswiYp6O+QKC34P75rlII= "]]}}, 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.662347996957409*^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:eJwl0gk0lVsbB3DkOPM5GslUpEyZznsTqfs+1MWNrhSKKHPK3CVJMlwyx0FR Slw0aTB9KCFDyEwqpCSOoch4SKZvn+5ea69n/dbea69n//eWtvM44ijAx8f3 Ek1e/W+Q4L8qDEXDecrRKhRwI7wv6ZxiQkahYftqGA0Or/q6WQ8wYfivydaH 3cLg1mxvGPeKCSnqZkd9CSLQd07Du/UKE9jpfWofqkTBL8nxHPs8E+pP+oQ5 JogBIeONsKwTE2TPUo5720qB50LxheYDTDAqstHz1NkGq3INkVl8TBjAg4Mr DGVBbsF/r80EA4K9ZZ4+MtsOZ/nrPlA+MaDQpTvvsZ0cJEhZKu8uZcAhuarL JHt5oOXU+5Q9YEDNuAJ221oBGqa/fNJOYcDx58m+UvuUIDjQUFTehwG9lTV6 mlt3gsKjGP4b9gyQe9Jt8G1hJzSKBBwgHWGA1pmQ3VvjVKAsa8c/QyoM8DnV 6SpwUBUq6bHRxyQZwHEWidL6qQoeGiW0RioDbn5r7D99Ww3ObZaJfDFCh2w5 j81VHDVwEKN+0H5Ph8sMVVVRBXUQPDVr/vIVHUqa+hSST6uDy5HDgQcL6bBs O29yLkMd8EdT9T3/0mH44rn4/nfq8Nbu6C0vNh0K394a2UJiwU+nMXNmEB0+ 1totL59kQQJrpdbGmg6KcX1FomEscOt6Qd5oRIdSk4d35+6xoD5Xw7h9Dx3U 2dtKO+pYwNLParqmQIdw3aTrnEEWcLZ65tqJ0sHQsiDJcZUFW/R1ZDWJdDCb ++50RQQD/shmo01zNOiIdim7oIzBVTnrg8uDNLDwYz87q4NB07MyrfE3NBjT JzPDj2IgL/ACG6qigaeGs8ScPQZq++78NZJHgw9Es/becxi86ze8PZNOg6FU yS/WQRioR9xRI8XTwMpCeCElBoPylrebdgTSQJEw55+TjMEbFUurQ+40sCl+ YJWbgfbPP6ZetqaBifv0vpcPMdCazlQuMaJB0oH1/VP5GPg6Fr1e1KaBt31X vcVzDMQGuzkGSjSY8yYSqZUYeLxnhN8Ro8Hvt0cSN9Wh+3Qt5q6SaSA1MmQe 04TB+dRbTmcWqNAQIW1+sR2Dr24ND3pHqBDM3Fo8/hYDb/J7f4suKuxeM02Z 7sag0lnhY18dFeyvhs58+IjB0Y3P9BlFVPBzfy5v1Y/yk4gf18+iwqbysuaA QQx2myo1xyRQodaVXnBwGIOn6ocHe4OowD75rbZxFIMYcUM1TQ8qZGizV0hj GMyKDj1Os6bCGKPYn/wdgy7bAAthIypQSuTS3kygfAYoqrF7qKDku0PLbQqD Yxd1FTcqUMGlyMz4wzQGq0nuevdFqNAV/C5ZZhYDgkHRZT0htB6duUGfi4H1 0mjzxCwFNMvbYozmUH9JZSpZAxTwR6dqzGOgo5qVat9BgV0y6yWIPzCQo3IZ KpUUwF5FtlYiY9m9IWtyKZC5KfmJ8wIG+7admRpIo0DTrdM9y8gBj+JNW2Mp sKfQX/mfnxgMPlK6++oSBWxaLtj+RM7YT/hc60KB1pC/xh0W0fuLtyx2WFKg QX2VVYOsOaoyM/onBRLM2k5uXsIg7v6LCqoWBT4G6kc6LvHyJJzQlKeAXIB6 9gPkHpeocg8RCsRVlMYNIe+KrB3IE6JAS2sFXXIZAy8t5ZoVLhn0IqqbDiFz bylZmHPIMG67z9UPWSZc6VZJJxlYyiI+6chpg+wr22vI4O9cdroambmlau3t AjJQbk7lfkE2fSaxe0smGeq7ZIpXkJUj+CdzEsjgvXqjXWQFgyMerbv2h5Dh fIyaiQryqCKHn+NFhtU59iVd5EfLGYfYtmT47OWaaIo898pxnYEJGbpKTFsd kNuC2H+Sdchg7yFi+zfyVYgY61Qjg7D3D98g5OLox7MPt5KB4HBbOwZ5vYbz iWhhMlw7VdJ6HfmztbTYeT4yCALJPB3ZeN1+RZdJElyrXPv1PnKE1t6ws59J ELay824ucpdRgKx3Gwlyuc0Zxcjbpf34I16SYCwtcrIMOaXm3sa7uSRwTwzJ rkbWlQuxbEknwZ7A/rZ6ZPd/w17zsUlwQYob1Iws3Y9Z7QsmQc2pwOp25P6r TFqoFwl04p9kvkU2cM3o6LQlgfyS0M5uZKsJ7weqR0jwcyPLqRc5/4JIVJIu CQ7vVWB+Qf76UNL40jYSRB9MiRlElit6umN1PQluuLW/H0LuKf/7e7QgCYyi VOdGkGWfUdNluETQ7ysZ+4r8ba34rmoOEap2WleMIT9cd+ie2zsi7Nce9PqO HJlhPCldR4QyLzx0CjlxLHL83n0ipAvGdUwjs/+42XHxJhEcj00wZpGxNlkN 8xgiLBz1Pj6HzJ00c9jhSYQY30zneeRGo25ZcTsiJCqOuv5AbpmMMxY1JULI d4r5T2S/wPOtSppE2ENZ3LuI7PolWUxXkQgOxEnxJeR2l3elNhJE8CnvrVlG 5qc/GyviI4L/nSudfKsY6OkPXKp9LQQPjY/E8iMLOleOS70QgsbeMF0BZO5R sdjAJ0IwWHkgcw3yOvKT/RaJQqAnGfqHEPKARSXR3VoIwuodNajIfLuE8kRn CHAmcKmF58nKoeopDgEelKY70ZCNnFRKOroIQA8wSaQjD19ex3pUTgDV0OcD TGT2b4IjL6IIkBBl5bwBuc960MBElgDbYn+XlEROjXfs2m4tCKrRX2X3Irty I93MkvhBN4BNskCuYi/6ly0s4/M9DGlPZKvfXoV5dS7gu69eZIchZ/dB4YlF Lh7TEx2cjFy94+AZmeppfPTB1g/ZyBF5Ty2CVCdx1hM2VoCsqff0vk3SN1yh OPRTGXJkT8rIk4phnIMZDtYi44n9nzNCBnFXykuDVuT325zU/C7345yQWv53 yEop03WXFD7h/bP3oRdZjKPSuvq4G5cQ3+jdj5zpY+8fPv4W90qWusdBlt1V kPZ4rAMX+PwlehRZu7r9lIV0G55kU3tvDNkkVMQgtaERN+EGRE0gB62tad+c XYfnFfJxp5CL2l/u365cg5sStAdmkO258Zm3nCpxQuJoPhd5050RCf+sMjx9 5ZDUPLIQfdYsP/QZ/m112fQHsgRbzPoUqQhf8agMWEA2VFKoTnudj689YWb9 E9nzy5H1kuq5uEgDF1tElowCo+PxObj4Kf5hnl2u2bxpDr+Hh8TE+C4hJwts d7Q9lIWfZM0u8axL7vJsMsrA540OOC/z/tth3YLQptt4rnt5Fc9OalTmsbIb +I9KI9EV5PSomX/2cK/hBdfLTvKs6dItkPM9Ae+cu5rCs96YsPKGLfG4481d jTy39G/g/1QRjZPWSnJ5FrH8UTc4E44PBP1PYhWZOpdbaukSihf6iezlOdi9 Tm0yORgP4jey4DmSYyeqOBWAy6jke/G878WroaZYP7woVSeMZzUBn9FWynm8 ItYjmef470xp0bpzuM6f89k8E6+/fPdD0RNPDYJ8nlvHmGNnxlzwBeO8Mp79 NlRtjm07jY+KdtXyvObnbIJ4uAMuX7it7pct9SxnXBzw6G/uv5xz7Pj6icMO OEFdsJ7nNNmcVBNxB7xoh+JrnkvbrzSW5dnjUyyHRp4raxOBzLHDQ5kTv1y1 Z3VEqcAOtwu40MRzvtQxrCPYDj8gEd7Mc698u0nWVjt80pTWwrNFPde3edIW V1eJ/2WDBuqgd6Ut3pOzrvWXXTbgNxJs8fL2xF9eGhZmaTnY4v8HPXOO+Q== "]]}}, 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.662347998910616*^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:eJwl0wc41tsfAHBbyHqnoiukjMwfiludbxRSZtKgMl6j12hQRiQzKipb6iYu mRWvJLllVUbJSgllXKPIH5FZ/c/73vM85/k+n+f7nO/Zcs4nbVx5uLi46nBn x/8aBf6LUnD0h6lcTioVpIO3ui19p8OIxur6CCM6rEuZ39M0Soe8R97GVR1r IUplsH+hhQ46NyNGlvbLwYjY8ZapDDpcXs/c//6IPPC6izr9iKfDtWW5ZF6G AvAashq/htFh/+GBjiR/RXi2Ts0y350OV7m9wDhYGbT+lQKWNh2MdWaCOnxU gFfyjLS7Ih0mlH5Ouh9WBdv9gc/oUnQwjIwojSKrgb6KTK7vTxpctiHT+wfU QJqx9Gb9FA0sGoLXDeWowx/U1piWQRqEJG7MEBfUBFHpE1fVG2jgpHTwucJh TeCmN5MGK2nQusT3Yz5HE37JVBekFdPAgRablfdNE3JGL72xyaRBwuv6OntN LXhSMzUqkUiDCBlBAvloQWPZUfuOKBrsNa9+GH1PC3TmdUNvBtBA9vqJfEVR bTDreLBq6zEaZLtb/zmnrw3ZdYkzYtY0eN+sfJfmog1vjeRTx41ooFcCoq9i tKGAsVapRY8Gzn4f7xoWaoNwTXVbuTINqJZmH+qatIGc6Pw0RwabISARMaoN 3YfmBTPEaaArKW5Uzk3ArEFMVwoPDRh9VQ+vryFg1/YLhulzVHjNPPDEUIMA iw8CjKwxKlzPC/68aEiAz5UAq9IeKuSb/1XTb0uAvknK2sYWbL5eIVlXAjYI 1bSP1FBBWE8nrtuXgK57Oy+JPKKCt6ifq3YYAVb8F3dtyaOCAFVdYU8cAaUb tlGZGVQIdDXPNUgjIN0giefveCr8iD1eJZdFQOHxU5ThMCqkNDewaIUE+NvZ WqudpcJFUZcFDRYB2oy6Z8EeVHCtIvOEVxIQK7XdqcOeCl/cip1la/B4H8tt WpZUMG/J01vzioAEVrtNiiEVpi0yBy68JmCpsyWXR48Ku9blhR1rI+A5Ktc/ p0yF8sSSb/XvCDC2c189LYPXP9Wa/U83AUmTSrK+ElRY3F2wY38fAR6zfadX eKnQrBJoe62fABH1Fd64eQp4+LudDBoi4Dp/TOeGcQpcTSiNlxnB53fSfLDu EwX4zldmBI8RYGt8T5XZTgG3NLHNd78S4DfScp/+kgLDWfM7Hk8QwDOeyFNf QYEJo29tSZMEfM/lfRtQSIHYVFoDmiJApnlAQu8vCvQKuUf/M02A0lJN/fJ1 Cshl/+816TsBtX1LI40RFBC6PD1qNEvAkWMVvpnnKJCoKVhkO0dAU8SyW+gJ CogZD0Sa/iCgt+9OlZsDBWLSVsfKzxPQH+Fyys6SAuGZ7t7D2LUjgyEWhhSo VTRgJS4QYG2R32+hSwG1PNMLGosErOgrJB1UooDUkez2SuwKPlaihzQFqIJn XLcs4ftzzvoYJobXv9q+LZftGW1GNjcFBI+I+QsvE6DolKz6ZpYM7608brpg R0jPqP4eJcN4VfphFvaDU+6O+j1keOGkzFrCLnAXqT/fQgY7X5NrBisE8OdP Wr+sIQN6YVbui70YtpZ/zSMyqPWkDuViJ19L/uCbRwaZjSLlndh2E46vujLI MHfMt+4ntufG4Nc7r5FBLH/jFfmfBLiGDg2WhZOh7FR+hBG2gkwCv+Y5Msyu LJOdsG2rwzRYJ8jAo7k7/Tx2T0Sh/Y6jZChuo7snYsfZrYpqsyJD53TOyzxs 3u03cr12kWFqYn75KfbhdPRUfCsZ1rt2ozfY8FP0WaUqGVbODzb0YhOm8/ne smRw/HvpxVfsfaoLgUpkMswshzrOs/M0ceVxATJ0JTQ+5/mF/6/wn6WPlkgQ lSfNJ4rdrn9BImaSBC82K+2jYWsJ9+5wHiSBTpwm6w/sPXcOb93VRYKQXl2r jdg5eb/n1ZpIsLMjUV8Ne3XY+wDZZyR4WDMXRmAHiEw+kirF9RdsDPSx3apt CtbkkkBX6MaZHdgRRjRr+ZsksJUv2WLEnq/VPEc7ngRcr7anmGAPtItnmYWT wOj00z2W2N8NboXdYJIgIG1tiQ12TGPWkepjOP9N4f0B7GqVgoYfNji/ZFF7 CLu3+8NLHRMSDLNEAu2xs8YPWAX9SYKPZGHuY9hC7WaOrzRIsEl56Jgj9sea kWnpDSSIvjRazMAuHnwX3bOaBHzD53LdsI1/JVw05sb7mxKO9MD239c0VjEn CQ9iQhY9sWv506ZZnyRBYURc+Ay2/s9VjjwlkuARzWXriy2Rs5ztkSsJqQN8 SX7YerKTO99lSAKtsmvxHPZZv6zgmmhJIClnnDmP7eSwPWLYQRKMBnbGRGA/ f/Px6iFhSehaLJBIwH47M1vbekICFravZBax97uQnh/pJwa24VLiPdjOz6NO B5oJwy7d04YCv/H7fmqiso0pADMPkt8qYbcHRf6ZacULEk4z/abY5akkI+5T XFDZaF7ign07yZp1UHwFhbiMbA7EXqVkJ9Sps4BiBgvl47DTNMbiUu7Mouor lfTb2Ioy+7NQ7TTyeFKRXIBd3N2ZmFE0ieofWMo/xlbKmoh9o/wVcYUwtWqw y1hFobNxI+jzwaKYJmxyxbkPQ8QQEvEIP9iOPZz3WrnItR/trnuQ+gH7jXrK mWDzXmT+jYv5CZuZVbm+SPYDUnVulRzCXrhr46jY0Ike5u9zHsXO+Ldhi+P9 NiT0RyzvOPbeyInBus0tSFz1yuA37ENbPfbKCjShL6ThG1PYA6fUTS+qvES0 cKGqGWx760oZHmod2l9CjZnF3uH/+OzzwGpESkdbfmDrn+HfJL6mCs0U7k6Y x/5Cn4iELRVo2CumYwE7Xt1MJPlWGVpnxde3iN23ab1tQFoJsm869HgJO+rv MvKSw33U4hoVuoxdGHP3ReTnfHTewlN2BVtr86v+k6a5KGTV5Vy25/vnKh1Q Npr3GdvwE9vZWe24j2cmoqA9sWzfSVQY8Le6hQ6avhthm9di7d49Mmnouke9 1i/sDSjp0PrqJPT7tvBptuWCmE0WX26glzt4C9k23+hY+j42HpV5N35k2zSq +38TlpcRWsgS+I3d3Vhs0eMSjVyCu5TYHjPd+alLOgLZpDNN2VYem/FIO38R 3Rs8yGDbr+9AhSU1GBm1uASyrauhukJRD0AmVoVX2Db4bXr15IAfKvIXv822 6cX44IGW0+jAN6KA7dZnR8wPsHyQmOBVjhn3n7Q9vuyDrjV/5rjP5JKwlbMP 8juuUci2dIl31ToDH6RyJ4jjdIedd9eTfFB5IG8R249M147W13qjG8KBHH/N G84y8PNGUtwfOf6+xj8jdKM34kvWLWa7S27x5o1uL1QqxeA4ZExn5OxLL2Rh nMuxUrN/txvLC+1BXPfZPl7D8mJkeiEUdYrjbWb5qT5xXujftyMc67SlWEcF eSFa59EHbF89IZ2Q5+6FKkO2PWTbXjA6imLohTYJ/8Xx4NfWTUc1cP3tixxL aOZplMp4IWdd8xK2ecpbk8jCXkivN41js9uilmHznkjiJLWU7fgd43fD2j2R cUc4x6HTwaOlxZ7oF980xyDx/hMl1hO1tzJYbL8VKdx8j+GJ+E0sy9gmKwdR HWQ80Zc1NRwnqssoXVhgIsdr6o/YrrE58extJxNl9f/mWPJX5FaueCbSuMUo Z/twjV7zJyYT2fm+4linlzY3bsJE/xgrP2abmantuEmRie7PXOa4J2mh4joP E/0foB6bHQ== "]]}}, 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.662348000895073*^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:eJw11Hc8Vu/7AHCrKCFKKYq0kJUjLVzH6GNmy3iMT8p6VAqFSqRkpcJHQoOQ WUJGVsrOikhSVjzL83iMrFC/c8739Tv/3K/36x7Xue/rel27XLwtXLk4ODgU ODk48PF/3xT8//i5S0nraeY0WApcK+GeZcPnD1WGEn6z4F15aI0OhQ2uovmO oxfmwE/EXPpVPxsCZy9uupm0ANt+N44h7WzoldILz1tYgh2n1RM/1LJBuF07 WClmBXQWrrfrv2FD+hHNzGE1TtQnItVgJZkNKbHqLrrDXGjHy6puwXtssGF/ WVrN50HnrxrIjIWwIUAx0s25jA9NTzj9ttqVDcpcrTtFF9ej7p9IBpU2bGDk 5FnvchFAfVsfV543wNYfVjUusBdGhVciZhcOsOFPkJVsfKcIarD9cE+vBBvW fvJRzrXbjMY37VDyE2AD/0md64Z5W1FvVf17PJOTYPhAcX/GHkm0JPJUwXz+ JPDcGn4w0LIfvTttwf187yTM3FTe8W+zDOrv2Sh1ZfMkvK9p4n5eK4t+OCN1 VId7EjLuZIZA8gF0IFNhuWOYBZMvRZvGAuVRQeq1+AedLLiALAmMmiqgFVNr ec1rWNBwlDy0+4ciCgff0LtTWDCaoudAiVdCVdaf0HsYxQJPcbudy5LKaK62 0RmHQBZUrHfcXXhRGbWaDVje68GCsry5Y3JVymhG6YG3s6dYcEk65k0Y10F0 ONn2UOMJFnTGfZZg6hxE/Vzkl56osuCjRqZ4avlB1MpSCXEWYcEZiWU/B+ZB tE+TFmjAyQKTDX/m/oiroAIb9uQem2LCQMBr4QY9FXRn5l0RZIgJD52aIpje KmhTzcdFlQ4mSFZ9vVP0nwp6Tls75Vg1E3YNJglalqqgTgefHDTMZ0KScpYK X68Kqie9wDydwoQZs6D+ddMqqLHZ7qWbUUyYlTb2jVmPoH27+IJyA5lwU4w/ tGEXgl5WbU4e8GDCFTfuhhE1BI159fqiqC0T5CsNlIUMEVSyLmqbrR4TXnrV H/cnYfPqJvnpakwIyXpsoeWFoEWW3zQW9zLB3a1pb2wggr75eH/IWpQJXZa1 9yLuIOiGy6ppVTxMUDmc6XIoHkGP9rpbNc9MwAzCik14iqD+WlODg8MT0JHw K64lG0G7nWONOTonwCNDMJZdiKAuDZey5KsngCnq4iJdgaAfh/0WzuRNgJYk Z93V91i8WUmdF0kT4KIn1iLYjKD83Yb3Z8MnAEUqHsy0I6h0fuKQ0ZUJ6FJu Uz3yGUG1b19WKzg7ARecNefZfQjaHAsJOywnIPGr2mOZ7wgadCFw5aHWBDy6 rabJHkLQY70lXuLKExC65vwP059Y/LXskdydE0C+F8dlTUHQ1BwrxxMCE6DX K/svHx1BOWUyBxjLDDigHu3jM4GghdardikMBhjtP78uhYWgyk67ek/1M8CM o4ASwUbQgXcnDXY0M8B3UstVexpBi29FF02WMuB6oTqjbQZBw8bJG1ozGdBI okQo/EJQEy5ri6L/GKCQHf7MYw5BrYMdgjNuMUBEhK4cOo+g2RKe99J8GDC6 NrQ1aAFB27QUruacZoAayTTRaRFBxWPlNarMGPC8pyJg/xKCPhxq6+oHBvBl zTO/Yf7+4R8VLiUGJJpIGV77jaA/kQ+2h3Yy4NX591Hrl7H8j6LgK8AA/q/x d6Mwt1eUD1Wu0GGTPT91BXM2jVtViEmHC6+iNp1ZQVDf2alD5wfowOv6wKEG c8x+jaHej3RwulTCEFpFUPu0JFmDCjqgHpIfbTHD/EvBphw6nJeRHEvCnETe EmGaRIdKw8vCPZgzbcj3RyPoMErOXM/3B0EP/D28MySADi+bI4LUMB9mrleU 8aDDmKbnqjNm0n2P8m82dBBs19schtnNr6c0UY8Oe8xpR15gDlvukHY+TIdl 0dW+Osyz4n1zyvvp8HVzo/kgZq86sqzAVjo8LE9fnsP8Q4e3bHYtHZJ67u3m /4ugAQJyj8bmafDWrWrDTsz3dpo3DlFoIFJ+YV4Rs+jv38fGv9BA5o+cqCZm Z/ls9q9GGgjFTCYbYW7v6f8hWEaDcoWuYhvM48c3/VXJokFkWmSYC+blom/G pxNpcHLgruo5zLwSZbVJ4TSoDl4z5If52SnE6rs/DTj9KlOv4679ziHrQQM/ 10vxtzBrbNeuDbalgTOi2haJuersQvSwPg28fp51vo/53VSso9FRGgR/rnP6 D/Nfcq5irSwNckvj+x9h/uUyuKi5nQZx958PPsZMN2wraVpPA1Ki1NVUzD2r W53tl6nw+umJsnTMlskek/MTVJCTiUp9gVmp0s3xyXcqlB5R2JyD+czkCzHz DiqEUMuM8zA/+IfPVaCWChwXb5q9xHxF3VOju5AKR9UC9xRgjt326vGzdCoY W+i1v8ZsyKgOupxAhXXad0yLMAfTgwcsw6nguGv1VTFmrpKh18cCqVBd9mvy DWaBrPZFWS8q8L8uFSnFnGFz7NUuRyr86ym3owyz2Bextl2mVHjsnbOxHPPQ GZeTclpU+B4Swcbtd1hY4ThCBYMby5VvMYcES7lZ7aXCEqvMvwLz84zY2ctb qTCicbeuEvNaNGeha5kCOT4xp6owR9U7nN0wSQEhx8hB3PNXHgibDlMgvemF QzVm9pD576RuChh40rpwiz57vY1VT4FvT8ShBvOen3Vk/TIKJHmtz8T9zSCD lpdDAROSqtM7zCWlyjbR9yhA6sgtwi1ZsE2T9yYFDvFactVilrU21Iz2pYDE MGcC7hGNrf55thRInGF9xe3NXvNCz4gCx63zt73H/Gdy7NuEBgXqk+vicJO2 H9Ux3k2B3XUlrbj9rbac59tCgVqNWK4P+PyJisg2PgpU6jWRcdNVjSPOTY6D /Zsh2TrMrs0jovwV40DuvHC1HvOLM042NRbj8OxtYRruUJ8zAQsnxsGuNKIZ 9znQ0zpydBza57eKNmA+P2ZZ/klyHJRP7M3FnR2rWuXAGgPZ+xeaGzGbmynt /RE5Bt87Cxm4pZ3eO7UFjUGCZdeGJswfXna++nBpDILrc81w83S036y3G4MP RylfcF++lBIsITcGLvKOg82Y2zpWtLa3/YS+KyPfP+L1/M7imbLoT8hbffu0 AzN/NvufyYoROETLovVi3mEqPRR9fBDk4gsrBvB8uTgW7RH4Br4v26RHMLfQ A/1TDnyBTqOqrHHMyu5GYMbshpm++hN0PJ8rPSNrLD/BUgZnOhN/nzaOacqn Vqheem7JxuspMy2mbbYJ+BtKkqYxJ1AVOtIL6qF1OLtsFq/n7T9sv2e8BzfD D+Q5vB+Ff4/N86qB6beklXnMU7nBQbaUt5C3m3xyEbOuQWJoxo5SyORzJS3h /cX4XTCypRhibGJP/cb7QcJV7hrd1/CtXkdyGfPetuD+BN58YIal9OM+H+qx f/lrFniLlV9ewdx3ePgk2S4Ddv3eN4P7murt7a1308BTi9ttFfNHhWm5O+NP QCy3pR43kkp2el2VBII0buE/mHevWkmUHXkIi11i9rgjTbhl+mfiQE2yPh73 NgfmkUyrB9DWIPwR98sMdf+etmgo5IleINy12HipMhxawtJ2/MWsmI62TgXe hrAb3pq4SU6OlQ0+N+GTlJAtbntZZ3mWwA3Iat/gjVvm4WJ831wgpBYshuHW HOFZe3HvFeCNcE7A3Znb0qUh5guhx90ycZ/1/+p9XfsiBHBlFOF+xokafF7x AjmfCcKFRQoicwNeQBXVKMZdcb7IT6PSCyhz/YRvD5qiWle94M2XtSW4HY/G TAr9JoMAmlWG+wRX6JoXHGS4LbpC2MxbdEBi1BNyIk3LcRdzCYWW1XlCMmmK 8KMD9Ks24Z5wqki8AreRu7mR5kZPYG7NJOzxbV9ZaK8HeFzcU4mbtVa9xizF A7KvbanCXR88+fK2jAe4+z4m7Pe3sFx6wh3M3shU444I0S12eOUOPNk+hH98 aShAHrjDVZlawlmbx7Ye9HEHEVGPGtz+Jx+6PlVzB2p6HeGF1kCawDZ3GOPZ +47YfyRJOH3ZDf6jLxE+F657R+O9G2j0uNXiviK1bK6d4QYNHJ8JyztQBc6F u0HFvxrvce9fpSZXkN2g7v5zwkMFfAkKJm4gvpLbiDvUZGn56KOz4K/f0EHc /5jF/aSa09CL6vTiHr+0RmNFxQl2J/gR1imWvGCX4gih9iKEUxM5n75e6wjP nz3pwS1s37r3crUD/KM79xk3j6ym+4yfA5AP3iWcWJxaPGfmAF+dNhMud82/ vE/BAaKqurtx+47HdZqKOwBJKYZwUnSr9+11DlC3Zr6LuP8lj7NCFBKMp1wh XHCgXarpAwnODJkQpv5qCch7RgKZWDfCpVu6TBTvkUCZI41wUZ7ckarrJCD7 cBHn1x8wU7T2IkGTwn3CJYKVckt2JEju1CD+b+6N4L5sfRKslf5N2LZbV8zl MLY/OIa4v8cN84b5zSTgf+T4hai3PTI+XdwkqLZW6CPysXFoqXTGHtIGqr/i /jvOCnr+yR4eT/j34+7wCjuX/c4e4oL2fMOtvvBp+9tX9uCrU0t4X8Hfmz1P 7CGQqTeAe7FBI+X3XXsYzRD4jvtx/6NmNy97eCFFIszToHNu1cAeHsUvEesP RWzQyuG1B8VVI8Kxj4MurjTYwernXYSjR52FTSLtwLa6nvCM60KS4kk7CK/a TZzHfbF40VLJDjqOhhKuUlLOfixsB085+gir7S+y4ftlC1vKd/4g6uvcIP3h F1uwsE4eJOr9l9DcMNsGoHOScFe3emhYmQ1US2gPEf0k1VmpINgGOEMohJ8M kqXuCtvA3wD5Edz9NQdv6B0/Bff9bxFu1h5VvMt1CnhL+glP0tYkn/5oDZW7 Q0Zxn4z9T/AtyRoM5jJ+Eu+X/z535rYVXLgjMYZbazq5oUbCCuKPRRK+98WC llRqCR1/Fwlvt+i1um5hCXLdAeO4uco07tlPWcC1qT+E9yTcFA9+aAEBn+Mo uF1uhsibalrA8ZNUwm4DT6Z7RSzAoUuXStRrY/Oveqo5LAjkE47a5CgkXm0O +wt30HD3xJdu74szh3WeyYStxA5OLXiYQ/lPaTpuMR8Zz+tgDndqvQjrpw69 +ClnDuadZYSVfBadzmwxB3a/HQN3Gn0lN49lBvqcrwlfvAVzvv1mYJ+9boJ4 377cfVYNZqC104GwoFHoYeNCM/iknkz4vG7+PocnZvDTqobwgfctp+oum4H5 6CnCInJrhCNNzEB9ve//3CJl9MDEBGyEyEzcup9CJz7JGoOoiygLd+OKC2/c RiPYHH6NcF/eiIQBlyGEbPyf09PEwySDDeCn4mnCoZ+zRaL/6AO3diNhkSvt 71Tf6cPKdo9J3EIOqtM3QvSBQ3ic8G3HcH4PbX2woO9jE/0pOu76JI8+eO95 Nk3EH373r/obXWjs2zxD9Cuo2xmlpwvVZrmEXcneZpt+6ICvLfcs7krZlu6c bB0IkXcnzGvXzW7x0wGT4i+EnRQ2Ui9p6YB9s8ov3JJky5RUQR1IsdMk5ofK hY9H62vDpcZswiYp6O+QKC34P75rlII= "]]}}, {{}, {}, {RGBColor[1, 0, 0], Dashing[{Small, Small}], LineBox[CompressedData[" 1:eJwl0gk0lVsbB3DkOPM5GslUpEyZznsTqfs+1MWNrhSKKHPK3CVJMlwyx0FR Slw0aTB9KCFDyEwqpCSOoch4SKZvn+5ea69n/dbea69n//eWtvM44ijAx8f3 Ek1e/W+Q4L8qDEXDecrRKhRwI7wv6ZxiQkahYftqGA0Or/q6WQ8wYfivydaH 3cLg1mxvGPeKCSnqZkd9CSLQd07Du/UKE9jpfWofqkTBL8nxHPs8E+pP+oQ5 JogBIeONsKwTE2TPUo5720qB50LxheYDTDAqstHz1NkGq3INkVl8TBjAg4Mr DGVBbsF/r80EA4K9ZZ4+MtsOZ/nrPlA+MaDQpTvvsZ0cJEhZKu8uZcAhuarL JHt5oOXU+5Q9YEDNuAJ221oBGqa/fNJOYcDx58m+UvuUIDjQUFTehwG9lTV6 mlt3gsKjGP4b9gyQe9Jt8G1hJzSKBBwgHWGA1pmQ3VvjVKAsa8c/QyoM8DnV 6SpwUBUq6bHRxyQZwHEWidL6qQoeGiW0RioDbn5r7D99Ww3ObZaJfDFCh2w5 j81VHDVwEKN+0H5Ph8sMVVVRBXUQPDVr/vIVHUqa+hSST6uDy5HDgQcL6bBs O29yLkMd8EdT9T3/0mH44rn4/nfq8Nbu6C0vNh0K394a2UJiwU+nMXNmEB0+ 1totL59kQQJrpdbGmg6KcX1FomEscOt6Qd5oRIdSk4d35+6xoD5Xw7h9Dx3U 2dtKO+pYwNLParqmQIdw3aTrnEEWcLZ65tqJ0sHQsiDJcZUFW/R1ZDWJdDCb ++50RQQD/shmo01zNOiIdim7oIzBVTnrg8uDNLDwYz87q4NB07MyrfE3NBjT JzPDj2IgL/ACG6qigaeGs8ScPQZq++78NZJHgw9Es/becxi86ze8PZNOg6FU yS/WQRioR9xRI8XTwMpCeCElBoPylrebdgTSQJEw55+TjMEbFUurQ+40sCl+ YJWbgfbPP6ZetqaBifv0vpcPMdCazlQuMaJB0oH1/VP5GPg6Fr1e1KaBt31X vcVzDMQGuzkGSjSY8yYSqZUYeLxnhN8Ro8Hvt0cSN9Wh+3Qt5q6SaSA1MmQe 04TB+dRbTmcWqNAQIW1+sR2Dr24ND3pHqBDM3Fo8/hYDb/J7f4suKuxeM02Z 7sag0lnhY18dFeyvhs58+IjB0Y3P9BlFVPBzfy5v1Y/yk4gf18+iwqbysuaA QQx2myo1xyRQodaVXnBwGIOn6ocHe4OowD75rbZxFIMYcUM1TQ8qZGizV0hj GMyKDj1Os6bCGKPYn/wdgy7bAAthIypQSuTS3kygfAYoqrF7qKDku0PLbQqD Yxd1FTcqUMGlyMz4wzQGq0nuevdFqNAV/C5ZZhYDgkHRZT0htB6duUGfi4H1 0mjzxCwFNMvbYozmUH9JZSpZAxTwR6dqzGOgo5qVat9BgV0y6yWIPzCQo3IZ KpUUwF5FtlYiY9m9IWtyKZC5KfmJ8wIG+7admRpIo0DTrdM9y8gBj+JNW2Mp sKfQX/mfnxgMPlK6++oSBWxaLtj+RM7YT/hc60KB1pC/xh0W0fuLtyx2WFKg QX2VVYOsOaoyM/onBRLM2k5uXsIg7v6LCqoWBT4G6kc6LvHyJJzQlKeAXIB6 9gPkHpeocg8RCsRVlMYNIe+KrB3IE6JAS2sFXXIZAy8t5ZoVLhn0IqqbDiFz bylZmHPIMG67z9UPWSZc6VZJJxlYyiI+6chpg+wr22vI4O9cdroambmlau3t AjJQbk7lfkE2fSaxe0smGeq7ZIpXkJUj+CdzEsjgvXqjXWQFgyMerbv2h5Dh fIyaiQryqCKHn+NFhtU59iVd5EfLGYfYtmT47OWaaIo898pxnYEJGbpKTFsd kNuC2H+Sdchg7yFi+zfyVYgY61Qjg7D3D98g5OLox7MPt5KB4HBbOwZ5vYbz iWhhMlw7VdJ6HfmztbTYeT4yCALJPB3ZeN1+RZdJElyrXPv1PnKE1t6ws59J ELay824ucpdRgKx3Gwlyuc0Zxcjbpf34I16SYCwtcrIMOaXm3sa7uSRwTwzJ rkbWlQuxbEknwZ7A/rZ6ZPd/w17zsUlwQYob1Iws3Y9Z7QsmQc2pwOp25P6r TFqoFwl04p9kvkU2cM3o6LQlgfyS0M5uZKsJ7weqR0jwcyPLqRc5/4JIVJIu CQ7vVWB+Qf76UNL40jYSRB9MiRlElit6umN1PQluuLW/H0LuKf/7e7QgCYyi VOdGkGWfUdNluETQ7ysZ+4r8ba34rmoOEap2WleMIT9cd+ie2zsi7Nce9PqO HJlhPCldR4QyLzx0CjlxLHL83n0ipAvGdUwjs/+42XHxJhEcj00wZpGxNlkN 8xgiLBz1Pj6HzJ00c9jhSYQY30zneeRGo25ZcTsiJCqOuv5AbpmMMxY1JULI d4r5T2S/wPOtSppE2ENZ3LuI7PolWUxXkQgOxEnxJeR2l3elNhJE8CnvrVlG 5qc/GyviI4L/nSudfKsY6OkPXKp9LQQPjY/E8iMLOleOS70QgsbeMF0BZO5R sdjAJ0IwWHkgcw3yOvKT/RaJQqAnGfqHEPKARSXR3VoIwuodNajIfLuE8kRn CHAmcKmF58nKoeopDgEelKY70ZCNnFRKOroIQA8wSaQjD19ex3pUTgDV0OcD TGT2b4IjL6IIkBBl5bwBuc960MBElgDbYn+XlEROjXfs2m4tCKrRX2X3Irty I93MkvhBN4BNskCuYi/6ly0s4/M9DGlPZKvfXoV5dS7gu69eZIchZ/dB4YlF Lh7TEx2cjFy94+AZmeppfPTB1g/ZyBF5Ty2CVCdx1hM2VoCsqff0vk3SN1yh OPRTGXJkT8rIk4phnIMZDtYi44n9nzNCBnFXykuDVuT325zU/C7345yQWv53 yEop03WXFD7h/bP3oRdZjKPSuvq4G5cQ3+jdj5zpY+8fPv4W90qWusdBlt1V kPZ4rAMX+PwlehRZu7r9lIV0G55kU3tvDNkkVMQgtaERN+EGRE0gB62tad+c XYfnFfJxp5CL2l/u365cg5sStAdmkO258Zm3nCpxQuJoPhd5050RCf+sMjx9 5ZDUPLIQfdYsP/QZ/m112fQHsgRbzPoUqQhf8agMWEA2VFKoTnudj689YWb9 E9nzy5H1kuq5uEgDF1tElowCo+PxObj4Kf5hnl2u2bxpDr+Hh8TE+C4hJwts d7Q9lIWfZM0u8axL7vJsMsrA540OOC/z/tth3YLQptt4rnt5Fc9OalTmsbIb +I9KI9EV5PSomX/2cK/hBdfLTvKs6dItkPM9Ae+cu5rCs96YsPKGLfG4481d jTy39G/g/1QRjZPWSnJ5FrH8UTc4E44PBP1PYhWZOpdbaukSihf6iezlOdi9 Tm0yORgP4jey4DmSYyeqOBWAy6jke/G878WroaZYP7woVSeMZzUBn9FWynm8 ItYjmef470xp0bpzuM6f89k8E6+/fPdD0RNPDYJ8nlvHmGNnxlzwBeO8Mp79 NlRtjm07jY+KdtXyvObnbIJ4uAMuX7it7pct9SxnXBzw6G/uv5xz7Pj6icMO OEFdsJ7nNNmcVBNxB7xoh+JrnkvbrzSW5dnjUyyHRp4raxOBzLHDQ5kTv1y1 Z3VEqcAOtwu40MRzvtQxrCPYDj8gEd7Mc698u0nWVjt80pTWwrNFPde3edIW V1eJ/2WDBuqgd6Ut3pOzrvWXXTbgNxJs8fL2xF9eGhZmaTnY4v8HPXOO+Q== "]]}}, {{}, {}, {RGBColor[0, 1, 0], LineBox[CompressedData[" 1:eJwl0wc41tsfAHBbyHqnoiukjMwfiludbxRSZtKgMl6j12hQRiQzKipb6iYu mRWvJLllVUbJSgllXKPIH5FZ/c/73vM85/k+n+f7nO/Zcs4nbVx5uLi46nBn x/8aBf6LUnD0h6lcTioVpIO3ui19p8OIxur6CCM6rEuZ39M0Soe8R97GVR1r IUplsH+hhQ46NyNGlvbLwYjY8ZapDDpcXs/c//6IPPC6izr9iKfDtWW5ZF6G AvAashq/htFh/+GBjiR/RXi2Ts0y350OV7m9wDhYGbT+lQKWNh2MdWaCOnxU gFfyjLS7Ih0mlH5Ouh9WBdv9gc/oUnQwjIwojSKrgb6KTK7vTxpctiHT+wfU QJqx9Gb9FA0sGoLXDeWowx/U1piWQRqEJG7MEBfUBFHpE1fVG2jgpHTwucJh TeCmN5MGK2nQusT3Yz5HE37JVBekFdPAgRablfdNE3JGL72xyaRBwuv6OntN LXhSMzUqkUiDCBlBAvloQWPZUfuOKBrsNa9+GH1PC3TmdUNvBtBA9vqJfEVR bTDreLBq6zEaZLtb/zmnrw3ZdYkzYtY0eN+sfJfmog1vjeRTx41ooFcCoq9i tKGAsVapRY8Gzn4f7xoWaoNwTXVbuTINqJZmH+qatIGc6Pw0RwabISARMaoN 3YfmBTPEaaArKW5Uzk3ArEFMVwoPDRh9VQ+vryFg1/YLhulzVHjNPPDEUIMA iw8CjKwxKlzPC/68aEiAz5UAq9IeKuSb/1XTb0uAvknK2sYWbL5eIVlXAjYI 1bSP1FBBWE8nrtuXgK57Oy+JPKKCt6ifq3YYAVb8F3dtyaOCAFVdYU8cAaUb tlGZGVQIdDXPNUgjIN0giefveCr8iD1eJZdFQOHxU5ThMCqkNDewaIUE+NvZ WqudpcJFUZcFDRYB2oy6Z8EeVHCtIvOEVxIQK7XdqcOeCl/cip1la/B4H8tt WpZUMG/J01vzioAEVrtNiiEVpi0yBy68JmCpsyWXR48Ku9blhR1rI+A5Ktc/ p0yF8sSSb/XvCDC2c189LYPXP9Wa/U83AUmTSrK+ElRY3F2wY38fAR6zfadX eKnQrBJoe62fABH1Fd64eQp4+LudDBoi4Dp/TOeGcQpcTSiNlxnB53fSfLDu EwX4zldmBI8RYGt8T5XZTgG3NLHNd78S4DfScp/+kgLDWfM7Hk8QwDOeyFNf QYEJo29tSZMEfM/lfRtQSIHYVFoDmiJApnlAQu8vCvQKuUf/M02A0lJN/fJ1 Cshl/+816TsBtX1LI40RFBC6PD1qNEvAkWMVvpnnKJCoKVhkO0dAU8SyW+gJ CogZD0Sa/iCgt+9OlZsDBWLSVsfKzxPQH+Fyys6SAuGZ7t7D2LUjgyEWhhSo VTRgJS4QYG2R32+hSwG1PNMLGosErOgrJB1UooDUkez2SuwKPlaihzQFqIJn XLcs4ftzzvoYJobXv9q+LZftGW1GNjcFBI+I+QsvE6DolKz6ZpYM7608brpg R0jPqP4eJcN4VfphFvaDU+6O+j1keOGkzFrCLnAXqT/fQgY7X5NrBisE8OdP Wr+sIQN6YVbui70YtpZ/zSMyqPWkDuViJ19L/uCbRwaZjSLlndh2E46vujLI MHfMt+4ntufG4Nc7r5FBLH/jFfmfBLiGDg2WhZOh7FR+hBG2gkwCv+Y5Msyu LJOdsG2rwzRYJ8jAo7k7/Tx2T0Sh/Y6jZChuo7snYsfZrYpqsyJD53TOyzxs 3u03cr12kWFqYn75KfbhdPRUfCsZ1rt2ozfY8FP0WaUqGVbODzb0YhOm8/ne smRw/HvpxVfsfaoLgUpkMswshzrOs/M0ceVxATJ0JTQ+5/mF/6/wn6WPlkgQ lSfNJ4rdrn9BImaSBC82K+2jYWsJ9+5wHiSBTpwm6w/sPXcOb93VRYKQXl2r jdg5eb/n1ZpIsLMjUV8Ne3XY+wDZZyR4WDMXRmAHiEw+kirF9RdsDPSx3apt CtbkkkBX6MaZHdgRRjRr+ZsksJUv2WLEnq/VPEc7ngRcr7anmGAPtItnmYWT wOj00z2W2N8NboXdYJIgIG1tiQ12TGPWkepjOP9N4f0B7GqVgoYfNji/ZFF7 CLu3+8NLHRMSDLNEAu2xs8YPWAX9SYKPZGHuY9hC7WaOrzRIsEl56Jgj9sea kWnpDSSIvjRazMAuHnwX3bOaBHzD53LdsI1/JVw05sb7mxKO9MD239c0VjEn CQ9iQhY9sWv506ZZnyRBYURc+Ay2/s9VjjwlkuARzWXriy2Rs5ztkSsJqQN8 SX7YerKTO99lSAKtsmvxHPZZv6zgmmhJIClnnDmP7eSwPWLYQRKMBnbGRGA/ f/Px6iFhSehaLJBIwH47M1vbekICFravZBax97uQnh/pJwa24VLiPdjOz6NO B5oJwy7d04YCv/H7fmqiso0pADMPkt8qYbcHRf6ZacULEk4z/abY5akkI+5T XFDZaF7ign07yZp1UHwFhbiMbA7EXqVkJ9Sps4BiBgvl47DTNMbiUu7Mouor lfTb2Ioy+7NQ7TTyeFKRXIBd3N2ZmFE0ieofWMo/xlbKmoh9o/wVcYUwtWqw y1hFobNxI+jzwaKYJmxyxbkPQ8QQEvEIP9iOPZz3WrnItR/trnuQ+gH7jXrK mWDzXmT+jYv5CZuZVbm+SPYDUnVulRzCXrhr46jY0Ike5u9zHsXO+Ldhi+P9 NiT0RyzvOPbeyInBus0tSFz1yuA37ENbPfbKCjShL6ThG1PYA6fUTS+qvES0 cKGqGWx760oZHmod2l9CjZnF3uH/+OzzwGpESkdbfmDrn+HfJL6mCs0U7k6Y x/5Cn4iELRVo2CumYwE7Xt1MJPlWGVpnxde3iN23ab1tQFoJsm869HgJO+rv MvKSw33U4hoVuoxdGHP3ReTnfHTewlN2BVtr86v+k6a5KGTV5Vy25/vnKh1Q Npr3GdvwE9vZWe24j2cmoqA9sWzfSVQY8Le6hQ6avhthm9di7d49Mmnouke9 1i/sDSjp0PrqJPT7tvBptuWCmE0WX26glzt4C9k23+hY+j42HpV5N35k2zSq +38TlpcRWsgS+I3d3Vhs0eMSjVyCu5TYHjPd+alLOgLZpDNN2VYem/FIO38R 3Rs8yGDbr+9AhSU1GBm1uASyrauhukJRD0AmVoVX2Db4bXr15IAfKvIXv822 6cX44IGW0+jAN6KA7dZnR8wPsHyQmOBVjhn3n7Q9vuyDrjV/5rjP5JKwlbMP 8juuUci2dIl31ToDH6RyJ4jjdIedd9eTfFB5IG8R249M147W13qjG8KBHH/N G84y8PNGUtwfOf6+xj8jdKM34kvWLWa7S27x5o1uL1QqxeA4ZExn5OxLL2Rh nMuxUrN/txvLC+1BXPfZPl7D8mJkeiEUdYrjbWb5qT5xXujftyMc67SlWEcF eSFa59EHbF89IZ2Q5+6FKkO2PWTbXjA6imLohTYJ/8Xx4NfWTUc1cP3tixxL aOZplMp4IWdd8xK2ecpbk8jCXkivN41js9uilmHznkjiJLWU7fgd43fD2j2R cUc4x6HTwaOlxZ7oF980xyDx/hMl1hO1tzJYbL8VKdx8j+GJ+E0sy9gmKwdR HWQ80Zc1NRwnqssoXVhgIsdr6o/YrrE58extJxNl9f/mWPJX5FaueCbSuMUo Z/twjV7zJyYT2fm+4linlzY3bsJE/xgrP2abmantuEmRie7PXOa4J2mh4joP E/0foB6bHQ== "]]}}}, 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.6623480032233157`*^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, 1970, 59, 138, "Input"], Cell[2552, 83, 421, 12, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[3010, 100, 1103, 33, 52, "Input"], Cell[4116, 135, 2120, 65, 98, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[6273, 205, 184, 4, 31, "Input"], Cell[6460, 211, 368, 10, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[6865, 226, 189, 3, 31, "Input"], Cell[7057, 231, 215, 5, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[7309, 241, 319, 6, 31, "Input"], Cell[7631, 249, 396, 12, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[8064, 266, 184, 3, 31, "Input"], Cell[8251, 271, 253, 7, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[8541, 283, 1922, 54, 52, "Input"], Cell[10466, 339, 6151, 107, 71, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[16654, 451, 1826, 53, 52, "Input"], Cell[18483, 506, 4276, 77, 121, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[22796, 588, 1793, 51, 52, "Input"], Cell[24592, 641, 4259, 77, 154, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[28888, 723, 167, 3, 31, "Input"], Cell[29058, 728, 13912, 238, 186, "Output"] }, Open ]] } ] *) (* End of internal cache information *)