(* 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[ 38136, 916] NotebookOptionsPosition[ 36800, 864] NotebookOutlinePosition[ 37144, 879] CellTagsIndexPosition[ 37101, 876] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"f", "[", RowBox[{"n_", ",", "p_"}], "]"}], ":=", RowBox[{ RowBox[{"4", "n", RowBox[{"(", RowBox[{"10", "-", "n"}], ")"}]}], "-", FractionBox[ RowBox[{"5", "n", "*", "p"}], RowBox[{"3", "+", "n", "+", RowBox[{"1.5", SuperscriptBox["n", "2"]}]}]]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"g", "[", RowBox[{"n_", ",", "p_"}], "]"}], ":=", RowBox[{ FractionBox[ RowBox[{"12", "n", "*", "p"}], RowBox[{"1", "+", RowBox[{"2", "n"}], "+", SuperscriptBox["n", "2"]}]], "-", RowBox[{"7", "p"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Solve", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"4", "n", RowBox[{"(", RowBox[{"10", "-", "n"}], ")"}]}], "-", FractionBox[ RowBox[{"5", "n", "*", "p"}], RowBox[{"3", "+", "n", "+", RowBox[{"1.5", SuperscriptBox["n", "2"]}]}]]}], "\[Equal]", "0"}], "&&", RowBox[{ RowBox[{ FractionBox[ RowBox[{"12", "n", "*", "p"}], RowBox[{"1", "+", RowBox[{"2", "n"}], "+", SuperscriptBox["n", "2"]}]], "-", RowBox[{"7", "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}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"n", "\[Rule]", "0"}], ",", RowBox[{"p", "\[Rule]", "0"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"n", "\[Rule]", "10.000000000000002`"}], ",", RowBox[{"p", "\[Rule]", RowBox[{"-", "8.526512829121203`*^-14"}]}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.6623091301711726`*^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[{"4", " ", RowBox[{"(", RowBox[{"10", "-", "n"}], ")"}]}], "-", RowBox[{"4", " ", "n"}], "+", FractionBox[ RowBox[{"5", " ", "n", " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"3.`", " ", "n"}]}], ")"}], " ", "p"}], SuperscriptBox[ RowBox[{"(", RowBox[{"3", "+", "n", "+", RowBox[{"1.5`", " ", SuperscriptBox["n", "2"]}]}], ")"}], "2"]], "-", FractionBox[ RowBox[{"5", " ", "p"}], RowBox[{"3", "+", "n", "+", RowBox[{"1.5`", " ", SuperscriptBox["n", "2"]}]}]]}], RowBox[{"-", FractionBox[ RowBox[{"5", " ", "n"}], RowBox[{"3", "+", "n", "+", RowBox[{"1.5`", " ", SuperscriptBox["n", "2"]}]}]]}]}, { RowBox[{ RowBox[{"-", FractionBox[ RowBox[{"12", " ", "n", " ", RowBox[{"(", RowBox[{"2", "+", RowBox[{"2", " ", "n"}]}], ")"}], " ", "p"}], SuperscriptBox[ RowBox[{"(", RowBox[{"1", "+", RowBox[{"2", " ", "n"}], "+", SuperscriptBox["n", "2"]}], ")"}], "2"]]}], "+", FractionBox[ RowBox[{"12", " ", "p"}], RowBox[{"1", "+", RowBox[{"2", " ", "n"}], "+", SuperscriptBox["n", "2"]}]]}], RowBox[{ RowBox[{"-", "7"}], "+", FractionBox[ RowBox[{"12", " ", "n"}], RowBox[{"1", "+", RowBox[{"2", " ", "n"}], "+", 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.6623091302179728`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ss1", "=", RowBox[{"jacobianMatrix", "[", RowBox[{"0", ",", "0"}], "]"}]}]], "Input", CellChangeTimes->{{3.6623091012643223`*^9, 3.662309101669923*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"40.`", ",", "0.`"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", RowBox[{"-", "7"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.6623091302491727`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Eigenvalues", "[", "ss1", "]"}]], "Input", CellChangeTimes->{{3.6621746301969433`*^9, 3.6621746534039087`*^9}, { 3.662309105413929*^9, 3.66230910581953*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"40.`", ",", RowBox[{"-", "7.`"}]}], "}"}]], "Output", CellChangeTimes->{3.662309130264773*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ss2", "=", RowBox[{"jacobianMatrix", "[", RowBox[{"10.000000000000002`", ",", RowBox[{"-", "8.526512829121203`*^-14"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.6621638160807314`*^9, 3.662163836993928*^9}, { 3.6621639928428416`*^9, 3.662164005888588*^9}, {3.6623091090019355`*^9, 3.662309109360736*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "40.000000000000014`"}], ",", RowBox[{"-", "0.3067484662576686`"}]}], "}"}], ",", RowBox[{"{", RowBox[{"6.9185829116836185`*^-15", ",", RowBox[{"-", "6.008264462809917`"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.6623091302803726`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Eigenvalues", "[", "ss2", "]"}]], "Input", CellChangeTimes->{{3.6621747699894943`*^9, 3.662174776784913*^9}, { 3.6621748483928537`*^9, 3.6621748579668093`*^9}, {3.662309112433942*^9, 3.662309112870742*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"-", "40.000000000000014`"}], ",", RowBox[{"-", "6.008264462809917`"}]}], "}"}]], "Output", CellChangeTimes->{3.662309130295973*^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]", "2"}]}], "}"}], ",", 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}], Cell[BoxData[ GraphicsBox[{{}, {}, {RGBColor[0, 0, 1], Thickness[Large], LineBox[CompressedData[" 1:eJw91Xk4VdsbB3BDRAMaCVexjRmSkF+c/a6LSIYrFBJRuUVKcXOiqIS6QsjQ lQbTJeUSotKuSCkNjiGuMcIxy3hIzjm/zTn37OfZz16f513r+661/tnyh3zt PAX4+Pheke/Cl/OMAXeA1vpup9ZWjMMlV+ttK1LZ4NP40iC+cBIeiJ3UGTJg Q92SqjmnumnQ3eJ+htLIgs4MO18JpVm40RTrlu/PAnG3b0H78+bAb8a/bE6C BQwUXH7qBBOqpaPi+/KYYFpzMejKJ37UMaDlLEafh5vivl+6wwRRfOQ2r42X 5qG2TypCz1kIvfv7vmaIzDzQEoIlYxxEkavLJRs7m5/g5JDpf1xjFdLgv70L y/8BhWdiHc8fXI3uKc2OSJv/gDWvTM/REtcgv8YaBdOvswBmmkJLBdejaj3v 7qcSs4CeoAlLPUk0F8IWscidgfs552Wfe0ghr+QKzRjTGUjPYlKt0qSRitxT 0cxzDKha+eCT1CMZpG/bV5wuxYDbTjSLs8WyqGpntUtRyTTY1l9P2h8vh+zD uzstGFNQtPRB2dvjG1H/0b2xA8lTED7zT7G5zCaUvfT8/2p2TMF4cfvlfbRN iGX/+8HIsEm4cySM6bhJHiXRKG+91Sch0lLUPs9LHs1nB9pm10/AaLlBRs5D eXRs2lfZOHgC8Hjl25kD8ig/jnHISHUC+MbcrfnlFVDiHz66SQ3j0ED/xUPe XgGxH1bM2YSOg02hwSH9+wrIytjrQlPXGPilbg0d+qSAHF7at6TfGANLv6XN XaMKqHKiNKbSbAwSRzJSfVdgaEPZy0I/xnf4Hh59ulcZQ0R5mfFPj+8wc+NY UQyOoSUF6lIytFHYGhhm/cgeQwWKYmLVaBRKKKWh0b9jKGrCIY1dNAL57GhD NyqG/GbLhIrURqDNUHevcwSG2sWF9FvThoExHUrLv0F6iTqfr+wwuGv++SP6 Lobs7r2Q8PxrCILqJPpkcjHEfzqXVSQ1BBrriM/BRRhalXDe3i5lEPTiR7c2 lGEIszoary83CCqbH2sbvcZQhr44cswcgEm943K09xhaE+NPz9UYAAnKuEl2 DYb2OXbKqJb2Q+TriTuNDRjy7NUJf2PSD1JLjgyeasbQx8enmwPr+mDWNO1p RDuGPriv6jM83AcFN8eZyl0YYjfKxgsw6FCitL3HrQdDbqzIlspIOgSq+mvs 6MOQBZvZGihPB23GtsOvBzCktl7vo9SzXgh0DV63dRhDT89vXqdk1QttmteL 7UcxlPNwrfyBzh6Q2X3aQncMQ94Bhb8WBfTAd9oXsdZxDEW0iw5sFO8BS59X TJtJDLnOOlck3+8GocsffJKmMHTP+EKKhFk35O1MJ0qmMTQv9SbgYs83MM8b +lzMwFDzQP2F9rBvAKO73iXMYOgM5eRaORXSQU89988u9Hv7xOBDF8y5+7aJ /MBQ5mAlpni6C9ZFnBLJIF0TzNzWJNUFtw7VHlafw9B+rxeMXRWdILY8Si6T NEV0c8apE52g7pOpIP4TQ+n691J2SneC3vYUmZOkD7yZ+BIS9xVEdbXCKkiL 3zgX82C8A6ot2tJXzmNoz9mjCvx7O8CnBUu1JW0e8SPJt6wdKtA9ryjSjS6H 9/QptsNhg8/D5aRzfpWacIprA0uh/p5x0tlxRoGP+dqA2ry0xoRJnm9ZQsBU XwtcCV/W/Dtp/5OuVs8OtsCXG5kxEaRfTNqVb21pBhEP+X0ZpEXctmhZOzbD jItTOEGa5TJkym78F0zs9Xy+LNTp69t2OP8LxvEBVkOkPx+r0ppsb4Lm65Se VSwMvVRNqMkebYQVu2yYCqQNVoqWXAtqhD5+eqMO6fe6MyfLRBoh63OC/W+k A7zcFWnqX0DTsU77AOkQq9Ybd141wMQGY9NjpEU03zy65tgAoWkvXUNIr+vP d0u+Vg9WqQ0RV0l/T6xanaZaD66hWerxpLM+hLTcraqDxLdN0pmkFUsUlu9d Xge+ZmN/vSYd7Pl+OCGHBgXpXs+rSXf0FCIjBxqIbWm6U0u6kX9/Tgw/DYIS 7Uo6SP9zyX5YqeIzKDk8S2SQluVbYuTr+RHkJvxNNrIxJJb3Vbti2ztQWbU2 DCO9x23umXF5FSgI5eapkDY8mz7WaVsFpdS8F1tID40nxNoEvIWL3ncHcdJ4 5qozXf9Wwlm69Fdn0tMOfS/Fp8vh8b3EhjDSOlre119kPYPLI8N6z0kX1VLv TOzLA6vihNlB0rk3JXoCdePhCeRLzpEOUKOGHjjzB4iNKeiySH/Meu/NJ+wM RPz0BTZp+mdq3LVSQ6gTqGldcJHTT5ktJzaCdpTe9IKV8w8yu44Igin76/yC owOGn7c/6MVfRecxF1w94h4mLNOJq/U6shacwh8wbjXShJfvHV70VIRzwq3K WnyGz569YMdIhl32kWrc2yZqsZ7+sPK2xjkCbze6tug1IL+l9kQBvkv34+L8 pwp27nyn0/BukYFFRyZ4NWtR4/GYczWLbi1oiS71CMOZZimLDujNsSztCMB/ zk8vmt5Z4tBBP4rLlv1YtMHleo+uIkdc9htz0Yr3/x6Vx83wC8sXfvOK6NVt vtHxDbq4uZLwooWL8MikSC0cXDh1MZ3Vy14fUsbzdwgsuvRK69neAC18uT57 MS/6ouiOqN0quGoSx89Uppo/BKnguk851kh+4CG2VQ6fK+DYjSEZbSK+DvdZ x83vfVDlPymI19A59Sc6GX3hvT0UtRCOhYP27qa6txKG3H7XZEaDK8qHiBJb juWu5AgZvRsimtZy8ho9qNbJ/l3EuRiOa9mVahtmO4gVyRzbpDfNDa3vJXI9 OX5dZ3XeOmWYGFXiWHm7d/hdnzHCks3Jt/HKkw7+ZYRIW8qpO1VaXIz06SS6 V3LscX9T22+uDcQZMY5ZY2+HRFi1xG3u/SaUSgZWPqwn4oQ5fvul6YNgeANx kZ/jTRHG1PaVnwhRFqffZZ35rKMtz4jV8xyfFowSYJncJhx+crzCq8uSLzyU 4s21XVf/YNnKJIoad77GmEoniouh2DA5RtbxkBrnSwRy85Ovui+bm7lJhHDP JxkZMrX/biYhy8e9H7RzVjX3b0KJa8kxio7w3rsEH9dCtrFWOXJRxCB3feyE i66N03Gii+urO3dZD5fvp3Ry7cR/1lXayYvSx/XHbO2jqlsOUsa5Vj2v+lKb ak385LraMzI6bqcvIcjt93y5X6HcnrPEUq4ZdPMyB79AQojrYrEQ12bXU4QA 19Td/W62Igd4+92dty1zslOPYHPzC5Qe3bKDbZT/fDPf0Ky71IjnY90R/0ys 1qT8t/6UpvFxQ4HNBD/X+RF/HDt21ZRnd5/wTSc0rHnOMKf9Wa9nyXOds75L TgDw9mNWT1Fws9bmee1H45gCTInnMOvSVJPHijxLrJ7PeF6vwrNnbOyuegdV nllIUtOlcRPP2kLaNrWEIM+qWo9G74IAz16GBt525kt4pvV890ymsJ//H/nA orI= "]]}}, Axes->True, AxesLabel->{ FormBox["\"N\"", TraditionalForm], FormBox["\"P\"", TraditionalForm]}, AxesOrigin->{0, 0}, Method->{}, PlotRange->{All, All}, PlotRangeClipping->True, PlotRangePadding->{Automatic, Automatic}]], "Output", CellChangeTimes->{3.662309130389573*^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:eJxF1nk0Vd0bB/BrnqcyJDIdMiVKRJyzNxWlkShKhhs/RYbqkhdJhjeKCOUV KoqiKCQlQ2QupV5JyhCuynhxzdPv3HN137PWWXt91nc/zzlr7/3HViR7W7my k0ikVvxljMyHFzJHURjD5y7LxscPM7Pp5k6RIpAiWMM7pC4I9XlTaTtEReB+ /lv5EtaikCtvU5CGvDCc4bj2ylNSAt7pLCj7iArCXT2/KwJ3SUI5m6e1N2oF oB1/Y0uKnxRcr1VzVeKgAHz8uqqvqlwaSrzz1mtw5Ycjzf8UrPu8Fj66O8mr P84HJcOOfDXrloH/U9R0bb7IB3WkyzZENK+DH8jHkztv80JK6ryXyDM5qEIe mDmnwwvDbacD90XJw71nFy0C3vDAwkzuIRE1BXj6kAy6cIQHxn/bNr0/XAEW qj7yYB/hhi9V9iRd/qoA3wkERyZEcMNyLEuTpqoI6S7tD57JccPoJmy1X74i tKTk5RbacMGALGoeNqwIfR9YNCdOcMLEjpDrlxElKIU4FQsmcEJdj1NF+dZK cIeJ+fm1Wzjh6keoktwlJUi+Yu9c2soBTX7ucNmQrQT3Hd8gSwvggJeTPnEY NinBS83iRQUKHHBIlHI5dVgJCnrDV6R6dkiJvUrL4UdgC3efdasPO1STTdMo VEbgTOCdvdoy7HDWLeLukhEC7x8NleepY4PJ3Hu8Ph7E89xBITsKG7z5PTXE /QQCT5bD9QoIGzyDtAhxUBAoUB5vceRfEjwWWt/yIRSB2Pd6OBtOgv7mcldI cQj8H3n9M5IBCbpY7HF8nYJAkxLLC81tyyBN98QL/SwEvvS+12G9YxlYpHve +PsJAoPDGm6GFSwBz9oT/M3FCAz5UZBviSwBgzJ7J+0KBEbegYbPbi6Ck646 NuU1CHz1RBreFVwEBxrFqhPeIlBXqaifP2IBFKjVr6tuRuC3plin38vzwNL2 YqPzZwR6yA3RdS7Mg/pZoamgrwicFDfl/bIwByqom+2VOxDotdpf4mPwHJhI Gity6UZg/D2jcnGOOYDNfm3Y3YvAfwTTXiZFzQJ36S0ZHVQEqoQHDx1YNQsa uwNj1X8hMCcQtdJNmwHNCsohegMIfHdeSMpAYwaUqvY38Q4h0LW+ONj65TTw /36+JHsYgfPPw5aDLaZBh3qhiNwoAvUSldXyOqbAgds7S91pCNz4Ub+k7ewU SGoJ2pc4hsCBistxk3xTgGtRJOv2OAKvHdr5Yj5jEuTqatrFTCDQd32kbQ86 CTh8FfxO0BF4gna94W47HSh3LTYikwjMSjEK3vgXHRzmfBvUiVvPhj34qhgd hFTybY6YQqD3W5879EcToPat3Z610wikuTmlOu6aALOo4Ogd3IeHL1hX9Y8D Ni+MJjWDr4+nOHXV5XFQNDgyFYa7RHxjrLn6OHDeVBOHziLwk3+q+dZzY+CK 6mhwDO6c6csfutaOASGuhY2tuLnb1VzNq2lgXb3Y+TVzCOyQXb/N1ZsGkJE3 kja4bZU8s3RkaUAPmS2Mxr0lgM/Z9MYoMP925UcFbmh3UuT45Aiwq8nrWDuP QIsXQvZL5cMg9RVbyHbcitTUw65qw8DDlOJwCveLqJzE2sQhIJr+kp6Lu++D G4+x3yCYGu5rV1hAYAJFy8bx6G/gI188aYh7gxCfOkf7L6Cvtb/XEvfNqkvV Ova/gLZp9t1g3JopHOZTbj/Bj4ctk024xeUTDn9KoYKorgGhbtw2aVsq3TdS gc0z9WtjuH17ZNi03vSByGyB0VWLCKyWNYjhmOgFziGB0gdxx5/97PGvew9I 0nOae4G7J5CiOOPSBWQVRgbslvD1DQ512f6+FWzLrqtwxp3MPlTlSG4FaHml 1yncfYY8CQEzn8HgLjMHf9y+Bcu5ERqfQXKn0j+JuDU4FVZP3v8XSMhLCdbh ltnkpPt5sBk8Jf1ql1pGINe9L990ZBtAA1hnFILbVPWnVk7mK9AUnn28DbfD Rq3cTNV0kGQv1knHXdlgmWNqHgZSnWDyAu75YfaS6EQ3ILg5SXgZt/8pKS4T OQvgwI+VMjxyW1uO97YWuO9tOMhw65mMiCCzVeD5yWLClKTYX+oDYsAxSG+I 4Wxy4kBOrBhg+7R5mGHfkF+SnN2iYHZ5yyjDgmnrQzU1RAB8Wkh4KtYyvCZa GBxh06Ex3N7Z0/jXuBCwmpIbY/j8E8mUqVpB4GYVR/gRfNiVricIHlUsEG7b G9y54YEAUMGqxxkeon8OeRDPDy5Z+UwwHEv96jkWyQcoq5cJi9qJv06M4gW6 q87QGX696PVbI5EHGD15QxgliTQfecgNDm0lTTKs6VgXzjXECex8c4j84cu/ 0VJLDrCFuobw2W6eXx9b2UB4ehHRf8K4Zd28PwmsMWN6k0+EQ9jDZazfi05Y ypB8jeaxhGV81SPqvyMz1VafFjF0ZzxhCR8uoVKjRYw3WpD4/nvV5lmZhwuY 4aUbhP0v7t3hKrWAuUkrTDFs32hgdiVqHos0yCJMNnRBgpbmsMQqxWmGryym LWv6zWEyj6MJn97hw3GVNovF39k2w7BM5sdu5eEZbLLxOeHM6t4nPRPTWJW1 zizDFwTWZWnLTWHhrwoJv020Zrc9Tce2B+2eY9ho9ZKx6eQ4VnewibCN5ptP YsljWHy3wTzDxe/ai1ttaJhIThzhxxK3lDOsRzAvFYUFhmcqJ3bmqQxgAY65 hK09nJTqIvsxkrP1IsOX/PgVYjx6MWqTyBLDFeO6Ww+UdmPha14RLlLpbUPa vmOBhXLLxH6O9dMHv7VhsZObidwfmkRoV77Fvj12J/IEuupfB6bKsSuz2YRt 09lTTYOfYkqalcR890WbuJdB6djgdCWRk2NVzxtaJmDO0oL4vUgZzgGKj9Xu COyQvCRhs/dcKli4L8Y3sEDMb/G5fvLdA1fMJHmKsG+9Mu89UVuMfOY9YVUV 82PBvQcx3iddhCVmRrSH3u/GnDM5iH5tNTtJ+aWm2IvnNcz/owrTHFYbYmxL M4RdskxO7X27CevGTzvDoc6yZ37ip9/hEbN+W2+dWrKuNEZNZ+Yp9NGkn7IC WFAT07YlHG12nHT07PeV+oX1Dbs3dKM8Mkzf27Pc4in5AZVvXyLcQ2vTD5Up R7EBpjOeRmSmRxSiuVuY8wd9v/FyqxSj5BtM7wPXPdUPv0bNepgWLqd+V2F/ gyYTl0xleGS/FdX6agXqMM/MFY4d02/gz0fPTTN931HL+0n+XZRYPtwHTU/X G1PiULUZpoeetYjvvBmN1s4xvXdPtEyjcwIas8j0Vn1Faq1vMsrYXYbPmdj8 lRd6C+1dMf3LSar/sTi0aMWfAsSKBpoD0bQV+0lcvX/7mTWavGJjo5gMGbJ1 2d0/9fGdfc9TjpcVrPhX/vbslP59ZQ0rfmp+diTgvAXateK81uSidCcXdGzF X+3uqWNRnujCikOFvfiAuie6uOI6pVynv685s+xXwtcoVbWX5YzIFl1ylwK6 tGJKcaVZ25B62R8/ltyfeOqWTtnyittkyE3JhrIsh5u2tkRxIegfdzTK3Tg+ aYAy7/3KUHj0SFGJjgnLnlbD5keFIct9biLGp3U3s6yWSj1pNPDDmGV7g8yU G6SyP34t2WN+eyM7yy/M+DyslPpL/5hs3xJd+nueVX+1ppykmy/C6t/iED4k ICzOct27DUeT94uyfHiQRLF6zsYyV/8av+KIcVa/rolE8wb6b5bDsthno44P spxmnp2Q82OU5eXcL03BMf85S+NxSfS5Ppa5tXe9OTtfx3J0x3vjgAv/2fe1 083Tdo0s08qaLgpYVRv/H+8cak8= "]]}}, Axes->True, AxesLabel->{ FormBox["\"N\"", TraditionalForm], FormBox["\"P\"", TraditionalForm]}, AxesOrigin->{0, 0}, Method->{}, PlotRange->{All, All}, PlotRangeClipping->True, PlotRangePadding->{Automatic, Automatic}]], "Output", CellChangeTimes->{3.6623091304207735`*^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:eJxF1Qk0VO0fB/ARMgkZxpKs89D7NlJvRMi990clabPlTWhR/FNka6Fea1JJ mzEqtKimRaGIlCZLKO1Cb4VQkSWKBsNY3suM+d9z7rnnc37P/X6fe8895+p5 BTh5T6FQKBnkOX4VHnQQXtXBV2flXM58FYjt4kg8SFWDxHCd5bc9VYEzsOTw wCI1oJoeuGJ8Wx1cu2a47g1WheeHjC1zn82ETKpffLuiKuy0fCKdU6cBHnpy n5ffUYHPS04ocBo1gQpHp9N66TD3kYzOoddawGNqpaYm0UFwr7CnMFsb9oXA LwdzOthd+IV3xuhAxLoG9uIGZbCoqXvJR7rg7hdote2gMozxZGuDInUBuzfz axVTGR7L0BZMqdEF79n5lPhqJWBFeQba++rBm2vTi9r+VAKnhUed2Bl68CnJ YYRVS4MTc36uUm3VA41MCxPWQRr8I6mZPEeTAXkW6i6txjQANcVao9UM+CuU zU74qgjuRlXD0WEMOOu49UgsWxG6pDmJB9MZkDg/eN5LO0VAhgy34+UMKCd+ apcOzYD19A7b/hYG7NLGWFE+M0CjvaJFURKBwhXc9lWNAiT4b21iaiGwPGsz L32ZAoQFxt5MWIjA93DxE16BPMx+/mpr2AoE7o1KasXz5OFUzNjIPA8EJjGV 34auyUG50RvnVj8EUyMGOGl6chDbXOFZfQCBbenmC2kXpkNGoLKfwVEEzKOK 7T1a0+Ew0lguw0ZwTPFFDOuiLOgJzOouXkSQejJ3XSiShQLmCiulmwiCs92J MzenQVUxuz78LgINQa/69wXT4JzKqr/GCsg+NzPuxkdUOBxHycouQtBgukm+ z44KdBxTvFaOIFr/Vv+lf2XATCmuqv85gue5B3e7bZeBssqKuffeIKhay0tS EUwF82bVTd3VCJLCqw+8PDkV0peHUO/9i2CeponXntlTQTVusI1Wh6Aw2uai TJE0RAxChfJnBDKd9Wn73aRBFw+98aQJQchj7c7yPinwM6xvZX5FUKxvMLeZ JQUv2zfscmxBMJ9eaVhmIgXaa4sv4d8ReCf//cynVhJSlrSU8NsQaNtZTC8J lQTGwIv0Qx0Iln7cT3uhJQnDF5xMOzoR9J4pYISXTYH2dyw5wy4EciqMtDd+ U8DgPJ2+qhvBJdYuyULVKVDf7nnE6hcCn/W4OREgAcOrLetpPQguWGter9SS gPq+YIdq0m2UKpu3ryhQv1Je5WAvAlVe6p+rIinASo6ehX4jkLpzMMnMmAJe MTkaeaRTmnaW0l6NEe9urU215CFwrl4/Pcx8jFh2/+pGgz4El/c9TfJQGyWu PU+2ZJO+32alIRs5QrSqrv4oIH0sGqwftg8TnqyVAXf7ye8tgCPV+lRA3Edu aMoAAl5Ers7SxQJiwdQz6qtJl6RHNvjeHSKU9rea1ZDeH1biUsMZJIaea2+n 8RF83CHzcgQNEueaiR570lvlz7sVXOUT13KzIu+Rnr34nXJm5gBxOe1krO8g AnvZsxEW9X2Eeo/Ri9Okt6QFrj3h30cwlqf6FZB2DflgrSzRR7xZV2crMYRA 4slmVvQ8HrEpjOYUSzqKUh38NL2XqLwRnrBJgGAn4zheWPqTUAmifNhH2uyI 2dduvZ/EgvP84ydJl3kcjXKI7Sb0bSMcuaTZpvL5nk5dROegzVOlYQRFq014 v2idhKbnsk83SJt0uviEKX0nnoXWX2GNkO8LMza5VtpE8La81b1ImnlnWZhH SBOxonDYJ4P07ZaHj44aNBGWztrOxaTtAAbcNRuJrNzsinbSVXGXIsob6wnT L14dJqPkvKFAWbfxA6GdyM/PIH166G4PtvsdYd/pZ7R0DMFr88uRbMcyohZb yy0nnX5Z97fuoptEx1bryz9JH7j14QrTNYE4zPEbGiI9vDIqZ7dLEGGygzg/ Sjon2XLxyAJnIqs688wY6YwfKmx6mSnxeVlv07jnttl1P6zXIPJ/dPWPm8UJ p4ZKjeIpAytGxh1U4l4ZZv0Fv3VqbMLrXeToBcvqcLWsQ6PjVnasKcnaWIM7 zmqccHyPVn58yGv8/umBCTcHUcOriypwU2/3sXHb+PAsc5gFeGpy04Tfqd6R Lc+8ib+e1jzhmGy+oSAgDY9o65zwNx/5WRTvk7hxzeCEG65YeZcHRuGf/nEk /0P6UGr7sMQ9PAjnSPZOzPewHr950O2FZzbET9gwtXvQ38QZ3/1a2Lexbo/8 YusluP+GsQmfv+73+4u9Mb4hlDKR9yk3L+HQkB5uPyycr1l0qoapqoHr8YT2 KnhfFJmkiLeL5rHebnuM3svgaQNCb8PsZ1JURrFeWWHe2+RVXVsUO7BtKkIv XCS1bm9LBaYwU2ht6nLXwcJibq2a0Lf9+Z2LDCq5h5SF5qfw25hKXG63nNDr 8zo22MTdxahThVZYrLa5XPAWq6UI3TJ74+a5HR8wu1Hhfkrr15VV5X3ApMeE NjzXsN0k/xXWLnLF9RijRYEPsVDR/YY8P8+V7EQsR+Sanzt8q36wuZki8xxv BzXrpHFjJud7G8NDdxzjYiKbvw/9Hdt3DBsQ5f/I+Zq+VIuDPRM5nnLKI3hN JvZE5I8prjt+SWVgTSLfv0LMWUVLw/pFTtEv34IFx2GSonwqY79+wqwlGFVk /raLo7eld3FlRD713V+qxcWfO7ne/BmzN13bgTsqyrNzijW/KfDGBkWm8Krt e6zDxH1JFsbrNjdFiffvvMWEIev3DyaYnGv6ypmx/LHJvMbkRikDD1fh45N9 lL6XNBqd90hC5Hj04NjSDmvupNfvMg0JCye4k+sDRxI3RDLo3DFR3r6Avd41 XzBs0jN2VmQz89aI/WBrX9cfPs5iu4R2bt79a7W435DKn3G6HRO7yPb1l+Iz +mI7bDvBzO/tezRprwWV57+bq4j3c0mHTlh0K4hdPCcxePumb1aTXhgYIbB8 qCrOM/DSC2zZh8RO3Cm54nGbgdhH/ohGeziaYvd//T1EXSEjdn7KEU31Bmmx wenzjKdXe8R9u/gbHQJ4/+//G+U9+Z9Wi9ih9CLfQPkmq/8AMGWvjQ== "]]}}, Axes->True, AxesLabel->{ FormBox["\"N\"", TraditionalForm], FormBox["\"P\"", TraditionalForm]}, AxesOrigin->{0, 0}, Method->{}, PlotRange->{All, All}, PlotRangeClipping->True, PlotRangePadding->{Automatic, Automatic}]], "Output", CellChangeTimes->{3.662309130451973*^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:eJw91Xk4VdsbB3BDRAMaCVexjRmSkF+c/a6LSIYrFBJRuUVKcXOiqIS6QsjQ lQbTJeUSotKuSCkNjiGuMcIxy3hIzjm/zTn37OfZz16f513r+661/tnyh3zt PAX4+Pheke/Cl/OMAXeA1vpup9ZWjMMlV+ttK1LZ4NP40iC+cBIeiJ3UGTJg Q92SqjmnumnQ3eJ+htLIgs4MO18JpVm40RTrlu/PAnG3b0H78+bAb8a/bE6C BQwUXH7qBBOqpaPi+/KYYFpzMejKJ37UMaDlLEafh5vivl+6wwRRfOQ2r42X 5qG2TypCz1kIvfv7vmaIzDzQEoIlYxxEkavLJRs7m5/g5JDpf1xjFdLgv70L y/8BhWdiHc8fXI3uKc2OSJv/gDWvTM/REtcgv8YaBdOvswBmmkJLBdejaj3v 7qcSs4CeoAlLPUk0F8IWscidgfs552Wfe0ghr+QKzRjTGUjPYlKt0qSRitxT 0cxzDKha+eCT1CMZpG/bV5wuxYDbTjSLs8WyqGpntUtRyTTY1l9P2h8vh+zD uzstGFNQtPRB2dvjG1H/0b2xA8lTED7zT7G5zCaUvfT8/2p2TMF4cfvlfbRN iGX/+8HIsEm4cySM6bhJHiXRKG+91Sch0lLUPs9LHs1nB9pm10/AaLlBRs5D eXRs2lfZOHgC8Hjl25kD8ig/jnHISHUC+MbcrfnlFVDiHz66SQ3j0ED/xUPe XgGxH1bM2YSOg02hwSH9+wrIytjrQlPXGPilbg0d+qSAHF7at6TfGANLv6XN XaMKqHKiNKbSbAwSRzJSfVdgaEPZy0I/xnf4Hh59ulcZQ0R5mfFPj+8wc+NY UQyOoSUF6lIytFHYGhhm/cgeQwWKYmLVaBRKKKWh0b9jKGrCIY1dNAL57GhD NyqG/GbLhIrURqDNUHevcwSG2sWF9FvThoExHUrLv0F6iTqfr+wwuGv++SP6 Lobs7r2Q8PxrCILqJPpkcjHEfzqXVSQ1BBrriM/BRRhalXDe3i5lEPTiR7c2 lGEIszoary83CCqbH2sbvcZQhr44cswcgEm943K09xhaE+NPz9UYAAnKuEl2 DYb2OXbKqJb2Q+TriTuNDRjy7NUJf2PSD1JLjgyeasbQx8enmwPr+mDWNO1p RDuGPriv6jM83AcFN8eZyl0YYjfKxgsw6FCitL3HrQdDbqzIlspIOgSq+mvs 6MOQBZvZGihPB23GtsOvBzCktl7vo9SzXgh0DV63dRhDT89vXqdk1QttmteL 7UcxlPNwrfyBzh6Q2X3aQncMQ94Bhb8WBfTAd9oXsdZxDEW0iw5sFO8BS59X TJtJDLnOOlck3+8GocsffJKmMHTP+EKKhFk35O1MJ0qmMTQv9SbgYs83MM8b +lzMwFDzQP2F9rBvAKO73iXMYOgM5eRaORXSQU89988u9Hv7xOBDF8y5+7aJ /MBQ5mAlpni6C9ZFnBLJIF0TzNzWJNUFtw7VHlafw9B+rxeMXRWdILY8Si6T NEV0c8apE52g7pOpIP4TQ+n691J2SneC3vYUmZOkD7yZ+BIS9xVEdbXCKkiL 3zgX82C8A6ot2tJXzmNoz9mjCvx7O8CnBUu1JW0e8SPJt6wdKtA9ryjSjS6H 9/QptsNhg8/D5aRzfpWacIprA0uh/p5x0tlxRoGP+dqA2ry0xoRJnm9ZQsBU XwtcCV/W/Dtp/5OuVs8OtsCXG5kxEaRfTNqVb21pBhEP+X0ZpEXctmhZOzbD jItTOEGa5TJkym78F0zs9Xy+LNTp69t2OP8LxvEBVkOkPx+r0ppsb4Lm65Se VSwMvVRNqMkebYQVu2yYCqQNVoqWXAtqhD5+eqMO6fe6MyfLRBoh63OC/W+k A7zcFWnqX0DTsU77AOkQq9Ybd141wMQGY9NjpEU03zy65tgAoWkvXUNIr+vP d0u+Vg9WqQ0RV0l/T6xanaZaD66hWerxpLM+hLTcraqDxLdN0pmkFUsUlu9d Xge+ZmN/vSYd7Pl+OCGHBgXpXs+rSXf0FCIjBxqIbWm6U0u6kX9/Tgw/DYIS 7Uo6SP9zyX5YqeIzKDk8S2SQluVbYuTr+RHkJvxNNrIxJJb3Vbti2ztQWbU2 DCO9x23umXF5FSgI5eapkDY8mz7WaVsFpdS8F1tID40nxNoEvIWL3ncHcdJ4 5qozXf9Wwlm69Fdn0tMOfS/Fp8vh8b3EhjDSOlre119kPYPLI8N6z0kX1VLv TOzLA6vihNlB0rk3JXoCdePhCeRLzpEOUKOGHjjzB4iNKeiySH/Meu/NJ+wM RPz0BTZp+mdq3LVSQ6gTqGldcJHTT5ktJzaCdpTe9IKV8w8yu44Igin76/yC owOGn7c/6MVfRecxF1w94h4mLNOJq/U6shacwh8wbjXShJfvHV70VIRzwq3K WnyGz569YMdIhl32kWrc2yZqsZ7+sPK2xjkCbze6tug1IL+l9kQBvkv34+L8 pwp27nyn0/BukYFFRyZ4NWtR4/GYczWLbi1oiS71CMOZZimLDujNsSztCMB/ zk8vmt5Z4tBBP4rLlv1YtMHleo+uIkdc9htz0Yr3/x6Vx83wC8sXfvOK6NVt vtHxDbq4uZLwooWL8MikSC0cXDh1MZ3Vy14fUsbzdwgsuvRK69neAC18uT57 MS/6ouiOqN0quGoSx89Uppo/BKnguk851kh+4CG2VQ6fK+DYjSEZbSK+DvdZ x83vfVDlPymI19A59Sc6GX3hvT0UtRCOhYP27qa6txKG3H7XZEaDK8qHiBJb juWu5AgZvRsimtZy8ho9qNbJ/l3EuRiOa9mVahtmO4gVyRzbpDfNDa3vJXI9 OX5dZ3XeOmWYGFXiWHm7d/hdnzHCks3Jt/HKkw7+ZYRIW8qpO1VaXIz06SS6 V3LscX9T22+uDcQZMY5ZY2+HRFi1xG3u/SaUSgZWPqwn4oQ5fvul6YNgeANx kZ/jTRHG1PaVnwhRFqffZZ35rKMtz4jV8xyfFowSYJncJhx+crzCq8uSLzyU 4s21XVf/YNnKJIoad77GmEoniouh2DA5RtbxkBrnSwRy85Ovui+bm7lJhHDP JxkZMrX/biYhy8e9H7RzVjX3b0KJa8kxio7w3rsEH9dCtrFWOXJRxCB3feyE i66N03Gii+urO3dZD5fvp3Ry7cR/1lXayYvSx/XHbO2jqlsOUsa5Vj2v+lKb ak385LraMzI6bqcvIcjt93y5X6HcnrPEUq4ZdPMyB79AQojrYrEQ12bXU4QA 19Td/W62Igd4+92dty1zslOPYHPzC5Qe3bKDbZT/fDPf0Ky71IjnY90R/0ys 1qT8t/6UpvFxQ4HNBD/X+RF/HDt21ZRnd5/wTSc0rHnOMKf9Wa9nyXOds75L TgDw9mNWT1Fws9bmee1H45gCTInnMOvSVJPHijxLrJ7PeF6vwrNnbOyuegdV nllIUtOlcRPP2kLaNrWEIM+qWo9G74IAz16GBt525kt4pvV890ymsJ//H/nA orI= "]]}}, {{}, {}, {RGBColor[1, 0, 0], Dashing[{Small, Small}], LineBox[CompressedData[" 1:eJxF1nk0Vd0bB/BrnqcyJDIdMiVKRJyzNxWlkShKhhs/RYbqkhdJhjeKCOUV KoqiKCQlQ2QupV5JyhCuynhxzdPv3HN137PWWXt91nc/zzlr7/3HViR7W7my k0ikVvxljMyHFzJHURjD5y7LxscPM7Pp5k6RIpAiWMM7pC4I9XlTaTtEReB+ /lv5EtaikCtvU5CGvDCc4bj2ylNSAt7pLCj7iArCXT2/KwJ3SUI5m6e1N2oF oB1/Y0uKnxRcr1VzVeKgAHz8uqqvqlwaSrzz1mtw5Ycjzf8UrPu8Fj66O8mr P84HJcOOfDXrloH/U9R0bb7IB3WkyzZENK+DH8jHkztv80JK6ryXyDM5qEIe mDmnwwvDbacD90XJw71nFy0C3vDAwkzuIRE1BXj6kAy6cIQHxn/bNr0/XAEW qj7yYB/hhi9V9iRd/qoA3wkERyZEcMNyLEuTpqoI6S7tD57JccPoJmy1X74i tKTk5RbacMGALGoeNqwIfR9YNCdOcMLEjpDrlxElKIU4FQsmcEJdj1NF+dZK cIeJ+fm1Wzjh6keoktwlJUi+Yu9c2soBTX7ucNmQrQT3Hd8gSwvggJeTPnEY NinBS83iRQUKHHBIlHI5dVgJCnrDV6R6dkiJvUrL4UdgC3efdasPO1STTdMo VEbgTOCdvdoy7HDWLeLukhEC7x8NleepY4PJ3Hu8Ph7E89xBITsKG7z5PTXE /QQCT5bD9QoIGzyDtAhxUBAoUB5vceRfEjwWWt/yIRSB2Pd6OBtOgv7mcldI cQj8H3n9M5IBCbpY7HF8nYJAkxLLC81tyyBN98QL/SwEvvS+12G9YxlYpHve +PsJAoPDGm6GFSwBz9oT/M3FCAz5UZBviSwBgzJ7J+0KBEbegYbPbi6Ck646 NuU1CHz1RBreFVwEBxrFqhPeIlBXqaifP2IBFKjVr6tuRuC3plin38vzwNL2 YqPzZwR6yA3RdS7Mg/pZoamgrwicFDfl/bIwByqom+2VOxDotdpf4mPwHJhI Gity6UZg/D2jcnGOOYDNfm3Y3YvAfwTTXiZFzQJ36S0ZHVQEqoQHDx1YNQsa uwNj1X8hMCcQtdJNmwHNCsohegMIfHdeSMpAYwaUqvY38Q4h0LW+ONj65TTw /36+JHsYgfPPw5aDLaZBh3qhiNwoAvUSldXyOqbAgds7S91pCNz4Ub+k7ewU SGoJ2pc4hsCBistxk3xTgGtRJOv2OAKvHdr5Yj5jEuTqatrFTCDQd32kbQ86 CTh8FfxO0BF4gna94W47HSh3LTYikwjMSjEK3vgXHRzmfBvUiVvPhj34qhgd hFTybY6YQqD3W5879EcToPat3Z610wikuTmlOu6aALOo4Ogd3IeHL1hX9Y8D Ni+MJjWDr4+nOHXV5XFQNDgyFYa7RHxjrLn6OHDeVBOHziLwk3+q+dZzY+CK 6mhwDO6c6csfutaOASGuhY2tuLnb1VzNq2lgXb3Y+TVzCOyQXb/N1ZsGkJE3 kja4bZU8s3RkaUAPmS2Mxr0lgM/Z9MYoMP925UcFbmh3UuT45Aiwq8nrWDuP QIsXQvZL5cMg9RVbyHbcitTUw65qw8DDlOJwCveLqJzE2sQhIJr+kp6Lu++D G4+x3yCYGu5rV1hAYAJFy8bx6G/gI188aYh7gxCfOkf7L6Cvtb/XEvfNqkvV Ova/gLZp9t1g3JopHOZTbj/Bj4ctk024xeUTDn9KoYKorgGhbtw2aVsq3TdS gc0z9WtjuH17ZNi03vSByGyB0VWLCKyWNYjhmOgFziGB0gdxx5/97PGvew9I 0nOae4G7J5CiOOPSBWQVRgbslvD1DQ512f6+FWzLrqtwxp3MPlTlSG4FaHml 1yncfYY8CQEzn8HgLjMHf9y+Bcu5ERqfQXKn0j+JuDU4FVZP3v8XSMhLCdbh ltnkpPt5sBk8Jf1ql1pGINe9L990ZBtAA1hnFILbVPWnVk7mK9AUnn28DbfD Rq3cTNV0kGQv1knHXdlgmWNqHgZSnWDyAu75YfaS6EQ3ILg5SXgZt/8pKS4T OQvgwI+VMjxyW1uO97YWuO9tOMhw65mMiCCzVeD5yWLClKTYX+oDYsAxSG+I 4Wxy4kBOrBhg+7R5mGHfkF+SnN2iYHZ5yyjDgmnrQzU1RAB8Wkh4KtYyvCZa GBxh06Ex3N7Z0/jXuBCwmpIbY/j8E8mUqVpB4GYVR/gRfNiVricIHlUsEG7b G9y54YEAUMGqxxkeon8OeRDPDy5Z+UwwHEv96jkWyQcoq5cJi9qJv06M4gW6 q87QGX696PVbI5EHGD15QxgliTQfecgNDm0lTTKs6VgXzjXECex8c4j84cu/ 0VJLDrCFuobw2W6eXx9b2UB4ehHRf8K4Zd28PwmsMWN6k0+EQ9jDZazfi05Y ypB8jeaxhGV81SPqvyMz1VafFjF0ZzxhCR8uoVKjRYw3WpD4/nvV5lmZhwuY 4aUbhP0v7t3hKrWAuUkrTDFs32hgdiVqHos0yCJMNnRBgpbmsMQqxWmGryym LWv6zWEyj6MJn97hw3GVNovF39k2w7BM5sdu5eEZbLLxOeHM6t4nPRPTWJW1 zizDFwTWZWnLTWHhrwoJv020Zrc9Tce2B+2eY9ho9ZKx6eQ4VnewibCN5ptP YsljWHy3wTzDxe/ai1ttaJhIThzhxxK3lDOsRzAvFYUFhmcqJ3bmqQxgAY65 hK09nJTqIvsxkrP1IsOX/PgVYjx6MWqTyBLDFeO6Ww+UdmPha14RLlLpbUPa vmOBhXLLxH6O9dMHv7VhsZObidwfmkRoV77Fvj12J/IEuupfB6bKsSuz2YRt 09lTTYOfYkqalcR890WbuJdB6djgdCWRk2NVzxtaJmDO0oL4vUgZzgGKj9Xu COyQvCRhs/dcKli4L8Y3sEDMb/G5fvLdA1fMJHmKsG+9Mu89UVuMfOY9YVUV 82PBvQcx3iddhCVmRrSH3u/GnDM5iH5tNTtJ+aWm2IvnNcz/owrTHFYbYmxL M4RdskxO7X27CevGTzvDoc6yZ37ip9/hEbN+W2+dWrKuNEZNZ+Yp9NGkn7IC WFAT07YlHG12nHT07PeV+oX1Dbs3dKM8Mkzf27Pc4in5AZVvXyLcQ2vTD5Up R7EBpjOeRmSmRxSiuVuY8wd9v/FyqxSj5BtM7wPXPdUPv0bNepgWLqd+V2F/ gyYTl0xleGS/FdX6agXqMM/MFY4d02/gz0fPTTN931HL+0n+XZRYPtwHTU/X G1PiULUZpoeetYjvvBmN1s4xvXdPtEyjcwIas8j0Vn1Faq1vMsrYXYbPmdj8 lRd6C+1dMf3LSar/sTi0aMWfAsSKBpoD0bQV+0lcvX/7mTWavGJjo5gMGbJ1 2d0/9fGdfc9TjpcVrPhX/vbslP59ZQ0rfmp+diTgvAXateK81uSidCcXdGzF X+3uqWNRnujCikOFvfiAuie6uOI6pVynv685s+xXwtcoVbWX5YzIFl1ylwK6 tGJKcaVZ25B62R8/ltyfeOqWTtnyittkyE3JhrIsh5u2tkRxIegfdzTK3Tg+ aYAy7/3KUHj0SFGJjgnLnlbD5keFIct9biLGp3U3s6yWSj1pNPDDmGV7g8yU G6SyP34t2WN+eyM7yy/M+DyslPpL/5hs3xJd+nueVX+1ppykmy/C6t/iED4k ICzOct27DUeT94uyfHiQRLF6zsYyV/8av+KIcVa/rolE8wb6b5bDsthno44P spxmnp2Q82OU5eXcL03BMf85S+NxSfS5Ppa5tXe9OTtfx3J0x3vjgAv/2fe1 083Tdo0s08qaLgpYVRv/H+8cak8= "]]}}, {{}, {}, {RGBColor[0, 1, 0], LineBox[CompressedData[" 1:eJxF1Qk0VO0fB/ARMgkZxpKs89D7NlJvRMi990clabPlTWhR/FNka6Fea1JJ mzEqtKimRaGIlCZLKO1Cb4VQkSWKBsNY3suM+d9z7rnnc37P/X6fe8895+p5 BTh5T6FQKBnkOX4VHnQQXtXBV2flXM58FYjt4kg8SFWDxHCd5bc9VYEzsOTw wCI1oJoeuGJ8Wx1cu2a47g1WheeHjC1zn82ETKpffLuiKuy0fCKdU6cBHnpy n5ffUYHPS04ocBo1gQpHp9N66TD3kYzOoddawGNqpaYm0UFwr7CnMFsb9oXA LwdzOthd+IV3xuhAxLoG9uIGZbCoqXvJR7rg7hdote2gMozxZGuDInUBuzfz axVTGR7L0BZMqdEF79n5lPhqJWBFeQba++rBm2vTi9r+VAKnhUed2Bl68CnJ YYRVS4MTc36uUm3VA41MCxPWQRr8I6mZPEeTAXkW6i6txjQANcVao9UM+CuU zU74qgjuRlXD0WEMOOu49UgsWxG6pDmJB9MZkDg/eN5LO0VAhgy34+UMKCd+ apcOzYD19A7b/hYG7NLGWFE+M0CjvaJFURKBwhXc9lWNAiT4b21iaiGwPGsz L32ZAoQFxt5MWIjA93DxE16BPMx+/mpr2AoE7o1KasXz5OFUzNjIPA8EJjGV 34auyUG50RvnVj8EUyMGOGl6chDbXOFZfQCBbenmC2kXpkNGoLKfwVEEzKOK 7T1a0+Ew0lguw0ZwTPFFDOuiLOgJzOouXkSQejJ3XSiShQLmCiulmwiCs92J MzenQVUxuz78LgINQa/69wXT4JzKqr/GCsg+NzPuxkdUOBxHycouQtBgukm+ z44KdBxTvFaOIFr/Vv+lf2XATCmuqv85gue5B3e7bZeBssqKuffeIKhay0tS EUwF82bVTd3VCJLCqw+8PDkV0peHUO/9i2CeponXntlTQTVusI1Wh6Aw2uai TJE0RAxChfJnBDKd9Wn73aRBFw+98aQJQchj7c7yPinwM6xvZX5FUKxvMLeZ JQUv2zfscmxBMJ9eaVhmIgXaa4sv4d8ReCf//cynVhJSlrSU8NsQaNtZTC8J lQTGwIv0Qx0Iln7cT3uhJQnDF5xMOzoR9J4pYISXTYH2dyw5wy4EciqMtDd+ U8DgPJ2+qhvBJdYuyULVKVDf7nnE6hcCn/W4OREgAcOrLetpPQguWGter9SS gPq+YIdq0m2UKpu3ryhQv1Je5WAvAlVe6p+rIinASo6ehX4jkLpzMMnMmAJe MTkaeaRTmnaW0l6NEe9urU215CFwrl4/Pcx8jFh2/+pGgz4El/c9TfJQGyWu PU+2ZJO+32alIRs5QrSqrv4oIH0sGqwftg8TnqyVAXf7ye8tgCPV+lRA3Edu aMoAAl5Ers7SxQJiwdQz6qtJl6RHNvjeHSKU9rea1ZDeH1biUsMZJIaea2+n 8RF83CHzcgQNEueaiR570lvlz7sVXOUT13KzIu+Rnr34nXJm5gBxOe1krO8g AnvZsxEW9X2Eeo/Ri9Okt6QFrj3h30cwlqf6FZB2DflgrSzRR7xZV2crMYRA 4slmVvQ8HrEpjOYUSzqKUh38NL2XqLwRnrBJgGAn4zheWPqTUAmifNhH2uyI 2dduvZ/EgvP84ydJl3kcjXKI7Sb0bSMcuaTZpvL5nk5dROegzVOlYQRFq014 v2idhKbnsk83SJt0uviEKX0nnoXWX2GNkO8LMza5VtpE8La81b1ImnlnWZhH SBOxonDYJ4P07ZaHj44aNBGWztrOxaTtAAbcNRuJrNzsinbSVXGXIsob6wnT L14dJqPkvKFAWbfxA6GdyM/PIH166G4PtvsdYd/pZ7R0DMFr88uRbMcyohZb yy0nnX5Z97fuoptEx1bryz9JH7j14QrTNYE4zPEbGiI9vDIqZ7dLEGGygzg/ Sjon2XLxyAJnIqs688wY6YwfKmx6mSnxeVlv07jnttl1P6zXIPJ/dPWPm8UJ p4ZKjeIpAytGxh1U4l4ZZv0Fv3VqbMLrXeToBcvqcLWsQ6PjVnasKcnaWIM7 zmqccHyPVn58yGv8/umBCTcHUcOriypwU2/3sXHb+PAsc5gFeGpy04Tfqd6R Lc+8ib+e1jzhmGy+oSAgDY9o65zwNx/5WRTvk7hxzeCEG65YeZcHRuGf/nEk /0P6UGr7sMQ9PAjnSPZOzPewHr950O2FZzbET9gwtXvQ38QZ3/1a2Lexbo/8 YusluP+GsQmfv+73+4u9Mb4hlDKR9yk3L+HQkB5uPyycr1l0qoapqoHr8YT2 KnhfFJmkiLeL5rHebnuM3svgaQNCb8PsZ1JURrFeWWHe2+RVXVsUO7BtKkIv XCS1bm9LBaYwU2ht6nLXwcJibq2a0Lf9+Z2LDCq5h5SF5qfw25hKXG63nNDr 8zo22MTdxahThVZYrLa5XPAWq6UI3TJ74+a5HR8wu1Hhfkrr15VV5X3ApMeE NjzXsN0k/xXWLnLF9RijRYEPsVDR/YY8P8+V7EQsR+Sanzt8q36wuZki8xxv BzXrpHFjJud7G8NDdxzjYiKbvw/9Hdt3DBsQ5f/I+Zq+VIuDPRM5nnLKI3hN JvZE5I8prjt+SWVgTSLfv0LMWUVLw/pFTtEv34IFx2GSonwqY79+wqwlGFVk /raLo7eld3FlRD713V+qxcWfO7ne/BmzN13bgTsqyrNzijW/KfDGBkWm8Krt e6zDxH1JFsbrNjdFiffvvMWEIev3DyaYnGv6ypmx/LHJvMbkRikDD1fh45N9 lL6XNBqd90hC5Hj04NjSDmvupNfvMg0JCye4k+sDRxI3RDLo3DFR3r6Avd41 XzBs0jN2VmQz89aI/WBrX9cfPs5iu4R2bt79a7W435DKn3G6HRO7yPb1l+Iz +mI7bDvBzO/tezRprwWV57+bq4j3c0mHTlh0K4hdPCcxePumb1aTXhgYIbB8 qCrOM/DSC2zZh8RO3Cm54nGbgdhH/ohGeziaYvd//T1EXSEjdn7KEU31Bmmx wenzjKdXe8R9u/gbHQJ4/+//G+U9+Z9Wi9ih9CLfQPkmq/8AMGWvjQ== "]]}}}, Axes->True, AxesLabel->{ FormBox["\"N\"", TraditionalForm], FormBox["\"P\"", TraditionalForm]}, AxesOrigin->{0, 0}, Method->{}, PlotRange->{All, All}, PlotRangeClipping->True, PlotRangePadding->{Automatic, Automatic}]], "Output", CellChangeTimes->{3.662309130467573*^9}] }, Open ]] }, WindowSize->{1226, 914}, WindowMargins->{{25, Automatic}, {Automatic, 4}}, 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, 1719, 52, 124, "Input"], Cell[2301, 76, 405, 12, 33, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[2743, 93, 1103, 33, 52, "Input"], Cell[3849, 128, 2229, 68, 90, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[6115, 201, 184, 4, 31, "Input"], Cell[6302, 207, 241, 8, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[6580, 220, 184, 3, 31, "Input"], Cell[6767, 225, 139, 4, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[6943, 234, 346, 7, 31, "Input"], Cell[7292, 243, 354, 10, 33, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[7683, 258, 236, 4, 31, "Input"], Cell[7922, 264, 188, 5, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[8147, 274, 1872, 53, 72, "Input"], Cell[10022, 329, 3954, 72, 119, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[14013, 406, 1826, 53, 72, "Input"], Cell[15842, 461, 4292, 77, 187, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[20171, 543, 1793, 51, 52, "Input"], Cell[21967, 596, 3500, 64, 252, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[25504, 665, 167, 3, 31, "Input"], Cell[25674, 670, 11110, 191, 284, "Output"] }, Open ]] } ] *) (* End of internal cache information *)