(* Content-type: application/mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 6.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 145, 7] NotebookDataLength[ 38686, 920] NotebookOptionsPosition[ 37574, 878] NotebookOutlinePosition[ 37914, 893] CellTagsIndexPosition[ 37871, 890] WindowFrame->Normal ContainsDynamic->False*) (* Beginning of Notebook Content *) Notebook[{ Cell["\<\ This sets up the basic equation(s) for the numerical ODE solver\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Eq", "=", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{"Prey", "'"}], "[", "t", "]"}], "\[Equal]", RowBox[{"RPrey", " ", RowBox[{"Prey", "[", "t", "]"}], RowBox[{"(", RowBox[{"1", "-", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"Prey", "[", "t", "]"}], "/", "KPrey"}], ")"}], "^", "ThetaPrey"}], "+", RowBox[{"AlphaPrey", " ", RowBox[{ RowBox[{"Pred", "[", "t", "]"}], "/", "KPred"}]}]}], ")"}]}]}], ",", RowBox[{ RowBox[{ RowBox[{"Pred", "'"}], "[", "t", "]"}], "\[Equal]", RowBox[{"RPred", " ", RowBox[{"Pred", "[", "t", "]"}], RowBox[{"(", RowBox[{"1", "-", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"Pred", "[", "t", "]"}], "/", "KPred"}], ")"}], "^", "ThetaPred"}], "+", RowBox[{"AlphaPred", " ", RowBox[{ RowBox[{"Prey", "[", "t", "]"}], "/", "KPrey"}]}]}], ")"}]}]}]}], "}"}]}]], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ SuperscriptBox["Prey", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", RowBox[{"RPrey", " ", RowBox[{"Prey", "[", "t", "]"}], " ", RowBox[{"(", RowBox[{"1", "+", FractionBox[ RowBox[{"AlphaPrey", " ", RowBox[{"Pred", "[", "t", "]"}]}], "KPred"], "-", SuperscriptBox[ RowBox[{"(", FractionBox[ RowBox[{"Prey", "[", "t", "]"}], "KPrey"], ")"}], "ThetaPrey"]}], ")"}]}]}], ",", RowBox[{ RowBox[{ SuperscriptBox["Pred", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", RowBox[{"RPred", " ", RowBox[{"Pred", "[", "t", "]"}], " ", RowBox[{"(", RowBox[{"1", "-", SuperscriptBox[ RowBox[{"(", FractionBox[ RowBox[{"Pred", "[", "t", "]"}], "KPred"], ")"}], "ThetaPred"], "+", FractionBox[ RowBox[{"AlphaPred", " ", RowBox[{"Prey", "[", "t", "]"}]}], "KPrey"]}], ")"}]}]}]}], "}"}]], "Output", CellChangeTimes->{3.448368347393125*^9}] }, Open ]], Cell["This gives values to all the necessary constants", "Text"], Cell[BoxData[ RowBox[{ RowBox[{"KPrey", "=", "1000"}], ";", RowBox[{"KPred", "=", "100"}], ";", RowBox[{"RPrey", "=", "0.3"}], ";", RowBox[{"RPred", "=", "0.1"}], ";", RowBox[{"AlphaPred", "=", "0.8"}], ";", RowBox[{"AlphaPrey", "=", RowBox[{"-", "0.6"}]}], ";", RowBox[{"ThetaPrey", "=", "1.0"}], ";", RowBox[{"ThetaPred", "=", "0.5"}], ";", RowBox[{"InitPrey", "=", "5"}], ";", RowBox[{"InitPred", "=", "1"}], ";"}]], "Input"], Cell["This provides a numerical solution that can be graphed, etc.", "Text"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"solution", "=", RowBox[{ RowBox[{"NDSolve", "[", RowBox[{ RowBox[{"Join", "[", RowBox[{"Eq", ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"Pred", "[", "0", "]"}], "\[Equal]", "InitPred"}], ",", RowBox[{ RowBox[{"Prey", "[", "0", "]"}], "\[Equal]", "InitPrey"}]}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"Pred", ",", "Prey"}], "}"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "100"}], "}"}]}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}]], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"Pred", "\[Rule]", TagBox[ RowBox[{"InterpolatingFunction", "[", RowBox[{ RowBox[{"{", RowBox[{"{", RowBox[{"0.`", ",", "100.`"}], "}"}], "}"}], ",", "\<\"<>\"\>"}], "]"}], False, Editable->False]}], ",", RowBox[{"Prey", "\[Rule]", TagBox[ RowBox[{"InterpolatingFunction", "[", RowBox[{ RowBox[{"{", RowBox[{"{", RowBox[{"0.`", ",", "100.`"}], "}"}], "}"}], ",", "\<\"<>\"\>"}], "]"}], False, Editable->False]}]}], "}"}]], "Output", CellChangeTimes->{3.448368353736875*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Plot", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"Pred", "[", "t", "]"}], "/.", "solution"}], ",", RowBox[{ RowBox[{"Prey", "[", "t", "]"}], "/.", "solution"}]}], "}"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "100"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"RGBColor", "[", RowBox[{"0", ",", "0", ",", "0"}], "]"}], ",", RowBox[{"RGBColor", "[", RowBox[{"1", ",", "0", ",", "0"}], "]"}]}], "}"}]}]}], "]"}]], "Input"], Cell[BoxData[ GraphicsBox[{{}, {}, {RGBColor[0, 0, 0], LineBox[CompressedData[" 1:eJwV13c8lX8UB/ArM1SokEpWVrJX1v2kyHafq0TWvSUrQsoWGSkVDYrICIVU FJVEVChFS35Kw0hUlLIqye/rr/t6v57zjPM853u+50pvD2TunEej0bI5aLS5 X5m0mcJcyUaTSYGDGjTaKP0mn9KqAJ88ulWkpJWHwCi9MfxKkaLPVTrd/7JM s+govVk92JvH5y59wb8vSpfWjtIfHn1hMeH9lF7cliXz0mWU/sggoWbY+wNd cPufrXk1o/TWE93BQ97f6Yd/pZ6yi/pBt6rjtfvkPUuH5N0cM/Exusbs9gnR 0Vn6psyvAv3yY/SjUh7tskY0/G5o2XREZ4y+d/vtfZsf0PD8baj5X+YYfeOn oAfVHRzwDeExZ6eO0ftG3rHDJzjhfJIj24R3nC41czNnRpcPF3KXKidwTNDb FujO53ITgoaeg+M6iSn6C2uRzPZoIWTdcj1ZrDxF7zr8TT4rRwjKn+LNVhpO 0T9yX9yo+lYIPOjRVnOdov+dEYtzchHGju9TGu9zp+jK335PljuL4F7c7fNB 8r/oh9rrPm5xXIK1SwT1ykx+003TzBou2osjMN/kTWXGND3Oqqr0xW5xyDZO xQyWTNPruWVPzRwTR4qcwxXlO9N0g2ial8MTcYgkSX563TdN1/KtFZi1WAZ/ 70L+d+p/6XIbNLY6mkqAtipiOO7pXzrP1IpvXDorcKPm6FCm2D/6TrU/CzQW SKGzpvGO8C4abLqTi8UlpWCQ/ijvWBgN2slLjWlqUqj7HVuwOJEGrg/q/k8Z UpAtm6dmco6G86nerQHpUjgWfbFq4TMaPox0HCxbLo1u5njxRT0OOJdfnZVV lsGP2oR1IovmwVZp5w/RTXIoS+H2HezlxJ8j+2rfOcpB4LFwp80oJ0q+JSUV ecnhdWelef0/TnBVXVimeVAOIbRjgU3LuXDbZGi9bZMcmnzbDns5ckF+866T CRtWg4dH7xBHOxdoscFao3R5NOwUbxJu4caV3gMzN+zlEReqtU3sP264bDzZ EuMhDw/Fr2JrBrlxY/51F4FYeRyOU3aJ4eWBf/p4gny9PMoj7so8teDB69Kw DlcjBXRKl0zot/PgWkfMvkf6iqiIEO6u/ciLkiTP7X2Wiqg+vlxC/RcvcnWt 7ae3KSKg67LyLQE+pGSKK62NUYQYd5EJlxYfdrhef3u8UREbpL2K0+L5sOTj 0AYnKyUkG3+mp62ej9AxB5FBF2X4Pp04LhPND/8ig9lZf2XcH/x7gTrJj+1b pIfF9yvDQfDPl8wSftjd/NZkla+MWvG/Z4Ne8UMh6nDElY/KcMqQ/2ahLoCu eXd79gWsAfcDuaCGLwIwFFG+yhWrAobrU/HfIQvQzrHtUVuaCtRfhmfwHlsA 1o/D/afzVVAi8iRf88ICJD77LKZ0TwUbdnGs+NBF4lNL42y510Jm3vmK0/SF YAkoMU8fWYt3pfljJSKLkMilOKGQqQp/U86VLj5CEB3fuuhHiSrOny6VehMq hJL+ZKXbNao4UvL1iF+SENruDbrZdKsiqi6Tr+E8iT9wsTlIUg3CVd+K8V4I pTPymTVFajhUPtJd4SQM7+mT54qd1HFEWvvyMlcRpNsN349gqePQwDPeN7tE 0FBg9sXWRx2eb61vl0eJQNzil+5kmDrY1I9jh7NF8DDd7Zn5GXW8e3+0Ufat CJRUFTmHXqlj4ZrghT9Yi/GVVeej7KABQ57fouURSyBeJZr2b5sG8vlMOr4c WYKNvEHVL7ZroF48x5qeuwQ5l2XmRe3RgLziuzdK95fA5k9y9uMTGnAP4M3m W7gUl08x2/2faSA6IKmZu2QpApsHta7aasIycHln24goGlC+S8VREyUpitwt XGIQqg0qLHXXxJ+hYMn3y8VQcfW3SFGgJm6F15R7WYnhZ6bAz8wTmtBQWZEw UCKG0F1qFfGvNMGfdObqnV3iiBYKU3F000J2p5Pl1/kSeJJi6PlypxaiBJ0U vstJYCU3LYexWwv8MJ9eCgnU/z4kYB2rhROPpJTbQkmf6D/7hV6ghRtvRhQ/ DUggobq+RGlACycj1E0GW5cjZRvP6hl/bRRIGMTZVK6Eibrg5Oe92lBT1mrq b1+JH9wiLZ3R2vCvqhHNHV4Jp8qVvhVHtZGgzO8XrSiJ1Xw6l3eUa6OH9Tu/ v0ASDdU7dFq/akOyo3kwMnsVJoUazc746eD+RO0P1isplH1qFk3co4PdpXIG sZ+k4HbnyWBQpA6G+jauq5ySwgPvrsNWKTp4WDVTwJCQxon6720zpTrg/SzI H8uShoq/pKPnZx0sErPILR2Vxo6WKC8NH10IjxUvrpSShaeRvN1koC520uTa uLRl4VX5TKc2TBcSB3ZnBGyShU/2am6zZF28zzqevn+3LHYHPi10uqCL2i2r VuTXySJSTLY39qMummKF7n/3kMNJ71aXdpYevg+0e2lXr0b6270bTvnowWHb gh7Lx6uRQa1a4xSkh1tdovtCe1cj03Dvn95YPQydd9ZaulAeuYsksyZy9aCr b9bs5yOPspvBnSve6QF3E4MTpRXQyLuM2uWsD9rg5GhvsSI+FapavWbro/n+ XauVpK8IYuOGTX76UIuP0g59qwin8EAd2Sh9BDz4wnNgsRK+DzUt687Rx3bT +Od/45WwsnVPn1WPPsTTPW1yvZURcfRJiLL3OlxbEMkUtlCBhnDc6c97DTDb WSjzWVsdt/n+U2yJNsCey8Wvt5ipw5SmWluUaIAP/0RDnm9Rx+bv3R880g1g eXSF70SoOkLb9JQ6rxtgq2n4yh816rhz+HvtvZ8GEKz9ZVRnqgGLeR692UGG yFFaX9jhoQmPMWMVO38jdDynNDc+0oagCO1Zx14jOEVHdNl90MZt9XshLjFG YCm7F/tPaGPJbvNan1Qj9BZyt/RKkzoYsrNKqDCC9xHha9eidKD+3sPn1rgR qqXX+45p6oL26ECRTIwxFF/wZQlV6CE/t2nF1DETqCTWhocMGCDbVeW/LadN cPBo++vBWQOcljh1oirXBGE3NnkHShji6GkWz56rJuAces17y94Q4anT34ef mUBBKHq7Wq0hqP0a9/sW06E07BZ3I8MI89zP+bZn0ZGw7om0pLMJrqUadL9g Ar8WLpplC65HyG1jxjmb9ZBYzi3oM7wR9a3FPw+7mMJh3Qz9RK4lfDatq/N1 2gDlxBV59G12yC7RqN7ltxEvnDwVs7ooHB+rOxAUZgbWf1xVBRGbMVC2L9Rn jzn2SaxxvtzliPvmJzPj4zehcKhnTZyJM2j2b/rS4izAcaJE9rGpC2ZxgyM/ whJWdfe67v5zxbaNg61dMZaYzHAPuMnthqpN4unC8Za461quViboBh+7SPmE FEtsS4o2OiDhhmcuJtZe2ZYomGoOmtBxQ0Foc7pKnSUSeO+JTO5yw4byToUa mhWmvJZ9cuxwQ7LYlO2LZCsoUnnVhafdcfLhw/57R6zwvkSgkivHHbkRZyOu p1lBMz3mELvAHTe6jYvTz1ihVCS9fWG5Oz7mJf51vGiFs0oRbXqN7livKHKp u9kKExlpZg++uGNafy3fALc1gpo1K54YeIDny79zr+ZbQ/YSs7aM7gGR7Gda zQus4WUTVp+40QPKMyEeF5daY/5Lj0oVOw84Nd6+4bvaGryOxb42bA9UW1rs /LbRGi9d62SfHPRA4LYd96YSrFHnbni/o90DgnlfeowOWeO8/M3yEy89UNYf PHvgqDUmY+V4rbo88Mk/zkgwwxoiFotEL/d6wD323A2pC9Yw6nzQbjTuAZui /8osH1rjaoajlKkgC4rfrE9lC9qg3Ojg1w+qLDRrvqzsEbJBUae35pQ6C55h 256tXmqDOsXT9oJaLOTRfBdUrLRByt5pIXU9FkQXH0x+oGqDxavfPHGjs8Cl 3xA1wrDBz7oQIQd7FnoOaHvSM2xg5sTPVx7AwmrRhalBWTYQtb+eHhbIgl/Z 4M2CczaoHNBTWR/MwtiLswJcF2yw/I5wXPteFnjkOK4/vGGDEydvdT2PYmFt cxsHs8sGGYJbDaMPsxDJ753nudwWV6Kh413Iwt1cPMpYZQuHq50qk0XkebQk xpplbWE9mkclXGDhmEubuZKKLdKPrlPILGUhv1zr24iRLWRvp0iXXyX529GM w9xtIWBeabW/loXFp7LeHC6whX7zHo2Xz1ngf7X4+Z5iW3g/jhkTeckCTSy1 xaXUFp42LiZUBwvDZ+Or1lbaYmqB9p5HneR653elPm+wRUfR0K6SbhbCrxmb LuuxxX+P60OkB1h4+7y3tETSDmfq+dfxTLHwYrFr/gkZO1wWWPxG/BcLD7d0 no6UtwPHzVXSyr9ZqH7dmmCjagd9+S9FFtMspPZedxs1tkPD1rs5+2ZZwI8k YX03OzBTrZYd52WjWEg5oiXbDv6KDvltomyofnL+xZVvh7R9i71axdi4UXs4 zLTIDmWNC282ibPR7PV5X105ud9yy9mbEmwM3SnZc73ODp/UnlkckWRDxVfB P++DHR716jaMrGbjWqOcR5isPRqXl8uWabFhcGbz+yoFe0xa3oiP1Gbjvn+i 28819ljPx3HCUoeNDrGPLgHa9jgq1Ha3X5eNyYBCpx1m9jA+G5zEZ8DGOgkZ pr23Papciy4sARsNe1aZKVyyR59MVquJNRt/Jr4fX3HVHnk1UVxvibUjGt4K X7eHT8USvXAbNkpi2Xv/3raHoErn4cu2bBw/Wlj4vNUe+8eFZwQYbHhcUJgX /cUevLfXdGZuZiNT6Zdt8Dd7HNP2oslvYeNF+cMsr5/22OHTrH6N2Py6rwb1 xx6bwLzY4siGWM8W9cXcDIynLuIcdGJjdviMx9ASBvl/1aL+y5WNHNU47n5R Bg6NP5AOdWNDP8jn0jtxBr4fWCs/RrxnTH/qxQoGRGyC/L+4s/Hp9+vUOjkG SuwWcbez2GjnWVF/UpsB6wyLZX6ebORK5a8w3sxAzg931rQfGwtTNB6udWTA 5PcRN5ddbMSO3QuRdGIg8X1q2m1iVvNA66wLA9c2j6SH+rMh478m4t4OBor3 fNMfCGDjws0brzaFMKB8bEHSuSA2RKUtDujvYyA7/4ncD+LklNcqSmEMnHhb x7UxmA0ft78J/FEMNPQsODRIrMxpqtkWz0C+VmjF6hA2rti1pTJPMjCx0Fo2 dh8bkrfcDTakMzB/pGpZC3Ga9OiA1mkGlK7puiwMZSNoXMR46VkGVovJ3jtL rHnW6et/BQxsduDVvBRG6nGgf5NbJQN57cZrL0WwIW+/b8z2OgMXkvjkhojP 3OLJM6lmQOF2YJBcJBvhR5QmV9UwoNVYJZhNbKARWNTXwEC13YGdcVFs1EX/ ofk8YyCswHONSgypt6ufRuNfMBBlZieyg1i9/0XPuQ4GdolPGWURF1leanjZ xUBbTBqLcz8bKaKucfReBt6XLp9sJ35gaRG0rZ+BOEW6EC2WjZlobda+AQaO GWrvVScO6hfEpc8MlFE5K9KIHSvqZkV/MqBRUMS1MY7UY3/pd41xBoLkM+ID iR+Jnv5gM0nef7Ev4yyxUczuu/F/GNj7oOX9CLGM1arYb/Mo7B3M7D56gI3v /XH/mhZTqHh7YrYsng1FsYBvPUspqI696nxMvN3K+f20GIUz7TL8w8SdFRr1 GisojPD0JCsnsHEnpi/mnByFx6dkJvKIJyvaA27JU+DraWusJVb7eNvtpSKF A9Neo53E561OGc9fS86/IBEqkMjGIbGNM3t1KMjNpukGEN+zUh9J06NQOV65 OpF4OmbFu7J15HotGtFniXd/HL/TY0zhNR/n9ibiErHe8mk6hRtRFT9fE/da teWImpLrvQgZGSHeXFkcbWNOIVR9XHZJEhupH0/4e1tQqDdcwJYnbhHb7xpv RcFNY5xfn9hgv6PRLTsKO/v6LzoTS1vz/pXdSkH4CHPZaWLtVxteSDhT+D26 5vN5YnOPuBJhFwpZC9ONrxL7hfzePOtOwfGeWWgzcfSMzppJFoWVHGUuz4nT kvdwjGynMGbxoqWbuCr76+U3XuS4WZHud+Lm1YqJz32Iz/sY/CLuuuq57aEf hThn+RbaQTa+rCtQv+tPIeeu3zs+4pn773hu7KbQrHg/RohYyE7iXXkQhc2f 75SIEct0OV4v3EPitwy4SxKbDz/1OBFKYSHv993KxE6hgrqHwilsvfnykRqx H81SMDaSQojcRKU2cXRKUt++aAqvPjxTWEectuTeLf/95PtFflxjTHw+91/q jjgKasmX60FcpWi4c1s8hQWnS3s2EDdfCzOkEilM6TueNCfuMqoStjhIgc2h 0WExl0/z6KDJIQq1cTcvWBH/Zayt10mhEGHhzW9DvLDbN13lKAXBL3Wzc5be ecFPNpXCwKR+si2x9vc+SBynEE5zLJizecQqMeGTFHQyT9rN2YnTdYQ3nUKa vX/a3Pl+xzLv/8ug8Osmy9N6Lj+xV1kTZyj0bh9psyROLRAOGs4i9dGb9mQT cf4aO/P+bAr/mppYZsTXq1NWvDlHQbRA8IjpXH70lp/P8kh9cfVZ0+fye8T5 qKWAzLUH7hcazuXngLz6Qgr/laqe0pvL7130vupiCncWFkhozeXnU2NdfpHU +/pMA9W5/H5OSBeWUlDqODGmOJdftOavrEsUMrYKmMvO5ccT2H78MgVvt3jN lXP5Hb9UlHyV1GfY3hrRufwkhiL3V1Lg1vXoXjT3vVTZCv7VFNI/L/k5Vz/n b52b2X6TgghH9+BcfVWZvnnpXEOhkPN45OhcvTk6xG2qo9CwqC78PfGXnjRH k7sUNI3UBzqIZ/yeqOg0UojJOjTaSiwTa/ZapomCr8qHgeq59TE//uqyFgrb HnU+K5tbH6fqk4QekXizHNc84l0X9TT/PaHQ98TaKXnu/k+VjrzuoHCxZtsw k/i4bGPxn07y/J90b22Yiw9zalj+mkKgOENce+5+q5In3N5RMOPffWbx3Hra /dGjb4CC+P6W/Day3v3uRUVyDlH4oTu9spbYXHRxhtwX8v6+dmiVEP+tW9/q 9Y2C1P0raw7MxS/I0/46ScEjeZeL6lw8W9de8DeF65Eml8WJpavbfNdOUzjl IZUzj7jL9W9u4CxZD+FmDq9IPzMrd54/zstEolmibiixlPWS93/FmegatrG5 QPrl39xLv1YuZ2JBmkxlCvF/P00X01cy8f7XjdrdxKlngy0OSDMx283/QId4 +nP7NW5lJsaGB67cJf2589Dh5IWGTCzLK3C+T/r30eZ/6lJuTOwM+ul5jOwf CP6ksdKDiZfHOKd8iceXt2tKsJkIMOrcZE7sGnxOe8lOJn6IOyj8I/vR2hVG +nwBTIg3nRnyI24PjqCPRjOx6S3/A/1o0j9Wjts25DCh5vljYwnZH0+FfPZj vWXCwVQ1/DPZnyXMI4Md3jPxeU3Irhri8+IC4eY9TETMT9A+THytTiVJ5SMT j/ulLRWJn/EG5f36yoTisFqi5x5S7zmTL49PMyHLqRL3iswHhx5wGjdKOEB1 WD0+m8wT0UslF8k4O8AzaW9nkxcbpUrcrXt7HOCWYdxeTuYlnnGOZY3um3H+ h77TI3M2LqpX7s8d3IyLWab8LDIPKm0vs5sJ2II/NXwcF+TY6Lasubbl6xYw ynyYdUvJPFYYuZYe4ojM868Ss3nY+HWdMpL56QijNL5tuWQenr21p/rg7q2w yXXxFPjGgrW6uVn0xFYw2517OvpYaHFT890V4oTglKimgtcsJM+0Fn6ddoLq S4UxGpnnN6cvm5gm3rR/aHb2GQtSa3zMBf86IVXgiuCcbztxf1YhPvxkk/wM 8fB1uupuYuvBauffxAy/qlujxCPh0Y0/iCX+O9c2MeMEgY5r6X3EnwKGJXn+ OcGuetH5XuLrXIZBosT1YcFXe4htNLtE9Ihp+Rta3xPHHhNxCidWi5GhdRN/ 3JDc94dYYUer/0viijevtARmneDusC/yBXFMsFzScmLIrT70nFg0r1HBiPij 96nCp8SWf6b9Y4ib9pW9eTx3/IRV/THiQL/ooVbiPoWsRbnE8XSHyUfE0Vt0 r9UTN6csEnlIbDGcyPmU+BNtYlUL8dKEl5s/EDO2fljbPHf+MpkL34mnj7Ub NhFfqQiamiX+UHTP8gHx/zHgbw4= "]]}, {RGBColor[1, 0, 0], LineBox[CompressedData[" 1:eJwV2Xc81H8cB3D5WZWkQiUrK0my53EviUT2LaWIcmfLSkUJGXdISEiolFJK 0iKiaBCVRJKSJKWSBmX1+/TP3eP58B2f73t9vvew3CvYZQc/Hx+fNvn49614 aPpUoVyDueaMRgsfnwSui6yUD+QUUf90vxFv05RAQ9TFEjXOJeon+aS5Ilsl cE9rJ1uIc5uqHK+rca9WAg9S221+sR9TB+7c/ZgYJ4mHJvE3P7PfUO8k++rw KS7GecqiffqX+qhmzietFrEWY7Rne1F42Vvqp6Net0+lL8Y+aZH+78XvqGXi V7ZbTC/G0aMO7JG0D9T03ys2H+pbgubDPTuH2CNUk7l7t+6rl0Zybs7DLO1v 1Nw+g+YfM9LQ4FA36V76RuWf8b7+x2wZQoUP7w0tG6UWHTTUvVW7DDNWenXf in9Q6TH37qNRBkvu7rH8kjZO5cbJl0V+lYNtrbDDIPsvtSCtMLOzbzm0/3r9 kvr2l6qz0Yz7Z54iUhU82pQofFhatVsjzVQR4V7VEbRGPjQYz22dk6uIdYMh jVc7ZuGZtqPAPoYS+r/0bov69R8UdAxs8UkZnem8qxnGAti4ap7FTzkVNK8x nn0uRgDuYuFYRVNB5c7sim5BQfQNBKqvqFfBgV+2f02khNC/c6R4T6EqFKav F0wbiICqnTl1LUINEoU7vknuFcGrLu1bBlfUIEJdtE7ztgjapCxbHUbVMBIb NLx1/Wyc+VNy1j54JeoEVEzqmXNwpbDzzMVQdWwRzeyKixKFetkGD/tcDbTO M5gtsEUcWVz/o8uvaqHdbmFuW7Q4onbx7Je/0sKLlK+qeQXiWFqzkLLnP20M CJau03wljqFdfztMXLUxNb04lrV5AWbpBUSJjmtD/eufsQtuC6FVfKvylJMu tDQ6D+7avRDdezqeHInRhYFfpcTavIU4cihQ6WuZLtYO+ul0vViI/Jbqhx1C enB73RM4i7UI1dpRvw/f1UNyW+0AnSGBMT/jD7nOBkgXzQ9TiJSAUMsRe7eD Bsi2jeQfPiIBdtzphPSbBii6p6kQ+1wCmZpdKzcrGeJaXdHmMpokRmPKb+ZO GeL9xQPtMy5SaNWPupR2zxhrD1nVlzougZ2X77ve22aIta061x60BE32xeYf +M1RJ6iUNZ22BO6njTbEWZvDJJrPx/XREsT83asT+tgcur41c//aLIXOYu+m TUNUKFtqMxlrpZHqZWlSvswCXjNFOLBNGtcuxS0TXmuB4pti6hdipTHEPBZ7 hWMBGa3PU7NuS2NHSX2o0HULSMqWniw3XYae/EqzUuZaCI3LfBXQl4Gh0acO n3OWWFfJ61pDkwE9VCDL8Lkl4gIn6jeFyYBvv2H0jlnrMPOuK6visgy+02XP 9Gxah/EnmcbumrLQHXDXtVpohaEykYNXVsiBHs9svpBljR1rJuZpz1PA7Zai F2KDG7CxJ+n0EjkFrLpw7+BDJVvoJUma8a1RwOCG8TAxL1sIvNEKeOykgF7Z DmvTPlucTGc3B2YrwOKRuKb5Ozu8+dKRWLZsOWot+Us6BRxwL89LLlNjOU74 C7+osnPARatvV3ebLYfdW9mjmlkO2Hd87nsbj+W4fNJU54+SI2Ts11p+OLkc RXt15A+qOsHtwqW/SuqK2BcYkCPq4QywzI7OJX1oU13+/ViSM9QEWjR/2CnC aFUU4iqc8Xvz4JY7gYoIHm+Zw/jPBTlzZW95XlZEfBn7qXq5C575caOOGykh f0/9m3kLabBfuWNUar0yGtJo9d+VmZjgRdT0MpSRQQv9xs9g4uzXgwdLfJQx Yrb+UHASEwJVZ5bqJCrj1ejO+Z7DTFSbD1nYNymDpm3Xef46C6o0/8x4SxVw 54mJuHttInnZqfuNqooIWeFNYf9txcW3B6avOaqiw+1j719sxeZ1mfdjPFTx 7YHqm/F9W3Ft9pXNc/eroresXkJjaisCsn/Gq9apgh7Y5nt50gPd53Z1uFNW gDL99Pl94W2o7IiJeGikBqGP/Jol1t44e3C7V/8GNRhYTnlVBXuj0MDOcXKT GuLf3Kmdk+cNbu6Slatj1JCebZLN/ewNb/crrzIa1GA97HO4J2c7JAaGLFm2 ZG6YZeoN/d6ByB+uCz9sVsdO997EsB8cBJSY/P0boA7B4hEOV8UXXvTln5fs U8d8HbEjn5i+cLj+tcm2WB3Kgum3Emt9sWJvyu6LA+ooref0rU71wwv+230R gatQfWT+Im3DAJguVL8ksF8DT2qrLU3vBqNt1qaHrYc0sOB+YCp1Ohieoynv coo10B6uZnR3WQgSnnxcvPKOBqgGu/YEMkLQln4u1l5wNdLjTCvsWkPgOXel Sw5vNZpfDuTMb9iJBAG1XytyNaE2N+bP5tthYE9mHj/N0oLtoE7rTZEoBN/7 oHvJXge/Qz1XDlnHoh4X/DUYOnAdOmlGd4uFeE3IqXNbdTDblbuk1j8WFZf+ LCwJ1sHCRb/1ozNi8T137vfcwzpgJsnnXe6ORaT/moq45zqwnp840+R3ANHi uzQYW3RBa/ykuygxDo+4ptuf7dDFDYp0vf/ROMgK8hU4Bemi22iVzY2zcaj7 kzzXbr8uwmV4yRotceB7l/+JekIXl6aiI6PnxyP+at3Zle918TL3JCPrSDy4 m4RUpgP0EO1ToRSXkwBzLdGxj+F6oN8qm8d3OgGjggvvd0broZUmrxN8JQGs y7K+Fal6MGvhv63wJAEqIvrl3hf0wM5XcGgTOYj6q976zcN6OH/au9Ji10GM iTdYHfXTh4TL+nPCNokoG7wnlRCqD+7TTzGJtERsufXoQ8gefXQNmHj+9ExE I/tFii1XH0HOH5cd352Iw3UjrdPn9FGb+1RneVkiNALkGNs/6mPv3nzrr8JJ 8L6/10ebY4CfG1Pf7b6WhO0UVYexYAMMRMgq59cnwefyE/2aXQZ49WtH3rnm JHCOqQhaJRlg3kNVlZzXSQgKfnyKdcYANR4XOk8LJmPPYqW3+wcM8Hl16LUT TsmITm19sO6zAd5XnTe75paMmFlRFbN/GqD+1aw11V7JiB1+tD/rP0O4hX7f kRGejMTbkXKlioZI8fetizyajEx28+Y2T0PkCoxt0nqZjOxX4ZZZHEOcDbuz xbw/GUec5VexQgxR7zFqQvmUjFzT8Im3+w1h8+rBWeGJZBTOl8v7VWgI+1PH nV4vSUHZ9Z2dMr2GkM8q78l0TUGD8FJnfzcj+F6Ybt5em4LBU5q23duMUNuV riDUmAJRrLNc72eEMKb3/JzmFLCigvWV9hqh0m/N1IGuFIwMNS3tKTCC/P67 H3nfUiDbHNpv22cERaWaJ5HyXKz1Se65MWQEXs2SE7NVuODwF3aojhrB6+57 6VR1LqpMHt7j5zfG5C95F099LuzOy5+vVjImfex8Z70tF7tTH4Wps41xU/aI peJOLgrV+gNyg41x5I208I0ILhobx3cIRRnjo++mTPM9XIhPK7L6k43h0iJp IxfPRWngbkp+mTGucwbV1LO5eDT7kL7IFWN0SZ5StMnl4vvpEs3IGmOohRqG Mwu4MH/9WMHlkTFqVo352JZw8dxBVXDOV2Po2urFx1/hYuKT6cyuMWMotprp aV7nQiHJefz9jDFEre63NldzEXA7+mPDfBPQ66v5XjZw8Z/Ws9bd2iaYk7HW 2eUxF9oLYnM+hpug8Gvzr5sfuKgW6VK7H22Cl5/kuZs/kXjxadaUJJigUXDi 0ehnLmgjPW88sk3w3691R6dGuYhsNVzZecUEwizhe6mTXPxtTK+5UmOCkOkx +vNpLpJvvbc/fNcE/Kk7+Rfx8ZB3Pit04zMTmGRdHwwU4OFWykjNne8mUMtg 0U6J8mB1wNqheIJcrykEOWI8tEUd74vhN4UFwqv2ifPwhm0naLzQFB5GX3y1 JHjgty51uKRtiqDTUXFLlvHAM5vpSzU2xWbZ8yo3ZXiQ0KeH+VmYgjdpWeMo x4OK8n+5Ks6mmBPy9+vW5TzY8Hu8PRZiiuevJ54KreDh6Z+rYbujTKGpWbeN pcbDplFRIWasKcK/ZlsXr+TBv69afWGGKWTWP8MiDR5+di2oHTlqCjXFwSDz 1TxEP+Y4thaZkn3CS8lTk4f0Oqnw5EumWD3zzSlJi4fF1wKFfK6bAnvK+1O1 eSgub8y1vG2KTzvVNVJ0eLhcEFo702YKptutxT56PJhkP3R81WmK/Pygdmt9 Hu7wFPpvvjaFlZpfhJwBDx172oQivppi7uXG4kpDHraEquS5jJlC5Mt2arAR D4O+0au0ZkzBV9L+XdGYh3GWutPwPArOCOQ1BpvwEOt0oP+BJAV6Izp/REx5 ELF5EX5GloL9nTLBecTShol521ZTME/W266YwsNJzd5VVH0KtGv5HkiZ8bBK Va9OxoyCZz0hlQnEVbI8p4l1FHgv7ZD8REyR7O/v2kiBi5TdtLU5D02ixhFX aRSIh42FHSN2EMgQznIn138ylvCBuHNyMC9kOwXVrQGmq6g8ePww03AIoCBR dPdJH2LRhXxPOsIpWLqfU3GcuFrrTtjmGAq6p5QkW4jZjgmL+w9SkB/l3PuV WCLIuoaTToFYcprkbPDQkCriMZJDwYf68OalxEHnm/kjiyjovHFKZDnxsubU M1OlFKQKXn8lR/xgyME2voKsh2Lnuog4QnjB19k3KVj0cypqhlxfUfXZ4YwG CtY3cDb1ET9ed0R/cTMFwiarBG4SR3szu4+3UyD0sy8umXhl3NIY5R4Kfiao vnMg7izuUTj/joL4uFxtUeL428cbtT9TsD3u2Z56Eg+t1x6cGz8p2CAd3uJP 3Du1XJQ6TUH/GyktMWLusoFLTYJm6OewrpaSeBuanHHdKGaGP7bP2UbEAyzO eLuUGX66z6XVk/wd3qV+zE3eDDYNiQepxOY5n837VphhAWeH8DWS79xnIYlf jMyQu/SuWgqpD6vvOurhFmYoCo15/Z7Uz3fxX60TG8xwsWzBAhPijQ67JUU2 myGCZr70Iak/vocHShRjzFDaOfg5jtRv+QdLm3MHzTAcfLG+QJf0k5DQ5zXp Zti03MzvIqn3KkuurlmRGRx9Vf1vkv7g1GXeYTaYISGV6RZM+kuyl+bz+qEZ XPaOztuwitT/pNScHe1m+LiG/4u0Og8yxsecQ9+ZoWNOrHMp6d8nV071pQqa g+l3rdFbkYeY9h0Ji8TMoeT3uPevAg/qoyvU8qXMsT/EMThbnocEzQvBpSvM YXc5WrGUzAujc1UzDRvMMfye9jJMivRnYZPMeJo5Zs+qb6wU4eGYu0YXPccc 4y5PY98J8ZAjnXW4qtAc1CH3pHmCPKTmeAqFXjJH4WfFbPtZPESlT458fmIO gRb19H1/uHDep323fxEVBi/O/Xd+iIuNlNwYCxkqVOaXKYYPcmEzMWNYrEyF 5/NrwoYDZP5HPrqwRZ+Ks+uNUfaGC/Ugds4LBhWydje01ndywb/1uG9bHhVd z6fMtt3hYnqZgPLqk1RUiq4Ubb7Nxe9uv9e8MiqCNoTMX13LxVe6keuGGirC /hi1vSX7RY/9M0rjKyrEDRb1CJVz8bJQQKR3PrBbb7hkXQ4XlekmPe0uwEM/ 57U3t3Ox7ZiYYDUDkOEvEirfRva3s/2aJzYBIToiOwq2chHUwI0L9gKWWarx B7HIen++XCkaClj11Br32nFx0m1vlNVhoDR3yHhEm4tM5VuSNx4DlQn6rdMT KQirNnM6vtECs6eLJvMPpKBq033+bicLlP9WW2IWk4KxCaerEnQLjPtuUO2J SkGU6Xbp1C0WWNrepDY3JAXRt7gDe4IsIPxEj2K6NQWJdV1RrAwL7KHZWr8x TkHunZ0nFnVY4MiHJYtoX5NR13z6e8rmtTjkGmP+bGMyOOuNa31Z5Hek44aO JXxJOHZW+6q/3zq8Vjc4uibyIDJ+1B4I2WUF1nqxjbPexuN9WUQkJ9Qac5Pn X3jiHYe71pm5cXHrcQHKFfOfxILP8WX/oVgbiDTcM6x/tg9/cW1W8e4NeOu+ YoWCdTSSFo/btyfZYv/+Yp3H9bsRvMn7zni8HX5+Fy/YcWYX+g7obace2YjV Cce695VEYFFW3suUE/aI+ha9ynMkDKfF1XffP+aA4S7FOb6HQlEfKm+14rwj 6SeR/mqbnShUKJYxozlhsqLlR7tVCGqjJ/g4T5zwI/nKtqWvgjB2afBbXLsT 7huM7vN8HgStd+19xzucYOV+P/xCWxBKNpyvf/bCCU2Xt71xaQgCV8o9lvrW CdFCjtu7zwSBUVH7V+q7E3a6irc2hwZh5F3sTNMiZwxus5fWEAvCcjvhKSWm M1IC3sW7uARC77llu7SbMzrOb/GXtg+EtUfs2QWbnWGR2Df+fn0g/ML+0P5u dYbTrnvjSWaBqDo2XP7Sxxl0g1i3cTVy/OfHHocjnXHbR/HFAT5yfFru3Zkj zggbjXxjWBmAqscred0dznjUdk6lVj4AGUoNpyc6nXFHZwP/rGUB8N/Fql/W 7Yzm+T9zbKQCoCif9GtLrzNanoa/6RcNwKGgAY/+986YrNt60XLCH37zivSG x5whkHLNvuq5PxTsJF5PLXGBx0TqW4tD/ki9N6OlsMUFscuHJ7JF/YGdg9qy Hi7oKviuMi7ij5/L2nSkt7kgXMV/sbugP9x3HteT2OGCJ4Pj1pozflgtQzES CXTBj8Yr536M+KFt527qt2gXMAZfLZvV4Qdx2Z/29QUu0Dwj+2FZoR+ywj76 eb5ygWPEpdWTxn6Qtt6z0/W1C27MlypcYOiHk0vmRln3uUAgX9ZdXc8PlbUa BzUGyP3pNTVemn54IhxS9HvYBXc3DLz7qugHsYKxZxmTLrhgFF8eJuqH5Mb/ zBqkXcHxpXervPVFtKTcfEU3V3zMPW46mOGL1suzN1E3u8JS6vcBl3RfyDn8 KnHf4op1wcLyDTxf3E58ZHJ0myuUM+YxTyf6gv/3Hh9RP1dk+UxnpUT7Iull V+3YHldsj87mq+f4IrPocEBLgSuKLp+rfrnWF+dWCjaH97ni7Z9HvVHTHLTL WEwG9ruiMWuPa+YkB1PzYzTYA65oivE9WP6HA4exn+luQ67o/NW19OMvDn7e GXA1/+aKDx4ZH8K/cgD3xl5hPhre+Bqsev2Ggxfp8aP58jSg2O3vyF0O+ONu K2Yvp2FKYPG4+R0OVkVMuqYp0XDXhcI9VM/Bvs1h1/avoOHzwNA+/VoOlNW8 925fQ8PChNHE7GscBDesFdSk0tBXMW1SeI5D9vlZSxu20rD49cfo7gwOflfk pg940rDspTG/ySEOhgPXCIh401B/g2tdkMbBkw/uIw5sGmZu2Qv7cDnI77ne 9CqYBp89w+vFEjjQvBsY+juWBrf0ww8rdnHAyOxuWXOKhk/tlX/1PTmwcQyx cD1NQwd9WuvRVg5MRYWvR5bSEBRkPO29hQOFRL2Ttedp2LhG3Ct3Ewefdh+K squi4dSRKTMJOnm+bVYq7CYaks9ssGLaclCqdXlf4Qca/B0i1P4z4CBIcvTP 9Y80DHvNcXyix4HehFbE02EacpYrvi3U5aD+boW/wDcadopNWUKbxJdRwfL7 TYNWre3HLA0ORGIu6RjMpqNR9VV3jBIHnIflg63qdAgYix6jLyTPd/HLtg8a dKwSuM9wWEDyl7m6l28NHfSf3z1sxDmIdS9/pqtLh8TsRUyIcZA7cqE+z5SO Z76nb1PncPBA8kL+jo10rNQvvLKXn4OVXmUO04F0PHIz/vH2BxsPjptUfA+m o+D4+WMLiH26mxcM7aTDPOGG/drvbJQ4DXe0R9AxtXGm5cw3NhSoGpvPxtDh WVIrF/2FDSmZi2zXNDqS1NVTmB/Y4O+8HHvuAh0amRXtj3rYOLFgbX/hRToa /ATZMsSwb7fMrqCjtL5jJOAlGzGN34X2V9FR/tr8lXg3G78rddNot8h6ahrn be9k4+uhq3l/W+jgcj6ZGz5lo2fDzUr6MB31Os6c1PtsrHUsYUV8ocMq+P23 z/fYKKMdmskeoUP+HfOAA/Eejx22HT9IPLmPr0k0sSEdsaDfZYqOawntW8vv sMEq9l3gNI8Bw0sBf8Tq2Kg/TbsePJ8Bfdp4yf5aNtTOU7ccWsCAwRFsHL1F 1ndV8lybJANBAsK8rho2jrY0wF6OgXfrNaou3WSjc2xpiO0aBtwGZOemX2XD fEpAyk+bgWLK6cVixGdmfatJ0WXg6JGWRYeq2IgUvSf80JCBO15furKusCGp GFq0HgwYSzvsOXuZDZp9c9s6ZwbKunv2T5ez0X5qz2pqGDn/TOJt51I2FvJd u2kSwcDIjucjP8+w4eo+amWwi4FvD4p6c4k7JDhbV+9l4LLxnpH+02x0HaQf WhbPAFdgdu/+EjZesbVGxzMZcNA5u+PVCTZkGv1jfmQzoHDaW4NHvEWhdPZI Drl/9ftSE+LXL2QVB/MZOJ0mXnOsmI2+DXNdO04ykLmgca1PERvvVw1WVVQy IJhwa4vMcTZUk5dbXKhiQP1Y9prOAlJfA+6tpdcYkKDtK88gHip49r6QvDD9 9nM7L0z8aV6DVNodBlj+o0XT+ST/I8eifNsZZK5q/xTJY0NzY5fAjg4G7G/1 WLfmshF8duFhz05yfdv1ZpnE3zxTzjFfknzo7FokR/zjaeRLq34G9oo+YJod Jfm74kxR/M7AxwMFHXlH2BA1Mmhb9JMB+jZdRV9ihVvSnoJjDCiVtAgaE9s0 vosf+sPAUiHttJfZbOR1RLRcnMXEjLRyhjKxya+8TaYLmOhJcJnVkcmGw+59 nzVIK36Mroq6QOw17bVPTpIJoxMv9x0k5gpqnJy1lIknCiX6xsQ9knUf7ysw sbMnFqWH2Yg26I+iaTFx0UTE53gGGxnV9+ZY6zDB6N6WtY+4xPx8gaEeE3KB a208iR9Zh9dLGzGxXvWhtQqxLFNI5C2ViYQ6C+7VQ2zo9HzKa7dg4tr7FR/y ia09Hq9qtGQiymfsfixxEDvXsXQ9E+9UYxY5ENftUj8a6MiEl8zlpq/ppN4m xVZ6ODOxlj9C/wXx4P4f1U6uTBhfuapyh1gs+dZrXSYTNcGSiUeJt+Y6qE5s ZWKesqqpDXGojO6NYU8mruiePaVPnFi82LbXiwnhWbtSlIgvnu0LrPdhYqrV q4WP+O+N0KuJQUxo8NvJNqSR+hSq5J0MYWLiXab6ZeJq2ui2ulAmtjxPu3+C OPxbsNhYJBMnH2ddTCB2Nr80sGA3E2835o9HEmumfq1evZfE945tmy/x0IpA 9o79TFD3Fys6EzdFXDA7cICJlyYWRVbEJ+8OLzoez8SBrrhSE+ItHn71HUkk HwVvWKrEJuXncr6lMHHYeIxPlnjJ5FCAaCoTJsUlOhLE7TnspesymCi4cltL gPjSwJkRj0wmKjtWTE2lspGqM9i0N5uJDWflHMeI/WJVCo7mMJExtHzlN+L1 bdtDr+SSfDleSvlErCxTYvM4n4lPR9eFvCfm93snN1zAhNtej1d9xG+uK/4S KmJCryqm9RXxLUGvFsUTTJydkrHsJs5zPXHC/BQT78VqKJ3EkSf6dm06zcSS 089vPiN2HZF3iCwl8dn96eZTYi0zD+XMc0zIW++iPCGexyucKD/PhFPTHzwm /vSi98nDciZ0Lo/cayN+oCpb+v4SE7/udrb+8+lw95hZlUzcv65B/3d83J1j rrJVTHRZB7v/u56neM9K42tMuK6SeffvfmZbpfnoN5jwl7889G890hfcOkOq mdjYez7g33rH/+ReSL3FhKo8N+jf83SsfxF3to4JtaKHX/497+Uji90a65nY v3jo4794pL9jrOm7w8Tmtd5e/+Llr50jONXIxHD3Dca/eNrsf96z+D6pn1yf eyPEqq0SlboPmXhj+KniF/F/y2jJji1MnAhtkPyXrz5O1lb/Vibm2/r95Cf5 rLvWrpf0mAnxUP+Nc4ijXJzf1j1jYq7HRR9pYnpxxvWXz5nIzpurpkSs8/Vx 2lgXqS/H354axF9SHEw0XzFRK33dEMTNXWnitq9JfktKOmyJz6q0Du7oY+LH EdYHOrFXg13W8QEmqo/aBgQSU+fzfG8OMiEiV9K7m1hmSzP1+RATLmMTdxKJ O3/bfBb9wsSfpK4P//rDTsvaKvoXE8HHH6X1EvM9fz4rfJzY/Fn0Z+Jru3fU +f9h4ov93dZJYsW7CYabp8n6mZ+zlpF+nKDfXWkiyMKDYJNSD+KKCZdBbWEW pl9vP7aT2Keo/+TK2SyITh7+Gv9vXgzNklk6j4XBF7EV54jL9kLstwQLzQdT cv8Qeyo8aR6RYmH+tRN7xMi8kWrySPqwhAXu0+c3lIjjxPbzdcmw4OArZ+NI vOlE7feryiw4Rs4+W0Y8+75JV5geC1ovZsuEknlX59+c5W/Agpr35uxU4nDx TU7eRmQ9jEmvUuI+t90PXShk/QuvzHpNfOPz9RrtdSzYrH1p4kjmLWehfvGI Kwvil0LtLMg8vr9lja9/GAuv8PDYNjLfrR8t9UuLYGHDo9wNKcRNJgL+l3ax cD3FyOky8d3F3QE/9rIQ+n2tDz/ZL24/jQvZm8DClv7Jc2eJr1t17uIdYeFN WIOuINl/SlfHJP77v5+BXNPOSbKfqRawk1pvsrDCrzJoNdnvTs9xSR6pYeGY b951D+JTH1S5uvUsMDQ7VzURF514klbzgJy/Lycr6xh5/5BUPtLSzUJgxSFR c7LfJk03nxqeZEFHxUn4PtmfadlLf00SXzj1/f4DYoVVHGvRKRYke04lNxNX swQ/ahAvCdWZ00b8+QpVM4g4X5NP7Dmxk1/VjW/EFo5blrwjlu463vprmoXy gqVGf4kHAz/LCc2wIP/r8AQfeT+4ImAaIkWcP7qklp94o86LhYbEYh0ua4WI 96ctZEURl2hE2s8jHrBM6p8gpuTv8VlGXPHyue7cvyz0DrxaKUscs1P54DLi RQ/tv8gRSxU1rKAQl2zeHK5IvGFiMiCGWOyoSvzKf38/bFuXRrywYMp6FXH/ irz5hcRh+v1zVhNH0w0q64gvKHZnahHbfE747zGxfsQnhg6xZPwz2hvifW5z lun9O3+p4pkR4jdtRm/0iS9WhIz/Jb7cEXHKkPh/O0nNew== "]]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesOrigin->{0, 0}, PlotRange->{{0, 100}, {0., 872.0161297598631}}, PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Scaled[0.02]}]], "Output", CellChangeTimes->{3.448368355424375*^9}] }, Open ]], Cell["\<\ This command will export the preceding graph into a JPEG file that can be \ read by Word/Excel/etc.\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Export", "[", RowBox[{"\"\\"", ",", "%"}], "]"}]], "Input"], Cell[BoxData["\<\"graph.jpg\"\>"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{"i", ",", RowBox[{ RowBox[{"Pred", "[", "i", "]"}], "/.", "solution"}], ",", RowBox[{ RowBox[{"Prey", "[", "i", "]"}], "/.", "solution"}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "0", ",", "100"}], "}"}]}], "]"}]], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "1.`", ",", "5.`"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "1.0944330403234594`", ",", "6.724839651682503`"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "1.1974113857237125`", ",", "9.037649313691434`"}], "}"}], ",", RowBox[{"{", RowBox[{"3", ",", "1.3097171652323845`", ",", "12.133752318696736`"}], "}"}], ",", RowBox[{"{", RowBox[{"4", ",", "1.4322337915891048`", ",", "16.269532968224496`"}], "}"}], ",", RowBox[{"{", RowBox[{"5", ",", "1.5659740217812157`", ",", "21.778627126285734`"}], "}"}], ",", RowBox[{"{", RowBox[{"6", ",", "1.712119162228486`", ",", "29.09022025223283`"}], "}"}], ",", RowBox[{"{", RowBox[{"7", ",", "1.8720730848203775`", ",", "38.747621560449616`"}], "}"}], ",", RowBox[{"{", RowBox[{"8", ",", "2.04753649607906`", ",", "51.423905147334075`"}], "}"}], ",", RowBox[{"{", RowBox[{"9", ",", "2.2406071917311636`", ",", "67.92817382927544`"}], "}"}], ",", RowBox[{"{", RowBox[{"10", ",", "2.4539111575911994`", ",", "89.1921080722443`"}], "}"}], ",", RowBox[{"{", RowBox[{"11", ",", "2.690770107337426`", ",", "116.22300828772778`"}], "}"}], ",", RowBox[{"{", RowBox[{"12", ",", "2.9554046506416083`", ",", "150.0066302313234`"}], "}"}], ",", RowBox[{"{", RowBox[{"13", ",", "3.253165909612729`", ",", "191.34737342268326`"}], "}"}], ",", RowBox[{"{", RowBox[{"14", ",", "3.590776359449059`", ",", "240.64856663589492`"}], "}"}], ",", RowBox[{"{", RowBox[{"15", ",", "3.9765466861610084`", ",", "297.66599733511407`"}], "}"}], ",", RowBox[{"{", RowBox[{"16", ",", "4.420528087814391`", ",", "361.3050466772271`"}], "}"}], ",", RowBox[{"{", RowBox[{"17", ",", "4.934549555850298`", ",", "429.56281214026717`"}], "}"}], ",", RowBox[{"{", RowBox[{"18", ",", "5.532126995274104`", ",", "499.6871639138506`"}], "}"}], ",", RowBox[{"{", RowBox[{"19", ",", "6.228239907751799`", ",", "568.5597519548836`"}], "}"}], ",", RowBox[{"{", RowBox[{"20", ",", "7.039022411364166`", ",", "633.1969346807349`"}], "}"}], ",", RowBox[{"{", RowBox[{"21", ",", "7.981424006266057`", ",", "691.2038739059901`"}], "}"}], ",", RowBox[{"{", RowBox[{"22", ",", "9.07288731574801`", ",", "741.0413966265966`"}], "}"}], ",", RowBox[{"{", RowBox[{"23", ",", "10.331059612487564`", ",", "782.0602233839204`"}], "}"}], ",", RowBox[{"{", RowBox[{"24", ",", "11.773522172365478`", ",", "814.3517132141089`"}], "}"}], ",", RowBox[{"{", RowBox[{"25", ",", "13.417502581589988`", ",", "838.5122409388217`"}], "}"}], ",", RowBox[{"{", RowBox[{"26", ",", "15.27953226465809`", ",", "855.4075223547995`"}], "}"}], ",", RowBox[{"{", RowBox[{"27", ",", "17.37502266795725`", ",", "865.9868642760865`"}], "}"}], ",", RowBox[{"{", RowBox[{"28", ",", "19.717745448736785`", ",", "871.1616619665604`"}], "}"}], ",", RowBox[{"{", RowBox[{"29", ",", "22.319230661295673`", ",", "871.7412385463069`"}], "}"}], ",", RowBox[{"{", RowBox[{"30", ",", "25.18809987854585`", ",", "868.4098597319843`"}], "}"}], ",", RowBox[{"{", RowBox[{"31", ",", "28.329353189345525`", ",", "861.7288539987572`"}], "}"}], ",", RowBox[{"{", RowBox[{"32", ",", "31.743698990807992`", ",", "852.1529582758902`"}], "}"}], ",", RowBox[{"{", RowBox[{"33", ",", "35.42690872747931`", ",", "840.0509687719253`"}], "}"}], ",", RowBox[{"{", RowBox[{"34", ",", "39.36932529819863`", ",", "825.7275934039181`"}], "}"}], ",", RowBox[{"{", RowBox[{"35", ",", "43.55553465502072`", ",", "809.4427874386072`"}], "}"}], ",", RowBox[{"{", RowBox[{"36", ",", "47.9642426782308`", ",", "791.4278245530151`"}], "}"}], ",", RowBox[{"{", RowBox[{"37", ",", "52.568458914063925`", ",", "771.8978259446254`"}], "}"}], ",", RowBox[{"{", RowBox[{"38", ",", "57.33589660332869`", ",", "751.0605288791025`"}], "}"}], ",", RowBox[{"{", RowBox[{"39", ",", "62.229716683430055`", ",", "729.1218670649092`"}], "}"}], ",", RowBox[{"{", RowBox[{"40", ",", "67.20947081521958`", ",", "706.2886453285402`"}], "}"}], ",", RowBox[{"{", RowBox[{"41", ",", "72.23226944984954`", ",", "682.768880165848`"}], "}"}], ",", RowBox[{"{", RowBox[{"42", ",", "77.25408310746869`", ",", "658.7703372429672`"}], "}"}], ",", RowBox[{"{", RowBox[{"43", ",", "82.23109212245704`", ",", "634.4977914847046`"}], "}"}], ",", RowBox[{"{", RowBox[{"44", ",", "87.1210023452244`", ",", "610.1496185270157`"}], "}"}], ",", RowBox[{"{", RowBox[{"45", ",", "91.88425665727087`", ",", "585.9141064967284`"}], "}"}], ",", RowBox[{"{", RowBox[{"46", ",", "96.48507784416087`", ",", "561.9659218252674`"}], "}"}], ",", RowBox[{"{", RowBox[{"47", ",", "100.8922657185529`", ",", "538.4631395693636`"}], "}"}], ",", RowBox[{"{", RowBox[{"48", ",", "105.0797465995039`", ",", "515.5449673957179`"}], "}"}], ",", RowBox[{"{", RowBox[{"49", ",", "109.0268710249764`", ",", "493.3302707017075`"}], "}"}], ",", RowBox[{"{", RowBox[{"50", ",", "112.71847439668524`", ",", "471.91691282870863`"}], "}"}], ",", RowBox[{"{", RowBox[{"51", ",", "116.14472703082326`", ",", "451.38185986732424`"}], "}"}], ",", RowBox[{"{", RowBox[{"52", ",", "119.30074151526452`", ",", "431.7822399507738`"}], "}"}], ",", RowBox[{"{", RowBox[{"53", ",", "122.18614254927216`", ",", "413.15656734522696`"}], "}"}], ",", RowBox[{"{", RowBox[{"54", ",", "124.80449444682735`", ",", "395.5265745061537`"}], "}"}], ",", RowBox[{"{", RowBox[{"55", ",", "127.16269505332941`", ",", "378.8992366314366`"}], "}"}], ",", RowBox[{"{", RowBox[{"56", ",", "129.27035962138007`", ",", "363.26890362731496`"}], "}"}], ",", RowBox[{"{", RowBox[{"57", ",", "131.1392238371565`", ",", "348.619424405366`"}], "}"}], ",", RowBox[{"{", RowBox[{"58", ",", "132.7825870054788`", ",", "334.9261810643746`"}], "}"}], ",", RowBox[{"{", RowBox[{"59", ",", "134.21480995703445`", ",", "322.15796778445`"}], "}"}], ",", RowBox[{"{", RowBox[{"60", ",", "135.45087589099563`", ",", "310.2786771846732`"}], "}"}], ",", RowBox[{"{", RowBox[{"61", ",", "136.50601697586043`", ",", "299.24877708845776`"}], "}"}], ",", RowBox[{"{", RowBox[{"62", ",", "137.3954059077029`", ",", "289.02656346160035`"}], "}"}], ",", RowBox[{"{", RowBox[{"63", ",", "138.1339088354608`", ",", "279.5692123172848`"}], "}"}], ",", RowBox[{"{", RowBox[{"64", ",", "138.73589458724823`", ",", "270.8336298996664`"}], "}"}], ",", RowBox[{"{", RowBox[{"65", ",", "139.21509303510436`", ",", "262.7771272173285`"}], "}"}], ",", RowBox[{"{", RowBox[{"66", ",", "139.58449576883376`", ",", "255.3579400192219`"}], "}"}], ",", RowBox[{"{", RowBox[{"67", ",", "139.85629383562934`", ",", "248.53561291319994`"}], "}"}], ",", RowBox[{"{", RowBox[{"68", ",", "140.04183706426318`", ",", "242.27128678415897`"}], "}"}], ",", RowBox[{"{", RowBox[{"69", ",", "140.15162483618488`", ",", "236.52788154562646`"}], "}"}], ",", RowBox[{"{", RowBox[{"70", ",", "140.19532133461271`", ",", "231.27018484430334`"}], "}"}], ",", RowBox[{"{", RowBox[{"71", ",", "140.18177560514135`", ",", "226.464900351519`"}], "}"}], ",", RowBox[{"{", RowBox[{"72", ",", "140.11904464912394`", ",", "222.08067114532605`"}], "}"}], ",", RowBox[{"{", RowBox[{"73", ",", "140.01442846152173`", ",", "218.08805825327255`"}], "}"}], ",", RowBox[{"{", RowBox[{"74", ",", "139.8745243227997`", ",", "214.4594545124398`"}], "}"}], ",", RowBox[{"{", RowBox[{"75", ",", "139.7052719059649`", ",", "211.16900863533243`"}], "}"}], ",", RowBox[{"{", RowBox[{"76", ",", "139.51199731715246`", ",", "208.1925439655472`"}], "}"}], ",", RowBox[{"{", RowBox[{"77", ",", "139.29945284917915`", ",", "205.50748326690163`"}], "}"}], ",", RowBox[{"{", RowBox[{"78", ",", "139.0718619917038`", ",", "203.09275356675752`"}], "}"}], ",", RowBox[{"{", RowBox[{"79", ",", "138.83297701027456`", ",", "200.9286481477497`"}], "}"}], ",", RowBox[{"{", RowBox[{"80", ",", "138.5861149633196`", ",", "198.9967447751001`"}], "}"}], ",", RowBox[{"{", RowBox[{"81", ",", "138.33418387089606`", ",", "197.27985305369242`"}], "}"}], ",", RowBox[{"{", RowBox[{"82", ",", "138.07972950242615`", ",", "195.7618970541213`"}], "}"}], ",", RowBox[{"{", RowBox[{"83", ",", "137.824974398206`", ",", "194.4278161267936`"}], "}"}], ",", RowBox[{"{", RowBox[{"84", ",", "137.57184398256484`", ",", "193.2635035175955`"}], "}"}], ",", RowBox[{"{", RowBox[{"85", ",", "137.32198818805807`", ",", "192.25576044958407`"}], "}"}], ",", RowBox[{"{", RowBox[{"86", ",", "137.07681750912172`", ",", "191.3922004475671`"}], "}"}], ",", RowBox[{"{", RowBox[{"87", ",", "136.83752792918025`", ",", "190.6611859810161`"}], "}"}], ",", RowBox[{"{", RowBox[{"88", ",", "136.60512405049028`", ",", "190.0517680973482`"}], "}"}], ",", RowBox[{"{", RowBox[{"89", ",", "136.38040781877862`", ",", "189.55376308727222`"}], "}"}], ",", RowBox[{"{", RowBox[{"90", ",", "136.16406897241248`", ",", "189.15742349976534`"}], "}"}], ",", RowBox[{"{", RowBox[{"91", ",", "135.9566108002886`", ",", "188.85376639808422`"}], "}"}], ",", RowBox[{"{", RowBox[{"92", ",", "135.75845122968914`", ",", "188.63417629579212`"}], "}"}], ",", RowBox[{"{", RowBox[{"93", ",", "135.5698619620409`", ",", "188.49070643391985`"}], "}"}], ",", RowBox[{"{", RowBox[{"94", ",", "135.39104547774446`", ",", "188.41574491258493`"}], "}"}], ",", RowBox[{"{", RowBox[{"95", ",", "135.2220910110908`", ",", "188.40227148749275`"}], "}"}], ",", RowBox[{"{", RowBox[{"96", ",", "135.06303212265527`", ",", "188.44357576410223`"}], "}"}], ",", RowBox[{"{", RowBox[{"97", ",", "134.91381598863748`", ",", "188.53347627278123`"}], "}"}], ",", RowBox[{"{", RowBox[{"98", ",", "134.77434425252235`", ",", "188.66609321872292`"}], "}"}], ",", RowBox[{"{", RowBox[{"99", ",", "134.64445923212637`", ",", "188.83598101828707`"}], "}"}], ",", RowBox[{"{", RowBox[{"100", ",", "134.52396071642178`", ",", "189.03804225302974`"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.44836836003375*^9}] }, Open ]], Cell["\<\ This will export the previous table into a CSV (comma-separated value) file \ that can be read by Excel\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Export", "[", RowBox[{"\"\\"", ",", "%"}], "]"}]], "Input"], Cell[BoxData["\<\"PredPrey.csv\"\>"], "Output"] }, Open ]] }, WindowSize->{628, 580}, WindowMargins->{{22, Automatic}, {Automatic, 6}}, FrontEndVersion->"6.0 for Microsoft Windows (32-bit) (June 19, 2007)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[568, 21, 87, 2, 29, "Text"], Cell[CellGroupData[{ Cell[680, 27, 1076, 35, 112, "Input"], Cell[1759, 64, 1147, 36, 86, "Output"] }, Open ]], Cell[2921, 103, 64, 0, 29, "Text"], Cell[2988, 105, 462, 12, 72, "Input"], Cell[3453, 119, 76, 0, 29, "Text"], Cell[CellGroupData[{ Cell[3554, 123, 596, 18, 72, "Input"], Cell[4153, 143, 649, 23, 50, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[4839, 171, 581, 17, 52, "Input"], Cell[5423, 190, 19985, 335, 242, "Output"] }, Open ]], Cell[25423, 528, 123, 3, 29, "Text"], Cell[CellGroupData[{ Cell[25571, 535, 98, 2, 31, "Input"], Cell[25672, 539, 44, 0, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[25753, 544, 339, 10, 31, "Input"], Cell[26095, 556, 11145, 305, 658, "Output"] }, Open ]], Cell[37255, 864, 127, 3, 29, "Text"], Cell[CellGroupData[{ Cell[37407, 871, 101, 2, 31, "Input"], Cell[37511, 875, 47, 0, 30, "Output"] }, Open ]] } ] *) (* End of internal cache information *)