(* 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[ 33698, 822] NotebookOptionsPosition[ 32521, 778] NotebookOutlinePosition[ 32866, 793] CellTagsIndexPosition[ 32823, 790] WindowFrame->Normal ContainsDynamic->False*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"DSolve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{"s1", "'"}], "[", "t", "]"}], "\[Equal]", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", RowBox[{ RowBox[{"s1", "[", "t", "]"}], "/", RowBox[{"(", RowBox[{"100", "+", "t"}], ")"}]}]}], "+", RowBox[{ RowBox[{"s2", "[", "t", "]"}], "/", "100"}]}]}], ",", RowBox[{ RowBox[{ RowBox[{"s2", "'"}], "[", "t", "]"}], "==", RowBox[{ RowBox[{"2", " ", RowBox[{ RowBox[{"s1", "[", "t", "]"}], "/", RowBox[{"(", RowBox[{"100", "+", "t"}], ")"}]}]}], "-", RowBox[{"2", " ", RowBox[{ RowBox[{"s2", "[", "t", "]"}], "/", "100"}]}], "+", "1"}]}], ",", RowBox[{ RowBox[{"s1", "[", "0", "]"}], "==", "0"}], ",", RowBox[{ RowBox[{"s2", "[", "0", "]"}], "==", "0"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"s1", "[", "t", "]"}], ",", RowBox[{"s2", "[", "t", "]"}]}], "}"}], ",", "t"}], "]"}]], "Input", CellChangeTimes->{{3.4149251778684053`*^9, 3.414925349643406*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"s1", "[", "t", "]"}], "\[Rule]", RowBox[{ FractionBox["1", RowBox[{"4", " ", SuperscriptBox[ RowBox[{"(", RowBox[{"100", "+", "t"}], ")"}], "2"]}]], RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "t"}], "/", "50"}]], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1000000"}], "+", RowBox[{"1000000", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"t", "/", "50"}]]}], "-", RowBox[{"5000", " ", "t"}], "-", RowBox[{"15000", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"t", "/", "50"}]], " ", "t"}], "+", RowBox[{"300", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"t", "/", "50"}]], " ", SuperscriptBox["t", "2"]}], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"t", "/", "50"}]], " ", SuperscriptBox["t", "3"]}]}], ")"}]}]}]}], ",", RowBox[{ RowBox[{"s2", "[", "t", "]"}], "\[Rule]", FractionBox[ RowBox[{"25", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "t"}], "/", "50"}]], " ", RowBox[{"(", RowBox[{"15000", "-", RowBox[{"15000", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"t", "/", "50"}]]}], "+", RowBox[{"100", " ", "t"}], "+", RowBox[{"600", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"t", "/", "50"}]], " ", "t"}], "+", RowBox[{"3", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"t", "/", "50"}]], " ", SuperscriptBox["t", "2"]}]}], ")"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{"100", "+", "t"}], ")"}], "2"]]}]}], "}"}], "}"}]], "Output", CellChangeTimes->{3.4149254222674055`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"{", RowBox[{"u1", ",", "u2"}], "}"}], "=", RowBox[{"{", RowBox[{ RowBox[{ FractionBox["1", RowBox[{"4", " ", SuperscriptBox[ RowBox[{"(", RowBox[{"100", "+", "t"}], ")"}], "2"]}]], SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "t"}], "/", "50"}]], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1000000"}], "+", RowBox[{"1000000", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"t", "/", "50"}]]}], "-", RowBox[{"5000", " ", "t"}], "-", RowBox[{"15000", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"t", "/", "50"}]], " ", "t"}], "+", RowBox[{"300", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"t", "/", "50"}]], " ", SuperscriptBox["t", "2"]}], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"t", "/", "50"}]], " ", SuperscriptBox["t", "3"]}]}], ")"}]}], ",", FractionBox[ RowBox[{"25", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "t"}], "/", "50"}]], " ", RowBox[{"(", RowBox[{"15000", "-", RowBox[{"15000", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"t", "/", "50"}]]}], "+", RowBox[{"100", " ", "t"}], "+", RowBox[{"600", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"t", "/", "50"}]], " ", "t"}], "+", RowBox[{"3", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"t", "/", "50"}]], " ", SuperscriptBox["t", "2"]}]}], ")"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{"100", "+", "t"}], ")"}], "2"]]}], "}"}]}]], "Input", CellChangeTimes->{{3.4149255582904053`*^9, 3.4149256044674053`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{ FractionBox["1", RowBox[{"4", " ", SuperscriptBox[ RowBox[{"(", RowBox[{"100", "+", "t"}], ")"}], "2"]}]], RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "t"}], "/", "50"}]], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1000000"}], "+", RowBox[{"1000000", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"t", "/", "50"}]]}], "-", RowBox[{"5000", " ", "t"}], "-", RowBox[{"15000", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"t", "/", "50"}]], " ", "t"}], "+", RowBox[{"300", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"t", "/", "50"}]], " ", SuperscriptBox["t", "2"]}], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"t", "/", "50"}]], " ", SuperscriptBox["t", "3"]}]}], ")"}]}]}], ",", FractionBox[ RowBox[{"25", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "t"}], "/", "50"}]], " ", RowBox[{"(", RowBox[{"15000", "-", RowBox[{"15000", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"t", "/", "50"}]]}], "+", RowBox[{"100", " ", "t"}], "+", RowBox[{"600", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"t", "/", "50"}]], " ", "t"}], "+", RowBox[{"3", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"t", "/", "50"}]], " ", SuperscriptBox["t", "2"]}]}], ")"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{"100", "+", "t"}], ")"}], "2"]]}], "}"}]], "Output", CellChangeTimes->{3.4149256059134054`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Plot", "[", RowBox[{ RowBox[{"{", RowBox[{"u1", ",", "u2"}], "}"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "500"}], "}"}]}], "]"}]], "Input", CellChangeTimes->{{3.4149256071844053`*^9, 3.4149256526864057`*^9}}], Cell[BoxData[ GraphicsBox[{{}, {}, {Hue[0.67, 0.6, 0.6], LineBox[CompressedData[" 1:eJwVyXk81Psex3GEY1/KlsmWZDDWUIdmvh/ZjrVDIkL2JVmGLFnPlaUpZZQb Wa6h0JE24tjLUpaOBmVrjHuocTiWmES2nJ8/3o/34/l4KfmEO/pzcXBwXMO2 9/PsbVuLABYpsvWoUnOEDzGy/Moxs9JOFHoxReZmpzyKE1wMOFj6HjFB48h2 DQkdN5N7QSUzkZA5LjXP2AEltFVwcJd+Qir7xF8N3fZEHz+SKhFlFnXiTIIt RgKR0fqIXQJ5ETmlxtk3K5JRP+V4PIOwgno+42r5zWLRTt5Vuw3aVzRAZMwU n0lGhAq6orTkGnqiL4HjUb2K3F/gVvUp39HolbTRlfxMlNUR2O24s4l+NNkJ TyVnoZaB2oII8g46nfxoR3SBihYmd0NvsXYRkzpe+GPwDrq1Vq44QuCEdCbJ +fPyXdTGzf76lbYPBCufnTwxVoxGBPOZXAxuWJPOeLtNp6Gl/aRucUleeNlW eKyaXobkla4XaFP4wDTOldWvVIkM8TpppE5+wA9+Pr0Q9Tuy1x4JtdsRgCzu 6jrdqUcoiXjYJIQsDOzAvBsWGc9QrlmPevwjEahN036PZz1H1TZhEhSWKPDl 4jXsHGrREoEZ9Z4gDtVhOSZCDvVI3ljNcoV2AMR3FQqJCy2oirTd/J0oAT/a eJb4atuQ4Sm6DgdDAkRmzhsTTV8he6toWVFJKZAbq4tMVetAyS7tiwSKDBRR pb+5wBskcD7XR//oQfAhPXd0PNuN7noGjhp3HoQQ3XpT54ge9MRfuN16Rxaq 6+64ujT0IWaUa24QWQ7KRvQ0rbLpKCiWwB8hIg96P90y7w8eQKvxu0mxj+Sh boDzpzGRQSSUWh6YwVKAiwFNG6neQ0i+J05uiKAEidLE/ib8MLIf6Vr+QlOG GyVqRRrGDJTyfw9XkTVlMMmaUOXDTaBns2vtBJsjwMnAvdnemUDiW/g7wd+O gNZ/8Iq4/kn0QeGmwSero5B+td3ig+UUcgs+Gz/CxoMTaeYulwML3Yhcml61 VIMHJLZI0RgLtSRk2hwoVoObJLxzsfcMkrvVeOhXS3UIKhMefhf7N/qrRu5l b6EG2EzHcte+mEOB26x9rWZaUHJ0Tr4hbAnl86RcYtzTAvexynejgl9Qj4jM 8MaSFpxrnxdFVV+QupJ1heE9bShwnXULi15G64ddjozF6MBTXKA6tXUFNVk2 v3F4qwvUr1srMiGrCHKuCppf1oeFSu9oA4NN1D0pReIq0oekUTl3MeomsidU Rbzs1IcYzV1Hn/lN5NE9MPzzAQOYjLz3Pej+ForfkaNp1RjAQ43jIssyO+hF cIO+zLIhEI3LV54KcQD+1KLHfMjP8NhoWJdlsA/Evzo/y/ElgXozvXX4kADI eyg5ZV4nweGHpS7jZwWA0D2/nvicBAbkD5682QLwS+FvKGiXBAf+m/0Pg1MQ kk2r3pEKEbDDKxPC5wRhrTd5LToYYOsjD1X4lTD4yg5JKJeaQMaDhoXFLjEg Ncc5JImZQYnYKbpipwTQotoeT5yzgs8PEmp1DOVh0oo5tGVhD2Gq89mIpQz+ NvGLPnQHCIxltTP91MCkaZ2bHuYE/RlLjb682iDYMMWBm3KG6StqQn7H9UBj xTLs5DVXCPkzlaJoaABGjbKviFzuULYkuakmfAKeWl1QeRPuCY3VR5rvMo2g VUDshzjFC+6WugZ5rxBBtqGC8azTGwSFHs5pXwQ4fHWL3CTjC94E8zS+bRPw 1OO/nRnoB5pN7YNGmaYgrmb3oKvEHzbrtq6vr5tBOPF0+CV2AKjE5UvVZVqA 2b1MQppWEGz9vpx+ftsSuvYnl/GmB0OHR7lHcLYVxEZoltkOX4Q/6tsjdGVt oFTyQC6X6CVg1JYxMii2EDtWT1L2DQXRYt6Acil76JqmMmppYaAgmO9tee00 mA0qc/GNhUObec1SV9yv0BDWFPEpNwIINYq2UicdgCKTPB6vR4YQy/MrL1cc QGeQ9j9LJhnq43vdN+ocYXmrBL+WFgmq9/3N2dFn4JJ10yN3pSj4bZcbF4x3 gtbKYze+90VB7x/c+NNzTtCH7/G1Db4MrxsNRfuKzsIK7mfOVO5osEmt/KBr 6wyLcvxbC9XRQKh6kPeF3wXKTVcXRa1jwJAwJRPa5AJ6cSRi0GoMxJ6R8tO8 cg4SBu6nx92JhZNuseGv1V2BWcC6fVE7DuxFXvMFjbmC5OWM7ZyxOMhL2Rw0 ynUDfvKpC+SYK/BSjs0tYnoehhI95DOU4mHhON8JLg53qNt0KkvsjYebpaoC fI3ucKYoUEEjIAGsu5+26ft7QLft+HcRsUTwMC24EaXgCWuKAW+rniSCodu4 zOKfnoAC6T546ySYUxkXo8VcAM+WPGcGOwl+kZ8WIkh7wQmDBFMLajIsKjdm 1F3ygtmW1qwZoxQ4VLDFNu71gqCNC+xNzPl11F0LzLMGXOdEjVNAX0FFyBHz 308slE9g7puwUQnCPEOjN1Iw53DmuORi/pw2NUM4mQJH2RvNC5inbHlMoogp 0HnZIq24zwu8KJUVmZiLzr6mPsT812sroSLMj9NMi2v3TLo12oU5sZ9Y37tn HekwKVIKDPGrz37DPCmhVtiIWUprysb+rRd4OrzleIe5g+10zhXz5M3QgGnM dTw9fn6Ymbw1uoIoBQKkqpPiMU98N+pxx7z/U8jTCszuBkxNMmb2wY/Nz/c6 OeVOOmZciWVPC2bGP52eTzDb5SpODWF2Uw3o6sDsyp21yMT80ZdPfRRzTcfa xuxep1Vlz2N2oXvxru71Cdtvu5ht1Pv272L+F1Th6DA= "]]}, {Hue[0.9060679774997897, 0.6, 0.6], LineBox[CompressedData[" 1:eJwdynk81Hkcx/ExjkeiHNFlS47sSFolG8K8hRFN+JVUrkJm8HANQq6piKxS 1OahtA1yLKXS8XBXiMpG2EKOdliiUo7c0n7tH5/H+/F8vD4q7gH7POk0Gs2e 3OJ+Gptnszj9JkkHrBmLG5RzQtc8s4bJ5j6wSU2vYYZLDXPWZLYydV6U9pen tDIjq3JpYpl9zLrCiB1up/uY39Ni984IxpkLMm798sfGmVViY+PjAlFsT1fK v8YQxfqdmpajghX4c2Kuo+OhPGze1o58FaiBGzcTkPReBfz3LoeXT6rhh7Zs Jk1DFXcHJ59u3qMOA/37G5i+qpCbY1zynlBHgoeRyJppVfytfF6vz0oDH/x0 h1jL1OHofSDi7RgDU6eaNB9ra4A73y9aab4FDIvkkNu7NwEpsVIWIdtRemfu 5tC0DuTGHe6meJgg0LCmY9ZiB0zKw6loWXOEpD7V4SQbQRBcdbvrkBXKFKLl hBVAj1V3yxzLBgNiYnSXkl3w3BMx7N5E4Vna6Bv6FXOYlk2JNfnbg1JW8CvT ZkGqREhTEjog0jKN/jnSElqjlv5GZw8jVMtWMql8NwxL1z4xpjvD96DZSGSN Fe5YHdlYF+CK0MSrXpwCa1QulV2QSzyKTMt6ra+Je7C2JLfzbo0bhD26pkHm bKjGzvHKVnsg/ArX8d4MG67bJFMTuMeQbrhlVCV9L+Q0996sveGJV8/ERZ2V bRBgbBvgO8YBv/hBsEy2DczTEzbHbfGCVVp+Qai8LWrlY7IkznjjMsvnmuJJ W4QFamex3/hgtrGqt0Zoi0zFFZfpMr44JyMnfXKtHcLaH5moefihtfpt1Qsz O9T2Xuy8L/CHfnGJzjquHcyb1ehL2gOwfMc6WckEO5T4lwX2XQ6E+Ldg+6+5 dkhcHdMRsY2HK1lO7seq7aDTLPjDspuHpKllqcJuO4zM3WBMxgVBt9spqGHS Dr7WZYXOKsEQiZS3cJKmUJmnmzT9MhixxolvFFQpvGQ892B7h2DOezR/dDuF USUDkdNix/HgVJ2uAYvC8DrJuc+3jqMwLOhG+UEKOWbfhmWsQ+GyL76FxqWw LdzE2OtbKK4zB6xjQilEvs4+E34pDF9MlQ4vjaPQfbU/1eeXcHDoE7SzFyko hsTPp7SHI+SD2hOz6xQkebuO8EJPIDtJFLx8Ci1RLuvjVSLgWPibf/t9Cg9n 7bOiXkRgPqUt/WgVhf0ZXGUtTiRexcZz8uop1LM7ppfLRsHTj53HaKYwuYHT UFAUBS1b51vnOigwuU3uDOtoKC80ZzsKKbhWpDl0jkWjWDquzW2Igr5epBnr YgzMJJb+FD1CYbCi8tyAIR/1OwWjf01R8Jo5MjZLLKzlNbwiHtSjH5LZyUf/ R9OcRuIPRSw1feJ7y3oOvSYeEDSVJhJXFIk/aSX+N044sNmIDz03+eR3xEK2 uGmwMR+PMlQ0B4mPJublJhBXP6+jDxH/88xKOoNYMdWn63+bJLfVEhfM377w adE6q/xXmvDxeKv65BfiHgXNa6XE9usbqyeIXakGWiOx4waXjMnFft6P00vc +fjj8SnibonirVJMPva3iWjOEHdNGz53Jg7yUr4wT+ys163NI/7dIdvr+2Ln 8S+dIbbq3bhrgbjzY41r0eJ/i8bED2LHnzm11cRffs1ppE1TeOexZFMbsYSB ar4IsaOg4MIn4qr266foi72LPfGDeNXK1U6ixP8B1XdvuA== "]]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesOrigin->{0, 0}, PlotRange->{{0, 500}, {0., 134.3748868893184}}, PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Scaled[0.02]}]], "Output", CellChangeTimes->{{3.4149256326624055`*^9, 3.4149256537714057`*^9}}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"euler2d", "[", RowBox[{"{", RowBox[{ "f_", ",", "g_", ",", "s1i_", ",", "s2i_", ",", "deltat_", ",", "n_"}], "}"}], "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{"k", ",", "s1", ",", "s2", ",", "f1", ",", "g1", ",", "w"}], "}"}], ",", RowBox[{ RowBox[{"k", "=", "0"}], ";", RowBox[{"w", "=", RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"k", " ", "deltat"}], ",", "s1i", ",", "s2i"}], "}"}], "}"}]}], ";", RowBox[{"t", "=", RowBox[{"k", " ", "deltat"}]}], ";", RowBox[{"s1", "=", "s1i"}], ";", RowBox[{"s2", "=", "s2i"}], ";", RowBox[{"While", "[", RowBox[{ RowBox[{"k", "<", "n"}], ",", RowBox[{ RowBox[{"k", "++"}], ";", RowBox[{"f1", "=", RowBox[{"f", "[", RowBox[{"{", RowBox[{"s1", ",", "s2", ",", "t"}], "}"}], "]"}]}], ";", RowBox[{"g1", "=", RowBox[{"g", "[", RowBox[{"{", RowBox[{"s1", ",", "s2", ",", "t"}], "}"}], "]"}]}], ";", RowBox[{"s1", "+=", RowBox[{"deltat", " ", "f1"}]}], ";", RowBox[{"s2", "+=", RowBox[{"deltat", " ", "g1"}]}], ";", RowBox[{"t", "+=", "deltat"}], ";", RowBox[{"AppendTo", "[", RowBox[{"w", ",", RowBox[{"{", RowBox[{ RowBox[{"k", " ", "deltat"}], ",", "s1", ",", "s2"}], "}"}]}], "]"}]}]}], "]"}], ";", "w"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.4149256996774054`*^9, 3.4149260041814055`*^9}, { 3.4149260668694057`*^9, 3.414926072263406*^9}, {3.4149262018424053`*^9, 3.4149262243454056`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"f", "=", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "3"}], " ", RowBox[{ RowBox[{"#", "[", RowBox[{"[", "1", "]"}], "]"}], "/", RowBox[{"(", RowBox[{"100", "+", RowBox[{"#", "[", RowBox[{"[", "3", "]"}], "]"}]}], ")"}]}]}], "+", RowBox[{ RowBox[{"#", "[", RowBox[{"[", "2", "]"}], "]"}], "/", "100"}]}], ")"}], "&"}]}]], "Input", CellChangeTimes->{{3.414926007150406*^9, 3.4149261132814054`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"-", FractionBox[ RowBox[{"3", " ", RowBox[{"#1", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}]}], RowBox[{"100", "+", RowBox[{ "#1", "\[LeftDoubleBracket]", "3", "\[RightDoubleBracket]"}]}]]}], "+", FractionBox[ RowBox[{"#1", "\[LeftDoubleBracket]", "2", "\[RightDoubleBracket]"}], "100"]}], "&"}]], "Output", CellChangeTimes->{3.414926114353406*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"g", "=", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"2", " ", RowBox[{ RowBox[{"#", "[", RowBox[{"[", "1", "]"}], "]"}], "/", RowBox[{"(", RowBox[{"100", "+", RowBox[{"#", "[", RowBox[{"[", "3", "]"}], "]"}]}], ")"}]}]}], "-", RowBox[{"2", " ", RowBox[{ RowBox[{"#", "[", RowBox[{"[", "2", "]"}], "]"}], "/", "100"}]}], "+", "1"}], ")"}], "&"}]}]], "Input", CellChangeTimes->{{3.4149261158694057`*^9, 3.4149261661894054`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ FractionBox[ RowBox[{"2", " ", RowBox[{"#1", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}]}], RowBox[{"100", "+", RowBox[{"#1", "\[LeftDoubleBracket]", "3", "\[RightDoubleBracket]"}]}]], "-", FractionBox[ RowBox[{"2", " ", RowBox[{"#1", "\[LeftDoubleBracket]", "2", "\[RightDoubleBracket]"}]}], "100"], "+", "1"}], "&"}]], "Output", CellChangeTimes->{3.4149261670974054`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"es", "=", RowBox[{"euler2d", "[", RowBox[{"{", RowBox[{"f", ",", "g", ",", "0", ",", "0", ",", "1.0", ",", "4"}], "}"}], "]"}]}]], "Input", CellChangeTimes->{{3.414926468294406*^9, 3.4149264813754053`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"1.`", ",", "0", ",", "1.`"}], "}"}], ",", RowBox[{"{", RowBox[{"2.`", ",", "0.01`", ",", "1.98`"}], "}"}], ",", RowBox[{"{", RowBox[{"3.`", ",", "0.02950588235294118`", ",", "2.9405960784313727`"}], "}"}], ",", RowBox[{"{", RowBox[{"4.`", ",", "0.05805244850561585`", ",", "3.882357086617171`"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.4149264833964057`*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"ee", "=", RowBox[{"euler2d", "[", RowBox[{"{", RowBox[{"f", ",", "g", ",", "0", ",", "0", ",", "1.0", ",", "500"}], "}"}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.414926179138406*^9, 3.4149261913354053`*^9}, { 3.414926235204406*^9, 3.4149262369344053`*^9}, {3.4149264118134055`*^9, 3.4149264120334053`*^9}, {3.4149264536024055`*^9, 3.414926457275406*^9}}], Cell[BoxData[""], "Input", CellChangeTimes->{{3.4149262463424053`*^9, 3.4149262900164056`*^9}, 3.4149264620264053`*^9}], Cell[BoxData[ RowBox[{ RowBox[{"eee", "=", RowBox[{"Map", "[", RowBox[{ RowBox[{ RowBox[{"Drop", "[", RowBox[{"#", ",", "1"}], "]"}], "&"}], ",", "ee"}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.414926293077406*^9, 3.4149263090504055`*^9}, 3.4149264232794056`*^9}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ListPlot", "[", RowBox[{"{", RowBox[{ RowBox[{"Map", "[", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"#", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"#", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], "&"}], ",", "ee"}], "]"}], ",", RowBox[{"Map", "[", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"#", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"#", "[", RowBox[{"[", "3", "]"}], "]"}]}], "}"}], "&"}], ",", "ee"}], "]"}]}], "}"}], "]"}]], "Input", CellChangeTimes->{{3.4149263280094056`*^9, 3.4149263896424055`*^9}}], Cell[BoxData[ GraphicsBox[{ {Hue[0.67, 0.6, 0.6], PointBox[CompressedData[" 1:eJxtz3k4lQnYx/GHJpUoLVKUTkKFUNEi9CNCmxNly3LsO8fh2A+nUtpDO41O q0KlrM2YHNWUNBWpmSTNibRIZaJoEm/vNffzz3u9559zfZ77uZ/7+s7wi3YK lGcY5rMcw/zv////6172fx4gZeJFuxfJGfR8JF4FfrOsX5hPVsHeg9kfla5e IE/E/lX2nzR9qsiToXO/SMnkey15KrIXpbiE19eROdhgsDFhv20DWQu/bH1x Uy73EVkb0y+1t2l6NZN1Yf79+6umHc/Js1G7OaO7QLmNrIcpYr+yTQodZAP4 +k9YtnbuG7Ihyl7d9ZATvyMb46bOqMfnBz6Q5yHre9u7Su4/5PnQ5c8qcVjY Q16AqU8uh5gafiaboOimzYOSRX1kUxw5uDbz27qv5IUolO/e4Z36jbwIfWrd efoV38mLsXzU339UDw6Rl+BMsufQVT05/GczdA6VmrV6ypOX4sHkOwvFB4eR zbG7UtXL+9FPZAvMcp1lsnCKAtkSFr8JTe74jyAvQ9XHD8dbSkeSAYO9IzxN Rin+ZzGwJlPr/jH/0TS3QtdR3V9k15VoboXTFhoz7+uOobk1diWen2u6byzN rRHz1ki+V6RC8+W4su2et9fCcTRfjqqalBO13WTGBro7J/+sdGE8zW0wqb24 Ty98As1tsVhjNkdRfyLNbdH4oMTmSBeZWYEl7zc53LykSvMVGLq1UiqKm0Rz O/zaM0pWZqZGczsMM94o8pCfTHN7jJg0bxjvLllsj0dfesIrD06huQO0jwyZ Ofuq09wB1xSS43UNNWi+EuvHp53QHiCLV8Im8ECF1d2pNF8F803xWoK8aTRf BVUztZrCCE2ar8bBnK77MsvpNF+Nb9lLB34axaH5GuyqlpsxJpYsXoMvTKNC Zis7XwuxxXnFcPsZNF8L3a5LpeWlZMYRG77GXfabrkVzR9yYt+Zp5C4yw8W3 cbVzmvrI4GLL4zE1ewNm0vtcCMqvtRxrJEu50HO1b/tqqU376xCSUrf+ZDEZ 65CrdKspR12H9tfhXdXwaXXbydJ1MF8Xrrisj8w4Ydr8N1m9gbq074Q+hfia jkdksRNeW027NtZmFu074WSxwr3oUjLjjJPph5cPmzmb9p2xXXh43+0cstgZ MW48k1/k59C+M3btnVT/REBm1v+43znAaSdjPVrseLuznPVofz1Ci0+t1LlJ lq7HHddWjVYTfdrfgAe5+z+XnyFjA24OjX1+YZIB7W/AgtS7765nkqUbwD0r s+npJzMuOHx6i6ZF2Fzad4FtTnp1fgtZ7IK+SVf5k9YY0r4LuD3luSevkRlX lNnNKbAyNqJ9VwR4ZYz9dIIsdsXu578bNI8ypn1XiIwTd0mtyYwbVoXGCy6m kOGGu8N6TE+WkcVuWPs20yjvPbvvhoc+iteO6M6jfXeU5/Jsc33IcEdui5vh 8SNksTt0xm9+dbaRLHWHcudwWYnifNr3QP2ijOzq5WR44HPZR+GdVLLYA4aD 4Up/lZOlHvh+0+DJyw/s/kY0z3bhfZq1gPY3YrL+c6MhHlm8ERrzox8r5pKl G3Eua1X7xCYy44mHCvM/TlUyoX1P+M9fnqNlSxZ7QrvjdYhOGlnqiQ7XqvE6 lWTGC0nzAubP6Gb3veDinhc5ZY4p7Xvh1Z6+CGU/stQLHRWvWgZyyYw39ncd DHndRIY36gxab/yhtJD2vbGiUv9EsS1Z6o0BRWHFtjQy44OZaZrVHpVk+CA2 /H3E7G523wftHz5v/Th7Ee374IzBzluXfMkMD87T9LpDcskcHiKfxt5SbyKD h/SgiWN/H734P/N42NlVlR1sQxbz4CFgJsuLyBIern2NjDxQTpbyUKG00Ufz A1nGw+OSzurjukvovi8WXxtwmexD5vji3KEZzPbDZPhiz5K9uz88IPN8sf2y ysNVI83ovi/OPvnp3HGQJb6w4Du8fptIlvri1A7TML3LZJkvFC32Kvm+JTN+ 6ChKLdozYynd94P36gydEncy/BCekuRRl03m+cHKOVH/zztksR9ab7QkN8uZ 030/7HtnYvBwCVnqh4gL1trSGLLMDxVb3tqfOk9m/OG4W39Tygsyxx8+9zjl 9lMs6L4/Iks+3lNcR+b5Y8mV38trt5PF/jhv2+AaKSVL/LH8jalEuZ8s9UdB 0NjUk0aWdN8fexx9m/WDyUwAXu1XLTifT+YEoNrp8yONP8kIwOYHJ503Ky+j +wHYl71/YqsNWRyA7UWKigapZEkA9syZMCu6lCwNwPAzyT5nO8myAMQOdjS9 UQHdD8Svk88kJxqQOYFYePmzprw9GYGw6ZtTvdWfzAtEffObNcPSyeJAfKvQ b0jOJUsCEalhvqyrnCwNhErW5cOujWRZIPLjo/+s7mLvB2F878l+jZFWdD8I 995Iv8XNJCMIzJctrbctybwgDO7rzJvoQRYHYeo0U9ONQrIkCHWclwV5WWRp EFIFFZ8eF5FlQYC4YKLibTITjNBj75UWt7H3gwGX/qc+39n7wXBtjUrbNNma 7gfDINqw79gCsjgYb+7b215ZS5YE4/Ut+dDaULI0GPdLH/rVZ5Blweic7LPg 3nEyEwKJhajlzi9kTggKmow8pY/JCMGTM2cqSrrZ+yHwdlzclTt6Od0PwbBv YYNpumRJCPh6bz94WpGlIdDZurXGxJMsC0H/xcUChQQyE4rXOpxRD7PJnFBM nCXadriYjFCU5ZZ2udwm80JxV8HBUqWNvR+KzdvbRTcG2Puh6G2yv8hXs6H7 obAy83qoNp8sC8Xv5Vrvq1aTmTDYDVTKrQ8mc8IweMh1wttNZIRhOeOqk3iM zAvDvVs6S+UqyeIwuFasdstoJEvCsHt4pFiui70fhvE7llUmKtjS/TBI3u1n OjlkJhyfOx15G5aSOeE4eDTq6S8byAiH0ZkhgTqfzAtH1gldw7idZHE45N3U VepOkyXhYDS0p02qIUvDoat53N27mb0fDrOHfXckPez9CMC2SNCqvILuRyBx 2wq3ibPJiIDJmKB0W2syLwIhc8q7+J5kcQRkIe9PH4onSyKgPPHP4oossjQC tseXj2osJMsiEKmjVtNxk8xEovelRVPvc/Z+JEZKS1YO9rP3I9FfE6snP8GO 7keiOO9gAjOXLI7EdCujRV9XkCWRqPf1j+rikaU/vi/nrfU0mSyLxF3F9W7X D5CZKFj+vVn1zEUyJwrbeox8NteREQWH4HNmHm3s/SisOYkLBgPs/SgU9ThX f1W1p/tR4C+wiaw1IkujcDIo6fctDmRZFA4E+tbDn8xEY1tG9Kb+VDInGgmW U7sLD5ERDZ22lxpuJWReNAytViky9WRxNMx9H0tPtrP3o2FhpGyL7+z9aGSZ iI89meRA96Ohvbbmj3BjMsNHXKB3y1cHsgof/vq99zb7kzl8tD6aUqggIhvz sbQ5TLj1EBl8BF7JMB26RObykZev0B17h8zjo8M/u7C9jcznQ131dtjaAbKY D/GTlQvLVFf+5yw+3lySjlc1Ikv4WG99b3iMPbmEj1EFqqp1vmQpH2rzl61Q TyE38GFt/TE/+ABZxoetyl3tkgvkbj746QmtPbfITAyyPh65P19GVomBhkrT 94ivZE4Mfi6r458cv4r6Y+BRqDG3SZ+MGMzYF7F4yIbMjcGK8x5HdL3JvBjU SLc4OSSQ+THw+rolIjiLLI6B2sXP78TnyVkxCF4T8+fB62RJDC4t3DPvbAu5 JAZ5R1rlr/SSpTHYzu11vKq8mvpj0FzP06zWJctiYP3yUuwvy8jdMRi8G76h zI3MCPD2lcad8zFkFQFSS20e5e4kcwQ4eTE4NfMU2ViA0CXDbkdXkyGAKD3n qvNjMleARO897gs+kHkCfFoTUjxmxBrqF8DN/kBlx3SyWIAtZWniqsXkLAFG cbOVMteRJQIsU5/vuy6MXCLA52U3tqhtIUsFeDriVFJzHrlBgJ01+vaHy8gy AWzOl/Zy75G7BbDq3J2h8IrMxOJ80vDvlYNklVjMUAn2D1RbS/2x+KTWVz3G mGwci9qfv4wpsycjFs1WJ703+JK5sZDTm3/5UxKZ92Pe0zx6dw6ZH4trBXWx M4rI4ljMK5794coNclYspFeni/GMLIlF0IP7hvW95JJYhMxyGnJUdqT+WBQZ 1H5p1CE3xCKBb6jOtSTLYlEqKYiqdyF3x2JnmulXRJOZOKQ9brlWmklWiUO7 4OhNLQmZE4faaSFj91aRjeOgsc/89OcGMuJwOExhu/tbMjcO1zWe/P2M4VJ/ HDgXxs4erUzmx2FAQf34oilkcRy8LgXY+emQs+LwOdtj5q55ZEkcxN/9Ta9Y kEvisObEq/S/HMjSOPTstB3+bQO5IQ5ztF/cnupHlsXhprFavXkUuTsOz/6w HrMxmcwIoRpefShhG1lFiNFNzYE5OWSOELLpH0RF+WRjId5dc3hxvZAMIXZE meU8qSBzhRhjMXp/13W2X4iHmT+9GrzP9gvB+zd8z9gWtl+I2j0ZuzRfs/1C tGfubdPvYfuFmJpy+8CiIbZfiDu2YaetRq+jfiFGjL80YaUauUGIzA13X3Jn kmVCJDX+M8XFiNwtxL4cp3L3pWQmHur/qF/daEdWiQeyhbM8ncmceFj575Tf 6EM2jsfTU2mr3cLJiMfVZ3Ej1ieQufFQ6NhtunYLmRcPubb2Jyv2kfnxKLc7 0GORRxb/2G8uz1xQQM6Kh2Kx1+FZpWRJPD7NzNdWryGXxGPY123zRt9l++Oh tVBH+u+fbH88Tphtefimje2PR9Sv58Ief2D747H6Zu5O6b9sfwKmXd1oUqjg RP0JWOr6LjRnPJmTAFnFirlJmmTjBBxdHSL21iMjAfWddoHWC8ncH/s3nrZq W5N5CSgPmdA1fC2Zn4DChhfZHe5kcQI61i9ovBFIzkrAmvJPRZIYsiQBRVNG aqeKyCUJ2DzBfanLDrI0AZ3H63oMD5IbEpA6Z5GjwgmyLAFWOzc7Pismdyfg svK2LyVVZCYRZQL1FRk32f5EXH6vCpcGtj8RZ/RWv9J9xvYnwsJxr9mX12x/ IoYGCi1v9rD9iYhzjurNGmL7E+Hz6YSX52hn6k/El29KqbpqZHEitNXi1n3U Imcl4mHM0b8rDMmSRJjeMJ8tMiOXJEJrrcZc6xVkaSJC2z9+HO5EbkjEgpcC QZ0XWZaIYBuPih2h5O5EqO5x/81BSGaS0P58xo6Rm8gqSdjU66txazeZkwSO zi3h5iNk4yS8D+o/anGajCTkrJPs6LtE5iahiuO/ouRXMi8JeyvaG4Nvs/1J cM6/Plezie1PQseY/I1Nz9n+JKx7P84rs5PtT0LAg8oFZl/Y/iTYDuq+eCe3 nvqToJjUG3RMmdyQhA19V26tmkKWJaH73Se5f7XJ3UmIzrRWLzAmM8m498vc cc7mZJVkTDCzfPPdjsxJxte0b5ICZ7JxMvbd7TTn+pCRjHaDfb/1hZG5ySgU hWrnx5N5yegM6Rcs30zmJ+N0RFXh6z1kcTI2PZt5b+dRclYyRC+vtBicIUuS cXBXy1/3SsglyYirGHcjsprtT8YMvdv5SnVsfzIu7g2PKGxi+5PRERMx1+5v tv9Hj4Npe1sn25+CeYlG+0Rf2P4U5K66Pl9NfgP1p6BsmPODS8pk4xRoNtkE 2U0hIwV2bf39rdpkbgqmON3fHmtM5qWgwFJLbaQ5mZ+CJ7WOZ/LsyOIUjP9y ycTQmZyVgtf9R2/VeJMlKah+7+fJDSOXpCCW8e77W0iWpmBiyLMjUZvIDSkQ hM61HthNlqWg0T2/d/sRcncKdp9IvjzxNJlJReW10cnHL5FVUnHd7qjjnF/Z /lSofQhZcOUW25+K38TXdM0esv2p6P3QolfbyvanomGlnJXdW7Y/FY/vRoT9 0cv2pyK70eMcl3Gh/lTMezV8oGk0OSsVk2vPBrmokSWpiHWIe/uXFrkkFfVH Tm11MyRLf3zve7D5kyXkhh893z+NcbUly3583z1u8BGX3J2KDCV1JSdPMiOC o7PKkvvBZBURGgo3bXOIJXNECH106PONNLKxCFpP4jMtdpIhwuwRdhYVB8lc EVYP6k02PEHmiZB2zmbKmWIyXwQ7qzIrjSq2XwRL1yM5WTfYfhHS7ysrD3/A 9oswrd3oSuJTtl+EbP3JW951sP0/7uc9F3v9w/aLwM84Wnx/gO0XIT/Pa8Sy ka7UL0Jfv0XOxQlkJg0LNFavnTadrJKG528Pm+/SI3PSUFVp5N1vSjZOwzhb TlmAFRlp0L4et6xhNZmbhsjGhXJL3ci8NAzuiPr3tD+Zn4amJC2DMdFkcRqm n3M9HJ9MzkrDscPjLJ5vJUt+vG/vpGmbTS5JQ7C5unnRMbI0DS7CkMMq58gN abBKMZsnLCXL0nDWc7di8zW2Pw1al920zevZ/nQU5h1Jz3/M9qfj+h/LpzEv 2P50lGZwB3ldbH86ag6WzqztY/vTwQ3i7eUMc8P/AET+9+Y= "]]}, {Hue[0.9060679774997897, 0.6, 0.6], PointBox[CompressedData[" 1:eJxt13k0FHrcx/GpZCllpLqSarK0UF2EQvRJSERIlCXG2NdZ7GY0IsqNVEja VLdoIa03t220KFdJivZLqzbltqrcenrO/X3/eM55/OO8zu985zPv/xgfkuAV 1p/D4Xzsx+H87+///6dn9v/9zUHNvIfpw2t+MKsi8u8ts3xzVfCfuRjuKbHR 1NBgHo5LA/Tyjay0mLUxcuilvE23RzLrYpbDhBu+4TrMPNxZoTL60aMxzHq4 gQEvfSt5zAY4YKd5ftRiPeYJsAyIzO/+qs88CcJt/fqpbjJkNoL6Bsf2zb9O ZJ6CqS46F++cmcQ8DW+0vtT0OBkxm0BiPv3YkMvGzKa4ELTaI8B+KrMZ6s9Y Rigfn8Y8Hc9PNpiUa5owmyNLu2Kb8kWyBW7fO7HQP8mU2RJ737WI9+mZMc+A NL2hfOBV8kz84Ks2yBKnM1tB9L44TH+UObM1gnLXCJRPkm3w8E3QZEs/C+ZZ uFu5Y3f9J7ItZhnJPf8osmS2Q53grYXhpBnMs7H2r9we7dNkYOj4USe2L5z5 n+XAik2D+A0dzJw5mHEoir8pzoq9z0F9Q/NNgy/MHHs81leWJ66wZu/2uNqd d3ODmg17nwuXafXvCwqZ5XNR2G1yOE5zFnt3QN+YgHLb9cxyB8x4Uf5WRdOW vTvi9Nn62vYCZrkjfudXlR9RtWPvTnj1tsioKotZ7oSQf4QVp3vpfR7+Gj62 379xs9n7PNT4BWsnPGTmOCO2PMm2cDLYuzN88sZdyOUzc+Zj/QJry91l9D4f emu/FHU307sLElaurhMozWHvLrj7pF022IqZ4wppn1vJ81h6d4VW2rod37fT +wI4OU3+1bOV3hdA6+vk448H2LN3N8hN2srrzJnlbrCaMX1aeyi9u2OAYE/N rGJ6d4e7INW55xy9L4TZY/7IDz30vhDO59T1XMfOZe8eEBivC+11YYYHtJq6 G74nM8s98EyQocffyazwQMtX/zl6V+neE8dan/7t8JnuPfEselBVE8+B3XtC fKjOpm4+s8IT4RvnGHPFzBwvJJvt+9q6iRleMCgZOl6pnu69YJHzTnVPF917 YdS0Eu7JIY7sfhG6lOfkO05nxiI0vCg7NW8Js3wRjm579Oi8lFmxCMMfuU89 VUH33mhQm11vdoHuvTEtZFXd5C6694ZGwLEJu9Sc2L03Pu100NhhzMxZjJ1m zxIM3ZixGMfkx23M4pnli/GuOSWgsZDuF2NmXePh1zV07wOZZtDIbc107wP9 0WY+Hd1074OJln/MPKw+j9374KHaUKm2MTPHF3Z2pu1a85nhi0PXDT7uDmeW +yJ8lP7ylmy690VF6uVhJRV0vwRr7tnO+3iK7pfAuV/IgXe36f6nK/75q/AD 3S9BS0KNQZOGM7tfin28leZ7jZixFA/cD8cbOTLLl6Iq6Rk8gpgVS1HzJnKY bhrd+6FRydph7Xq690ONnrnp/v1074exhlOmiS7QvR/G3ObmP7pP9/5w3aNW yPlI9/4wyyqWXVGfz+79MUJtV+k8A2aFP7K/HteU2DBzAlDd9e94Dy9mBMB1 buzHjki6D8Bv524d1V9O9wGIHK2Zo1tC94FIMPg3q2kf3Qei26/6sqmC7gNh cDMlzqON7gOhs21o1qSXdL8MvsWlGie/0/0ynD1wmqem5cLul2GNouUPrYnM imWY19Zy4541MycIshd5CSHuzAjCBeetWVV8uv/5Hr23ryaR7oNQt2z+teQ8 ug8G/4rpE045My8YTdlz9NwP0OcFo6ZqQ2zQGebgYPgkr9tp0UKfH4yuigEV 1x8yVwTDy0ff1vw97QVj+8N97suUXP9zZzA475/u9hjBzOFDf209V30CM4+P YVX3lq6zZAYfD/nrnR86MQfzwR+otaPPh1nOx9udkQZPwpkr+IiMyC8qS2ZW 8GG3O6pCJ5f2+di64cfgmBLaD4Hoot+q/N9pPwRj/Q5cSjtC+yHY+8I7zeoc 7Yfg09R/5v/VQvshuBId/nhKB+2HwLj0zWt+N+2HIN3eYHD0N9oPwYgFQzsd 1RawfQG2uRoP/jSSmSdAbrN0coYBMwTg+a5vazNlDhag2LC4Unk2s1yAL3+5 WQxdwFwhwI28mcO6lzArBDBJWXltVxhzpwD5eRP7zMS0H4oN9g8MNmfSfii4 YdY37+fTfiguzTu0vreU9kOxTWfLszc7aT8Uj985eitqaD8UZ6vfB8b9Sfuh aK18nf3pIu2H4o5LtkPAddoPg794jO72+7QfBtOqWSvPdtF+GK5Z5faefkf7 YTCY2PS6/F/aD0PXK41bPqpubD8MIw2Kjd4NY1aEIVn97pK4McydYWjtqL5+ eSIzJxynNKcOGGjGzAtHTsTVpXqzmBEO3R4BxjsxB4ej5J4xr78HszwcT82u Bl9cSvvh2C4URMYIaD8clt+yN3+Ipf1wJGU5BPKTaT8C56zmaBxeTvsRaFDa pdm1ivYjoJZoeJ6znvYjMPqYbtaPctqPwFV+RuXjXbQfAXXL20XVB2g/Ah37 JaWBx2g/ApqZttwPp2k/Et8173qLGmg/EoMbK+vammk/ErOsL6wff4v2IxHq rTl9cQftRyIoyvq1sIv2I3FVr25o0lvaj8TZ4Lsv+J9pPxIn+p9vmvGD9qMg 35zY77OyO9uPgi63tX37UGZEIUTv91KzkczBUcjp2pp/cAyzPAqfjfjfRhoy V0Rh2JkVJlFTmBVRuHS0ULhnOnNnFFrWqve7Zs3MicYe7+wJj+fQfjR2bJAZ PHam/WisCSiwvbaQ9qOhrudxpNKH9qNhsEbcEBtI+9GYbH30+NhQ2o9GeuKr c6ejaT8as5pVJ7uIaD8Gkj90JlxIof0YnPrh8ck4k/Zj8Ho+py0rh/Zj8Ey4 u+9iPu3HoOPB8z29RbQfg5CWHeo6G2k/Bu5ld8KMttJ+DLjnml9M3kX7sfC3 6G3W3kv7sdjspOr6uYb2YyE4+Xb9xaO0HwuNc513s/+k/Vi03/f0MVPQfixW Kjgzr1+k/ViMmfjxUHAT7ccivzT7Y2cL7ceherTHdK922o/D0y8em4/eo/04 +C28HqD6kPbjMKVGssftGe3H4dI+8z05r2g/DiuD4lZX99B+HK4/bZNe/kj7 cfh8Ie9A21faj8ctH45D2w/aj8dWPf34S0oL2X486qbOWnxAjTk4HpZz++vn DGWWx2Ob6SbOQi3mip/+sWHkEG1mRTxW/vJ5zRld5s543Iyy/k0wnpmTgH9z XWf2GTLzEpBRu/bPVUa0n4CEh09nqP1K+wm411zVJptO+wmo7ag98WwG7Seg XtdnoMMs2k9Aq7a4owS0n4AbLfzo+w60L0RqjkrjL/OZuULMy7syxtmNvo8Q GLQrP86T2UQIu79VTFctpu8nhNmvSmZlS5k9hFgy89earYH0fYVY8c/Bi2V8 ZqEQG4/Zla0Oo+8vRPijCY7CKOYiIV6InF8tiKMeIYYVjNo7TsRcK4SDdtWW 54nUJ8Qqq82PK1OZW4SovG2+OVBKvUJYjVrUqSZn7hGioVl07UA29Ysw/NX+ HKc86hfh5LNTo2/lU78IXIP+tcsKqV+EwgXSgAfrqF+EhqoWm0Ul1C/CJ7FL pKKM+kXwNbj/zWAL9Ytgf5rHzdpO/SJYqQaeubmT+kV4c0s+jreH+kU//6Kd tECwl/pFUGraErDtAPWLINzTHnn9IPWL0NsiW9V3mPpFeG/3qnXcceoXYVnf 5CCbOuoXY7BL79yFp6hfjITjlzf6n6V+McaYnRcGn6N+MQbyGp4EXqR+MbL5 EWrel6lfjMXLPr63b6J+MVaqptYZNVO/GP73D8UPvk79YoTf5016eoP6xTiw +seXP9qpXwzN24692XeoXwzVcfss59+nfjFKlHdeVumgfjFOCgX1Zx9S/8+9 oKfTRU+oX4xN7Z1mul3UL0HLyestihfUL8GkQ5/HBb2mfgnudvuYfX5D/RJM 06obtfof6pfA6Xt11/AP1C/Bn6Vrqss/Ub8EOuUSqc4X6pcgJNQlrPgb9UsQ daw+VfU79Uuw32qfwObn/+X/9UswyuL32oR+zLUSOOf4Gu/oz6yQ4Lj6kXvX BjC3SLDddk3jNyXmTgkmR2/vM1Bm7pHAUq9I6qrCzEnEk9UvPBJUmbmJOHtw bFaRGjMvEdUP84cfHMRskogijnhY02BmJOJIK2fFE3Vmj0T0ertFfxvCHJyI 64N8WzU0mIWJ0HBTqR/PZZYnInujqa2pJnNRIlLbNrjZDaP+RHCr23udtag/ EfLGAkfP4dSfiJkGM619R1B/Ij4uwAO/kdSfiPraPvOAX6g/ET0dB+f5a1N/ EgbfWDV6ySjqT8L5abuOe+lQfxIe5KqMcx1N/UmYVLHfb44u9Sch/ZA4wWIM 9SehbtdUwcSx1J+Erx6bbX4ZR/1JOHAu8rsSj/qT8CTe9FAPuSgJpxOXL707 nvqTcP9WU1+9HvUn4ZeTmRWV+tSfhNnlb+f/ZkD9SRiRHfY91pD6kxDn9sv5 BROoPwnLW423Gk2k/mQk6P+5TnkS9SdD+Wzj7k4yLxl7G7z+PjGZ+pOxxswO a42oPxl6A8NbBMbUn4wbspINllOoPxkfXVYXqkyl/mRwT7692E6WJyOqO9Lu 92nUn4yOfnl9Cb9SfzK0nZuVrE2oPxk/4o759zel/mRsqt49sJHckgzzr1+V C82oPxnvpiiFek6n/mRUmwwYr2VO/SmwmzrV/gaZm4LJh/c1rbOg/hTs69xY 725J/T/fS5QmDZpB/SnQLfrQ/wLZIwVe24J8pDOp/+e9wnTCdCvqT8HsRa7J z8nyFPQ3z3LZYk39KWh+vmWPuw31p+D1hYDCH+TaFIwYmjjw4CzqT8Eb3m8a gbbUn4LvxbNr1eyoPwWGwQavjpJ7UlD6pf78stnUnwp3aamtCqg/FSMaRwTU kHmp6F+1V897DvWnYtWLQ8W9ZKRCIy/n2GZ76k/FyPNxq2znUn8qHjvXcv8m C1Nx3nWlp8yB+lNxf7eu+2hH6k9F4actaifIFamIVejnLXKi/lQ05l5s7CYr UnHvZcHN3HnUn4qMy8v3jXWm/lSkeZW4HyP3pGLDNMV5l/nUnwaldS/UO8jc NBz3fTdF7EL9aeAaHpmg5Er9aYhM7+srJiMNo1tLDxssoP402NfZuxwhB6dh 5IsTDXCj/jS8C9pndJUsT8PAdEXqEnfqT8Neta1HHpEr0hCRcPdBzELqT0Pl ocFf35MVaRhSfGpQhgf1pyG3RjSsnyf1p6FfTpVWLrknDZ1lt7iDvKg/HcNW x6sXkLnpeDmmXVVjEfWnY5dAqFpENknHljNlQ7je1J8ON87vOmvJHunoyeg0 HbKY+tPxRKd0cT5ZmI4mofIqZR/qT0fItdimLHLRz+/j9318H7kiHXNHvSxI 8qX+dEQNWar5hqxIR+uc+JqwJdSfDtMOp4j75M50vLD9auu1lPrTwevdZH6J zMlAwTdTNxs/6s/A4M2XV9eQeRl4v33Za54/9Wdgpv/b5PVkZGDBcaHxgADq z0DD3BuDJeTgDLQt7K/ziCzMQEbPY1+PQOr/uT9i6aXT5KIMLC+zjTZaRv0Z +Nc10r6UXPtzb/B+z35B1J+BtMbWTTHklgx4V1TqtpE7M2DV8fqebTD1Z6By FP/ObjJHikqLMyOH8Klfin7pdSUSMk+KuVsGLLlDNpFiU7NzkF0I9UuRfNWk eifZQ4qh3cvmKguo/+fnTxDrR5GFUvxpOHBhE1kuxUvvmktTQqlfimsvfxQW kCuk2Pgpb083uVYKWWyLhlsY9Usxyizoyn5yixRZRscfqIVTvxQu9oscI8g9 Uvhl1CpdIHNkKOudO4YXQf0ySNTkRRlkngxOFjeC2skmMhQ5Va0xiaR+GTR8 wrXzyR4yHG4TKj0mB8sQuXmil00U9ctw5dWz/hvIchmqc1S0X5KLZEiq/Hs9 oqlfhjj9M4ml5FoZHgZ9OPOKrJCBF35FjBjqlyH4VcG6YnKnDJl+Mt3n5B4Z tj64N9wmlvozUR3RlVFA5mbiyum7Xh1kXiby5d3lJnHUn4lhFt6Ls8jIhIU7 sq+TPTLhZdFpOD6e+jNhtCjITkgWZmKd9rOrZ8jyTKiolLepJ1B/Jszqtvr4 kSsygVdcr0pybSaMs7mN78mKTHSVHT81W0j9mRBl6kz7jdyZiYZLHmPbyT2Z OO8VU8gTUf9yPHdPWx5N5i7HQMe8d0fIvOVoHrv9eR/ZZDl8pl0OcRRT/3JI LFRCCn76fwCwL70l "]]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, PlotRange->Automatic, PlotRangeClipping->True]], "Output", CellChangeTimes->{3.414926390591406*^9, 3.4149264330474052`*^9}] }, Open ]] }, WindowSize->{607, 750}, WindowMargins->{{447, Automatic}, {Automatic, 83}}, FrontEndVersion->"6.0 for Microsoft Windows (32-bit) (February 7, 2008)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[590, 23, 1194, 37, 72, "Input"], Cell[1787, 62, 2007, 57, 114, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[3831, 124, 1869, 54, 135, "Input"], Cell[5703, 180, 1759, 52, 114, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[7499, 237, 262, 7, 31, "Input"], Cell[7764, 246, 5102, 90, 242, "Output"] }, Open ]], Cell[12881, 339, 1750, 50, 132, "Input"], Cell[CellGroupData[{ Cell[14656, 393, 538, 18, 31, "Input"], Cell[15197, 413, 466, 15, 47, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[15700, 433, 557, 18, 31, "Input"], Cell[16260, 453, 467, 13, 47, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[16764, 471, 247, 6, 31, "Input"], Cell[17014, 479, 543, 15, 50, "Output"] }, Open ]], Cell[17572, 497, 424, 9, 31, "Input"], Cell[17999, 508, 124, 2, 31, "Input"], Cell[18126, 512, 309, 10, 31, "Input"], Cell[CellGroupData[{ Cell[18460, 526, 729, 24, 31, "Input"], Cell[19192, 552, 13313, 223, 242, "Output"] }, Open ]] } ] *) (* End of internal cache information *)