(* 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[     98904,       2306]
NotebookOptionsPosition[     93789,       2146]
NotebookOutlinePosition[     94197,       2162]
CellTagsIndexPosition[     94154,       2159]
WindowFrame->Normal
ContainsDynamic->False*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["\<\
Math 481A 
California State University, Northridge
Newton's Method
\[Copyright] 2007

version 6 Compatible\
\>", "Section",
 CellChangeTimes->{{3.397782792213933*^9, 3.397782800485338*^9}}],

Cell["\<\
\[Copyright]2007 BE Shapiro
\
\>", "Subsection",
 CellChangeTimes->{{3.397747927170925*^9, 3.397747928052795*^9}, {
  3.397748805779945*^9, 3.3977488403282413`*^9}}]
}, Closed]],

Cell[CellGroupData[{

Cell["Newton's Method Algorithm", "Section"],

Cell[BoxData[
 RowBox[{
  RowBox[{"NewtonsMethod", "[", 
   RowBox[{"f_", ",", " ", "p0_", ",", 
    RowBox[{"\[Epsilon]_:", "0.001"}], ",", 
    RowBox[{"Nmax_:", "10"}]}], "]"}], ":=", "\[IndentingNewLine]", " ", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"\[Delta]", ",", " ", "\[CapitalDelta]", ",", "i", ",", "p"}], 
     "}"}], ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{"\[CapitalDelta]", "[", "x_", "]"}], ":=", 
      RowBox[{
       RowBox[{"f", "[", "x", "]"}], "/", 
       RowBox[{
        RowBox[{"f", "'"}], "[", "x", "]"}]}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"p", "=", "p0"}], ";", " ", 
     RowBox[{"\[Delta]", "=", 
      RowBox[{"\[CapitalDelta]", "[", "p", "]"}]}], ";", " ", 
     RowBox[{"i", "=", "0"}], ";", "\[IndentingNewLine]", 
     RowBox[{"While", "[", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{"Abs", "[", "\[Delta]", "]"}], ">", "\[Epsilon]"}], " ", 
        "\[And]", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"i", "++"}], "<", "Nmax"}], ")"}]}], ",", 
       "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"p", "=", 
         RowBox[{"p", "-", "\[Delta]"}]}], ";", "\[IndentingNewLine]", 
        RowBox[{"\[Delta]", "=", 
         RowBox[{"\[CapitalDelta]", "[", "p", "]"}]}], ";", 
        "\[IndentingNewLine]", 
        RowBox[{"Print", "[", 
         RowBox[{
         "\"\<iteration:\>\"", ",", "i", ",", " ", "\"\< p=\>\"", ",", "p"}], 
         "]"}], ";"}]}], "\[IndentingNewLine]", "]"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"Return", "[", 
      RowBox[{"{", 
       RowBox[{"i", ",", "p"}], "}"}], "]"}], ";"}]}], " ", 
   "\[IndentingNewLine]", "]"}]}]], "Input"]
}, Closed]],

Cell[CellGroupData[{

Cell[TextData[{
 "Use Newton's Method to find ",
 Cell[BoxData[
  FormBox[
   SqrtBox["2"], TraditionalForm]]],
 " to 50 significant figures"
}], "Section"],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"f", "[", "x_", "]"}], ":=", 
   RowBox[{
    RowBox[{"x", "^", "2"}], "-", "2"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"NewtonsMethod", "[", 
  RowBox[{"f", ",", "1.5`53", ",", 
   SuperscriptBox["10", 
    RowBox[{"-", "50"}]]}], "]"}]}], "Input"],

Cell[CellGroupData[{

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "1", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   "1.4166666666666666666666666666666666666666666666666666666666666666667`52.\
66224719731129"}],
  SequenceForm[
  "iteration:", 1, " p=", 
   1.41666666666666666666666666666666666666666666666666667`52.66224719731129],
  
  Editable->False]], "Print",
 CellChangeTimes->{3.39778268709312*^9, 3.3982768778289967`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "2", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   "1.4142156862745098039215686274509803921568627450980392156862745098039`52.\
360089651515274"}],
  SequenceForm[
  "iteration:", 2, " p=", 
   1.41421568627450980392156862745098039215686274509803922`52.360089651515274]\
,
  Editable->False]], "Print",
 CellChangeTimes->{3.39778268709312*^9, 3.398276877831108*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "3", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   "1.4142135623746899106262955788901349101165596221157440445849050192001`52.\
0590586775016"}],
  SequenceForm[
  "iteration:", 3, " p=", 
   1.414213562374689910626295578890134910116559622115744`52.0590586775016],
  Editable->False]], "Print",
 CellChangeTimes->{3.39778268709312*^9, 3.398276877838221*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "4", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   "1.4142135623730950488016896235025302436149819257761974284982894986232`51.\
75802868183689"}],
  SequenceForm[
  "iteration:", 4, " p=", 
   1.4142135623730950488016896235025302436149819257761974`51.75802868183689],
  Editable->False]], "Print",
 CellChangeTimes->{3.39778268709312*^9, 3.398276877840528*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "5", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   "1.4142135623730950488016887242096980785696718753772340015610131331133`51.\
45699868617291"}],
  SequenceForm[
  "iteration:", 5, " p=", 
   1.414213562373095048801688724209698078569671875377234`51.45699868617291],
  Editable->False]], "Print",
 CellChangeTimes->{3.39778268709312*^9, 3.3982768778428173`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "6", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   "1.4142135623730950488016887242096980785696718753769480731766797379907`51.\
15596869050893"}],
  SequenceForm[
  "iteration:", 6, " p=", 
   1.4142135623730950488016887242096980785696718753769481`51.15596869050893],
  Editable->False]], "Print",
 CellChangeTimes->{3.39778268709312*^9, 3.398276877845152*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "6", ",", 
   "1.4142135623730950488016887242096980785696718753769480731766797379907`51.\
15596869050893"}], "}"}]], "Output",
 CellChangeTimes->{3.397782687374732*^9, 3.398276877847658*^9}]
}, Open  ]]
}, Closed]],

Cell[CellGroupData[{

Cell["Damped Newton's Method Algorithm", "Section",
 CellChangeTimes->{{3.397747720687557*^9, 3.397747721500311*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"DampedNewton", "[", 
   RowBox[{"f_", ",", " ", "p0_", ",", 
    RowBox[{"\[Epsilon]_:", "0.001"}], ",", 
    RowBox[{"Nmax_:", "10"}]}], "]"}], ":=", "\[IndentingNewLine]", " ", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
     "\[Delta]", ",", " ", "\[CapitalDelta]", ",", "i", ",", "p", ",", "j", 
      ",", "  ", "pnew", ",", " ", "fnew", ",", " ", "fold"}], "}"}], ",", 
    "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{"\[CapitalDelta]", "[", "x_", "]"}], ":=", 
      RowBox[{
       RowBox[{"f", "[", "x", "]"}], "/", 
       RowBox[{
        RowBox[{"f", "'"}], "[", "x", "]"}]}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"p", "=", 
      RowBox[{"1.0", "*", "p0"}]}], ";", " ", 
     RowBox[{"fnew", "=", 
      RowBox[{"f", "[", "p", "]"}]}], ";", " ", "\[IndentingNewLine]", 
     RowBox[{"\[Delta]", "=", 
      RowBox[{"\[CapitalDelta]", "[", "p", "]"}]}], ";", " ", 
     "\[IndentingNewLine]", 
     RowBox[{"i", "=", "0"}], ";", "\[IndentingNewLine]", 
     RowBox[{"While", "[", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{"Abs", "[", "\[Delta]", "]"}], ">", "\[Epsilon]"}], " ", 
        "\[And]", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"i", "++"}], "<", "Nmax"}], ")"}]}], ",", 
       "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"fold", "=", "fnew"}], ";", " ", "\[IndentingNewLine]", 
        RowBox[{"pnew", "=", 
         RowBox[{"p", "-", "\[Delta]"}]}], ";", "\[IndentingNewLine]", 
        RowBox[{"fnew", "=", 
         RowBox[{"f", "[", "pnew", "]"}]}], ";", " ", "\[IndentingNewLine]", 
        RowBox[{"j", "=", "0"}], ";", " ", "\[IndentingNewLine]", 
        RowBox[{"While", "[", 
         RowBox[{
          RowBox[{
           RowBox[{
            RowBox[{"Abs", "[", "fnew", "]"}], "\[GreaterEqual]", 
            RowBox[{"Abs", "[", " ", "fold", "]"}]}], " ", "\[And]", " ", 
           RowBox[{
            RowBox[{"j", "++"}], "<", "Nmax"}]}], ",", " ", 
          "\[IndentingNewLine]", 
          RowBox[{
           RowBox[{"\[Delta]", "=", 
            RowBox[{"\[Delta]", "/", "2"}]}], ";", " ", "\[IndentingNewLine]",
            
           RowBox[{"pnew", "=", 
            RowBox[{"p", "-", "\[Delta]"}]}], ";", " ", "\[IndentingNewLine]",
            
           RowBox[{"fnew", "=", 
            RowBox[{"f", "[", "pnew", "]"}]}], ";", "\[IndentingNewLine]", 
           RowBox[{"(*", " ", 
            RowBox[{"Print", "[", 
             RowBox[{
             "\"\<Revising \[Delta]=\>\"", ",", " ", "\[Delta]", ",", " ", 
              "\"\< \!\(\*SubscriptBox[\(p\), \(new\)]\)=\>\"", ",", " ", 
              "pnew"}], "]"}], "*)"}], ";"}]}], " ", "\[IndentingNewLine]", 
         "]"}], ";", "\[IndentingNewLine]", 
        RowBox[{"Print", "[", 
         RowBox[{
         "\"\<iteration:\>\"", ",", "i", ",", " ", "\"\< p=\>\"", ",", 
          "pnew"}], "]"}], ";", "\[IndentingNewLine]", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"Abs", "[", 
            RowBox[{"p", "-", "pnew"}], "]"}], "<", "\[Epsilon]"}], ",", " ", 
          
          RowBox[{"Break", "[", "]"}]}], "]"}], ";", " ", 
        "\[IndentingNewLine]", "\[IndentingNewLine]", 
        RowBox[{"p", "=", "pnew"}], ";", "  ", "\[IndentingNewLine]", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{
           RowBox[{
            RowBox[{"f", "'"}], "[", "p", "]"}], "\[Equal]", "0"}], ",", " ", 
          
          RowBox[{"Break", "[", "]"}]}], "]"}], ";", " ", 
        "\[IndentingNewLine]", 
        RowBox[{"\[Delta]", "=", 
         RowBox[{"\[CapitalDelta]", "[", "p", "]"}]}], ";"}]}], 
      "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", 
     RowBox[{"Return", "[", 
      RowBox[{"{", 
       RowBox[{"i", ",", "p"}], "}"}], "]"}], ";"}]}], " ", 
   "\[IndentingNewLine]", "]"}]}]], "Input",
 CellChangeTimes->{{3.397746800560541*^9, 3.397747103594524*^9}, {
   3.3977471748801823`*^9, 3.397747218537197*^9}, {3.397747265513595*^9, 
   3.397747309650807*^9}, {3.397747354940393*^9, 3.397747362265942*^9}, 
   3.397747400976905*^9, {3.397747466962101*^9, 3.3977474766165867`*^9}, {
   3.397747541035405*^9, 3.397747627845809*^9}, {3.3982798861714697`*^9, 
   3.3982798960046387`*^9}}]
}, Closed]],

