(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 11.1' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 158, 7] NotebookDataLength[ 47752, 1012] NotebookOptionsPosition[ 46295, 982] NotebookOutlinePosition[ 46645, 997] CellTagsIndexPosition[ 46602, 994] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[{ RowBox[{"ClearAll", ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"BackwardEuler", "=", RowBox[{"{", RowBox[{"\"\\"", ",", RowBox[{"Method", "\[Rule]", RowBox[{"{", RowBox[{"\"\\"", ",", RowBox[{ "\"\\"", "\[Rule]", "\"\\""}], ",", RowBox[{"\"\\"", "\[Rule]", "1"}], ",", RowBox[{"\"\\"", "\[Rule]", RowBox[{"{", RowBox[{"\"\\"", ",", RowBox[{"AccuracyGoal", "\[Rule]", "MachinePrecision"}], ",", RowBox[{"PrecisionGoal", "\[Rule]", "MachinePrecision"}], ",", RowBox[{"\"\\"", "\[Rule]", "1"}]}], "}"}]}]}], "}"}]}]}], "}"}]}], ";"}]}], "Input", CellChangeTimes->{{3.741965896356278*^9, 3.741965896357958*^9}, { 3.741965934719986*^9, 3.741965950254174*^9}, {3.741966143168725*^9, 3.7419661455969543`*^9}, 3.741966417382311*^9, {3.7419665314989023`*^9, 3.7419665373124647`*^9}, {3.741966787610839*^9, 3.7419667889026814`*^9}},ExpressionUUID->"33eab7b7-2a12-40bd-9558-\ 324e8e2e4bd0"], Cell[BoxData[""], "Input", CellChangeTimes->{{3.741966074622986*^9, 3.7419660787345867`*^9}, 3.7419665339914837`*^9},ExpressionUUID->"01917e4f-2166-4ab9-bb29-\ a631abc026bb"], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"solb1", "=", RowBox[{"NDSolve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{"y", "'"}], "[", "x", "]"}], "\[Equal]", RowBox[{ RowBox[{ RowBox[{"y", "[", "x", "]"}], "/", "x"}], "-", RowBox[{"1", "*", RowBox[{"Sqrt", "[", RowBox[{"1", "+", RowBox[{ RowBox[{ RowBox[{"y", "[", "x", "]"}], "^", "2"}], "/", RowBox[{"x", "^", "2"}]}]}], "]"}]}]}]}], ",", RowBox[{ RowBox[{"y", "[", "5", "]"}], "\[Equal]", "0"}]}], "}"}], ",", RowBox[{"y", "[", "x", "]"}], ",", RowBox[{"{", RowBox[{"x", ",", "5", ",", ".00001"}], "}"}], ",", RowBox[{"Method", "\[Rule]", "BackwardEuler"}], ",", RowBox[{"\"\\"", "\[Rule]", ".01"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"solb2", "=", RowBox[{"NDSolve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{"y", "'"}], "[", "x", "]"}], "\[Equal]", RowBox[{ RowBox[{ RowBox[{"y", "[", "x", "]"}], "/", "x"}], "-", RowBox[{ RowBox[{"(", RowBox[{"2", "/", "3"}], ")"}], "*", RowBox[{"Sqrt", "[", RowBox[{"1", "+", RowBox[{ RowBox[{ RowBox[{"y", "[", "x", "]"}], "^", "2"}], "/", RowBox[{"x", "^", "2"}]}]}], "]"}]}]}]}], ",", RowBox[{ RowBox[{"y", "[", "5", "]"}], "\[Equal]", "0"}]}], "}"}], ",", RowBox[{"y", "[", "x", "]"}], ",", RowBox[{"{", RowBox[{"x", ",", "5", ",", ".00001"}], "}"}], ",", RowBox[{"Method", "\[Rule]", "BackwardEuler"}], ",", RowBox[{"\"\\"", "\[Rule]", ".01"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"solb3", "=", RowBox[{"NDSolve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{"y", "'"}], "[", "x", "]"}], "\[Equal]", RowBox[{ RowBox[{ RowBox[{"y", "[", "x", "]"}], "/", "x"}], "-", RowBox[{ RowBox[{"(", RowBox[{"1", "/", "2"}], ")"}], "*", RowBox[{"Sqrt", "[", RowBox[{"1", "+", RowBox[{ RowBox[{ RowBox[{"y", "[", "x", "]"}], "^", "2"}], "/", RowBox[{"x", "^", "2"}]}]}], "]"}]}]}]}], ",", RowBox[{ RowBox[{"y", "[", "5", "]"}], "\[Equal]", "0"}]}], "}"}], ",", RowBox[{"y", "[", "x", "]"}], ",", RowBox[{"{", RowBox[{"x", ",", "5", ",", ".00001"}], "}"}], ",", RowBox[{"Method", "\[Rule]", "BackwardEuler"}], ",", RowBox[{"\"\\"", "\[Rule]", ".01"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"pathb1", "[", "x_", "]"}], ":=", RowBox[{ RowBox[{"{", RowBox[{"x", ",", RowBox[{"y", "[", "x", "]"}]}], "}"}], "/.", "solb1"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"pathb2", "[", "x_", "]"}], ":=", RowBox[{ RowBox[{"{", RowBox[{"x", ",", RowBox[{"y", "[", "x", "]"}]}], "}"}], "/.", "solb2"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"pathb3", "[", "x_", "]"}], ":=", RowBox[{ RowBox[{"{", RowBox[{"x", ",", RowBox[{"y", "[", "x", "]"}]}], "}"}], "/.", "solb3"}]}], ";"}]}], "Input", CellChangeTimes->{{3.741966501723892*^9, 3.741966528800671*^9}, { 3.741966571060986*^9, 3.741966572848772*^9}, 3.7419668953237677`*^9},ExpressionUUID->"292841af-fbb0-4901-b982-\ bd43d688dd7e"], Cell[BoxData[ TemplateBox[{ "NDSolve","ndcf", "\"Repeated convergence test failure at \\!\\(\\*RowBox[{\\\"x\\\"}]\\) == \ \\!\\(\\*RowBox[{\\\"0.010004999999999948`\\\"}]\\); unable to continue.\"",2, 5,1,23909664869109542353,"Local"}, "MessageTemplate"]], "Message", "MSG", CellChangeTimes->{ 3.7462030749367666`*^9},ExpressionUUID->"7d56a2cd-5bba-40e7-bb96-\ b3d811431c1e"], Cell[BoxData[ TemplateBox[{ "NDSolve","maxit", "\"Maximum number of iterations \\!\\(\\*RowBox[{\\\"32\\\"}]\\) reached \ at \\!\\(\\*RowBox[{\\\"x\\\"}]\\) = \ \\!\\(\\*RowBox[{\\\"0.010004999999999948`\\\"}]\\).\"",2,6,2, 23909664869109542353,"Local"}, "MessageTemplate"]], "Message", "MSG", CellChangeTimes->{ 3.746203074953768*^9},ExpressionUUID->"b120db42-8032-4754-9f6f-\ 0ccbad6ab90f"], Cell[BoxData[ TemplateBox[{ "NDSolve","ndcf", "\"Repeated convergence test failure at \\!\\(\\*RowBox[{\\\"x\\\"}]\\) == \ \\!\\(\\*RowBox[{\\\"0.010004999999999948`\\\"}]\\); unable to continue.\"",2, 6,3,23909664869109542353,"Local"}, "MessageTemplate"]], "Message", "MSG", CellChangeTimes->{ 3.7462030749697685`*^9},ExpressionUUID->"be4d9328-c3a6-4702-b76c-\ aca14fc589e0"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ParametricPlot", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"pathb1", "[", "x", "]"}], ",", RowBox[{"pathb2", "[", "x", "]"}], ",", RowBox[{"pathb3", "[", "x", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"x", ",", " ", "0", ",", "5"}], "}"}], ",", RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], ",", RowBox[{"AxesLabel", "\[Rule]", RowBox[{"{", RowBox[{"\"\\"", ",", "\"\\""}], "}"}]}], ",", RowBox[{"PlotLegends", "\[Rule]", RowBox[{"{", RowBox[{"\"\\"", ",", "\"\\"", ",", "\"\\""}], "}"}]}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "6"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "6"}], "}"}]}], "}"}]}]}], "]"}]], "Input", CellChangeTimes->{{3.741966546889476*^9, 3.741966550747978*^9}},ExpressionUUID->"6dd662b4-c2ac-4ef0-a31b-\ ed2f6439eb1b"], Cell[BoxData[ TemplateBox[{ "InterpolatingFunction","dmval", "\"Input value \\!\\(\\*RowBox[{\\\"{\\\", \ \\\"0.00010204081632653062`\\\", \\\"}\\\"}]\\) lies outside the range of \ data in the interpolating function. Extrapolation will be used.\"",2,10,4, 23909664869109542353,"Local"}, "MessageTemplate"]], "Message", "MSG", CellChangeTimes->{ 3.746203075108777*^9},ExpressionUUID->"9d487384-ba43-4ded-bb01-\ 711285ad336f"], Cell[BoxData[ TemplateBox[{ "InterpolatingFunction","dmval", "\"Input value \\!\\(\\*RowBox[{\\\"{\\\", \ \\\"0.00010204081632653062`\\\", \\\"}\\\"}]\\) lies outside the range of \ data in the interpolating function. Extrapolation will be used.\"",2,10,5, 23909664869109542353,"Local"}, "MessageTemplate"]], "Message", "MSG", CellChangeTimes->{ 3.746203075127778*^9},ExpressionUUID->"ad5979c9-f74e-4482-8c87-\ 892c8f3dcaed"], Cell[BoxData[ TemplateBox[{GraphicsBox[{{{}, {}, TagBox[{ Directive[ Opacity[1.], AbsoluteThickness[1.6], FaceForm[ Opacity[0.3]], RGBColor[0.368417, 0.506779, 0.709798]], LineBox[CompressedData[" 1:eJxF1nk0lN8fB/AxT6Rsw6P9G7JUtn6+FKn43G8klVCR/FIia/oqvqlEm7Jl iaQka4aKFjL1Q8pYCoOxjKUoCYUxhsGMNX5P53fM7znnOc95nfvc+/m877l/ 3HXOZw66kkkkUjrx/v4uPDLIS5k/gaHtRQ5vBxTtAfF7D//fnlAjNl//288i gi2G23xhDRZySCAc94fUi2aPf5vRr6fR634VKviSsQLh/0EgEM9q/T1/XJl3 gzMTDDL8So9x4fxwcI4IVR8jHHSgavSnSSRs3tEmNypc7zZE0W6wRwiL0+tz R3bFwvl0+/PDwvXvwMxLtcAhwtZ22/uHPO+CahStlC2sFw86E61FfcJ6D4C2 29esm3BnsJxdf14iPCSxXb4J6yfBcFnCmi+E16+9Ev2jJRm86w7qfBL2kwqP 966OYhHWbx935gWmgVVqiVy9sL90sPjLpJxBuD0ze25Y+hHwo/xdPgr7fQQp 3frFpYSv+JxI5KZlQHKSVUyxsH8q7F6VkPVG2H8muAzJeWYT/ihe2zhYkQVu FKPLVMIJ3CzFwajHsGa1sVuKMN9TsGhczowhLJlmsGSg/ykw/lQUhAnzZgM1 QAldF+Z9Bp2MoXW+wnwvwXJpjrUVYUv3xcU+K3PBQ/TYGjNh3lyYMA2e2CHM mwf1fhEz6sJ8r0Ci0mR4ToCh+DnfMe8d+eBmw1YeFyzkzYcUG5LcgGAhLw0O SjINGwULeV9D4jIzvVTBQp5C0BK/e+JfhKPDfM99jSwEepl+rALhrsNh2o1T hVBwQ81KinDIWH5qAasIsoJ3nhvkY6hJS+JmaGgxPNp8RzyD8KnUQgu14RLY 682/IEq4801fkOQmOvQYrmeNjBN5UnwdQ1zokKHpktlB+OjpsJWXGunwOYWs 8pKwyBLaLaecUpjfaH7dlvD+nRJndBzLgWHn5Bs7hqHz6vf2ZceXQ8+1FWUX CKdR1m1UrS0HjwQfpeOERzv1v68wrAC0N7VLnXBCgPOhebkPINv9bm3xKIZ6 aYUGzI8fIUW3/VwzD0OBGzxIpzcxoL66437NMIYKblXl7XFiQPVwvU0OYf7Q xpMb7jJgSt0/6xZhb9rAh+4pBhg/lDXdTfjEX14R9h9qoOfRJVYJF0OmR72X 73aoA/VdqaoZQ8R+7d62XHZXA9h8UqIYDmJI1Z/NvXq8Ae4H331NIVyWnVjJ vdAAxvFes31sDM1LzVyszW4ATPtaxF3CF5vfdoRSGoFC14viDBD762SUPtfR CBWH3hnf7ifOTwDS4viyIN/gsW/GDwz9/Oq6vyOCBR/W83j/EL6CIrwZVBbk BaaomBJ+Idqa+6SVBbFBZUXdvRiSifHa4rqtGbhU9jIFwg3Ue8adIi1wSq+m LawbQweZHOuGmFawyrXZsekbhtg6cr4lT1vhtHkSS9BJnO84g7gXZa1wbKpj fwnhV0eCWiLHWyF9zY/FVoTxnhX/3nOkDdhMk2unvmKoeWLnyTLFTwDR2mb3 OzBku+6BH+3FZ/Cqsjajt2GoKtH8Ib36Mwynh5wOILxNfpJe2/sZ7hjIUvUJ Ky22k/yxuh3OMLmqz1qJfjjyGctC28GM7dER20LkLYhpOO/YAecHT507wMLQ E6swTUPKV/ilXFRwn4mhufdvRzq3dsFU3nT63goMFbdXTV7e1QXhjoeOTJVj 6JKghbT2YBfwliXveUxYoD1COerVBVcuqU6TCXOTVP9sS+6COA07TkEphr5d ivKpJ38HuvbVpOUlGKLrO46W1H4H61nWFWohsT8vyONpJ3qgNzlt4/xzDB32 jwzDan+CuYOypm0shlxPFmue2cAGSZMDGdSzGMqUNv1j5Z4hGDMx68MtMfRw 7+J6HbsRMFHCOBtVifPqnaA/zuEBl7Hc4OskGU2nB6oVXxoDKqnA8SaDjJRM yMdsNPiw0s9u+nMcGXVW6G2XeCeAaAs19/cOZGRKE3O+4zMJW0qz34X/QUbm +VoUBYVp0J3SfDf5SQT13plzqXo+A0kKc39w4kXQxi4fsonjL2h7Xs6tNhdB i6Rjzb6MzUHv+mu/smZJqKL5yYCyLQnd0/SazMgkIdvsV2P2WiJIK1ucl7WP hHTPyh1Q4ougkBKxhB1989DDKWtiMsnoQvRK0j+K89BlX/1YMhFDTzfPv45x mINnlv50C6dFaLV4ovqL+F9gMe3/65yhKGri5TervJ6FrfcaorQwMTTEBYcT HTMw2X9hZ+1XMaTQMvbS9ec0hJZXbfOmLUba55k6tZwpcI4+E5IaKI5IvG+q vTOTIEtxFxTYLkGubcc31A1NgFQQrdBJaSka1TKKOzYggLaAHLpk31Lk8d6u aVc3H/y0xgZT30sg2UF5v+XF43CiLuZUWagkarjSctYrfQwU1a5K1TlJoVS3 4WWybqOg8Xc55q4tjULUec2GBjxYZJRnhPGl0YcAl9mLYiPwqnz4YWidDAqJ KvU0cBuCw2vLsur2UVBrUsEboy9sMKLRYvcVUpDSp/+M9C7pBwmZwBSGqiyq fB235KbuD5h34trIhsuislxpgzcB3XCMEr62li+LnKol+gx4ndD9hnnb8qgc SnirVrNG5zNYrhtQaXovhzLl47CC0yxQ7dNkVGvgSOrGrDVZjgFKmatWh9/G 0ce5jIbIyrcQaZzH3xuLo8m+4+YJ9m9hos28UTIOR7nzHtbPBouAufRiWMw9 HKk0Snt+kCmCwLOtgvvJOKpXal+x1a4APm2PY2Xl4Oj7omLbW700iGmSjKqo xBHrnliC1HAOzHpRPUKqcbS65QHeRckBd9EdpuY1OPK6XnMxXDcbjLd6zdQw cdS6RU/Dwe8JcJKrPVktOMpPttSrmaaCuWeIWXcPjhrkHQQR5CSgkRWUqT9w xO1Qz3hukgiKSa9/ufbhqNpnlXPhzQQQMHtpA2wcmckm+3FE44G62USFx8OR rr2IxZRkJFCY7XOvxnAUNNl9+2lkGAS4+7af4+PIsi5r32WJYDiYmB47OYmj y5QN7FJKILzXM/y7aBpH9Pm0AQlpP1CvazAPnMXRqaYM0chlp+Gum4eq8RyO 3oVF2+ivcoD/3f7k0cI98L+N6uFF "]]}, Annotation[#, "Charting`Private`Tag$1322#1"]& ], TagBox[{ Directive[ Opacity[1.], AbsoluteThickness[1.6], FaceForm[ Opacity[0.3]], RGBColor[0.880722, 0.611041, 0.142051]], LineBox[CompressedData[" 1:eJw12Hdcje//B/ASStqHFKUkSpEVTV5HhShZISUNIxktScgoaQkNiVAkWYkm EaVCS7uMktLep3Puu3NOy+/6PH59zz/n8Xyce7yu9b6v+8xzdN1xaJKAgMAp QQGB/77/95F6K604PloOg3f73ncp7cXT3z+0uxT+Z2fUCvevudRRhqSrV8wH vnvgg/2blFsb/vf7GegFf5ssIFOB4s6VGq1OFxHw6OKfr+kVE8f7wXLG48eh lypBqQxe7h25gmQWi6F7tGri/GDEb5HtSDtcDb/thex241AYMAvru47WTFzv BgJ/ZyXnH6uFSG75a9b6cEToxypo7q+buH4EnI+r/52m8x3b9hh09jnfhFzQ 9OxNPd8n7hcF3bXlT3Kv/pi43x1Qn+IXJP78icYrMns6U2LwO0dlhdH5XxP3 v4c45YTFknL1WKh44Xpb7X1YVjQoKz6un8gTh1J/wwx1lQas/kU5Dvo8gOnr aWbNcQ0T+R4iQz4x+p3S74l88bD5KZZlNLsRF9ztY/ofPIL0pjfQvNE4kTcB qYsGm54MN07kfYyUQypdz4v+TOR9BvPGH0Xir5og9kBnWlfnM2zZWdrtVt80 kf85fMRmVoiINE/kT4JPZny3mmPzRN5XsIlOzFoq8nciTyqOBkrcN1zVgqhx D46LYRquRu+b3WfVMpEvDbXfd4Qp+rRM5EtHWbPAgFlOy0S+DLjkb3EoYrZO 5MuCt0lJiNfyNlwP8vD8HZqF5U5P5Rq3tqFpd9CSSn4WBKYH/Xp2og0BnLS4 t9XvELdLNor3rA1Vi6f7BwZmw9OiV9FNsR1H47LMFwzkwHbywTvC3HY0Znb4 iWnlwm9GZf926Q6IxHrYBRzMhYSe60tNzQ7YHA+SO1uZC1HFD5NlbTsgOC09 xOHFJyRlLZLbnduBLUbTXZfZ5cP9893IzLOd8Fp0y+x5VD6WG0XFaId34oHU PHXV0nx4mCUufv+kE+zG1c2z9ArwtfwapV3TidvnHHf+k/kMlyl/MEWjC63p WTplX74gPTyzb2Z5F3zUjggc1yrG4rhF1wfHuvE2pDBlk0Mxir/WRVdJ9IDu Uz+gdrMYls/MrK4q9cAlvevzX34xspJOM/SZPbBfd+zq3s8lONJ8w2+Lbw9M bFxkN+77Bs+xzfoPx3sQsFFfVnp9Bfa1KzVVdvVC9Ux3/8X9FViwZeOZSbxe 5D2P+dp/ugI980by5Kf24Z/4iHfp8wq4BLS9yJnXB++a9/WBUpXwq5p/Psuq D0cd1jwcr6+E/gILZ4v8PlicYy7u9ahGqu43K3ZoP9p/H9pSf7Ua9zaasftu 9+MC86pLcUI1Fm8Oz/2Z0I/kKXWvn9ZVo2onI3Nndj8kw46tOqRfA9/Prt9a u/tRkXBrbaNgLZh5crzU9QPYUda7rSKsDsu/yu7wZQ+ge5mMR86zOpi3MGYL jA3AL1InMjmvDilLBmU/CLGQauVXG0rVgXbvfK0owQKjZZb1JqvvmNLda3hV hYUartGBPKUfEBGqWRmzmYVd8+6cSk/+ichbQu5zolkojDG9m1v0E/6Tj8rN v8eC/gxebmnrT+wrK+nQfMiCsvAesbbZv2BTujZ+yQsWuntnPJoZ+AuH518d qf3IwoW3YRVedvVQ61zPcmhj4enWIE09qd/YnHNmiNIaxPjH96xG3SYk+izx 7MsYRPavQt759U0QTAxweZQ1iLNDtQKKO5qQN+N6zdYPgxhawpKyOdaEZrWT 9OOCQfTfU13+/X4TPiZ8+c2sHsSfs9fcyyc1w/SiimLJwCByV9uxc0qbsTTO yLF8IRt+yZOoB/YtsM70z3kdysbLWssLFi4teL9ZrYV/g40fo4kio+dacL3q +1HDSDaWmpkp7oluwW+x2+Ypd9ho6IjcIF5Gjs9xPhmZyIbOvAV3zhi0QkBM XyUrh43eyI1rd8i2YZdqvP4cFhu7z4QGCZW2Iz7qtfTaTRyM68rruP9sh+i2 7YW1Zhwk8h63Nba3I3abcNhxCw64p3OM3gl2YHLD5PTrOzm47cUeddPpwPrU trVpthz8Omnl3hjfAc4Uvvxddw72u6haZ53pxAAvx6ssmoNDB7I1XdW60abh 8XpdIwcyfhfK6ZXd8JwSEby8iYOcB8yTPsxuTL19YbfCXw7kGwveBe/txpuV DDdOGwffrL5tSrjajZQayz83+zhYZdHo9HOgG2VZjFsXRzmYoieQYPy2B6O3 KqfoyVNI3ZO3sbigBzd80m/yZlOw8/Lv2VbZgwdrGAvTFShkpYms3N/dg+DI +2dVlSmc0JLO81bohWtEf9EkNQq181WaXl7qBaNg9doj2hQeS5goyG3qQ1zJ JSt1Cwq5umZGL3f3wV32eW/9Vgq/HHccMTrYBxUDBfHQ7RQk3tinH7/YB5mg h/rNlhS87c6bf8rog+yBXXZBNhQ2v870OarSj6uljSe9j1Do377od/ZIP946 JmzT8aUg4rNMaMe0AQhEO1qk+VFQSdRZ1CE7AH4cz3CJP4U9w+tPSa8YwHGt SSdVA0meeEfxw0cGwJrb+4QfSiGCc3eNZO0AjIdOfXa4TeHuZuHyZXtY0Pr1 tVI8mYLGjccNgbYsbHi3rmv9Kwrvqo27Gw+wwP/03v7Ca9KefZemXHNj4cqn 4q89qRTkXPkGXcEs+Eu6UklvKNyM7Hn6MJuF7f7uweO5FEIbyv1kVAbBX/Cr 8mUVBYV5LmHO6oPwiGj4/aGawotDYrG5ZF19CpeLKauhUNxvmuViMIi8L9S7 vjrSXsG8gWLLQUQxIiiZegr+C9L3XQ4cRIjZw6NLWyj4uNxeTfUOYuYqr1A/ Num/fzt3a3IGofdmL9uJQ/KFSXo58gfRnsQwt6AopKcEZFRMYQPx+7vkhyj0 Up7aL+eywR190HmHT2Hfue0rDm9jw6tA3dpVgIZBkKjWjzQ2bl+8/dFcnIaI /Jctku/YiJ19kCshQaPume+JDbls2NeX2VUSe5TwktJL2fB1vlO7R4rGc4kO zfA2su6E9S/uYNCYHZW/aPMsDrwOCoXKy9MYfuizIPssB1Xlb+M759Ok/lgH nPUl83jZC7UYVRpC5bodukFkndRssDFdQENslHqaEcXBvaI6TuJCGkq7Tyx+ 9ZoDh5A1l3YtomEiun9lfDsHy0P9F5zTolGobHjTnqyL2EJ7K4WlNLbozKbn UhzsUPk7K5t498G6zLuCZB4383rZy2gc+WihH0XWAUuwzc16JY1rHlgXtINC iVL3ml4dGjLBivEb9lJIqbgh4K1L41bcyKQp9hSUjTOLJuvReFD6psD3BIXZ SbOnK+rTSF24zPRcEAWjsa8zdQxJf/1S3urykayjhCWKYNLkvEm2lho0Ep37 /T030lCbdk43juQINo3t+EGsVc5hdK+m4SS3IdPQlIahTVvxRWMaq2YFaPwj tjpZqPtiHw2BLrkTXptphD26PmPSDRrJsiUDhltoRDsLs8yjaNzcdVwinDh2 6aWS6Ls05LtW7/9LnPTe3W/JUxo5VVrMAAvSX9WWLKtPZJxU8fPzVhr/hOaU vuLQ4D8biGDsoDG1OPLJMJ/Gn/1LpK2JxcLELq8XGIKM8RHPB8SzFcb16sWG IOyrGLtoJ43V2s1Ppi4cgltu9VpNSxprhvde3r54CFImx08eIzbJrdp/b8UQ wkwW7npOvMO8YOYKDGHFh/OS83bRcDn45LKt1RCyS6SyRHfTOKWhbPd0/xA+ zdFuMSL2Yd3W5xwcgruLbqA3cYhPyGCQ+xD2vG1AM3HizRN26SFDCHAdrH21 h0ZjwUqD6R+G4PS5YFB+L417UuV7fucP4ZrK/f4NxHv3HfV8VTwE6agrBp7E NZwHL3f+GEJB3pOsEuISFQnle5whuFQ7NDhZ0whyeWboMjyEc8maDteJN7wz 2csU5ILPCtROI87bfi6iVYKLs2f56jTxpfuyrzJncqElkP1vpg3pj66UkiAF Lpk3suo6xG8vdU5eokGODzJ18yL2Kr08798yLtYqnQmPIF4pp7S2UoeLbQfD DiURJ7/adfrUei6SKBP7v8THR1iRG825OElH148SL9oY+lp+Jxefb2kkzyTz pD1C7VvPXi5ytd1+aRInNOZ1fbDnYoZ2WhaIHTX2Tw1z4kLz5dRblsRKXnwV RxcuVp5y/+xM3PDpJrRPcZFdv8D8AnGM+LJ9U324sFy2Wyuc2GpvifcPPy5C 8rTdHhLLPj4c9TyYi0PuWodeEYcb3i+ziOYiUFQ8+huxRZBuj3IsFwm+j1Y2 EIvVVAtzErjQMWhe2kVcrOSq+vkFFxWM1+Ec4qBjouuiU7nY+PRj4gjxhjeP bZ2zuLCqq02dbEvqhdC6swa5XMg+P8uQIM61aLgl/pW0f+fGClniCzGn0/58 4+KpUB1LkdiwXaYipYYLoe5cHxVi/vLk3sv1XNj3ymarEmee3zRt918ujkWf HFAj9ixqXaDeRcaDmumtSbxi5iWj4QEu9M+ss9MiHrCfY1c6xMWne3rJS4mT kjLPxY6R/qtd5/Lf70d522+7TeaBo97RrEGsbtKXbjSdh1y3iA3qxG03gipn yPBw8aSe4ELi+Pr5/e1yPASI5inMJ7ZXyxHNUuJhkaRMwlziuSet1a4u5CHb WjJ+FnHDR9rYdgkPZ/F353/t37N78XlBAx7o5ybdY6S/ZsR/vVO9jgdxuX/G NHFVn2PmY1Me0vYcF+8hDtMbrzq9lYcFveI7G/8bnyt3Bjbt5qFoKN+8/L/x qdQWU7DlwX/Yd/kn4iKFCvX+Azw00w9PpBEHHDm2PvcoD62++pKPiU3SpzpG uPMQpuOoGkXMNM8JM/HmwVUqWO4ysWHr6ZyhCzy0U2cKXIi1Z3Qp2ITywM6w XmlKvCzpoZl4JMk3ln9Hm3ixifXZnDs8rN/LjlEiVm+Qeeb+gIeu+ed3ihCr epZ8n/+Eh4z1t4L6yfxXFvOfWveSR95fw47UEiskGK4KSudhWhOdn008o+Zl RO8nHqKFNQ4FE0sdP/wptpCHvzYSHceJxSYrsbaV87BsWe+HLcRTtMO2ZDTw 8DtYpUqKeFKpqY9TCw/dy2fZsMn6Hj8g+EK+m4zXvY43VcRDkR4iF7g8qM+s vx1K3M6xzN8gwwd7v+7XEVJP/l4VZ/Pk+HDrPNldQ9w4/4vyCyU+smMo8xfE dTt1Lkgu4eP+qsU6VsRf0uT1fmzkI2StoPhzK1JfzKqcQiz4oMM+hp0l/tgS cstwFx/LradTZsSZjFHOA0c+guaWd/SR+pd4sjHZ+TwfZ9cGDysRx0+P/j3H n48U38nTe0m9jH20VawshI9d+ca2b4hvVec6r7jNR5dWp7UpceDKR6ojqXyo xqTscSD12b/EZufLLD5mKVwbWkB86cAMP7tcPiwi7dd2knruHXnlT/43Pmyq RZ86Ex/hOMWEdvJRFl325wh5PpimLZaaO3cYIt5FOau3k/YNX8zapzoMptvF O+xtNNYbVTve1RjG+5ubVyUTG1WeyZDTGcZka52xecSGA1+sGduGYRCjEzNK nl9LNR0ei/gNI7suPSLSnNRHj3SLjUHDaFZpXmVMvPidMO/KdWLH1JcsMzJ/ NiVvFro7jHe+1VKbiVWcRgbG0oaxz/XgOH8TqV+PovQ5bcNoU+V5LCXP45vd nS3Le4cRHfOEqtpAg7HC8Jobexg5qsa5XsRSn/7+6RsfxogJFZG9nobon6VX OmeNIMTiyHQ9Expjs4vLf28awf5HySqy62i0RowfLHw5gss/ZhiYkP2CvL3n 1uvpIzhY2NXwmewnLJZ061m+H4HtIxMbU+K3X2slmgpH8LlgRZI52X+Ejia9 5baMYE6x5DsTsk9YddhGTG32KHbpS06fsYKsR723aVcCRqE21yVVkuyP3k/V ijW7Noq7C47eO6NOg1X9KEj65igWGchz/6rRsHa5YXv/4Sh4qmblr8n+SuvR YeGM7FH8865cYkz2Y9/FZtq0skchsynt0jol0n9N7pOM7cYQE5h2/vEMGk37 /JyiDo2htdHzviBx9M+I0o5jY7jTWam4j+z/plan3Qr1HoPlI7NeMRnS/s+0 Rl34GHTMaGkXSRpxL87sdC4Yw+7rtz1FREn/n/ZNuLFoHOIXjzvkjlOYLBG+ oYEzDmr05O5l7RQ8pAp6M4bHMWfut9xtbRSaZLgRNwT/wa5C/pdrK4XsWbZ/ jCT/IWBps/HTvxQ85y3yfqbxD4XPBrRm/aHQpv3pmZfDP5xblWdR9J3CF5tB MenyfxjQ/Dj2qJBCQc3TLpVdAsyrbqJvZJ9S2PU8lbN3sSDzmbSVbb4dhRVu MtuVaUGm6PVe+6syFFp686rKyiYxV6wpM4xN56Bpb9ETsRghZuXY1I6X5P0w yeJMrrnDZKbFrOpXxr1smA+fGfPUm8Lcvv1sCteDDd1bFdcWC01l8pUrGBfI vp3Xedqo9PdUpujexhCpE4MIzC/Ud0kXZsrpHfEP+cWC43XXgDgfEaaIccZZ m80sSEs5Db3dNY35KKU1dNhsAOJ+6VkOyqLMFa5lBY6i/fh+7kWuWIcoU1TC KvnBp16cWszpifs4nem7R7RrjkUP7L+FHc0LFGMyD0q01bZ1QWnBRfFvDuJM +RzdtanbOqFxIl/IaYkEk3nnuQRV2I7Ja1LWCNESzN2DdmN+hm1IzR+4G/hN krnlkDDzx6UW7FbMS/xmJsWkJByFNOqasSY9PdwsS4oZeHPjvBDqD6ZL+sQW q0ozG9d9DtZc2YB/Dv2W0sHSTNlJcQcfxvyArVSwYiktzZypE8wUuFuDv5ll NyxsZJixG933Zn6tgMW8rvlVH2WYU4S3xN9uKYJqh2ZxkQaDGdLC3RZi/wnK j+VnB99gMPn+jQ5F/qkIXZtCbw5nMC1DYxwuGKSC+920UiySwQyLMlLXY6eg TNQ7KOwWg2m3sPNegl0KfNzqhqLvM5gC8ydvHdZ7jR8GkdWJLxjMVt2gDs36 FwirErtW8JXB1FmYG2ceF4/RYwlHAooYzNjuNc1Lp8fDaYqhiWkJg/k01DCS 6/UQa3WPjZSUMZi+RqemLrd4gN77Rc7VtQymdNDbgXHeXZg6B2z428JgOgZN ursnKQLpk+aqJLQxmCe6o4Zc/MKhdC9j7FAHg/k2T+3ZaaswDJW1pnd1M5hP pK5nNAmHIkHbeP7gIIOZ59pWTV6QIVX2azyVQ9p7RZNbYHsJ55w8fnnSDGan 1bmkrZI+2BHzMJzHI+1ZnORQm+GBjyv1TrwbZjBvv28dfStzHIu+VZj6jDKY 1/NPrxHOcsTNw0dU146TfAKblc+834H///d8BvN//6P/H7Hyzuc= "]]}, Annotation[#, "Charting`Private`Tag$1322#2"]& ], TagBox[{ Directive[ Opacity[1.], AbsoluteThickness[1.6], FaceForm[ Opacity[0.3]], RGBColor[0.560181, 0.691569, 0.194885]], LineBox[CompressedData[" 1:eJw12Hk4lO/3B/BBlpAw2lOUFlKpFKW8n6RIJZUlUSlJyJoWoVSEUlEoqSja 02ILKQqF7GvZyz5mnuzMjOV3f66f7/wz1+uameece577PufMKB513msjyGAw HAQYjP+e//cQWCRjvmVaArTSLD91zjfHhSd2iVm/4idsh3N0J8NJKQFvrvvt /FftBoW4hW4hbf97vwfuXH6h5T0lCfkda1RabC+iX1Y8xMc/eeL9l2EvWPOU pZOC/gU9V9h8P0wfeBknviJt4vOBmHdLom6lRjou78ntbdsShFj34GMvd36Z uN4ttCvZbLJekQmxzOL33VtDcMh8kXKXyteJ69/Gz5d7psoxv8HITKuDYxcK dMQHCjd8m4gXhvphqXUPI7Im4kVAyad38YW5OWjwkzXr+HAf/95edGfX5EzE f4D8mtRk+5DvE/Gj8I063rap+sdEvCeY3GNyocAgHxdcre7T0TFg34x/WNOQ PxE/FuumeSf/sPs5Ef8pirVTpdacK5iI/xJx6WN3dK4XQTJaY3Jnx0soZSVa uwoVT+TzCubL7S0VvIon8nmD1V833v+sUTKRzzuMKH4/Mz6ndCKfeKzcO0t7 e1YZwsbc+pw2JsDao91+2mDZRH4JyJsdnma0tHwiv0RceCavohxYPpFfEoIK axxst1ZM5JcKxkHJoIQXlbgZ4OZeH5QKnVrmQu2qSjSZBiwv5aZCPvSKloxQ Fa72JUSllKdBadWwWtjBKpSpSvj6+6eDvVAq1VqyGvZRqTsX/ctA1Ho933zz X9ilI+GsdjgLJmqfEzen1uCMcviOV2FZcA/IXSZYX4NoacWlSgVZEPytXFDI qEVvw7o/M9ZnI5vf7R+rV4t7nkf3jcvmYEZX8Fqt0lq0JKZqFH3/Dg8bT4Gv tXXwWnKCcXJFPoqUFySioAEp13I/bD+SD5szah+VOhowwFlqvSQ0H95Z7EVl Qo1wSuzM+cvNx6/QqQcVNjTCarPDdfOcn/hbMktr/7NG6Fo4TdezLETt4gBb sY1NuKq3YbrM1hKwJGxP1Qv+gaEnpcp2K0frkn//Yh3+oq3eZlft9XJMLeuo ljnzFxeo6075seXwULcRPOXzF2+Fq96/qCqH5fEvqf9C/2JqsMNamw0VUM9S EIn+/BclseHaDQKVWEJTG7dMacbeIrZRSXAVFi9jNXU8bYaJYsTpxLe/0cu0 uH4oswW59/UjM/N+41aYcppbXgs2yA1nFrT8hujO7R8DylqgIGom2Tq7Bg/2 sBoiW1rAYsvFTPOvwef9a8bWiLXiQkpwyZnDtfj0SkVIzbAVL3YHLFsvXY+V K7pLHlS0YuzLp+4GzSbc0dhrJ1PUhvSa3GHvrU1I1h6Zza9ow/nBSob83ibI iRntrqxtw+DybmkLhybE/DJ19etsA/1AaVX1wybopUvdyZvUjsbzN1yLyfeW nKF3U2tDOzLXHe7NKPgDrf5+3azH7bj8VrA/2qoZEs84OqesOxBXaXzB0Ims u+z4ULpdB36NPBMb8WxGxEmNj3znDqzcsUPe7G4zrBoqThzy6kBd+51tU4qa MVnBZNmMsA5oKC6K8NBqwXkXq6WXcjrAvqOnvXd6K2za3iUvWtgJU4+gAKGC Nhi4zdF/WtqJMc1ZGq6/2/DB8UL4t6pOPBt+2trQ1gZ5Sak9DbWdGDqboZMm 0A4lzZZzgm2duHemd8RFox17bOVUpvE6UXNqv2vDk3ZsSHIYD1Nk4ZCT0oFU jw5Ibn99us6RBRvr9GXOS1jYkq/2R2+EBdnLF4oH1rAwdQkjVovRhYxo6pQX xUJZlV7W4kldmNWQnRZozkJCeWXouEQXCvcXbo+9zsJ2GZhHzunCWsMG29// WLB8bFOltaELwusZsVtSutBj4O0bdroLT6V0587czoHosefSfi1dyNTcoRNn ykGCZ0b8vI4u1Bzde0LnGAcm/ZKH47u6IPXRKvHkRQ5miSXc+97bhXOHvXd+ TeJAPD3Tp5PBhsH7ZC/7BTQq+moCg+ayQe9Rrk/n03jwOGB9wx42Ig1Ei9XM uvFgTaCtayIbKree1vkf7Ebk81ymykc20sq3sBqsu9GdXn28OZWNGksf4Rsu 3fA8Zs7ancHGTGeuVmdgNz7qc0+M5LEReqfrxeP0bjSKPIlhNLIRVFd8WXZB D0rGWYNWohx4Od1b18/uwaw4/4jhPRwYjO8zXdbXg2n2dcbzjDmYGTz1zFFu DxJ/W6ZuJetO/HA1qUS4FzVG5g9vH+CA3e+uHjevF3NzHZuFrDmw9Nyz+rhR L7I+0FoPTnGgFSC+4lcCeb3lvvvxOxyIzfq+a2paL44k2pUbhnFQ9fKS47bM XpxKdvTTuMuB28/hN4kFvSgNoHIlIjl4JdW+LKS1F6Yv7Z1Dn3AwOyxL2WBG H853xg8af+CA99hrUfr5Pgg9bLkSUMgh5/HA1fOX+uBl5G1nVcyBULFmu2ZA H4aexV9ZXcqB5Ej/i6SwPkyv+BBUU8HBfFNH1Xfv+3DGNllcpo4DXfFDa560 9cHY+Z2KPIuDG27YHLC3H6X3dTakTKIhGyj/ZJt5PyS+HFI2FKERHsUXFLbq h4vW0uZWURrRBR+zLzn2Y3C/hDhTgkb8YjV9z4B+cMKZ07bJ0KiqUdjt9KUf mxqC1Q7L01DYInjQWGUAT5m+xkx1Gksme2pGqQ1Ap/FrjPlaGiuK+5isdQPw embc/WgdjY0WrfkXtwzAVp7pp7Cexv5TuZqvLQeg7qxbKaFNIzjmppzgrQGc cL8c9ECPxrjQnIJ3fQPwmbJw90VzGiL5d57zuAMwO3TMOuQADclgyStbGYPQ SD6tH2NBY/bcsfW1koNYWi2XlXmQxjr1P89FFg9C4NWJoIYjNJyOPb9ycP8g 9IbmRNXa0WjIXqMl8XkQ77O4pU/OkX0sXWxWnzWIAV6VcJgHDXNLe/d3+YMo 7vLJ8T//376Pjtv3axCsA+eqbb1o/FwgpfCgbxBi6itlRH1opPh0TFquMoS9 h/+ljFylcabgiuK42hDsUeJd6k9jzcz52qUaQ2hZZ6D9JIDG23cmZ09vHcLa QlaB3jUasQ3fOj9bDeGeu669xw0aIRsfFhneHcLvRstjN+/QMAzQ7FJ4NATJ iMkhO0LJ+ivKRftih6AyuqFdOIxGgIP45rvxQ/icHHzeMZzGhftnExoLh6Au NEVEM4KG/fCeey6ThpFzIkLy5CMaS3U5iToSw4gqr1shEUWj9VZAqZzsMOr8 5NpfEVstyRBPnT+M2pgbDq3RNMxMVb0FtIZhMJQprRVDQzdR5Oht12EcnDwp //FzGtTOjGDdc+TzNVUCC1+Q+91yNmPwwjBWxakgmlhdrnOuRdAwynw3H3r0 koaS+8/qhc+HUd4lme37muwnSV+RqrhhKJ3kxHKJ58ZuXBuQOAy9rAZL5zc0 5CribrO/DiOtkq1iFkdDWD14V1LdMEr8uPGi72gIFuh72TYPo9/m6tZTxGPW Aq9nsYbhnb+0sI548I6b2IWhYXBu1dAx72m09RlnbZPlwlz1naZKPI3vCbPW /9Ljom5vlFhzIo1vO8psrxly0f9u+4BaEo0vzdfCN5pwMV1wha0PcTJzpC/6 KBfx1i218sk0np1qeGvnzcUN6/up6z/SeCJxt36OLxed5yLaLhI/itktWXSN C6en53KzicPLM+1W3+PieWNhkm4KDf81MUr8eC6MRlZ2qKfS8P1psS8ulYs3 ZtsUnIh9rOUuH87kwnRu99YY4nN3/BqzCrlI3e8zXSyNhvsybakzFVz4/675 qUXskjW4cWktF5EJlrediE/02d4P6uDiRGiZy0/iY9cV8rT/ccGu8eX3EVst /D3UPUCu71tZJP+JhuWnkMUxI1w05Hw11Sfev8/AxESIB+/EufJuxMZdgr6i 4jyo2KS+jiA2uvIpPlWah3q2fdBnYv0EVel583jYyTV9NUr8nXcx1VKJBzOl L0Vz02ls1Sk/GqnCg42zp78WcXbgYskaNR4wqfPZfmKdUo+kmRo8xEkeGXcj zpxZeMhsEw/n4zdZBxLDSkEsfAsPGw2cjzwi3vjv+wGmEQ8xkXn2P4nT1s2e tNeUh/mut+f9IV5/wTEu2JKHq2LmH/uJU3IyTYuP8tBWt2uy8GdSL6bIMaTs eLA6zY+QIU40tn2505mH7BM6ffOIVz9I23v9NA8m+dYrlxO/b54ykufJw4iQ lZIW8cplR56KXebh8UWNF3rEb90SDfUCeFB9zgreQ6yaJjrsd5MH/xaPRDPi VwIWj7NDeWjWC8+2IF66/a2BUCQP0zcbGBwhfhYs0L/5MQ8Rf/+ssCVe9Mv4 oc9zHk5uWWJwkjh2/ottGXE8hJxvOeVMvMCW/280gYfJUc1X/nP0W8OIjWk8 eJ3fFudAPH/wsY5nJg/Beu4rThA/2jTQlfqdh8IVofE2xHP99MOGC8j6ytYw rYnvF0Rqa5TzYFDRpHyYeJbcv/bTv3korlpasp/4noVOSGIjiecV0WhIPD0m bENfKw+HDzXs1CEOZXU0r2LzkBbJUdQgZq7eeMOll4fXu/afVSW+7XFr3bth Hjwe2S1VJJb++reRM8aDXcpJaTnim2LrAlWF+cjMeSM+iXiKUeBqBwk+lAbv ynWT+yfeuNKvYwYfOzpudOQSByy+smLJPD6kxvh/k4hFnaqqbZT4sLmZlvqY +GqS8qVYFT5ajxdsv048adRLpVmNjyUm6i6uxJd1S8oVNfhYZPvyijExI2ih t9UmPlLLf1dpEo/Ozi+u387HF1s5CBJ7HZX3mGvEh/Jff5lWsr95L10WWJjy oZK4SDGHeFBzxulfR/lQ/HDgvRexu4/9vBl2fDg66BmZE/f9+PzDxJkPYa1L MuuIu02PzS735CPC97BHFzmPTo8+Zslc5qMiIafpKzGnVdzRKIBPzq3Jn1Bi lvuHjMJQPjKuz47eRNxye+xYbhwf80MXll0i9WGWlfvum4kkH+45H31iw+Ws 9caf+JD8phgnTZzyo1KqKZePHyWP66JI/QkaeZMy1Ey+L9OT3o9JffqatyDm M4sPq2nv5x8hHgy/d+NKD1n/pfqL84mtVvkenTrOJ7+7hKpvk/q39riF5JLZ I+D/WN3uQOqj/ZrSQbbCCK4url05nziaofcnfskIMh9UehaR+ioeuTpZe+0I CvQDnZSJG4onW5kajWB3M3dLNanHV9enJPhdHcHPLb8s2kj9/iSy4tGOGyM4 q7Gq5wpxd3lMgEzoCJwWX7qpQHzA6dbBh49HMBMtB/e8JfNHzHHRpPQRtDcc q3lC+ke15DSLlt4RaFmttBAi/Uey5trWV9wR3JYrlbxK+pPOc4aaC2MUa/Yt GZhM/GYze9LIlFFoJc9VlyX9zedM1lum8ih6XF+Jij0l57XJVXDL4VGInXQw Syf9scnysm2YzSiZM3uNNYjv/r5d0O4winizz77vSD8VKU8IDzo3isjFk2Pv kf7bkjOgUhUyisLZlnL2D2hEvfbYZ5c9CtUHESK5d8l5Onsp9pbyGBq3mtm8 IfNAUV/I5L8rx+DQbmwnSOzn8sRJfd0YNh+NkzcNotFvl6X5W2cMHQJHavvJ PFFmKVy4wHIMSnnbHFeReeOGTuBA0s0xxLhKHrC7TPa/VMi2ur4xXN1Xfkzw DA036Wx2Em8MVgWbWhacJuuRHbp9S2AcgdFn5292p5E+42CjztRxZPX8jD7l RvazovK5lyrjEN+veKDSicwP6l9fnjkyjg1ZYuG7bEm9t+iRlCkex2hpl/ce M9L/DyklsCrH8bw8jKVrSvqjlZl5dt04Ijw+zllnQuOizeenZ1njaCqQPTVn H5nHPurO/yfIoHI/1H8pNST95qt0isk8BlU16f0BkW2kX1S86FxgwqAunfjn fXYNmf9GSnIu7mdQKkNmLdtXk/63kPu4zoJBlfRMOTJ3FQ15t+0Hwo8yqLZN d/6kryDz2VRW/mQXBqXfWOvSoUzup8GyN93XGJSfoqfLFAVyHjLinL5kMqig jhxVBzIPm7yK7zNXFaA05sW9etXAwUnlpVO8VwpQVfRh7KjnwPfFw8XRqwUo RY0KVkctB/HPAszbNAWoa1VP0xV/cyAdcyjDbasA9a5GpPBWOQcFkeLXrx8S oG73F7euyiXz+42jCz4HC1CuLQdVvr7nYLWL7B6FAQFqg0sku+sCB83sb2VF RYLUeSbH6bswB03mec8l7wtRFYvFgsX82Hhj6JG588gkSjeGp/JxuAs7eR6j 7uuFqRTlmV4b93RBM7zkhqqQCKXrr0svfM7CcMdZnYJ6EWrFw8lrNzV1wj8r d4NToigVPM80Mk6qE0dvOl+N8hKjJPZ69TZod0BG2nYwxWQytef7sjlORu2Y cjkx9YiCOEW9f7vpj00bqj1fZ0q2i1MbrU1kRU+34rRqX1fUFwnKMIBx7bt1 C6wKg+2/+UtSDLu5Wq3mzZi/6OKUwiNTqEdrs60mr/gLFccsIdvlUpRpT7jR YYk/mLTpwyahASmq+6OAK5XXiPisf5H+hVOp+27cbEv1epjKf3tWuEOaatvX crLkYw02JSaG7EiVpl7erJnX010Nialej/KVZKjtK2ctWbCgEuNHaGOZQBlq LfVmka5rGQ5KB8oXDMhQD6xvfz7VWoS/yUW3DC1kqeU21+p5ynkwVOxcWPZF ltrcfffwqF0WlNqX5eepMKkQH9kZx6TToPB01uzAW0zK0hssyR+vEaT9YcAg hEnxX42UyZu/xlC1fqnkHSZ1uFbLWaPrFYrEzwUEhzOp910/3CymvoKXS9Xg 3YdMSqKrtLrJ9AV+ad0pf/aaSRmbJDjvbYlFcJnkjewfTKojRz6yk47EiEPs iat5TKpU2tM4QzoStsIbdfV/MqlmkbYoh9X3oa3pwP9ZxKReZKRO2nL6HtgP 8+zKK5mUDTub8YkXCn27q9v+NjMpq1+y9j6CQUgUnLcgtpVJPX3ixInZcg3z HySN2rQzKZEPRY4ffAMwWNSS2MliUgFVOkvbhf0Qq75lYU8Pk/qeNV11QNIL 0kU1Y/F9TMpX2DLpWdA5eNq61bgPMClKzVLDU+I09t5/HDI8zKSamM6OGdIn 8WXNesc0HpO6tOudsbiUDZQLS/S9RpjUZe1Z7temWSL0+Akl7TEmZR0xembt rJ34/3+n5aj//U/9f3Hnhng= "]]}, Annotation[#, "Charting`Private`Tag$1322#3"]& ]}}, { DisplayFunction -> Identity, Axes -> {True, True}, AxesLabel -> { FormBox["\"x(t) - miles\"", TraditionalForm], FormBox["\"y(t) - miles\"", TraditionalForm]}, AxesOrigin -> {0, 0}, DisplayFunction :> Identity, FrameLabel -> {{None, None}, {None, None}}, FrameTicks -> {{Automatic, Automatic}, {Automatic, Automatic}}, GridLinesStyle -> Directive[ GrayLevel[0.5, 0.4]], Method -> {"ScalingFunctions" -> None}, PlotRange -> {{0, 6}, {0, 6}}, PlotRangeClipping -> True, PlotRangePadding -> {{0, 0}, {0, 0}}, Ticks -> {Automatic, Automatic}}], FormBox[ FormBox[ TemplateBox[{"\"k=1\"", "\"k=2/3\"", "\"k=1/2\""}, "LineLegend", DisplayFunction -> (FormBox[ StyleBox[ StyleBox[ PaneBox[ TagBox[ GridBox[{{ TagBox[ GridBox[{{ GraphicsBox[{{ Directive[ PointSize[0.5], EdgeForm[None], Opacity[1.], AbsoluteThickness[1.6], FaceForm[ Opacity[0.3]], RGBColor[0.368417, 0.506779, 0.709798]], { LineBox[{{0, 10}, {20, 10}}]}}, { Directive[ PointSize[0.5], EdgeForm[None], Opacity[1.], AbsoluteThickness[1.6], FaceForm[ Opacity[0.3]], RGBColor[0.368417, 0.506779, 0.709798]], {}}}, AspectRatio -> Full, ImageSize -> {20, 10}, PlotRangePadding -> None, ImagePadding -> Automatic, BaselinePosition -> (Scaled[0.1] -> Baseline)], #}, { GraphicsBox[{{ Directive[ PointSize[0.5], EdgeForm[None], Opacity[1.], AbsoluteThickness[1.6], FaceForm[ Opacity[0.3]], RGBColor[0.880722, 0.611041, 0.142051]], { LineBox[{{0, 10}, {20, 10}}]}}, { Directive[ PointSize[0.5], EdgeForm[None], Opacity[1.], AbsoluteThickness[1.6], FaceForm[ Opacity[0.3]], RGBColor[0.880722, 0.611041, 0.142051]], {}}}, AspectRatio -> Full, ImageSize -> {20, 10}, PlotRangePadding -> None, ImagePadding -> Automatic, BaselinePosition -> (Scaled[0.1] -> Baseline)], #2}, { GraphicsBox[{{ Directive[ PointSize[0.5], EdgeForm[None], Opacity[1.], AbsoluteThickness[1.6], FaceForm[ Opacity[0.3]], RGBColor[0.560181, 0.691569, 0.194885]], { LineBox[{{0, 10}, {20, 10}}]}}, { Directive[ PointSize[0.5], EdgeForm[None], Opacity[1.], AbsoluteThickness[1.6], FaceForm[ Opacity[0.3]], RGBColor[0.560181, 0.691569, 0.194885]], {}}}, AspectRatio -> Full, ImageSize -> {20, 10}, PlotRangePadding -> None, ImagePadding -> Automatic, BaselinePosition -> (Scaled[0.1] -> Baseline)], #3}}, GridBoxAlignment -> { "Columns" -> {Center, Left}, "Rows" -> {{Baseline}}}, AutoDelete -> False, GridBoxDividers -> { "Columns" -> {{False}}, "Rows" -> {{False}}}, GridBoxItemSize -> {"Columns" -> {{All}}, "Rows" -> {{All}}}, GridBoxSpacings -> { "Columns" -> {{0.5}}, "Rows" -> {{0.8}}}], "Grid"]}}, GridBoxAlignment -> {"Columns" -> {{Left}}, "Rows" -> {{Top}}}, AutoDelete -> False, GridBoxItemSize -> { "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, GridBoxSpacings -> {"Columns" -> {{1}}, "Rows" -> {{0}}}], "Grid"], Alignment -> Left, AppearanceElements -> None, ImageMargins -> {{5, 5}, {5, 5}}, ImageSizeAction -> "ResizeToFit"], LineIndent -> 0, StripOnInput -> False], { FontFamily -> "Arial"}, Background -> Automatic, StripOnInput -> False], TraditionalForm]& ), InterpretationFunction :> (RowBox[{"LineLegend", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Directive", "[", RowBox[{ RowBox[{"EdgeForm", "[", "None", "]"}], ",", RowBox[{"Opacity", "[", "1.`", "]"}], ",", RowBox[{"AbsoluteThickness", "[", "1.6`", "]"}], ",", RowBox[{"FaceForm", "[", RowBox[{"Opacity", "[", "0.3`", "]"}], "]"}], ",", InterpretationBox[ ButtonBox[ TooltipBox[ GraphicsBox[{{ GrayLevel[0], RectangleBox[{0, 0}]}, { GrayLevel[0], RectangleBox[{1, -1}]}, { RGBColor[0.368417, 0.506779, 0.709798], RectangleBox[{0, -1}, {2, 1}]}}, DefaultBaseStyle -> "ColorSwatchGraphics", AspectRatio -> 1, Frame -> True, FrameStyle -> RGBColor[ 0.24561133333333335`, 0.3378526666666667, 0.4731986666666667], FrameTicks -> None, PlotRangePadding -> None, ImageSize -> Dynamic[{ Automatic, 1.35 CurrentValue["FontCapHeight"]/ AbsoluteCurrentValue[Magnification]}]], "RGBColor[0.368417, 0.506779, 0.709798]"], Appearance -> None, BaseStyle -> {}, BaselinePosition -> Baseline, DefaultBaseStyle -> {}, ButtonFunction :> With[{Typeset`box$ = EvaluationBox[]}, If[ Not[ AbsoluteCurrentValue["Deployed"]], SelectionMove[Typeset`box$, All, Expression]; FrontEnd`Private`$ColorSelectorInitialAlpha = 1; FrontEnd`Private`$ColorSelectorInitialColor = RGBColor[0.368417, 0.506779, 0.709798]; FrontEnd`Private`$ColorSelectorUseMakeBoxes = True; MathLink`CallFrontEnd[ FrontEnd`AttachCell[Typeset`box$, FrontEndResource["RGBColorValueSelector"], { 0, {Left, Bottom}}, {Left, Top}, "ClosingActions" -> { "SelectionDeparture", "ParentChanged", "EvaluatorQuit"}]]]], BaseStyle -> Inherited, Evaluator -> Automatic, Method -> "Preemptive"], RGBColor[0.368417, 0.506779, 0.709798], Editable -> False, Selectable -> False]}], "]"}], ",", RowBox[{"Directive", "[", RowBox[{ RowBox[{"EdgeForm", "[", "None", "]"}], ",", RowBox[{"Opacity", "[", "1.`", "]"}], ",", RowBox[{"AbsoluteThickness", "[", "1.6`", "]"}], ",", RowBox[{"FaceForm", "[", RowBox[{"Opacity", "[", "0.3`", "]"}], "]"}], ",", InterpretationBox[ ButtonBox[ TooltipBox[ GraphicsBox[{{ GrayLevel[0], RectangleBox[{0, 0}]}, { GrayLevel[0], RectangleBox[{1, -1}]}, { RGBColor[0.880722, 0.611041, 0.142051], RectangleBox[{0, -1}, {2, 1}]}}, DefaultBaseStyle -> "ColorSwatchGraphics", AspectRatio -> 1, Frame -> True, FrameStyle -> RGBColor[ 0.587148, 0.40736066666666665`, 0.09470066666666668], FrameTicks -> None, PlotRangePadding -> None, ImageSize -> Dynamic[{ Automatic, 1.35 CurrentValue["FontCapHeight"]/ AbsoluteCurrentValue[Magnification]}]], "RGBColor[0.880722, 0.611041, 0.142051]"], Appearance -> None, BaseStyle -> {}, BaselinePosition -> Baseline, DefaultBaseStyle -> {}, ButtonFunction :> With[{Typeset`box$ = EvaluationBox[]}, If[ Not[ AbsoluteCurrentValue["Deployed"]], SelectionMove[Typeset`box$, All, Expression]; FrontEnd`Private`$ColorSelectorInitialAlpha = 1; FrontEnd`Private`$ColorSelectorInitialColor = RGBColor[0.880722, 0.611041, 0.142051]; FrontEnd`Private`$ColorSelectorUseMakeBoxes = True; MathLink`CallFrontEnd[ FrontEnd`AttachCell[Typeset`box$, FrontEndResource["RGBColorValueSelector"], { 0, {Left, Bottom}}, {Left, Top}, "ClosingActions" -> { "SelectionDeparture", "ParentChanged", "EvaluatorQuit"}]]]], BaseStyle -> Inherited, Evaluator -> Automatic, Method -> "Preemptive"], RGBColor[0.880722, 0.611041, 0.142051], Editable -> False, Selectable -> False]}], "]"}], ",", RowBox[{"Directive", "[", RowBox[{ RowBox[{"EdgeForm", "[", "None", "]"}], ",", RowBox[{"Opacity", "[", "1.`", "]"}], ",", RowBox[{"AbsoluteThickness", "[", "1.6`", "]"}], ",", RowBox[{"FaceForm", "[", RowBox[{"Opacity", "[", "0.3`", "]"}], "]"}], ",", InterpretationBox[ ButtonBox[ TooltipBox[ GraphicsBox[{{ GrayLevel[0], RectangleBox[{0, 0}]}, { GrayLevel[0], RectangleBox[{1, -1}]}, { RGBColor[0.560181, 0.691569, 0.194885], RectangleBox[{0, -1}, {2, 1}]}}, DefaultBaseStyle -> "ColorSwatchGraphics", AspectRatio -> 1, Frame -> True, FrameStyle -> RGBColor[ 0.37345400000000006`, 0.461046, 0.12992333333333334`], FrameTicks -> None, PlotRangePadding -> None, ImageSize -> Dynamic[{ Automatic, 1.35 CurrentValue["FontCapHeight"]/ AbsoluteCurrentValue[Magnification]}]], "RGBColor[0.560181, 0.691569, 0.194885]"], Appearance -> None, BaseStyle -> {}, BaselinePosition -> Baseline, DefaultBaseStyle -> {}, ButtonFunction :> With[{Typeset`box$ = EvaluationBox[]}, If[ Not[ AbsoluteCurrentValue["Deployed"]], SelectionMove[Typeset`box$, All, Expression]; FrontEnd`Private`$ColorSelectorInitialAlpha = 1; FrontEnd`Private`$ColorSelectorInitialColor = RGBColor[0.560181, 0.691569, 0.194885]; FrontEnd`Private`$ColorSelectorUseMakeBoxes = True; MathLink`CallFrontEnd[ FrontEnd`AttachCell[Typeset`box$, FrontEndResource["RGBColorValueSelector"], { 0, {Left, Bottom}}, {Left, Top}, "ClosingActions" -> { "SelectionDeparture", "ParentChanged", "EvaluatorQuit"}]]]], BaseStyle -> Inherited, Evaluator -> Automatic, Method -> "Preemptive"], RGBColor[0.560181, 0.691569, 0.194885], Editable -> False, Selectable -> False]}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{#, ",", #2, ",", #3}], "}"}], ",", RowBox[{"LegendMarkers", "\[Rule]", "None"}], ",", RowBox[{"LabelStyle", "\[Rule]", RowBox[{"{", "}"}]}], ",", RowBox[{"LegendLayout", "\[Rule]", "\"Column\""}]}], "]"}]& ), Editable -> True], TraditionalForm], TraditionalForm]}, "Legended", DisplayFunction->(GridBox[{{ TagBox[ ItemBox[ PaneBox[ TagBox[#, "SkipImageSizeLevel"], Alignment -> {Center, Baseline}, BaselinePosition -> Baseline], DefaultBaseStyle -> "Labeled"], "SkipImageSizeLevel"], ItemBox[#2, DefaultBaseStyle -> "LabeledLabel"]}}, GridBoxAlignment -> {"Columns" -> {{Center}}, "Rows" -> {{Center}}}, AutoDelete -> False, GridBoxItemSize -> Automatic, BaselinePosition -> {1, 1}]& ), Editable->True, InterpretationFunction->(RowBox[{"Legended", "[", RowBox[{#, ",", RowBox[{"Placed", "[", RowBox[{#2, ",", "After"}], "]"}]}], "]"}]& )]], "Output", CellChangeTimes->{3.741966605920697*^9, 3.741966792886639*^9, 3.7419669050422573`*^9, 3.741966993891625*^9, 3.746203075241784*^9},ExpressionUUID->"a12a28fd-8f5b-4acf-85de-\ bebdf67baec1"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Manipulate", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"solution", "[", "x_", "]"}], "=", RowBox[{ RowBox[{"NDSolve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{"y", "'"}], "[", "x", "]"}], "\[Equal]", RowBox[{ RowBox[{ RowBox[{"y", "[", "x", "]"}], "/", "x"}], "-", RowBox[{"k", "*", RowBox[{"Sqrt", "[", RowBox[{"1", "+", RowBox[{ RowBox[{ RowBox[{"y", "[", "x", "]"}], "^", "2"}], "/", RowBox[{"x", "^", "2"}]}]}], "]"}]}]}]}], ",", RowBox[{ RowBox[{"y", "[", "5", "]"}], "\[Equal]", "0"}]}], "}"}], ",", RowBox[{"y", "[", "x", "]"}], ",", RowBox[{"{", RowBox[{"x", ",", "5", ",", ".00001"}], "}"}], ",", RowBox[{"Method", "\[Rule]", "BackwardEuler"}], ",", RowBox[{"\"\\"", "\[Rule]", ".01"}]}], "]"}], "[", RowBox[{"[", RowBox[{"1", ",", "All", ",", "2"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"Plot", "[", RowBox[{ RowBox[{"solution", "[", "x", "]"}], ",", RowBox[{"{", RowBox[{"x", ",", "0", ",", "5"}], "}"}], ",", RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], ",", RowBox[{"AxesLabel", "\[Rule]", RowBox[{"{", RowBox[{"\"\\"", ",", "\"\\""}], "}"}]}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "6"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "6"}], "}"}]}], "}"}]}]}], "]"}]}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"k", ",", "1"}], "}"}], ",", "0", ",", "1.5", ",", RowBox[{"Appearance", "\[Rule]", "\"\\""}]}], "}"}], ",", RowBox[{"TrackedSymbols", "\[RuleDelayed]", RowBox[{"{", "k", "}"}]}]}], "]"}]], "Input", CellChangeTimes->{{3.741966986702653*^9, 3.741966988901229*^9}},ExpressionUUID->"f704fadf-2098-418f-8237-\ 45609089c550"], Cell[BoxData[ TagBox[ StyleBox[ DynamicModuleBox[{$CellContext`k$$ = 1, Typeset`show$$ = True, Typeset`bookmarkList$$ = {}, Typeset`bookmarkMode$$ = "Menu", Typeset`animator$$, Typeset`animvar$$ = 1, Typeset`name$$ = "\"untitled\"", Typeset`specs$$ = {{{ Hold[$CellContext`k$$], 1}, 0, 1.5}}, Typeset`size$$ = { 360., {153., 159.}}, Typeset`update$$ = 0, Typeset`initDone$$, Typeset`skipInitDone$$ = True, $CellContext`k$1861$$ = 0}, DynamicBox[Manipulate`ManipulateBoxes[ 1, StandardForm, "Variables" :> {$CellContext`k$$ = 1}, "ControllerVariables" :> { Hold[$CellContext`k$$, $CellContext`k$1861$$, 0]}, "OtherVariables" :> { Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$, Typeset`animator$$, Typeset`animvar$$, Typeset`name$$, Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$, Typeset`skipInitDone$$}, "Body" :> ($CellContext`solution[ Pattern[$CellContext`x$, Blank[]]] = Part[ NDSolve[{ Derivative[ 1][$CellContext`y][$CellContext`x$] == \ $CellContext`y[$CellContext`x$]/$CellContext`x$ - $CellContext`k$$ Sqrt[1 + $CellContext`y[$CellContext`x$]^2/$CellContext`x$^2], \ $CellContext`y[5] == 0}, $CellContext`y[$CellContext`x$], {$CellContext`x$, 5, 0.00001}, Method -> $CellContext`BackwardEuler, "StartingStepSize" -> 0.01], 1, All, 2]; Plot[ $CellContext`solution[$CellContext`x], {$CellContext`x, 0, 5}, AspectRatio -> Automatic, AxesLabel -> {"x(t) - miles", "y(t) - miles"}, PlotRange -> {{0, 6}, {0, 6}}]), "Specifications" :> {{{$CellContext`k$$, 1}, 0, 1.5, Appearance -> "Labeled"}}, "Options" :> {TrackedSymbols :> {$CellContext`k$$}}, "DefaultOptions" :> {}], ImageSizeCache->{411., {202., 209.}}, SingleEvaluation->True], Deinitialization:>None, DynamicModuleValues:>{}, SynchronousInitialization->True, UndoTrackedVariables:>{Typeset`show$$, Typeset`bookmarkMode$$}, UnsavedVariables:>{Typeset`initDone$$}, UntrackedVariables:>{Typeset`size$$}], "Manipulate", Deployed->True, StripOnInput->False], Manipulate`InterpretManipulate[1]]], "Output", CellChangeTimes->{3.741966994979238*^9, 3.7462030778639345`*^9},ExpressionUUID->"e005cb13-0ca7-4dc5-933c-\ 721a0e02d4fb"] }, Open ]] }, WindowSize->{1139, 1198}, WindowMargins->{{573, Automatic}, {71, Automatic}}, FrontEndVersion->"11.2 for Microsoft Windows (64-bit) (September 10, 2017)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[558, 20, 1223, 25, 105, "Input",ExpressionUUID->"33eab7b7-2a12-40bd-9558-324e8e2e4bd0"], Cell[1784, 47, 179, 3, 28, "Input",ExpressionUUID->"01917e4f-2166-4ab9-bb29-a631abc026bb"], Cell[CellGroupData[{ Cell[1988, 54, 3754, 113, 124, "Input",ExpressionUUID->"292841af-fbb0-4901-b982-bd43d688dd7e"], Cell[5745, 169, 390, 9, 21, "Message",ExpressionUUID->"7d56a2cd-5bba-40e7-bb96-b3d811431c1e"], Cell[6138, 180, 406, 10, 21, "Message",ExpressionUUID->"b120db42-8032-4754-9f6f-0ccbad6ab90f"], Cell[6547, 192, 390, 9, 21, "Message",ExpressionUUID->"be4d9328-c3a6-4702-b76c-aca14fc589e0"] }, Open ]], Cell[CellGroupData[{ Cell[6974, 206, 1006, 28, 48, "Input",ExpressionUUID->"6dd662b4-c2ac-4ef0-a31b-ed2f6439eb1b"], Cell[7983, 236, 439, 10, 21, "Message",ExpressionUUID->"9d487384-ba43-4ded-bb01-711285ad336f"], Cell[8425, 248, 439, 10, 21, "Message",ExpressionUUID->"ad5979c9-f74e-4482-8c87-892c8f3dcaed"], Cell[8867, 260, 32693, 600, 331, "Output",ExpressionUUID->"a12a28fd-8f5b-4acf-85de-bebdf67baec1"] }, Open ]], Cell[CellGroupData[{ Cell[41597, 865, 2226, 61, 105, "Input",ExpressionUUID->"f704fadf-2098-418f-8237-45609089c550"], Cell[43826, 928, 2453, 51, 430, "Output",ExpressionUUID->"e005cb13-0ca7-4dc5-933c-721a0e02d4fb"] }, Open ]] } ] *)