(* Content-type: application/mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 6.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 145, 7] NotebookDataLength[ 39829, 1003] NotebookOptionsPosition[ 37310, 915] NotebookOutlinePosition[ 37736, 932] CellTagsIndexPosition[ 37693, 929] WindowFrame->Normal ContainsDynamic->False*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell["\<\ Math 481A California State University, Northridge Horner's Method\ \>", "Section"], Cell["\[Copyright]2005 BE Shapiro", "Subsection"] }, Open ]], Cell[CellGroupData[{ Cell["Parabolic Fit", "Section"], Cell[BoxData[ RowBox[{ RowBox[{"ParabolicFit", "[", RowBox[{ "p_", ",", "q_", ",", "r_", ",", " ", "fp_", ",", "fq_", ",", "fr_"}], "]"}], ":=", " ", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{ "deltaf1", ",", "deltaf2", ",", "denominator", ",", "qmp", ",", "rmp", ",", " ", "a", ",", " ", "b", ",", " ", "c"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"c", "=", "fp"}], ";", "\[IndentingNewLine]", RowBox[{"qmp", "=", RowBox[{"q", "-", "p"}]}], ";", "\[IndentingNewLine]", RowBox[{"rmp", "=", RowBox[{"r", "-", "p"}]}], ";", "\[IndentingNewLine]", RowBox[{"deltaf1", "=", RowBox[{"fq", "-", "fp"}]}], ";", "\[IndentingNewLine]", RowBox[{"deltaf2", "=", RowBox[{"fr", "-", "fp"}]}], ";", "\[IndentingNewLine]", RowBox[{"denominator", "=", RowBox[{"qmp", "*", "rmp", "*", RowBox[{"(", RowBox[{"q", "-", "r"}], ")"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"a", "=", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"rmp", "*", "deltaf1"}], "-", RowBox[{"qmp", "*", "deltaf2"}]}], ")"}], "/", "denominator"}]}], ";", "\[IndentingNewLine]", RowBox[{"b", "=", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"qmp", "^", "2"}], ")"}], "deltaf2"}], "-", RowBox[{ RowBox[{"(", RowBox[{"rmp", "^", "2"}], ")"}], "deltaf1"}]}], ")"}], "/", "denominator"}]}], ";", "\[IndentingNewLine]", RowBox[{"Return", "[", RowBox[{"{", RowBox[{"a", ",", "b", ",", "c"}], "}"}], "]"}], ";"}]}], " ", "\[IndentingNewLine]", "]"}]}]], "Input"] }, Open ]], Cell[CellGroupData[{ Cell["Do an example and Plot to Show it works", "Section"], Cell["", "Text"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ParabolicFit", "[", RowBox[{"1", ",", "5", ",", "7", ",", " ", "2", ",", "3", ",", RowBox[{"-", "2"}]}], "]"}]], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"-", FractionBox["11", "24"]}], ",", FractionBox["25", "12"], ",", "2"}], "}"}]], "Output", CellChangeTimes->{3.39860154129207*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"dots", "=", " ", RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "2"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"5", ",", "3"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"7", ",", RowBox[{"-", "2"}]}], "}"}]}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", " ", RowBox[{"{", RowBox[{ RowBox[{"PointSize", "[", ".03", "]"}], ",", RowBox[{"Hue", "[", "0", "]"}]}], "}"}]}], ",", " ", RowBox[{"DisplayFunction", "\[Rule]", " ", "Identity"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"curve", "=", RowBox[{"Plot", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "11"}], "/", "24"}], ")"}], RowBox[{ RowBox[{"(", RowBox[{"x", "-", "1"}], ")"}], "^", "2"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"25", "/", "12"}], ")"}], RowBox[{"(", RowBox[{"x", "-", "1"}], ")"}]}], "+", "2"}], ",", " ", RowBox[{"{", RowBox[{"x", ",", "0", ",", "7"}], "}"}], ",", " ", RowBox[{"DisplayFunction", "\[Rule]", " ", "Identity"}]}], "]"}]}], ";", RowBox[{"Show", "[", RowBox[{"curve", ",", "dots", ",", " ", RowBox[{"DisplayFunction", "\[Rule]", " ", "$DisplayFunction"}]}], "]"}]}]}], "Input", CellChangeTimes->{3.39860155566932*^9}], Cell[BoxData[ GraphicsBox[{{{}, {}, {Hue[0.67, 0.6, 0.6], LineBox[CompressedData[" 1:eJwV13k8VOsbAPCxDDNnZjjHlpC9sieuZMk8KUuyRoWuci1l35JUtiwRrrKE UpIUZU1ayYkiiXCzptyxFD9taLFcmt/rr/l8P+e873ned97zPM9R8Aje481L oVAYPBTK6u815/IHKRoJJswvJsFOe0fJkiWJ+Ou5f7GVD0kNXZAfJW3SqKrR uaHssg1sOc4nDln8JmrZN/c0e/meF5yL55Adx7ceds/NZIeUsyuq3f8lR3dJ FBzIvcaW160oVJZ+R37WSdJzyb3DrvFOmnW+0keGljbwOuc2sUOqQgr7K5tI La3Ki7WO/7DLgmo/3ZzqYauTsyExhsPsiAd9RR1vJ9g8yd3+u5dG2cVtqtnR C9/Yt3L26FPvf2Tflv6RTdb9ZOdG2xWPyn9ii8QPt/k0LbHJjogb/tPf2O8z dynE/4/L9n+iutBpO8cOue4WOeHIC5LxewPnY3+w82+1azz05ofxE7lVxzm/ 2JdlJZ/5RwmA/b857LOSi2zzhOafOlY00P0R+CPQ7D92uikkDHyig57xn8du XVthl/TVN9WfZcDQrKB5RB+XPXJWolF8kgm2gV6z1jI88JO3zTvfQAj8DGod LC15oVQ2XGGHmzDsbxhxYbrzQUCZYE+UMg5t+Xlr2r34gH3lZtTeDTgYHrPv S/blg+8uF0Y0VHCQ1Wq25zvKBxR1j5YBdRw+Xr1hsZTEB7fVH8nK6uIQkRCg N1XOB1D6JyVwOw65Vkv483k+iEg+ESvqhkPvoHjbqXP8wP8ZC0w4j4OCdlDZ VDY/9F4cvP8tE4eg5NaUvfn8YHEzcNQ1Gwf6luO7Nl3jh3XDXwY1c3EwyR5s H6vlh5g6c4fWAhzKbAo6d/Wj6xtcjzTdxOFUs/ybNeuo8IKfpsqpR+tZe+Ju ogIVVJOb0lSe4CAR2pM9u54Km3xMu4IbcaiRi3d6pYWum19aWHqKw0TUeF8M UKE3tLiXrxUHW/2bgx88qcCfIlb7bxcOihXqI3dvU+Gf6ljep2M49GvKBMRW o/GVDWPL4zicrWYuWdVRQTvoSejWDzjM1H4RH3tCBVsj7HL5JA6Nj6pt8B4q tPP0KMR/xsHlhe6TgHkqPGOLzn37icP5MaPL63cKwONYvxuvBQnY6aWhNmsp ACLf5Q6O0giY/yDzsMFGAKTULrbM0Qk4+L+VN477BcBgvL+IYBKgMUMy4v0E 4OHp2BQ2TkDbyo5TI+eRg5Nk/pQkgCtp7Zr3XgAahL32528kQMQ+McN6TAA2 J2RVR6kQoJzc0MwzKQCm76ubD6kSYPlTQz1gRgCMzB6sVVQnIKuHtQx8gpA0 sZCap0XA+tSuK9MqgnDc3eKkmR4Bu//bM7otXBCGbIj0yu0oPp008e+RgjCe I2Pqa0pAiO+zXWXRgpDaLHFXaQcBOQO6tSLJglDRM9CSs5OAkbtiCZMXBcGK TK86YkFAaEC/ciaJxg80pXdbE5D33sVnAqNB9m/nqKp9BCSe20xaCtNg+qpc g/5+AsK20yUqRWlgW7LCIZFtbzx8dkyGBkb93KEOZwL4A9fICWrRoGRE5uyQ K7p/ubdP1YEGz08/Gn5ykAD3qgqNjL00GHH8FaB3iAAb98SEORcaCJy1Gi9H Vnmuu7negwaRMgaxF9wJ+DctK906nAazk1rhBzwIsJa23xGcRwM1InpPrTcB hp0bL70poAFNTFtP8jABG2O5M/pFNJg7sbQchcw7XlVIuUWDT0F9RjuOEPDo Nmsp8zENQlauRzT7ELDB8FVN3XsUX8V1Tpo/AWKfiwXXjtHAi9MhxkHmKTx5 MPojDUSOs7R0Awh4x6PGNP9Gg6zOXuF+5OyXKUcGKHRQlos6LxREAMXFXHZJ iQ5ZTppnrEIIcNK7r/hEhQ4e3v/Lz0IuIzZsjNOkQ/uw+KEhZPt2qraAPh1c F0mqZygBV41atotY0SF13rnhcBgBc2v0zPts6SCkl+9Xgmz+o8Qq35EOaoMW HRzkrxWJjrJudFC0W7i19ygBJrI7vdVC6KBEKXy9KRydp6W7vl/C0fx+7pKe yB/7lYJqTqDxhtniF5AzzvEd35JAhzTMSu0n8pj/0VOLyXRYx3i9XekYAVss x2Mb0ulgHtrMtEd+T3mWYpqL5gtabriBrP1e529qAR0sqabFr5ETHxVntl1F 6+NZMPyFrBEWf8n2Fh2y1e+d2R5BQNT49grXBjr8YRqe3IbcTd6pWddEh17q NtoksvJlhXucFjpgmyKM+Y8T0OHE03i4iw7GtQakIbL85tBm1V46bGE+Zzkh h7NGWz8P0uHFq0bhAGTp1qddYWN0qDhF1cpHDirW7tWbpANvgINVBXJzTNHg Asrrn5tYUiSyxAH8ff0MHYKiDa90I/vpx43G/KRDUuHAGw5yo+jMh+1LKN71 H558QxaZOTTNz6VDJOF5YAX5cEfX1xd8GDh5OlTQIwl4XMb+nkrDoFWqpkYM WSipet6GhcFr/VgfWWSPv+SWcREMAh43vtmAfH/bOUqvBAbFxUFzmsiYFJc/ TxqDoJmM57rIB38F0V3lMei+Lme1Fbn2nxHWuvUYvGyQTzBCFqi2FeGoYuBi lBW8DdklrVHiuhYGMXLhwibIFUe0pA/rYmDs2+qxep2ys1BOdSsGGZJnfVfH O8kLKX82xiBSs0Fpdf6y5WiV6u0YLJQeylx9/vLgF40wcwyEUo7WrsZnf89t s95uDB53zp9Zjb8ks1NvwQ6D35GfhVbXtxC4zbDeCYPBsw5Wq+vf2GL9l6YL BrXOXtzV/dkn82fKVTcMLAXPZ/1G+5d01L+a8MDg8oVevjnkuvaT/QmHMajm 03GcQB5XSF356YfW41AW17e6/ycuKvsEY8CM1U9vQQ7d+DDMOhKDgyeKthUj F8W8uNgYhYGy+cnxDOTuvv6n2qcxKJsLOHISWSvxp5B4KgZT05fptshuw/xb kjMwKNnyduMW5HQdMbfFLAxkPbSV1yFPc3TK3xdgYOW88cHq+Vy71fQfuyIM lhXbbF4hW55zWGwqwSCw83RzJfLNbSEWNysxONPENg9G7suJDZKsxUCLsXOv LTL/54wLqfcxuPDHX6YayJ6XKseDSAwubfh4fwK9L1mzDdjoMwwMv9mYkMjN lh2bHdswyMt7fT0fWX5+Ola/B4OzyVJ0S+R3jioylHEM7CfeGV9E7y+jXH9H 2CQG724z7vkhG/Ja+E18wmCDgyPDCDm/xvth2w8MZo019wyifLGXdd0pU4AB fpmMXVTkRK/aU3wMBkzs5antRPmlrr6p+JgwA6SWhBdzkEX8ODMukgxoHD2u II/8unXd34pqDJjWDBtWQflrZZ1mXY4WA+xntPSmUH7TOGY8LKDLgKNZ/IE3 kdOUDqh+MmLA7LGlCDlky7j8lrs2DOh+hP/gCyagyUD0985Q9Dy9yaxclF/t zf1eFR1jAK4rFmqGPLKnKX/5BANyTYq2fkf5eck/WLcungGuMX0Ju5F1Cl/5 KOUwoK6oDvvmS0Axb+IbngfIIWZveVD+1xYeLnKrZ8AzC0/bq6g+NErrBD0i GeClwSk3Qn77B4cW2obmk51XDEX1ReSIsQlniAGSNy7ydHkSkPDqR1njMgM4 2+7E7UP1SWRwd4QUDxNoX8nUMVS/iiaKd0RQmdBcrhkTiNywYj+iKcSEabsN CnGo/v3YVCl6RY4JLKP0/fl/EuCV4x1zypQJhZfVzuSg+mnm1u+4NYUJX9Lv TNfbof8vp+X1XDoTxqWcLSWRtV7V7arMZELYPtGDR21R/jLIZisWMMGr/8Wd 9TYETInZqzOrmJC4qLJ00oqA1I6XvJw3TGjveCU0hep/l9GTO8lyLHAtzbbo 2ErA87AKDVNlFtyc6irl6hPw8FZB6bIKC9IH5J9uRi5ac/JKqA4L1Bopa7NQ /xH2fcvZA2bo/qqH/jt0UL0sr3Hf5M+CIhn6BR/Uv7hIlQgP3GdBgYbGfIMM qucdOpVQz4J6orRjQBrlh+gmq9skCyz5zyfOSqHzODqSFN3GgiqiokBhLQH9 ZWuXld6yoLiTeSRcnIA9+uemQlZY8Ed7SNgEC9V/p5NN2E4hmOtrPvB1BQcZ AfqhcEshGPDhxnxYxuHTg7zl99ZCoEq/nTj8Hw6p0vf07+wVAsmwXfrPF1G/ Ovatcv8RIZD/UGOQgvpBs1DvSyWpaLyd4t23qF80ybA/atIjBGP65yJvvcUh 3O94o0C/EAS43yj3HMLhtnkhveutENw6sOuczCAO4txPVw9NCIFhu6FPah8O n4OTO+LmheB+8jP3fd04FDg0rn+2ThikEgMCqltwWBDTGLTwE4bCeJmv3tU4 DHUmFJbx4BCT++PSRAwOJRJGeTz8OCTWsIL3RKN+/dDcOVcBFOfhAW7jKRx4 Z91PMxk4UG5E7c+JxEFDzMQrRAwH62ChfzeF4RDruqC6dSMO5vp2dmreOKyf DLj3whr14TlmdRpWOIRS9nV8zMPhWBL1/BEch7CMNbZ0ZXSugk88NtMSBkPf kXieBnQuI0raErgs8L+awxjZJQIpFmklOSQTTpze7vatVwQSn5a6L3owwK7m +1/9/qKQoWTQlIzq+N+GFx5fWxCF3zJiTy3XoD6SMG41zRQDjd1bbWQVBIB3 5z/r2+TEgZZsVvr3Bn6wVTxkG18uDu5QZHVFhBc8bGZVRXZIQNyAvaauE/rO ehPmcqZbAmqyP17VSVtiK5YuGnX7roGnlfd6urGf7B0bjk/aUiSBUozHxpV8 ZXdSJj7eSZeEmX0HN457j7PlVPZdnFZeC3E3NleFObWzh24/uCp4dy1QEtrZ +0dekKFDwsEuJlLwtEo8IKlnjOx2CpyN65MC2stwX5GaGTJf/fGLdA9pSDFu Xj7/fJ7caMd0CPpPGuIMVD/kNXHJwefX8uKQJbXfuWU2c8lUQ/2RLGSOQ9NA 6jMu+XW9h/8D5IcSil3RLVzy/vKDJMqyNLTd2Nbh+ZJLmt/2fJyFTBl6x93U wyUPU+uVHqxIQ352vfVLDpeUPGXv9xK56MvdT82jXLJ95kPNMPJg2EhawxiX 1HqHm1B+o3h+v+yqnuCSv2oPO+9Cpkj/5507xSWT3EX+HkYG35ZXnjNcckt/ 6ZsvyFPhnGi3WS45tXubFIUrDd3lmpv3z3FJa32fUmXkwUsHC6x+cEkxocam QGSVGKmYzQtcsjXBkRaH/HSs1UB9kUtGLk7ZZiFHht34pbzEJd9NiL67jwyi i0cll7lkhustxZfI2spBf4iscEnoNvEdRnZ2F//F+M0l58x6q78gF734+pDK 5ZIl9b6/uMgzLktRXOT/A44qIaQ= "]]}}, {Hue[0], PointSize[0.03], PointBox[{{1., 2.}, {5., 3.}, {7., -2.}}]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesOrigin->{0, 0}, DisplayFunction->Identity, PlotRange->{{0, 7}, {-1.9999995119047718`, 4.367424241507136}}, PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Scaled[0.02]}]], "Output", CellChangeTimes->{3.398601557625267*^9}] }, Open ]], Cell["", "Text"] }, Open ]], Cell[CellGroupData[{ Cell["Muller's Method Implementation and Examples", "Section"], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"Horner3", "[", RowBox[{ RowBox[{"A_", "?", "ListQ"}], ",", "x0_"}], "]"}], ":=", " ", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{"y", ",", "a"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"a", "=", "A"}], ";", RowBox[{"y", "=", RowBox[{"Last", "[", "a", "]"}]}], ";", RowBox[{"a", "=", RowBox[{"Most", "[", "a", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"While", "[", RowBox[{ RowBox[{ RowBox[{"Length", "[", "a", "]"}], ">", "1"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "=", RowBox[{ RowBox[{"x0", "*", "y"}], "+", RowBox[{"Last", "[", "a", "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"a", "=", RowBox[{"Most", "[", "a", "]"}]}], ";"}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", RowBox[{"y", "=", RowBox[{ RowBox[{"x0", "*", "y"}], "+", RowBox[{"Last", "[", "a", "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"Return", "[", "y", "]"}]}]}], "\[IndentingNewLine]", "]"}]}], ";"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"Muller", "[", RowBox[{ RowBox[{"A_", "?", "ListQ"}], ",", " ", "x0_", ",", "x1_", ",", " ", "x2_", ",", " ", RowBox[{"\[Epsilon]_:", "0.001"}], ",", RowBox[{"Nmax_:", "10"}]}], "]"}], ":=", " ", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{ "p", ",", " ", "q", ",", " ", "r", ",", "i", ",", " ", "fp", ",", "fq", ",", "fr", ",", "a", ",", "b", ",", "c", ",", "x", ",", "discr"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"p", "=", "x2"}], ";", " ", RowBox[{"q", "=", "x1"}], ";", " ", RowBox[{"r", "=", "x0"}], ";", "\[IndentingNewLine]", RowBox[{"i", "=", "1"}], ";", "\[IndentingNewLine]", RowBox[{"\[CapitalDelta]", "=", "\[Infinity]"}], ";", "\[IndentingNewLine]", RowBox[{"fp", "=", RowBox[{"Horner3", "[", RowBox[{"A", ",", "p"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"fq", "=", RowBox[{"Horner3", "[", RowBox[{"A", ",", "q"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"fr", "=", RowBox[{"Horner3", "[", RowBox[{"A", ",", "r"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"Print", "[", RowBox[{"i", ",", "\"\< \>\"", ",", "p", ",", "\"\< \>\"", ",", "fp"}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"While", "[", RowBox[{ RowBox[{ RowBox[{"i", "<", "Nmax"}], " ", "\[And]", " ", RowBox[{ RowBox[{"Abs", "[", "\[CapitalDelta]", "]"}], " ", ">", "\[Epsilon]"}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{"a", ",", "b", ",", "c"}], "}"}], "=", RowBox[{"ParabolicFit", "[", RowBox[{"p", ",", "q", ",", "r", ",", "fp", ",", "fq", ",", "fr"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"discr", "=", SqrtBox[ RowBox[{ SuperscriptBox["b", "2"], "-", RowBox[{"4", "*", "a", "*", "c"}]}]]}], ";", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Abs", "[", RowBox[{"b", "-", "discr"}], "]"}], "<", RowBox[{"Abs", "[", RowBox[{"b", "+", "discr"}], "]"}]}], ",", "\[IndentingNewLine]", RowBox[{"\[CapitalDelta]", "=", RowBox[{ RowBox[{"-", "2"}], RowBox[{"c", "/", RowBox[{"(", RowBox[{"b", "+", "discr"}], ")"}]}]}]}], ",", " ", "\[IndentingNewLine]", RowBox[{"\[CapitalDelta]", "=", RowBox[{ RowBox[{"-", "2"}], RowBox[{"c", "/", RowBox[{"(", RowBox[{"b", "-", "discr"}], ")"}]}]}]}]}], "]"}], ";", " ", "\[IndentingNewLine]", RowBox[{"x", "=", RowBox[{"p", "+", "\[CapitalDelta]"}]}], ";", "\[IndentingNewLine]", RowBox[{"r", "=", "q"}], ";", RowBox[{"fr", "=", "fq"}], ";", "\[IndentingNewLine]", RowBox[{"q", "=", "p"}], ";", RowBox[{"fq", "=", "fp"}], ";", "\[IndentingNewLine]", RowBox[{"p", "=", "x"}], ";", RowBox[{"fp", "=", RowBox[{"Horner3", "[", RowBox[{"A", ",", "p"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"i", "=", RowBox[{"i", "+", "1"}]}], ";", "\[IndentingNewLine]", RowBox[{"Print", "[", RowBox[{"i", ",", "\"\< \>\"", ",", "x"}], "]"}], ";"}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", RowBox[{"Return", "[", "p", "]"}], ";"}]}], "\[IndentingNewLine]", "]"}]}]], "Input"], Cell["", "Text"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Muller", "[", RowBox[{ RowBox[{"{", RowBox[{"1", ",", RowBox[{"-", "10"}], ",", "1"}], "}"}], ",", " ", RowBox[{"-", "5.`25"}], ",", " ", "5", ",", " ", "7", ",", RowBox[{"10", "^", RowBox[{"-", "25"}]}]}], "]"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ InterpretationBox[ RowBox[{ "1", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "7", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", RowBox[{"-", "20"}]}], SequenceForm[1, " ", 7, " ", -20], Editable->False]], "Print", CellChangeTimes->{3.398601571965*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "2", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "9.89897948556635619639456814941178278405`25.422454610311895"}], SequenceForm[2, " ", 9.8989794855663561963945681`25.422454610311895], Editable->False]], "Print", CellChangeTimes->{3.398601571968337*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "3", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "9.89897948556635619639456814940658940414`25.11696980693702"}], SequenceForm[3, " ", 9.8989794855663561963945681`25.11696980693702], Editable->False]], "Print", CellChangeTimes->{3.398601571972378*^9}] }, Open ]], Cell[BoxData["9.89897948556635619639456814940658940414`25.11696980693702"], \ "Output", CellChangeTimes->{3.398601571975727*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"Muller", "[", RowBox[{ RowBox[{"{", RowBox[{"1", ",", RowBox[{"-", "10"}], ",", "1.`25"}], "}"}], ",", " ", RowBox[{"-", "15"}], ",", " ", RowBox[{"-", "5"}], ",", " ", RowBox[{"-", "1"}], ",", RowBox[{"10", "^", RowBox[{"-", "25"}]}]}], "]"}], "\[IndentingNewLine]"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ InterpretationBox[ RowBox[{"1", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", RowBox[{"-", "1"}], "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "12.`26.079181246047625"}], SequenceForm[1, " ", -1, " ", 12.`26.079181246047625], Editable->False]], "Print", CellChangeTimes->{3.398601574106414*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "2", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "0.10102051443364380360543185058821721604`23.675947728839844"}], SequenceForm[2, " ", 0.1010205144336438036054319`23.675947728839844], Editable->False]], "Print", CellChangeTimes->{3.398601574109618*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "3", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "0.10102051443364380360543185058078882138`23.37035785790146"}], SequenceForm[3, " ", 0.1010205144336438036054319`23.37035785790146], Editable->False]], "Print", CellChangeTimes->{3.3986015741135798`*^9}] }, Open ]], Cell[BoxData["0.10102051443364380360543185058078882138`23.37035785790146"], \ "Output", CellChangeTimes->{3.398601574117743*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Muller", "[", RowBox[{ RowBox[{"{", RowBox[{"6", ",", "20", ",", "5", ",", RowBox[{"-", "40"}], ",", "16"}], "}"}], ",", "2.5", ",", "2", ",", "2.25", ",", RowBox[{"10", "^", RowBox[{"-", "5"}]}], ",", "25"}], "]"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ InterpretationBox[ RowBox[{ "1", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "2.25`", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "30.75`"}], SequenceForm[1, " ", 2.25, " ", 30.75], Editable->False]], "Print", CellChangeTimes->{3.3986015770334253`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "2", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.960592288479808`"}], SequenceForm[2, " ", 1.960592288479808], Editable->False]], "Print", CellChangeTimes->{3.3986015770363894`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "3", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.9705636020344424`"}], SequenceForm[3, " ", 1.9705636020344424`], Editable->False]], "Print", CellChangeTimes->{3.398601577039824*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "4", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.9704465385226368`"}], SequenceForm[4, " ", 1.9704465385226368`], Editable->False]], "Print", CellChangeTimes->{3.3986015770430927`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "5", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.9704460787306015`"}], SequenceForm[5, " ", 1.9704460787306015`], Editable->False]], "Print", CellChangeTimes->{3.3986015770463037`*^9}] }, Open ]], Cell[BoxData["1.9704460787306015`"], "Output", CellChangeTimes->{3.3986015770494757`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Muller", "[", RowBox[{ RowBox[{"{", RowBox[{"6", ",", "20", ",", "5", ",", RowBox[{"-", "40"}], ",", "16"}], "}"}], ",", ".5", ",", " ", RowBox[{"-", ".5"}], ",", " ", "0", ",", RowBox[{"10", "^", RowBox[{"-", "5"}]}], ",", "25"}], "]"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ InterpretationBox[ RowBox[{ "1", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "0", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "6"}], SequenceForm[1, " ", 0, " ", 6], Editable->False]], "Print", CellChangeTimes->{3.3986015798588037`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"2", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", RowBox[{ RowBox[{"-", "0.5555555555555556`"}], "-", RowBox[{"0.5983516452371671`", " ", "\[ImaginaryI]"}]}]}], SequenceForm[2, " ", Complex[-0.5555555555555556, -0.5983516452371671]], Editable->False]], "Print", CellChangeTimes->{3.39860157986287*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"3", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", RowBox[{ RowBox[{"-", "0.4354502836022195`"}], "-", RowBox[{"0.10210124885056898`", " ", "\[ImaginaryI]"}]}]}], SequenceForm[3, " ", Complex[-0.4354502836022195, -0.10210124885056898`]], Editable->False]], "Print", CellChangeTimes->{3.39860157986624*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"4", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", RowBox[{ RowBox[{"-", "0.3906314607027883`"}], "-", RowBox[{"0.14185223238975053`", " ", "\[ImaginaryI]"}]}]}], SequenceForm[4, " ", Complex[-0.3906314607027883, -0.14185223238975053`]], Editable->False]], "Print", CellChangeTimes->{3.3986015798694477`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"5", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", RowBox[{ RowBox[{"-", "0.35769842880905456`"}], "-", RowBox[{"0.1699262688971378`", " ", "\[ImaginaryI]"}]}]}], SequenceForm[5, " ", Complex[-0.35769842880905456`, -0.1699262688971378]], Editable->False]], "Print", CellChangeTimes->{3.3986015798726473`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"6", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", RowBox[{ RowBox[{"-", "0.356050666609772`"}], "-", RowBox[{"0.16285601376462905`", " ", "\[ImaginaryI]"}]}]}], SequenceForm[6, " ", Complex[-0.356050666609772, -0.16285601376462905`]], Editable->False]], "Print", CellChangeTimes->{3.398601579875893*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"7", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", RowBox[{ RowBox[{"-", "0.3560617022464755`"}], "-", RowBox[{"0.16275830739471184`", " ", "\[ImaginaryI]"}]}]}], SequenceForm[7, " ", Complex[-0.3560617022464755, -0.16275830739471184`]], Editable->False]], "Print", CellChangeTimes->{3.398601579879154*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"8", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", RowBox[{ RowBox[{"-", "0.3560617617471714`"}], "-", RowBox[{"0.16275838285153657`", " ", "\[ImaginaryI]"}]}]}], SequenceForm[8, " ", Complex[-0.3560617617471714, -0.16275838285153657`]], Editable->False]], "Print", CellChangeTimes->{3.398601579882408*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"-", "0.3560617617471714`"}], "-", RowBox[{"0.16275838285153657`", " ", "\[ImaginaryI]"}]}]], "Output", CellChangeTimes->{3.3986015798856573`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Plot", "[", RowBox[{ RowBox[{"6", "+", RowBox[{"20", "x"}], "+", RowBox[{"5", RowBox[{"x", "^", "2"}]}], "-", RowBox[{"40", RowBox[{"x", "^", "3"}]}], "+", RowBox[{"16", RowBox[{"x", "^", "4"}]}]}], ",", " ", RowBox[{"{", RowBox[{"x", ",", RowBox[{"-", "1"}], ",", "2.5"}], "}"}]}], "]"}]], "Input"], Cell[BoxData[ GraphicsBox[{{}, {}, {Hue[0.67, 0.6, 0.6], LineBox[CompressedData[" 1:eJw1mnc4le//wI8RkmTvcbKdCvnYyX0biYZRElpGKcoohUpllFSysyKSEklG Qxn3g5JC9l7nHMfe2znW7/n+8fuH63W95/1e57guu529TlxiJhAI5/Ef//ut VhH1qYx9BqXYW28TWzWBzblZ15KHp5H5NT+d2jcmkP/q+y2N+ilk97ujTp7d COpI8NzOrZpCk3oVPDFzhvB8vd+87I8pdP23806eXkOYo3p4iC9rCl0wMeA1 LDCEcHHo72zgFApVLXX1PWsIr92Tjc/VmEKj5ktXFL9C+DPy1V65lEnELZem 2TlpAH0/x9nzX51A0gkOZ9Os9eCQV9nzvc4TqB0LrFHYrwdt9gxXHLKfQL6n zbA/PHpQPUOb5H94Au0prN/t3KgLpyO7Gb2yEyhb5ZNtj5UudL1KTHnXO44u a5oP0Gx1oK1Mbr+uxTgquemVvuWrBX/1tfKdPDSOVE/2r5LttaBG0obpNf1x 5EkZsO/V14J8PJafXpHGkbBxy7bdLFrw38bcAxb2cUQ79egbOUoTmnZp765H Y+jMyoKoaoEG1I6qcnHcP4Zyw17yfWH9D46Szx59pzSGdtU8e2U0oQ6T9y+r T0qPoSqC+q+RJnW43qzE7M89hmIuffrqn6YOMcGI9MiJURQlZ7FSdEAdHnlp 11f+dhT9WxE4F+q/H559N2UrKTaK3uwSSbZgUYM7V0INXHhH0dO354tvdqnC 8sNEhWyOUWT7nbnEP1cV7h47saSxMoLizo8NbTupCof3FMceax1Bl+0b7qhn qECvguDGuxEjiOcibaHCbB8MLhU+0s08gh6+9UiS/0qCl4sD253pw2g+SmyV 7RkJHvs85jwxM4x0O3liRB1JUDC35O567zD68VN8ZzMnCb5/eeGj1LdhJPro p+ZzZ2VYfyeLx+XqMHJUuiMoKK4ERXS1OyZahtCN75wXCz/Jww2NdJebf4eQ 4SXJkcdP5SFVbfvsOjaEnhpWOjy7JA8/KPVwcOcNoXMVSx/EJeShgcgDPbWw IbRZVkS88EwOuqxUp97UH0KVRadLHW7Iwrwvthc3Mmno6jVv7etOu2FvSeAH xRQaEvqM1U4Z74aclTnz1rE0tOV59GWkwm7o+m8z8F0QDeVJ8yUbTRKh9Mi7 VMvzNGTgf7XL5TYRRomstL8WpqEnnqe4WDSloc/dBHPTp4OI6s11Q3tGHGYE VUR5Bw0iZZOiN2pl4rDx8URHsv8gmjOe1fd6Jg73xgHXGddB5EN7x1aqJA5p uSPBicaDKEbVOiXSVQza9muXjq9T0ZM/WiLr0yJQx7BTJdKLilSufHNp3y8E k+ebxNYvURFb6y7SV1YhuP6mls3tLBUJnFl+T+kQhIgN9RuZUxFnhsgO8QeC 0LTubcSyDBVVVu7a5G0SgCdP+0yfa6cgmYNihUPB/PAzh0f33zoK4vrTwXBz 4IdCP1yrtaso6A6ry6SWOj/slnB4xZtPQZzmr9lzqHzQkQotfj2hIN30t00z pnzQ04M7b+9BCpL+4uvTIMELG6U4kpP+o6CcpaFXv1d5oHojUygbiYIYXQVe 29p44JL60jmyEAXJyasbWUfwwIDVnp1xs2QkcMjMR46VBz4JyfZYe0NGs+rx g1L0nTAzyWTfX04yUr42EX7GgBOakfPqXJjIKDmsQW2LnRN2HlfvZMwMII/r ehWSzdvh5Ss9p5z+DiCo7Wf+9/J2GJq6z1rlwQCacrB060jkgFXszYdrRvvR M8NUdIuXHXK+eBe61taPMl2XrNOobPCEzN1fKlX9SOHfXWeDz2yQqi9n/CK1 H2F1qW977Ngg4YavgfOJfrTdkGDxJGsbNOgV01wv6UPtjUnQ1YEVhrrN+Khm 9yH3jiuPkTorrF+uKnSO70ONa5OLYTtY4VleD7U/3jgP2M65l7PAu6ZoT7x8 H+osD9rvrMwCi/NdZNQie5GbarVSKz8zJBjoOLkE9KIS0212HAtM0KyWKz3e rRcVll7wyGpmgh1DXyQ3jHtR9oDjhlwsE1wW4xD9u9qDTqReTREQZYL/Pcrd ddG5B/3tuuZ9T4cAuRVZfnZb9aCF+ZphujQBjtbY+1uDHnSMLjPyj50AX3Gx UwwkcI4yWWUu2AIccU6Fou3dSGfjjSfL+iagaha7Rv3sRhZfrO74Nm6Csg5u cfaibmSQWeUb/3YT3BAvDVmM7EZJ9gvqPdaboDdD6FSDOa4/32DrVrgBvpp4 bDfV6UZTVbyhmuEbIHq4qqxUoRuZZ4Y/y3HdAKak6wo5LN2om+k293PJDVBQ ULvysKwL2eU35MvHroPwkzK5ax+6UBDns+OkG+vAdcnf8UZyFyqoTs5rsV4H EroKfy74daEA/tPf+PjXQSgWmKy3vwuFqWiU/5e6BhydOywKpLuQ55no01ZB a+AAqwqLEncX4vrOlrvTdQ3MHu65KjjRiW7O73t4Zf8aONOgeXA2sxP9Pvhx PKqBAbSuh8+5xnaiTmeRU2+LGYCHf/BtX1AnKmXmuGObwQDVtlHctec70Zve 0qVwPwZQ6x8feCvSiX6Gba/TV2YAzkDDOAn2TpTzTn23sjADDO1ONItd6kCd rEs6edsYIPnSoYLA5g4UNkb5Gkijg5vsKZdWsA4UpVNArG+lA8vseVHPTx1o 1M7FJv8XHbBOpQWfCe9AlFgtGd1sOuiPWNZqvtOBrobyVvS9pIPvascnzNw6 kArD6D/hSDrwvMk4qWXagZw5k1NM/ejAXMia46NGBzoQEZF04BodyBVnlcrK dqBNJ79bVU500LVmI89L6EBGZbDjiQUdiC0ekd051Y4Clqa6eQ7RwdlJuHt7 dzvSlE/WM9Kng1c0Leltv9vRvsqx05IadEDu3SvJ9Lkd7ZIriHu9lw5k2mTE N9Lb0ZWHiyf+ydOBS72IKP15O/qwObKYIU0Hmb+4hZfutCPTqJleKTE6GC5j FZy73I4eiKR5GArSgdJXBt+UTTuKt/uxnZOXDtzzZnnGDNvR+wUNAf+ddJD7 bph7SKUdlXqqzz7lpIOpV71cFPF2VKvGxmTEQQcqCc2cfRztKD/4v5wUNjrw jqzh6FpqQ3qRHpbJ2+ig8HE5Wxu1DYm4JNjo4bz44DNrU0Mb6p+skgzAWcs/ h7m+tA0JbWjuuoDb+3mnE/5ktyFtU8fnVHY6KL4Sv/kzvg1Fnx1a34bHZziG r2Mhbejq14TKGi460LcPZpR6t6HxD7eNSTx0cN/af7X4XBtS9iujKQvQAWbu ufz5SBvaq/Frx28ROmAyuriYr92GOskyUkxSdGCk5zCfK9eGCOcSgzpk6eCh utXse942dN5dOMWcRAfVJNPpzM1WlJJb0WK7nw44ZPUn0yda0Yns0YhNHTo4 Iq4+ntLZimZ+79SAhnTwjF9pNPFXK+qwpiqLH6EDblYBWlRaKyrWsniTfI4O rNa3U8PDWxFXhPbSwSt0EL24NRB2uxWdN7+p4e9DBwJDEz2BJ1vR96Rk0U/P 6MC2j9IVAFuRm96pyJxEOkhs6+jw34ez7bYYvXd0IF5d1eLN3ooy65YHBKvw +Sn/3nRtsQV1GIatX2jG5+frp4YrlBZ04MonM30qPi9ZL2sdS1pQepGHQDYL AyiF+VRZebWgOxKTrMbHGcA90K3i2NkWhH5JKnE7McAH/wvIzLwFdY9I/3G9 xQAqbkdLoGwLqi+HquVpDKB5RLZIraMZzR6QMFxiMICvkWjB3p/NaJo5h02C fw0U6+36pFTQjPofHsis3bMG9Pes5RCfNaMhKV6RqvNrwIirJYMHNKOQ3tgK 5to1ICR6VMhqXzMC24nht4fWwIR81ZNI8WakmtWX5kJYB3GgyJt7tQmJCSZV R2qvg+EbsQZcBU2IL6CY78CHdfCs62Q3u2wT6i3fR/LN2gCOw7XHD/M2oeFq yDReswE0FowrQrcaUWPkJ+zP+Abo59J4v623EdXSXORrVTeBGhTwZXnRiNzq b/6uKd8Ere9aebfYGtEnu5UZ89ktYH05ifPeeAM69fLehAwXATICn8ytVTeg 3vY5jZNEAjz+2Q3RHzSg9p+vaObmBDgvRjqzNPcPbdaF+D1JJcADozkxk231 6J7ybI60JRP8F5zL1JNaixw88zcl25hh4Nz5r1x3a9GdlIfXfaeZoboj71UD u1pUbXfr0So7C3xx0Lf1NV8tIo24VfgfYIFnV0HWlcd/UW/rbapLJgsc92g5 tuz9B0n1y/3+do8VstkzEnhNfqOc36N3546xQaBqpmI2XommW49/Lz/BCcuP 6Po+/FmJLrirh/f4cUIDV1I59qoSSUW3fhVLxTmVy0L3ZCUSFhAibo3gvKPJ g1RagTKWWV/bP9gBD47Yf+SKwNBaZk4K8TMX9DXosMtwKUUCHn2GpzR2QQ8p qVOf5wtR7ozq1HFDPpgvwL7ea3IL3RsWfVD3QhjWtpoKNfwsAB++85/sTpSE 5rXRz42MKkHFLcfb5fm74Rfq0j9nQi1YMXL7j0VdDppUl5ZaKdSCaCf+qqP6 crA1OyTH4FgteHCU932EqRxc9OYNFUusBbIeotYsZ+Sg+tbeg82qdSB6QUIv K0QO5olf/GB4oR60untZN7fIwfc2zY+J5Q3A5a+AyH0veWjhcpo4c7AZzIku PQmJUYAhWPHrDMtmcFfzmkt0sgIslhSTtXVqBn+tPuckZihAmc5e+bKHzeCL Rvm12EIFuHTMec+z2maQRByU121WgClaHlqK9i2gXvLhZ8ijCMe3hxy7cLMV EIyzvcQfK0Lpy7R6vtBWcD3MuehYhCK0+XnIsjqhFdAvF5X5vlCE5fc5Tuwr wfWHnPi+vFGE0QvPT68xtYEdvB3q3UgR6vQlOidEtoGy3TuEKpcV4aP8PP+G nHag+aIlcOS8EjzxzTt/o6wdaERN1X+8pASJ5eqje5rawfHv4zPu15Rgae3X 02Er7cCvxSblz20luDhcrmVo0gFSa/PWheOUoLNkw2JhXwdQ02eRWqxWgvDp jHc8Txew/1Gu/0heGXJHF7z/KdsFOG0tOsT3KMPeBB/yvFYXYA2pTH+vpgz9 3q5YWp7rAh6d2h8zDyjDPGxLhSOnC3xyqFNUsFaGkiu7pm4bd4NCh/TYI3eU IcNFze2cbw+4cdR2B+OnMjSBvHrqT3pAeNeDdN0/yvC5xDwne0oPuOahpu1V rwyJbUW5+RU9YM6KYV/apgxNTbVnWXb2ggDug1z9Q8owRsnAL+dNL7CRrZw+ wEqCvazSZg++9oI+ZkFKFzsJylO2RGz+9ALeU13JHjtIsDip4sfGdC+Iv8yw 9eUjwX7OQ5tWen2g4WBZDBeRBElTRx+tNPYBfjpRw1yXBG/+2WtbN9gHWF54 7zTSJ8GytzsVXy/3gY3yvttqgAStzjfUHJHoBz8Nfm1NmZCgb8MJrleX+0Fe B8N9wpIEsQ//9d240w8OBMOwqBMkuD1MIO/w834w5ZZB3HOKBFNgu+VcYT84 /F7xG3QgwcoC+1iTzX7gb+yzZ8qFBLki9C6K8g4AuemPTQddSfCUu7jmtOwA +Djr0xtyhQRHZfraE8wHwIfxaoHZayTIHecoNhE3AAqf1/8yu0WCdt6GEyhr ANBCHxwy8SPB18dkSuN+DICj51MENG+ToMY22jkD8gDo/LsQvxxAgvtKihpr BMlgePry8GoICXJIu9yPliKDu7NH3wU/IkFaMN9eB0Uy4E95yM7ymASTj1x/ PKFDBl2Bx96NPMHfn0fU/GxIBkdH+chH8b93rfkaqQFHyMDBNqXoXTjur1vF gPsszknxRBhJgoMH+yfaL5LB+ftOX+9HkWD56+dJaR5koPD0TX9RNAnecptc VHtABvF7rw8wxZEgdTX7w8kMMii87HJ9byLej7P29hIfyGD6N91NOokEEzEO 9qEiMpjqgiMcySRoGXbZ2fcXGdDyrdzQS7y/k0I84B8ZHN7T6x6eQoJsVtVl 7B1kYFqZOm6ZSoKlwvKiSaNkcCLL//yPVySYcLe12mmODF6E+1k6pZHgjYGQ myQGGXSW+hVv4mxh/J/MAjMFGI4vJ8Skk6ByFrWhZAcFaMj6TUq8JkHWHTH3 HgpQgHnz/c+pOJM9Dfccl6SA0+WJ0wIZJFjSPNspqEABlzhtkoJxjtdKD+1X oYByTZOiYZyvJ1tqZGlTQI840wHDNyR4bHOT4gUp4O6igmoMzkrOeZE65hTw 6KRtRBfOLNXnDjKdoABlZY3Twpn4fCvvnPjjQAG/KtzCjuL843lpYowLBaz+ yiT64vxi7qrpmWsUEK+XIJWAs/cp8UXZWxSQG7gSmIfz0e9/X0/eo4AHpwJN S3FWkLxj+SWUAqKbZz0wnJmDlDfuRVJAogfX8g+c+2idOaaJeD6ryX0fcC42 C7Pb9ZoCHr++IhuLc1yuNltnNgWszR+q9sbZi2ekKL0Ql7ew/zLB+cjNeCe3 EgpYevRYahfO8p2Hdqn/pICz2RlNDfj7CPpLpYw6CmDSPdAXinNvWqZ7VRsF XAzXhho4f2OxEQnvpwD908H0Drx+MZdZqm1GKCAom5PjBs6etYU+krMUcG6m yp0ZZ7lY3oY8JipIbclVYsP7Q1jBAvw4qeCucFfQbbyfPQ7eJMhPBW1JqnsG 8X7HyDQ8apKjAueXD92T8PnwCL3/X/I+Kmh+nbE5hM+P2fg+irMWFWgvKQ8o 4bxVEK6/eJgKmh4PBkXj8/fO/dqhO9ZUkKHK5J2Jz6eF7DELwhlcnuU//wGf 31dxOxy5PKkAvjmUn5hAgqbHJ67E+FFBcEjM2L14EpzeVntdJIgKVKYVi06/ IEEDv6chcnFUYOlybmE0hgSHVd3Dc1KpQHRT7ko6vj8Ro+Yv1LKoYMH4odsx fL/67Ldn6f+gAm6dl9zBz0kwQD/srw2ZCgLWbtW64Psst3S5pXuMCkLWxCMQ vu91Hw/3Oi5QQb24aMeuhyQoLs0+fY1tEDydc6c/D8TniSmU99HeQSBcs+FB wO+J849LYju0BkGoT5U1Ab83nD6HZKPBIOiz2P1t8SYJ2tNYNVNPDIL1+GLW Ery0q9Uhdl/8BwFjYE9a2WUSTH/g4nQgeBB0W10IcbqE11PH2L3i2SCoezk1 RHfG9zGbOaD+1SD4FzEtx3GeBLXCg9KGfg4CKXpqSPJJ/F5bPxgW5KWBmCt2 6fH4vZbgvDDzUowG2JgfnQzG7/nPSoPV3XI0cDRJKMRZC2+txuZ2FW0a+EUq P01QI8FCoXv7TM/RwHOZ6u4CGRKc7blz0zebBmx/FPtPbsP79TSp0qGQBtqP /xtyZCbBXJ1iHlBCA+/dh47VbipDtxdLuez1NBBy6gVn8IoypFp4DyXO0oDg +GVLv1Fl2FZxybZEewh8bTjoFlSjDDW8H2amwSHgzROZewr/PIyRejMfYj4E Eo+MzUljytDiLjni2JkhkFrt8Cj+mzKs0Tjzu+8+Lmd+Td3xThn+yLLSIVQP gV1vdXeLBClDEVuvUNq/IVBzPK7vWIAy9GWNaK3pGALf/WPf+/kpw/+car2j x4ZAu6+p5ycPZZgrZpots3MYnLg415JlrwxfPT8gdshmGIT7L8yXqirDEB+F 9SfUYWCgXfszsFUJpoe3xhtODAMlNZGRhH/494u3wfvpC8Ng4rtn+tsa/PtF R7/rlW0jIPb8TF9mqRK8qJ/QZKo0At517/monqkETbdxZLF4joCMvK658OtK kCNhzPoefQS0WpxMObxNCcrlJ0xqMI+CvZpHP2tsKkL459DjSc5R8OYJf5LI iiL0X3tdekZiFPz8Mtb4c1QRjl1wUNQDoyA80HX/p1pF+Fepdn354SggfzwV VB+pCMN/5GZ584yBhvbCRMCrCG8qPUsdEh0D94vyC05zKMJz8W6xDrJjwNPk +JFLWwpQ5YbCAxOtMcDXIszuNKUAm5TTbUXOjoH6AdXtPTUKUCgpdhuWNQZK z0QI9N5TgOl+d1x2gXHgzPkuk0SWh0UaZtJ5HhMgU5ezpiFUDp5IvOQQ5DsB DLkMeZUC5OD8WvALmwcToFxidqf/dTmoVlW2Yy1qAsQ51qpt4d9nc63/Wz38 eQL4Sm6G5qjJwXeeUk0UxgRQipXv9O+UhYnvF4MFwybBmpTxTh4ZWXhP4vVw QMYUECr5ZSP8ajcsfMl+0JplFnB9q/90bVwcUgi5ifWqcyDXaUz71yg/pEio pvCZzQPHjLzJoR526N1cuo3FfwFgu5UTXhJGQaOBsnBp7CIIjNvekrRvCU0y h/zoz14Cjd3D0j7ynFiC40hVbscyCLz6onHtuwA2ei1Xir93GfDEFNR5tgpg erev375NXgZqCtumG6YFsP7oNRXTsWXgPRame1ZWEFOo2pU8wFgG+caxMsnP BLGv8jrefFIrwDv9YlycnRDWNh4m4X9xBWC/XUfqBoUxhZXjfv1XVsDs/jax 5jVhzJ+Fv9nEA5d7eO2t5BfBJCRSH/P6rgDoQb/vbCyCXTxeMJ/zeAUQtbXm 0l+LYIufumr6cnD7m0LXNuxFMf5bpFvGcysgPdJr/dZ3Mazuw+N7HMsrgCzL GL7SKIY9pNIe1TNWQNSRqSWDETFs0TIt/jTrKnBMHpyPFRDHWkkCxVeFV0Eg 16GQPg9xLJ68wYg9uAqIxp+zGsQkMJFjTYG0J6sA7tdtUj8niTUGqzzJjlgF s8P3Hpt4S2Jh359Fe8augnzR8i6NEEmMrmCasZqC279aLst/L4l1M5VUceWv AoI3xy3eeUns5bdMNs12nP1pC6fuS2FSsv7hj2ToAK7YuA3elcYKW29+oyvQ QbqTpvOTJ9LYodDrVI89dKCWQ9vNnSCNXRtz17HVoAPHZgH/ogJp7MenczQF UzqwWr3N7T4kjZ3SN9avcaMDO5+HqZqyRGx0ClzW96SD0ai/7itKRCwgTT8m /wYdKGnteZWlQsQymLVGE+/SweqdwyvNukRstkYpzu05ru+pW3jSkoiFn+Ke 5CygA56q41ozvkSMyL5D+MEXOni/h60v8i4RKypmN1r4TgfkmpkRyUAi1i3O lNhTSQezY3pjTE+ImBJ1wTi3lQ4Cn042jiQRsdLYWU9iF27/5mnjXCoRszw0 lRTXRwf5aTkFo6+JmO/74ZmAYTrQWcPaU7KJWJVn18vjq3RQo3BnQKKYiNkS 26sr1ukgKuSOVHgJERtvap7TZGIAf1uFkpFyIsajWXdYipMBrMIP1rr8ImLn 18oXpsUZIArE3fJoImLzuSVSLkQG4EjRazrYSsQenS8275BjAMJIxDZGOxHL rShIQ/sY4Mqo+xfNXiLGeJx5NArgXMV7rniIiE1Oj0pMGzMAZhQ60z1CxAZO 7Zs+asYAxQKnMsbGiNhP2a9R7NYM4Bh0V7l2ioh9fcpwunSKATrPbetPnSFi 7+cM/quyZwAyj4XP2Tki9hz9brvvzAAiGdaE1EUiFqjAldXrygCzHyU6pZaJ 2I3nVv56V/H3TFRbha8QsYuLceaJXgyQ/ve4w9AqXo8zXWLLPgxQczuHpswg YmaVkpMn/XH7waYB+zUidkDZuawggAF4Cj8fuLVOxFSi3kXsCmIAnSNg9O4G 3s+V8Qsej/B6NNq3eWwSMb7zqvtrn+D+5Ubnj2wRsW2/fJiVI/D6Df9R4yfs xlb3FLeExuD+bnTfr8Z5PGY9kxbPAJA8W3eRaTfWS4e+Ri8ZQImDzDmN8z/H R4fT0xjA7L6vkjPzbgz7/Udk8w2uH3dfsALnIhXu8TPv8XrZl1dwsuzG3r04 UfI9F6938ooEwDlxPT5cuIABwigjQudxfurSc+7WF7wevmZpl3EO+Cut2vId r+e5f2/P4Oy5/yJhfzneHweSmB7OTonvmyIqGWB1kmOSGWebrcmMyWoGaLTm 2/EVj2/quv/mkVoGUHs47noSZ536W4feN+DzpL9/phvPf4/GDyG2Vpxng2KP 4Sz1cnPEpRPnFh/j9/j7eZiNv1f04vHbkqan8Hoxuz1+Kk3B/b0LfiiJ82JD 7Zl7Q3j/PxUvauH1Hdbi2dczxgD54pUqenj9u1JtNnWmGSCwkSijiPenjjWp IX4et79xonwD72f51b70xWV8Pp0b58vo+L7quhrnb+HyS0FUBj4fcek5Atys ayBwOjrAb4mIPWafGbrKsQaIhdnX+haI2LU2vzBF3jXguHqH69Isvg/6pfaP BNdAeu9Mcug0EbN6Q9gzKIrbl1P1oyeJmOb1J/WvZNaAt/NTf6dRfH8761+t K6wB7HrnrMIwERMDfN4Oe9YAFO3QbRskYptcL/mENNaA1cwYY6Ifvx8+A4M+ Ori+7S1rG3x/BrtlvzTprwGemjTNt11ErOZ97unnh3D/56NdCPj+RZuUp7Ce xvMd1Lyf/Bvf57Zk8YozayBfbf+u9SoiZnTZL/meI+5v7303Qwyv5xO1xGW3 NdA4+1jW73/3oCEjdjRgDah9vWoRnkXE8hwf8L0NWgOEvyq1298QMfb5M9FO ofh7S9Jc3F8RsWIBwcjuSNze2DyqNQ5/j8Pjp3UZa4C8mXHzCH6/bky4bH+S 9b/8ht377xCx2gAYdigXz8/xR73Nrf/dT/qj8i84vxYMGnXD9592NSi/Bs9X x04z25qIZXpZ+cdN4/kqfPQNECdi64R9K1YLeD4nprr2ChCxUzHbfXeu4vZX hA+XcxExti+VPqFM64AY6FzLviaNuTE0vPwF1gH2vilQtUQaU3kk6npWbx0Q lGf8UlSlsVChZZoIWAfpYpJXt0tKYwNZzS5txuuALM72WpNTGov6+8zJ4vg6 CLRQkK4dlMLmd22ehY7rwJG1/cpQrBT2LZl6Ui4Uj0dY3PtpRBIbf3d+KfIp zgn66YHNkphUUXf8WgTOabf+Ey6TxEJrm7uaE/D4sbdlY6MlMdv1ygtB2esA dghzVWpLYkvn3lztr8Pjqf1e7PaXwJTcpXceacLlbcTvkxcksLO+L/O+tOHx RS/AClMJ7GdE7Fx4P65vs3OrgF8Ci0Mhfgdm8fzN2nfl5Yhj/+2+GJLIvwGw pPHl6ToxzHUfWY5VZAOQje+siuWLYcm6Z6u9JDYAbJixZY8Vw5hP2Gw/LI9z Zksk0U4Maw42iVzS2gDE7jufrQZEses0uZcnHDaA48a/qkcUEexT1lAh1+sN EDi+uqr2WQg70k9n6ni7AdJ/fC/ZeiaE0QS4rV/n4DzfsXHPRQgTCdaa0fyM 63+N0AzjFcICzz7ec+E3ns8hylC9myBmxaucWTC1AQg/b3+f4xTA5u9ce2Gr twmIG0+lNbh5sfCCQBoRbILA2XNYZA8PpjAa99+EMS7/J3Iz7j0P5nCqrOnB 8U2AlQcYJxryYFWq3NzZjpsAHoyMZR/gxuJon0LXQnG5i7/A7Vs7MC2LBb/0 FtyfbjpzqBgrpqpn6WHTicvt6kzo6SyYksIHZ46+TUCYmdAoVmDBxDadj3sP 4/J/6mHCaszY5sdmGUjH7e1E2q30CFj1zsL6AaktEFiZeCDbbhWdrveWlXbf AoRMgfGHUr3I6nudSIsnzp5MS/EJHejIWyXuxz44c5LOpXG2IP175NXpAJxN 0CWusipE3Gf1rzwS5/y+KxceVYDRcNXb57/g7Cr1hTV+BOgr+I1YEAgwsDWo X12IAN913fQJZCZAQrkb8cVnAuR5fmOrgJUAsfWzTB9OMkHagoew4HZc31hz Yz6WGT6vuHi4hxfX31J035DYBvvPnsi6LItz80iZmTUnDIrddznQFLfnCGCK zuCDk6Z7FgrMCBCS7+y33uCDtgylB4NHcH/naz972/FDkpNcgqklrh9CvPJv lwBsVhH/zWWH+5sqcd0eIgjl/nIoJrnh9sw629/eFYHVBNpwQTgutxTpyPsq BV3+fPzDF4lz3g2ziwrSkBDjl+sTjev/0/ONipeGB2R33NCMx+OpG64q6BNh /iH1jeI0AiT/MKMWEHfD5KeB/FghHl8mRlx4SQZ6CUjCf50EmC45FjlVqwBF lGyTxuWYIJZGjbfWUYWyeV47ZRWZIPmkFyfpripU0XgSdEaZCcKhR0i/XBWa wFK3un24foNAdpiMGvSyk9HL02KC6Z0p46U5arA6bLr7uhnOX64TP2H74a2x UAn6VSZIULkZxseiAZtzvqWxF+H+K5MsXtF0YPiOpWGNL7h8Qtd6hEsXHr6m ruL8jQnmK/U7W2rqwtJ9H0tLS5ig9z8nWkCoLnybn9F1/ScTdNwl8N2CpAf9 vz7n721jgo0KW7c6/Q5AqUqXx/mrTDCQfpS/h2QA3bt2edkbMENyXYj2R7oh tLMU3fkVMsP8osNCi/xG8NAvmQ98xswQ+0MUslQxgsQCzZHaw8ww7GJLkqGL EWwPO3PBwJoZFss3x4v+M4JGOlmWsheZYdTDssSULGP4wqZd92wgM/z//+f6 Pybbuys= "]]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesOrigin->{0, 0}, PlotRange->{{-1, 2.5}, {-8.506649949641684, 41.93829409095926}}, PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Scaled[0.02]}]], "Output", CellChangeTimes->{3.3986015824260893`*^9}] }, Open ]] }, Open ]] }, WindowSize->{879, 1145}, WindowMargins->{{190, Automatic}, {Automatic, 0}}, Magnification->1., FrontEndVersion->"6.0 for Mac OS X x86 (32-bit) (June 19, 2007)", StyleDefinitions->FrontEnd`FileName[{"Creative"}, "PastelColor.nb", CharacterEncoding -> "UTF-8"] ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[590, 23, 92, 4, 112, "Section"], Cell[685, 29, 49, 0, 34, "Subsection"] }, Open ]], Cell[CellGroupData[{ Cell[771, 34, 32, 0, 72, "Section"], Cell[806, 36, 1762, 48, 188, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[2605, 89, 58, 0, 72, "Section"], Cell[2666, 91, 16, 0, 26, "Text"], Cell[CellGroupData[{ Cell[2707, 95, 152, 3, 37, "Input"], Cell[2862, 100, 192, 6, 54, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[3091, 111, 1495, 46, 83, "Input"], Cell[4589, 159, 5908, 103, 244, "Output"] }, Open ]], Cell[10512, 265, 16, 0, 26, "Text"] }, Open ]], Cell[CellGroupData[{ Cell[10565, 270, 62, 0, 72, "Section"], Cell[10630, 272, 1230, 34, 158, "Input"], Cell[11863, 308, 3720, 93, 438, "Input"], Cell[15586, 403, 16, 0, 26, "Text"], Cell[CellGroupData[{ Cell[15627, 407, 277, 8, 37, "Input"], Cell[CellGroupData[{ Cell[15929, 419, 303, 8, 20, "Print"], Cell[16235, 429, 319, 7, 20, "Print"], Cell[16557, 438, 317, 7, 20, "Print"] }, Open ]], Cell[16889, 448, 129, 2, 37, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[17055, 455, 361, 11, 53, "Input"], Cell[CellGroupData[{ Cell[17441, 470, 343, 7, 20, "Print"], Cell[17787, 479, 319, 7, 20, "Print"], Cell[18109, 488, 319, 7, 20, "Print"] }, Open ]], Cell[18443, 498, 129, 2, 37, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[18609, 505, 280, 8, 37, "Input"], Cell[CellGroupData[{ Cell[18914, 517, 302, 7, 20, "Print"], Cell[19219, 526, 251, 7, 20, "Print"], Cell[19473, 535, 252, 7, 20, "Print"], Cell[19728, 544, 254, 7, 20, "Print"], Cell[19985, 553, 254, 7, 20, "Print"] }, Open ]], Cell[20254, 563, 90, 1, 37, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[20381, 569, 302, 8, 37, "Input"], Cell[CellGroupData[{ Cell[20708, 581, 286, 7, 20, "Print"], Cell[20997, 590, 377, 9, 20, "Print"], Cell[21377, 601, 380, 9, 20, "Print"], Cell[21760, 612, 383, 9, 20, "Print"], Cell[22146, 623, 383, 9, 20, "Print"], Cell[22532, 634, 379, 9, 20, "Print"], Cell[22914, 645, 381, 9, 20, "Print"], Cell[23298, 656, 381, 9, 20, "Print"] }, Open ]], Cell[23694, 668, 184, 4, 37, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[23915, 677, 382, 13, 37, "Input"], Cell[24300, 692, 12982, 219, 246, "Output"] }, Open ]] }, Open ]] } ] *) (* End of internal cache information *)