Cell[CellGroupData[{

Cell["\<\
Use DampedNewton to Find Minimum of a function without a root\
\>", "Section",
 CellChangeTimes->{{3.397747760977425*^9, 3.397747764366465*^9}, {
  3.397747874539793*^9, 3.3977478791772623`*^9}, {3.397747912155539*^9, 
  3.3977479153224783`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"f", "[", "x_", "]"}], ":=", " ", 
  RowBox[{
   RowBox[{
    RowBox[{"(", 
     RowBox[{"x", "+", "3"}], ")"}], "^", "2"}], "+", "3"}]}]], "Input",
 CellChangeTimes->{{3.3977467353828993`*^9, 3.3977467405968237`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"f", "[", "x", "]"}], ",", " ", 
   RowBox[{"{", 
    RowBox[{"x", ",", 
     RowBox[{"-", "6"}], ",", " ", "0"}], "}"}], ",", " ", 
   RowBox[{"PlotRange", "\[Rule]", " ", 
    RowBox[{"{", 
     RowBox[{"0", ",", "10"}], "}"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.397747735597315*^9, 3.39774775499949*^9}}],

Cell[BoxData[
 GraphicsBox[{{}, {}, 
   {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJw11wk0VdsbAHBd7j373HtN1xyZypBIvVep9wwfQpShEpJ6iAyJRObQJGNS
lKQ8wkO9KA0a76dCg6EB6ZXMURpQLyH5n/9/rf9Z66yzfmudtdf59v6+b++j
4R2yzpclIiLyhbn/+7TWD6smOfIo8r9rHvCHyv/KTpVH/bC0Za19c8G3bH7R
gEAeZ6/dWxNYqQml2YIb0llyeN0o8OCFnRqg5rrw+aCCHH7z8+FKu6hDolU1
fTJbFs0lYmK++qnCz9hndpUysjh14ihX1VEFriTPywg5JYO33L0SzkTOhvW1
XR+4SjJ4e+g22zJXEfSeNNk7FQjwCulSu/FNHlR1XK5P6ghwTOvA3YRoOZgo
t7FILJVGsylYvX5KBsRFgxorDaSxoKnayfyyAEjyRNiNCikUU/IfGw6VBo0N
Gtp6elLom6x2T8pFCvpHxi8d8pHEd33ld1+3SoBzZoJzqKcE/rw7UrxEThwS
V0quLvWQwMIvbcF6kuJwbuKM+Ss3CWyxS3s6lxYHls+dhVZrJfBZ0DYvjWk+
VC2fJkoWEuif8KbCZoAP/N6Y2zhXAt+8aeuHq3yoWxKhLTUojiNJPV8H3Pmw
7J+g7xd2iGNBbNDByQs8qO0YbDgaII7T/tsu+pfzYM0L7xMRvsz7dLL8q7M8
8Gp1W2a6WRxztlmkN+byIL15ZXjTGnHsTz2m1buPB713VUaGF4hjasJfkh2u
PDh8rnFQ9z0fpxWP134S44FShc01/ls+8qOTBP0zXCguu5s00sPHqvQz+Z2T
XLheUqNV85KP15qtPDpHuNBfULzV5iEfH93emSLZyYXfsuO6fMv4+N1m+SnT
q1x4G6f/omgbH+P4h/MbtnPhkv2P5pPefIScDsFGPy7EqzbWH9nCx9V9Nqc+
eXNBDrdfjXfh44qW7F4tdy6sFD2X427Nx91vB+3qVnGhMGW+s0CbjzsEYfsa
tLngcVz76d4BHm7xnJPxcoAGXb9vD6J6eKhxVGTVxx4avhjVY0gnDxekEDn2
GxpSX/pWbW7joVJD0gC003BdpSRzRR0PN+wz9nxXT4P82bkOo8U8tP+Wfmu6
jIanleqPPX14aNpaK7VkJw0OCpOC0548rORLZovtoKEx/rn7Sw8etrM2FnUE
0FBvn/R+7QYertTt+yNtKw23hodpS2selro+2aLsSkOZzrVVWro8jF8lETRg
SkNCgX39u2EuRu5SvK8hTcMPjo6E9hAXC7ymS7dK0BATLOLi3c/FPRJdsuU8
GsJNqgf+ec1F98LcLFMODYGvFNmNTVw8ea5fInuKgItCv+WFSi7qretWfPeW
gEFmTO2uMC6eGPQN/B0JcMGsd3MIF8Ow2uL2bQJvR0RFbbdzcei7vpXlTQJ/
rstYqbaVi38mFz7YdJWAQP7PhsfruGg299jDO+cJfM+va5r3CxetrHKr1+QR
aLVP/SS5kIs1CvkLpXMJXPzpIDk1n4t1to+fv8whEOjZ4fRUg4vnGs/lR2YR
6Jw7/HyPNBfL7XfrDiQTuFsh9U/7CI2Bxnq6yVEECja1TdZ+oHE36bfKjCAQ
x89T/nuIRqUlDzPywgksDZ67eX83jQsr5PJv7yRQtnhZt+ETGlOkr11aGkAg
o2bT2+RKGtlNyg8+uhNwqy8d+y2YRoXxoGN7zJl4xXdv3R5I4/VAuzhNILB7
g2XrqW00Hu05wnlsSuBof9eVH1to7O5wm9E2JtDIUoq+40jj+MfWvVwjAmZm
aTPmi2mUurTnUYIBgeakjTt3GdA4a+Pwwy36BDyadXqK5tOoI3o30HwBgajN
9++KatIo8Ko4Ij2fwKXY6YN1Ahq9LP12dM4joF0Twrf9SnC4MCtw7hwCl2dM
9kSPENSLy0w0ViFgacP/VP6BYMHGb0YblQl4tpe1cAcIaqnOpJ5WIpD3tedo
YxvBC00/atfKE9A1rhKdfkowu9ai75AcgWv748MNmglW1oyW1Moy6yVQdjlc
T/BAVaokyBAQX+ys5HSNYFl9WZi/FIH8KM3UxGqCS1uX/XlTkoAejkxWVRL8
3m8QLGBs7ZjxWrqM4Gwp10+PxAkk7qgveJ5L8FV6ToIfj4DElWwpsRyCkpe3
X3vBJXD6h/feJVkEo99UHFnN+EbajHdOCsH3xj46pjSBr+XLtV1jCNqo3xtw
oAjcv7dJmBdB0CTojucbDoFjnfFub3YRnHk27RbK+BfpulTf7QT7RlbVFrIJ
zFowNLfcjyCnzHmxGeMnK3m3P2wl+CW6WLVbjEBI1NrPYR4EDYOlQvUZmx4N
T77mRjD+iH77C1Hme8+f0JhyJpjblF6TxPh8V+f6ffYE9XOV/v3EYvJvQuRj
nS3BUi1rxzLGq2XmJdHWBD80XdX2ZfzeOrDmiCnBjmAf+3ezmPg8M9a2/kYw
wa8wv4pxSkzVewUjgiMxmgGxjN2yn+/f9CtBg/LBS7aMdS98UykwZOZ7bCRa
mfF4g9LV3gUEKRfTOyMiBOp7jB21dQmWtLYkPGR8fOqPoYB5BCWCS+6UMPaV
27/3b3WC5prCmIOMlxiWzh5VIQijKpf9GYvZPqxeokSQfnUn0JHxc+8Pa6Lk
CJ7pKj6zgnFRnOTbm9IEf/5sdtRhvOv4Lwkz4gR1V5jEKzI2r9qgaMklqHr4
83xxxlKPoi4mcQh2TQ+sYzPu7jtl94hFMDBZfUKEcdX0nT7xGQprDQpU/usE
hd64tVMU9n70us1i7LCYLZ8zTmHLA/9XhLHqat3Kji8UHrpTHSFg/NFn9SqV
EQrpFotMNca34oN7/vhAoeuknNYixum5WTFnhyjcCQbLVzLedOmyzGA/hesL
kx9sYqzX+OK8Xg+FrDm67RGMJwcmrYI7KYyr5nnlMH44M6fr4ksK0efXgGuM
c5XMo/5to/DJooJPrxn7/+ojveIZhRUKjqNsZj2M7A9VxDVT6KRkGfoL4/bE
pteiDRT+CBofymFckjey2+YehaLCa92NjMMvy0imCSlsW3BrA8Xkh8zQRnNB
DYUvN2RdS2LcO2vPPxsuU8iX8S5+zPiS8p9hJ6uYeIcjJGSZ/HNyfFuiXk7h
tuEg90rG6gG0mU8JhfdlHIxZTP5+3qff8VchhZ9cwk65MT58dRfXMI/CIoPU
p5JM/j+eMxNskkZhpso7lTKmXiYTxXsMD1GYf7k+Sp2pL91+5fWa+yn8VfXT
53zGSRXLjahYCid9Zy8qIgTAaNfMk0AKvzcMjXcy9Rqclxh6bxuFDh9kF0Yy
9X16+nDfFW8KDXpTE2T5zPj3zjXkuVOYbOEcvonpB5edBjJ97Cgsb/Momc30
k97qL7NcrSlclfTiSitjaQVWuK0FhaEkv+eoNDN+p6rbwt8o3LyXlabM9KP5
gW7q3+cz8VhZaW1h+tmZ/Y+r0mgK+xWKIsuY/tj49qVmPJtCuZ7rxsWqBKZs
h7J3zmLWL1Z5TrEaU19S7OgNExzs1rplWKXBzP9pU3O1dxzcsSYsf0KLQPLV
i08vPeCgWV/NsbWGBK4qoWXJfQ6K1u0Qz1xEYCCu+coJ5KBjVFDRk8UELCyH
T8bVcHD/ehtVzyUEfrTM87Yu46D7612bq1YQCH13YuzlIQ7eMiqckLMisK6v
8rL0AQ7ufa8UHGvN9KvOhgjbRA5WnVOI6LchMPZkfLImioNrNhZBrR2TLzWu
orkBHIz0+ioscSIQcUhB1mU1B09rpfnN3szs53sN2zNsOPjzjrPs9S1MPsba
5NZZcrBh21nnTZ7MfhwcqbLUmIO5PS2ZFVsJRLu8mCdrwEHDA74jYYFM/9I6
sfSZBAcX7NnvfonZXz3UKsdpHgd1pk+bpcYQMFFquG5OcTAtzaTGN47AT/64
ycUZNh59zvPRSWTq96uL9ZERNh6TGihpO0Rg3z15V4dnbPQzTX3WeYKZT6/j
UY9y2OgYUNCSf4PAg+KTMhez2Bj/sChMkTlfUEP5F05ksDF8k/TqXCGTb8Fn
+30OsrFadlF28X0CB2KrnGbC2Rj3OizzczOBvTmPdJc4s9G4VlbPsZ9AzMOf
HacFbNzTdfPzEgENN/is8AMSbIyM0pplJEfDhCNbcjuXjcstJ61MFGmIaudZ
LWexkRc8EOCiSkNkn+LFZ6NizPn27uybesz5afqXFOqpGErfz1Z3WEnDjsV+
K0IzxXAiuOBWXgwNn6yFPwLTxHA0qvSCIJ6GYA8F9Dkkhg7pi/4+vJeGkEMN
1m4JYri0panlaDINoZ06zmYhYni9Tv/c4+M07E4ZChZ3FMNTu+NuzFTTsKcn
oLhcQgzrur1/Dx2lISMrSLI3QxQv9K+uj47iApo0a2qkimKXyfL0v+K4MPbO
cJlnkijmR3S/6EjkgovFF4838aJYneQjYp3CBfUv0RX/hIiiu8bvx+3ymPO0
8yGr52tFMRs0237c4kK7QlHsfTlRrNTuKGCzeaB25sVg6WkW7s/q9S8r4EHV
sXmt2SdZ2Oi2Zr5WCQ8sUkJxXw4L24wMJUoreOAbzju5JYOF+iaeVpev8OCc
nZmd/B4WvvcJCWA18sBovPT8QQ8WFupaJPtO8MDRKSLUR5mFM4nX1W1c+ZDA
kpvSODkLI3d/vhWpJg5OMjVSoyUiOJbVM5ZZKwE3/NUiZUNmhEVSN4dv20pB
qHuYRXXAtPBBcbXmB19p4JUcFUwmTQn7KpoPJ6cLwLzAxe3v2AlhemORP/eD
DLiPDRSRY+PC9W0B/LB1cvDO943C9+p/hYITW4XqzfIglxLLcj7/RWhXcbPn
bIoizGP3j7o2jAqNl/GvOFrNhjNevX8fCfssHKBmVR+wVoHuiebGk9uHhXUk
aa3+H6owaLnu37GyQeGPZ58jU5erg+5HVVorp08YV7nOocVPAzYfEbN3yOgS
HtZwXDSWpQmWD0IWn5XvEPJ/G4u2vT8XNNLtpCWq24X//z/+D5wumZI=
     "]]}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesOrigin->{0, 0},
  PlotRange->{{-6, 0}, {0, 10}},
  PlotRangeClipping->True,
  PlotRangePadding->{
    Scaled[0.02], Automatic}]], "Output",
 CellChangeTimes->{{3.397747746228539*^9, 3.397747755947188*^9}, 
   3.397782687519788*^9, 3.3982768780084248`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"DampedNewton", "[", 
  RowBox[{"f", ",", "32", ",", "0.0000001", ",", "100"}], "]"}]], "Input",
 CellChangeTimes->{{3.397746747738326*^9, 3.3977467895511913`*^9}, {
  3.3977470821960907`*^9, 3.3977470877322817`*^9}, {3.397747235140584*^9, 
  3.397747235517539*^9}, {3.397747388576626*^9, 3.3977473891066647`*^9}, {
  3.397747499353001*^9, 3.397747503636622*^9}, {3.397747636415607*^9, 
  3.397747645387024*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "1", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   "14.457142857142859`"}],
  SequenceForm["iteration:", 1, " p=", 14.457142857142859`],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.397746790041247*^9, {3.3977470881158447`*^9, 3.397747108319227*^9}, {
   3.397747223419484*^9, 3.397747236311907*^9}, {3.397747293357253*^9, 
   3.397747312592518*^9}, {3.397747366127379*^9, 3.397747404509283*^9}, {
   3.3977474816462193`*^9, 3.397747504036358*^9}, 3.3977475644588337`*^9, {
   3.39774763082836*^9, 3.397747645761821*^9}, 3.398276878034697*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "2", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   "5.642646714987142`"}],
  SequenceForm["iteration:", 2, " p=", 5.642646714987142],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.397746790041247*^9, {3.3977470881158447`*^9, 3.397747108319227*^9}, {
   3.397747223419484*^9, 3.397747236311907*^9}, {3.397747293357253*^9, 
   3.397747312592518*^9}, {3.397747366127379*^9, 3.397747404509283*^9}, {
   3.3977474816462193`*^9, 3.397747504036358*^9}, 3.3977475644588337`*^9, {
   3.39774763082836*^9, 3.397747645761821*^9}, 3.398276878043386*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "3", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   "1.1477654128654677`"}],
  SequenceForm["iteration:", 3, " p=", 1.1477654128654677`],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.397746790041247*^9, {3.3977470881158447`*^9, 3.397747108319227*^9}, {
   3.397747223419484*^9, 3.397747236311907*^9}, {3.397747293357253*^9, 
   3.397747312592518*^9}, {3.397747366127379*^9, 3.397747404509283*^9}, {
   3.3977474816462193`*^9, 3.397747504036358*^9}, 3.3977475644588337`*^9, {
   3.39774763082836*^9, 3.397747645761821*^9}, 3.39827687804615*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "4", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "1.287757803743499`"}]}],
  SequenceForm["iteration:", 4, " p=", -1.287757803743499],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.397746790041247*^9, {3.3977470881158447`*^9, 3.397747108319227*^9}, {
   3.397747223419484*^9, 3.397747236311907*^9}, {3.397747293357253*^9, 
   3.397747312592518*^9}, {3.397747366127379*^9, 3.397747404509283*^9}, {
   3.3977474816462193`*^9, 3.397747504036358*^9}, 3.3977475644588337`*^9, {
   3.39774763082836*^9, 3.397747645761821*^9}, 3.398276878048829*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "5", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "3.0199231923812757`"}]}],
  SequenceForm["iteration:", 5, " p=", -3.0199231923812757`],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.397746790041247*^9, {3.3977470881158447`*^9, 3.397747108319227*^9}, {
   3.397747223419484*^9, 3.397747236311907*^9}, {3.397747293357253*^9, 
   3.397747312592518*^9}, {3.397747366127379*^9, 3.397747404509283*^9}, {
   3.3977474816462193`*^9, 3.397747504036358*^9}, 3.3977475644588337`*^9, {
   3.39774763082836*^9, 3.397747645761821*^9}, 3.3982768780529203`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "6", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "2.9831560534309776`"}]}],
  SequenceForm["iteration:", 6, " p=", -2.9831560534309776`],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.397746790041247*^9, {3.3977470881158447`*^9, 3.397747108319227*^9}, {
   3.397747223419484*^9, 3.397747236311907*^9}, {3.397747293357253*^9, 
   3.397747312592518*^9}, {3.397747366127379*^9, 3.397747404509283*^9}, {
   3.3977474816462193`*^9, 3.397747504036358*^9}, 3.3977475644588337`*^9, {
   3.39774763082836*^9, 3.397747645761821*^9}, 3.398276878055881*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "7", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "3.004899507057452`"}]}],
  SequenceForm["iteration:", 7, " p=", -3.004899507057452],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.397746790041247*^9, {3.3977470881158447`*^9, 3.397747108319227*^9}, {
   3.397747223419484*^9, 3.397747236311907*^9}, {3.397747293357253*^9, 
   3.397747312592518*^9}, {3.397747366127379*^9, 3.397747404509283*^9}, {
   3.3977474816462193`*^9, 3.397747504036358*^9}, 3.3977475644588337`*^9, {
   3.39774763082836*^9, 3.397747645761821*^9}, 3.398276878058606*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "8", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "2.9955563766280977`"}]}],
  SequenceForm["iteration:", 8, " p=", -2.9955563766280977`],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.397746790041247*^9, {3.3977470881158447`*^9, 3.397747108319227*^9}, {
   3.397747223419484*^9, 3.397747236311907*^9}, {3.397747293357253*^9, 
   3.397747312592518*^9}, {3.397747366127379*^9, 3.397747404509283*^9}, {
   3.3977474816462193`*^9, 3.397747504036358*^9}, 3.3977475644588337`*^9, {
   3.39774763082836*^9, 3.397747645761821*^9}, 3.398276878061301*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "9", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "3.000707203403132`"}]}],
  SequenceForm["iteration:", 9, " p=", -3.000707203403132],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.397746790041247*^9, {3.3977470881158447`*^9, 3.397747108319227*^9}, {
   3.397747223419484*^9, 3.397747236311907*^9}, {3.397747293357253*^9, 
   3.397747312592518*^9}, {3.397747366127379*^9, 3.397747404509283*^9}, {
   3.3977474816462193`*^9, 3.397747504036358*^9}, 3.3977475644588337`*^9, {
   3.39774763082836*^9, 3.397747645761821*^9}, 3.398276878064001*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "10", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "2.9996958170701773`"}]}],
  SequenceForm["iteration:", 10, " p=", -2.9996958170701773`],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.397746790041247*^9, {3.3977470881158447`*^9, 3.397747108319227*^9}, {
   3.397747223419484*^9, 3.397747236311907*^9}, {3.397747293357253*^9, 
   3.397747312592518*^9}, {3.397747366127379*^9, 3.397747404509283*^9}, {
   3.3977474816462193`*^9, 3.397747504036358*^9}, 3.3977475644588337`*^9, {
   3.39774763082836*^9, 3.397747645761821*^9}, 3.39827687806668*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "11", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "3.0002836670851187`"}]}],
  SequenceForm["iteration:", 11, " p=", -3.0002836670851187`],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.397746790041247*^9, {3.3977470881158447`*^9, 3.397747108319227*^9}, {
   3.397747223419484*^9, 3.397747236311907*^9}, {3.397747293357253*^9, 
   3.397747312592518*^9}, {3.397747366127379*^9, 3.397747404509283*^9}, {
   3.3977474816462193`*^9, 3.397747504036358*^9}, 3.3977475644588337`*^9, {
   3.39774763082836*^9, 3.397747645761821*^9}, 3.398276878069365*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "12", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "2.999968484343618`"}]}],
  SequenceForm["iteration:", 12, " p=", -2.999968484343618],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.397746790041247*^9, {3.3977470881158447`*^9, 3.397747108319227*^9}, {
   3.397747223419484*^9, 3.397747236311907*^9}, {3.397747293357253*^9, 
   3.397747312592518*^9}, {3.397747366127379*^9, 3.397747404509283*^9}, {
   3.3977474816462193`*^9, 3.397747504036358*^9}, 3.3977475644588337`*^9, {
   3.39774763082836*^9, 3.397747645761821*^9}, 3.398276878072195*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "13", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "3.0000128110060698`"}]}],
  SequenceForm["iteration:", 13, " p=", -3.0000128110060698`],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.397746790041247*^9, {3.3977470881158447`*^9, 3.397747108319227*^9}, {
   3.397747223419484*^9, 3.397747236311907*^9}, {3.397747293357253*^9, 
   3.397747312592518*^9}, {3.397747366127379*^9, 3.397747404509283*^9}, {
   3.3977474816462193`*^9, 3.397747504036358*^9}, 3.3977475644588337`*^9, {
   3.39774763082836*^9, 3.397747645761821*^9}, 3.398276878074957*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "14", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "2.999999180305889`"}]}],
  SequenceForm["iteration:", 14, " p=", -2.999999180305889],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.397746790041247*^9, {3.3977470881158447`*^9, 3.397747108319227*^9}, {
   3.397747223419484*^9, 3.397747236311907*^9}, {3.397747293357253*^9, 
   3.397747312592518*^9}, {3.397747366127379*^9, 3.397747404509283*^9}, {
   3.3977474816462193`*^9, 3.397747504036358*^9}, 3.3977475644588337`*^9, {
   3.39774763082836*^9, 3.397747645761821*^9}, 3.398276878077772*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "15", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "3.0000000124712267`"}]}],
  SequenceForm["iteration:", 15, " p=", -3.0000000124712267`],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.397746790041247*^9, {3.3977470881158447`*^9, 3.397747108319227*^9}, {
   3.397747223419484*^9, 3.397747236311907*^9}, {3.397747293357253*^9, 
   3.397747312592518*^9}, {3.397747366127379*^9, 3.397747404509283*^9}, {
   3.3977474816462193`*^9, 3.397747504036358*^9}, 3.3977475644588337`*^9, {
   3.39774763082836*^9, 3.397747645761821*^9}, 3.398276878080484*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "16", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "3.0000000124712267`"}]}],
  SequenceForm["iteration:", 16, " p=", -3.0000000124712267`],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.397746790041247*^9, {3.3977470881158447`*^9, 3.397747108319227*^9}, {
   3.397747223419484*^9, 3.397747236311907*^9}, {3.397747293357253*^9, 
   3.397747312592518*^9}, {3.397747366127379*^9, 3.397747404509283*^9}, {
   3.3977474816462193`*^9, 3.397747504036358*^9}, 3.3977475644588337`*^9, {
   3.39774763082836*^9, 3.397747645761821*^9}, 3.3982768780831947`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"16", ",", 
   RowBox[{"-", "3.0000000124712267`"}]}], "}"}]], "Output",
 CellChangeTimes->{3.397782687602878*^9, 3.398276878085896*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "16", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "3.0000000124712267`"}]}],
  SequenceForm["iteration:", 16, " p=", -3.0000000124712267`],
  Editable->False]], "Print",
 GeneratedCell->False,
 CellAutoOverwrite->False,
 CellChangeTimes->{
  3.397746790041247*^9, {3.3977470881158447`*^9, 3.397747108319227*^9}, {
   3.397747223419484*^9, 3.397747236311907*^9}, {3.397747293357253*^9, 
   3.397747312592518*^9}, {3.397747366127379*^9, 3.397747404509283*^9}, {
   3.3977474816462193`*^9, 3.397747504036358*^9}, 3.3977475644588337`*^9, {
   3.39774763082836*^9, 3.397747649220955*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"16", ",", 
   RowBox[{"-", "3.0000000124712267`"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.397746758028408*^9, 3.3977467907583303`*^9}, {
   3.397747088286936*^9, 3.397747108791889*^9}, {3.3977472234247828`*^9, 
   3.397747240797082*^9}, {3.397747294683403*^9, 3.397747313699307*^9}, {
   3.397747366867434*^9, 3.397747415033289*^9}, {3.397747482482562*^9, 
   3.397747507870969*^9}, 3.397747564461244*^9, {3.397747631287345*^9, 
   3.397747645785097*^9}}]
}, Open  ]]
}, Closed]],

