Mathematica

(* 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[ 100332, 2147] NotebookOptionsPosition[ 98468, 2084] NotebookOutlinePosition[ 98811, 2099] CellTagsIndexPosition[ 98768, 2096] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[TextData[{ StyleBox[" Lagrange polynomial. ", FontWeight->"Bold", FontColor->RGBColor[0, 0, 1]], "To construct the Lagrange polynomial \n\n\t", Cell[BoxData[ RowBox[{ RowBox[{"P", RowBox[{"(", "x", ")"}]}], "=", RowBox[{ StyleBox[ UnderoverscriptBox["\[Sum]", StyleBox[ RowBox[{"k", "=", "0"}], FontSize->10], StyleBox["n", FontSize->10]], FontSize->18], RowBox[{ SubscriptBox["y", "k"], " ", SubscriptBox["L", RowBox[{"n", ",", "k"}]], RowBox[{"(", "x", ")"}]}]}]}]]], " \n\t\nof degree ", Cell[BoxData["n"]], ", based on the ", Cell[BoxData[ RowBox[{"n", "+", "1"}]]], " points ", Cell[BoxData[ RowBox[{"(", RowBox[{ SubscriptBox["x", "k"], StyleBox[",", FontFamily->"Times New Roman", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontVariations->{"StrikeThrough"->False, "Underline"->False}, FontColor->GrayLevel[0]], SubscriptBox["y", "k"]}], ")"}]]], " for ", Cell[BoxData[ RowBox[{ RowBox[{"k", "=", "0"}], StyleBox[",", FontFamily->"Times New Roman", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontVariations->{"StrikeThrough"->False, "Underline"->False}, FontColor->GrayLevel[0]], "1", StyleBox[",", FontFamily->"Times New Roman", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontVariations->{"StrikeThrough"->False, "Underline"->False}, FontColor->GrayLevel[0]], StyleBox["...", FontFamily->"Times New Roman", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontVariations->{"StrikeThrough"->False, "Underline"->False}, FontColor->GrayLevel[0]], StyleBox[",", FontFamily->"Times New Roman", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontVariations->{"StrikeThrough"->False, "Underline"->False}, FontColor->GrayLevel[0]], "n"}]]], ". The Lagrange fundamental polynomials ", Cell[BoxData[ RowBox[{ SubscriptBox["L", RowBox[{"n", ",", "k"}]], RowBox[{"(", "x", ")"}]}]]], " for degree ", Cell[BoxData["n"]], " are: \n\n\t", Cell[BoxData[ RowBox[{ RowBox[{ SubscriptBox["L", RowBox[{"n", ",", "k"}]], RowBox[{"(", "x", ")"}]}], " ", "=", " ", FractionBox[ RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"x", " ", "-", " ", SubscriptBox["x", "0"]}], ")"}], " ", "..."}], RowBox[{"(", RowBox[{"x", " ", "-", " ", SubscriptBox["x", RowBox[{"k", "-", "1"}]]}], ")"}], " ", RowBox[{"(", RowBox[{"x", " ", "-", " ", SubscriptBox["x", RowBox[{"k", "+", "1"}]]}], ")"}]}], " ", "..."}], " ", RowBox[{"(", RowBox[{"x", " ", "-", " ", SubscriptBox["x", "n"]}], ")"}]}], RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ SubscriptBox["x", "k"], "-", SubscriptBox["x", "0"]}], ")"}], " ", "..."}], RowBox[{"(", RowBox[{ SubscriptBox["x", "k"], "-", SubscriptBox["x", RowBox[{"k", "-", "1"}]]}], ")"}], " ", RowBox[{"(", RowBox[{ SubscriptBox["x", "k"], "-", SubscriptBox["x", RowBox[{"k", "+", "1"}]]}], ")"}]}], " ", "..."}], " ", RowBox[{"(", RowBox[{ SubscriptBox["x", "k"], "-", SubscriptBox["x", "n"]}], ")"}]}]]}]]], " for ", Cell[BoxData[ RowBox[{ RowBox[{"k", "=", "0"}], StyleBox[",", FontFamily->"Times New Roman", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontVariations->{"StrikeThrough"->False, "Underline"->False}, FontColor->GrayLevel[0]], "1", StyleBox[",", FontFamily->"Times New Roman", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontVariations->{"StrikeThrough"->False, "Underline"->False}, FontColor->GrayLevel[0]], StyleBox["...", FontFamily->"Times New Roman", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontVariations->{"StrikeThrough"->False, "Underline"->False}, FontColor->GrayLevel[0]], StyleBox[",", FontFamily->"Times New Roman", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontVariations->{"StrikeThrough"->False, "Underline"->False}, FontColor->GrayLevel[0]], "n"}]]], ". " }], "Text", CellChangeTimes->{{3.5884331150297766`*^9, 3.588433134627365*^9}, { 3.588433901117075*^9, 3.5884339293098936`*^9}}], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{ RowBox[{"Here", "'"}], "s", " ", "the", " ", "implementation", " ", "of", " ", "Lagrange", " ", "interpolation"}], ",", " ", RowBox[{ "which", " ", "takes", " ", "both", " ", "x", " ", "and", " ", "y", " ", "coordinates", " ", "of", " ", "points", " ", "as", " ", "input", " ", RowBox[{"list", ".", " ", "Fundamental"}], " ", "Lagrange", " ", "polynomials", " ", RowBox[{"Lag", "[", RowBox[{"k", ",", "x_"}], "]"}], " ", "are", " ", RowBox[{"plotted", "."}]}]}], " ", "*)"}]], "Input", CellChangeTimes->{{3.588615210991409*^9, 3.58861525733857*^9}, { 3.588615361929764*^9, 3.588615395730233*^9}, {3.6198896419932528`*^9, 3.6198896580140533`*^9}, {3.745921773781764*^9, 3.7459217762516546`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"llist", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "3"}], "}"}], ",", RowBox[{"{", RowBox[{"2", ",", "5"}], "}"}], ",", RowBox[{"{", RowBox[{"4", ",", "2"}], "}"}], ",", RowBox[{"{", RowBox[{"6", ",", "7"}], "}"}]}], "}"}]}], ";"}], "\n", RowBox[{ RowBox[{"plot1", "=", RowBox[{"ListPlot", "[", RowBox[{"llist", ",", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{"Red", ",", RowBox[{"PointSize", "[", "0.05", "]"}]}], "}"}]}], ",", RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]}], ";"}]}], "Input", AspectRatioFixed->True], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"n", "=", RowBox[{ RowBox[{"Length", "[", "llist", "]"}], "-", "1"}]}], ";"}], "\n", RowBox[{"Do", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"x", "[", "k", "]"}], "=", RowBox[{"llist", "\[LeftDoubleBracket]", RowBox[{ RowBox[{"k", "+", "1"}], ",", "1"}], "\[RightDoubleBracket]"}]}], ";", RowBox[{ RowBox[{"y", "[", "k", "]"}], "=", RowBox[{"llist", "\[LeftDoubleBracket]", RowBox[{ RowBox[{"k", "+", "1"}], ",", "2"}], "\[RightDoubleBracket]"}]}]}], ",", RowBox[{"{", RowBox[{"k", ",", "0", ",", "n"}], "}"}]}], "]"}], "\n", RowBox[{"Do", "[", RowBox[{ RowBox[{ RowBox[{"Lag", "[", RowBox[{"k", ",", "x_"}], "]"}], "=", FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox["\[Product]", RowBox[{"j", "=", "0"}], RowBox[{"k", "-", "1"}]], RowBox[{"(", RowBox[{"x", "-", RowBox[{"x", "[", "j", "]"}]}], ")"}]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox["\[Product]", RowBox[{"j", "=", RowBox[{"k", "+", "1"}]}], "n"], RowBox[{"(", RowBox[{"x", "-", RowBox[{"x", "[", "j", "]"}]}], ")"}]}], ")"}]}], RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox["\[Product]", RowBox[{"j", "=", "0"}], RowBox[{"k", "-", "1"}]], RowBox[{"(", RowBox[{ RowBox[{"x", "[", "k", "]"}], "-", RowBox[{"x", "[", "j", "]"}]}], ")"}]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox["\[Product]", RowBox[{"j", "=", RowBox[{"k", "+", "1"}]}], "n"], RowBox[{"(", RowBox[{ RowBox[{"x", "[", "k", "]"}], "-", RowBox[{"x", "[", "j", "]"}]}], ")"}]}], ")"}]}]]}], ",", RowBox[{"{", RowBox[{"k", ",", "0", ",", "n"}], "}"}]}], "]"}], "\n", RowBox[{ RowBox[{"a", "=", RowBox[{ RowBox[{"x", "[", "0", "]"}], "-", "1"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"b", "=", RowBox[{ RowBox[{"x", "[", "n", "]"}], "+", "1"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Plot", "[", RowBox[{ RowBox[{"Evaluate", "[", RowBox[{"Table", "[", RowBox[{ RowBox[{"Lag", "[", RowBox[{"k", ",", "x"}], "]"}], ",", RowBox[{"{", RowBox[{"k", ",", "0", ",", "n"}], "}"}]}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"x", ",", "a", ",", "b"}], "}"}], ",", "\n", "\t ", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"Table", "[", RowBox[{ RowBox[{"RGBColor", "[", RowBox[{ FractionBox["k", "n"], ",", "0", ",", RowBox[{"1", "-", FractionBox["k", "n"]}]}], "]"}], ",", RowBox[{"{", RowBox[{"k", ",", "0", ",", "n"}], "}"}]}], "]"}]}]}], "]"}]}], "Input", CellChangeTimes->{ 3.5886146754875093`*^9, {3.7459417556599874`*^9, 3.7459417766151714`*^9}, { 3.7459418787159576`*^9, 3.7459419132733574`*^9}, {3.7459419489199905`*^9, 3.745941978035424*^9}}, AspectRatioFixed->True], Cell[BoxData[ GraphicsBox[{{{}, {}, {RGBColor[0, 0, 1], AbsoluteThickness[1.6], Opacity[1.], LineBox[CompressedData[" 1:eJwVl3k4VV8Xx433IknXvfiRIVNEKDRQ9kIJSQqVUlFSMhRpUJSkVGQeUxIS oswU2hnvPScKkYxFCBWpkCnvef86z+fZ+9l7re9a33XOWXn09J7jPFxcXA3c XFz/fz7a/7T0lsZ1w2GDg06JiYKQPisemBbnhOQCpBKLQwVhZwi/mn+cF4JZ mWdtVwUh9b3fvGvcNRQf2eyx0kUQGi5sdHGMi0Tce5fofV0nCH3m4kkH4x6h xQDCTb1RAL6vu6FnH5ePLF94V0fxCYDXkwqe/XFVqCQgxBZfoYGmZm5igU0L 6krWTPFx5gV1PHHmin4XOs32c35TzgXcwU1uO2b70JSsjt/cs3mUFbNnA3/J EFKrhtS+P9Mozn9Xap/8N6TomHJ6ydRvhBvOP3YbHUeN9e5F4aI/kVul2t9G q1/otxRsarEeQZKBdh7TV/8gmw1/j4w9/YK++MY9u/B5CqlHrVQuju5G1p9i 0G3JGcRop91N/fAe6fzx+OOxbQ4pzTg5nujlINtcr/5eqzlktCpXmZbGQT4u 55p27Z9D9wNyGhpPcFDRR7+ctW5zSGRZUWz7BBvpVoY4T0bModLlh368p1N8 M6vVr3sO7XImr7gY1CFdycGikLPzyPPn8+Cs11VIb7PDuaxHC2jXcdUtRYEl SPhbcubF7AX0Vke0316lBPUn9nVtL1xAuU/qFY3JYhQx7WI8VLuAJP4Zbepk FKPRAq9lisMLSKfnS6PKk0KUrBacdV/zH4qXccxS6M5DdImCnvDyf+jYoPaS g2MZqGOCbnq+bRE5PnM/ot5zHvlNXmrl+bSIjl3O87cr9ELyMz+Ohg8vokr+ hXGVFa7Ihav1aubcImoTCfbTuGOMfOlNDj9EuUD0fUhj2oQX5he/Hye3gQsK DJsqxt6GYhkdXaGg61xwdte8jWZhCh4r8thtcYsLpDvnZzxMH+HX658kiN7l goSr6bHyHY+ws77UqgdxXDCpIuLwaT4VZxtxbS3J4oL+zVbD7C2P8XrrBv/h d1wgpHdf3O1xFrbycJ6wXMEN1cdTt8XS8vGV9MHeeXluEH4xI3bENB8/73Jp yFHmBmJhx52vN/KxqIXrE2EtbpC1NKnS5SvALSqehxqMuOFy5OrE6oUCbPf5 ArHjBDf4njCbqhwowg42IY8sCrjhdZBKZE1gGQ69Ixw+W8INTsmWpndeluHK qrt+2eXcYKDlv57xqwzLakfsE6rjhtVPjki5OL7An4Vjl5Lt3KB2o11bXv8l Plaf7Gu+wA1yqYsPxvvK8alNBbvNzHhAKT0uPYQX4/PfWpyeWPLAY97Cn8e1 MQ588NuLtpsHDL8VfV5yCONEbr2oWnseqLXePP6yGGMOUdoCbjzwW2dZo826 11j5QKXNpjAe8M192u/c/Rp/vsSxW93KA81+7mttZKrxd42R47c/8sCm0u+1 d7Sr8d9ewfPD3TzA++ypXpZJNV5msiMuY5AHqnftarrnWo0Nl7z9oDDNA1t7 lrXuLqrG95Pe75OW4oXs0dvO4ttq8L6KXnthR14IWcmloGtbizkJ8RIk5eMT HXOzhcdrsf4567ZgV17wxYn+EhdqsaxmtTXvWV4gzhzKvp1Yi4cePt4+e4MX UL9Hg3RPLT5/3V1v+CkvvGz+d4PrSB0ePqL8+3EeL1wWov3FnnX4wObevGPF vNDQpfvq2JU6bDi5S+PTK2p/arbX/gd1mHZCR/FDMy8MGDg+nO6ow3EWs6K1 01S87shV3KoeC6wqfHttnheMGzJXyDvU40u87qGImw9ITfxW5FQ9PlzRQy8X 4gOjfTH2STfqsYpm1b98GT74+9xc17u8Hhcvv/X9oQkfvNitsVtegY1bP7I4 l8P5wOqzVg2ni41XantmDkfzwRlHxaCmQTb2DK6/ZZfAB6FPaoaqxtlYcP0F c61HfHCyuKjQiZeDDaM/kv0FfNTcXTwcp8bBmTuTGs0/8MGa6z4RXV4cPJn+ O7ekkw8sfJT9rl3iYOP5HWGKn/hAv109Zvl1Du7Jnrda+MoHu2S+uP2M5mCG wOGmvBkq/rbUnKwiDr5cLf9eQoYf9CxTBBi/OJjzn29h0Ep+mDmYSBye4WBx r+boCWV+cLsVYRTHReA8uUDbN5r8wDiq6FsnQuABvy9tV4AfTpoEziivJrDV hoyPg8f4YZ344dYpBwInhf0r23OSH8x+qbF+HyXw8ODeROzOD8Ox9+M7TxL4 egz9QOI5fvA7LB/k5kPgsomTXZbB/PBDznws9g6BFXLUewuz+eFn+f77zEIC f1izwv3qc35w2uJVp1FG4NvPhWctivgho1t5w5pKAv8s+MHqr+SH/Bx/+a/1 BH714vlO0WZ+IKQeIuggsLdBSldXGz98ZVj3J/UQWLkywvVJJz90JDQrf+4j cOhrrxuGA/xgdNOxXGGUwPZsnUr3aX6QP33cLnKGwMLmSjs2zvPDxKE9OjoL BH5NMjt4uWlw32Kr0CsuEqu8nfxzT4gG3na734XRSdy5ayjQZRkNUsN937cK kfhu8wfRdUwaZC6+nucVIfGf1lINQoYGhi/YozJMEmfuzXwZo0CDRtmYOgEJ Eh/8mGDmuIoGur8SJ7r/I3F1l6/ztDYN7kq5uxrJkTii3+C+8lYauMz0Cc+r knirs8bqCTMaVNitzVinTuLpwRVlFTtpMBo+NrxnDYkPjyy8t9lHg/dtDXHW a0nMcBtzknOgQfjCrISmDonrvveOjzrSQMkg33NSl8QaP/GSwFM06NZq4hht JPFnr7zEnadpwMtr2l+/icTRv1NW/edDAy6lwwt6BiQ2PRdZPHCRBh3ycsrh m0k8O3XNJM+fijc43LllC4mfXfRuvhxIA5mh4vpFQxI7zR49sj2YBjfJhL0S QGKmn80PRigNyGpDuf+MSMxZMLncG0HF7/BMnd+YxH5XdQWzY2nwVWXsZhfF WtzK8efu0aDyGd/m+yYk7g9kKRs9pMHpwL9227eSOI6PViicToN/Wk2fuim2 uDkFHzNpcHJb1IDDNhIv0L++TculQZ71ttNsivNutzucLqD0+TR+Z4UpiZ2X cEb1S2nQEpBg7ECx5N2yi7QKGuRWbXt0k+I3Ilm0lteU/tr/8pMovhKRGPOg jgazBzhX/8/rGHcUXEkaeLzLZv1//2D0pTzddzQol865cpDiBJabIVcrDUYa W+qkKbaMP9jw5iMNJs4qT9VT8SxKWh6I76HBx6RsqUMUM6yDwiz7acCZLKr8 TOWjFFxRzf2VBk52pzSsKV7/6s9UyTcaXB3a5PWU0sNsUkPd/ScN5lqj4n5R eh3QOH5k5SQNSiwHMpQp9jj2IPrDDA1u+z7I2UrpHdW8dB546fD1EW+1KSJx uoCp9hSdDhqTYUNqVL1K0RXnp8J0mM4xXD1D1bc7d6yBJU6HxBn2Rgd9Eo8N qHC/kaKDytBVyZ9Uv3CtOKIXIEeHnCFJozMbSKx8592DUVU6LHufHadH9dfG KnrLwzV0eCwe+PHSOqo+fxHNbh0d0mSmQ3O0SezpkueJDehgpnzhQpsGiQMe DKf6AB0y/B2rG1ZT/dcq3662jQ7dZulZeVT/lxlHohgrOsDX6fOGSiQmfYmz 5jZUPIwhr+GVJO7J4878t48ONke+rbtK+YdHzlvU1YkOU8v+ap+VIvGOuT19 W3zocCrO2Oj8Msov60JYvy/SIfkhIzJYmMRnXGvMM/3pMKN2dvaGIIlj2nUK GMF0SPBKNbHgJXFvIfP610Q6+Cf0v5+YpObRqGXp/WQ6/Dt/hOnxi8C8Cje+ 7U6jQ5DeUFXbGIFVIyZtKnLoIMnz9u7prwT2cv+gFImp+x9rLa1vJ3BQmoi9 aS0ddKyVOnPfEziu0/TuHIcOA5vun7j+jsAvzUonj7fQwcO35cc0Ne/4lBPq 9QfpIHfu2ht6MYHje+xPDggJQDjW+vHmLnVe+FpstkwA4oUU34zdouahkaB4 rpgA7KtW2M8TRM3vx2U151YIQNaN1Y/mfKnzPCTk6JoCIF0e0b3hOLV/vrVN bbcAoJ6jR3s3EtjxWY5GmJ0ArL9bpiaqQ+CdjkHXf9kLwKXPGat01lD51eqs LT8qAFFXtEcsVxL4U0hUqKWPAFz+q3GmSYDAltLWJqfjBcD807lp2TYO1m9c de99kgCcOTQhf/ItB6+6uvhzQ4oAXCyOFEhjczDPl2fJXFkC0MSlpDv2goNf ZC+djXwpAGrOtSOdDzhYRf9NXlGPAMTZP9zIOsbBXPamsrOKgiDymVVg9oWN bfVKFCpVBWHJ3i1Zip1snLlcZVXAGkEwWHNJZryJja1Jfm3aBkEw3+Gcue8V Gz80qDNiWAhCyAHL82cS2NhQduvx1WcE4YrrvDNpzsZ+X4xyDlQIQtgJJwW3 R/X4r8cW/XJbIRhtP1FVoVKHV9VZOq2xF4IszP7+RaoO713hcOvhISEIz7b/ vihSh4vISx+uuwiBial+iPh0LfZaVeZteVEIoioWtuD6Wjz6ed3TniQhqOql qdKca3G3jeoKri9C8JnzLfJJXA2u2iT2b6vXEnDLs+pS6q3C2w59sNl4SxjE qgd3tOlUYHup9GXtJUtBeDRZyEe6EBuGWZ81bBaB5WfHwt4y0nFH4/XkTG5R aN9UcbjVajPy4trbMBQvCs075+eclNKRd5iElaDScghzKZd9xShE+q69gdwV y8FiWd2BW3oVyO1hzJJecwa8JO89MB2rQn/2dFkLWTIgYkq4Gs1XIT+aQtx6 KwZoF5hv0haqRqEez+XD9zCASbdanFKuRjkGpB4cZEDTqTxe6UPVaKx98Uiq OwMM9axwMVGNvEXdilzCGdDw7Mrd+aQaNFuTPxMVyQAzvTVKTZk1KPDCX0Mc zQD1ykXaveIaFN17g5BIYMDhRYMI0Xc1qCgnpZeTwoC6CPQmmLsWTZl9EFAv YACHZ7Ao3rkW+V4zOjTeyoDE+U39vIrU/5aF7YxFOwPu2WdCqFYd+sk4EZvR wYDx5hPz9M11yCU9tPFQLwN2Fn116rCtQ7bs9i0NXxlQlmyB+2/UIc2lnrJZ swzq/R4fd3OwDg0k3Pt0VF4MjJJVdxH36lGKU+7lVwpiMN1oQ8t6XI8Orn4t IaUsBpcktry5lFePWsoHdzWriUFN/D7bP3X1CPdqv0Y6YsCVs8TQ8Gc9uqfE TllhKgaEpJuxkjEb7cr77fTBTQwynvdmjHey0cX1e8MFPcXgzPglgaNf2Cil oqxi8xkxcM69p8D+xka/OP4SaT5ikNC557DrPBvF9dEbT/uLwUZfNf6lMhzU y1ixUTBcDE6/DrlPHOQgeuKV45sjxUBqp8WtY8c4SFuuL+p0tBj0r9uWPn6K gwLVM360xYvBLei16PXlIJWt2mmpKWKw7Yu8MzOegzzPbRXZXCAGzxGfTfdb DoqfyzA4XUTdx/HeGtvGQa+vCbqmlohBY9CHHtTNQYywtzUC5WLQ3eQ25DrC QSUZ+33basQgW2e3oygvgf61uw94tomBlozi9Yc6BFp1+N3y1HYxiLJ+N+m5 iUDWA2tRW4cYTESp3dJCBEr9OZVo0Evlf/nDyHULApkKXtsl8FUMNhyILbjo SKC7+rEvH/0Vg5dmNYN3bhNotK3GkG+Out80d3dTGLXf61e1y4IYjBooPuKN oe5/suuNOg8TDjfvwDrJBPJkCXUXL2HCC4lKTlk+gci8jY6SIkwwM9Bd6V5C IBXLEwOXRJkw2n+StrScQL2Btd8RiwlzQWk7JGoJZPXzygIhywS7wJBIyzYC ZYfkXtNYyYRA3RXybh0Eoq3q5g9XZMLwdNLA+R4CvTq0ScRWlQl/LnIP2A0S SOvNb9netUwQbIkM1f5DoFAXhTTQZQK/H+de6TSBhrl3r0pbzwTHWrlb6nME Stn4TOukARNqpjitHdwkEs04Cb+2MmGXtQ5tlwiJ3I3ia223M8H04P0eg+Uk 4nTXmZWaM6E6ct8WFpNEAWKKu/2smNC9/dCpu/+RqPvZ7rZeaybUrv6zUmMF iTZaBNgb2VDnL55bXypLoomAHifafiY8MH1vEaVIIqsVwkMnDzChSnpOt1+Z RNml+qfeODAhPCvDWlaVREfH4r0jnJhQWbxy4oAGiV7drp/6dYwJJ+O3OB7U JJGU8uQlOxcmWP8917Zdm0QtB/dcl3Jjgvu/1owBHRJpTgfQ/T0ovWlmP2L0 SHQn6nnIp9NM2Gh3RmTtBhIZE8Ix6T5M6LjtW6yhT6JkZwNJ+gUmnH+foRxq QKLZRdf7rr4Uqylsa99MooL17HTNK0zQrUpy0UYkEmmZVI0MYIKBy9AeQyDR KQ+l3N+B1Pm9EjUbjEhUL2Czdu8NJrQc2JuzwphECunXisuCmZC9lMP7g+Ir KG+T9B0qH9XoV5kmJOrs7K30D2VCXndP866tJFp/fqnx5zAmKDm8XddPcfTy zfXGkUzw+Hjnw5FtJHI49bLLL5oJa0o0jldTrFKzcaIkllrnN2xgmJLop3Qp bSKeCeayIzy7KX7po7dC/R4TPu7fz+VPcVBj4drj96l4l5aUx1G8U2Xd9ofJ TLgaraGRTLH41TyHjhQmNJgMmUVT/Lld01ssjQlUXILnKc7Wzg3e+ZjSL7Dh mCnFPrfVHwQ/YcLbqp3WfBQb9mcVVGUxQSekvO45FR/dQJUz95QJMhug3Izi puiMHr1nTHgqQ1N9R+V377vS79N5VD/HbeA1odh5W5pAdgETWHP/TNIpfTST V8oOFDFBJfVmH/X9iv5OPdSRLWWCycDfFk2Kq3fJmu9/wYSb/wUs20/pH5p5 /3BUOeWHe6ZhnlR97LilfRoqmTBSed7Um6qf7IHE27TXVP169LWOGZJouEDi IVQzQda6EhltIZGfM5Msqqf8GZFRVEv1h2ll1KcxDhPUlh2ln9pE+UV8+aTq GybsNrpxeJ7qpzT2UvkH75gwuTJn6qsu5R/5UL32ZiZkBPCpGFH9uN5XaMfy Vmo+/PZYH7yWRORq+vkbH5lw5oJMe8caEsVcvxGCO6l6bSyzG1In0eFu3kcz 3Uz47rE06pMa5Z+7XA0efUzgTYiQT6D88nLoSt+TL0zI2izmtI/yUxBamOob ZMKNkhkD3pUkkpyYWbl3lAlb7yQ9kKf82G9+cUPEd2pexTxlRFF+zUmdsiTH KH1bOf/GxUmEbH9fMPzNhEbbVkkvyu8upd8bVeapfv/X5BDMTyJtUbcvjv8o /fd/jzrIQ/nj5Mjfe1wsuJN5YkZ6kZqfUkOKy/hZwLi7StfpL4GK/D/5Ti9l wT+6p+v8CIF4trYoc+RYcKhbQdKTINDlCz4i0ytZcOCjZ7xEHYH+ZItPKyux QKrZIij3NYG+Lj/Iua7Kgl1Vx28klxKo4dOXU2gtC+zbhTfHpBMo9vLk8xJj FpzLM1/HvkQgkefxCYNbWVDvE5Gmc45At/o3XWNup9b36SyEnCbQJfMre7x3 sMAxAMy4nAl0WII+ucaWBTJJh2yndhBItVBS//FxFmCyYU5UkkCPhl4qtJ5g gZxX05KNywkkJXVoCe8pFlzYN3/FQohAS6896nb0ZEHg96sea+c56JfV6qsr LrCg0Vs6kPaJg8pHDWqjb7GgrL4R6lM4yErhiFXgUxZ0Ic9GJyYHzZ5qvPs5 lwUJip0e4kIclFFo0GiYx4I12DU6f5GNFrdKWs4VseCWjY9J0igb5Z1oNvd5 xYJvJ8uWGmA2Ess12ebSzIImPoHF1GNs9HoyP6j2PQuS84qnXu5jI3dD+VqF DyygfwsmynawUe27eeNPnSyqX6U+HNdho/O/imH/AAv8vl1wCuJmo48bVDdb TLNg9I/CJJFQj4IC4i5nzlD5tMql64fUI22Cr5w2zwI+KeOgCP96dPtA36Ya LnFQFrqQ+82pHun739uwWUgcBGYULx1Wq0cPqoV1NGXEofJdbJ5GYR06unNC jWEiDl22G9atLKhFfjbPZ7dvEwdX9SWupqm1KM7e443/dnE4y3SQtYmqReTx EfeRHeJgqD3cp+Bdi9Zd6Xv+2lYc3ETedd7TrkXcz5v1TruIw9+jvPyD1Pdf imiBUcNtcTA56FL7+GY1eil+hsETKg5i+xNFss9Wo9YVml82hInDgZm+6hjH aiSglh2UFiUOtw1WbZfRr0ZnjFLZl5LEYa3hk4aZb1UIeUftVMsVh4v6WZs7 LKpQ73tv+5tN4pA0q+YzOoKR+vh8wLcWcbheczcv7wVGF4WCM63bxKFK6o7F xtsYMYySpqU6qfueOaCQVRiZPquJff5FHOa4umq3OL1Cz28xWz5OiYP9EoKn mlOBAraUWKjLSEDqmSaX3BMvUON+8I6Qk4B0DYVXPlovkJQPmTi5UgKEzcIK Z6fKUFF27zBWkQCt9jp54xtlaFiCHmyrLQG5PxQfVCSXot2/9tX4m0jArwLV Psu6YqTwZMagyVUC2AdeFsDnfCT74Mn61e4S0OwCSz6l5CPpaLu1QZ4SUJv8 qszAKR8xr+arbDwrAdeOuWOtvjxE2+e6PMVPArRHWKe4Kp6jUVrH0Olwar34 upRT3FNU4FIWuaxYAjrr6amXhlKRicqFr1ZckvCv0uNy5RJl/KzD52wAjyTE bltcr3PMDkve9V7M55OEFQOHH2pIOePvvz0kWIKSEBZ70X3LtbM4psp5e9dy SdiRH/tmfsd1POSw58kJRUn4+Hp9ML9jNL4VveZEgKkkuOv8OWnGScONXAND +aGSsNu9ISKhogifInIJRrgkjHY2iW6RK8YCURdyzkZKwtDiuIBWYDE2Vlzi rRdHxbczqXTOrASXbVu3UPZQEnJVz7jqd5bi1DsBYq8LJMFtk/PkkSXl+DxT Bt5+lISpMwHaXq2vsZzq3sRRpf8g7FbB0spRNu7ILn1IL/wP7FIehLd4v8de HctO2xtKwax4c3jM7y7cZOsxEdAmBf9O/uCKNv2CE9RfskOPSoN1N31129Aw /jjiE8NzXBocRJPZV7lHsGSmltPFE9IwtmO19sYVIzhR6fHsMXdpeLBNxLNl 9wi+JxOxxuC8NOiud/mZUzmCk5a5RI/clobyjOPZ8/GjOPn38iOm+dJgmCDX IL3vO+7Nb1CvKJSGOq7guzfOfseyZ4L/ri2RBv4mLd+F8O/44ff5SJlyaUgT 0eqSIL7jlMGvtX9qpcF1cZvMtMEP/Ki9cnXaR2mINR/UPao0hvtjL05LdklD SXXl4nOjMaxgq1Mb1iMNcnv8lYSPjOHUpsxDl/qloVEuKXg8kWIiOmL3d2kY dC2kHxMdx/3BVofYY9LQfktlUlFzHCuYCq7eMiENw7E2tr93jOO06is1qlPS sH+o9mxl8DgeCNgUkfxXGrI/LWwqezyOldAfB+acNPxZ1+lfVTOOnReeqd1Z kAYDA1vjjr5xnF7uOrW4KA3+Ni53FhfH8f8Av2+jfQ== "]]}, {RGBColor[ NCache[ Rational[1, 3], 0.3333333333333333], 0, NCache[ Rational[2, 3], 0.6666666666666666]], AbsoluteThickness[1.6], Opacity[ 1.], LineBox[CompressedData[" 1:eJwV1nk8VF8bAPCZQZZkm2vIIHuWX5YIUZ1Hi6KyFJF2ZaeNKFHSRilrSHZZ UlISFToVLUTZdzOjJCGmsq/vff+az/cz9zz33Gc59yo6n9jlQqNQKOVUCuX/ v+8EFuN3STWhreqvbtzImcda+M/JCyZdqPj7v1ctPyYw9Xq91/aZXnR/1NXq tzMXP4jbZcRX0o9eMwvshAz7cHywdWavwhBql1U0zvzYhHGtf7bX4Cj6d/aP 7v6NftirQmOqzuovunp0PNGaUo8ueC+kxDj8RatenlutblyPophNmxwO/0Xu Bf4B3SfqUcm5oEjOqb/oRY6vcjurHlEMv6r9i/uLaq3zObJlDSj+sZ+9dNdf 9Enydsc37yb0Lg0XObv9Q3loh3ji81YkHWrvM3lxDL0ICZ8+0tONJq86HT8X PoaiNSNM9or1oNbwQydmY8bQvQ8jDPVNPehOtMepxewxlFAtV7MjrwdJZASf 4a8dQ2oP1A2enWShpW+yg6Wkx9E4WvFzbpqNFubHbxk9GUc/8+dVXAJ60fdz 8Y8DOBMoqMTl+hu9PnTWt3PM4dcEWmxJuTe9uQ+J+MibGv+dQENNLtpKjn3I 5HDOxymeSaS9RpBjdqEPRZmXcM6tnES5hGN/WXUf2kBvlQg6PolsQgeJ1AM/ UNIjyYCQuUl0buRyntG5fqSbs7fiMN8Uqq1ZdtH9Rj/6kJbCYyYyhcTatjiE 3+tHf2JUo6gKUyjVQkXuXkU/2nZuzYPQTVNITd7+1CLlJxrfYt915cYUKsBO DV+u/EQ27DgULj2NHrxKkM8OGUD3N7k3rlOcRmqV1QnRkQNoMtfUhasxjQLl ZHRPpA6g1JPfbjiYTqPd1dfq5ssG0DBVt1X14DSaFMhd6TExgMJVPnu/y5pG mqFlL8+7/0LvPKhJs9ozSHLVpNNDs0HE+NL8X6HRDPIxfrK60moQeazOw84w g0IyRbNr9w0isVmrH9W2M+hE3GqP4jOD6GBEsm6i3wx6NbOiqy1vEM08Nvpo 8GoGKWvZfjglPIT0x3zGfLbMIvrwVyxcPYTsCk59Y1nNomcvkzTsW4aQn+uZ emvHWaRxoWD1rd4hVNwe9EjPaxZ5Kx4zrJ0eQgYVN4+NR82iOAu129c0hpHB tQfNQd2ziOf7pa/cq6SlfxTf9J1DgtbFcq9X/0Z2DQOZc0FzaDfNYIFn/W/k d2M4yufaHLrt9J5psvU3Kp775219dw4l6LQqBuz7jQx6aap0PIcs+958N7hM Ol/hTpLQPKo32RBZ+fU3WrNu/5kHGfOI9+mEgM2BESQ8lJp3Nn8exehrFku7 jKBvd3u7tj6bR2vD/HSavUdQ1KTrxv6qeSQTsztHOWgEDRadElUemEcboo7x qSaNoFSN6w+StRdQmL+c0POmEeTXXt3tZbSA/CJPJDzqHEHbrwuLmcICylLm GiT0jqCpvmj/DtsFpJ4X+mTL6AjalZ68iXFmAZWEM2QEhEYRv1RRT2TZAooY 3Wystn4Usd6PiR2qWkBleroPbDaNomI/o83adQtoX+IRz+MWo8i5sTy/lrWA LCuSf1+2H0Xltz4GCNIWUU9dipuUzyg6ydsjfnnbIjpuyXj4/e4o6vjDb+7f sohuPDcJXEWeG0Hjgc009iLSN31v9t/oKFKY/u0cObCImLI2L2THRpErpfli 3uwiEt369UPj/Cg6x1+//7cYBZob1eYsBLiIj5Ecv8KIAjxNIjGly7moNaGl 54cpBaxFx9LeM7kob7mo6iOggIy65uNaOS7aLhf6zNiSAt9kKlreK3JRtIp7 ve0BCmwxGko8psFFcvoGQlcuUyC6q0FzjxEXjRT72FqGkfEWe0+MG3PRG8Pc RLFbFPB+TGfdMuGiYyYyK1PiKRDQKQr567ko34yyueQBBbrutR+J3sRFQZVr bwY9JuPtcDaa2cxFO7f4Nm58RgFJzhnL/eZcxN3Wf/hrOQXWdR3XF7bgIkOb 2uCBr+R+7SQdmVZcJNDI9/5xM7m/nfuOWFlzUcduJHymgwKPQ2LLA23I+A5F SdTvFOC+FfHFu7jo3cHEUpkJCphxUlZO7+GiOE7DAmeGAod7vvFPOnCRi/NS 89xFCmh7h1qNOJLxXS806wtSIcyTb3m1ExdZ+Rz7s0OWCuayYYM6B7nowv0f rDkFKqzROZc9Trqwy7X2kSoVTJpTpZ4f4iIxS49cYR0qvLTS0VI8wkVwaTCu XJ8K1r0+DXWkT73wCvU2poJuSJixnzMXNaodP1BrRgXBA9O7Hh/lItrBUcsg cyrIOSov3XiMi/TvnDT+bzsVnnbU3/pKOo7Hlx5hR4XYwQNzbS5cVGUyRlm3 lwqzhuK/d7ly0dipMyNDB6iw/Mro04+k7TkB1dvdyPVJzW9T3bjomtR0yawX FTYrFopT3bmoxCrw/sOTVKhU2o8OkO6/Ohu97wwVKoK37HxGmlERdHFpIGmx lvU8HlxkPjbvXXaBCrtf+dF3kg7QuujkdZnc/97b9VGk85wp25hhVDj5+WbA F9Ltdy+t+RxBBYLbKrTEk8xvA035fDQVQrM7IoxIrxW4IqYVT4VtNfOLzqQ9 Ed9CZxIV/q0rcA0jneR/behGGhXEu20rc0nXFPB3mNyngvY1e+m3pGf7wj4M 5lGhY7mSezPp/2SFipMKqJB2RvF5L+n9u29mWBZRYdnZCuov0hE3hCNnSsh8 LB6yGSRd8fZWUH4Zma+h01k/SI9MiXg6vaFCg/Th+U7S8rpRDkLvqeDqFn24 hrSVm/iWV9VUcHt+8usz0hdTY1Z7fqHCkV63HQmkn7TQFWSaqPCtmtPhR5oj fGdZTRsVzuzcEriDtNhmxuy5bio8tFzUkSdtdj5hQKOXCiLpZ+d/kfk6XSTd 2vGDfL7/mD8KSWf+ulsZPkiFhDqbXydINyown64dpUKr82VhTdI0x+TUX/+o EFgja8si63X0Q+o5i3kqqN7U2raGdNz8CrdpKg3MM/wW28j6Vxlk2D1YQgMa b22vH2mVrCwdQTEanNfYqZdC9otdp4rcS4IG9g8ss9RJXxXPEfJYToPF0XmL x2S/9V/M+/FJiQbj21ehPLIfGaUaTWdX0uCgfHCMLOmtI/lv1P+jQa2epspN sp/z9hfcC1tDg5E7m4QdyH73XFtku20bDXA+v34+OS/+Q41HcnfQ4LWajuc8 OU+hKf9OLbGlQYr7tifbSd+lromp2kuDykmVUx37uehTdWkjeNFAxtFh7uVe Lmo+3/4t7QQNwqVYvT/I+WWvmv674EsDFzFKvzDpiWgTekUQDVa5D9ruIOdf 1ali99rbNDhapTEaRZ4PusKso4kxNHh4dB8l3paL1r1e8J2Mp8FwIfmtRZ4n dkoQ9zyNBqqCR2zDyPMndPBds95TGnQ2do2vtiTrG/jJXrOZBhWuNTNZwEXD //1yCW8n1ycZdO1HXDTFEvQf6KaB7lFDQfENXCS6aXt8zg8alHm02HuZctGG pV9alSZpsOPoroGfa7go+V6TA1OGBz48rM+jaZL52znmFijPA8H9GTtPqXNR 8SJxtkOJB5yTDxzqUuOiz0f3JCZo8cBFfY+GDGUumtHqaKev5wFP5Uy1KVku cihn7RU+zANnC5R31IqQ+UtMkKo5xgPNn4eO1ghzkckZm5brHuT1ffKoUojs d+13Njy+PHClpdIiewlZ77TsrTNXeUCY83fIfm4U+V/2XjPwkAeexq9+Ftw+ iuItZ8SqJnlg8NqzWIPwUSSw8tmXS3M8IMr7Y/vQ5VEUyOMdgai8cDi8WjPx wig6WN7DXybEC4LWuSebfUeRmvbbhadyvFBolJDXdmAUPRcPG07bxAuOAnbq y3RHUXO75KfzkbzwaDLlkEPdCFLUPZ43EMsLCwpaHts/jqDj1z+E2SfyQtMR LQ+jtyNI0DDAQieDFygjVfsni0fQhtj2mm9FvJAeIXNFKnmEzOe9OotWXlhd tfFIt/sIOv9OoUlKjg/81eTfx839RkqPtFjP8vng5JKXhVnLfqPWVbLeFwv5 4Hib3O6lfL9ReKHwjGUxH+gQ9i8854YRt+i35LcKPljD/NktNjSMXr8s3CnW wAfSbtYrZj8Oo70f9Su8J/lgm0GvyI6QYRT1zTRZdfMSKI67+YFnaAgtSu9w SuhZAvv/CPEOPxhECT173fuEBCC59KVuSd9PRNlrLj+jLAimb4DiutiHpnzW m5TZCcFErM7NjIRe9HYtfWHzqaWQ8tOI7urfjbYcaN1tHCYMzx/hoKK8FrRX 5r5oW8kyuBadF1u5vxZtuG3ju6FBBMSrk/OKZ8pRR93l1DyqGBRt6PNb1r8b n6Lsqe1PEIOSpU4Bmx+V4dO3pawEVcTBvEHUt7S4Gpt4sEKp5eJgprjUgZ7X iOkbqePGr8XB0tYNW5U34mEZFbeTb8SB/+DDhgv1jTil1mM7u0ocqoVQU/lU I6bojtMr6sTBKZXjUritCX+YWJp9li0OSnTvLLf+JrzrytqPozQJMDuctmtI ugV7pcUtZVlIgAXXZrm/Yxse29VlI7RDApaf93bY4tGGg5YoxRtaSUD7iPsr /sA2HOFTqBC5SwKihNfXu9xrw49Ma9bAPgkQDek5eq2nDY+0LR7K9JaAT+B0 YtehdnxazKvYNVICHMdkTDiOHXim8ul0TLQEvHEs1OJz78ChAVMbcKwEpIQM bJcN6MCxrKvVUolkvPhL0wp3OnDxo3TWp3QJ6Dbs1Ltc34EntrUKaBVJQK5F yL77mzvxuUtmB0abJcDlYsBzZaUubGBpN23ZJgFWEv3DGbpdmCvhdienQwIM RBPLRVAXdr0fUXeAJQE5wr7zz/Z1YbuPbetrf0rAqouv/2jFdWHtZcflH8xI wPMwf6FqajfuS0xiOyvQwerQe03fr904/UjB+ddKdFjNY3vwaVc33qf5RkpG lQ7p6Npo989u3Fj2w7pBgw4LW5N45ha7MWbpvkH6dPBbVGQ76vTgJJWP6bLm dIi0Pl995UYPtn7y70irFx1aswvvHl7NwmcN90QKHqfD4ZdX9kWasnB6+Yvy dSfp8NnrWNOjzSz891OwVJYfHfIou98/28PC8b38dSeCyXh6jo3ccyzMkpA1 FoykQ9Kv/XEKFSzMf/eCy7poOrQ/zp0yfM/Cuit6Y07E0mH9YZNY0zoWDtXK +d2SQIcPnebbRXtYWG2zblZmOh3E/d1+smdZ+PiZzSLriuiw9KJik7whGyfM 5pieKCZdU9JDrGfjN5cEPTJL6HDKqiF1ahMbS9z+UilQRod3e7JVwm3ZuCTH 8VxLJR3uPzI4VOfFxgtt3n3HW+gwkWMZ25XMxisPfhXPbKND1fbFrQ8y2dim Tw+1dNDBNMBf/1geG2dyJ+6asujg/bzsdsYzNjYXvGQt8JMOMoxlAWur2fiW yZ1XGVN00MpUlWjjsvFgS+UG3lk6dG2s7G0dJ68/9fed6zwdHg68XPJmhrx/ rvVnLRoBv19TiJU0Dj4uKdT9fCkBX1fH2P+3lINrnhgflhYhwOyoE6/tMg5W 2+HWFyhGgI+N/BpPUQ5mhVYNI0kCli6AQCCdg624F+ar5QkY/nOtX4TJwfk3 Cy79p0jAh93vGG2yHLxkZTdfpDIBVx5tS4+V5+DXB9aK2KkTUMnuyf+lyMEy 024xJZoEfAl7o3hZmYP94+IZy1cR0PT+To+YKgfrfP4nz9IjgJWIRQXVOTjC VSkLDAg4Pv/zvK8GBw9QbVdmGRKwLDZfulGTg9ONH+u4mxIQPpL91WsVB4vl uMPfzQTQi4SyxVZzsLdZQpXdVgLyBu/yy+pz8Kfu99tKLQg4c1viNNOAg0Po yrZBVgR8LnaznFjDwd2PbVtYNuT9+aQfNxhysLFlyF6z3QTAob18GUYc/Cek 58gSRwLmPxoGMteS+ZAV7nd3IuBHds6tD6TzS008P+8nQFvDOsTFhMzHbo+R VYcIkFsrs3OKtPNIwumoIwTMlP8ZuWBK5if8w8TfowTwJ+PD06RlVMcD7V0J uFnhk+K2jszPG+XFUndyfwL9adWkG/ftuizjRYCwi8TRFes5WHsyhD/Yh4C4 l9297qRvxBTeZJ8gIPe3tlQO6f5VLNGNp8l4A3/nW0lvrBaOu+9HwNVoiehZ 0qnHTKX5Awjo/Brwmb6Bg2cWPZI9zpH9cUM4bwXpPfcSFWrPE/DwwSvF/7vI 8ON97QsEPKV760uQFmkcV48OIWDlE6H2KTKep49Kwb9QAna5hgg1kf4gsFtv z1UCzsqVvk8nrXT/0vMX18n8lifyOpO+gJ6sZd4goFpL9B2DdGcnqyI4goBg W97p1+TzGvov28i5TcBGSef7TqRjxdd92BhNgNNBCfyLzNd+z1ddQbEERIXL TfqRVqs0/lNyh4C3kl1642S+uczSJX8SCFjbreblSfqV3xpZrSQCtpR0pzaT 9bpS90zPJZmA1cHtn/RJ71RbvTUtlYAdQhN91405mNOmfZqeRYBRN993EbL+ +boF13dmE5DRVFhuRvaHX7hWyvVcAjwnBgM8yP7hN1X/NPuQAHPlK6EJZL/V x+b0rHlMwBj7bnUK2Y9Jwyr/TjwhQG/Vk95EPbJ+qYryfcUE3GDt9PPR4eCp iTR9+VKy38Yj+s21OfidtbyF40sCSqnrlkuS/W1PZfrVVhAgc2RvYSTZ/0HH iJriD+Q8LCRrzqtwsHlFDHvkEwFq193+CyLnS4whPq7+mYCj9SF1o+T8ZX1c ppDylYBAAdGeYnI+azT5/a+2k/Vpyq54LsXBcZev3sSd5Lwr/d3SIcnBB7t5 Mqa7CVCZ9drxh5z/P7cotT69BMgXOVQukOeD9J9pxT2D5LxPyvy4zc/B3yzO GkUNk/N2kD1oz8fBjzIndtSMEPBc5slOcR4ORnb/Ajb8IyDk92f9mjk2di0d rlObI+ef7+BTyi821hXz+n54gQCvh5XZTd/ZeMb911QSRRKSIxfFI1nk+SfT ryzKJwk6rHZd3MTGxcHsc5PLJKFSViOdp4KNaZsbVT+tkAQLlzk+xk02Ph/g JzKpKAla8yfQ7StsPJbPmFRVkYSr9w7tGw1m45/i+z5dViddGcj0PsnGtezv nkhPEnpypb5Z2bHxnfPjhSUbJWF+JrO4SIqNRQoTEn9sJu9nvtHnuSgbh31b e4nYKgkXp4JupvOzcaDFhV2nt0tCg0tEh84kCx+U4h9fZScJHKsbvXotLKz+ TNok20USXsW+rC2PYOGM/ldKzW7k+puBhhsvs7CMzIGlPJ6SgLdNhBWS77dl lzK6Dx+XhPihH5/MXMn3oZXmRdkASXj/UrZz0wYWLhs0rYoNkwSzp2k+pb96 sJXSIavQh5Lg0FJcoqrdg2c8625xCiSBIcRuva7Qg3OemdZteCIJjgc0z9ZK 9ODFzdI7ZosloSh6B2N+ohs/cWuw8HstCXox7Q9mKroxvWDTFtcGSagb6/A/ aN6N243U11lOkv/bd46EQhe+EhJ/Pm9aEpzL0vsVdbqwbjVv2ZI5SShvrFbI lOvC4U69ayspDFi/1dvAfKYTmwQnGa0TYgA7d+k7ZlEnTnknrK8txwCm5UX8 j9mJnXf+0ZDYxADFlTHh4u3tOGh34czWLQzoFvGWXV3VjuP3+nwO3sqAEbVt TYZP2nGNyy/vX9sZYFqdGTF8vR2vvtBb+MaOAToLKtsFDNsxtbBhzQlXBpxz 8++djmjD6WJFZrXhDJBaKVUVodqKXzFOStAiGLBS+YWakkgrbpbV/m50mwEe jib5SRMtWEAj/0pWDAMyrS/8M/jYgk+aZX4MvMeA7MStj5a5tWB0OmanRgED fjVOpW/LaMasptN7r9UzoKct7efDxUasNToXMtTIgJNiOQY23xrxWaHreTYt DDD2HHdvqGrEEmb3JmU6GZAva3TQM7wRmz+uvFP4nQFfy8qcxsUacWEY0dg+ wQDJDt8E1+UNOGR9iaWWnBRkXWRNd8l/wXWOcDpqhRRo3I8OMOyqwzJ+NXfH FaWAiW10tyXU4eJ81gBWkwIsEkqJFqnDA1L81+10peBmzL3O59Ofse1fh8rg TVJASC94m7yrxkq506b1HlJw9rzC+1Tmeyyfkmuo6S0FfpY3VynUV2FmrL3e leNSIPglbMLxShUmLj5VM/aVAhUj4/3coUq8xMFDPD1ICiieVnbFpe/w4JKO /hORUnBJ7nr/RbM3uMj1RbTocykI04TlJqYleJNawE8rijSs72h9+eCBMXrc 4ecbQpOG3rcZwQGiHkj61unFp7zScJm3HlveCkDD/3ykJAWl4bZDs1JP+zUU 9/bY1i5xadhy71IRX2I86t+/K9dNWRqKdv5ZEZWfg8JiV7mFmEuDy5N+8zSZ UlRH6et/GiEN+1SErJ8qf0Ke1QXVEpHS8M3k+82WO5+QQEzAI99oadD3evty p0A12qi89PSaeGkY1Nq21mSkGr3Ysnr+RZo0aL9MnuN7/Rll3gihvymShumd VusDfb4if0IOvrRLw8k+mQC9a01ohfqeu4Mqy4GbHHffOrcLdeSXpvE/Ww5x OnaHrzr1o1Mdoif2bpCBi1LmJ/X/jaJ6O58/IS0y8LDDcU3BwDhK1Hr1McKZ CdyeUqNkjTnU/ssvjubChBRVxpOVW+eQdJ7OkbNuTCimSCa9PTaH7qpkzxz1 ZoK1+YovCulzKEkuapWpPxOWf4/mvSU1j+6Jusb+CmdCDY+V7ArBBZT6T/yQ +VMm3P7JCF0/uYhYT2u1yp8xodYna7vVUgrIn7w+pVfCBNUlpaZ8TAqkDc9F y5UxwSuo3zbchALpP35WjVUxwcgjKHjwHAUy2io0s9pJx+kLScxQ4Nuds5PS XUz4yuWPMxaigpKdftXtHibwxRfcdZWhQmZ93oHAb0wYUaHJd5iQro6Nsh1m At7QrycdRIVv160OfBxhQkCN09vACHK9uaDm+j9MmPL93NaXTIWsdxcq1SeY MHnDMO1rBRX6QtZGpU6R6wtmd9l8oYIKGttPzDKhx3FXbDuLCsfmH2vcmGcC 9WfaEfdRKtwv85hYXCSft8gCLy5S4X/t7FqM "]]}, {RGBColor[ NCache[ Rational[2, 3], 0.6666666666666666], 0, NCache[ Rational[1, 3], 0.3333333333333333]], AbsoluteThickness[1.6], Opacity[ 1.], LineBox[CompressedData[" 1:eJwV13k4VO8XAHD7lm3u2JIKKSWVpEJyXlEoX0pFikjZ10QSlUQJJSlKZcua fR2JXlFIthZFQiSKMNc2dr/7+2uez/PMvHOec88597xyth6mdhxsbGyF7Gxs //9MOpHFCFW+oS1xXjlweXkMUuYkgp7HnIEZ1ueO0r4x+C+ce9OVmPMQcSAq P+DdGCR/DlhwirkOITo5DTvCx6DRV93eJiYKyrsEHEwkx6DXUOLJqZgk2Kd4 ekRDdRT+qYbstIgpgD3Y89F7r39wPr2C40TMG3jGk5KhveYvbN2a87jw6CcQ mJe7wp/cB5sx6XlVsxOcMu3fhf3+DvYBEZ/srTpBLOPNw1VV3yFZQ1HNJLAT zhw6/ett3HeQKrZkra3thF12sb02h78DT2bdlTdHfgBDKbW9sqIDft5/Fsbl 1AUNCrQNvHHtEGNvmBIW2wMJhNbvQvuvwH6r1eXQXC9EVjlsMdL/BOqRS6aa Qn2gTdN+umn9J3CPVdbcJNsH4tYRx+Q5P8GP9FA+3gN9oMotfMsOf4TSepRa HdUHCmw0TVWNj+AiUNCtsekXXG8/v6Jneyu03Yk6svFEPwQkso3kuTdC5gPT 3dylA+Bxu/h2uvw7MApdysivG4Dn536eNhx7C0z/F9KWHQOw7dGp75IVb0H9 LMdS4cIAGLtseKxp9hZqtxe+tdEdBO+LJ/liImqgv5lmWtE6CKm1SfCEvRpk +T+6Xhj+A+kcJYaKK1/D24WA7jWLf2BLPlHnXFkJjsyNhxuE/8LIXp8eLttK yPsWuENO9S8QVeLhWtkVoJ2mMtd86S+411t5MfRewSnde7eUeIZg/a621PqQ Moi5YpLcKzsMA0Ovat/oF8EnacsGXrVhODywIPR3rBBEyhzHt+gPg3e60E+v R4UQOn5932W3YdB5OViaM1QA/g7FfUT5MBxvvrqtMjofbI6sXKd39B9ElAY8 MmNmgdL6X8/TQkag84O1Ix89Geyrxz40PhoBUf+6M5H+SfDcemFiPGsEbB7b PvzSnwgyceJ68GkE6iE8YfuFBBARMej/tmYUeHwD5bnb42ByJluBv2wUzJvW 3zX7EwW48WKqy9AY8CzNtW4TMIVLkGyvvjgGRVMBs/9uy4NUx0A0Bx8TruVv tdR4fhC3HiKOt0gzQXj6fV2Ukj3W2e70zRExYVOOM6/g8ctYfkHix5NwJgT3 KGstDYTj31Fev9nlSfB/u7hXPfsJ3neI74DjZhK0Hn18bRPwFMdzxac1q5GQ SX/ySu3QM2zm+97hiT4Jz7FZ/c/+eFxrKTuk5kZ9/26RnhJXEs5QbB51KCOh rrS5ioMvBbtUbpppMh4H+f3TnIX3M/FV16Vn983HwWhawenwz0x8b9VnXXOb cTAU8jay2fICl/oFRP48Pw70rtmdbe9eYLZdLRsmHoxDq9gRFSlmFo7J9T4u 1TkOAh4xCYpKuTjTynD+x69xUHkekUa65+IKwTVJSf/G4ddiRdO2olzc61z3 T2lpHFbHmq7n18zDyorSwXvlJmD1zmsH42TzcXUCLrR1mADHuz8z30sX4Dbj BycUPSfgoYdQZId6Af6z6Lg0fGkCzuida0sxK8BClsRBn7AJGGf73OEWVYBP SNn13sqZgHjjcJFGzkI8ck9ANGdiAiYCVz370FmI2VBPyfmFCWjbWu/9ZKoQ E2NFp3ZxT0JuftaSokgRVv/PKh1LTIK748krvjpFOIg/X/uTxiSkqGzxHk8u wlJBx91Y1yZh2NGy1NmiGLNCTrr73Z6EAs3NMz2exfjrbWuP+fuT4FbEnsob WowfRjmdX06dhKeXd2a5lxRjIumKD2/jJPiE1FnEC5ZgMiXo4u22SUgu7vuV KVeCWzNu+a7omQQHo995jrtK8N38KD+R8UmIiBH1Y7MuwSuqUq9ISk3BbLlw 1ZPsEjxU8+LqY7kpSD7mbaiHS3BDXd61VZun4JZU0nDixxIc2vLy+lrtKWhe xr8vTpVgrp6mEMVzU7DaSz7pkUYpXlqcurM7fwpu0ieKG4pKcRf7/N2XL6fA NnPwzJs3pbiCm+3enpopeP3Pts2tpRRfFhK4D1+noM7Syv/P31LMklkTo78w BW2854scVjHwV9l1se+5p8F524tgEUUGLlXY+OiQyDTI1xpOuKkysLeyapyJ /DTUySmLmRgwMHPPgXhzg2koYEQcS/Vk4BY4lNBxZBqIq3zs3JcZOFf3cOKp U9Og7tMgJ3ODgd0OnUy2dp+G9K/3Np9/wMBGJtbPe32nwWz7fwfznjKw8tFz KWevT8PTZ9ggNYWBh066pzk8mIaK8F2HXxYx8PvTF9L/PpuGpcXxu9/KGTjD 9lKGS/o0ZNi2SeW+YWB756AXHuXTwN3v7xfczMD73W9lkTVU/CQrJPgLAyt4 RWRfaJqGix/5Dfd+Z+BffjG5vj+nIe7f6dMd/Qx86cL3SfO/09DK8YTp8peB hd3W7FEfn4aWz7Kvq0YYOMXe9rrU/DRcr2748Y1kYE2btLoZThZ82X7BtmCK gVsthoQ6BFlQEDx57r9ZBrY7uvXYS3EW8K7ZNJe1wMD3DpT+9FNkwYc8cbUC jjK8Ac1tOKnCArp7ENcJ7jJcoaHtpqnBgpClFfgNbxk23RFUJL2PBYJvI/xG +cvwoHLt7NxBFkTL9Rh2ryjDARsEUOdRFghc7tofJVSGCVnjm68sWXBSxfYG r0gZzlh5v/GJHQs80/WldETLsDb9KxHgzgL1fY7S2rQy/EVQ2sLSlwVJZ+MT lig785xO0ApkwQP3t5VBRBlmY0v+LXObBUvvisMaKD+c/b15MYoFjY2GMl8p K01s8uqKY8HQb6vQVMpV/9zKKp+zQM/7Q786ZbOBgqVn2Sx4991V/x51/r+e Kb2rJSywDJCpzaLiCerQCD/9mgVTr4ovhlPxSn2+8lG7jgWXFkQ8tgmX4dzG N5JrW1nwNZHvXaxgGdat5T693M6C2RV+d98KlOF2bJjS08uClsQ9vSV8Zdj9 5Z0hPMSCfdVKn114yjBX0UeVxAkWoM+ylwc4y3Bctrhv4AILSnbNDK9nL8Mq aRaVNtwzEGYSrq20xMC1Cc84dYRnIC+w+NbEHAOT99ffY5edgTiG6vG2CQa+ GeH0tXfjDOy8ocU9OsbAMjdzZKq3z8Cn8f2cjcMMbOC3MzNIdwbKDd+g9j4G 7vbyG7M1moGe6cEsnm4GvuBauVP3+AyUfvatmGtn4Hjr/dWcDjOgmKG1Xo2q 16n9xzuDw2Zg7M5ev9ESBg6Dx3J20TPgHXbvglEeA8tqdDnsfzoD/M1rqz0z qP5Rtpvkzp0B2VLLIsE4Bk4mfIRufZyBD52qBx9Q/acu+PKow/cZqOw7eFDR i4GbuRcf6/+aAXkDr/+uOzHw7EzwBr6pGWArP1EcZs7Ah3sewG2pWWhc4eA4 pULVu67jJy25WbDhq2uY3sDArPQ9dsxNsxDGqSRaLUPF79kXZr5nFn5lHI0p 5WXgf+wqX9efngXrmPzWlZ2lGNlzOnXYz4Kd+/NwV2r+RDd8nY/wmIXrWcbN t2tKsWb0FbnJwFnYy7dpr0JWKb6t8MG1+vkseKk6snf6luIfofHLPtmzMLjy sWawcylWGTl/f1PJLIw27k6csSzF7SVSZfdqZ2HRYKWurk4p3qBvz2n9dxba rw0q2vKW4mon9rj5rXPgtCb0fPvtEizR/EU5b/ccJFftktL2K8FOqhnYFs2B z3ylmqdjCRadN/79/sgc0Ebz5TQPlODTEU9VHnnPAZnzfIXDcjGey91dp1Y+ B/rN0et2OxbjHZNuk2775yEenQhCYkX4WM75vm7jeRBpiHNynC/E3vY+rSYn 5kH3cW+OWV8hLm4PyN7uMg9pijmlN/MKsVpl+Lmpe/Pg3BmduVuf8s3MLwE/ 5uGk+WnJSZcCrCb1uzj8wgIw6suPrJHMw8c+/kleCFiACeFBA5dvudg77N89 t5sLsK9xy51Dsbm4eGHC1eTxAgT8dW2pFc/Far0c6+l4AWItz8znCOVgtRey D+MEFmFXTOK0698XeKeWpU9m0iKs7z6btto0DQsOx2dcerEIrPU3/WwF0nDf 495O/aJFuDEcbm1cnYrvsez3DbxdhNeFm0VktqfiocLzIuv+LIJyqGxmtEAK jt90K/Pp1iVgX2bfvDInCfNKFnZFvloCgAvLh88+xh0k74GLbcswfTTgJH2d Kw6YuvyFo2cZfA91+3H+sMOysyO2kX+WoaLuaKPTGmtsz/blWsb8Mghdk+Q8 uG439uNttRwRZUMTnTb1oSX2wC3xNGbtbjaUHavkz6oMgtU71ASCb7ChnMXx nJaJGDB2O0caybCjKJ3YO0dWvABnjcIjBgYcqCnLOoEmUgnmFd0WgjacKO1x S83J37XwpV283j+SC9nEKB93t24F+ezN3UUvuFGK7I9wf+02WJYyOhnbxYPO rJwdTCv7DrFdFo79AnzoWI7EPueVP4HN4sCauXX86L4fqP6X3gczbns1Xx0T QOK/eNwDrH/DGw36kt75FUjhc+HswMdB2G/19ah6qCCytzt0Q0RyCCykU0S+ lQqhvM07TwZb/gPtu4cvaH8URpoTETtrdUbB29n3Nc9XYZRSmnfmguEovDgQ z9/yXRi1qvJmrToyCuLLwwnW/cKodmvGLSsb6v7lcasxkCWMnm5Y5Eu7OgpP jrxeX7NaBJ1IOrhVv3wUZsSU2/WdRZD3GhWJwi1jsJU0XSfqIYLkdG4ZVqiN wbkmP/f2CyLo+GCB8ps9Y9ASUsfldEUE5ce2TZYbjkHatO32sHsi6JFfkrK7 3Rgcb38c1lgqgtKakrVan4xBR9ON+Ax2UcTFH3xpJwcTUiT2xLJziaLQmnqP o9xMcLcejzzJI4oYsb+S3am9m4O0uS64QhR5PVxX+lSYCcpi2uc8xUTRM9mQ s9+oPfzayZlN6oqiqFP/nsc7VSasH3QtqTMSRX3ObXOXbJjA3KaQK2siinaf fXNwz1kmlF/qTPM7IopaGDXJc3ZMMBY4+GizmShafTX6nKcLEy5tUfSPtBZF AS8lOdQvMuHDhV5k7iWKLCq9NXdRe/15NrPGgVhRFBgWafyrkAkveKzMJONE UdBEXI50CRP6BM/91H8qilyfJRw0YTDBdKXXZEaiKDJZOKaR/4oJqqp3ZVwy RVH0z4K3h98yYfxsretYuSgSEddzW9HGBK+6XUIzXaJooNLgq9UkdX7T3tiN P0XR7KEEJ81p6vzPenIWfVS+rEhRsRnq/J+mO8sHRJHbUYOY1/PU+bPuVgFj oigymNOcyUHCuFJ6Dhs7Dek1vw57J0qC111JY34FGrKPrAvfokRCV+bQVqcN NJTYeJ9ooO4lBu8qRd5vpCFp1Y0ztltIWLNw9mPoFhq6VbncfluFhA9O+cf4 dtOQaWNIddYuEhT2G57iPUhD035Cr7X3kRBpI6PlYERDLcnpjom6JMz5j8nU GdNQA7HffEmPhNbCh903j9JQp0j7bC51z7ki23eGx4qGrngHrf9lRMLX+cuO 3J40JPSIL/GhGQk6ksaGdl409CiOs7rcnIRsVTmld940pMi5R+jHCRKCnOqG g/1oiNs4R4p+ioRt3+geXDdoqLmq5oupNQlx4wMm50Jo6PLfdzQLGxJ4hMtV 3t6ioTjeJOdTZ0j4oWczfiOChkZcVJyOnSUhtDDbhzOGhv40fHOWdCBhouma 2dlHNKT58J35EuXTf01318TR0Pn+jUd7HUlQk52dCUqgIV2LePcnziT03tkf wJFJQ8H3Hfl/Uvc2o0wpK9ssGmK/vGFPujsJZW+H91bn0NAXn2k3Zw8S7s5H LV8vpKFt49Kt/Z4kaDp1B7FX0ND7icfmmRdIoO9jn1J/TUMJRl/t93uT8E9a wcGzisqXpexYF+VnjU6Het7SUE+iix3XRRLYVKbolU00NF9WVb3hEgkdfFI3 J1toqKvWyDuNcmGv5szmT5TDOB7I+5FwNjrwR9xXGsJqflKil0monV6ReqmH hj4J+ohU+ZOQ0LJVMr+XhhzqVQTXBZBwKePI7cFfVLwaqq5BlDedjHUz+0ND PzQeHN51hQSOHeU9d4aoenm/51MY5e8ruo68+0dDTCFGQSfliEr5XWokDQUW lwZ7XSXBLmZ/hssEDblJBPmUUdb2cJR+PkXVo3Le6znKEgbhEd9ZNLR4sfe0 xjUSxmRzl2hzNGTnU3HMm3L9bKun4QINVaT/is+inPhpoi9wiYb678zt6aFs GqxRN8ZBoB0+IU57AknYbGWpochNoNUXEessZa5d17JO8xLonJLfp1DKP4ST V8fwE+gP7fNsJuWSwbeRTSsIdCebdKmjfLdqkJ1bmEA7L19S7KVs/1jAW0uU QDVzAptZlMFry8AFgkDSuUY+/NdJkDp0+ESWGGWxXi4pyuS6Cw19EgR6ke3S Lk+5YeGhlvRKAmWyIkc3UU5uK8s9sopA897k/i2U/XM7ZW+vJlBQyo4uZcrH bi3dr1pLIDmeucKNlJVt5Lhn5Aj0auWfd7KUuTX0fLcpEEjVOklMnHI3zeGv /QYCOZm9f8ZNmTF0+1T8RgLpXOO0GafivVeT3dSmRKAEiwmLTsqOT1tAaAuB 7Oy0wqoo6/iMF+htI5CPfNBYEmVpY3GFgO0Eyhd1DL1GeXyDekzRDgIt5EUd taD8Yfkk3/BOAmkEJptso5zSfuWyvDqBjNi1/dkpXylI/GehSSDxIt62Fur5 mIXVnI7SIpCv+PNTcZS3nR1ordcmkHlGkdAZynxa/LpsOgR6MPzj7zrKvWLK Jbt1CUS3rh7to+qjfMRY0WM/gRhBvKviKbskPFjRbUigyQ7tfj7Kk6adhwWM CPQo7NTqcqr+AnjkY3YZE+idLX+Cw//r0S1PNtKUQHMGrIsMqp7F5KbtXh0j 0LPStAYLyk+/aGUNmhGIfUXtmRmqH7L3NOxEpwhkrLjjuBJltTFRf1crAn00 ZuZWUv1TkWxe9ciaQCYyCwf+o9zE//sgeZZAXo7dymeofhv9tmyd7Eqg/uX+ rFO+JFwMP5Da7E6g2qitX1qp/l3SvjM050kgwiBYZR9lkTRpn6M+BFLoVQ1c 6UOCivfOMK6rBMpbLjKP86LmiWJAi0oggdzHfBaY50lAndViVkEEEghi/6BL +ci+wwklNwnEY3LsZxc1b7xEXYrtIwmUUb1c98uVms81BbP3o6h8cG6IXUU5 yHdGG0cTaIv95xATFxKiu0PeSz6i8m2cXJHlREJxdmJ3fSKB6vsVXXfak6Bl M7huKplA1c7xdGM7Et7StzrJpRLo4j2LVttzJLRdfjXhl0mg/ZcDggNtSZg2 +Mq3uZBA/BwRRRGnSbi6IGNsXkyg01XavSFW1HzOPxt9o5RA3U8nVvpbkiAp Sa7+UU49f4996RYnSVD/vWLHnRoCnS0NFxo8ToLfdR2rsS8Ean12868O9f5Q O3hs9uA3AuH2h9t6DpLAJBwepnUQaK3ns2JfQ6pfUyKarLoJtK2N1yfuANVP dd/2Ng4S6Li1Y1IaImGrkPuazDmqX1+YaLTvIGGoLbCca5FABrf1XmxUJSE1 PtrMZplAjQZz4E29D2VUXt6V5KKj+tvSnxeUSeA/yrUcIkRHNnNGzp/Xk9D/ KK7HVpaOXBo9Rx5IUPPpTI7/a3k6mr/z+XeMGAmnlKokpdfTUU0R/hNNkPDp 1W+Tj5voaNQojfeaMAm4W6UKdtBRwNMP/it5qPeZQl2izAE6ahtwa/7/fmCS P3HmqwsdOXU2pEg1UfvLLrNIfnc6Ggn2+JfawITEirIKLU86umjvWbK1nto3 6q9IPvemI1kBtVrVGibE9PI2eVyho3Lra5VcL5nQTcio80dS/58s0NH8nNrH fPSEtQrpqDpQ4sK4DxNi59P2eBTT0f7ITaYjXkyous7vlFxKRxpxHKG/PJhA 3G2u4XtFRzNdReNvnJhQmnbCr62Gjt5ktMQfsGLC0jfXfvc2Kj8vWhXddZlw R/NhedIMHekcX5WtLMik8l2jzTVPRz4tnzt0qP3vwPnxavtFOqKJrD51hIv6 fbrJh80cYojPf8fz1zNj4C4u8KNkhRgS6w2pdGsfA2Pm1cX3a8SQMim9kPpw DETTHNG4nhgqeXm/3553DFx1Yt8e0xdD0RLH769dGoX6H+8MGIZiKG+09HLj 5CgE0tcdCTAWQ1nN9SMcfaNABnad4TkhhsQDdyz1vhqFT6dMb0i7iKHGPTGt dm6jEE3Tqt0XJYYU8lZfc60bAUvn8s6AaDH0xVeczapyBDbUqJOlD8XQjVPK PfuKRqDce6fM5jgxJGObnT0aPwI/v231oj8XQ4Jj6btrfUZga7zcmv5iMcS1 Vu65l9wINCjxXgxpF0PyURUBOZ7/gEPv0/r6teLoZvizBNHRIfD39RZmyYkj nvaZAzd6h2DyhQRrvYI44t5iJjj8ZQgGaafqb2wUR5K/cgcevxqCxp5fzrBd HInLs0V5hw7BQ/+pvNJ94shm6bTPKtkh2FgkpZlqJ45s66cPrNb/C8by1sZB WeJo57nPltYXBmHOuenOzxxxNMbL/oRhOwhpRXuatPPFka7vSn0O00FY1pMy mi8WR4HvqzvdVQYh3+GjofdrcdRVsPw45d8A0HN099t/FEfX9toHXbMdgPbd G7UOssTR66ec7zp1f4Ptf+QmQlcCRVk6y3v/6IPuz14WN1sl0MYvwv7joZ0g nz67p9VJEplF7x3+o9MEuht8B43ZpNBVLwnc4IlwE1v/QEGEFJJvPmZevaEF r91o9nhIYSUSy7Dq1TXoxR0vGAm8RSvROg4OweTtw/h8h4iHhbY04pK54Z3J Gsetx9zIwDZpVPb4bITCuVlcLwlLHSqrkP7G8rCwtEX8P1nWX2Y= "]]}, {RGBColor[1, 0, 0], AbsoluteThickness[1.6], Opacity[1.], LineBox[CompressedData[" 1:eJwVlXk8lF8bxo19xi6kRpQIJUlKKPexJJVs8Strqy1LyJ6QJbLNKBEK2Y1d WtEJIUUUoYWQJVHZZ2x55/3r+Xw/zznPfe7rXPf1bLtwxcyelYWFpYnAwvL/ 58MzxU+jFcO1yjvUlVJTW3HuslhYTvJ5COyrKWa50YpPxnIoXE/2BLbUM2U+ Tq04uyto1Tn5Bpw/Jn0gWq0Vt/kddDiXnAiVmU8WjXre4KFjYunWyQ9hYnTE 85zYGzylErnfMrkS/szas1zPasaeBbWsZ5LrIWbw8r+Jr41YSak0terUR5Bq sbd8L1iHvRCp9aTfR/D3EzJ2nKrFT8wcGT/TPsJ/QaG3xd/UYi2/baelhj9C 7dH9PxtDarHRq7sicR5dUKbdJXp4uga7mQVT7OO7wSG/uOLRpxe4xNc4Uryl BxQsZNPDq57hXXjGI1jjK6idWVjQrKrGDkFxHx1sv8LTuY7zJ5Oqcba6nKpx 6Fewmtl4x9u3GotX29Clmr9CVBzbbnbNasxZ1HK93vQbnLz7gSLX9AgP3n4Q w+7cD3E3Yx8E9FfhZIdjuTEp32EyJVZTSrISE6I6XU4sDwHlT6OMXG0xPkj5 Z6bBNwx7D7ife3CzGLunKGoobB0GdtbPirqmxfhbQTQ3l/4wLN00PrhlnIaf vEF5DYnDUD1pbUMUpWEXUuWAusIP0DPs+P0ooBB/ik80lT8zAqdaif3D9nmY Nxkf3OgyAmz+dDvFA3lYJ+O3FGfwCOyuzPLO4czDFWXH/vzIGYGs7ZHkzYW5 OKaDEJv5dwSKyvamyf3OwSDk9VosahQ23lq88jgkGxclmalxPBmDGdu9BlZx mdgw+l9hRcsY9I55p1YezsTT12ibbT6PQYR+yRfzvxn44EXWf1WrY7BrILIp +FQGbt5b9fqc7jgov+hwspV6gEfeC5nVdo7DdKj4kT1NaXgr8YPr1cmf4Klg auxhkIxfrwYNSK79hDqVoIgNbMnYaVre5C3/BECVqPt63V1c3hu6b5vKBAw4 bqdXqN7FWvnKy+/9J2CakcTrL5uErXWpUTs5f4EX/67cEJHbOPm6cfbQ1kno bJ8khBjH4Y+bbd5yqU6C08WeaO+2WCzwzGl299FJ2OTfrM53PBZHz97QCXSb BP9xB2G5ozH4mmP1sPCLSXht6eeyYBCNz5lu2q53agp0XY5vb7wcgXfK/sjJ j/wNpS4VL4qPB2KHhr/v2u79hpD8PaoXXANwztnVudni36Dj9CTOiOKPJdJE 9eDjb+gX0mhr7vPFAgIGI72Sf6B7/omkXOBVPM8okSE++wNc7v21jlPOGLf5 5rn8+gsQH8Z6c1UB+0O2w8G1v5CwNK7sWbX8Uvzz2B1W7mlY/m9XBEVxO3Se ELbo2DwNORl39A+PAGjvde51QtPwnCXiyplJM5BeFfuWHjsNSedihV3XL8Jo otcoQXoG2htivd/VXAWdE9z6TrtmoLab/WKSkTdksGfkv1edgUvVD4TNBr3h P79Wx/SjM7A+F1FQyuILzTZbf6m6zcBvkQiDOVV/KJR7/8fx2QzA8KbLsq5B 4FKnwGg3moXRMdmH+rxhEOz678Ht07PQ8vul2X3tMKCSu3RPn5uF0DOa4RG+ YfAkIIgy6DkLVP1Nu/i+hwHLgY4dc0mzwM/CKclRHA7JZd4W4l9nwaD+ZnuN UiQU2R5b+fZjFva1Hd5qahsJtbySDx9OzQLb38bNQbGRMHS5ZWrnv1kwd/b/ cGcsEhTlNkcc3jYHw8WS7R0pN6EhE1ddcJyD3dp6GvljUfDJKOmMnMcctFK7 a0v5o+HnmtO/Sf85OF0iJWh0IBr4bISP+8TMgbbqtjnLiGg4I24/FFU6ByGs 7iq55Fvwm0oSLJ2bA1lP4S4ltRhgQd8fe67OwRZ9NZGvZ2JA+O8j6wMc88CT zmq8PTAGDp60LcBi8zBW2XItuiYGwogVWh/V5yGAv1C9WD0WxMMs3Ogh8/DB Oi88eHsc0COt3ANuzcPqQ/kuY+046Ll19srK7XloVGLXKbeLg7uJzp7refPM vh+EbLoXB8IPr/twtc1DadKxqGyOeJjJDfO99WkeVg6EDkxKxUNnYZQfz/d5 cKyrP9OkHg8JFYkBArPzcKFJ65mCazzwvMq7vlF8AdQmhJTxu3j41UgLTt22 AL1d1BLqcDy8bSkPIe9aACnUbjzPiIfojuc3pLQWwHbqNO9VmQRg/94eKXdp AdZf5F228E2AH8Mfbxa5LcC2L4Sd+FYCNIz1Ru3yWwAP8fcb6u4nwI0/Q7f2 xCyAuTbVybc+Af6tLcSrVSxAnEvJ9jEOCvQTVhKeP1+AS8WBugxRCtRysFA1 GxdgqsAlMUuWAoF8pNvQswCfZCnPq3UpYCkkcKf++wKkGZWJbDOjwEFRkSTd iQXoDGvNlDlHAbqEZPLR1QWwt2k/OhtIgZ6t21NaORZBuih7tDaKAk9k5O+d EFiEjSYnH0slUcBbUSXNWHoRXFZWJsuKKTCtqZ9x2mARxP6WNy5/okAHnMj8 bLoIP1W3rcV/p0CZrkmWtfUiBHhtdq35SQG3E1bZZ90XIX178NToEgUMjc/m DPktgl7mZTzDQgXFU5dyL95YhH8d+jIv2anwy8o93zFpEZIa1rRf8VCh1e5q wcSDRdjWZ6S5yE+Fwgv+hS4Fi2By8LDpshAVohyuF/2uWIRlXvbcNhEqOFwO o115sQj0kIFTbhupcMQ9qnimcRG21nDd+LaJCjJecSVX2xfhxM8XJzZLUIHd N7F0oWcRzPfu694hSYUfAcllfoOLIPSqWIN1KxX8r36ZPz2xCLuvdr9P3EYF fjdJzYOzi2ATY9CyKE2FXIcLN8RXFkHqmLbWXhkqaJzLb2Gw0YE+Sz+pLUuF TstffJ956eD4Lk9QbgcV7E8pmT8XpcOMaETWMJOXDb3SUiXpICrWttFbjsrM gSeDAXJ0OCPwlPqVyTvQ8g4rZTqoHwmW3SRPhVp1LTcNdTocWrGaVmSy2b6w R5t16DDsfZNnI5PHFZuXlo/ToWRBPbGXuT9oBwl9PUUHyfqydHcmC281ullj Q4dnmyWPfGPWL9x0uy3dng5rR+oLtjJZa0OPcJA7Ha7c/NCvyTx/N+9mSxs/ OlCkoglKzP4uc9plHgqlw4C55G46s38WluxRiVt0GHRuD0xh6nN3aXTXWiId xLPfswsw9ds5p+DVn0YHdm3rUVumvq+m3J7V5dBBLqbkQART/6nvC3rBj5nn UwhLOi1OhbDP6rF2L+mwwBHwiV2MCuJd1z9otdCZ+bFp+tYGKug2c9it99Fh deUQWYLphz58LPf7EPM8RirnFJh+cX8e/wv/YtZzejLOz02FtBJRv9BVOuxe MZa7SKCCcr5l3TkOBlxmlMvlrlCgOfMBmzY/A0RzCk7vm6fAzG1ZKmErA8z4 xRMcRihwM865Z0ieAZPffNZI3yggcbNUomEvA/Z6RN4930UBg4D9RWG6DOhg iRDaUU+BAa+AvxcMGaDEe8Yu6SkFrrrW7de1YMAJRkHe41IKZJw90sDmyIDm hMgFoVQKLByx+BoRwwD5gJf/qblSIAZSt9nfYUC7EWOfLnNet6r3Ox65zwDD wzrlbObM+VG0n+coY8BoaKdpjgYFsoV9+KI+MGDBwZbtIjsFTL4nwS3xJWj5 4L4qFJMAubpOHw9tW4Kni0Le3QEJQC/QtJ9WWAIrZd5MU+cEyPAYjjmtuQSV NiXdQUcTYIqg3CNrtwTZi/I77rMkwC2Zd64NOUtguzE4Q8EhHr5FZ6z7lCzB kRkRnzLTeFD+7Xlb4fESfHl3X3TmUDz0PRZ/Rm1eAtI/2cAK4XjYcdSB7ezE EsQ8Kj64pSYOGpwJaStKy/D6+N+Wz+uxIPa+W7FcbRlcUsUdzcZjwVmlEF9A y4C728Zi38eC4IrRaKvpMpyPCzqpcz8W7OLuK9/zXoYUm5tSDaqxsFym1qL6 YhnuOZUJq1vFwL55t3m3Iyvg5lYgbucfDealnsMDRivQeMUiNtkiGrwdfDqN z6zA0v5RRqxKNFT3BZXsdVkB0xmp+rLJKFCti720QF0BDmvN31+smHyzqDvo 2wooPO76zb77JqiKj1bHXl2FfMs+O8vscDD/8DN7NWgVpub7fIv9wsE7Zorq dnMVXrhfdmswDIfq1TlX49RV0Bp6PAuLYaA6xCq7Aa/C8SW5I6H6TKZtvZtG WgOBCh3dk99CYf8hG5+ih2sgREajuXXXgHcyo9CftgYiO0WSyeevwXDq0Nej j9ZAq6mCQ539GlDpDjpjr9fgdIF3jfOxQPhV5Smw/eca0I2Dn5V2+EOGQlTR faV/8POz2aNNbT7AtbGqn1LzDwTF7rNl8l+BzzNc+r6f1qHkyeNw4m8ERm6X ZgwlCMh0sTxNelswvqxeZWpgwIqC3/OX5opT8enaAUvec2yoU9Jxo8Ghe7i7 T/TNNQo7opREJDedz8LSJbsGHtE40KUg/onXhFzcs1vCNaScAxVr6U06C+bi W+W8y8erOdAH+ko4u1Qunq76LTpcx4E6llsPbD6Ui18+Lz8p+IEDedPbPjb5 5GLLln11rnQOtNel20VgPBdThzXvy+pxImlNMVvGqzysd0lx54wBJ7qlECyj /D4P00clntWe5ERt7rs9rb/mYbuJta5TpznRZJHGmbiFPKw4jXnCLnOiyOTV eAuFfPxmTffaAJUTLf0Ml/JLyMfr4oZWKf2cyMOzXTTcqAALm0QkGA5zIgu9 lS+fLAuwTFRtA2GcE7Ec/Z622b4AGywo7nKd5kTR3UptAdcK8O0PfKuIjQud GCjTii8owLIxHQ9+yXOhPvKxyuXVAnxixWzosDcXsnmeejjsfiG2U4kVnfPn Qlz8vpyLeYXYw7nxWOF1LnTQUvmTZXkhTurdVyUcxYUE23KOzNQX4oFHIuHj qVzo5oxBre1YIfZ07ZFJxFwI07KNpBSLcEq/pdMIiRvVmRz176UV4QjKXmwg wI2ybqqpVlQWYS9toljpBm50tirEIuBZETbKe9boI8GN/v6ne/dbUxFmd9so xaXEjSSaWvJGvzPXr3Z/UjDlRvE1elOLgjR8rqxEMcGCG3FNvOXuFKPhk+ci wmctuVHQi//wPQkaln+9b2/NBW40VE1rYcjR8PfY23GG3tzoFOPq0MxhGjYk m+heSeFG1cBDrnOgYY12ubSudG40rs/645YLDcuFrE+rZXEj1aySgCMeNMz6 oyyDpYi5vlp/MDqAhp/T+JYTX3CjOYtXShBLw/k2oyaLL7lR+8AN7hsUGk7i ryuwauRGPdSTJyrv0LCHl6vF9jZutPZV6fO3dBreofGuorqfG51ZN6iyK6Zh kalsrk3D3OhiQ6yHcDkNEzIC7a6PcaPSyD0q1VU0/I2wk1f/LzeK82+Kq39O w3daox17WYhImsVak7WFhkOvnXupyUFEPx9/qpB5S8Puuw+KZhGJaF2H+7Vq Ow0fTxxvcBImonyjV9QtXTTMYqkvubydiNIlSs4KDtCw+f4n0nXyRNRvK53c PUjDhUI75EJ3E1ES6fbTsB80bPKWQ5lTjYhA9W7g/Z80nJvvo9qqSUQ7/F+v cE3SMCNs9GAcIiL7xzJytr9pOFOzSVv4OBH1FpVrvp6h4dmN+/U/GRHRipwv e88cDevP5x6/d4qIbHfnHfiwQMN/SiJOSdoSUbO1d3L4Eg3r3Jo/PXSeiPCJ ihyNFRpOtr9kk+tARF9N4j5/WaVhLUk9+50eROQsbxPcs07Dt5cfOf/2JqK2 +sm2vYRiPNaz3b0igIg0/5mx+bIWY41Hd7yuBhORUv1ZoRy2YpxAYfM7EE5E W/58H3vGXoyHXa5eW4oiokinh1efchTjAwY/QmrjiOjbBq+CDM5iHCNzKiIk kYgyu6S83LmKcT9LY7ROMhHpxvo37uAuxsr9KvEc6UR0RPJ4bAuTI55nJ77J ZH7f48JzE2Ix7rsrnBybS0T6lyiHG5is6BWWZlRERGoDGVwSpGIcajSbIVRG RC71dlx2TO7eeSGnu4qpz8x9xWgmy3N9LEh5SkQRpnKOaUwO+qFdYlVLRFZN 7/KTmNyJKyu21BPRGeUzX32ZLHN/2+PBJiJa8kmbR0z29098nvOWiJ4EnBue ZdZvMye8dOggIi4p/6Q4Jm/d69mg0M3U81j5Mh+TvfmGmqf6iMjkWxe7H7Of NxMm78r7iUij5lF2C7N/cvOrDq9hIspqlXi9ztTHPVu5e/84ES3/+nZGiskN wVl9jEkiamBtMNnB1FPMWrC/ZpqIPOZT00WYel9WCx0KXiCisBwVhV/M+3m5 YXpUe5mIXOcsJ/KZ9yc8ffYX+zoRlbT0vzvOUoxfFMJcDDcJkWJVnh5do2H+ yHL6ST4SktTOzMxi+uPCealVQWESciFmOQ8y/UPavM6eQiahVLu750UWadhu 0Z1otZWE3Kpkr3DP03DVxwG+LbIkdC7/xpYfTH9axr4Uy1EiofjVLrIh078l jkpkh30k9FTW4cSXX8x50cuQUjhIQhNBK1+MmP4vXL0uX65NQjq6jp9/DTP9 73ZYo8achDYJvYmf6mHmR5Ph+d2WJPSF+ocwwpy3/yRsojNtSUjZL7b8dScN V78N7Al3IKExG/cwbea8eso98zL0J6H8D/V35WppOCu4JfVlEAnteRdxyfEZ DXd+6nmlfIOEfi4cY6NW07BSxAK/aAwJqb/eGpReQsO/BlWK+9NJKOLWxEQq M182HdT5aJzFPL9pQYtkCg0bUEyX6nNJ6InoHE/cbWZeHfY4ml9KQs6+gXek btHwxbTSH+6YhJq+792t4MPMm1PyEiw/SMjjWjxrjQEN8xSr6XqNM/X8bk1+ pcPMS9ajl0cmSailu6Ol5BAN36uwf/Zmnln/9J/8Q8o0bMGXY57IyYNS2P7y 7mfm8/vmLfHSO3mQgdSxv0H9RXhty+7qJCUeNOP1UGThUxFW9Dn0lXMfD+qd erJ86n0Rjt1urTCpyYMKD57tbnpZhA1C7zU9OsmDKD4Jhv4ZRbhefcM/PU8e dL20soXXsgib6F9+l+XDg9oLW+bFTYrwgFn9vdUAHmT5N2QD79EivOxyZV91 GA8ampefz1ctwioZ75y2J/EgE8K9sVL+IpzNGtFFeMqDFsmNlQt1hVhZ4GuW bQ0PuuMbv8PlUSF+SVZxf455UKIfam4uLMRfVAe5Pd/woEi2KnPlO4VY2PGQ 1uBnZv0L/0zfORTi8HfzhS9XedA+TdFpUSLzfd8J380EXtRcsIR+Mv+/WSPZ ur4cvGhSlG8xc7oA166ZDOzm50VF8cLWuLcAz+8p3fBAihfFkLqqa3ML8KUk ++BrOrxI5QBba756AT5i23PqYDQvmlf7GTZklI81kprez8bxIjfT/fLzh/Ox 0rvqY6WJvEj0ZkHBqGI+FlO/A9LpvMi8I27hEikf/xQx2cVbxlzfbCHU1ZiH Y9paWQe7eFGa+zOns8p5uEOzrjJKig/1DgXcbf2bg197lSjqyPChBJvJpntf cvCzovSCVXk+9N6e/5NeUw7O2hj4wFOFD0HVsvG+tBzsNXfglvURPvRpo+Oq o04OFimuOLfHhQ9dOOxgLB2XjS035wr0PuFDfxSffnNhf4jl21RKUQ0fGj/j OeAxlYUZ1+uP0zAfKjc1plh0Z+F7QwOR19/wob1cv97W5WThnsJNq9u/8CFr wQbErZ2FzdQoPz3W+NAmsteruV8Z2NA8sJ6kx4++xht6TMakYwlO4llvA35U W+BrWG6Wjiefpqz2G/KjIRVJoz2b03EM+bFapQU/KjndZnWxMA2/Gf5betqR Hykw5lru1afiI572abkx/Cg+YOKi9HgK1kowuar1gR9556keWF27g70v+73k 7OFHmxhFnRlP72Cafgax4ws/StBxmvztcQeLrk9mnh3hR1MngiUTh2/jqStR baF0fpRF0h8Zr0/E6aYvZRu3CCBz8bEr3kMJmCGi2Hf0sgAKDb2aJi4RiZVm zLYLXhFACVmqL/4FROBL7QHufVcFkGSIR/7+3nDcEdnC7nxdANUmfupPooTh /MULe2OoAigl7Mb+7dMh2KIvNabtiQCyHpl5ddfYB39uD88oJAiiWXq0U9WP PZArpplCYBdECdpdX0eP64L72VmKFacg0kRPrSYLjIB15twNXh5B9MN8/4m7 9VagKKJ1yUNEEJXyeJQ8/+sEIVYMhYNygqiw8oh17zF/kB13fdxiKIiK6DKD 5cHR4MnyX9tYiiCSnIvTDixPBRqn7X8b0wSR9O4ZrSub0mCY99Lg0fuCSK1R 5vvT8DQw2+Q1X5gliE5fqnfVs0wHFZUECZciQSQsoRd8kvsBzF5sdv37QhAJ qXE2btbMAq+WA3yMfkG08cabBs7WHPBK2GhElBFCQrE9ohWOxdBf9EvJeYcQ erucsjpeUgwGTXUCrfJC6FjF02bz2WKQXL34IXq3EPLOtjjSfr0E3jlXmHOr CaGu+6OlmimlIHPkmDXXcSE0XSWrKdFdDj0rgU4cHkIowtLMTkmzCjScB8II tUJISf1N1B7lp+CSmcQzcEwY7dcXvGr38RUE3NC2/dstjEo2hB5Sxc1gXDF3 vsdlA6p5peDAu+09xGvcffGQsQG1Fvb8sVv5CHeEDjXrJIqg/bM5Scebe4FV 76PsGylRxD9Y/uLw229gJH3WKKxYFAX77uJJlxyCCydnFIR1xdDboU4L3osj MNDlZXmzUwxF6TQ6Sr0cB+mCJc1O542own3bOfsXk6C7w2/ciEUcdWwt/LaL MA3tLCNjlXHi6B9ZjLZLcRak5P9L/SWzCQ2uQVHe03n4THuayfVoE7I1fbmk Y0AHz88CVyy1NiO2/8wH55eXoNPcbSb002YU532j95XkGtzb9aIl7gIZzcp1 zb3dyoL6JryTWO3JKOJJ7g1OeRYkXrjnvL8jGS2d2Cmht4cFpcrkLV90JaM8 tyv+TYdZUNoW6m5NXzLqelH1adCKBaULONyZuEVGy7P55wPvsqCMOaGz+pVk ZJSQ+EiDl4AGKtt21T4io/MhtheebSAgSY8oxt4nZKQd3bbjIJmAMqdWE7fU kFFCY+XM4Z0ElDU6/nr+NRkJCFxOuWBAQA9763bm9JER8iCfF4wkoOG7/nTx r2RUapsaWx5HQNLm+14n9JORzeTeTpMkAsruLLQNHCaj3NbOwtQcJrfeoZpO kVHJq+RMvQbm/igj25Y/ZMT9/UX2fCtzvz5x5+EZMho5v78h/wMB5TQEN8ov kpG634nTQkMENBKqTs1gkFHjsk97608CkoF5G5EVMhp7LXI+YpqALq2VKcSs kZGS0RVhbQYB5dY4L66vk5G0QPvQ+joB/Q+Ozrsd "]]}}, {}, {}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->{True, True}, AxesLabel->{None, None}, AxesOrigin->{0, 0}, DisplayFunction->Identity, Frame->{{False, False}, {False, False}}, FrameLabel->{{None, None}, {None, None}}, FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}}, GridLines->{None, None}, GridLinesStyle->Directive[ GrayLevel[0.5, 0.4]], ImagePadding->All, Method->{ "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> AbsolutePointSize[6], "ScalingFunctions" -> None, "CoordinatesToolOptions" -> {"DisplayFunction" -> ({ (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[ Part[#, 1]], (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[ Part[#, 2]]}& ), "CopiedValueFunction" -> ({ (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[ Part[#, 1]], (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[ Part[#, 2]]}& )}}, PlotRange->{{0, 7}, {-1.853159467604671, 3.1999995809523973`}}, PlotRangeClipping->True, PlotRangePadding->{{ Scaled[0.02], Scaled[0.02]}, { Scaled[0.05], Scaled[0.05]}}, Ticks->{Automatic, Automatic}]], "Output", CellChangeTimes->{ 3.6198896140618563`*^9, 3.7459415494882455`*^9, {3.7459417576832657`*^9, 3.7459417775316563`*^9}, 3.7459418803402133`*^9, 3.745941914536869*^9, 3.7459419832599325`*^9}] }, Open ]], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ "this", " ", "is", " ", "the", " ", "unique", " ", "degree", " ", "3", " ", "polynomial", " ", "that", " ", "interpolates", " ", "our", " ", "4", " ", "points"}], " ", "*)"}]], "Input", CellChangeTimes->{{3.619889703967351*^9, 3.619889735308918*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"p3", "=", RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"k", "=", "0"}], "n"], RowBox[{ RowBox[{"y", "[", "k", "]"}], " ", RowBox[{"Lag", "[", RowBox[{"k", ",", "x"}], "]"}]}]}]}]], "Input", CellChangeTimes->{{3.58861471522016*^9, 3.588614716124977*^9}, { 3.588615424760027*^9, 3.588615424869229*^9}}, AspectRatioFixed->True], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"-", FractionBox["1", "5"]}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "6"}], "+", "x"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "4"}], "+", "x"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "2"}], "+", "x"}], ")"}]}], "+", RowBox[{ FractionBox["5", "8"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "6"}], "+", "x"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "4"}], "+", "x"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "x"}], ")"}]}], "-", RowBox[{ FractionBox["1", "6"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "6"}], "+", "x"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "2"}], "+", "x"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "x"}], ")"}]}], "+", RowBox[{ FractionBox["7", "40"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "4"}], "+", "x"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "2"}], "+", "x"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "x"}], ")"}]}]}]], "Output", CellChangeTimes->{3.6198897566599855`*^9, 3.745942010254305*^9}] }, Open ]], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ "the", " ", "graph", " ", "below", " ", "shows", " ", "the", " ", "data", " ", "points", " ", "and", " ", "the", " ", "unique", " ", "polynomial"}], " ", "*)"}]], "Input", CellChangeTimes->{{3.6198897489596004`*^9, 3.61988978134122*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"plot2", "=", RowBox[{"Plot", "[", RowBox[{"p3", ",", RowBox[{"{", RowBox[{"x", ",", "a", ",", "b"}], "}"}]}], "]"}]}], ";"}], "\n", RowBox[{"Show", "[", RowBox[{"plot1", ",", "plot2"}], "]"}]}], "Input", CellChangeTimes->{ 3.588614428686929*^9, {3.588614617578576*^9, 3.588614618041622*^9}, { 3.588670360724497*^9, 3.5886703615978727`*^9}, {3.7459420247286344`*^9, 3.7459420333412237`*^9}}, AspectRatioFixed->True], Cell[BoxData[ GraphicsBox[{{{}, {{}, {RGBColor[1, 0, 0], PointSize[0.05], AbsoluteThickness[1.6], PointBox[{{1., 3.}, {2., 5.}, {4., 2.}, {6., 7.}}]}, {}}, {}, {}, {{}, {}}}, {{{}, {}, {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[ 1.], LineBox[CompressedData[" 1:eJwVl3lcTV0Xx++9zfM95zRJkcpQkR4NSmmvBlERDUoUKVRSmjOElMjUREVI QkSRJESKIlIhTZ7SRLOHRtU9d3j3+9f9fD/n7rXX3mv91l5rgc9+590sBoPx iMlg/P/3xpb7TxOXxls4xG+sXmFAVd7iKMbdzNiJcssSjGrnUJUbzopoH8kI RVZ804+BDKoy92sMNyDjOKp8YOhVX09W1kWb7PHOSEUaX1VGJALIyh47xSvb Mm4gxXkjun43iMpfKxKMPDIeIT3vm7p5c9iVoXdesrZkvEZhPjsG7m2UrNTT K7xc7NKIBpZmneIs4lXoVoyFHF3VjrKTj9it8v1RwTz1OdCB04M0P/LrjKo6 UP5F55Uipf3IZPVlhfal0yjjyMbcHvURZKpT5neFxYJGFc9aMcMR5PXa5baL KgvknvmPL1s7goxEzy3mGLEgcfy41aGgEbS0uECKDGDBYb+SXrJsBD2xOvwp qIEF3k5zNG1cfqHVVcW639KFQGfhj5t5Cf+hc9pPA7KVRKCiLup24PAfxKiJ SO8YEIMDKHePCe8P+uD/3LCPKwbK3/ovsMRHUeSHM5++E+Lw2YHc/EllFKU+ Kha/biYOlv8EtPrDKMqJXy7YmSQOGlzFjitnR9Hy3TmF/+lLQF9qWB9TYwx5 fTqQYbpXEgLLtWfqHcdRyocv75MfSMPRffxrae7jqJ13viL1lTSkzP1q7e49 jsb9VR+dbJCG0oMxyd2h48jceGea3W9pYBh/WjRxcRzZjvoWGiyTgYwHEZuV 28eRAUfJ5updGXhzvaLYx28C+Sw8O9ySLgvNjhe3LA6ZQP9ah323uCULgzx/ /siBCdSkV2qSVSwLMp6kfeSZCSRRLmKs+0kWtijv7jlVOIEO867eLxWTg/9S JNmFExMoZlHU8agoOVCO2xw0fWwShf0iBrYsZMN0wtbgg6cnUWx8m63TUja0 nN6xn06bRPyBqWNWBmxITw0IFdyeRBMf/pyUs2QDeeNIpFjdJNr9jqPl6skG qcrbR5SUp9BLW5szp1PZwOdNnV9ZNIU6k7dCyV82fGfSSc+fT6FwP3+ROi4b XoowUsyqplBdau5QB4uAQzKSaahlCoUluykPyxIwrTovYy13CoFQm0PxIgJG zWyz3df9RVZhSpO0KwE/DmY8iO7+i5wurU6JLyDgQPi/k+5Df9EWgz6vvY8I kA2aZ2Yy/hdFqDwMcCglYJV3Xs2M0DTqnK2y4lYQkGJb2n1w8TRaZJK/dM5X AiyoFjImeBqVnpobZDlDQFaBQnQsdxrtVg2W3WJBgn6eR7m3yAzqFUgTZ6xI eHf9mpCl7AySybx+ptSWhLG0hSlM9RnkK1M2xHckYd1Bo/w46xlU7Hr49IYd JEyt2dx+4swMUtIc6vhxlIRNXRfRaeVZdNQfLJllJNyy9m80XzCL2opEhsPL SZi+Y7Z7VHsW1a9NiO6pJCE7pPeMu9kseuSxa21hDQm/mPotC7fPohMtB4bE mkk4rfVx35ubmI88Kqz9TcKbAGYWrcdBm2upX6nzKFBsaFr6cCUHbVJpD9i1 gIKAFXcrfICDTA/9ZK7QooBNO/Z9cOKgXzuDT73SpmD7uav6lyI4aF6WkdI1 Qwo4D1bWGJZxUNPftZWhdhQYTAZNBq2hUe89vd8n9lPgWhja2+lIo8hqzTm/ QymI2BP5eeMWGnkWuVY4RVBQ0hZT8E8gjb5zpxLED1JgWH5211QKjc5skTQx jMN8Mr8ppoNGpy/PkT15AdtDhW/+66OR2wWD/zLTsb2ZoqLtf2gEbs8DbmVi e3ufnQMWF2kaXbQuvILXb6yxEV7CRXzbbaVHb2JW7is5G85FYqsaojSKsb0v g7ncGC66Giu9rOsxtnfmV0rQSS5SOdWfm/EE2+NO7Nt4mYsOf5ltnnqG1/ew FlIVXFTpEJoeWYH5nnp6liQPWW623hLzEdvz1YqXkuehZZtvi/fXYXuqS0Jj 1Hjoo9TXLPsGbC9p+Ybt+jxUk6UWJPwFrw+3ENFw46H7y61FzFooMDL3jMy/ wUNzSvb3PumiQHok++6Bezy0rfCCzMtuCnov97SvfYz/P3yjp7yHgpTpPVb9 1Tx0sSj4bMkPCoaLQ+U0B3ko5IevYvAABZU7S6zGR3lIKDriudMgBZns6cjX szz0c70z9c8QBTbBMR07pPjomYXfeO8wBdnap/Kv6vHRwZKUJOo39r/tQ0fg Sj4KbLqt3IzZ4ZQ02wz4aMQw0C/tDwUzP1Ojvjnx0WvjGQ/GGAXOOVetFSP5 KOvpk6r8CQqWbOyK6jvCR5ObB27aTVIg4C24V3KSj+YbFLj1YS7YdoftcomP ZLV7VrH/UiCmVPw9+QUfLTuSpmc4Q0Hn20n2jmo+4p2d+VqMuSRipY1ePR/1 tDZ5LJulwKfx5b26Tj5i5fXIKHIoMDku+H5lgI8aB86gE5jl9K2IwFE++tp3 ZctvzC/P10RLsAQoM71/bSlNQZq55P02SQEaFr28QJ6L839kfecdSoBSH98a CsKsaNdoY7tQgDirnzso8Cj4NS1/QEFPgJzaVw35YK7Kc7//01iA5hmIxRRg DhH+TsSvEyDT1NBEAz4Fto/nr3F2EqA57okioZjVfHwOLNgqQElu7UfvYZ5k 374/6iNAPn8W/e3CXFsx0FkRKEDbjy8OIgQU5ATrkMkRArR23+U+C8zRakFr th8RoH45053+mB3rHh5YdlKAJOre9JzHrHV4/D43SYDCGRP+DzFztI26PmYK 0Ij0yZk6zF/aoskrOQJk7mSS1o/57qmyNXvzBWjZsncruZiPGvMOmBYLUALj 17AMQx4296EC8RcCROjuKZyLeenFuK7WKgG6ID0WuxAzy/oteadOgNqE7fx0 MX8bE7ONahagW+E6Psswx0wdamJ1CdAd3dWq//+uPvufT/Ig/t4Uj7QwV3G9 x+aOC5CIn2GDMuY9jKZjd2kBWt5+olYcs3Cc5qksFgPOa983nsT+tQh5NXiL MSDB/rt8+//9T8hQWCzNwH47+ZVjPij22fM/NgO6+sz1r2C2Py1x67ECA+Sd WkIjMKtKWY8cVGHA8RAXPTvMv8/FrID5DEgxH/abg7lCtvSgqBYD8gs/q/fh eKSk/KmsW8KAFS/NtxVg9iG1xS8sY8BF711z9mMWUbyaMX8ltrf+lfYAjn9L ZvP3PjMGdJzwjL2K+e4cuYUFwIBHhkm+GzA7qMU9NrFnwM6O+LYbOL9Ur7/g 8B0ZIPX2cfAazL/VpyzfujCgQXRj1k+cn6la/p+dvBhgFn29TAGzT94NJWUf BjzJqHx8B+ez4ZL27Z17GPBzVMvBCHOrruN/e0MYILEy2csa60PNwFDyRDwD Xhl9cdSdxvZLgpzsExng4DYRnYX1Vml85xL7PAPqI8SdRDDvWqWy+FoGA9zO /5jbgPV5z5JhU5rPgDfuNs5LxymIqTI9G/OAAe4X+kwOYr1vWBPeaPWYARop 2j2vRykYXdfv/eklA6YfhB6zxfXBeFPdkcFPDBgtGNJX+0WBeKPI2wdNDDit luizboSCby5IOvIbAy6QB/324/oT416cxfzBgNxAnckiXK/ebL/0VOUvA4RE it1YfTjfg3aNrVdlQv7ypOXdHThfb/V1ctWZEGY5XNnQTsHD9j11BQuZsFXT O/j5v/j9sg+4I72cCcN20wvj2yhoXBTsVWfJhHVVu6wGv1KwuTv6g4MfE6RD PgjTHyg4qTRbSgcyQfS/aueC9xSUOh66dT+ECYNCzp0eNbgelMcckzrEhM9j R0vzqylou3zc6OM5Jvjuf3JbFb8Hni5nb9gXM+E+Z4GrE35Pzp2RTuaUMiHx mGlUzSMKyl+fj7n3ggnK/lunTYsomKef4i75lgmv3BTnyxdS0C2dLlPbygRX B3F0LY8C33fZB+14TKjI2pkufJmCvabFTuvWsSBqrjMYH6YgaqRx5531LHBZ 9uDEdvx+xl2bCBV1YkGXo3Z7fDQFl5lGadUeLPioZ/v2bTgF7z88bYRAFuT6 XC1ftI+ChVvLXUyTWLDzem3Zdk+8/6H3m3WaWDCS/qCmeBWud0uHdp9uY4Fr /ZHLsSa43ndKRA12sOB73NPHDsa4vlo7ZOT1sUClxeD9t38osJBqaNGYZkG3 WmNO8xIKrl756j5XRQjKVj45t0GRAveXnR7S3kIQXp6+3P4PCe8vZSrV7hKC Xv17Ged/kbAqclPzqQAhiDuRINUwRMI8vTebhMKF4Gma6WLbnyT0X7+9lpMg BF+fGvpS30iIit9nNHhfCDae9WnUfENChj2HXT0tBHZnp0t000hoalN4fzhZ GBYkd7/dpEOCRoFu5+N7IvBQalAqfCMBAuX1WzO/i4JRRdPBGBs2ZH738P8p KQ5BL5zZLk6ywPCwncfRlIDRivRC/mJpmAlaveqFqyRoppj2xL2TgNemFN8m VApsdfpnky6IwRqvFheTRGlwTdLLO3RMBDxUbsm1lspAlKtqtj4+p0XSpnCL L7KQORzRbnedCRF7o1+JtsjC0PoqofmXmHDPNlvi07+ycLJ1pmEmhQkKgpHr O37KQldm146yOCb82n+qLnZaFmq0xAYj9zDhitOrhVVqcrA3aE5kox4TZuSX tq3dKwd5+RrXoIoBemPOmuz9cqC1z+1VMdblrvqDwW3hcjAymMvWKWXAp4Qa 4YAjcpBuvX7rEqzzvL8+/5xJkQOJ7w5832QGbG67fKauVA4isj32tXsy4Ft9 fPZdJhs46bLTX6dxHVc0y2QKs2H2yN+j638JUPCO8eStomygPzarfO4WINaY 93FpKTb8vqFWNPNBgJbKW+wKkWfD0L615p+vCtCxrTPaJovZULuq0DjHRoAW Dux7UrOeDZtT5VL5V/lodLnWA/WNbDCR6+gdSuOjsgPteQed2BC14tnmnkQ+ cpS0v6Trxob+EKvKyQg+OrBs8eHkHdg/YfmLDzbw0cfwHnAPY0P8ms1pkUw+ CmW41fVnssFUKPvd2TAeuifq5aaUxYbE0eqEa3t5qFd6V/faq2xoyGt9UebD Q85zwibv5uA56I3UXSUXHlqxIkk1MJ8NPyVfIU9DHhr3fbfvTxkb6na+3aU4 y0U6e+v/zi9nw466mLe9Y1zkE9IUu6mCDU3OGyRfDHNRY0zvxUdVbFAUT9iW 3MFFxen8l+F1bFDfvb29ppKLwmqMZWa+s0EuzMDu6Fkuule/OnNJNxuOJW6a J53ARb1fbRZ49LLhqUyJTd5RLnLudjYq62eD1fPLOYJQLloxG+wV84cNH/p+ u6z14KJxnTuFDCYBJ9TMbujrYPtJSo4SWgTc3rA8K6sZ9935w3oBeI46tHIq 0v0zjda9LZf7sISAyK1pnzQ/0mge1/dL4jICXt9cMH+kkkYfA4pcxVcSUHQP t9iFNDJOiDfyNyXgffXenuC7NMrNcVN8b0ZAeY2U77ObNDrQym09BQRodB/9 eiiLRlpr7LaJ2RPw+XWT56tEGiV7q5r7rSeg+/RWnegTNOIc/qNa40hA4vAk A2Jp9Lk4vfOkCwFB9ca/mAdodES9d6eoFwF5sVkP1vrTqIU+5C8SQsC7JvPY XgcaWSo52u0OI+B76JaJ6+toVLBigc7bCAJmRN9x/j+XxAXUjJw4SAA7dM+Q GaLR8lZqv3A8AdX3/UxzDGiUNd6/cVcCgXX2yuqHPo1EZcv0q08RUNsusmil Ho06bLzH488RkMYyvkZo0yixuCBSKIMAk4l3lmnzaTRRf8zN9xIBWp6Osl5q NNo+5LyyKouAY6q7OSZzaWSoPjsTd52A7bScq5YSja6b1X3rvUFAw1uprwYK NJJ0v15mdQufH+256ErRqOf8mhhWPgEvL4kya+VotD5f2cvnPgFqWiXl82Rp 9Kx6ZPWbQnwfpcldJ6RplESnCo4XE8C0y4lMkKDRrOLu7p4SHL9Jc3q+OI12 rTB5bfmUgOOCFvE6URqtCuiMY74k4Ofpxs4twjSirJhTJq8I+H21Ys16IRr9 UtHyC6kkQGaOM7GRRaPqCdtvd94QQKzIXreTSaNrdQEOXdUEvJr7bDyOQaOo 2+fKFWsI6JcoUSwRcJDj0YfLHT8Q4Kyf92KSz0GL3RtvJHzE8X+TO2iDmaE/ RZXXE6AzVX7rFo+Dvokrn5z8hJmWmqYwF/esmtFtJCBX/EZ/GpeDzpR57fVt IsDTKfXQAsy+F2I7sloIWCQx9ryC5iDzfTcdG9twPvg23Q/CLL/mXaVEOwHS 6RvddTH/pza0wvI7/v45qmqGw0Hv/krdPtBFQJn99ulmzNc/6SkV9RDwwliF 8wbzgbtOpwd+EHDx/avGCsybjkfQ8/oJkDT0TqzDrL01M8htkIDJAlJjCDPL oKzr/DAB8n7DeQp4v3+lvju9/UUAZNBqzpgf/xRUcX8TYBq8NSkb87lyDWPD MQKaNeYKOJh3Z6y5GzhBwL8DTtF++HwW+/1Vbk4RoDBG8n9gVlx39ty/0wS4 hYVlheD7+aP+gE9wCNh1/YSLDL7P97OfQ+y4BMgWeOmXYc5pnOiN5ROgXkeY RON4OJ8wrfnDIqFQwfPnAhw/XS9P08UiJJTnqWTK4fgKGx+7v12MBH8VkVQZ HP8nA9XJ9VIkFB94ELVKBOdX5QBTRJaEX7XpIQE4f/ZclowwZ5Pwu2G2/q4Y jZQdNm25L0/CoReJn7ZI0mhMM7y2V5GEXKmO4x+kaFTLTTdXmUNCZQf3sb0M jQ4/aFc/rUaC/Jfpx0fZNOok/Ib2LCIhtkB8qBDr4+nw6W3ZS7B/D3g7bs6h UUpVQX0zfsfdPu9Mzsd6sowcf2SznIS38N5gEOvvVtuRQxom2L9Ya4/lWJ9H HuX88lhFwoYlAx73dWnkdqZqe6o5Pt/5Kw4mWM/i5hLWDEsS+DtOzUvA+g+8 flGq044ExpecdSMWNJp0bt8kuZ4EJ5n+ld2WNIoR1cgwdiQhpsj7dp8Njc4F PVRPdiahU+GmggauPwVmtUawDZ/nnaKozxas9z/sw/u8SDiT+/pG+zYavcx1 r7y0g4SaMnp09w4a1Uv02Y/5kqB8ura6YA+NfrcKduTuI8FR9DxXKArr6azt 7YZgvL/Ovj9Ch2jEtzg/zAkhQXqH4Uf5ozSSy1OJdIkkQbNv6ZaAkzTSjzA6 I3yUBEHdgbjOTBqFsQNL9iRjfzdv/M0vx/W16tFsWioJPXxWCrsK18PoGYuK CyS0+bgkGLyn0YXOhA9Kl0gYt4lOKGykUUlBTuf7HBLU2xMywgZp9Hddi7hu MQnWl1O/2Cpx0VGuqqN7Cc6P3a71IWpcJFrkeyG+lISUDL9TBZpcpKQ0ptZR RoIse6mThz4XmfRJGZyvIsHw48BUoj0XHTxu6fWniQTnAp6HbRwXGdq7ztq3 kpCd88+g3hkuGiX90vNwnyf6W/jb0jQu2nPrXL1XJwmWcdnCAblc5FrTurpu gIQ0Sj/sRxUX6ckEz8vnkCA27j13UJKHhptjy4R5JIR487/oyPPQ7ewLbt4C Es51J4XHq/GQqv7zJCVhCqSEsy9H6fOQhIuwIEGGgiXbxLx3u/PQz0tZXT7q eK73ybLfl89DOTsLD7/SoKCoqeMBUcJD23QqlVQW4r6X2HOy/hUPNb7o2/hF G/flrTsfJTbxUEWnfiUywHPS3NQ+Gdx/ZGnV5KjaUrDqso/Ztu18tLFoYmdL IJ6LVoRZrtYRoAPGbskSwRS4rHeX+mYsQDkvn700D6HgR+eKDaetBWj8/RGl mxF4Tus+7LrYS4AyesTq9x+h4MwRFYOtqQLUSaqaSCRT8Gf3Y4sbPNyPRdrI mhfjuUxmY+x4GAMy6Tyz/SUU2KoMXUo6xIDK4xIBuaUUdOibeBrFMYBMaqgS f0HB2dK60dxUBpTmbTnYXEWB17NHxfwiBvBb9/0MbqYg1iQ258UYnstXpZfd mKHA1LLJ0/cwnrOaqyyE8ZzqORMupJbABNvQ8Td78Nx7SuTakx9JTODf2fhR lyUPrsQyx0u5TAhWkOx4IiUPGoeKP2bXMsFx9Cjvwzx5GB5+9VVLjQXsPH+c i/LQpGo5YVDHgn2WmdWua+VhXdzh9MxWFrzveLvuqZ08BD9z/Cz+gwWxlKZT jKM8KJ/cuUCew4Kx2O87RbfIg8WKvIKGJULQuM05XiVQHnSWKFu3JgrBBcL8 nVWqPGhtZzpsdhEGz71l7TEX5CHi9ZLc9B3CsKjKZKw0XR7UM/sdewOFoSzC SFU3Sx7sa2Of5pwQhu5WvTDqpjxU3C7ULi4VBr3sBfN+lsjDVFf7Ex01EajV EYtKaJOHyJb5589OiADLpnHh+/kK0Bg+ZtD8TgwOR0fITi9QgKVmdqmzzWIw eU9xeqGWArge8sjT6RODAWLb+/glClB0S0fksZA41HX92Iv+UYCTp60mPS3F If3w1MNSKwXIv9h1XrlCHJY8Vl51e7cCVD89uc22WgIcNXY4xt1XABV/u3em f6WAs7f+fHehApT8OdZmLCkNeY/N6i2KFCDKeVO39TxpENgor6dLFODuBVbv eVtpKPL7YhfxSgE4gYmLejOlgSq0XrPniwKoP5w8lLxaBtpWLjG3n1YADVft xLnpsuCzYUybtFaEoG7/yCncF3d+DfM4+VkRRJ6by519R4HGnVmzzwFKIFK1 KybghhJYL4oecGQog8mgQ1EESxXqGT/7H51Thq3F3f4NLuowf4nb5WGtOYA0 BdJ1OpogrhW6dZfFHODt+PcPr1gT/gf231iK "]]}}, {}, {}}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->{True, True}, AxesLabel->{None, None}, AxesOrigin->{0.921875, 0}, DisplayFunction->Identity, Frame->{{False, False}, {False, False}}, FrameLabel->{{None, None}, {None, None}}, FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}}, GridLines->{None, None}, GridLinesStyle->Directive[ GrayLevel[0.5, 0.4]], ImagePadding->All, Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({ (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[ Part[#, 1]], (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[ Part[#, 2]]}& ), "CopiedValueFunction" -> ({ (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[ Part[#, 1]], (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[ Part[#, 2]]}& )}}, PlotRange->{{1., 6.}, {0, 7.}}, PlotRangeClipping->True, PlotRangePadding->{{ Scaled[0.02], Scaled[0.02]}, { Scaled[0.02], Scaled[0.05]}}, Ticks->{Automatic, Automatic}]], "Output", CellChangeTimes->{ 3.619889758890097*^9, {3.745942012440316*^9, 3.745942035010105*^9}}] }, Open ]], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{ RowBox[{"Here", "'"}], "s", " ", "another", " ", "implementation", " ", "of", " ", "Lagrange", " ", "interpolation"}], ",", " ", RowBox[{"used", " ", "to", " ", "interpolate", " ", "a", " ", RowBox[{"FUNCTION", "."}]}]}], " ", "*)"}]], "Input", CellChangeTimes->{{3.588615210991409*^9, 3.58861525733857*^9}, { 3.588615361929764*^9, 3.588615395730233*^9}, {3.6198896419932528`*^9, 3.6198896580140533`*^9}, {3.745921773781764*^9, 3.7459217762516546`*^9}, { 3.745921816073552*^9, 3.7459218729593234`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"LagrangeInterpolatingFunctions", "[", RowBox[{ RowBox[{"{", "xj__", "}"}], ",", "x_"}], "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{ "i", ",", "n", ",", "xi", ",", "xjc", ",", "L", ",", "xgrid", ",", "num", ",", "den"}], "}"}], ",", RowBox[{ RowBox[{"xgrid", "=", RowBox[{"{", "xj", "}"}]}], ";", "\[IndentingNewLine]", RowBox[{"n", "=", RowBox[{"Length", "[", "xgrid", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"L", "=", RowBox[{"{", "}"}]}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", "n"}], ",", RowBox[{"i", "++"}], ",", RowBox[{ RowBox[{"xi", "=", RowBox[{"xgrid", "[", RowBox[{"[", "i", "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"xjc", "=", RowBox[{"Complement", "[", RowBox[{"xgrid", ",", RowBox[{"{", "xi", "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"den", "=", RowBox[{"Times", "@@", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"xi", "-", "#"}], ")"}], "&"}], "/@", "xjc"}], ")"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"num", "=", RowBox[{"Times", "@@", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"x", "-", "#"}], ")"}], "&"}], "/@", "xjc"}], ")"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"L", "=", RowBox[{"Append", "[", RowBox[{"L", ",", RowBox[{"num", "/", "den"}]}], "]"}]}], ";"}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"Return", "[", "L", "]"}], ";"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.58843061679911*^9, 3.5884306168021107`*^9}}], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ "next", " ", "call", " ", "generates", " ", "fundamental", " ", "Lagrange", " ", "polynomials", " ", "of", " ", "degree", " ", "3", " ", RowBox[{"(", RowBox[{"based", " ", "on", " ", "4", " ", "points"}], ")"}]}], " ", "*)"}]], "Input", CellChangeTimes->{{3.5884325909131465`*^9, 3.5884326411989546`*^9}, 3.5884328124170513`*^9}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"LagrangeInterpolatingFunctions", "[", RowBox[{ RowBox[{"{", RowBox[{"x1", ",", "x2", ",", "x3", ",", "x4"}], "}"}], ",", "x"}], "]"}]], "Input", CellChangeTimes->{{3.588430638692415*^9, 3.588430638692415*^9}, { 3.5884319763569508`*^9, 3.588431977236038*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ FractionBox[ RowBox[{ RowBox[{"(", RowBox[{"x", "-", "x2"}], ")"}], " ", RowBox[{"(", RowBox[{"x", "-", "x3"}], ")"}], " ", RowBox[{"(", RowBox[{"x", "-", "x4"}], ")"}]}], RowBox[{ RowBox[{"(", RowBox[{"x1", "-", "x2"}], ")"}], " ", RowBox[{"(", RowBox[{"x1", "-", "x3"}], ")"}], " ", RowBox[{"(", RowBox[{"x1", "-", "x4"}], ")"}]}]], ",", FractionBox[ RowBox[{ RowBox[{"(", RowBox[{"x", "-", "x1"}], ")"}], " ", RowBox[{"(", RowBox[{"x", "-", "x3"}], ")"}], " ", RowBox[{"(", RowBox[{"x", "-", "x4"}], ")"}]}], RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "x1"}], "+", "x2"}], ")"}], " ", RowBox[{"(", RowBox[{"x2", "-", "x3"}], ")"}], " ", RowBox[{"(", RowBox[{"x2", "-", "x4"}], ")"}]}]], ",", FractionBox[ RowBox[{ RowBox[{"(", RowBox[{"x", "-", "x1"}], ")"}], " ", RowBox[{"(", RowBox[{"x", "-", "x2"}], ")"}], " ", RowBox[{"(", RowBox[{"x", "-", "x4"}], ")"}]}], RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "x1"}], "+", "x3"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "x2"}], "+", "x3"}], ")"}], " ", RowBox[{"(", RowBox[{"x3", "-", "x4"}], ")"}]}]], ",", FractionBox[ RowBox[{ RowBox[{"(", RowBox[{"x", "-", "x1"}], ")"}], " ", RowBox[{"(", RowBox[{"x", "-", "x2"}], ")"}], " ", RowBox[{"(", RowBox[{"x", "-", "x3"}], ")"}]}], RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "x1"}], "+", "x4"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "x2"}], "+", "x4"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "x3"}], "+", "x4"}], ")"}]}]]}], "}"}]], "Output", CellChangeTimes->{3.619889385099207*^9, 3.745942108004756*^9}] }, Open ]], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"f", "[", "x_", "]"}], "=", RowBox[{"Cos", "[", "x", "]"}]}], ";"}], " "}], "\[IndentingNewLine]", RowBox[{ RowBox[{"points", "=", RowBox[{"{", RowBox[{"0.0", ",", "0.6", ",", "0.9", ",", "1.2"}], "}"}]}], ";"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.5884308686807957`*^9, 3.5884309011966624`*^9}, { 3.588431087068935*^9, 3.5884311327957015`*^9}, {3.5884312858621244`*^9, 3.588431287166255*^9}, {3.5884313649681005`*^9, 3.5884314859766254`*^9}, { 3.5884315613612556`*^9, 3.588431568195939*^9}, {3.5884319703333483`*^9, 3.588431991438819*^9}, 3.5884326555494957`*^9}], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{"next", ",", " ", RowBox[{ RowBox[{ "Lagrange", " ", "fundamental", " ", "polynomials", " ", "calculated", " ", "above", " ", "are", " ", "multiplied", " ", "by", " ", "y"}], "-", RowBox[{ "coordinates", " ", "of", " ", "the", " ", "points", " ", "and", " ", "the", " ", "products", " ", "are", " ", "summed", " ", "up"}]}], ",", " ", RowBox[{ "resulting", " ", "in", " ", "the", " ", "final", " ", "polynomial"}]}], "*)"}]], "Input", CellChangeTimes->{{3.5884326744893894`*^9, 3.588432756433976*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"poly", "=", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"f", "/@", "points"}], ")"}], ".", RowBox[{"LagrangeInterpolatingFunctions", "[", RowBox[{"points", ",", "x"}], "]"}]}], "//", "Expand"}]}]], "Input", CellChangeTimes->{{3.5884314999610233`*^9, 3.5884315266618214`*^9}, { 3.5884315842675457`*^9, 3.58843159455448*^9}, {3.588431906658903*^9, 3.5884319280030375`*^9}}], Cell[BoxData[ RowBox[{"0.9999999999999999`", "\[VeryThinSpace]", "+", RowBox[{"0.02271748399367368`", " ", "x"}], "-", RowBox[{"0.58434428998366`", " ", SuperscriptBox["x", "2"]}], "+", RowBox[{"0.10217161531292529`", " ", SuperscriptBox["x", "3"]}]}]], "Output", CellChangeTimes->{3.6198894202788553`*^9, 3.745942111761964*^9}] }, Open ]], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ "this", " ", "is", " ", "the", " ", "comparison", " ", "of", " ", "a", " ", "function", " ", "and", " ", "a", " ", "polynomial"}], " ", "*)"}]], "Input", CellChangeTimes->{{3.5884327684925294`*^9, 3.588432804381879*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"Plot", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"f", "[", "x", "]"}], ",", "poly"}], "}"}], ",", RowBox[{"{", RowBox[{"x", ",", "0.0", ",", "1.2"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{"Magenta", ",", "Blue"}], "}"}]}]}], "]"}], " ", "\n"}]], "Input", CellChangeTimes->{ 3.5884286943944497`*^9, {3.588428818404659*^9, 3.5884288202128396`*^9}, { 3.588428959268738*^9, 3.588428959960758*^9}, {3.588429011049819*^9, 3.5884290275341544`*^9}, {3.588429457279316*^9, 3.5884294733917327`*^9}, { 3.588431671495558*^9, 3.5884317373169475`*^9}, {3.58843191213245*^9, 3.588431932639454*^9}, {3.5884320042990775`*^9, 3.5884320055391026`*^9}, { 3.7459218870535326`*^9, 3.7459218879962416`*^9}}, AspectRatioFixed->False], Cell[BoxData[ GraphicsBox[{{{}, {}, {RGBColor[1, 0, 1], AbsoluteThickness[1.6], Opacity[1.], LineBox[CompressedData[" 1:eJwV1Hk41VkYB3ClGDTJMpRpbEkqWZO9b2KaajLV7VaWhrJkC4OkUpZUTCNE 6mYtlC3Z9ylEkTZr13ot2bcf92JkndN5nvOc5/Pf+7zv9z0yVq4029VcXFyx 5H5/OV1iynm/m+1dWvl+KMQGSEmvk9KCXXxC+3e7bi0MYkoehOZcsMt3Bxll +omK02AjUFa0TGwyu1r/paQJgtu0SpaIFVJOfbOVt8CU3bz7IvGiYrMZr7gV rhZwf50n7lu7UtQuYYv8d3T+b8QfWQriWZL2aHo6PjRLXFBI87wp64SMo598 ponjwq41msi7oOzD7Psp4srdjQyxyL+wY4td3QSxpsxQNZe4OwpNpYJHidN/ XJodYXigzktqZZBYal5IvlnCE3/fstveR8zbqHszVfIyTCM+J7YTXys7lhvx +Ao+PaM4X4gn0217r8t6I7bu9EI9ccvN0H00eR9IR3KpVBMfcUty1UvxhflB gz8riLtfnDkqEukHLsNag1JierSbz5LYDWQdyFd5QRwitBjJZt6Ai7O4bjJx TeDtjEFGAObEXi/GE+tdjO6ol7iFqd1TT8OI5Yzf6DyTDAJPeFGoC7FF5VFa dFcQtMvP0m2JGdptDmGP/8YCD63OnHidPMW4KvsPTmpsmPiNeHp5439/yIfA L3qpfDOxkmfieqPBEDQLTkGI2H50l7x2SijemytdWks84fQkTSgyDA4Dt38Y X6ZgZVLQuiAWDqXOnKoiYt8Ib9tLXuGolZ5bSiOO+bRvapIZDge5oIkY4i9G 738YYERAfixoxpf4sGq3Zp1EJAblhKL3E6vz8z9IkmRA7b65UckSqT+oJ1Ev kAHHS3dtnhEv8hRnN1EM2DcZK4cTV3Hbf1xT8QjeHVer7InpS2/W2FpHw04v eFyI+OKkv8fWlDic2lKuQ18k83Ex9X+5IR71q2a4dYn1xlVCT16JR5T8v14y xKtnXJv3b32MUfpqy/EFCrVTklPus4+xl+Op6E9sPnpNoZGRgE2PdCoZ82Re afbKUh8TUPg+X+MqcbkDfY/TqkSc5lTTzInlhnYacTsm4sqFnODNxGN9bZbq ukkwNEpdifpG8sTSehjR+RQbzjv1Bs6RftVPr6XLpoJPghkrM0vh5hO1VoUT qUjOGZFZmKFwwe2v50sBqSiNu2feRKwvNEZL7k+FdHsUdYuYdexr/HxKGlKo pKj+aZLnunrtJ8rP8U5sZPEBh0LCpxfOE3qZiLSYnK6YpHAnbmxfpXMm+jx2 KIYTe7jsEGXEZSLPc80mK2JZ3hOd0lxZOLzySGQVMfvZp4Dac1nwiTgdrEtR CB94U/fL1my4GdwIThin0GSb51iVngMpg/5arREKsy1RWnmVOfjR1uwLF7H4 EX+epPYcoO5cVM0wBTO1PxICBHJx6+L2I6eJu5YHW/c75YK3O/WU2xCFkYc/ H6rYkQcqrZuXMUCBq/bGtlep+TAUb2Bl9pL69B1mMiryESketuRCbJh1tDK2 NR8Wx6aadhHffrDZ8jpfAY7nxmek9VAQsC54oOdQgHOm//I86aYgtji8tlSh EOm8l038WBR2KR3vK0wuwqZCFR++NgrSav8KZ1YUwVtwxbm4lYLonm0Gz9qL QPcZ5XMgntdfir0vWIydlMyB6hYKb41TT7t5FYPy3lzny6Rw1nnV+x2/laCK n2dnTxPp1/Ps7NiBUszffbva9jPpzw5h34CtZbjYnuHsVknBSPOZSJl+GRKU DimKEIcb6qTMnyyDuwBfQt5rsp9nrOrdbpdB8O23nNkKCnZ3c7acHSzDZelC R69yki+v7TZx9HLwMazUnF6S/6Dt3Q/eYhXIN0xR2FxIwSCen7Y76jVi2l37 ElMp7DOUnAnvfQMNjP51N4jC0JEtMnE/v0NtVKfuoBWZ97lDwnt++YANOq5Q 1aWw+94q0RLLz+iu81gnL0Bh3NbZ2VajHsftOPxPWycgWkD7ySi3AVdcAlbX JE3ARMT/DeN4E2xDDvdGnZ+A8Sxa1JuaMbDk8nCTwgQsBn7jGHkx0TB3Scun axxGf9fQngu3QrBYtdkkYhzt2a2iZuFtsLT3iuQ9OI6QQfWu62odGD3VUZU9 NQZ2guS6oYJOXKzXlvgQMwZ3J5pBjXkXTrzq6lLQG4O1h++fvte6cUWuQzap ZRT9nNj3+vo9CG5amsj3GsXMfLe+5vpeCAk8rNLiHwXL6ELPifpeBPlIp4/H j0B1SD05OuwrlDV8eebURyAsLO0mebYP/H53eAtLhvHBh36hSqYfMzTxc2W/ DkPUPcPpKacfHy13+ve/HsLbV/vvZlcMwDzn+JjV3iEc/0O6u9t/EOs0YjT2 lQ3C0UZW83faEHp0bp95pDoI85cXylpEhxGy309tIHUAha4ch3uNw1h/1pKd LTOAJylxG28ljKDlM3+G851+WMuFz0VbjyJ54L6j5HwfcuXTOsZ3j2E4Mqrs qFkfOpz42VZzY0g8kpXL/forXPmTt3MXjaPv4GHtoe1fEaZsb8sKnAArUMas 178Xdx66NXQcogCtadVjvT3gr/wsI7B+EoNvK7h91XqwZvCmk77nJO73H1wU DOiGoyLL2L95EpqtjcYLxV0Qzz+U06s/BbZ+uv3bIyyoWJpmnoyfAk13KMuI 04FljUB91vIU0pX2/+Qd0w6RY9XGN8+zYZ576YSuXBsi/0mf1a1mo3L2m9iB ihb4WQQ0jtawQbf6ZWZdaQscVc2zomvZWBPTFdyY14K9TH7HhY9ssOiBw5Yp LRiSdewsaWJDnYdu6hnaAp0ShSrtXjbcDc/LRvzZAtbQ03uaS2wkJuV/ifrG RE3pdZfBZTZ2JYgImHKYyAk5+ftDLg6iyw8uiI0zEbh77dr/uDmojvUVDO9m QsXP5nIBPwcGucJ8t6qZuCG+xVJjIwc24XbHLCKZuDAyr9u/iYM9duGHN4Yy ceplw8bInzloS8te1RDExE7rGw3TkhwYj4kmG11joulFz695Wzm47G9VJ2/N xCv/YlmbbRxoTcnlsc4wkUq/tyKynYPrrxxPPjjFhM/8vmIPRQ4k/quuW3OY CfuPGx9sUeKAaxt3eqkhE7THk+6NyhyE6vMZuuszoedRczRAlYMByTGGgiYT 2w48VlRX56Az8+VzlgoT/wN9iSfY "]]}, {RGBColor[0, 0, 1], AbsoluteThickness[1.6], Opacity[1.], LineBox[CompressedData[" 1:eJwV1nk0lWkcB3BFJWlDMSVrGZRUlMry1YSispVGM0Qhsg4GU5J9q6nsbraQ LEV2WSbipkS2qHstWcMV7uVeySjNPPeP97znc85znvN7vs/vd95X+rKbqd1K Hh6eV+ThvjmDW5XLTv2mZa4uRzSL1GBJKUHJwyjOOVbIT+y261kETeIklvqH VDYTmy+s1HwuYY73/KXPpYlbB+RFiyQcsE5luEubuOKZqVeIjBM0y8yUTxKn Rd3oMpdzhVnLRWsjYjVpxmseUQ+USKsdtyRe06UekifxFwI11AV9iG/UGZfG pl/D8aP3RPyIZ5/YjfjJ+GI4KYIdREwPuadtKncTylOeyneJT7tnuWnk+mP7 Bse7scRDTy2MhOMDYL20nkohPpfsfnN5axBc4+0LMonvbv4ez6YFwctH2iaH uCk8rGCCEgyjyGjmE2KNP5P7O7eF4t/MlLAyYu/Pu+Zf9YVi9YB6bCVxkXXR un9SwnCp8L7XP8Q7zzQezZaIQJDv/Roq8UWqkWnyYARC8iD5mphypPdqVHok TLrzzZuJ3xXaBoZeuoXQLaP2rcSCcizKdZnboIZ9Ot1BrJdyrcht9DZyjpWu 6SIOEOJrss36Gx1WJpT3xPM/xL4ayt2F2x4+9V7ivV4PN+hM3IXshh1n+okd ppTkjuTeA6Vq/cEB4s59/RKJqlFwPklnDxIznTIeb46Pgnb3zbBh4nU5Vw79 /SUKhd68syPE8iO7G1afj4aFvoPyJ2LdHXNnAiuiMWZXqDNGfNm8oufb1hgw xnqVx4n9Y33tvH1iwDs6Pct1Spv23CwtBh+8J8IniKvWrvZzOhwLnpq381x/ 0GnhH6fEQvJNmhqDmOMfFWf9byy6q62Mud5UYybVdyEONXmbwbXSwrZ8s+o4 eBZW8HJtsH9IrWNbPFLGjNO4+9k7P6Ia+Mbjg93IBq5DchyNGvviUafvYsat J2NEuQ8aCZh/zHHn1l+748uV6pQEbHnoY8s9X595NVt1OQENp77tGSVejPW/ WWiZiKTagBZuPlvbdQQUaxNxZIuA9hCxioBAQpYEBceu3r/FzZcSMfxQI5yC 61ROADf/76urirtZFKi5LOty78c6NKrO2fw+rgtX9dCIX/I6tPLV34e+jYIG 9z7lg9CXopCEvLUmLu+I/+YRnVSNSUJS6QGXdm7/Ljfy2dkkw/DV3p4m4krf VKHllmSIRRrqNhKLL/0pFa+agnBelYB64k8LshqNq1Lh0HXuchXxSa9v+pau qRgJSlxbTpzPeffrlw+pEEnL8y8i/nM20HNXbhoKWvTo2dz5cr0Q+HzTA+w4 TKvO4M7DzL57Ztce4NeXx5xTiFd+cXv/y650qM4WqkcTN89JzHkspEM0SP3O X9zzFwks1ipmQK9nYtyd259/LPwQsMrA5crgjU7EMcxWwYevM+DcENdqQfz7 1A35LkomDNhXljW58/nYQVmyNRPGXr1tB4lfXD13yGnFQ9gyTN2UuPPG2K3D 6/gQb0wtd4gTT3/qtVJRz8IKAz0L5n8sSGS9uuLvlgWGgtBsL7GxTYlLy8Ms HKgWknlNXDYc6Wsj+Aibzr22TiO+MXA4MfbjI7B04/RPEj9N3Zk2uDkbA4HJ 7fuJhyw2Pdqtl41rWivntxHr9E2UUJ9mozmDV2z6BwuC9IQ2jn8O3k06tUYS p3TOrzonk4ddJwLZJcsshGQc6JE/m4ey8PyTFGJn9z/yl4PzsHiKz8CPWHPz tGnOWB6G3hoY6hEPGI8+WMp9jPUFBZbd31mQ7Og8kqGcD8/WpOShbyzwp29c 722dDznR+9kviGfdzgwZROcj+f5ay3TiFxvfhM2z87GuwvTtRWJroxddJyoK 0Kr/XxltiYXMtqcuTI1CxFGjS6r+ZeFW2rQ21aUQCva13nHEnq6KIpS0Qtxe zBx3JZZZc/ajFE8R5k0evZEhZme3BTdfKoJGP39q2CILMeONHTt2FQMHt/tp fWWh267M8eWTEuzeO63tNs/CAj3pcBm1BLOtgl7axKKnA1dn9ZVglUex+Wbi 3w4YZgavKwXzQOLeEg4Lgz8men5xKkXurRNbZtksfE7crl+vWIYO+a5OizkW eJqDfq7NK0cl3XTn0gypT/Pql4L6cgRpMYzriY8XGVFTe8rx65pQpQjisARx K7+1FWAqjS2IEK+zqUjQuFqBF/UlzYrTLGz9PrmqRv4Z1iZvu37iMwtKe00+ PcupBOttj4juOAtSB/4RKqyvhOse19TvYyyIHPr5WHZfJda4BzeVEi9pLqfG bayCwwPd1dLEr87k/eruUwVngxqpr6Mkf5cVLYonqsHMWvEhepjklV9cnDpe g41NuVcC+0k+ikL+wbvqIJbbpWvQSfpLLVu4TrMOe4aslmc7yPrjR3OXzOpw rDPdNpF4r8XlTvewOsgxTH4ZbWfB/k6JrPVEHfiVGEd92liY9lGwTTv3AgKq D3ZGtbBwsfcNv+/Wenho8S+GNrJw7IGAqWpSA0QDDO8EVbLQMZPw1iynAU6T XxUFia01ZE/4lDWguTDOKf4ZCwE96hrVbQ04ylv7PruC9Jewi5w2LxV3ou0N G8tY0I5oXzrlTEV5cff72SJi97iHNlovEaIkYLMyj/i4xJeYkUbM3PTN3kRh gXFaVjpt+xs8Gy+8k+lB+umSvtChHW/xh31DeLsBC6rRK0SqrdohqhJ46MI2 FmbsXFzsDnbi+MtbMYcZTIhUmG7RKX0HhZ/OykWXMGEuHNhIMelG9obiv6w8 mTizALpK93u8TF1R1arCxMXxExwdHxq8j7gEzU/OQCeyyTRfqAfiG2TFI1Nm 0FfcI/JbTC+OCjP8VxrP4O6EyqDfgX6I83+zv7g4DXamhCCj4iPufLlTfyJz Gh5Opseafh+EMN3ntr32NGw8/S39bwxhfZ/45dq+KYxxUls0NYcheaVZrPva FL4sDWmqbRjBT1qUZqt1UxjQcR4+2zmC8dM5+/TTP2M/QyUnOWoUnD30gByV zxASknKXsP4Er4ONkm+rJvH25jnnl9JjYCoHhdnoTELEo8DpEWcMk4WvvPbX M/Cq9pc7xfXjsMpXjL2uwYCJodTQUOAEaL0e7TrPJ+BoK6N2ypSBe40XlHKV J/D7c+c6usgkoi2W7v2UO45nbpyr0V2TqPquKi0lNY6M3DSx0MzPUAiZ1D8Q OQabnTGLyTZTWD571Uf3308olXvcP6M6jVENn4jZC5/Q7yTAvkxy2phOrZVs GIWbQI4Cb+UM1H32uQUqjiJK2cFuIJwJWqO338rgEdxKdH/Xr89C+U1m1cTo MASo7dLrNpD/ns8bhyJVh8E3EeKk6UW+Q+f7bb1Dh+C4Z+BM4PtZGE+Vi2vV DUK0XL9kRHMOg2Lq9jAdwD6rC4VmD+YQ8F3RIXmxHz8OhmsO/JjDovkpR5OM Pggbvz4TcoWNKMrXCAH5XsTffrKg/pqNbU/Hf0ui0hFwMbhrqomN6ZrVFy88 p8Nx/+9Fyc1seHhf2ij2jA4tmoDjt1Y2DPkXbsQ9poMh4/ixupuN61pjIZHR dBytln95ZIQNB5a8qa0VHQOMR9Fqy2xUqRWoflyioanGz3XiB1k/PLclZp6G krtmpxJ5OOBjTmbrMWkIV1216isvB4J6svSnwzTsC7D9q0KAg7hovmu+TTQE icpaHRTjYCSqrXZVAg3On5fUx37iIGRzW3r5PRrOP38nFr+dg1kjTRW7SBp2 2wS9m5fggC2edJV6g4bup8O6Zbs46EhI2HHDhobawCoZ2585kHlkGaJgSUPe uej/hBU4qEpjJ384T8PNJe0qzz0cZDuWf1Q2oMGhVSxBdi8HFk1La/qP02Ca PuvRpcyBlLzmaLgmDRqeTUbB+znINYpwU1Gj4We99D0qKhyk8zHyB/bR8D88 5T3Z "]]}}, {}, {}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->{True, True}, AxesLabel->{None, None}, AxesOrigin->{0, 0.36235777730212054`}, DisplayFunction->Identity, Frame->{{False, False}, {False, False}}, FrameLabel->{{None, None}, {None, None}}, FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}}, GridLines->{None, None}, GridLinesStyle->Directive[ GrayLevel[0.5, 0.4]], ImagePadding->All, Method->{ "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> AbsolutePointSize[6], "ScalingFunctions" -> None, "CoordinatesToolOptions" -> {"DisplayFunction" -> ({ (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[ Part[#, 1]], (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[ Part[#, 2]]}& ), "CopiedValueFunction" -> ({ (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[ Part[#, 1]], (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[ Part[#, 2]]}& )}}, PlotRange->{{0., 1.2}, {0.36235777730212054`, 1.000221551857682}}, PlotRangeClipping->True, PlotRangePadding->{{ Scaled[0.02], Scaled[0.02]}, { Scaled[0.05], Scaled[0.05]}}, Ticks->{Automatic, Automatic}]], "Output", CellChangeTimes->{3.619889431398744*^9, 3.7459421160582438`*^9}] }, Open ]], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{"this", " ", "is", " ", "the", " ", "error"}], " ", "*)"}]], "Input", CellChangeTimes->{{3.5884323984402046`*^9, 3.588432429188205*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"Plot", "[", RowBox[{ RowBox[{"Abs", "[", RowBox[{ RowBox[{"f", "[", "x", "]"}], "-", "poly"}], "]"}], ",", RowBox[{"{", RowBox[{"x", ",", "0.0", ",", "1.2"}], "}"}]}], "]"}], " "}]], "Input", CellChangeTimes->{{3.5884323656187735`*^9, 3.588432381752536*^9}}], Cell[BoxData[ GraphicsBox[{{{}, {}, {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[ 1.], LineBox[CompressedData[" 1:eJwd2Hk8lU8XAPCbZE9CSCX7zrXdLHEPUYlIsmW7RJJdKSU/VLZEiYrsJFsl IaFIRSTKVlH2JBUltMnynnn/qb6fmTlz5sx5nkfEDvhbHmSiUCgF+Af5e3ZI gFppak/Hfzbc2rSSnnl2syjXZi2gUP7dGT4zQveXuh/7VsQYKBGDLmmhP+h2 v5j06kTscLzZzDRgBbQPygqWiXgCRZ8t9Z3tWmDt3hpZLHICKIn2fiYyYiBp 1qRTIBILlGmZ7q3NKqDOwXE1XyQV5/91NOTUh9TYkeu6MWjZbZ7HTPRhgaXm bs93NI9Q5ZFz+tC40rOd+fE1oDgWMV6xGIDVYhPzQbd0oPTob+dbNoCg6dNH pYqygMLWQXv70RAcvobKdqfmAWVSyWwwfidkdM6tshIvBkpAZZHldnOIzFXr k92H1rBiTQ8yB5/AgFuLZ9HTla01181Bb+2kZeHHYnK+RyEr9sCgxYfs+aIS oKT+FOd6sAc2d3Rq51JvAcXS11Uv2QLyXpb6ftO9A5Sl59tOfd8LcVmT+k99 0XHjMleYLeGonzx/ahb61iq99PWWIM66b0CUUgYUqxsCPoaWMFPw8myrKzr0 0gGxq5aQNN7UsUnqLlAyPqmOaO6DnoOVXo03y4Ei9G71pI8V/OpN06p8in7W 6a4XbgWCu0+z5L9Hxxtmh16yAns187yznBVA6Qp933jPCoaWPvVt80Zvs976 36IVfEnZsOuxfCVQGv6NJZ+zBkrrGZn64ntA+XJmxfnLNiCud/jn7cfod2Zl Bfk2YFi252lmH3rGV+ZupQ1EX93I+I+9Cig6P+OSemyA063qqu5hdH3HZAyf LQgsfF71QPY+UB5Kbze8YAtKynvH7hdW4/kVLJWC7UBU7SHvncfoce1B52g7 4N8iY1DwHu3+kfrfFTuY11vMvLymBuvhsjOiwg6emRXbBgaj/2To5X6zAxff FS/kd9Zi/TR+W7vuh6Rbd+9mjj8AinxQY7OWPUSVbRy+THkIlJILC3nb7eFE ZQx3vDDakbHWx9IenB86eYeYoynJs83e9qDQxi5lW4Vu+B5XmGkPjV9dU3mi 64Birxu1eskefsnzhp+VegQU42umZXcdwEizgO+RHprfTSnvgQMkGeoUzVuj s85To5ocQNnxQGdgNFrZqEykzwEOJZRLuHxCi64/kLPsAJPBcu5ZVg1AYdI6 ObbLEZzfPWc7JfAYKCE/H2V1O0LUMF93qjh6yb398XtHuDXulFmljL7M0/3m gyPMz/xQndmBtmgYeTPrCCmcwg6Hg9Hvqp5R+ZygQ9frtm0vesEl1GSPExhk c1hqpD3B+hnX7a3D8amrbdaFaE+Tz7sancBFV2JncCV6+E+S2gsniOjbqlv7 Eu0TEPG61wka+Hyl9Vc+xfN/Kzkw6wT6sa/mTX3QIu6BuVLO0PHG4ZjvSbTk c6sbCs7gIjXx/UI0est+7TRVZ4h4QvnYmYNWTpe21XOGhgXVV7av0RyFx0Ss nEE/8PJ1N3ojUIRXUL6EYbxHoiJRpugZWtOKSIzHfTu1wA5dbbya8xzGu9mU 8PkImm3LjpkkjPfxV7B/ITp0tS5/Acaz37/7FE8Tvk8eat1+jvGKPj7L2IT2 Cpsabsd4vwMN6uXRHSw3WLow3uW4LSu2o0OldGXfYbyXD0RjTqJbrCcaP2M8 Q5GfSaPogHBRbxYGeNi+FAv+jvZz0KlgZ0CCd5i5wwJapXfvNy4GvEseLJQQ eIbvE/MUPT4GHBvLsq/chY7mMTXYzICMP+YxqTboxAFRcXEGPOFargh1Qy9Y U39JMoCHxli9/T90fSnfCXkGlERtbnhdhlb2bv6gwYCOa68ma+rQOhwqDE0G /L4dvj6rFc38k/OlNgOM3gwd8RhDhwzqRNMZMCSdI/VbsBn754++604GsGy1 sHwvie75s+bkLgYo7qGEP1JFM52KjzRlQEiwS2+MKbqScvP4Hgbwt4ieFwpH 510xbrFhgE5/x/2F82gumnmYHQNcpyPGhlPRAeJbZe0ZUCo0oldcji69/cnW iQGvFS95JTxCl5TcG3BmwIK+QUpgG/pdRYytCwN2Hc6d1h5HT1+1knFjQMB/ ezeJzKLLnjj9586Aq5dWmDBRWrBfs2KaDzJgrMb1eut6tPHv/VqeDOB4ubaj VBrNU8PDOMwA1dHHC0nqaLZK5hAvBoRxiNs4mKHz2hrifRiQL9J1BuyJxVIi fRnwQu3MHYlD6FSR5gA/Bgg5jLJNRqDnfCykAxgA/km0jgR06MT972iPs9sO VKahdVyu3A5kQMXNvAehleiQ6JcsR/H+H1lOuDxG9+/VykNTepjWbX+JNtqn oRbEAJmJcgO592jmK8+r0OYLB/xWT6AbNi9TjzEgiIcv/ccc2W/+fSY6XfJp 8+sVz/F7In+Ichz7RevoXA03uqCi0BY9sVtCLGsDujcl/zpa9+OniHwZ9NV2 92F04n+3hkvU0ak1bnzBWL91gfp3Ac3NVqiL1iql5dw3RY+6aTii43fML9fZ oruOrz2CHh6sZzS6oXd8NgxDawSffdTqj/Y50BmBPrfGeHPnKbRu7KOT6MFC rvC3MWiXZb7DaDX9zsGBZLSRd485Orr3Cn0sGy3ix6KAfhdgn/XlJrq/6M4i 5qvMvnlp+j5a8sHrZ+gzuR+cfj9FBxiGRaPfaBfVLb5CTwyW6qLlu3w2Mfej Q+R8P2O9wrxU/+OYIPWoq49HdzH96ueZI+tpxdJo6fRaXUFKK1A8ZNXuY/1P qYdnbOJCVysx6OiOF4YLEkJotve0h3h/ku5sjvKS6JbGByrotsuXNmjqooWu fV3E+xdVsjmlZ4zWLcixQQc1Cb83tELz0DYUYL9s/HU9zcIHLb8mWM4f+/uC 57zNCbTAUTcn7LdGaSV7p0j02DmFWOxHX9uq9V4Z6MpZ9yfeeN/fQ04GFKEt rEc6sZ8FYqHveCXZP2ZPL/Z7fXVL6tk2dIaXYOshBnALvxfMXCD52H03wOfp QHl28HW2F0A543WZ/QADqkzc3xbzoyNMbVvw+WOcmrpapYje4bpbAZ/Psv7l dR2OaMrX+Bh8vpmPNR5744lmOl/DZssAu9XnXvcHoYOKqKetsb/pvFc+x6Pt NnDbWTJgb7Yk/8o69HSvQha+T2Zcd/Fu2dQGlCf2BYm6DLizzUJcXx49sfuB ow4DfCRs1Uw00ZN3f4toMWB8zN3SeS+xLFO0OvbHoYik6Ci0SfXGITmsh+99 3rdTaG5J+6f8DAg3rxcfmUcPXxb4uRb7n9qk9pW1Hb8nTIc2rcHzTndZLouh l+4P7sX3dcnRb0kyNui40I8mi86QfFKK70Q9+tdG+bMfnMHCXlHizAs05dEu sWFn4N6qrh7fi+5Q9qzqd4bYBf19OTNo/ukHTa+dITTcMblF+iV+/0eDTzY7 g3tUMt/6i2iK2sqtxc6gcWkFfy3jFX4vGoWtPJxB5NdG23Fv9AXDAfkDzsDu oJXGewK9ZU3ovJMzDEr6i3onok+4uUVaO8O56n7FjY/R+WlehkbO8P37e3FY 24Hnf/vhiKgzGAne35p7BC3/6/SDbieoeW7kdP8EWvpV+/2XTqAc2hXWHoae o6Xdfu4E60emHv+NQ+sOtkQ9coJvJZLGltfRrcUcn0ucIJWeZLWyB13qGxke 4QRTB319D9I6sT9rlNfIYLx/XBOOW9FPTgUeEHMCv8RbB6wM0Pau3GUbnGC6 9qutoRn624rDOjxOMLPGa5vYQfSgnjjXX0f4VeMhOHAZnTF0o/C5I57D5bHV HDqJWjfo4Qj8VZbrjCq6gOJsnb0zyQHs+E43pe7twf6ltCue2w9mv6BXvec1 /nwJUoaetuA8vnPWKPgtfg9OTyZPWoHRuRbLW7x9eP8lowkulvD+bh+/fdI7 fF/dfFL8dQ9c+KQ+9J9aP1AaG0Piy01gJk+Ea6JqAH/+TLTI/LYdjnhbGrQ4 DOH32MOldF4f3I6GO4WHDgMlp3LVfKcWfJzNfKGnN4Lf73h38x8q8HN+WE+T exSfJxU+lwBx2CjrX2a4Fq3S2iWnIw5G1oviFvxol3FmbmZxuHxnPdthYTRl 5INUqhiouVt2pUmjh0+dvvxYFALan3os0cl4UuRnNRH4kl1wsckfnWN1+Het IPC2azB1HSH7iVoaWwuC9vyToMFjZP2USsd3AYi1Gtr/+xSZP7U6SkoAZNiF JOVi0foy3yOT+cH9yLnq+Bx0g7/Ll5NrYdDIZ2RfJ4mXUuE6wgLebmcpe3vI +L7xMhkW+HM6bbP5W7K+X8HXbxXwPWpxMh4g+TFZKi2uBGNtqXe6n4nlql5u XgEVSoNdUis+oFekW8/+ouvv/vlDnJl49MethJ/0l15ca0VZ0REvD47IztEn CnX2CK8m41n7R11/0EXEU15wryfOKEp6+Il+U7/0C+dGYt0fG5fH6FqMJnb2 zcTPtd0NR+n7MmZ3rpQiTttpk9dHH67lOESRJeYdsaH10P36xKIX5Ykpe4pS 2un/fmvdmFci5o/JfPGEHitg0fhbhTjn6KbVZfR1tEMf5tSJIxwFHl6h5+0L Y5rZQlx2SmGglq5y5IrYd23ise4ikRZ6XeIt/UldYmO3AbdOusmdp4zPQOwj +7HyDb23/V3Y+DbifkFvpn76DAd73bAxsV9I8K4P9HA50f4BU+IdXT+/fKRz GWv+e2f+//EV5XETdJmog9qvrYiP3x5tnKI7Lj++99zl//U8xM07R1edUC9M T0Q3cD85LMYMvakO5X3JxNeyotqYIWLX2TrBFFJvOtuG4FXQUdLVnZxJ/Fu5 vJUF/HwDKfElJN507EVvdiiZKd0f2ojWD/s7l8wN4otyXE5/0KI0lSDgh9bb ewXT/5H1e01i4/kh0PmkeN8SWX+Sxa6PHxoetWhZrxrD9Xo75I+sA+czhw6a 8aJdHqy9liMAaaw36vUU0TlSydemhYCXVzRQxAUdYVn8x3ETNH3gYpdyQze8 6RjN2AQn7v3NUfAg663qE/s3waBdd4eWDxnf917TQQSKs6NV9p0g4013y603 g77S1PeYRDRl7+ub2aLgb1zr96OB7Hc9Jv+IOIgLF7L8eUr263rUdUEcer4m Zy09Qye63xwtEQedi74vOdvRFt9TC0fFgeWNqLJ0HzpA4l6ThQRkucVM2f9A lx3fZSkjCW1hVj6NYh9xXMWyuE4Kro3uneaUQufsvinQKwUeO/YE7ZNF87g7 usxI4fvOJHRUGU3hao6VlgZaJsRTdNGJfo1KCdKQ9UD+lq4NWtSpqmWfDAT8 oUxWxqE7Pp3Z+UwW9ByX/BcS0A09s/r9ssDR8G/W8BJ6OO+NwIws5Mf8+teV grYQ/GO7SQ76BCY5f+SjI7b3Uo/IwTbaWwWlR2iVr7ma6+SB/8ht7xuzxEsm VfoKsLNZPWPgF9pFRnR6rwKEbKxtWzePLvs1x+qmAEPPnilFU8Zx/H5RZ6QC lGwY/u6xGk0xPdf+TAH0m3iDZGXQZWKvDxspwlHha/ku8mieXPaofYpQ4L/5 daoSWlTDP+yAInAKK25h1xgnvy96zxmhCG/8tv/5rE/m70u7VKMIvkInQm/a Ewt9MZVWghzfpVsfnNANNqOm6krQ9SRyYIMretp7p7y+Emj5XqLHHyL5HPU8 s18JmJ+UUHyD0BbXQ3LilCDNeyBS+QLx3YKUcSVob3C753GJ7M/hE/FDCX8u +/Ix6zLZb/vYzgUlONjwa8eadPTwDsrRtcqgso6HfbqQeL2vjo4yPKvflnD3 MbqDEZAUpQzMtAsm5Y1k/G9UZ4IybLvZx1rRTOJFTMxeUYZHKQFnKttJflU3 Rm8oQ3VAVvD9PnSEQahnozL8Gv+sUd2Pzvkv6MkLZdBwos1UD5H8tAeWu5Sh bFebT+1HMj+eV2pEGYrF513rfhC/urB+QRkmUrdvrp8j65kTLjJRQXrNpf76 3+Q+tMQ+s1Eh75+MbcMiye+byM51VEjrsd79lP0Txj/xp1SRCr0mueyNXJ/I 94a6TY0KAo8nnzWuQfNck6jRpELy7bMGz9ahy3LsLLdRIT767pbnYmh9/Ug1 Gyq0LizMPpdEW9TTxuypwHbU+G6rDHr4V/tJBhUiGUMKbUrojs5mC08q/Ke1 WuyVNlrUP7j0BBUeltoNvtJFU9rnwkOpMC+Zn94BJP7mSFoEFY6v3SrQtZ2M v/q2O4YKfl88OV/vJd52SekyFW65VLa8tkKr/NIRTaHClzfLUW9s0S6PJH5f o8Khp1cpvU7ogLLne7OpUKA9WtfrQsZPJb3OpcLYHaVTfW7oiNy7evlUcMlo /PnuMMkvleNBERWyeHkq3vuQ8xlbdJRQoT/WIaDfH93Ay9x6iwrCywVKA0dI vvMKRaVUsDs282XgGHqa7aVXGRWuftUrGjxB8t07vbacCj2u5w4OnSL5fE3N rKDCXnPR4eHTJP71KweqqFBiE3pOPIrUT1vx3n0qMDF61dxjyfnSEj9VU8H+ kEb/jfNk/4etf2uoUO6fGPXpAtkvf/RrLRU4TkwqyyURM4brHlDBLcK41+sK OQ9Lm/9DrHds/ulbqWS/mtsr6qjAf4mi8C2dxCtIOIr2uebYQ80m+60KakI3 5Vb/F5hH6iNw+DdapIRfpuIGsfQx9nq8n/KAjrkicp7zOf9w/GVt28ktt0j8 gtl2tMxTWYkTd8j5R49HoCNeRLbVlJN6VarxoXu7h4/9u0f6L4Majfmp9Otu 1qsh/SJ1vA/zPzeW2hL2kPTjGPdq9OjkXGDDI5LvcSYxPK/OT4sNTE/JeaP3 r8N6JC3eajR8Rvq7ZeMXrN+XVex+Uc9J//61ysb6buM+KNjcRvIf5aDdo0K6 wOMGtg5yXpudN/E+5kQ2eZl0k/WqfEx4X7tlTvLFvyH94hisjfe5oKnqsWaA nC9Q0QTv30o/Yc3eYVJf54sS2B+3jT9XJ30g58kuHSikgvP+XE6BL+T8KTvm 8qhQ5bpYaTtF1ud1WuRQYY3Xfudr0yT+Lb0LmVRoCFl7d+Nvsn/fwE3sZ6Gz vvbO82R+vmQ89nvA+ecrcxZJvnWBZpeoIJZx2kaCeQLPn2Z0NI4KIfkDy+6s 6ID9A2+jqdB1S7u4gAMdYZ+36SwVztT9+Ce3ljj27e5T2P+DrrkqIhOk3xd1 vKig+6nO5IgY2kWSteQgFa58Xz9XIYlWaYmbd6HCDkrXDk0FEn+lrL4tFQrF tk3qaaF5jLqc8f3g5S6uaWpJ4lk3z/BiPyQY+lpbo3O0Wae4qKB23/06w47M f721nYUKf9kL1wQ5k3gxN2BeGWLvKExkeKETP/pMDSlDwbzGtW9n0A1wLKhA GTgkbV79iUIP3xToylIGP7PgVSvPkXwLtQVSlIGWUxMoeJH4oqxljDI0bqeb 6KejKbeLr3kow2jizn+XKsh5Esz9RZVhR62nakbVBPn9qp22oDKUfDh3qKCG 5MN58stqZQjc0tZd+4iMC4kz/cPvRb/FrQ8v0B0L4mw9SiAiY+9EGyP1P+kS fFoJHOp8HvXyf8Z8Gc6yrYog2SuQ1y2IDnijKV2vCN9mGiJfCqMTvfetLleE CNl1Jo2i6OnQEynX8HuZXNdzR5E4KTHjkCLMHOb+Em2E5pF4cmRZAc4JlPHT gtAqfmNXRRTAUtX+NzWYxHt4q2ONAmzczfxOPgQ93ET5SVGA0tN22aIRZJwr eHFUHnomKfJc8WjKqhfcBfIg+nQvfLiO1mfP0JWVh/v+s4cvdaMtaAW8+PNE 2sZSrfY3ZPzta00OOQh77snK/o7Ed/Xc/lsWdkgM5kcMo8vO/1jfKQs9b1qG /KfIeR/RSs7KwoxelrU56xfst6WyoDEZUOTaZcC5Fa2y68HGK9KQW5QlFJWH Hq7lV6uSgMTjIkHnbqBzRu6cz5aACKOsVwlFZP6aLT2xEuA8nBl9tRTtYrhV zUEChIUy5wpr0dPfyiSWxSEpJr2jtYusF8ueMBSHMx6psTwrv2L+6bnO+P8t N8mkP+luxNZHDopsAJ9Icf0RD7TLSNfimDAcHyuPkfZC5+xvOH5TGGLzuwXu BqCHy0Y4tIThtuQ62rP/0A3J2Q171sNvydTA6RR0hJ5l1ElBOC+V9cWojcTL q7p0nw8uRymrxb0invGMC+GDrI/1J191kfn9LZ56fFB2Y4jNvo/sF6Dz9gkv vJYSk/EfJ+OaMeFta2GzdL7bNcok7uemuqNqDVRIl/RPaRAbvGlpZ4WWwuWy f5pofXuNnTas0C9rHcW+ldhDKHeQBVgUKMpSBsQWlye+rQI7qk24oxlZr7jn 5WpmWNJkEn9xEB3x4rnGdgrw1tj+6vVEU2jHN6Qs02V0breOe5PxwwkJ15fo Frp2R5mOkHGzm5UPFujX9UsbtcOIOyZFPv+hG++yP1R0lcw/enJc6wfdsfXO 1qprxJ3FIwbT9EDTVTyNGWT+o/7dV77R08zKqofyiFulRe5/oU/tZeEQvEMs //ndz1E6U7fDoGQ5cXvveeYRuoDV3XK1e8QFM+VSg3R9G0cH8wfE1WrZjm/o Vm/vUh3riRUpP4930Q/bsTJ7PSa+nLpBrJ0e1ufYG9z4/3wpvZZN9CT78ltR zcTDwsYsD+gF71lPJ7cSi1rythTQHzg6Wee2E0dc4mK+RO8YKJe700EcK9w9 UUEfc2ZbethNHH+69Plj+t8hp67WN8QaSaY5z+ncrhUFvX3E9tObHnbQNd2c 98wNEfvZr7jaS989ViHB9OH/8Ws/ne2nuxxk/7NmnJh1fov2MD3uUGWOwiTx Id4i2TF6vzfHzIE/xL77+Dgn6TaDE0dc/xG/X3FqforeYdE8y1gi/swj/eU7 vYkW+dOReQrrX9wl+vEHvXR58Y8ND7HhdNWzn/TTyT+WTWSnSD2Tdl1cpEs/ 6F2tbjdFfh7QmE5lAX+OQrmV1cStTQ+l+SCR6nlwMOYbvk+G9tGYxCEuJbCr f9d3fF8cyg8/IQdGe7ZUS5mhVd4+80mRgyWWf5l+FmRc5bBplRwcOR55mGKL nm4u75+VAwfrFIqkOzpn9p+pnzwo8tUpe4ehy144b3NRgFcX2M7/LSfjGmcr NZWA4+krMU7uaaD8EQ4ctVGBtJRfXwrXoC22ana4qYC8z6ZKo7VoISNqRYAK mKzz3hHGj2ZTDLaLU4Fzh1i8fwijE/dVitSpABuXXuVbGXQRr3G8qCqwWpfs yDdAx2ZohQ2qwlW5zjUGhiT+6cXZz6ogvfS7d8AIzWO95PhTFXYUbvcWNEZX CwexcalB9N/hi3HmZH5C4FMtNWDJEuwLcECfGKt4cFENmD9FeusdQ4vGpRup q4Pwbu5T9cfRnscrHHTVgXo3JQ5OkPm23w9uVwf7UyXFBqfQxu8y99mqw901 HZ+2n0anhrgmhaiDs9YGd/MEtKzx20v16hCUmX/01QUyP/J7zzN1OMekfNYi Ea3FacPxSh0q2/TzLJPJ/uq9pkPqwOHqMWxzjdRnSlN2WR3un7vryLiBdom7 rKOrAW3fdHyGC9AtHfQMQw0Y3dd4yrXo/+dlfDfRAG6Rt2luN9EBu18c3q8B 7uULvYfuoicye24f0wCe9zutA+rQEVuqPG5qgLR+p/uPepJvrn3fXQ3YesM+ 6EgDOkc1TbdaAzz8fZODnqIlyyt6GjWgbmVy54nnZFzu5a4BDfBSHDQ7/Rrt bhy7lZMGCe/8u8Xfkvgy4dJraVAWu2J/Yy96cqXtkgANfo5JurP2o/tLDwSJ 0yAi0zskYRQd+eF3sBYN8kwXl5XH0GPvmFv0aND090JUx0d0w2F3ZkMacNiU J/J9JvfR221gRoMr3H8Lrn0n+zGFt7vQoPphnOLWH2jdNwphB2nwzmtjef8M ib+QLuxFA9FmqNv8C73bMUzgKA1KwqO7byyQ/vLcJxRJg3Ylof07l0j8CXGp WBpMvy8enFhGG6VKbYqnAU2r/bPCyh943xqhL5Np0DDNRylnR/faH9XIo8GH rBtR+zjRHT4jqTdowGKmyfmTC+1y58enIhqYFtsLaPGgc1aVbb9DAz/byfS+ tejdZdpW5TRIXBUmdooP3Viz0eweDd645irWC6A3/rixWEuDv2vUyxlCZL80 mXt1NNhY36i5Qhh9K63GpoEGrsIT24w2oWP9mmybaDDZpby/VQLN5nFqz0sa rDndMOgtRfJ5G2bdQQM1qqX7ahm07pdloy4aBJ8/5m8hjx7erzr8mgZp2qw/ fyige4Yfnn9Lg4efUkOSldDVTqGifTRgMnoY9UaFxPM0W9VPA8kZM84TamjH VTrWAzTYmTOUuF4DHVrFc2GQBl7mgQIPaOQ8/+rvDNEgfoEpw1GT1MtcrXaY BndKLostaaFF5e1uj9Cgy066MFsHLcskcm6UBnMs1YoGuujLxoFmH2ggeG9X +ageyccb5tENs3phnkDyHTeOHqOBp5qayTd9dPyqikU0T6C0QNA2Uq/rfi4f sV/uCI/+NUT3p/mUol2+cZeGb0c/fJz5Cc2utDJk1U7iyRnOcRqUe//eft6Y 3N+U4wa0fcnXtWtNSLzktnVops9DA1dN0ZE31f7h+psyPcUbzdARy+Ev0FYe LcfyzNEqtjei0Qv5Dw1kLUg/BScrovM/lK0u3UvOJ7G1DvM1E7/Rp74PzbUh QQv9y+XajRorMr4hJAvPm5WdEAg2aOPhmW9Ynx2Dp/WabEm+YhNy6O8bj7Ob 7kdXXjHZg/VMdfB63WGPnvjOwcB666c559o4kvsfWb8f72Oi19K334nUW8hT D+8rUXCn9gEGOZ/+OCfe58hlaqfvAXQ+O+3QOxrEdUtkzrqR+v4cn+3FfuIV OnzyIOnfw+e9sF/OXKSsiPJEB3j4reuhgfzLuTYOL5Kv9IVd2G/dXJ9TE71J PQYOHMJ+lIzrVMnwI/d5dorRRoO2lqYFsQBS76uROq00OMZa21IYSPwpaamZ Bs8i81wqgsjz4Pd221N8np5eVdQ+jm6o/PYEnwcBpvN/6oNJ/ragWE+DQ+FH E1tD0FbV5rXV2A/1hxwtQkk9Ar8N4vNWveAg++Y/dFv4uUl8HtlCjBpGIki+ x44+uYXvi6B103+iyXkM3GyzaWBZwf4wLBb9R17+SToN/v1YjGGOQ+uP+K9L pcFu//HNPAmknw20jyXSYOrw/T0yyehJcW6hCBpcLbq54fZlUv+365pO0YD+ KfuT2lXyPlDfbB9MgwvuseH0a2gtsT9b/Wigwth/xzqb1IeZc8CeBucFkjq7 c9D8u7aMWtNgvL11dm8euuxZ3ysLGmTq6mqZ3yD9wmzuvgPfn8KbH++4iR5T NA1QocFYz1iX5n100NcEVgoNIH7Tz3vVpH9eGiT81YA0QxtBjVq0hWH4/IwG 7C1vdlCpI+efsD//UQPqL5aMyT4l9RqcjnyuASkmAb/Wv0JniPq1xmvAzIpi odQO8v4JvvM4UgPMa0Z0BLrI8xZ2LjtUA5hl94Xzvib9tOkFt48GBKzawsb5 Hp0Y0l+zSwNMGv6tXxgn+aR/16BowBItRm9wCa1RLGm1Xx08trS8DqXM4Pdp KoTdTB1ebmH328CEtnpanKWvDtma57PsVqH7T2zzlFEHfe2Ly12caJ/rWsdn 1eCsbsrjZiG0sfCsd5QaTOr27vcQRlM4L68LVgMrvfUzzBvRzKc7Mz3VQJKe Lr5tM7o3qWuXqRo8g+yzD6XQ+U8bVq9RAzbDwu131dAPTw6pJ6pCgOHEwB4N sv/OubpwVegzlDv+jYY2yvWQ8leFEqObBQra6JZOvngzVTDdcYe1ANAZ8c5V bKpwYdf952m7yXl4N/ueUAE+i2azSA90kRzrrZPKwDTipcjqiY7NGU61VIbp QG7OuMNooT+rtisoQ1uydctFH3REw51375Ug8u2HbelH0NM8Q7+1lWDOmaJV HobOuXyt/7MC9Phqiw9fRasMs3oIysGTpX6KWypZn//+3ydZKLsYMfTxGjpg 6eXualmIv9uc8TWD7NeVsNlOFox+Wgv+ziPxCh5GXpaBytAjXNx3SD2nbFQo 0nDl/M1fW5vR1bUHN+eJQ4Tz2e6vLeiGNpOdbuLgpepQlt5K4ndz80uKA/0t h9e/drJ/+rH0fDGYEPcaqO1Baz0555grCjq1so3aozPk94vyC7ybYHDixiXN RbSF6skhKX5oefCf36cldMfch9uJfFB+wdo0hTKL/fWGt3aeF2I0Vq36vRJt EakTh/+/U4lwP1HFgXZp3brUzg1nBCUYNCF0w6uhtilm8Pkyv/XjeuJdx5pL VoJNXZfQlQ3EeUGqHkyg4Hama04ErS+3O7dsmd5TOrK9UgpNsVlltu8nvf50 jbi7DDpiXpG7coZebHVpmU+OrOf+k7t2mh42r19zVJHMf1ZzPPoj3bNd6KqE MnFpxynlYbplzvSRbirxpS/qOa/pukdb9pxVJY501brdRJfZkaOork5MiXCV /6j3P3+tstM= "]]}}, {}, {}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->{True, True}, AxesLabel->{None, None}, AxesOrigin->{0, 0}, DisplayFunction->Identity, Frame->{{False, False}, {False, False}}, FrameLabel->{{None, None}, {None, None}}, FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}}, GridLines->{None, None}, GridLinesStyle->Directive[ GrayLevel[0.5, 0.4]], ImagePadding->All, Method->{ "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> AbsolutePointSize[6], "ScalingFunctions" -> None, "CoordinatesToolOptions" -> {"DisplayFunction" -> ({ (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[ Part[#, 1]], (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[ Part[#, 2]]}& ), "CopiedValueFunction" -> ({ (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[ Part[#, 1]], (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[ Part[#, 2]]}& )}}, PlotRange->{{0., 1.2}, {0., 0.0019205816256888664`}}, PlotRangeClipping->True, PlotRangePadding->{{ Scaled[0.02], Scaled[0.02]}, { Scaled[0.05], Scaled[0.05]}}, Ticks->{Automatic, Automatic}]], "Output", CellChangeTimes->{3.619889438348675*^9, 3.7459421185653176`*^9}] }, Open ]], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{ RowBox[{"Exercise", " ", RowBox[{"(", "practice", ")"}]}], ":", " ", RowBox[{ "repeat", " ", "construction", " ", "of", " ", "the", " ", "polynomial", " ", "for", " ", "6", " ", "points", " ", RowBox[{"(", RowBox[{"degree", " ", "5"}], ")"}], " ", "and", " ", "compare"}]}], " ", ";", " ", RowBox[{"also", " ", "compare", " ", "errors"}]}], " ", "*)"}]], "Input", CellChangeTimes->{{3.588432844659224*^9, 3.5884329388739357`*^9}}], Cell[BoxData[ RowBox[{"(*", RowBox[{"The", " ", "polynomial", " ", "obtained", " ", "with", " ", RowBox[{ StyleBox["Mathematica", FontSlant->"Italic"], "'"}], "s", " ", StyleBox["InterpolatingPolynomial", FontFamily->"Courier New"], " ", "procedure", " ", "is", " ", "the", " ", "nested", " ", "form", " ", "of", " ", "the", " ", "Newton", " ", RowBox[{"polynomial", ".", " ", "Below"}], " ", "it", " ", "is", " ", "shown", " ", "that", " ", "our", " ", "constructed", " ", "polynomial", " ", "coincides", " ", "with", " ", "one", " ", "calculated", " ", "by", " ", "Mathematica"}], "*)"}]], "Input", CellChangeTimes->{{3.5884335873454294`*^9, 3.5884336088853054`*^9}, { 3.5884337133863897`*^9, 3.58843375243436*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"mpoly", "=", RowBox[{"InterpolatingPolynomial", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0.0", ",", RowBox[{"Cos", "[", "0.0", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"0.6", ",", RowBox[{"Cos", "[", "0.6", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"0.9", ",", RowBox[{"Cos", "[", "0.9", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"1.2", ",", RowBox[{"Cos", "[", "1.2", "]"}]}], "}"}]}], "}"}], ",", "x"}], "]"}]}]], "Input", CellChangeTimes->{{3.5884334021382203`*^9, 3.5884334929420204`*^9}}], Cell[BoxData[ RowBox[{"0.3623577544766736`", "\[VeryThinSpace]", "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "1.2`"}], "+", "x"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "0.5313685379361054`"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "0.4004353824203932`"}], "+", RowBox[{"0.10217161531292468`", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "0.6`"}], "+", "x"}], ")"}]}]}], ")"}], " ", RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "+", "x"}], ")"}]}]}], ")"}]}]}]], "Output", CellChangeTimes->{3.619889471918339*^9, 3.7459421232214785`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Simplify", "[", RowBox[{"poly", "-", "mpoly"}], "]"}]], "Input", CellChangeTimes->{{3.5884336304216337`*^9, 3.588433633259918*^9}, { 3.5884337014511967`*^9, 3.5884337068577375`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"-", "2.220446049250313`*^-16"}], "+", RowBox[{"1.4432899320127035`*^-15", " ", "x"}], "-", RowBox[{"2.3314683517128287`*^-15", " ", SuperscriptBox["x", "2"]}], "+", RowBox[{"6.106226635438361`*^-16", " ", SuperscriptBox["x", "3"]}]}]], "Output", CellChangeTimes->{3.6198894824082336`*^9, 3.745942124188427*^9}] }, Open ]], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{ "Notice", " ", "that", " ", "the", " ", "coefficients", " ", "of", " ", "the", " ", "difference", " ", "are", " ", "zero", " ", "to", " ", "working", " ", "precision"}], ",", " ", RowBox[{"which", " ", "is", " ", "15", " ", "decimals"}], ",", " ", RowBox[{ "thus", " ", "indeed", " ", "the", " ", "two", " ", "polynomials", " ", "coincide"}]}], " ", "*)"}]], "Input", CellChangeTimes->{{3.6198894970480876`*^9, 3.6198895843903723`*^9}}] }, WindowSize->{947, 810}, WindowMargins->{{Automatic, 159}, {39, Automatic}}, FrontEndVersion->"11.0 for Microsoft Windows (64-bit) (July 28, 2016)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[557, 20, 4670, 173, 181, "Text"], Cell[5230, 195, 805, 16, 69, "Input"], Cell[6038, 213, 693, 22, 50, "Input"], Cell[CellGroupData[{ Cell[6756, 239, 3190, 98, 226, "Input"], Cell[9949, 339, 39427, 669, 234, "Output"] }, Open ]], Cell[49391, 1011, 307, 6, 30, "Input"], Cell[CellGroupData[{ Cell[49723, 1021, 383, 11, 59, "Input"], Cell[50109, 1034, 1257, 47, 53, "Output"] }, Open ]], Cell[51381, 1084, 293, 6, 30, "Input"], Cell[CellGroupData[{ Cell[51699, 1094, 488, 13, 50, "Input"], Cell[52190, 1109, 10439, 184, 245, "Output"] }, Open ]], Cell[62644, 1296, 587, 11, 30, "Input"], Cell[63234, 1309, 1971, 53, 183, "Input"], Cell[65208, 1364, 398, 9, 30, "Input"], Cell[CellGroupData[{ Cell[65631, 1377, 299, 7, 30, "Input"], Cell[65933, 1386, 1971, 69, 107, "Output"] }, Open ]], Cell[67919, 1458, 680, 15, 69, "Input"], Cell[68602, 1475, 594, 14, 69, "Input"], Cell[CellGroupData[{ Cell[69221, 1493, 423, 10, 30, "Input"], Cell[69647, 1505, 347, 7, 30, "Output"] }, Open ]], Cell[70009, 1515, 276, 6, 30, "Input"], Cell[CellGroupData[{ Cell[70310, 1525, 846, 20, 50, "Input"], Cell[71159, 1547, 8785, 159, 236, "Output"] }, Open ]], Cell[79959, 1709, 182, 3, 30, "Input"], Cell[CellGroupData[{ Cell[80166, 1716, 325, 9, 30, "Input"], Cell[80494, 1727, 14136, 245, 229, "Output"] }, Open ]], Cell[94645, 1975, 535, 13, 50, "Input"], Cell[95183, 1990, 765, 14, 69, "Input"], Cell[CellGroupData[{ Cell[95973, 2008, 652, 19, 30, "Input"], Cell[96628, 2029, 682, 20, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[97347, 2054, 212, 4, 30, "Input"], Cell[97562, 2060, 363, 8, 63, "Output"] }, Open ]], Cell[97940, 2071, 524, 11, 50, "Input"] } ] *)