(* 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[ 20954, 454] NotebookOptionsPosition[ 20474, 433] NotebookOutlinePosition[ 20820, 448] CellTagsIndexPosition[ 20777, 445] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ RowBox[{ RowBox[{"implicitEulerLogistic", "[", "n_", "]"}], ":=", RowBox[{"(", "\[IndentingNewLine]", RowBox[{ RowBox[{"h", "=", RowBox[{"1", "/", "n"}]}], ";", "\[IndentingNewLine]", RowBox[{"y", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"n", "+", "1"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "1", "]"}], "]"}], "=", "0.01"}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "<", RowBox[{"n", "+", "1"}]}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"initialGuess", "=", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}], "+", RowBox[{"10", "h", " ", RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}], RowBox[{"(", RowBox[{"1", "-", RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}]}], ")"}]}]}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "[", RowBox[{"[", RowBox[{"i", "+", "1"}], "]"}], "]"}], "=", RowBox[{"yNew", "/.", RowBox[{"FindRoot", "[", RowBox[{ RowBox[{"yNew", "==", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}], "+", RowBox[{"10", "h", " ", "yNew", RowBox[{"(", RowBox[{"1", "-", "yNew"}], ")"}]}]}]}], ",", RowBox[{"{", RowBox[{"yNew", ",", "initialGuess"}], "}"}]}], "]"}]}]}]}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", "y"}], "\[IndentingNewLine]", ")"}]}]], "Input", CellChangeTimes->{{3.6647340662923274`*^9, 3.664734125098843*^9}, { 3.6647364349355927`*^9, 3.6647365369836917`*^9}, {3.6647366851973925`*^9, 3.6647366883065925`*^9}, {3.664736948504409*^9, 3.6647369636200886`*^9}, { 3.665042948715598*^9, 3.6650430133673797`*^9}, {3.6650434236718383`*^9, 3.665043615040267*^9}, 3.6650436827271633`*^9, 3.6650437171525264`*^9, { 3.666772834086974*^9, 3.666772925019039*^9}, {3.6667731078824453`*^9, 3.6667731392497787`*^9}, 3.666773309872548*^9, 3.6667733851237755`*^9, { 3.666773417124422*^9, 3.6667734496214476`*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"n", "=", "50"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "=", RowBox[{"implicitEulerLogistic", "[", "n", "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"plotAppr", "=", RowBox[{"ListPlot", "[", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"i", " ", "h"}], ",", " ", RowBox[{"y", "[", RowBox[{"[", RowBox[{"i", "+", "1"}], "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "0", ",", "n"}], "}"}]}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"plotExact", "=", RowBox[{"Plot", "[", RowBox[{ FractionBox["0.01", RowBox[{"0.01", "+", RowBox[{ RowBox[{"(", RowBox[{"1", "-", "0.01"}], ")"}], RowBox[{"E", "^", RowBox[{"(", RowBox[{ RowBox[{"-", "10"}], "t"}], ")"}]}]}]}]], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "1"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", "Red"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{"plotAppr", ",", "plotExact"}], "]"}]}], "Input", CellChangeTimes->{{3.666773329803655*^9, 3.666773333085955*^9}, 3.6667734693353815`*^9, {3.666773519837705*^9, 3.666773600336671*^9}}], Cell[BoxData[ GraphicsBox[{{{}, {RGBColor[0.24720000000000014`, 0.24, 0.6], PointBox[CompressedData[" 1:eJw1UWtIU1EAHlI+JoShf1KbKSla9MtgpK5Pk1ChwC2zsoWVKBFhmik+MAyc FSaL0tIVG+jAEh+IJmpijwmznKbL4eZ7zjm3u+06YSuUJO/d6YPD4by+14m8 fV+U78PhcOL3BzP/R1VId5qxshbeuREBGnciJM/IWglZSrSSo5ZikCehhIPd iLmaqHCrXpHzfrhv9KYmS5sha2EwBFfQOYtvhZzcH4Unqd8VsdIKu5DZ+YoY /fCSRdRO3o/h7jHf9otfOsGdaysScNV4erDi06mOXsL3HdRYU6hW8RFxYuk3 d6wG27WGJydyhgn/FP7kbaS8dX2Gcv/1XNs0kkuzNEaBiujNgJu5N1lyQI14 Fr8wVv7II8mcIPqz0Nzrq8kImobKHbuvoEMo/5BPwoCW+JlD1qRWVvhSBxEr qIdAax/RteiJPwNe3Hmtrx9dgLGSIZzHmcCCDMPWMvG7iM7aB+H8ICOKGXtF S8iPTj0+K18j/pfRtVgdf+XoOry/sYo6+Uo5r85M8qwi8nxCk79uA+leAaQd CTNU+W2SfEbsqgWZ6eFWMGpi6Ro6Gk03i4NtJK8J1YazQ5TDBm8cE8S5H7au dVEk/zqK837/nbpkRwRjh2fGQlzYQ5nJTvoww7X7wxqR60ABiw0IrzefHhp3 kH4sECreRLVGOcGwhXRbUFaQt00VOklfm0hU1I+W9TjhYejcm+D5NQhos5P0 ZwW4O1U/g2mw9QhsmL7Ar7mVQJM+bSjhZweqcmiw9iQUlAbRurWUJv3a4ZkI yHj/nMakhoEdj1cvp/YqaNK3A3qOZPhdDw02bogTA20NJ6kRmvTvxI58Zc9/ nAZLJ6YhSpqpFM7Q5D+2MO87mH3YQOMfrQCbnA== "]]}, {}}, {{}, {}, {RGBColor[1, 0, 0], LineBox[CompressedData[" 1:eJwV1Hk8VPsbB3DbjVtRQjWWK+m2iEkq7X3qqihKli4VspYkrZarhBqSLUlI 1iyFjLIkS2SdKLImZEkZM2NmvslFxeV3fn+c13m9X+ec53yf7/M8XzX786ZO YiIiIi3U9f/7QSdOy2vu0V2ZLFrYZ28GJPQUL/TTNqLGYt9oazADG5Ji4/Np e/GjyTqxOY4Bw8bOvASaORzlL/z4VsKAhePkq0iaI0Q+GT6vm2Kgd90tvTDa FfjJJ7g3+QYgWcE4MJjGwKcY0dzIu4Hwl8pfdZsWhcMOlUeX9QVBQ/9uUN6F KMhe7p+VHQtCa8B5bg8rCkZ9TQoKUrexQlwzm+5xH5Y7sl4Z69xG3X9p9PbW aJRvcQuJvXUb88eiN6mGPkBKNMlnbApGbK+3XuFMIhreSKyiJ4XA8lhK7rfR NChoN7xYnhMO0mM+95FiOpzrCsW1a8MRaCPlZKaXjiK9C4PGveHIdziv+CIq HbX3voaypO9Axm1nwNXNGUiNG1GOPX8HtTe6js/xeYzn3nv2OClFQOep7Bxl qSy03fn1gM+NwLxZfyv9P3KxOcX497H0SMSesbDR3J+LDaIpXjsKIrGiXdNW 1i0XXtMF6RFVkdiV2enQ8yoX+na6lsf6I1GcGV1YRH8G1W9bf7rS7kFgxi8z bHoG9SZ52bQ793D0Sczby9J5eL3Tr+s2Iwp/mgg51SEF2F0peUkyMhqKvvrb TqcU4GVEh2pRWjQW5iSHzC0qwPZl4LsXReOXlBndZLAA7kNBC9R6o9H8uuhy 35ZCEHftN/Q1MfDW9p/5MVQIDQNlG52aGDQtkJen7y5Ct0l1Nn3OA3g2bUXs v8WIO9GaV1r2EDL0wqF6qRJM+aqslW59iLQw7dAp5RJ8v3g2wm34IVqMVnXZ 7CsBg9592G5RPNa+lb/8Z3QJ9ihIfZA8E49+FsnI21yK3xI1pu1pCdCvTJdp 9C6DVMJFxqPARCwpkO0TFa1AjtauNd8bkmEVOmX+dXEFdl2/v6n7YzIeOQ69 rdOsgInfFbtGdjI0FYqLQywrMBhdv6VDLAV/edjeV3hWgXU8buPZ7SnQjbfI iFr3GmZmDl9KmSmws8mwftL3GoExxUfUEx/hxaDeu/c7qmAbxMmdTknDBwP/ 0TD9KvwR8OxmSl4aJpjli41Mq5BjY7rBrDoNm69utXtzugodYtsP9n5Nw0u5 dROvI6uwN9/UPGx1Oor3Ki7L41ThS798p2NeOsoyvl2OiqpGS/wnP53mDFS5 JCgd49cA6RdCT+pmIp979qzNRA0qOC2GGgcykea8rdRBpBa7VsqbilllIuBU 53E3+VqcCDtj3u2fCQN72bibO2pxzt+WtrYpE++PBSxlhtbiabnlj7/PZOGT wTkFca06OPUstpPLysb4yp0Lcs6xsLnnfoeEBRO7Bv9iBXuxELt0c5HBeSYC Ewx8nW+y4D2woC7xFhOL5c2J+gMWjvwnGuf/kgld0bPv42pYeJcoba+uRPVx T+ydIKU3OPmksfjzl1xM3/l3gWP9GzTFXTNqqnwGsZ/ZC5X+bMD8FPcNlr55 0J21fWHGfIf5wmSPCZ0XeMfY2FSn8R6el728DD1LoBox8jNT2IzjqdP7HrMr MDLZ39eh3Yr1lbrXY4erkMYrbZAJb0NSsv/6w3F1kFF5TL5+bUdOad5Ln+oG hJ4pPJW84gM0taxdB5a9x8Tu7K/Ovp0I8j1yrpiKc0WnpmBx9UewlgcZ6va1 o59ZE9Ki1g3XKNazyKlOTLm0P1zj3YMjtJmVSzf1YO0eyZlTpZ9w9JiKvpFx H/QfqZ1+vqQPlVp/H43rHoCnysnln636UW+RoH3VdRCa/KPzOg8N4Fe65akZ 3hfILfgmVfNrAGROqvpL9yFIVDzMVY//jIPMo943Gtk4u7G063fDQXB9nBz5 yznIPmCfXD4xiHcs7j6pI1yMjyguZEd/wdClvCt1B3goNNMXj9r7FaNDkoes N4/A1yAgqYD9FQxtydTgyRGo9R78mBowBEtZGuNAHB/bNjOiaDpsvFjbo8/W FsDqj/c/RdvYEIavVil/LIAM7VJhaQcb6i5iO2IyBaiUU7jo/pGNhuTd4uez BVj5+wkOp5eNJUFLwhVzBSBj7A/vOWxUPLbPtnshwI36mfyEGTbirJwLn9UI 8OQK3W2rxjDuWjkxLAYEGG8IG7zgP4y5p3g7x+SE+DSvw3n65jA0uv+9a68g RJWRMrkVOAzpBp2h5sVCRDRlTSWGDOOzXG5IFk0IjVaWfOP9YRxyP801URXC pktUf00WFa9p6bmLGkLUDbtnf24dRvAWy4JJCBErYXPFdAUHe+k2fnuchZhx WLKpZSUHH8ju9efPCOFY3TxuvIYDmc/BMvEuQmjf0PM8ROdgz6ztkTFXId7M rvY22MJBSob9bMxFISZ/jfnuMuIgfu/payXeQvw9ejtE4woHHm93XjkdKoRc X+EjsRoODq2zZV7LESJ8vo9PUR0Hsd9WvgpnCjF3+15L13oOnD3CJ5NyhRCL bZX+0MRB2+Go0YrnQnw3IV6ZXRxM6BTPnSgUorl2tfERwsHCDBPT7eVChOY8 nEpQ5OKMNVOL3iTEbz7+5lsvcpH/rZnVNiJEwSj/Z8VlLhjeMmoZfCq/U5ZJ +z24+NKb2u8pEKLGeB3P7CoX36a9w2lECMbyXl+3AC72/W1cfvS7EBJvtmSn PuAiTcuj8vkPan1y30Rlqrj4Pvlz2ZQEgUimTe7gIh72tqkMcpQI6vidbkkK PCgeH/0rQ5kgVNuEbrWUh/QS0XkOKgRLXurldKjwcHXycWv3HwR01ursN6t5 WOTyyKhajcBq6HsGcxcPAw40Y69VBMVqtxK9XXjgWa5IcNchuBSXG7aoiger J4eH2/cTFE3H3aqq4eGoivMTV32CaZtA/0ssHuSuzRBxA4Jb6tbure94aFsv L7/+AEH807nWkZ08LLYeLw4wJKgtd9JcJOChofbIu6VHCJZ+UXoru3QE6rXt 1ziWBK/WBknJuo1AfHNh6+4zBLrtSxojL1BzJG1yr4Dys2uP78pfHoHExvzp VS4EjxrrFJd6jeDkPFaD9FmCoPO/aaneGIF2/7wTHa4Epvk3TLSiRyDd+s8N 8wsEw9t84g6UU3OqMx0j70kge/CS5g0ZPpQ91uY/uUmQXLlg3jZZPuw0F4Ut ZBBob83hjsrx4SQaIetF+fBqToYdjY/7o7ZK+wMIgufYqO1ZwceYYdnWvkAC ieoDi0W28aE72dkzc5tgYruayHUnPhR23/QVjyD4pNXc9k8ZH67vRCZXPSRI GpH1Ea/gY+H15x0elO0zzVaGVfLRzfTwrKXMXdHplVzHR2L/uIR9PBVPqU+F 1czHiQfVRrEJ1Prn8k/LD/ERs0mu6GcSwX625DRTWgC96Odr76ZR+5WIFV9t BMjKm9OllEMw5j0maWcnwFCzT8gJyroWj0d6HQQIeOov/ZBy6YIF+Z3OAsj2 VfjRmFR/+fXteXtJgDIrsdXyuQQ99j4n8wIFcP6cEyvynGDOqtI4P6YA6OCP vyggsM7VlVWZEeAQaeEllxK8OP5DplmEmmP3p+LdlGXmlMy/KS6EyPVdgkVl BK+tdkpxpYTw5i/fwKC8fO7e2QLqXDKIUQp3fEXAdjQRGK0RQrflU4pyBYGb omv9NXMhcrxV9c5WEfgFpPh/yhLieF6+evAbAmlPyYN21Nw3cze4FVKOcz63 iE3NuabYipgByvmGW1NJgRA3XUoYuvVUfNmWGrEKIa6KB7cMUDZMFJFa00ad i/W1+uveEiwuOhnuPiWE4no3s4JGghyOykMZI4LUzOxFjq1UvyUE2PYdJkgJ n9cTQtnRVPAn04Qg4hrXL5+yclnZs0MWBDVl953E2qh+uXOiLtSeYNyZfT+R srPug9G5/xB0jL3Lbm6n8mfIG0hmUHW9yb6k3Ekws+WadOcT6jvr28t3U+4S fGnNyKbqWE7Pd6AcYZFvvZ+qw+G7I55ZlGc1TC8HUPsYw5ON1P1I9VvLnURx 6j+qfnzVA10ELwMnHdo6CJ5vkeC7UI7afnJNKvXe+MnyuFDKhun0gr/6CNSn rye/p1zi1VjvzyU4/Z4Ym3cTRGttumPCJ/jjR9cGd8oXB+PN1QjBl7P0qfuU Vxu59r/+l8r3ntfGTsoSou1pEZMEDL2kjEnKA4XbXWx/Uc83RU4v6SEoc0ld p/0fde4UnNHdQjlGdd747CzB/Imd5paU/wcrpKuc "]]}}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesLabel->{None, None}, AxesOrigin->{0, 0}, Method->{}, PlotRange->{{0, 1.}, {0, 0.9952475255706432}}, PlotRangeClipping->True, PlotRangePadding->{{0.02, 0.02}, {0.019904950511412865`, 0.019904950511412865`}}]], "Output", CellChangeTimes->{3.6667734548891487`*^9, 3.6667734949641004`*^9, 3.666773559860033*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"n", "=", "250"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "=", RowBox[{"implicitEulerLogistic", "[", "n", "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"plotAppr", "=", RowBox[{"ListPlot", "[", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"i", " ", "h"}], ",", " ", RowBox[{"y", "[", RowBox[{"[", RowBox[{"i", "+", "1"}], "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "0", ",", "n"}], "}"}]}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"plotExact", "=", RowBox[{"Plot", "[", RowBox[{ FractionBox["0.01", RowBox[{"0.01", "+", RowBox[{ RowBox[{"(", RowBox[{"1", "-", "0.01"}], ")"}], RowBox[{"E", "^", RowBox[{"(", RowBox[{ RowBox[{"-", "10"}], "t"}], ")"}]}]}]}]], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "1"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", "Red"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{"plotAppr", ",", "plotExact"}], "]"}]}], "Input", CellChangeTimes->{3.666773576403733*^9}], Cell[BoxData[ GraphicsBox[{{{}, {RGBColor[0.24720000000000014`, 0.24, 0.6], PointBox[CompressedData[" 1:eJw1lgk0VXsbxrtSuuijPk1SoVJdVCJyK480ap4uXUWDqDQXpYE0SJOpkDRQ JBpUGkiROUTm4ZxjOJzBGfem3FuS+s7e/t9ey7LOwrvf4Xl+D6Nt+9d6qPXr 169H9cV8//9zQi9lcevxc/iR3Fm1YtsByOcc1Pq6NpB8DkBHVvFbruV5fPf3 eud97QIOL+f+Eo0JIj8PwxLNEr7PkAvoqxOBb+l+Ytl/LpLfj0bAAqVV8vBL OCY1sRj0+BbqCnUfjpxymfx9HL6obf6RsOgKZresCr6+7x6G29dxMvcHk3oJ 8JsZpFGfEIIMR20106BE+G6PPr1THErqJ6FsB2/F0Zgw2L4sOpIV+xDL0r/o +iwJJ+97jKNhtXMO9IYjfWygfE16CvbsTW198uYqef8zRBl/8/Pxu4Zb+YbP Du18DnXLP3SGOEaQflLxRjvLuMwgErxu5oUvsD9HKf72LZL09xI9oYdni5ui MHpa04SBp1/hvOW4XOnH66Tf1wg6ZkW75EXDZTtTMQ16eqdHReTfIP2nw6tn 2KXP5TGIucE8b8D9N+RynegmmScDz4ZW+sYPvA3Op6/rBPffIr3J5OP7GXfI fO8wfVHspEivWIxSd/qPeWomUqRFM5c4x5F5sxDvrtcY9isOG9iC7zHt9Q47 t5S7ZP5sdP45PVvufg/LbS7O2yLKxueuI95HDOPJPnIwseCH4VxhPBRrmIo5 6D/kmH3A0wSyn1y0vBf0P3PmPi7vae9J7sxFhFHq2vOuiWRfedDtVRRX2T+A adBw1Qh54DxXz0gxSyL7y4ezxuLx7kbJ+Hh34esvvfnQ6Wl49ue4h2SfBfCJ HZZ62uQRmG2NVC/EU7nF0PM2j8l+C2Fjp77y+ton0KyPPzBX8wMSV5TtGOib Qvb9AYV/tzRMTHqKh8y4ukVIrDmVnXDwGdl/EY60pTgK9J5jKVNueDG03vpl d2c+J/coRtfU6R3/3ZsKdt0GJZio/ucUN+MX5D4lUPDMZgxofgFmexXGH+Es 1bOxjntJ7vURsfSMtLG7XmHKptDcfyaXIvS12kHFrNfkfqUo1PVfXqaThmJV Nf1pZdiQETyfVqaRe5ZhQW7MkoDqdOwMV57EzE+IEliZpmS/Iff9BMPp+ccT X2VA1Zxq5HJ0Z9jtiEp9S+5djvFhTbEP095hveD+zXzDCpjPKZ+mU5BJ7l+B Bm/pBpqbhQTV9urjK1BssHmPT/d7oocKzO3IWe0Rlo2u3hpfqUklfp37y+fy 5Byij0o0dozabl+QgwWsACpxKkB5Jcozl+ilErVPrI/GaudBdTzVQFVYmP4r aXdaHtFPFRr8dwZr7sgH296zKhg0BbpcGl1A9FQFry9F/36tKYAl+1Sj8ZtW s9u1QqKvahw64DGw2OkDzmaXqBRUjZBKUb7duCKit2qk8DxTCxVFqFZdf4Nt DaScqDCP98VEfzXg/Nozbez1EoxnC9Yg4I31us7DH4kea7Db6sDNpr9KcZg9 QC1Gn73l1TKnjOizFt991pV3TP6EvH8mqy5cC+2wd/t+1y8neq3D9vvOs3QG VWAo22AdFoQ2WIcHVhD91iEw7dzTMwMqwXT3sqgOdvPCl3CDKome6yE8bb08 WLsKqbfpB4XL6nHr5smE0GtVRN/16P5+vKHaoBpq7MD1COUZai1LqiZ6bwA3 NaZHOLMGa1nBNGBkrodGWEEN0X8DIkcZTrV1rsU9dtwGLNU1HFEqqyV+4KAw jjvSMqAOn5lyGzjY62jg4DyinviDA466nofhs3o4MPLncfDlrekEn6UNxC9c fJ5/v99scQNY+bpxces25eV0lkP8w0VExc2ah0ZctB5nDsJFiU9Hnn4Ol/iJ B++WZknQVh4sWIHzYNXoZlmt1kj8xcMQz++yxoRGnI6cr5KM6rPV7g/nFjcR vzXikpmR5kNZEyqNmQM3IsY0KNo0pJn4rxFmpzjDJRYtMFKpOVzZiJvSzf2z aluIH5vwtVspKZjDx0HGHgea4E/59Y8O4RN/NiHVoTJ5AZ+PHEZ+X5qg81jm /Nailfi1Ge/NwjPkZ1rBnvdIM47rT3lXUN1K/NuMVSZjw20ntIGdtrsZUa9K 8my824ifW2DjVoLneW1g2zvZAhONB5uuDxUQf7dAQ62qnbdFgL40V/WxvdH8 YIqA+J0P706B2tIeAZIKl6kczUdRh2nJ34uFxP980H68kxFXhVitcuPQs3ws b2p8SjUKCQ/4cAjiRbuaiPBtHVOQj3XhW4bW7xMRPvDhrKkTu+a1CHH6jCFa 8ehJ/z35P0SEF60w7+e6dpKDGEv6Doxg78/1h86LCT9aYRTtMC2uWIyOB0yD rbh69q51slY74UkrBnCfH7mwvB3R+xgHtcI20NLQ9ko74Usrtgyz+j2tpB0q caks1YbJtk+Sf2pICG/akGhyZ8TgBRKwNFnRhuuaT9T4/hLCnza0PnR96pMu AaOWTaFtmOQm6yztkBAeteHEqtLJoklS9OGwDVTmL+ssVynhUxsMg3T9N1yV gi2nusvm2u/1LwukhFcCxNrLw2q+SsHiaJ0AnSOlW95MlhF+CWCk8Pbw2CAD i5NIAXjJE0bzzssIzwTQGLbIecxLGfpwIICLb+i8aXwZ4ZsAg7iiFXpacpxi gSOElgWsK6zkhHdCXNmR2OXpKsck1kBCOM73zGw8Jyf8E6IpdnP3zEdylDN2 uinE3ewQ+6MVcsJDIc625q6/3yXHUdYOQpyJvmRQMEJB+ChEu7bzHZ6tAuMY OY8VYfP5+h6li4LwUgTZppC9v51Q4MMydoP4SxBrNzpGQfgpwtUAXpF9ugL7 WUCqdFTx88DRWgXhqQiX1s2Iz+9UYAS7PhGMJFnl5oOVhK9iCHV6dqRPUoJt b4IYspx9vF3zlIS3YhRsM4tc7KKEJ/uIMXiVa7fTISXhrxi8Bf3cb19UQmVm FWHFCGmqNzKMUxIei2FiFeQrfKUEQ0dfqRgt+dwYukRJ+NwO9QuzHy9qUcKN LdiOxPlLs6jPSsLrdvgPvluuGEBhYB+gEPg1vWThSIrwux1q4s5zP6dQYK6h l9KODxEZmaNnU4Tn7XjErRtzbxkFJ7bBdoQvzLWO2kgRvkswx9r8zA8vCj8Z u1lIcOX51a2cYxThvQTmZqIxFhcoJDJyPiTB4MTV7RpRFOG/BJITurVb4yms ZAeWIMHTWeTwjCJ5IIHUYklvyjsK/zLn+EcCh4JeSXIRRfJBimqzuxttaijc Yca1Uel8yxXapYUieSFF1w1fk5EyCov6Fogz+mm7DnVRJD+kGLf+sObenxQo Jh4ypDA52jtPaxBN8kQKLV+DGauH0IhittcjxdMUT5GdPk3yRYbf7eNWcY1p sHieK4N6+NiVE01pkjcyjNCZnmRsSYPFq78M2klBOrV/0iR/ZBA4PTKEAw0W f9kyHLLxcHd3pEkeyRB/23rZgtV0H17U5HidVmojdKJJPsnRu/yIz0JXGqx9 58shsoyZ6+VOk7ySQ3P81LT1u2iw9giUY8dF868D99Mkv+QINNL97aw3DVZ+ hSrffseHkmM0yTM57N6GmPL9adQy5x2kgGdr8q/8szTJNwU4tcEDTlyg4cd2 qMCk2K90/2Ca5J0CLUNqt20NpzGRFaACU11efImIpEn+KbAmx/SPuBs0ykqZ R4GBkSsDzt2mSR4q0KOv7T7/Lg0fFqhK3PG7OaM5gSb5qETTvRYHpyQaY1hB K+HIL+Y9fkSTvFQi50yGrSCFRoEhWxBGxi3B35/TJD+VyHwb8vPLSxp7WUAr sSEs/kFFGk3yVImp+VUx4Rk0WFzpUZiV8xvfJpMm+Upho+sNn/z3NFgcOFE4 /qFu9axcmuQthVvQd4vIp8Ha7TqFv2MtT3MKaZK/FPLSH0cOKqbR9+8Ohdyr Q/Ybf6RJHlO4ccmaM6mMBisXle66tp2MHlVOk3ymsdFxqNe3ChrsOTbR2M2b qV5QRZO8phE58+DwgBoa6n0GxuFh9Y5/1NEkv2lo6c+dl19Poy/uaVw3876x kkOTPO+An8EpaTGXxv8AasgHMw== "]]}, {}}, {{}, {}, {RGBColor[1, 0, 0], LineBox[CompressedData[" 1:eJwV1Hk8VPsbB3DbjVtRQjWWK+m2iEkq7X3qqihKli4VspYkrZarhBqSLUlI 1iyFjLIkS2SdKLImZEkZM2NmvslFxeV3fn+c13m9X+ec53yf7/M8XzX786ZO YiIiIi3U9f/7QSdOy2vu0V2ZLFrYZ28GJPQUL/TTNqLGYt9oazADG5Ji4/Np e/GjyTqxOY4Bw8bOvASaORzlL/z4VsKAhePkq0iaI0Q+GT6vm2Kgd90tvTDa FfjJJ7g3+QYgWcE4MJjGwKcY0dzIu4Hwl8pfdZsWhcMOlUeX9QVBQ/9uUN6F KMhe7p+VHQtCa8B5bg8rCkZ9TQoKUrexQlwzm+5xH5Y7sl4Z69xG3X9p9PbW aJRvcQuJvXUb88eiN6mGPkBKNMlnbApGbK+3XuFMIhreSKyiJ4XA8lhK7rfR NChoN7xYnhMO0mM+95FiOpzrCsW1a8MRaCPlZKaXjiK9C4PGveHIdziv+CIq HbX3voaypO9Axm1nwNXNGUiNG1GOPX8HtTe6js/xeYzn3nv2OClFQOep7Bxl qSy03fn1gM+NwLxZfyv9P3KxOcX497H0SMSesbDR3J+LDaIpXjsKIrGiXdNW 1i0XXtMF6RFVkdiV2enQ8yoX+na6lsf6I1GcGV1YRH8G1W9bf7rS7kFgxi8z bHoG9SZ52bQ793D0Sczby9J5eL3Tr+s2Iwp/mgg51SEF2F0peUkyMhqKvvrb TqcU4GVEh2pRWjQW5iSHzC0qwPZl4LsXReOXlBndZLAA7kNBC9R6o9H8uuhy 35ZCEHftN/Q1MfDW9p/5MVQIDQNlG52aGDQtkJen7y5Ct0l1Nn3OA3g2bUXs v8WIO9GaV1r2EDL0wqF6qRJM+aqslW59iLQw7dAp5RJ8v3g2wm34IVqMVnXZ 7CsBg9592G5RPNa+lb/8Z3QJ9ihIfZA8E49+FsnI21yK3xI1pu1pCdCvTJdp 9C6DVMJFxqPARCwpkO0TFa1AjtauNd8bkmEVOmX+dXEFdl2/v6n7YzIeOQ69 rdOsgInfFbtGdjI0FYqLQywrMBhdv6VDLAV/edjeV3hWgXU8buPZ7SnQjbfI iFr3GmZmDl9KmSmws8mwftL3GoExxUfUEx/hxaDeu/c7qmAbxMmdTknDBwP/ 0TD9KvwR8OxmSl4aJpjli41Mq5BjY7rBrDoNm69utXtzugodYtsP9n5Nw0u5 dROvI6uwN9/UPGx1Oor3Ki7L41ThS798p2NeOsoyvl2OiqpGS/wnP53mDFS5 JCgd49cA6RdCT+pmIp979qzNRA0qOC2GGgcykea8rdRBpBa7VsqbilllIuBU 53E3+VqcCDtj3u2fCQN72bibO2pxzt+WtrYpE++PBSxlhtbiabnlj7/PZOGT wTkFca06OPUstpPLysb4yp0Lcs6xsLnnfoeEBRO7Bv9iBXuxELt0c5HBeSYC Ewx8nW+y4D2woC7xFhOL5c2J+gMWjvwnGuf/kgld0bPv42pYeJcoba+uRPVx T+ydIKU3OPmksfjzl1xM3/l3gWP9GzTFXTNqqnwGsZ/ZC5X+bMD8FPcNlr55 0J21fWHGfIf5wmSPCZ0XeMfY2FSn8R6el728DD1LoBox8jNT2IzjqdP7HrMr MDLZ39eh3Yr1lbrXY4erkMYrbZAJb0NSsv/6w3F1kFF5TL5+bUdOad5Ln+oG hJ4pPJW84gM0taxdB5a9x8Tu7K/Ovp0I8j1yrpiKc0WnpmBx9UewlgcZ6va1 o59ZE9Ki1g3XKNazyKlOTLm0P1zj3YMjtJmVSzf1YO0eyZlTpZ9w9JiKvpFx H/QfqZ1+vqQPlVp/H43rHoCnysnln636UW+RoH3VdRCa/KPzOg8N4Fe65akZ 3hfILfgmVfNrAGROqvpL9yFIVDzMVY//jIPMo943Gtk4u7G063fDQXB9nBz5 yznIPmCfXD4xiHcs7j6pI1yMjyguZEd/wdClvCt1B3goNNMXj9r7FaNDkoes N4/A1yAgqYD9FQxtydTgyRGo9R78mBowBEtZGuNAHB/bNjOiaDpsvFjbo8/W FsDqj/c/RdvYEIavVil/LIAM7VJhaQcb6i5iO2IyBaiUU7jo/pGNhuTd4uez BVj5+wkOp5eNJUFLwhVzBSBj7A/vOWxUPLbPtnshwI36mfyEGTbirJwLn9UI 8OQK3W2rxjDuWjkxLAYEGG8IG7zgP4y5p3g7x+SE+DSvw3n65jA0uv+9a68g RJWRMrkVOAzpBp2h5sVCRDRlTSWGDOOzXG5IFk0IjVaWfOP9YRxyP801URXC pktUf00WFa9p6bmLGkLUDbtnf24dRvAWy4JJCBErYXPFdAUHe+k2fnuchZhx WLKpZSUHH8ju9efPCOFY3TxuvIYDmc/BMvEuQmjf0PM8ROdgz6ztkTFXId7M rvY22MJBSob9bMxFISZ/jfnuMuIgfu/payXeQvw9ejtE4woHHm93XjkdKoRc X+EjsRoODq2zZV7LESJ8vo9PUR0Hsd9WvgpnCjF3+15L13oOnD3CJ5NyhRCL bZX+0MRB2+Go0YrnQnw3IV6ZXRxM6BTPnSgUorl2tfERwsHCDBPT7eVChOY8 nEpQ5OKMNVOL3iTEbz7+5lsvcpH/rZnVNiJEwSj/Z8VlLhjeMmoZfCq/U5ZJ +z24+NKb2u8pEKLGeB3P7CoX36a9w2lECMbyXl+3AC72/W1cfvS7EBJvtmSn PuAiTcuj8vkPan1y30Rlqrj4Pvlz2ZQEgUimTe7gIh72tqkMcpQI6vidbkkK PCgeH/0rQ5kgVNuEbrWUh/QS0XkOKgRLXurldKjwcHXycWv3HwR01ursN6t5 WOTyyKhajcBq6HsGcxcPAw40Y69VBMVqtxK9XXjgWa5IcNchuBSXG7aoiger J4eH2/cTFE3H3aqq4eGoivMTV32CaZtA/0ssHuSuzRBxA4Jb6tbure94aFsv L7/+AEH807nWkZ08LLYeLw4wJKgtd9JcJOChofbIu6VHCJZ+UXoru3QE6rXt 1ziWBK/WBknJuo1AfHNh6+4zBLrtSxojL1BzJG1yr4Dys2uP78pfHoHExvzp VS4EjxrrFJd6jeDkPFaD9FmCoPO/aaneGIF2/7wTHa4Epvk3TLSiRyDd+s8N 8wsEw9t84g6UU3OqMx0j70kge/CS5g0ZPpQ91uY/uUmQXLlg3jZZPuw0F4Ut ZBBob83hjsrx4SQaIetF+fBqToYdjY/7o7ZK+wMIgufYqO1ZwceYYdnWvkAC ieoDi0W28aE72dkzc5tgYruayHUnPhR23/QVjyD4pNXc9k8ZH67vRCZXPSRI GpH1Ea/gY+H15x0elO0zzVaGVfLRzfTwrKXMXdHplVzHR2L/uIR9PBVPqU+F 1czHiQfVRrEJ1Prn8k/LD/ERs0mu6GcSwX625DRTWgC96Odr76ZR+5WIFV9t BMjKm9OllEMw5j0maWcnwFCzT8gJyroWj0d6HQQIeOov/ZBy6YIF+Z3OAsj2 VfjRmFR/+fXteXtJgDIrsdXyuQQ99j4n8wIFcP6cEyvynGDOqtI4P6YA6OCP vyggsM7VlVWZEeAQaeEllxK8OP5DplmEmmP3p+LdlGXmlMy/KS6EyPVdgkVl BK+tdkpxpYTw5i/fwKC8fO7e2QLqXDKIUQp3fEXAdjQRGK0RQrflU4pyBYGb omv9NXMhcrxV9c5WEfgFpPh/yhLieF6+evAbAmlPyYN21Nw3cze4FVKOcz63 iE3NuabYipgByvmGW1NJgRA3XUoYuvVUfNmWGrEKIa6KB7cMUDZMFJFa00ad i/W1+uveEiwuOhnuPiWE4no3s4JGghyOykMZI4LUzOxFjq1UvyUE2PYdJkgJ n9cTQtnRVPAn04Qg4hrXL5+yclnZs0MWBDVl953E2qh+uXOiLtSeYNyZfT+R srPug9G5/xB0jL3Lbm6n8mfIG0hmUHW9yb6k3Ekws+WadOcT6jvr28t3U+4S fGnNyKbqWE7Pd6AcYZFvvZ+qw+G7I55ZlGc1TC8HUPsYw5ON1P1I9VvLnURx 6j+qfnzVA10ELwMnHdo6CJ5vkeC7UI7afnJNKvXe+MnyuFDKhun0gr/6CNSn rye/p1zi1VjvzyU4/Z4Ym3cTRGttumPCJ/jjR9cGd8oXB+PN1QjBl7P0qfuU Vxu59r/+l8r3ntfGTsoSou1pEZMEDL2kjEnKA4XbXWx/Uc83RU4v6SEoc0ld p/0fde4UnNHdQjlGdd747CzB/Imd5paU/wcrpKuc "]]}}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesLabel->{None, None}, AxesOrigin->{0, 0}, Method->{}, PlotRange->{{0, 1.}, {0, 0.9954576073149749}}, PlotRangeClipping->True, PlotRangePadding->{{0.02, 0.02}, {0.0199091521462995, 0.0199091521462995}}]], "Output", CellChangeTimes->{3.6667735773344193`*^9}] }, Open ]] }, WindowSize->{1292, 599}, WindowMargins->{{-17, Automatic}, {12, Automatic}}, 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[557, 20, 2426, 60, 212, "Input"], Cell[CellGroupData[{ Cell[3008, 84, 1355, 41, 130, "Input"], Cell[4366, 127, 5814, 102, 246, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[10217, 234, 1258, 40, 130, "Input"], Cell[11478, 276, 8980, 154, 246, "Output"] }, Open ]] } ] *) (* End of internal cache information *)