Cell[CellGroupData[{

Cell["DampledNewtonPlot - uses Version 6", "Section",
 CellChangeTimes->{{3.397747990707266*^9, 3.397747992244878*^9}, {
  3.397748869777343*^9, 3.3977488731774607`*^9}, {3.398277503532343*^9, 
  3.398277504644359*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"DampedNewtonPlot", "[", 
   RowBox[{"f_", ",", " ", "p0_", ",", 
    RowBox[{"\[Epsilon]_:", "0.001"}], ",", 
    RowBox[{"Nmax_:", "10"}]}], "]"}], ":=", "\[IndentingNewLine]", " ", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
     "\[Delta]", ",", " ", "\[CapitalDelta]", ",", "i", ",", "p", ",", "j", 
      ",", "  ", "pnew", ",", " ", "fnew", ",", " ", "fold", ",", " ", 
      "pdata", ",", " ", "p1", ",", " ", "p2", ",", " ", "urange", ",", " ", 
      "umin", ",", " ", "umax"}], "}"}], ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{"\[CapitalDelta]", "[", "x_", "]"}], ":=", 
      RowBox[{
       RowBox[{"f", "[", "x", "]"}], "/", 
       RowBox[{
        RowBox[{"f", "'"}], "[", "x", "]"}]}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"p", "=", 
      RowBox[{"1.0", "*", "p0"}]}], ";", " ", 
     RowBox[{"fnew", "=", 
      RowBox[{"f", "[", "p", "]"}]}], ";", " ", "\[IndentingNewLine]", 
     RowBox[{"\[Delta]", "=", 
      RowBox[{"\[CapitalDelta]", "[", "p", "]"}]}], ";", " ", 
     "\[IndentingNewLine]", 
     RowBox[{"i", "=", "0"}], ";", "\[IndentingNewLine]", 
     RowBox[{"pdata", "=", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"p", ",", " ", "0"}], "}"}], ",", " ", 
        RowBox[{"{", 
         RowBox[{"p", ",", " ", 
          RowBox[{"f", "[", "p", "]"}]}], "}"}]}], "}"}]}], ";", " ", 
     "\[IndentingNewLine]", 
     RowBox[{"Print", "[", 
      RowBox[{"\[Delta]", ",", " ", "\"\< \>\"", ",", " ", "i"}], "]"}], ";", 
     " ", "\[IndentingNewLine]", 
     RowBox[{"While", "[", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{"Abs", "[", "\[Delta]", "]"}], ">", "\[Epsilon]"}], " ", 
        "\[And]", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"i", "++"}], "<", "Nmax"}], ")"}]}], ",", 
       "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"fold", "=", "fnew"}], ";", " ", "\[IndentingNewLine]", 
        RowBox[{"pnew", "=", 
         RowBox[{"p", "-", "\[Delta]"}]}], ";", "\[IndentingNewLine]", 
        RowBox[{"fnew", "=", 
         RowBox[{"f", "[", "pnew", "]"}]}], ";", " ", "\[IndentingNewLine]", 
        RowBox[{"j", "=", "0"}], ";", " ", "\[IndentingNewLine]", 
        RowBox[{"While", "[", 
         RowBox[{
          RowBox[{
           RowBox[{
            RowBox[{"Abs", "[", "fnew", "]"}], "\[GreaterEqual]", " ", 
            RowBox[{"Abs", "[", "fold", "]"}]}], " ", "\[And]", " ", 
           RowBox[{
            RowBox[{"j", "++"}], "<", "Nmax"}]}], ",", " ", 
          "\[IndentingNewLine]", 
          RowBox[{
           RowBox[{"\[Delta]", "=", 
            RowBox[{"\[Delta]", "/", "2"}]}], ";", " ", "\[IndentingNewLine]",
            
           RowBox[{"pnew", "=", 
            RowBox[{"p", "-", "\[Delta]"}]}], ";", " ", "\[IndentingNewLine]",
            
           RowBox[{"fnew", "=", 
            RowBox[{"f", "[", "pnew", "]"}]}], ";", "\[IndentingNewLine]", 
           RowBox[{"(*", " ", 
            RowBox[{"Print", "[", 
             RowBox[{
             "\"\<Revising \[Delta]=\>\"", ",", " ", "\[Delta]", ",", " ", 
              "\"\< \!\(\*SubscriptBox[\(p\), \(new\)]\)=\>\"", ",", " ", 
              "pnew"}], "]"}], "*)"}], ";"}]}], " ", "\[IndentingNewLine]", 
         "]"}], ";", "\[IndentingNewLine]", 
        RowBox[{"Print", "[", 
         RowBox[{
         "\"\<iteration:\>\"", ",", "i", ",", " ", "\"\< p=\>\"", ",", 
          "pnew"}], "]"}], ";", "\[IndentingNewLine]", 
        RowBox[{"AppendTo", "[", 
         RowBox[{"pdata", ",", " ", 
          RowBox[{"{", 
           RowBox[{"pnew", ",", " ", "0"}], "}"}]}], "]"}], ";", " ", 
        "\[IndentingNewLine]", 
        RowBox[{"AppendTo", "[", 
         RowBox[{"pdata", ",", " ", 
          RowBox[{"{", 
           RowBox[{"pnew", ",", " ", "fnew"}], "}"}]}], "]"}], ";", " ", 
        "\[IndentingNewLine]", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"Abs", "[", 
            RowBox[{"p", "-", "pnew"}], "]"}], "<", "\[Epsilon]"}], ",", " ", 
          
          RowBox[{"Break", "[", "]"}]}], "]"}], ";", " ", 
        "\[IndentingNewLine]", "\[IndentingNewLine]", 
        RowBox[{"p", "=", "pnew"}], ";", "  ", "\[IndentingNewLine]", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{
           RowBox[{
            RowBox[{"f", "'"}], "[", "p", "]"}], "\[Equal]", "0"}], ",", " ", 
          
          RowBox[{"Break", "[", "]"}]}], "]"}], ";", " ", 
        "\[IndentingNewLine]", 
        RowBox[{"\[Delta]", "=", 
         RowBox[{"\[CapitalDelta]", "[", "p", "]"}]}], ";"}]}], 
      "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", 
     "\[IndentingNewLine]", 
     RowBox[{"p1", "=", 
      RowBox[{"ListPlot", "[", 
       RowBox[{"pdata", ",", " ", 
        RowBox[{"PlotJoined", "\[Rule]", " ", "True"}], ",", " ", 
        RowBox[{"PlotRange", "\[Rule]", " ", "All"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", "\[IndentingNewLine]", 
     RowBox[{"urange", " ", "=", " ", 
      RowBox[{"First", "/@", "pdata"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"umin", "=", 
      RowBox[{"Min", "[", "urange", "]"}]}], ";", " ", 
     RowBox[{"umax", "=", 
      RowBox[{"Max", "[", "urange", "]"}]}], ";", " ", 
     RowBox[{"urange", "=", 
      RowBox[{"umax", "-", "umin"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"umin", " ", "=", " ", 
      RowBox[{"umin", "-", 
       RowBox[{"urange", "/", "4"}]}]}], ";", " ", 
     RowBox[{"umax", " ", "=", " ", 
      RowBox[{"umax", " ", "+", " ", 
       RowBox[{"urange", "/", "4"}]}]}], ";", " ", "\[IndentingNewLine]", 
     RowBox[{"p2", "=", 
      RowBox[{"Plot", "[", 
       RowBox[{
        RowBox[{"f", "[", "u", "]"}], ",", " ", 
        RowBox[{"{", 
         RowBox[{"u", ",", " ", "umin", ",", " ", "umax"}], "}"}]}], "]"}]}], 
     ";", " ", "\[IndentingNewLine]", 
     RowBox[{"Print", "[", 
      RowBox[{"Show", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"p1", ",", "p2"}], "}"}], ",", " ", 
        RowBox[{"PlotRange", "\[Rule]", " ", "All"}]}], "]"}], "]"}]}]}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", "]"}]}]], "Input",
 CellChangeTimes->{{3.397746800560541*^9, 3.397747103594524*^9}, {
   3.3977471748801823`*^9, 3.397747218537197*^9}, {3.397747265513595*^9, 
   3.397747309650807*^9}, {3.397747354940393*^9, 3.397747362265942*^9}, 
   3.397747400976905*^9, {3.397747466962101*^9, 3.3977474766165867`*^9}, {
   3.397747541035405*^9, 3.397747627845809*^9}, {3.397747994890922*^9, 
   3.397748073553388*^9}, {3.397748150894238*^9, 3.397748303145853*^9}, {
   3.397748374854624*^9, 3.397748380274761*^9}, {3.3977484512273483`*^9, 
   3.3977485952101173`*^9}, {3.397748727147669*^9, 3.3977487351522503`*^9}, {
   3.398277007570841*^9, 3.3982770109993477`*^9}, {3.39827707642953*^9, 
   3.39827710777302*^9}, {3.398277144234272*^9, 3.398277155160398*^9}, {
   3.398277466261368*^9, 3.398277468009128*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"f", "[", "x_", "]"}], ":=", " ", 
  RowBox[{
   RowBox[{
    RowBox[{"(", 
     RowBox[{"x", "+", "3"}], ")"}], "^", "2"}], "+", 
   "3"}]}], "\[IndentingNewLine]", 
 RowBox[{"DampedNewtonPlot", "[", 
  RowBox[{"f", ",", "1.5"}], "]"}]}], "Input",
 CellChangeTimes->{{3.397748098501985*^9, 3.3977481159294977`*^9}, {
  3.397748331958654*^9, 3.3977483365861273`*^9}, {3.3982774731924057`*^9, 
  3.398277474261052*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 InterpretationBox[
  RowBox[{
  "2.583333333333333`", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]",
    "0"}],
  SequenceForm[2.583333333333333, " ", 0],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.397748101211018*^9, {3.397748284925549*^9, 3.39774830595297*^9}, 
   3.397748337334689*^9, 3.397748384259513*^9, {3.3977484588305693`*^9, 
   3.397748508080845*^9}, {3.3977485505151777`*^9, 3.3977485988749123`*^9}, 
   3.397782687677505*^9, 3.398276878166772*^9, 3.398277112162354*^9, 
   3.39827747462904*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "1", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "1.083333333333333`"}]}],
  SequenceForm["iteration:", 1, " p=", -1.083333333333333],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.397748101211018*^9, {3.397748284925549*^9, 3.39774830595297*^9}, 
   3.397748337334689*^9, 3.397748384259513*^9, {3.3977484588305693`*^9, 
   3.397748508080845*^9}, {3.3977485505151777`*^9, 3.3977485988749123`*^9}, 
   3.397782687677505*^9, 3.398276878166772*^9, 3.398277112162354*^9, 
   3.398277474631856*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "2", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "2.82427536231884`"}]}],
  SequenceForm["iteration:", 2, " p=", -2.82427536231884],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.397748101211018*^9, {3.397748284925549*^9, 3.39774830595297*^9}, 
   3.397748337334689*^9, 3.397748384259513*^9, {3.3977484588305693`*^9, 
   3.397748508080845*^9}, {3.3977485505151777`*^9, 3.3977485988749123`*^9}, 
   3.397782687677505*^9, 3.398276878166772*^9, 3.398277112162354*^9, 
   3.3982774746347027`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "3", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "3.0937736371021947`"}]}],
  SequenceForm["iteration:", 3, " p=", -3.0937736371021947`],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.397748101211018*^9, {3.397748284925549*^9, 3.39774830595297*^9}, 
   3.397748337334689*^9, 3.397748384259513*^9, {3.3977484588305693`*^9, 
   3.397748508080845*^9}, {3.3977485505151777`*^9, 3.3977485988749123`*^9}, 
   3.397782687677505*^9, 3.398276878166772*^9, 3.398277112162354*^9, 
   3.398277474637488*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "4", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "2.968438842024392`"}]}],
  SequenceForm["iteration:", 4, " p=", -2.968438842024392],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.397748101211018*^9, {3.397748284925549*^9, 3.39774830595297*^9}, 
   3.397748337334689*^9, 3.397748384259513*^9, {3.3977484588305693`*^9, 
   3.397748508080845*^9}, {3.3977485505151777`*^9, 3.3977485988749123`*^9}, 
   3.397782687677505*^9, 3.398276878166772*^9, 3.398277112162354*^9, 
   3.398277474640183*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "5", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "3.0148671172908936`"}]}],
  SequenceForm["iteration:", 5, " p=", -3.0148671172908936`],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.397748101211018*^9, {3.397748284925549*^9, 3.39774830595297*^9}, 
   3.397748337334689*^9, 3.397748384259513*^9, {3.3977484588305693`*^9, 
   3.397748508080845*^9}, {3.3977485505151777`*^9, 3.3977485988749123`*^9}, 
   3.397782687677505*^9, 3.398276878166772*^9, 3.398277112162354*^9, 
   3.398277474642929*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "6", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "2.916330753669381`"}]}],
  SequenceForm["iteration:", 6, " p=", -2.916330753669381],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.397748101211018*^9, {3.397748284925549*^9, 3.39774830595297*^9}, 
   3.397748337334689*^9, 3.397748384259513*^9, {3.3977484588305693`*^9, 
   3.397748508080845*^9}, {3.3977485505151777`*^9, 3.3977485988749123`*^9}, 
   3.397782687677505*^9, 3.398276878166772*^9, 3.398277112162354*^9, 
   3.398277474645718*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "7", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "3.056718009266565`"}]}],
  SequenceForm["iteration:", 7, " p=", -3.056718009266565],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.397748101211018*^9, {3.397748284925549*^9, 3.39774830595297*^9}, 
   3.397748337334689*^9, 3.397748384259513*^9, {3.3977484588305693`*^9, 
   3.397748508080845*^9}, {3.3977485505151777`*^9, 3.3977485988749123`*^9}, 
   3.397782687677505*^9, 3.398276878166772*^9, 3.398277112162354*^9, 
   3.398277474648431*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "8", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "2.9533000976140547`"}]}],
  SequenceForm["iteration:", 8, " p=", -2.9533000976140547`],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.397748101211018*^9, {3.397748284925549*^9, 3.39774830595297*^9}, 
   3.397748337334689*^9, 3.397748384259513*^9, {3.3977484588305693`*^9, 
   3.397748508080845*^9}, {3.3977485505151777`*^9, 3.3977485988749123`*^9}, 
   3.397782687677505*^9, 3.398276878166772*^9, 3.398277112162354*^9, 
   3.3982774746511593`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "9", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "3.016080041825876`"}]}],
  SequenceForm["iteration:", 9, " p=", -3.016080041825876],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.397748101211018*^9, {3.397748284925549*^9, 3.39774830595297*^9}, 
   3.397748337334689*^9, 3.397748384259513*^9, {3.3977484588305693`*^9, 
   3.397748508080845*^9}, {3.3977485505151777`*^9, 3.3977485988749123`*^9}, 
   3.397782687677505*^9, 3.398276878166772*^9, 3.398277112162354*^9, 
   3.398277474653879*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "10", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "2.924975179062127`"}]}],
  SequenceForm["iteration:", 10, " p=", -2.924975179062127],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.397748101211018*^9, {3.397748284925549*^9, 3.39774830595297*^9}, 
   3.397748337334689*^9, 3.397748384259513*^9, {3.3977484588305693`*^9, 
   3.397748508080845*^9}, {3.3977485505151777`*^9, 3.3977485988749123`*^9}, 
   3.397782687677505*^9, 3.398276878166772*^9, 3.398277112162354*^9, 
   3.398277474656664*^9}],

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    {Hue[0.67, 0.6, 0.6], 
     LineBox[{{1.5, 0.}, {1.5, 23.25}, {-1.083333333333333, 
      0.}, {-1.083333333333333, 6.6736111111111125`}, {-2.82427536231884, 
      0.}, {-2.82427536231884, 3.030879148288175}, {-3.0937736371021947`, 
      0.}, {-3.0937736371021947`, 3.0087934950153743`}, {-2.968438842024392, 
      0.}, {-2.968438842024392, 3.0009961066927615`}, {-3.0148671172908936`, 
      0.}, {-3.0148671172908936`, 3.000221031176541}, {-2.916330753669381, 
      0.}, {-2.916330753669381, 3.0070005427815336`}, {-3.056718009266565, 
      0.}, {-3.056718009266565, 3.003216932575162}, {-2.9533000976140547`, 
      0.}, {-2.9533000976140547`, 3.002180880882857}, {-3.016080041825876, 
      0.}, {-3.016080041825876, 3.000258567745122}, {-2.924975179062127, 
      0.}, {-2.924975179062127, 3.00562872375676}}]}}, {{}, {}, 
    {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJwV13lYjO0XB/C2mbmfmSnTRqn8tKBQKrSp5hQlhBIhVCIlSouUJKJIWols
SYUiLZK16LRbkpBooaKdXtHbJvK737/m+lzPM895rvs+53vPqLrtWeMuJiIi
woiKiPz3KTfmdpczKsBnj/Nnqs+VgZDEscKPQwKs4kv+KJ8pA21aCYV3BgX4
Udgi46EqA7ecHt9x+S5ApzVvNxdOloElxVMK7ncJMEFTg7tLVAYCD9XkerwT
YN95u+c276WhgWWQVV0gwLUJ9W2VEdJwXoa5eNxbgKyaqB1VXQIIjxBYlOVO
QgX+hmOPXSeBcnLky01npTDmm0j26d2SMBsECwSnpHDsyDWVo+6SYNR74UJF
rBRWvJPWDnKWhHWL8j20I6UwIN9CzdtOEuJam0T/BkihvsmEe8ICSRDR1DVI
t5dCw9S+NYV/+NDxoOlyD18Kt+YllPmd5kN20zzfwKOSOLdHeqyshgdetg8z
asIkcXeyMH5jFQ9mP7F8rx4iiZEBjaeGSniQlbbO/I2vJBoU2H82u8OD654H
+POcJTFr2dr0SRd4cGW4KqvHSBJL29Ku3/fkQZKMc/umf/hYm+eMilweOER0
y93p5eOnMWkJZwkeyAz72nA7+Vjhq3c3c4ILiY1H8x4081Gus9BlxSAX4q7c
OCj/jI9l3S/fvW3hwgmdIcVXGXx0XGdj3ZvPhbDlMWssN1K/djaauYUL3GSl
6b1r+WiwbKmTzwYuJH+52R9vx8fEGxsNihy4kB/6LOrjUj4Wqv9W3racC59z
2SXBBny8+8fw6ldDLljLhs/Jk+Vjis3JC62yXJBqCZRQruUhz+TKTIdXDFzQ
ZL0pe8ZD3dLKWP0XDMwKTErdWcnDV7ryeZOrGQCpOyb3i3kY+kxvev8TBvxh
wHdNNg8T2rYWv8lloOGq18eoKB6mXktpORzPQIq3y70hSx5uPHm02HMtA01H
+u2E5jx0SZOsP23HgMK5A1+jjHlYXdgaXGHLQFJZ8nRlXR4+fLvJy8CKgegp
ddGWKjzU+Hnn22YDBvaVgkv8CBen3f7iU6jIwGp5NaJ1i4tphhY7YzsIxMzO
T/fP5CKOPdfyaifwXGhuVpzOReHDS0G2nwhYezn5rz7PRVddhSLVDwRMS061
BB7j4uy/Pu1SNQQ0d4rfLnPl4udpl7IHCwmIPu7YuFmei207/l3YHU1ALPPF
v0UCLiYMBo8fjiIgkVgQr8TnIqztnjr9GAG2++HKJjEuJieGe/iEE+BLKutt
HGBw3u9swxXBBKZscSCOLxgUT/xwh+VBQGGpScbdKgbPbPk1ZWg7gal6quby
ZQzOmGA19LoRUGZ9D6h/wKB5UZPLF2cCajnRrWuuM7hpRsjNyesJaP8pvbc6
nEGDhwlKYTYEdLqz7PNCGbxRuXpauTUB3dfx36SCGZx6fTRE0oqA/rUtarU+
DA5/5lgWWBAwWjkWa7uZwbO2yhzTRQRMDNs0b61n8NeoVmO2MYFFqtXlPAda
b78qS9WIgPlQ0tjzZQwqDL/0UlpIYHGKrvsyQwZDWpW97eYRWHJ8isgNfQZr
5+yX+qpN19tv4gLRYXC6bQlzci4BG6uaumoNBp8QRqpZi8Cqfg9TaxkGP/ud
WlmkQWCjWars4n6CteIHs0KV6H6Ixdj87CE4w1VR9NBUAsNVwQfTOgiaVGWL
HlMkkGS3pvtvM8FLeQW/0qYQqN3GKip+TtByxKZGWo7AJc2f33dVEYxMO7Pd
WJbArv5PGkpl1Os+XHWXIUCCHsTtf0iwpEmoUiMgYBm9y80gi6C8U3J/nSQB
weoNyZ0ZBK3WJFWqULfKWtUkpRLUX3vcag+fQGjKNIPBswSXBK42m84jsNyN
tyv9FMGdSRp3D3EJKM4aTbWPI2j7+Hv9F4bAvfzXTH4kweVzPfmPCIGIwCfm
LuEE3f3lNeZSO5hkB0gdJGhe9uBLOofAQHnEx917CWodane6zCYw573hba1t
BNu8izRfShD4dVGj64MzwW5RxcPbqZ+6SitFORG8lbr7kAj1jr6+iC57gntG
hwKtxQksyHv/8MxKgrvvzPEeECMgvrfinyXLCM4KXs9Ppb7yJ2VDBhBsmJZo
yaX2KYuOXWNK8JlIamOVKAGz40FlokYEQ76lK0RR8223j+TPJ9jbfunvSuom
gf1c13l0/dpj4hSos96ZbZ00h+DfXv/KbhECQRdmn30yk+D136uuFlFbu0x5
4a1GkKuoOieJWk5DQkRlGkFtYe9GP+rPPQMLahQJ8vdk6jpQ3875uPOAPMEb
WZvyjKkP+z+/PFuaoMQ3VoMG9SrD+28b+QSVja9flaNW/p1BThCCP+NNlRjq
r5hgZiRB93/gqbEY9aPIg/7dfznY4LSCnpYEopZ7ZZ4d52BHbbnPf9cdJ61v
sRrhYL6tfuR/359Rv1h66CcHjevPLpan/vecrvXVfzh42P3HjRnUZVtUDjj0
cTBCxPKBCXWiGjdfrIuD1teifNZSu3YPd9xu52ClQ2WNP7XOrS+KWz9yUIw/
+uYM9R/fulWCRg7+eTk9/DF16ESRnVs9B8204EU/9e+TmWsKX3FwPftbnDpd
31DF02vZLzi4zF2jxpn69/Uwxw1VHJSz7j2S+t/1BV4bbpZy8EGWyZ1O6vHS
dU6/izm4MFZhix7d318tc52v3OHgw+b9Ko3UIV4Krj9zOVhq9Wf+QtovYyPi
bktucjBNT7LiLPWoTLN7zxUOfot7Y+dF+yv4SpWHySX6PiJRK1uoR7ULdsYk
c/DC4LNyBxaBEZsT3rpxHLzlMqjsQPt3X8PePUdOcDDkQ9PfJurhba5+9REc
VH27arMn7ffhMMPA4AMcHFvZEn6azsM+vnrQ830cnF07tFyPzs/Qean9yv70
/esuxtZT/1vYGYqeHKzfGHtoDp2/vZavw2S2c9Cz/YNOK/Xgq+LD213o/vbl
rz9H5/Vn3+kI4shB06PsyUpSBH5Mt4yxW8xBvedLrdykCfjmaselm3Mw1nUE
LWg+DCxSTPjXmK5H8Px7M2l+DDgOnD6ny8FCHXtJUXkC/8RcvtimwsGKfDWf
EQU6L1OjU/QVOdi1zuCeGM2r/szA1Ag5DraHNgTJ0jz7VmabocnjoJtdsYeV
Cu3H0V83fUfZGCs57VefKgGvY123ygbZuK/0jJSaOoE+2Te5ct/ZaDRy/pIz
zcc+nRsFDzrZaNWzdEbPTAI929c/EnnLxpNtxr2mNF876+48S8hhY4nWvkX9
hjTvja+5699gY+G228KtNM/D0s6K1l9l4yp1T6bJhJ4//vuNp1xio/d8m/5G
MwK2ssIbKdFs3Jz8ZlhmCYGCdc+jbnqwUcRHIO1lT+Dv4yIN221sNBST8J7t
QO+fmYP9zmyUUQ6I+76W1h+OH9V1pPXaylWjNtDnn1vneX8JG2/mzIv+6ULg
UFPb0gpVNm5y1HW7tYdAjeWbL+4qbHy7pEX+oR+9P7v8EEeRjQMPxOOfB9D6
B67fXy5g42T9ytfjQQS6lHfPej3BQhlPmfPJh2l9lxH2p2YW1n1+p1N4iuZt
h2TVyBkWusiKBoU/IjBbdYNbXiILld59D1F8TM835/SJHbEszJCukLtXQmDz
BwPDhgj6vHFl7p8KApdrXLIKA1gYEJ8wXvqKgMbd29F+9iw0eRI2NquLzuux
dau/8ln4ZX7Q1u2TGTAvT/2aTlj4a5W8Wvx/v0dE+o47SbCwueNc8BNlBnxD
wvDZuAROHLbrnq3OwB2fTL2sPgmsDAl/YaLLgNH6MVn3pxK4aW9frsEKBiw1
Uz58OiqBVw59aL8cwYCETpeF8SEJdP2aOpYTxUDV/HnZp0MkcFHOpqCyGAaW
CUsP2vhJoOX8/KvjSQzYOXaqFThL4Cmbl1W51xjYEqHtHWksgWbn3msqPGUg
uO2x6Nzv4nj/xdztIwIu5J37NHu/kzguOdp4v7uAC+b7tXXUHcXRJM/ywOB9
LtRsDNV7aS+O+nkff4g/5kLf1KlGqsvEUXPOjVH9ai7MTFln9cyQWvbs7mfN
XLicVuOiKC+OkRnV/e9ZPIi5+Sjp4SsxzK08se7AFh54Fp+dGFsihoWToh81
T+VDlntHrTWIoZsH1/iCKh96pPRTTy8SwzoZKzNnTT54bn0p1NYXw97tvQ3D
C6nZEuGu/xPDKcP9q7zs+eBh5ydePSaKrjVd9k4n+LCjYwVJyhXFyzbJGc70
/4A7X1RGR0EUaw8NBo4OSULOVIugiva/JZ1sD/7cSQLQPOr9cuzRn5JHvLp+
FUNpmLtkIv2M+3jJRv/b3JE0GXi6K/TAhOFYSejds3I/leSg/6vz25fNwyXo
FpS9qU0efr3MXp8S/W+JzWnHk5+rp8CVyV2fjLf+LLGzyfL3vqsItfcjOk2Z
gZLJ1e/2VAwrgeYh1U+KBV9LtAPNhvcbT4PINVLLH3l0l8iFne8QWTEdBhvK
C2fh55Jb4/ckijNVQd0+a7lzXUtJQoyB4j0RdYhwdNLi7HhbUqMa4tJzUgN4
1bsCtq6uKHkY1tPsP2MmHMz13dwZda4k3tMvLv/eLLDN3WjZs7hEeEQt4+vr
MC3IWWq+a9XJOuHQ3QSPdaZzYMYlU7l0xUah9c2vfXHi2mBWXd+qrdcu7GpV
GDNmzwOjgPy8gchOoTrRYPWH6oKqT/nhLsM+ocayzrG1Z/Wg6eL4CssL/wgb
BVmxkpn60H5kXuuqgB/C3h/L/5QWzQfb1t8Kk58OClfapQyu/7kA/veuZbVg
aEiY9cN6G2uaAez038C5ajcq5PtNyg70NoQXOTtPGab9El61T848X2oEZ47E
Tp149lvoP1AJccQEGiZs1Z1V/wrlzPQCLgUuAmF0fFJwlgg0xqTPX1NvCv/b
vli07JEoZFUnnty50ByKpK10NnSIQcWbg36fnwiBx69r5PwSB76CeE7nRYDJ
TVlWErNY8NjuTGeItwW0LcoGTeqJfXybJOobKTmLbKmvvE++mUNturVAL4la
2ubrnjZqt54iFQ1NFiR3Roxb+VhA3nDt0GItFhjIHp0svccCbGSHrh2dwwI0
yt6U6WsBMoEjVzKpt2ffQKRubhi7+ILaa0PxjCZqnwsTibJzWXB567QBvp8F
JKmSsAzqW7rvIv2p2+cpOZZrs2BzqeRdM38LuJmoYt9Nrdu0Wmk99d7B/9ny
dFjgqX093Jea3NewdKD+eyR+VQa1jrmOzhfqN31re0mABYym6mpx5rHgSYOR
nRp1qeh8jTnUkn907y+iXltpODWAumF0S6QP9bRZJvLJ1KtqT307Tt0dZSoo
op79qdkhjTq/z5zXSu2paVz0iDrE1oItrssCbs5NtXrq/wPe8B3T
      "]]}}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesOrigin->{0, 0},
  PlotRange->All,
  PlotRangeClipping->True,
  PlotRangePadding->{Automatic, Automatic}]], "Print",
 CellChangeTimes->{
  3.397748101211018*^9, {3.397748284925549*^9, 3.39774830595297*^9}, 
   3.397748337334689*^9, 3.397748384259513*^9, {3.3977484588305693`*^9, 
   3.397748508080845*^9}, {3.3977485505151777`*^9, 3.3977485988749123`*^9}, 
   3.397782687677505*^9, 3.398276878166772*^9, 3.398277112162354*^9, 
   3.398277474659729*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"f", "[", "x_", "]"}], ":=", 
   RowBox[{
    RowBox[{"1", "/", "x"}], "-", "3"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"DampedNewtonPlot", "[", 
  RowBox[{"f", ",", "2"}], "]"}]}], "Input",
 CellChangeTimes->{{3.3977486445314903`*^9, 3.397748665435896*^9}, {
  3.397782716895771*^9, 3.397782728436283*^9}, {3.3982771203509207`*^9, 
  3.398277121355528*^9}, {3.398277477416404*^9, 3.398277478477113*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 InterpretationBox[
  RowBox[{
  "10.`", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "0"}],
  SequenceForm[10., " ", 0],
  Editable->False]], "Print",
 CellChangeTimes->{{3.397782723720703*^9, 3.397782729114873*^9}, 
   3.397782762638966*^9, 3.398276878277772*^9, 3.3982770433322783`*^9, 
   3.398277121837738*^9, 3.3982774791247587`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "1", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", "0.75`"}],
  SequenceForm["iteration:", 1, " p=", 0.75],
  Editable->False]], "Print",
 CellChangeTimes->{{3.397782723720703*^9, 3.397782729114873*^9}, 
   3.397782762638966*^9, 3.398276878277772*^9, 3.3982770433322783`*^9, 
   3.398277121837738*^9, 3.39827747912748*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "2", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", "0.28125`"}],
  SequenceForm["iteration:", 2, " p=", 0.28125],
  Editable->False]], "Print",
 CellChangeTimes->{{3.397782723720703*^9, 3.397782729114873*^9}, 
   3.397782762638966*^9, 3.398276878277772*^9, 3.3982770433322783`*^9, 
   3.398277121837738*^9, 3.398277479130583*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "3", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", "0.3251953125`"}],
  
  SequenceForm["iteration:", 3, " p=", 0.3251953125],
  Editable->False]], "Print",
 CellChangeTimes->{{3.397782723720703*^9, 3.397782729114873*^9}, 
   3.397782762638966*^9, 3.398276878277772*^9, 3.3982770433322783`*^9, 
   3.398277121837738*^9, 3.398277479133595*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "4", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   "0.33313465118408203`"}],
  SequenceForm["iteration:", 4, " p=", 0.33313465118408203`],
  Editable->False]], "Print",
 CellChangeTimes->{{3.397782723720703*^9, 3.397782729114873*^9}, 
   3.397782762638966*^9, 3.398276878277772*^9, 3.3982770433322783`*^9, 
   3.398277121837738*^9, 3.3982774791365843`*^9}],

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    {Hue[0.67, 0.6, 0.6], 
     LineBox[{{2., 0.}, {2., -2.5}, {0.75, 0.}, {
      0.75, -1.6666666666666667`}, {0.28125, 0.}, {0.28125, 
      0.5555555555555554}, {0.3251953125, 0.}, {0.3251953125, 
      0.07507507507507505}, {0.33313465118408203`, 0.}, {0.33313465118408203`,
       0.0017892057930191108`}}]}}, {{}, {}, 
    {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJwVyX040/sfx/HNpiQlSojDKpWK5FQknPfKTaVQkhTWoRElHKrDHEWHk7uQ
QzfmTDqU6FBUIsfnS35RlpmbmdyObUiYmXtyvr8/XtfzelyvjV6BTt5yBAIh
Et//OxZxcbBQgw2faKaVNOEM1DpmS8MobLBIPWOg0zgDqsy/dxhYsMF12xLn
ZewM5BnnZaSEsKGXYnWQMT8NLbSXjLN9bKD/eHUlfWAKDMreHxhGn4HyN9/5
RJsMOq8Mlq8O58DMXjWazepx2Fb2VYrd48CCEtXPvFUCweSRHSHFHBDq6W5K
z5SAAnM8o22IA64ccVKggQQ20CMsTxg0QlCDgt/R7FGwvfvHD5rMRhi8zpgO
dx4G1vCDrvxQLsi0x9QC9ouAqhX2Jj2GC5T4XT+T+oTQb3c2OTKVC+4e+0vV
E4Wgn6950OU5F9JmTxuXCfqh2Dcjh9jLhTUPE+5UpfbB/0TMy66Hm6CzQOkW
aWUvfOthzcqvbwb9L6PHrzJ5kKR8s0myqRnyKopC7yvywBjOF3QYNUNQrkVx
BKMVrrN0PV4ebQYsabvRHfcWWHJ/hLlHNAN/E8slaXsTqH7Jvl0ibIaUtr35
z+vrwawlR82zpAUozGt2t0NLIce0U+hV1QKE16UtdK/XsIa59hWd0wKRUYf8
VBxLQOx5y8l3GP+DCuGWZyGkjp1PDtzcCgQxMo88lQlfFbVW3PizFSJdXY3U
jmYj5sG7CxkhPKAGmJGA8h4ty62rz4ziQYrpjuvOpBoUrEBgspJ5QMhiSBxF
NegIJ8DscQEPHuXzthCffEAy9+PXn/Xz4IXKQG70hk/oeOhySempNug1PnCh
17cRLRZG9Lfs4QNBK+l2SzUP0bXTAgepfEhx8vBeq9OG2HH58/P2fIhkbVon
CW1DmfS2tZt9+fCW+9ktxpCPLDYYW/+SyYcg28XE3ant6LcYYc4q+XZo7Bdr
rLLrRAtux3xsee1AXhZf1v2uFy0oaA6VhnaAf1f1nGWICOX2J2TExXQAxemb
y3S6CDlUfj/mltoBFRc0/FLeilBWiLDwe0EHxF6b9vJaFKFD3UVXbXo6wFdk
/nt8jBjFldgSmm06IfHXxY0qKQNoPe2q+qhqF7zVH3yiljyEjEs4Nnr/dAO2
y8lJ0XIEPa56OKL0rhskb72xg24jaF3jhfTJum5I0cnv9QgbQVPfpoUfhN3w
IOFIvfXrEVSxlRLtp90DZCBbcXeOoiMZQdWFCT1AZx17R1o3hl5gsVna+3qh
fNng1i5fCap7ul376VoB/GlrrZxAlqKlPTcc9DYIIDT967lhZSkywZoisykC
OB3VuNxKS4py+eEipqEA1OgBz/nGUnRrRUNhymEBbCFbGbFoUmR+OfggI1wA
SVI2LeONFBXueudj3y+AcI32l/d+nkDpr4+/nCjug0y3Qrcglgy96GDTFsr6
YKLU37wsT4bYRHsl+ao+uOS6mz5XLENkB/uL6zl9MGvnU+JdK0PBYvsfzIb7
4KdL6/YLxmTIQd0x7sbmfrA68othCEyiZWEnPVek9UOsic/ggfZJ9KuFi4ru
NSHsDtYInpqZQiGXR7SVjcQwqqi6lDg8g1Z98a4xqR8AouAvjkv4HKqoqjGx
jB8Cwp2VzJmJeaQwRLYJ3zkMhOwX5i3bFpFf3E6dB2PfoL+2XdXa4TsScXse
hxWMghdD3fecwxL6/rzc86OeBDKv7aw9oErA0v8JPWMjPw7OD609PzMI2Lz9
Nq0u3jiUwSxdv5mABVRn7j2TI4VAZ+79nzYRMf+QmBLq7xOwOJ2jURRExPbX
0Ai2p2RQVW9r1FRCxJJPFhSX750Exl9Lvh+miNhEcWW5iuoUGJ0tXe24Rw7b
+qo6Jbp3CuaiTZ1yL8ph/eNmbB9sGqTnj9V6ZcphgZzsVvrdGWAFf8QOfJLD
FKIa99rQZkHnUsQlB5kcZhAYPJtmMQce51PPKeuSMJ1FciJBaR58LnZ6hFmT
MGatYnTHt3mw1a98nHaJhA0YE+2aKhegLL/ioXIiCaMGVXoXJC5Cpjv52boC
ErZ5eFQ8cfk7sC0FL3o+kbBW+TUmrhZLoGGi1DouJmGxu5wm1ZQJ1DaDs9Iq
EhkzrLk5J4kgUB3ynBsSdMhY2bNc8VkRgarz77Otx03IWIJhTAzjNJGKfdzH
ZzuQMWuix8rtlURqhqYeI51OxkigUVGnJ0c9xeWFLYSSMQWtO7qUNDnqiRwH
M0oSGZPZzvnJL8hR7R2JAXOPyNgWlY0XH/uTqCGTOU0xr8hYkr1vVxHu4D/2
acXinoorOvUvbuOR6gvxuOvkfqLycStZtMiScPvLzmmuvkKiSnoa1z/AXdyW
Vh+Ge5vtj+fycVNZCrtPBpCoPuXp3Q24n31xfELDDWuqtnJxq6rf1/bHXWAx
FNCMW5i8ZcVt3Gc0dy+14b4ddaivAndZ331dAe4G+m9p+oEkKu0d5/wEbtPs
94omuCMj6p5O4n7UpRhlhTtrfeXYNO5gl4dXaLiHP2fdXMCtduSNbRpuxknL
LPnXZOxG9GJFNm5TtvbActxizHpPEe4tOrO7FHG/NWvW/YS78HBe5WrcG69v
uNeG+4xu+DIV3PHFnkoi3HYNdg5rcU+M5t2S4k5xVr+nhtt9p2RmCbdhcW+X
Ou7/AO6w/4U=
      "]], 
     LineBox[{{-0.00029788278179667445`, -4.980441107437789}, \
{-0.0002969806530481864, 1.5557544164966424`}}]}}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesOrigin->{0.5, 0},
  PlotRange->All,
  PlotRangeClipping->True,
  PlotRangePadding->{Automatic, Automatic}]], "Print",
 CellChangeTimes->{{3.397782723720703*^9, 3.397782729114873*^9}, 
   3.397782762638966*^9, 3.398276878277772*^9, 3.3982770433322783`*^9, 
   3.398277121837738*^9, 3.3982774791399317`*^9}]
}, Open  ]]
}, Open  ]]
}, Closed]],

Cell[CellGroupData[{

Cell["UndampedNewtonPlot", "Section",
 CellChangeTimes->{{3.398277577635768*^9, 3.3982775832946568`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"NP", "[", 
   RowBox[{"f_", ",", " ", "p0_", ",", 
    RowBox[{"\[Epsilon]_:", "0.001"}], ",", 
    RowBox[{"Nmax_:", "10"}]}], "]"}], ":=", "\[IndentingNewLine]", " ", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
     "\[Delta]", ",", " ", "\[CapitalDelta]", ",", "i", ",", "p", ",", " ", 
      "results", ",", " ", "p1", ",", "p2", ",", "  ", "u", ",", " ", "xmin", 
      ",", " ", "xmax", ",", " ", "xrange"}], "}"}], ",", 
    "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{"\[CapitalDelta]", "[", "x_", "]"}], ":=", 
      RowBox[{
       RowBox[{"f", "[", "x", "]"}], "/", 
       RowBox[{
        RowBox[{"f", "'"}], "[", "x", "]"}]}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"p", "=", "p0"}], ";", " ", 
     RowBox[{"\[Delta]", "=", 
      RowBox[{"\[CapitalDelta]", "[", "p", "]"}]}], ";", " ", 
     RowBox[{"i", "=", "0"}], ";", "\[IndentingNewLine]", 
     RowBox[{"results", "=", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"p0", ",", " ", "0"}], "}"}], ",", " ", 
        RowBox[{"{", 
         RowBox[{"p0", ",", " ", 
          RowBox[{"f", "[", "p0", "]"}]}], "}"}]}], "}"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"While", "[", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{"Abs", "[", "\[Delta]", "]"}], ">", "\[Epsilon]"}], " ", 
        "\[And]", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"i", "++"}], "<", "Nmax"}], ")"}]}], ",", 
       "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"p", "=", 
         RowBox[{"p", "-", "\[Delta]"}]}], ";", "\[IndentingNewLine]", 
        RowBox[{"\[Delta]", "=", 
         RowBox[{"\[CapitalDelta]", "[", "p", "]"}]}], ";", 
        "\[IndentingNewLine]", 
        RowBox[{"Print", "[", 
         RowBox[{
         "\"\<iteration:\>\"", ",", "i", ",", " ", "\"\< p=\>\"", ",", "p"}], 
         "]"}], ";", "\[IndentingNewLine]", 
        RowBox[{"AppendTo", "[", 
         RowBox[{"results", ",", " ", 
          RowBox[{"{", 
           RowBox[{"p", ",", " ", "0"}], "}"}]}], "]"}], ";", " ", 
        "\[IndentingNewLine]", 
        RowBox[{"AppendTo", "[", 
         RowBox[{"results", ",", " ", 
          RowBox[{"{", 
           RowBox[{"p", ",", " ", 
            RowBox[{"f", "[", "p", "]"}]}], "}"}]}], "]"}], ";"}]}], " ", 
      "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", 
     RowBox[{"p1", "=", 
      RowBox[{"Graphics", "[", 
       RowBox[{"{", 
        RowBox[{"Brown", ",", " ", 
         RowBox[{"Line", "[", "results", "]"}]}], "}"}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"xmin", "=", 
      RowBox[{"Min", "[", 
       RowBox[{"First", "/@", "results"}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"xmax", "=", 
      RowBox[{"Max", "[", 
       RowBox[{"First", "/@", "results"}], "]"}]}], ";", " ", 
     "\[IndentingNewLine]", 
     RowBox[{"xrange", "=", 
      RowBox[{"xmax", "-", "xmin"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"p2", "=", 
      RowBox[{"Plot", "[", 
       RowBox[{
        RowBox[{"f", "[", "u", "]"}], ",", " ", 
        RowBox[{"{", 
         RowBox[{"u", ",", " ", 
          RowBox[{"xmin", "-", 
           RowBox[{"0.25", "*", "xrange"}]}], ",", " ", 
          RowBox[{"xmax", "+", 
           RowBox[{"0.25", "*", "xrange"}]}]}], "}"}], ",", 
        RowBox[{"PlotStyle", "\[Rule]", "  ", "Blue"}]}], "]"}]}], ";", " ", 
     "\[IndentingNewLine]", 
     RowBox[{"Show", "[", 
      RowBox[{"p1", ",", "p2", ",", " ", 
       RowBox[{"Frame", "\[Rule]", " ", "True"}], ",", " ", 
       RowBox[{"Axes", "\[Rule]", " ", "True"}]}], "]"}]}]}], 
   "\[IndentingNewLine]", "]"}]}]], "Input",
 CellChangeTimes->{{3.398277598561446*^9, 3.3982776456003933`*^9}, {
  3.3982777376868467`*^9, 3.398277999221881*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"f", "[", "x_", "]"}], ":=", " ", 
   RowBox[{"x", "/", 
    RowBox[{"Sqrt", "[", 
     RowBox[{"Sqrt", "[", 
      RowBox[{"x", "^", "2"}], "]"}], "]"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"NP", "[", 
  RowBox[{"f", ",", "1.5"}], "]"}]}], "Input",
 CellChangeTimes->{{3.398277647243602*^9, 3.398277650956357*^9}, {
  3.39827801642795*^9, 3.398278017407426*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "1", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "1.5`"}]}],
  SequenceForm["iteration:", 1, " p=", -1.5],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.3982776514103727`*^9, {3.398277728091199*^9, 3.398277763869632*^9}, 
   3.398277800875503*^9, {3.398277878243208*^9, 3.3982778942325897`*^9}, {
   3.398277926760487*^9, 3.398278018194806*^9}, 3.398278112892748*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "2", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", "1.5`"}],
  SequenceForm["iteration:", 2, " p=", 1.5],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.3982776514103727`*^9, {3.398277728091199*^9, 3.398277763869632*^9}, 
   3.398277800875503*^9, {3.398277878243208*^9, 3.3982778942325897`*^9}, {
   3.398277926760487*^9, 3.398278018194806*^9}, 3.398278112895454*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "3", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "1.5`"}]}],
  SequenceForm["iteration:", 3, " p=", -1.5],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.3982776514103727`*^9, {3.398277728091199*^9, 3.398277763869632*^9}, 
   3.398277800875503*^9, {3.398277878243208*^9, 3.3982778942325897`*^9}, {
   3.398277926760487*^9, 3.398278018194806*^9}, 3.398278112898254*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "4", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", "1.5`"}],
  SequenceForm["iteration:", 4, " p=", 1.5],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.3982776514103727`*^9, {3.398277728091199*^9, 3.398277763869632*^9}, 
   3.398277800875503*^9, {3.398277878243208*^9, 3.3982778942325897`*^9}, {
   3.398277926760487*^9, 3.398278018194806*^9}, 3.398278112900935*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "5", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "1.5`"}]}],
  SequenceForm["iteration:", 5, " p=", -1.5],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.3982776514103727`*^9, {3.398277728091199*^9, 3.398277763869632*^9}, 
   3.398277800875503*^9, {3.398277878243208*^9, 3.3982778942325897`*^9}, {
   3.398277926760487*^9, 3.398278018194806*^9}, 3.3982781129036217`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "6", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", "1.5`"}],
  SequenceForm["iteration:", 6, " p=", 1.5],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.3982776514103727`*^9, {3.398277728091199*^9, 3.398277763869632*^9}, 
   3.398277800875503*^9, {3.398277878243208*^9, 3.3982778942325897`*^9}, {
   3.398277926760487*^9, 3.398278018194806*^9}, 3.398278112906241*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "7", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "1.5`"}]}],
  SequenceForm["iteration:", 7, " p=", -1.5],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.3982776514103727`*^9, {3.398277728091199*^9, 3.398277763869632*^9}, 
   3.398277800875503*^9, {3.398277878243208*^9, 3.3982778942325897`*^9}, {
   3.398277926760487*^9, 3.398278018194806*^9}, 3.3982781129087543`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "8", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", "1.5`"}],
  SequenceForm["iteration:", 8, " p=", 1.5],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.3982776514103727`*^9, {3.398277728091199*^9, 3.398277763869632*^9}, 
   3.398277800875503*^9, {3.398277878243208*^9, 3.3982778942325897`*^9}, {
   3.398277926760487*^9, 3.398278018194806*^9}, 3.398278112911231*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "9", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "1.5`"}]}],
  SequenceForm["iteration:", 9, " p=", -1.5],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.3982776514103727`*^9, {3.398277728091199*^9, 3.398277763869632*^9}, 
   3.398277800875503*^9, {3.398277878243208*^9, 3.3982778942325897`*^9}, {
   3.398277926760487*^9, 3.398278018194806*^9}, 3.3982781129137506`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "10", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", "1.5`"}],
  SequenceForm["iteration:", 10, " p=", 1.5],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.3982776514103727`*^9, {3.398277728091199*^9, 3.398277763869632*^9}, 
   3.398277800875503*^9, {3.398277878243208*^9, 3.3982778942325897`*^9}, {
   3.398277926760487*^9, 3.398278018194806*^9}, 3.39827811293997*^9}]
}, Open  ]],

Cell[BoxData[
 GraphicsBox[{
   {RGBColor[0.6, 0.4, 0.2], 
    LineBox[{{1.5, 0}, {1.5, 
      1.2247448713915892`}, {-1.5, 0}, {-1.5, -1.2247448713915892`}, {
      1.5, 0}, {1.5, 
      1.2247448713915892`}, {-1.5, 0}, {-1.5, -1.2247448713915892`}, {
      1.5, 0}, {1.5, 
      1.2247448713915892`}, {-1.5, 0}, {-1.5, -1.2247448713915892`}, {
      1.5, 0}, {1.5, 
      1.2247448713915892`}, {-1.5, 0}, {-1.5, -1.2247448713915892`}, {
      1.5, 0}, {1.5, 
      1.2247448713915892`}, {-1.5, 0}, {-1.5, -1.2247448713915892`}, {
      1.5, 0}, {1.5, 1.2247448713915892`}}]}, {{}, {}, 
    {RGBColor[0, 0, 1], LineBox[CompressedData["
1:eJwt1vk31IsfBnCTlGiRKxSVohJZKpUtb3GpKFrcLEmWSpQlS5akSJa6IhFp
kRgSamTfeg9CZLIMJsYaISXDGIMM38895/vDc57zOuf5A54t9m6nLi7h4+Mr
IvJf37uTxV5cJFEp5p+4i4tcNNGOtd7FI1GtbKRub1zg4l+cG9UWcyRqTsDd
b5rzXHx2yTg+Z4pEDVO/sc97hovvjo5qOvwgUQ3cJ0RHWFxsXSMfXNNCokac
kU2t6efi5mcpIpFkEtW6Pbf3UhUX8/ISlTYYk6jJSRspu0O5WJSe9hcvlI+q
Tr7QY7w4jbJFvOnypEV8qRNuFxA4jcqzAUMU9QVstb/jbs7j4Hj5AeV9ufMY
Q5mNanLhoJXAfb/Cw38wBp4HXxicwhtN6zd+Lp/FVX72TGvjKZzy9OlbpTWD
w923FJ5S2CiR6JJu0jKNU/R2ersUG1n6STUyxhy0F+731fKeRE93px+sbjYa
DQd92dc6gSbeS3axfCZRs8B3JEV1Ant57TW32SyMf7/TbcGHhWqi4fWbRMcx
pur9OpXwcVyn2lFow/cLfz/Wey89PobpC7WD8/MjGMLr/LD28C+MbBmWjFca
wk2dembfUkaxNzuQd3rLAEqklcxNlo/g3QeGP1Ye7UOX9y9uC7cO4a4xIeFm
x04sOq3qxIkcRMqmj5V/W7Vi3e87pAzhb5h8v+vaW/1WfPvxYP0r0jeMmZmS
kVRqxfsX1P885vajZ8u2Wz8X6ag1ruHuOtCP+0PDtGLIdIxlxiR8LenH0jGj
vF5WC47PrpR4eLkfa8qayf5hzShn67DMEPuwKP/DWjsnGn4Yum6Zl92DInzc
xEvqNIyg6Dtuiu3By0YqcleX05C2PX7Ow78HJXuT9vuQG/Dc4iZGp2EP+goG
nY3s/4xy38bVk3q6UeOsPrnYsh7nHaj91Su6sYTvk7qoUS2GbkvxST/GRFFj
vkqJ9bUoHfv96YAKE53j1I03jtSgapq/KukvJm5QyDgnH1qDKcyUzb1fO/HG
6XvBOpXV+Lr7fvMGh07UTjvW4Kz5ETdyxNJUrnVguXHz+SqFCpRNNAiXd2eg
c1pAhiCPilpKyiJdJxgoSZJnH2+kouVzv7MXdzPQoyAw9KsnFePfHg+ZnmhH
Uz0LK6Yf4h7xYsFF93aUG1vFUXlRhlM/oxWnL7bhrVcXwkv0ylDclP9Yhn4b
dpqXShkMl2JqodMHuS1tGFV5Wc9KtRTrih8YuDNbcS6+Kiqkshiz2/3DNY63
YqOen0LHUAHapXU5usrRUWGmsdzhfgFa/NjBXr3QgiHZ20/+VinAL2UmFnaM
FtSUbPXh98tHZofbcrGIFkwdU65WEs7DzmNhpddHmtEvYfB8sEoOyohfIwnH
NWHrcS32ylYKJvoM8N5dbkJl/pjQeF8KcrYP1zO1mrCWGWhleecdZousCHGr
aMRnoWqc7PQsPMzyO7nsyBc83JWkYM5KQ5pdxKFuic8on6vWz38wDW0CyRZy
2fUodK8unhJBxisd4eRi3Xr8coC9VEg2FQXypnN87OvQ7JFhb9k/yeji63jl
amQt7ndmxjm9eokrDqtetBCtReXgvoRdPUl40Vw8ixxbg3ndez17RZ7jzluZ
45nR1UiN7dxhcD0eWygzMktsq/CO3tK4F2seY5u1k2JqXSUaspSWzLyOxRu9
LL0rypVIMw7qesN8iMLX9UKPD1MxevbN0aXe0Zht/0+mpCYV00wyXEqlH6D4
XW2/HCHEXIETW7cEh+N8U6T/AVIZatqlzNPz7iKd8mEp16oEqWXT7aHDwdi1
ez3ZO6EIDSWNcjQ23MZi/qjNW3ML8MmwpkTHE38s9HDTuJich2kv+j/rLHqh
YIP8JU2z93gxU4R0tM4Vh6/9s0rO4R0GrrdxFot1QCmRfHbx8gxUXWOejuUn
0CO4sLfh3Ev0OqUumfZ+PwT9OrT74epH4KzumSsqYQ3cvKgL+/eSQVr5Kqtc
1wn481yETVyygGQbpT+r6QGN0e8klMRyYOr1yINHbT5gorUrVS0sF0YmDnUq
ud8Eqn7kntHcfHB9ThO9cDcIxE5apuRnF8LwIX/1nFMhkOL3726328VgO7Td
ZmFzGEgxQv5M7i2Fznv0O8ZjEaCYF6FRRyuDL/Rdjd/DomBMy0SupApBruyq
/xOFh7AwNemZFUMF29X6q3QfxsCW0YhwR+0K6Hw/vifSJg7Om5rlpN6sBAmB
muq91Y8hwEvwlqV4FZiZP7PoVEyAhaxuEZc3VUD7cyRwx2wixBip+WnUfATb
rc2DfUdfgrsVSfRlaw0EGmsrq/FegmrO3rHZE7Xw3CvdJ4ySDAHqCby6+lro
rLkppCyRAtZJfOvsCz+B2RUF1RuDZDj3qW13m089JLT40ATS00CcXrfZ9Vc9
dGlUO0c7pUN79RF38fOfwUHQlkz+/RooVn0n2ToN4EGO29A0kwl150hnRX/S
4JPb7eCYZTmwpIOrCUFN4Cy7o5nhmAPrUgpdtlKaYCWDtnljXQ6UaH81vdPT
BKY6UuXp99+Dzi2+iC8azdC+soBbJpIH5+OXmU+ONsP3jJ9XhzcUwgHtPHdZ
TTqEnosp3RVQCGoJAtIe9nSQX6sh5NFdCGRvuxSXe3S46huaPp9UBMulaqU8
O+gwZbjl29ptJcB4KXxt1LMVlgyeMddWKYfrDsZBv562gcia/pefkisgzfVv
224qA/61qAhTbKwAt3dh8ta9DFj+Ktn1wXwFRHu9UsrjMWBBzV7bzLwSxHRA
bVLjK4xZDTB6V1bB25ZjZ/MpX6E+7ftqrs9HoH8I3lYb3wEhOqMB20xqYcP6
401CpkyYdWFbBM81gJZSCe2NZw8sGDaI5SrSoGBPnYVvZA8skSE3DVjTINPt
jqJqeg+sbDlzxABpIPJ0ddD5zh6Q2V9yQDDkC4Qdfp1ZCb0gd0tXuHZZE0RJ
t+104O+DY4eXjx6LbobpgM+6AVV98JQRm26R0grDhU8qJxS/AeOxV0puVid4
2O2NzH09CAXb/I4I7+iDf7Zm2uwsG4K4GerqRLUBMPc3iL1UPwIyNQzartVD
MKV9yIsXPQqz304ct+sYgTYrZs3gnl/gOHNK2rP/JySsKdr5sGcMggTlBsPz
f4Nv1I9lay3H4Ykg8g4OsCCAI2Mv7s4Cw3v0u8UnJ2HQJnaqedMEPN27/rdA
MRu+Pmi0eF41AZqWsgtUHQ6cDh0p7HabBJX5WKHAN9OQxDljJifMBjepzx8y
1GYgcF+f1pK3bDhjv+d0WuEsOCYuLWYdmYJH20mDqTJ/wMvjO+dy2xRIOcZK
i76eh998vrXpThyYkYs6GbRuATblRGSFsjhwU6j0UsWNRZBcwzTKujYNOwra
1qrf5tP9/vtPZtzcNJSQg1yj95F0D+2vGw/14YIHfc1mmTmSrvOciuS6eS7s
MbR/cYDwCgbroxxhdlHeRlPCB73Tr6kR9nphKR1IOH+aU3+asLdTyvouwqtE
JG/GEPZd3C+W8Ieku/Ste+8aHhcCFW1WiPBIuoGyXs+EFrig84ISvoNwouXY
4Q2EF0X4BYFww2lr9k7Ct6fTl7kSromVMjpKOKhigr+BsI7xZm4Y4RDzuwth
CyTdxqdapgKLXPi7nnEzifAH6bVzYoQFDirwCggvd+hKlSMcurXxz3fCv8Bs
Tp9w2JjkrP4iSdfgYF1qMOEjtld8zxKW/XnMNIawIL2c60F41qZmNplweJH9
dDJhZkKiScV/e8V872LCIuc5M02EV7xYzmki/OCnQUof4XoRK68Rwo46/x5n
Ef7/39b9/9+G/wF8mCX5
      "]]}}},
  Axes->True,
  Frame->True]], "Output",
 CellChangeTimes->{{3.398277886367362*^9, 3.3982778942487373`*^9}, {
   3.398277927296095*^9, 3.398278018208499*^9}, 3.398278112997156*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"f", "[", "x_", "]"}], ":=", " ", 
   RowBox[{"x", "/", 
    RowBox[{"(", 
     RowBox[{"1", "+", 
      RowBox[{"x", "^", "2"}]}], ")"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"NP", "[", 
  RowBox[{"f", ",", 
   RowBox[{"Sqrt", "[", 
    RowBox[{"1", "/", "3"}], "]"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.3982780337621202`*^9, 3.398278085240901*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "1", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", 
    FractionBox["1", 
     SqrtBox["3"]]}]}],
  SequenceForm["iteration:", 1, " p=", -3^Rational[-1, 2]],
  Editable->False]], "Print",
 CellChangeTimes->{{3.3982780486893578`*^9, 3.398278085676861*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "2", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   FractionBox["1", 
    SqrtBox["3"]]}],
  SequenceForm["iteration:", 2, " p=", 3^Rational[-1, 2]],
  Editable->False]], "Print",
 CellChangeTimes->{{3.3982780486893578`*^9, 3.398278085679923*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "3", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", 
    FractionBox["1", 
     SqrtBox["3"]]}]}],
  SequenceForm["iteration:", 3, " p=", -3^Rational[-1, 2]],
  Editable->False]], "Print",
 CellChangeTimes->{{3.3982780486893578`*^9, 3.3982780856834393`*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "4", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   FractionBox["1", 
    SqrtBox["3"]]}],
  SequenceForm["iteration:", 4, " p=", 3^Rational[-1, 2]],
  Editable->False]], "Print",
 CellChangeTimes->{{3.3982780486893578`*^9, 3.398278085684988*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "5", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", 
    FractionBox["1", 
     SqrtBox["3"]]}]}],
  SequenceForm["iteration:", 5, " p=", -3^Rational[-1, 2]],
  Editable->False]], "Print",
 CellChangeTimes->{{3.3982780486893578`*^9, 3.398278085686578*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "6", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   FractionBox["1", 
    SqrtBox["3"]]}],
  SequenceForm["iteration:", 6, " p=", 3^Rational[-1, 2]],
  Editable->False]], "Print",
 CellChangeTimes->{{3.3982780486893578`*^9, 3.3982780856879797`*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "7", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", 
    FractionBox["1", 
     SqrtBox["3"]]}]}],
  SequenceForm["iteration:", 7, " p=", -3^Rational[-1, 2]],
  Editable->False]], "Print",
 CellChangeTimes->{{3.3982780486893578`*^9, 3.398278085689472*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "8", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   FractionBox["1", 
    SqrtBox["3"]]}],
  SequenceForm["iteration:", 8, " p=", 3^Rational[-1, 2]],
  Editable->False]], "Print",
 CellChangeTimes->{{3.3982780486893578`*^9, 3.3982780856908627`*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "9", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", 
    FractionBox["1", 
     SqrtBox["3"]]}]}],
  SequenceForm["iteration:", 9, " p=", -3^Rational[-1, 2]],
  Editable->False]], "Print",
 CellChangeTimes->{{3.3982780486893578`*^9, 3.3982780856922073`*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "10", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   FractionBox["1", 
    SqrtBox["3"]]}],
  SequenceForm["iteration:", 10, " p=", 3^Rational[-1, 2]],
  Editable->False]], "Print",
 CellChangeTimes->{{3.3982780486893578`*^9, 3.3982780856934566`*^9}}]
}, Open  ]],

Cell[BoxData[
 GraphicsBox[{
   {RGBColor[0.6, 0.4, 0.2], 
    LineBox[NCache[{{3^Rational[-1, 2], 0}, {
       3^Rational[-1, 2], Rational[1, 4] 
        3^Rational[1, 2]}, {-3^Rational[-1, 2], 0}, {-3^Rational[-1, 2], 
        Rational[-1, 4] 3^Rational[1, 2]}, {3^Rational[-1, 2], 0}, {
       3^Rational[-1, 2], Rational[1, 4] 
        3^Rational[1, 2]}, {-3^Rational[-1, 2], 0}, {-3^Rational[-1, 2], 
        Rational[-1, 4] 3^Rational[1, 2]}, {3^Rational[-1, 2], 0}, {
       3^Rational[-1, 2], Rational[1, 4] 
        3^Rational[1, 2]}, {-3^Rational[-1, 2], 0}, {-3^Rational[-1, 2], 
        Rational[-1, 4] 3^Rational[1, 2]}, {3^Rational[-1, 2], 0}, {
       3^Rational[-1, 2], Rational[1, 4] 
        3^Rational[1, 2]}, {-3^Rational[-1, 2], 0}, {-3^Rational[-1, 2], 
        Rational[-1, 4] 3^Rational[1, 2]}, {3^Rational[-1, 2], 0}, {
       3^Rational[-1, 2], Rational[1, 4] 
        3^Rational[1, 2]}, {-3^Rational[-1, 2], 0}, {-3^Rational[-1, 2], 
        Rational[-1, 4] 3^Rational[1, 2]}, {3^Rational[-1, 2], 0}, {
       3^Rational[-1, 2], Rational[1, 4] 3^Rational[1, 2]}}, {{
       0.5773502691896258, 0}, {0.5773502691896258, 
       0.4330127018922193}, {-0.5773502691896258, 
        0}, {-0.5773502691896258, -0.4330127018922193}, {
       0.5773502691896258, 0}, {0.5773502691896258, 
       0.4330127018922193}, {-0.5773502691896258, 
        0}, {-0.5773502691896258, -0.4330127018922193}, {
       0.5773502691896258, 0}, {0.5773502691896258, 
       0.4330127018922193}, {-0.5773502691896258, 
        0}, {-0.5773502691896258, -0.4330127018922193}, {
       0.5773502691896258, 0}, {0.5773502691896258, 
       0.4330127018922193}, {-0.5773502691896258, 
        0}, {-0.5773502691896258, -0.4330127018922193}, {
       0.5773502691896258, 0}, {0.5773502691896258, 
       0.4330127018922193}, {-0.5773502691896258, 
        0}, {-0.5773502691896258, -0.4330127018922193}, {
       0.5773502691896258, 0}, {0.5773502691896258, 
       0.4330127018922193}}]]}, {{}, {}, 
    {RGBColor[0, 0, 1], LineBox[CompressedData["
1:eJwt1mk41XkbB3BZSo0We9SgjCXkjFBp6DZKydJxLI0xmig5WRqy0zPykEIU
6tFiKWt2sk7xuD2SZMm+ZPn/Qjjnf0o005Bpmf9zXfPivu7r8+q+v+++2076
2p0WFhISimLm/9tp68+D4bUCHLolrL1STFDkuXKHSqUA61MCRw8wLrsw2dRa
KsCDnw6WJhQRFJ12L5HKEyCId9RvKSRYUegdXXRdgPsN5gNY+QTX7AnTG/UT
4MLRYImd9wjWcm4kGmkLsGX6w6OI6wTdRI5FE3UB4lMxKEkhKFEjHxazXYDc
qrzxkWSCJzenufcqCHAN3/0UK4ngBpL1nae4ABtXDFJ7Ewie9inn35qhcX7y
uhUdQ1Dmcrv58l0ap07FtHqFEOS/T7hkmkYj5w+NIo9ggo3u7KexqTQWSRkE
uAUR5JoNWCgk0iic86bKPoDgo08TVvvCabxikvNKy5f5L3CR8y8HGqu+eV6Q
48Hkd9v8s9BaGpXmimfrHQhaGHsEi/vz8RzKHffbSXD2uGHTvbN8nCngpezV
IRhzQXSdkScfY7pS6r9oEWxuzs7wdOVjRobsq3hNgt8dIS3tR/k4auYSe0OV
oO4xJ+lEbT461on3BSowefwsKyRneKiotnotV5QgydGlFZ14mCqvP/t5iMKG
86vtn9nzkNobnBM7SOEde6o+hM3DFmeRI5IDFDqIXE0cOMRDI7rMSrmXwvaT
b769asjDvEN66rodFFZvKwlZJc3DkCqxOTmkMP7uDjFe1xzmbzWxlcqn0DDt
G5WaA3OoPNYSlulHoejEEbkp+VnUWGy+bfN6Ag1709iSQjP45aFe89mfJpBs
/PpIf880mmuv6Y3NG0cLqkP06/tTKCnRcwlGxlDmxH0Rz+RJ5CbEBud+GUUn
rcrmw+4vsfe1mJ/a8gv8zO6P46RSeE+jLOWH2RHs8MqQUvUdR0WD4DhdHMY4
xdYtigqjuGrl8Ngl1yGUYFkUywqGkbNl283yjwMYOaJZv/vqIMbnfXkp4tKP
AVnubK5xP1KzIhxubi9qd3/wTJfqxQIbJ7Ool91YqDq2dj6gC3+v7BfvWO7A
S/5B1iVSzzDhv5lev6q14SFB57Ua+xbM3rspvbr/MYodqqmrd0QMnTOetSON
eN1it6dSbi2WZliYbA+pxU6TU3t1fylCmcdJebtcitBom2X25YJkDNy16Q7r
fjIafx856zCaBB+XW/t2vEgCdlFFtY9aMUxq16UvShYDvbrY7cnRWrgcGPlX
vWot8N2juftbEULLz7AHLRFCspdKG2JbIMD+fO2IaAu8bQvfwqHaQHRy4kTT
n08h5olp7J3kLrDbv0tIRaoTLD7uCDy70gPdzY902A+7gR08ovxZpR+kYGG9
x7VecP9J08ozeBCyMw8oL1v2Q/WY0kJ3zjAoWfmUm3YNgPTTWhlJoVHY89u7
MA3OEJTavK+0MRiHo5WeRD5zGNqqRzrDj1NwA0ZtVZ+NQNhjXf3/Ob8E53Us
Jdu3L0A7vKWh2X8SnE6rfmUwOwqnZ7Uu6WVPgVzriNi9gTEgzacSN3dOQ1um
18L57HFYeu3L+TT5Cm4V9WUlW05Aaqe1I3vjLOx8wMnfyZuA/g8fZe9qzsHJ
gg6p26cpsBBdMpzXmQP9WIfVx7kUNG5YdDTRm4Pgst1dyp4UFKnO/GfMaA4O
J9oPZflQEGXTJbvZag4ch8XJzQAKdmVlyCafnYO7UdrS1v+mIMViv2x05Rxk
cC9aD6dRYHvzggzXmAd5mq7fXumiIP1VKF/6ex5sDDetUe+mgKfn39hkzoO4
3gcLTT0URHa6n1Fk88Bg6k3CYj8FFassG7rceHDF45jHgVEKJH1k3A0u80BX
TftD4SwFfaaFVav6eGBmdOVCzGcKHOk+u3QuH7TcdFYStAiMXESlD958+Oo7
852PtQk4K5fQjn585p6507IOAVeHmKiNYXxQOHj05gkWAa/GPZVR8XyQbi9q
3mZAICIlXdKzlA8XO/nyQSYE8vd59Bj+zgenAhcxZzYBzUG79JQlPvBcrgz7
2hIo9oUzb//iQ9u7TSUXOQQqcjcLFYrSsH67z/FiewKPNnSwtsrRkLTgS7/+
gcDzKdY1YSMa4j3+fG/qSmAp/oNNdwQN6jUK0iw/ArZnZp7HRdMQ7fIkdN85
AgXmveyDsTRUasQ2HfRn8ggVch4m05Co/27JMZBAQ5CTY04uDXq8UknvUAJR
rnUuIe00yNDGVq6RBCR2B3krywlga7HS/bBrBFQmF2J9SwXwtqzIOqKIwPrU
7bXXGbN1ww5cZLxi6TBdx/iIuMS+eMYD1XUmQmUCePw8Qj2V8eXYiHdJjPt/
u/WxlPE8S8KlslwAmbczMyYYN0RosN4/EEA9e6nSuJhAof6PLgqVAuAG2GWb
MU7lxceZMPZa459swdjPbn46hjEr+vMvDoxV1WtuyVUJoGXcWsOHcXyXmfDe
agHoVGBCGuOQqECWC+NxrZOhWYzd9+S7RDJuUXl26j7j/Vlr69oYv458ZVTF
eDGwx/vHGgHkqXGm2xlTO4Rv/8p4XaJ8Vw/jdkq/NYuxdFB27RDjXIubKjTj
rmdicVOMkz+12WyoFcCLPdQ5HuOIypXwXYxhx3nnecbeXJ2CY4zLM8fN/mD8
Tx+Af/oA/A0ff1iW
      "]]}}},
  Axes->True,
  Frame->True]], "Output",
 CellChangeTimes->{{3.3982780497760553`*^9, 3.3982780856956387`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"f", "[", "x_", "]"}], ":=", " ", 
   RowBox[{
    RowBox[{"x", "^", "2"}], "+", "5"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"Solve", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"2", "x"}], "\[Equal]", 
    RowBox[{
     RowBox[{"f", "[", "x", "]"}], "/", 
     RowBox[{
      RowBox[{"f", "'"}], "[", "x", "]"}]}]}], ",", " ", "x"}], 
  "]"}]}], "Input",
 CellChangeTimes->{{3.398278911415411*^9, 3.398278933074584*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"x", "\[Rule]", 
     RowBox[{"-", 
      SqrtBox[
       FractionBox["5", "3"]]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"x", "\[Rule]", 
     SqrtBox[
      FractionBox["5", "3"]]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.39827893503913*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{
   RowBox[{"NP", "[", 
    RowBox[{"f", ",", " ", 
     RowBox[{"Sqrt", "[", 
      RowBox[{"5", "/", "3"}], "]"}]}], "]"}], ",", " ", 
   RowBox[{"AspectRatio", "\[Rule]", " ", 
    RowBox[{"1", "/", "2"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.398278939578898*^9, 3.3982789840697947`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "1", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", 
    SqrtBox[
     FractionBox["5", "3"]]}]}],
  SequenceForm["iteration:", 1, " p=", -Rational[5, 3]^Rational[1, 2]],
  Editable->False]], "Print",
 CellChangeTimes->{{3.398278969332817*^9, 3.3982789845353203`*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "2", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   SqrtBox[
    FractionBox["5", "3"]]}],
  SequenceForm["iteration:", 2, " p=", Rational[5, 3]^Rational[1, 2]],
  Editable->False]], "Print",
 CellChangeTimes->{{3.398278969332817*^9, 3.3982789845374107`*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "3", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", 
    SqrtBox[
     FractionBox["5", "3"]]}]}],
  SequenceForm["iteration:", 3, " p=", -Rational[5, 3]^Rational[1, 2]],
  Editable->False]], "Print",
 CellChangeTimes->{{3.398278969332817*^9, 3.398278984539987*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "4", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   SqrtBox[
    FractionBox["5", "3"]]}],
  SequenceForm["iteration:", 4, " p=", Rational[5, 3]^Rational[1, 2]],
  Editable->False]], "Print",
 CellChangeTimes->{{3.398278969332817*^9, 3.398278984541194*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "5", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", 
    SqrtBox[
     FractionBox["5", "3"]]}]}],
  SequenceForm["iteration:", 5, " p=", -Rational[5, 3]^Rational[1, 2]],
  Editable->False]], "Print",
 CellChangeTimes->{{3.398278969332817*^9, 3.3982789845424137`*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "6", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   SqrtBox[
    FractionBox["5", "3"]]}],
  SequenceForm["iteration:", 6, " p=", Rational[5, 3]^Rational[1, 2]],
  Editable->False]], "Print",
 CellChangeTimes->{{3.398278969332817*^9, 3.3982789845436172`*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "7", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", 
    SqrtBox[
     FractionBox["5", "3"]]}]}],
  SequenceForm["iteration:", 7, " p=", -Rational[5, 3]^Rational[1, 2]],
  Editable->False]], "Print",
 CellChangeTimes->{{3.398278969332817*^9, 3.398278984544882*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "8", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   SqrtBox[
    FractionBox["5", "3"]]}],
  SequenceForm["iteration:", 8, " p=", Rational[5, 3]^Rational[1, 2]],
  Editable->False]], "Print",
 CellChangeTimes->{{3.398278969332817*^9, 3.398278984546088*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "9", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", 
    SqrtBox[
     FractionBox["5", "3"]]}]}],
  SequenceForm["iteration:", 9, " p=", -Rational[5, 3]^Rational[1, 2]],
  Editable->False]], "Print",
 CellChangeTimes->{{3.398278969332817*^9, 3.398278984547414*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"iteration:\"\>", "\[InvisibleSpace]", "10", 
   "\[InvisibleSpace]", "\<\" p=\"\>", "\[InvisibleSpace]", 
   SqrtBox[
    FractionBox["5", "3"]]}],
  SequenceForm["iteration:", 10, " p=", Rational[5, 3]^Rational[1, 2]],
  Editable->False]], "Print",
 CellChangeTimes->{{3.398278969332817*^9, 3.3982789845486803`*^9}}]
}, Open  ]],

Cell[BoxData[
 GraphicsBox[{
   {RGBColor[0.6, 0.4, 0.2], 
    LineBox[NCache[{{Rational[5, 3]^Rational[1, 2], 0}, {
       Rational[5, 3]^Rational[1, 2], Rational[20, 3]}, {-
        Rational[5, 3]^Rational[1, 2], 0}, {-Rational[5, 3]^Rational[1, 2], 
        Rational[20, 3]}, {Rational[5, 3]^Rational[1, 2], 0}, {
       Rational[5, 3]^Rational[1, 2], Rational[20, 3]}, {-
        Rational[5, 3]^Rational[1, 2], 0}, {-Rational[5, 3]^Rational[1, 2], 
        Rational[20, 3]}, {Rational[5, 3]^Rational[1, 2], 0}, {
       Rational[5, 3]^Rational[1, 2], Rational[20, 3]}, {-
        Rational[5, 3]^Rational[1, 2], 0}, {-Rational[5, 3]^Rational[1, 2], 
        Rational[20, 3]}, {Rational[5, 3]^Rational[1, 2], 0}, {
       Rational[5, 3]^Rational[1, 2], Rational[20, 3]}, {-
        Rational[5, 3]^Rational[1, 2], 0}, {-Rational[5, 3]^Rational[1, 2], 
        Rational[20, 3]}, {Rational[5, 3]^Rational[1, 2], 0}, {
       Rational[5, 3]^Rational[1, 2], Rational[20, 3]}, {-
        Rational[5, 3]^Rational[1, 2], 0}, {-Rational[5, 3]^Rational[1, 2], 
        Rational[20, 3]}, {Rational[5, 3]^Rational[1, 2], 0}, {
       Rational[5, 3]^Rational[1, 2], Rational[20, 3]}}, {{
       1.2909944487358056`, 0}, {1.2909944487358056`, 
       6.666666666666667}, {-1.2909944487358056`, 0}, {-1.2909944487358056`, 
       6.666666666666667}, {1.2909944487358056`, 0}, {1.2909944487358056`, 
       6.666666666666667}, {-1.2909944487358056`, 0}, {-1.2909944487358056`, 
       6.666666666666667}, {1.2909944487358056`, 0}, {1.2909944487358056`, 
       6.666666666666667}, {-1.2909944487358056`, 0}, {-1.2909944487358056`, 
       6.666666666666667}, {1.2909944487358056`, 0}, {1.2909944487358056`, 
       6.666666666666667}, {-1.2909944487358056`, 0}, {-1.2909944487358056`, 
       6.666666666666667}, {1.2909944487358056`, 0}, {1.2909944487358056`, 
       6.666666666666667}, {-1.2909944487358056`, 0}, {-1.2909944487358056`, 
       6.666666666666667}, {1.2909944487358056`, 0}, {1.2909944487358056`, 
       6.666666666666667}}]]}, {{}, {}, 
    {RGBColor[0, 0, 1], LineBox[CompressedData["
1:eJw12Hc8V9//AHBbZYVrzzeVBkplJDknM0r2SMgelVHJyB4ZlZns7IyPyEqy
jh3ZkRFFUmaEUETf0+Px+/11H8/Hfdx7H+ec+xrnkCyddWwoyMjIWMnJyP5d
2Us4Uj9t76Isx+r1v34kyH9MjpV8cxfV3j4q8syHBA/mmj8QXt1FiQ7qr9S9
SPB0aoGH/ewuihNNXot3I0HtsLP6q4O7KJXWc/yUIwk+Mr/GSFO8ixaYTzv7
GZMgOXOev5j5Lgp3T080lSTBxVtSNt4tO6g15OK25KIgvFrpyDB3+w/ils/n
L9IRhFTzJbEV0tvo9wT/yInHArB0tni9mGILxTXfUBKP5YebDq3SHuW/EMeo
mHNMGB90/zBjkHFvE/0s1D3F+JQXRncL6SeIbKCF+6ZhWoU8sPghp0Tdh58o
m6gbzezkhjRnuOfP56wh/rOsL0VpueGm/enAJ7qr6IKbtcItwAWFwow3VblW
UJLiqrRHDCfMmXaL03m7jO5UW3xLX+WAUvUOv8rEvqMC+5IKS10OKOfmyt5f
N4/ynqxL8XSyQ9dTJauPuGaR3T6DrS4ldvgyxOJBbsBXdLBnczukhw1erzo6
7sP+BbkxReYb3WSD2ssHFlfdJpHKgEkWGxMbzJXWeCXiOIlYZxsDBhjY4B+/
u4EmVpNohmpHMJYe39/fxtmmOYnyDdbTmPfh+xJ2ykmHJ1Go8VIiFzUb1JzT
OedXMYEavpDf19gm4J+1q1mcXp8Q15/ic8qzBDTc63Tz4p5xZK9EknBtIaCo
UMkP5rUxNK6/fuBiMwHJzq66jnwcQ4wUtFtCTQTMv+nmY1M+hkqPflF9jwi4
1esb6Wc2hhZrfW8r1RAwLT6qpLziA5L9qkCpUUbAaeHSnzwWo+ik8l8W2kwC
VsmteUxdHEUPikPEt9IJGKEvuZMvNYqyyi4eW0ojoFToa2op+lHE+DCyfSSV
gKHzjexar0bQuPdmTFUiAY+WvpMJZhhB2tYG5lnRBHSR/+m9+HoIJShcWOAP
IOA6Jcsy27Mh1NCREHnSn4D3Oo5bgOghZKRJRafqR8AgvRvKsbZDaFD/dNRt
HwLGX59ikGEdQmRaP1oGPAlYE9+XFnTzPTrye/hr/W0C0iw/b+AWGES1/LY7
0jYEfFDx9qTSvkHUcke/xdaagEz3ZnMc1weQ3BsL63grAnJSHwhv6BxAHg79
x39ZEPAYT4q2rccAOphC/7PVjIBaquFTpf3v0OXyVYkHRgRMSbehUgvuR5OV
PkEvLxGQOKBnPeLSj365j6TRYT/KV2ixM+1HTyZF3C0vEtC7VCA4RKofWebn
mLGoE9CseZSyZbYPnWdh3+ulSkDBGQ1KoNGHODVHOp0VCKhmepPqumgPEnX8
sGdGhoCytmotYdQ9iEzgbZMe9lHnQ8F5n7qRkUfxmWZpAtIHTFJ+je5GqSR9
sywpAvbk6FGar3chh8vcQQ6nCaizKEdhUN+J+JiDZXlOENDEi4FMQasDzcm+
N6c/TMCMuO4U7kMdqPH2V8dQEbz+RRHSa9vt6Phe4gol9o1JBuecvHbUqd3V
sXMQr48y4yeav28QuU4BsSNMwEQmptrOojbkyMimdVCQgOOHew1zgtpQynE6
UCyAx6cQteZ9pQ2ZvhCgl8HOc2U6Jk7dhrLfZsPL/ASsHGVKjjZpRWIunKz3
eQk4kL3fQ5+uBUH3clNeLgKy1/Wxin9uRk0BUeLlnAQ0Hop+QfOqGY3Oh39U
x57awzzzyrIZtY+MvvfmIOCKI7MBd3UTyhX+dnyBDa+3DMvpCftGZH63JHyF
hYBPEkUm+x/Vod499IM5DARcaJvv4ZavQ9/8o4JNsM+vF9VZLdci3eDyDQL7
u/bplHXtWiQu9lEgjJ6AynTnDbg4atAszav77nT4//Ux7jLPrEKOH4LT7+0l
4MXnfDUFOlVokpb5ujx25ofJglXKKpTfcG6DAltD2i402PYVkhn00I3cQ8Bn
P+4o5B+rRGZBd0nFtATUt4qoWn5ZjmrIw+lYaAhYGKOVJ2NXjsKm+X58piYg
eQNrfABnOSrU0Y0vwy7iTb7D6lWGVvvSRQywqYdyxaVhKQrWR6a5VAS8SnWd
z3+lBO3GGml6Y5eeFKPvyCpBvVdjKXSxTaPK54ypS5BL1EwRBfb1GH+1nbgi
VOXLyXmDEs//Eo1I1OZzNMleya6K3aQeQUUyfo7MZH8OCWNzUyUjJYFCdFVY
J3WCgoBv3colHxXkIyOfOFlzbLeBM6z89PlI3EVKRQFb6ETDjxdOeejKkgfr
Qex7c12FA6dyUcmjidVFcgIeMflG4qnPRtbap1rvYxvxqzu1CGejqUIRPSfs
kMmiasfwLKT++fRTQ+wp67u6DXqZ6MN4kqko9n6R0XSH6gxEnr0yyIEtPye3
yCKYgViULGmosCXH9uzrHHqKfrJHJnwiI2Dstb4svqxU1HpBbbMLe+lLwlkX
xxTkp8/HXIudu3jIkY0qCVX8ZbqTik1xa4navjcBiT0Vr4nAvrb+Mq06OR7F
sNnW+WFz7ir1XZOIQ5wJAj9ssF0D6R3K/sQi2YxMQWPsfppBcur2GKQaLEen
iS32MCXZ8HE0Aoo/Xihhswq02erVRSLBj7WMZ7G7RPQTf+s+QovaWSInsYNP
THekzYejrMyMrSPY587c2VYMCEWKXa9ChbA3zlOIzXHcRz39Mz082C/UY8wi
iwORYtmJPjZse13B6FPK/rhfiH24H/t8lVlQ6SsvtEjJQE6PXdbh/W172g0J
3M04sQf7wFiymgrLbSTfeImbGvvo6K0rJL/rSGWWuZ4Cm+FPkUDlqDk6M7/M
Q47d4JWfcKz/EmJ/8+0U2b/vrXSnUYSfAp98N6j/OU7qLA+HmBGIYhaI/Wc7
D4s3GzU24Jif+ci/54P1D5JbWjiD0jfVI5TYWkmuT22DXIHQgthjGuyJUANZ
XgNP4DP3mnYftpObzHD/YV/Q1GQmxYjte5KhutknAKy48QqwYrce8Q3jnwkC
dDRrjZzY9KRlA0+tEMDkMsUvgK3HaX5w8HUY2KmYlzyEncLUvyYu/BAMDdLR
imNP0Sg0hT+KAEm9yrFS/8a/WxY9vR4FVJ4lDgNsE7e37BrWseDa57BoPWyy
VQrfrcTHoFP1EJU5drbj2a953XFAKHxc4ib2vPXzcgrpBBCR+rgmBDtycpq7
5EYiyHaK43iCfdKEL9A0Iwlk8RSL5WB76kZqVe1NBeY0omEt2Ly9ba9s5J8C
fqXorvf/5lv9Lz/rnTTwxoShcwabilS1IeWUAYb5dTYY8P9+kfNIa8xwFlj7
Sk1jhd1HnNlL7pcNMm2OcHljGzCrabgcygEZKdq08djm+xzeX777DMQdLXHt
wnb7k/+VjjUf6A0+e6yI43X7V9URr+p8EEOz7GuF7b/e7jhvUQAEluWUg7Ef
Ls2ut5f8B05++WbZjp0xeZgm5HIREDe8M2WM88WhjzLq6+tFYOvs1LMA7MLR
C5HWT4sB57C+UgF25Tt7dsXFF0DtxEWWbWxbSz1WMuVSEOjktJqP85OXXpxw
3VQ50OJOaBbD+bFeJ7XVVaoCyFprS1tik2nn2Ik+qABu1vtDErBDNCoKUyRe
gunQ1AoqnF+jVQZP3fOvBAnG3/XnsHNlCCVp/tfgyvFbBWM4P89K8X5buvUa
WHxy6+bbh/sByQNhua2vQa1iQJc5donE6S42x2pgJ5FzdQ675qie3s+aGvA4
UaieCteHPt4469Ir9eBCPHmcNSMBlQSuaPbn1QMJW2X6F9hVJH7ZlfV64LKk
pLiFnXEon0kiFoEVEffvsUwElOE0MzWyaQBb8vknu/cTcK/Q7T0jTxuB6I6c
vw0rfl/24M5MSyNw05GbrsH2PSC9trnQCNjMN2lYCQKuHvrzkUO2CRjkZIQ1
Y384FlJuONQEJvuQtxg7Af+TTDYbZmwBbNUnDh/E9XW68o/ujGQLqI9bcA/F
5pe5prZpgi1iGTKPHSt74DRHYQtAboETZdx4/kHxXkPVVvC2KirwEq7Xl9Sa
Kob82gB59/jDbFzfv5vM7xv60Q5WJc8XTh4i4GC3wwFRng6gDTg4bXA/US0/
dy5QuQMwrf2ynscOE5x1EU/uALlXJtI2cT9yYPrrUIjiW6AiLPOY9xiOrxuf
s6TjO4FiPpdv9HGcDz1HZJPO9gA5z1QTwTMELJ831Fu26QEuPTNvG7CTrw47
Kkf3AIoeeS5LWZx/zg1lrEz3AKUuJoO8szg+yAdp1CN7gVTrOS05eTzeFjsH
H5c+MJbfoRKtiOvvRStJ8eZ+ICidKTd5GedzQ1W9jeF+MEFZdSZJE/c3lsfu
1C/2g9jbiFFXi4CW99ZKLrO/A8OSGbBDG+ejgiBRZ4d34JNuQTPSw+OjzREu
YRoAw99HVd4bE7C2eZr5pMkg8E+qOJ9lh+uhnO2S5PoQyO3V3NsShNcvzFMg
cc8w+Eb/eSYrmICN7x9pbfEMA6vmy4GB93F8OpeX1SsMA5ZUNK8Yivv3LDJ3
1ehhoP5nI7z3AQF59qb8NTo2AkgeteT7YnF/OtS739tiFNiZZg3R4X6exuXM
qabuMaD7aZRNsZGAz09fdXszNQaGZ5br7PB+Qee39+uuzTGw5f9TIQLvJ9IC
G8AwaRywW3Z8+dCK+/0nqpcW746D8olMf/+3OD6r9W3Y+T8Cz3YGkfUBArZR
3U647vQJyGj1HwqeIWDzpfc+4aGTwFuGdlWGmQ0yUW6IuG5PAbaCJB1dezao
I5zJTm71FaSuRLhbd7BBK3ke1R+7M+B58F3GbsAOybx86khO80Dmjf/ZyGp2
eKOs+O5L5u8AvFc0oNLggBZ07Xcehi2Dc04BycZTHFBkRnIs4PcPkHWEeBbh
zwmpVZz5li+sgmcSjGlXznDBgpepgTsBa0C90fmh1CoXrMgq3XPu/U+Qx7mQ
Tt/ODWta6Eo1hDZASVBWV14eDwzbvOL912gTSNxftboTzAtr86Lspwp+gdQZ
pUbRAD5Yaf/a3Xb6NzBcOX843ZMfFpodLqjh3AbXT6R/nfEQgLTvfLfP2/0B
RsD4qL+GIDxi0OHikL4D/lO31aQYE4Rlhts0TSu7oFaGMYYH78+JDk3PeezS
ai+B49h3ZXMWWFZ3gUkXRbECtgyfRp8V9kCaaacDdv1UWhL12i5QJVOhrcLu
dFQQVfu5C5Rffg7S8ybB6aBw7b6NXZA0oxQd4UOCKusfm39hXzCfJmVi59me
lCJt7oLn9k/LK7Cvq41x38Z+cthkeAz7B6PYNOuvXbBkTS941JcEd5L63Ax/
74KqXunKNmwzuoNzfthNii5qH7AbvD2v5mMHqzaNf8cOviYEf2OL7JZRsvmR
4L6DrntTtnbBmsgbHSvsm/HtXs3YRVxd39ywu2n5lhaweVPm7z3Ajp5vfSe3
ja3hml2KvXqVW9kGu8hnQroVW7fb6VUE9to5664R7Ar55iOV2MsRFBaL2P93
PgL+/3zkf94buek=
      "]]}}},
  AspectRatio->NCache[
    Rational[1, 2], 0.5],
  Axes->True,
  Frame->True]], "Output",
 CellChangeTimes->{{3.398278957075845*^9, 3.3982789846055202`*^9}}]
}, Open  ]]
}, Closed]]
},
WindowSize->{506, 748},
WindowMargins->{{195, Automatic}, {Automatic, 140}},
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, 199, 8, 172, "Section"],
Cell[792, 33, 175, 5, 49, "Subsection"]
}, Closed]],
Cell[CellGroupData[{
Cell[1004, 43, 44, 0, 42, "Section"],
Cell[1051, 45, 1743, 46, 188, "Input"]
}, Closed]],
Cell[CellGroupData[{
Cell[2831, 96, 156, 6, 70, "Section"],
Cell[CellGroupData[{
Cell[3012, 106, 301, 9, 58, "Input"],
Cell[CellGroupData[{
Cell[3338, 119, 471, 11, 51, "Print"],
Cell[3812, 132, 470, 11, 51, "Print"],
Cell[4285, 145, 462, 10, 36, "Print"],
Cell[4750, 157, 465, 10, 36, "Print"],
Cell[5218, 169, 466, 10, 36, "Print"],
Cell[5687, 181, 465, 10, 36, "Print"]
}, Open  ]],
Cell[6167, 194, 232, 6, 68, "Output"]
}, Open  ]]
}, Closed]],
Cell[CellGroupData[{
Cell[6448, 206, 117, 1, 42, "Section"],
Cell[6568, 209, 4341, 103, 413, "Input"]
}, Closed]],
Cell[CellGroupData[{
Cell[10946, 317, 256, 5, 62, "Section"],
Cell[11205, 324, 252, 7, 37, "Input"],
Cell[CellGroupData[{
Cell[11482, 335, 374, 10, 37, "Input"],
Cell[11859, 347, 5435, 94, 267, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[17331, 446, 435, 7, 37, "Input"],
Cell[CellGroupData[{
Cell[17791, 457, 658, 12, 20, "Print"],
Cell[18452, 471, 655, 12, 20, "Print"],
Cell[19110, 485, 657, 12, 20, "Print"],
Cell[19770, 499, 671, 12, 20, "Print"],
Cell[20444, 513, 676, 12, 20, "Print"],
Cell[21123, 527, 674, 12, 20, "Print"],
Cell[21800, 541, 671, 12, 20, "Print"],
Cell[22474, 555, 674, 12, 20, "Print"],
Cell[23151, 569, 671, 12, 20, "Print"],
Cell[23825, 583, 675, 12, 20, "Print"],
Cell[24503, 597, 676, 12, 20, "Print"],
Cell[25182, 611, 673, 12, 20, "Print"],
Cell[25858, 625, 676, 12, 20, "Print"],
Cell[26537, 639, 673, 12, 20, "Print"],
Cell[27213, 653, 676, 12, 20, "Print"],
Cell[27892, 667, 678, 12, 20, "Print"]
}, Open  ]],
Cell[28585, 682, 175, 4, 37, "Output"],
Cell[28763, 688, 704, 14, 20, "Print"],
Cell[29470, 704, 499, 9, 37, "Output"]
}, Open  ]]
}, Closed]],
Cell[CellGroupData[{
Cell[30018, 719, 219, 3, 42, "Section"],
Cell[30240, 724, 7034, 162, 638, "Input"],
Cell[CellGroupData[{
Cell[37299, 890, 453, 12, 53, "Input"],
Cell[CellGroupData[{
Cell[37777, 906, 542, 12, 20, "Print"],
Cell[38322, 920, 621, 12, 20, "Print"],
Cell[38946, 934, 621, 12, 20, "Print"],
Cell[39570, 948, 624, 12, 20, "Print"],
Cell[40197, 962, 621, 12, 20, "Print"],
Cell[40821, 976, 624, 12, 20, "Print"],
Cell[41448, 990, 621, 12, 20, "Print"],
Cell[42072, 1004, 621, 12, 20, "Print"],
Cell[42696, 1018, 626, 12, 20, "Print"],
Cell[43325, 1032, 621, 12, 20, "Print"],
Cell[43949, 1046, 623, 12, 20, "Print"],
Cell[44575, 1060, 6895, 116, 238, "Print"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[51519, 1182, 456, 10, 53, "Input"],
Cell[CellGroupData[{
Cell[52000, 1196, 363, 8, 20, "Print"],
Cell[52366, 1206, 424, 8, 20, "Print"],
Cell[52793, 1216, 431, 8, 20, "Print"],
Cell[53227, 1226, 444, 9, 20, "Print"],
Cell[53674, 1237, 462, 9, 20, "Print"],
Cell[54139, 1248, 3628, 65, 231, "Print"]
}, Open  ]]
}, Open  ]]
}, Closed]],
Cell[CellGroupData[{
Cell[57828, 1320, 105, 1, 42, "Section"],
Cell[57936, 1323, 3856, 99, 338, "Input"],
Cell[CellGroupData[{
Cell[61817, 1426, 421, 12, 53, "Input"],
Cell[CellGroupData[{
Cell[62263, 1442, 496, 10, 20, "Print"],
Cell[62762, 1454, 476, 9, 20, "Print"],
Cell[63241, 1465, 496, 10, 20, "Print"],
Cell[63740, 1477, 476, 9, 20, "Print"],
Cell[64219, 1488, 498, 10, 20, "Print"],
Cell[64720, 1500, 476, 9, 20, "Print"],
Cell[65199, 1511, 498, 10, 20, "Print"],
Cell[65700, 1523, 476, 9, 20, "Print"],
Cell[66179, 1534, 498, 10, 20, "Print"],
Cell[66680, 1546, 477, 9, 20, "Print"]
}, Open  ]],
Cell[67172, 1558, 4697, 83, 268, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[71906, 1646, 411, 12, 53, "Input"],
Cell[CellGroupData[{
Cell[72342, 1662, 374, 9, 42, "Print"],
Cell[72719, 1673, 352, 8, 42, "Print"],
Cell[73074, 1683, 376, 9, 42, "Print"],
Cell[73453, 1694, 352, 8, 42, "Print"],
Cell[73808, 1704, 374, 9, 42, "Print"],
Cell[74185, 1715, 354, 8, 42, "Print"],
Cell[74542, 1725, 374, 9, 42, "Print"],
Cell[74919, 1736, 354, 8, 42, "Print"],
Cell[75276, 1746, 376, 9, 42, "Print"],
Cell[75655, 1757, 356, 8, 42, "Print"]
}, Open  ]],
Cell[76026, 1768, 4875, 86, 237, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[80938, 1859, 470, 15, 53, "Input"],
Cell[81411, 1876, 327, 12, 63, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[81775, 1893, 344, 9, 37, "Input"],
Cell[CellGroupData[{
Cell[82144, 1906, 386, 9, 46, "Print"],
Cell[82533, 1917, 364, 8, 46, "Print"],
Cell[82900, 1927, 384, 9, 46, "Print"],
Cell[83287, 1938, 362, 8, 46, "Print"],
Cell[83652, 1948, 386, 9, 46, "Print"],
Cell[84041, 1959, 364, 8, 46, "Print"],
Cell[84408, 1969, 384, 9, 46, "Print"],
Cell[84795, 1980, 362, 8, 46, "Print"],
Cell[85160, 1990, 384, 9, 46, "Print"],
Cell[85547, 2001, 366, 8, 46, "Print"]
}, Open  ]],
Cell[85928, 2012, 7833, 130, 215, "Output"]
}, Open  ]]
}, Closed]]
}
]
*)

(* End of internal cache information *)