(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 9.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 157, 7] NotebookDataLength[ 21933, 492] NotebookOptionsPosition[ 21457, 471] NotebookOutlinePosition[ 21800, 486] CellTagsIndexPosition[ 21757, 483] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ RowBox[{ RowBox[{"RK4", "[", RowBox[{"f_", ",", "u0_", ",", "h_", ",", "t0_", ",", "T_"}], "]"}], ":=", RowBox[{"(", "\[IndentingNewLine]", RowBox[{ RowBox[{"n", "=", RowBox[{"Ceiling", "[", RowBox[{ RowBox[{"(", RowBox[{"T", "-", "t0"}], ")"}], "/", "h"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"t", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"i", " ", "h"}], ",", RowBox[{"{", RowBox[{"i", ",", "0", ",", "n"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"y", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"n", "+", "1"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "1", "]"}], "]"}], "=", "u0"}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", "n"}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"k1", "=", RowBox[{"h", " ", RowBox[{"f", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"k2", "=", RowBox[{"h", " ", RowBox[{"f", "[", RowBox[{ RowBox[{ RowBox[{"t", "[", RowBox[{"[", "i", "]"}], "]"}], "+", RowBox[{"h", "/", "2"}]}], ",", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}], "+", RowBox[{"k1", "/", "2"}]}]}], "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"k3", "=", RowBox[{"h", " ", RowBox[{"f", "[", RowBox[{ RowBox[{ RowBox[{"t", "[", RowBox[{"[", "i", "]"}], "]"}], "+", RowBox[{"h", "/", "2"}]}], ",", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}], "+", RowBox[{"k2", "/", "2"}]}]}], "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"k4", "=", RowBox[{"h", " ", RowBox[{"f", "[", RowBox[{ RowBox[{ RowBox[{"t", "[", RowBox[{"[", "i", "]"}], "]"}], "+", "h"}], ",", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}], "+", "k3"}]}], "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "[", RowBox[{"[", RowBox[{"i", "+", "1"}], "]"}], "]"}], "=", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}], "+", RowBox[{ FractionBox["1", "6"], RowBox[{"(", RowBox[{"k1", "+", RowBox[{"2", "k2"}], "+", RowBox[{"2", "k3"}], "+", "k4"}], ")"}]}]}]}]}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", "y"}], "\[IndentingNewLine]", ")"}]}]], "Input", CellChangeTimes->{{3.664889161346874*^9, 3.6648892186583633`*^9}, { 3.6691041921064377`*^9, 3.669104263666531*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"\[Lambda]", "=", "1000000"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"f", "[", RowBox[{"t_", ",", "y_"}], "]"}], ":=", RowBox[{"\[Lambda]", RowBox[{"(", RowBox[{ RowBox[{"-", "y"}], "+", RowBox[{"Sin", "[", "t", "]"}]}], ")"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "=", RowBox[{"RK4", "[", RowBox[{"f", ",", "0", ",", "0.0001", ",", "0", ",", "10"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"plotAppr", "=", RowBox[{"ListLinePlot", "[", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", "+", "1"}], "]"}], "]"}], ",", " ", RowBox[{"y", "[", RowBox[{"[", RowBox[{"i", "+", "1"}], "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "0", ",", "n"}], "}"}]}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"c", "=", FractionBox["\[Lambda]", RowBox[{"1", "+", SuperscriptBox["\[Lambda]", "2"]}]]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"plotExact", "=", RowBox[{"Plot", "[", RowBox[{ RowBox[{ RowBox[{"c", " ", SuperscriptBox["E", RowBox[{ RowBox[{"-", "\[Lambda]"}], " ", "t"}]]}], "+", RowBox[{ FractionBox[ SuperscriptBox["\[Lambda]", "2"], RowBox[{"1", "+", SuperscriptBox["\[Lambda]", "2"]}]], RowBox[{"Sin", "[", "t", "]"}]}], "-", RowBox[{ FractionBox["\[Lambda]", RowBox[{"1", "+", SuperscriptBox["\[Lambda]", "2"]}]], RowBox[{"Cos", "[", "t", "]"}]}]}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "10"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", "Red"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Show", "[", RowBox[{"plotAppr", ",", "plotExact"}], "]"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->CompressedData[" 1:eJxTTMoPSmViYGAQA2IQ3ecXIG7g99rRhOkqmL7mLq4Eop0MfhiBaLvH+8xA dBVnkBOIVpxu7AWi2e6dA9Mpj9cEgOg583aD6RuLrUJB9KZnpyJAtA6nbgKI 5soySgbRJfXVGWCatx5M35cpzQHrm2OfB6LlTnzcA6J/xk64BqJ5jiVrGIL0 P7yhDaLvnAgKAdE6y1hjQPSeexwa5kBapNpBB0Rv0o00BNFBeY1gWs/v9VwL IF0zk3k+iO6aG7kWRNcqGW4D0XIVibZloa8dVfY52YPoJNY+PxBtlvADTCdP PJgAohXcfoBpCYPwy05hrx2fnSu9DqJTQn9yuwDpSUquIiB6zo6amElRrx1P uKungmgAOGCaCQ== "]], Cell[BoxData[ GraphicsBox[{{{}, {}, {Hue[0.67, 0.6, 0.6], LineBox[CompressedData[" 1:eJxTTMoPSmViYGDQBGIQjQ50nWVePzKTslcwN90hxyrgAOFr2Tuz5zKsmJ7n mGQUaqC10tg+ec1kTQazs04QeSv7zxyvVwU80Xb5s/LjJd8kB3vdPMNnDxZ2 uULUO9tvv544mSHvpdvxXTt62T642c8XDEpg9/DwgOj3sk9asGZKiP5yz8mn 6z321/raR1QJ7fJTZvOGmBdg79ph4pGsnOqjn7eYcQ9rkP0cvu0hbPpHfCHm B9v/NdrstslV2X/qN7b4GT6h9vdPr+l9ntYUALEvzP5i95VonUkPA3/VZe0p mRxhH55ekMxz3CEYYn+UvdZZ1+kSrAtCEtjPSQbeirYXu3309tnp/0Ih7om1 L9zn7a9iFht+tN+wXFcx3r7+wWG94ie7IyDuS7DfsY/zNuMiqSixHq9XLCaJ 9usM7LgV8iqjIe5NspfekZaj6nEjxiPgj0Tx9SR78e9proX65nEQ9yfbb/ij KJ+lPC2+WmSd+8OqFPvtsZMLapW/JkD8k2r/yGqTn7x+SNK6G/Fl/nJp9jPT P8mfct2cDPFfmr3Lfcvgf2lCqQ/nCC7dezDd3m6WkpTTpMI0iH8z7F0XJqkp Hr+QLpJ4+LJ2aqY933fTLB1Wg0yI/7Psf2n4b/88vS/LTbWUaRZHtv2qboNS D7O32ZDwyLaXzTu+aeoT79zKl2qGHGty7Cdc57qgtGhVHiR8cu2rqhyb7fM4 C9asBTkwz96td0GFs0dGISS88u2XnHP51qt/vOh+YVffk0/59vmpwS2tymol U+02TAbxT9llPE1bY10CABimGFw= "]]}}, {{}, {}, {RGBColor[1, 0, 0], LineBox[CompressedData[" 1:eJwVlnk8VO8Xx20pWxhL2WZGhFSUVi3OiUqhiKRUKpWUvfhGaVUKyRJZEiUi 2ZIWwkUKSRJZsmRfhplL1rH+7u+fe1/v1733eZ7z+XzOmVG2dTE/y8fDw7OB uvz//vS6reDaMwF6PKqZLaz+g3olNUwTK+Z5GFBsbthKPw/zXwt8mEwvWCM0 UyTE8IKnhYGMBUw/KH2ZHmPJ8APH7R5OE4womDvfrSxO3Wum7cOGGK9gayax f4vSK1jY1KHOZuSAom73WX+ZHHj+Lc5tgFEOdocPG6SMl0F6hl+To9Yf8GMm koIejXAOqqMGL3aCzOa7rvO7OoAnfjy3KaAT6nifstysOyBSQLH5W0IndEbJ Y6NzB5SV29Ff1XWC099CY9+IDtA4OBN/dksXtDwMeHC4rwP67NVTW3m74Zub NW/NvU64EHqt4GdwD4SmOj2uyO8C5+4Vndnp/ZA+LvBTn78X2mwaFxWV98PC a/1/F8j2gkXDfa3Krn54vTLM7aNGL2yu6PXqlmdBvm1j18C+XuB/kyghe48F I22LIv0jeiHSW1nv8okBSHB3J56o90GRlFyUrgQbMhvfCGZv6Yf1gWUFu1ay YYnrnNV+4354KejZdWA3G57xbY6oO9oPD7h12he82eCQjpmZ3v1w6G/Y1+g+ NriWqC7wz+8HVorE6HQhB858W7iidCsLhJXvqP/dPATd/W/vK2kOQNm2EO3r u4YARE0+OG4aAN/DsZuUzIeAccZ1PGXnAPAFfzQ86jAEqnNud0dtBmBqbvBc /dMhyNq6ie9T6AAMNFsmV/ENw+1lSZ6LJwbge6T6isLvw/CE3hUTlTUI/tnr 155oHIYjNrzvPuUPwp6fO3Tnuoch4/l3VkXZIJQsPLZ32/wwVAqQj963DkLe f8HnP679B/QFp+jZQmxIO8h9lfn4H1xpfrvH6BgbHkpUaD4/OQI/vYpY70fY kLNJ6d1+5xHwujEksm+WDV02LjBzdQSkFYP0KwU5sDVN6qBVxAjo7ztbGyDH gb69x66J/RiBI9I833n0OGDgw67y2joKbu+Xm5Tc5oBzClir7x0FE9NDLyf9 ORBdHdJVe2gUajWd02UfcWCIsZGrfXEUUm+8yhR/wYHYvOvLepJHYc2vLU2j RRyYHBP3MJcdg4ZnE2zfaQ6oKNnOz6uMwXqXiNsWfCTs35ntl7Z2DNqFup2E hEhIDD0cu2jfGNTazaxfK0uCufbz0gKfMUhc1Dbvok3C9UMjBxxDxmA6Va/E dwMJr67tapaLGwPHw205vltJ4PneP+SeOwYJfpfU9Q1JSLPXkV85PAbsG9rh L4+T0BB0J6FhbgwekfccBU+TIPChTstXdBwe7eu0MLUnwXrBVYN29XHQrESD mIskLIr/7BhhMw7cR2HblvmSsK5cZmKn4zhUsR6s/OZPgs3QuVv/vMahqM5T 7GgQCe/0RB/vCx+HS+c7YzZFkmDbdJDg/z4O9sWfgpuSSRB+0u7k1jgOmddP hX1MJSHL2lnxbw+1/lPrPTcySeD/c88rl3cC0vjfvi//QEJCQ+46t00TwL/j lJ35FxJMIg07WndOQPBIZKJ3GQmjVrXBJuYTcGJD+K3gChJ21bPZ6k4TQByG hR7VJPT8Zia3xk+AgbefxNpmEh6Gpx0yyZwAq9SXM0mtJGy03LIgN38C7pSz PBa1k3Cv1sI2vGEC9DfsEvLtJmFNWJsEf88EDJWsZb7opfSzcCJcRyZg7dfh /Ff9JGjU+CqaLJ6EbTm7719ik/AzVKoiR2ESHgStydpIkuBp/sxLfcUknIyP ON45RAKTtlojfOMkEHP7/Lz+kVBenVPHt5Ni/w3ruSMkXAzZfdf1wCTsF8Tj tmMkyB+oWddqMwksERf+d+MkXPg5GJzjNQnH7IM/y3JJoAV7gfq9SVgFv5uW T5GQayrICQubBK3nnU70aUp/8UcxfPGTcOGyqBfPDAkiVQxj14xJEPTLnKmg OPthKrclbxLeX8XJm7MkHNuvm2z8bRLUarc5MOaofCz+eiinfhIOjG+wTKI4 tdJ8gXr3JPxc/v3d0nkSDgb+fRv2bxLkfDMeXKJ4xsTRlo+HCwXyzxo/UJwg ypVwFeNCZKlTTDfF+77fJVrkudBrNlk7Q/F4AM3ZWIML549q+8xSHGscp5iz gQt2jrSkXop3i6yqUDPgwtv+TPxEMefbR68wMy4kymmaeVL82H+XBp8NF8aX pf2mU6xn9KvOxYELaaJe1enUeXuETtxt8eQCJ7/ZQIPioPKBdca+XMib3bg6 gKp3s59nx8dHXCi+OB5YT+nRtmdBiNpzLszGR9mJUXx/USiEpXPh9Ej6h1WU nmvK6BzePC7cTKq9tYHSu+He6xiXci6ItX/8ok75cdNws3FLHXX+OLk7/JNU PhZ+4Rp1cUGaKMkvo/z7+fVA8sdhLmzgRly6TPnr6dt6SG2eC387br4UH6Xy sdthQZjoFFwZ9LYKofLh9uWOrYv6FPicHdl0gMoTvWu34Pi6KfCoMDmSTuWt gk8oxRunoGWFYv/EAAnLMXDY78gUpCQdPWzWR8IvG9NwCbsp2L9KaMuJHhJu XJPUjbg4BbcuJj207iKhMTf8ZkLAFIy3Z4fLtpHg22ilujJyCh6aFBnVt1D9 PilX9iZhCsq4k5fvNJEQuCFWnMifAgOzZEyoI2FH5svYRnIKXPc0vONS/cf5 Ya9/cmYKYg7kHaKXk/CErdnTs2gaREeue2t/JWFMM2P1qPI0CKQ3H19SSM2v xA95iw9Og8+QXtpsFgmHSrxOhp+cBqaIUIpvBjUPOrcKKDpNQ9+qU1tnXpNw nFlovMJ3GlLDPFTTEkiQfFL2x+DjNBB7/pqqhZNQkON//VvJNNiH/7dGOoQE hwaTZQeqp2GT8K6PnAckfJGpPm/DmoZFd/JLz98hwSu4cdJTcQaq32it8aTm Xftdlmz6zRn4vsI5dtKImhcJqbnrA2fgq47J9qxdJGz97GzzKWoG5BrC3hxG EsJ5RhLLsmZg5cYgdUdqHu+9Or2+s2sGtnc2JqXRSXjjKmqxdO8sFPV+koll c+B7tJRX2qFZ+P4pXDO6hwM9JfJx+mdm4SU3furGXw7Iy60YcLwxC0bNbW2L qjlwp2iXT9G7WQhwWXlE4i0HDtFuvL2wbA421DmWTl3kgOs238Z57Tk4lPRy afgFDvjbBc6HbZ8DvX+XrsvZcoDIfWJMHJ6DwTiTZ21mHNA4/bFTKmgOlNP4 Bu1Xc2Dq7bBU/jTFYo6m5zrYIN06qWsuNA8DEnX9iQ1s0FrEc7JXdh52Bs/x Vv1gw+lji1MldeYh22eitC2XDd8FVu60s58H+5AOd/5QNsRannEX/z0P0uoD 9cbb2PD8qgNrjRUPPhJ1WMy8NQhrCoUm7x3nQW2bBVHBHoNQKJC84O9pHvyP +8OAPD8IbYFdzEBXHmyXX7fksvkg0J8dt+r340EeyXt33iwbhOgS05LneTzo Eqm6tIEYgFCx9bG0Zbw4ZZh5VIXNgttxMwdGB3nxr6iAY5BMP5it+ei2coQX RwIVpJcK9gO96FKILZcXhVyvPw8Z74PcDtbPnwv40MP0dYhVfR/8W95omkbn w2ZL/Kkd2Qe2qe/22ZnxoT2jXkN+SR/syHHe2/CWDxlCv1YpCPcCT00b5l3h x+upE58lS7tAzd3O86DmAlwyRTsE7q3QNRZxSCRfEPuVospetP4Co+4L10Ld FuGqBne1tYLvYX9X8Eo6XRgfJep4DLoQBOuz7c2yNBGcnHvZViVYS2hJ1ikZ nBBD60U87l9WtRE9gT23w8+K4ZITdPlQ/TbiqfBEb6+DGA5E/Yk8eLiNEBVY mvXAUwzffGnoz/VpI1jjR3bXhYjhX/6FJ0T+tBGJzS3O50vE0L35VG3GnXZC MbmbCFqxGN+ZGW5Pqe4ghGHsVPPIYswULWI/NuwmPHdoLHo/tRjtrHtnFxzr JnoMjqYH8YrjwTd/z5x17SaK9xRN6YuLY+3q0OnRyG7iinngo1ea4jgkHZlq 2d9NDJxd/uW/U+KoKnxsj8S9HuL7A0sNySpxbE08sV7wfS+xJej+D9ZvcZzv E5M6U9ZLJId8ci9pFscKnre3Mv/0EnceKxddZonjH5ZBo/R8L0H/Ob6K5JNA q3CeWSHDPqKusfiLJV0C30ddfxpd00cYsq0nlllKoMa4yOqZ1n4iVtqn5sZh Cexj1m0s4/QTo1tfZzQflUCf7LTBW3P9xDP/6XOPbSUw7oLSmgIlFjGlHtMg 5CqBq7T4Og4fZRGpti05Q/4SmKQw7y1WwyIkGk54FxRKoJBkTLxm9gBhN3/v sEKJBGpKmyVGFA0QeWqZ6z1LJfBvsprw8I8B4rwHD3vtDwkU8XxQa98/QBTT nh9PbJLA9ofe5xyUBgn3fR16D8YlMDtoIDD71iDRUHyG13qVJHrk1vDLb2UT Qzv+yF/TlsSTHQdPG+5iEwuLTNc/06Ge7/1eZGvKJjYSW871bJbEaYkK29On 2UTYJ4nKi7skMX4yyy/An02YZudFBthIYqCc7IaLv9mE3TqdrPRTkuj7zyrL u5VNXM9Kqqg+I4mkg8VP9142kZoZOrfEQRJ7jdwaVnPZhHCa/ZmEy5LYZBF1 q1iBQ3xNlF6THyyJrwId49yOc4hWVf+9bY8kMcxL7rXFWQ4x9mLelj9CEmeL gwRUnDiEajwrfO9TSXSoSQ308+YQt2MLp3+/kkSN0Rl3j2gOoRfhWMoplsQn gbeflfziEIdkO9okv0rixjuhp7waOYRTuNXU+nJJPPJvSkOhjUPEPNJffbWK qqcgKHc5m0NMBS19tLBZEiPuu3kdEiSJd/dLbJhjkpgiXiWpsoEkko0qlt6b lETBuEcnFbeSRLTor1/saUm8lxggP48kcTP47+5PfDT0zmUSN01IwiRySuuQ BA0/K+ZbuduSBFjz9udJ0XC7od+Dj+dIQkdx0QuVJTRc/NJQtMeRJJY8k1ky rETDmZY3LTP/kURH0pr5gJU0/N3kflrOnyR+n9+U80+LhlvcIv2THpJE2Uq9 S0d0aPgp7V4b/RFJpGUY96rp0lChhyNS/oQkPD/YVRXtpuHlvLvCIqkk4eDl 5KdhRMNtdzcf6MggCZut7gZB+6j912emx78lCQPi1odjB2louW/5i9Zcklhc GhM3cYqGbzbs0TQoIwlevxdHbM7SUF/0Z6ttBUmMGqVIfbGn9l+55bHDD5Jo /PHhXqgLDY/dmhfYXEsS34MLdnAv0nDKfsXH6TqSKDD/Mn3iPxq6zf86m9xI Egl1v1xWX6Ph3tfc5LxWkoiIbFwRdpOGP59826LaThL+1m2dUz40XD5aTLh3 koRLK9uqzJ+GOS33Y372koTts1FJ7Yc03B+rPdzUTxIHbacrwkNo+EI9U+fn AEls6VmEpyNp+Nct+b9LJEmsThafKn9Cw8I/1/9bNkwSzAuy2WviKD9sBo/n /iMJ2iol54h4GtaIfVu5dZQkFnBUNOYSafg1nfdvwhhJTGZodpx5RcO23MD/ JsZJguW2NqYilYY6p7yGdSZJomXd5kM6mdT7vm37rbgk8XNMTyLqLQ13lJKB Z6ZI4vOHXd/m39PwpkdXmvU0Sbz3Mrljl0tDZoxouu4Mla+tFnqV+TTkm/sa yDtLEk9mj0yuK6JhkayPSRbFD4mTWdElNKQpFPUbz1F5u3XOkbeMhqu4kmd/ UHzJwFnNvoJ6fmRJge48Sdgt8Gj78YOGtTE+o4EUyzPK66eraahZfCCrlOIf m5WqNH7TMCW6bSeb4tvmbl8tG2j47fLS1DmKNzp+yb/dREPXU5zeaYpZd+Xe ZbTScELmzGQXxbFxTqnN7dR6wvdbcig2zyl6IdRN+XneI+IKxYI1Mk829tGw W/qgmjrFuYPnQ08P0LA+4bA/QZ3XWbDAL5hDQ+lzJZ8NKF7GpN3KH6bhgxPd v95S9dbp2nmyRmnIMZT9KE6xv0Wuy5JJGl4PKrhkReml57T43M5pGpZamwsG UHr+87W1cZuj4SEfc9cUSu/EZ+8tY3ml8EKESdY7yo/DucL7KgSkMKrxz/dU yi/RWpudkwul0BUTCoMmSMJ94cJ15oul8Fn9541ylN8aykc1b0hK4SKPpFxi hCSat2Qop0pLIUi0y5tTednpbCUhqECtT8xc2kLla/Le64Xr6FLIX1ntHMwm idTn83MnlKXQ41vH3l9UHqV/J7Fz1KVQV0ozUa6PJLq2Tn5zWi+FCmZatA9/ SSLS0qQ4epMUbrl5YY1jCzU/XJ7llG6Rwk8rxHXFmkgiO35PsvIOKfxWFyku RfXTHaGou7/3SSH7cfzcaqofN6uwvfkOSOEQqyBGn+rXwW073LUPSmHd2Wgd gy9Uf7j22/pZS2HsFs19iwiSUK3Xxe3npNDzV4ZH/BuSaBgK3HThghROnHYR WpJOEg+EO7QinKSwqMYx8XIK1e/b/ZWGL0nhKj6z+Zl4kihJ+DOVcFMK59IK ri2l5s/pi1ffiUZLYWN3Uk+YC5UfndZ3K55KYU+dv3PnBZIQ/ofvdz+j6jFZ xaHbkUSmm+CHmy+lsCnfoO3YMZKYcQ35OJolhZZWnp2jhiTx2CXpU3OFFG5X NV4RqUAS57WF87g/pLB2c/WzDzIksY10zJP9JYW7dNMNSsRJotNZJ9+sQQpf v+uUS+QniTXO+QUlXVIodI3rcpfFIcoda4pSZ6WQWDi/vCSbmv+rNhaX80hj QbSd4sc0DuE6GFncwy+NiRuWX4t6ySFkHW0+M4Wl0UV35KFqJIewdegvCZOV xl+fj7gcu8IhZs7Pl3prS+M+t/2qVls5hPa5lT9MTkmjWa9lRnQGm2g9z3BJ OyONdtmhD70S2cQDRymJxfbSaL6m28jgCZvod5s+UOUsjV9iuJcjfNnEC++K 3we8pdH+o/Mlg2NsYkmoQ8uhCGnU3yzzNJOfTcznvR48WSmNCTdldJYaDBLp RFxg0U9p9LevIE9tHCSOFT/SWlYrjQO9+4MjVgwSOaVXXTv/SGNWw+7Yr+KD xKVfxqNn+6RRpLdmaH3TANHbOzB9gV8GcyQq3lk7DxBVtFUiHroySOYGVIj7 sqj+cup/sk0G+6JoRwcvs4j00vTSYpDBxzsirn48zyIeeq+9I7FbBi//6A9l 7GMRJj0bZ1+by+A2CUVpTSkW8S1nB9nhIIPveaT2bY7pJ76ctKoxi5XBP+sI wavxfcSn9NtPVvHLYn/Rxk8lt3oIQRNL+QhBWRRREHAxde4hDvSrR/EKy+Jz z9av5dY9RK9K5eM6CVm8GZnzJHBdDyEVtST0Jl0WLU/HOFh0dRMOPqn3a3Vl Uc+YYdWl300oHK7z8HaVRTNzaVOf0U7iKq+m2fcWWTzlvGOV3Op2Qrj5yh7n 7CW4M+FszZnjdUT3S+WwOO+laOq7+OX9sBIiTENOIMdSDoOnHop1N72C4Gi7 0lNMedz+LWhmr1cVWIcd/yHaK4+a+lkFLXZNMLkzkRtXoIAXinYvEXbrBO9P cT4lRQp4pLDu4MydTuDRiRbrL1FAwZOldzsiO2EhI0hZp0IBJS5nBPsXdoIs 13NvSb0Cbgu5NWoo3gXr0kyi+oYUUNbW6JhTShc4y4xuXquiiMvtTO6K13fD vwBOsaWaIlayeCbs+7vBg6/f5MoKRZzvEerInu4Gb7Ll5GdtRfyiNTquwuyB gPLS+5bbFPH4xEz7cfseSL72pMHLUhFDBSwWbhrvgc5ufc/ie4poVCn8Smdh H+jH/YjXDFDEsdkoEU25Pnh22Loy9KEi3l1e9098ZR/YVLgtOx2uiLPhFxWS 9vfBn8xn3/lfKGI155cjb3gf/Lo6x9xVoIgS06oT1vR+KJbMLS8dVcQ8z4lg YxUWKFfsGtWeVMTXLk+OTeuw4OadanrktCK2RHxRi9Rngd5E3yV7PiUUeu3o F3mKBbnNS+hCEkrIU7N//9tYFmQleVw0WqmELd8ubjSUGYD47WsVKk8p4ZrI dediRwageL+wz/QZJWyP2urJnR+AjhOdrBX2Shj8el+NvuggqNwOz/V1VsI9 snf3PFMdhISv3CN4VQm3r7hw+O7BQUg0LY58G6aE8cscvIczB+HrySc8HRFK uPG9ZvzqvEHocXO3l3iihAFDfdFWpYOgFqa22ek5xY2r+91bBuFlg3+9WroS 7ioc31C7iA3Jpyxko0uVsHXzpTjGMTaUXVx1veybEtZXCREb7NjQ57OgZ7xS Ca+tp+3WdWWDxssP7yxqlZC5XMtt0R02vGIpWIq1K+G2vYzArylsSLnU9ejm lBI+PqYelDfKhrS7HpLnVtPRD8zKzl3hwOC247fvrqFj4NGnA5O3OLBqdOfI i3V0rIl/efuyHwdSbKV/t+nSsWQlwd0WyYFkeBtpvZuOpmpbg/qyOZDAHWLs P0FH9fXLHxewONCV2RDsaEtHhxCdAdVhDqjYF/IGnKUj7/Kzrp4THIivC+os daDjir8d0z38JMS91Ura4UlHWvrJ9HYFEqKdHLU2hNBxE136pNEeEv6oHoyz CKNj6m1+/tv7SJBr3ipxMYKOyTeGnr40JyHSSORf+lM63p81yXh/jIRw9ZR3 Gil0XNh7bGDYhYTa1hA1wzQ6hk2v5wa7kyD92CvibCYdtwcY9yp6kRAqsNcr /j1Vn9baI9O3SAhu792m+JmOcW2We3+FkvAzsipty1c6Lq0yXNzwmARxsw/0 I+V0TFn7PPZbNAkPC3x5HlfRUcxhZsl/8SQ8eLK8RLyZjjcaeFf3ZJLw3Vxs g9ZfOrq3+ouaZ5MgIjyWaNJBx8rQoMSUDyT4e5bc8+uj9IwYnVYuIOGe5Wlj /jE69lXJCaz4RkK8afOxs5N07MpZ8mZjJQl5ey2dS6fp2JH/X+zqnyT8224Y EsDHwMrfg2vrf5Mgurkwnr2AgXZz4zV+DSSo6+hmmwoxcJfTnvTlTSQcU1tZ LyXBwEJvawG5NhIuMxP6PKQY6Pgw2Nu1g9JDXmmqXpaBfBWTy990kVC6WFwp RomBsTJzCzh9JHQsuq81x2TgU7lHuiwWCTN8vHhSlYExG40iqwZJ0Jn4Z6u6 koHLDIurLYdI2Dfs4O6rxcDNDbOvxodJsB/outu3loFImCXfHCHhdvfxCKMN DCxTe14xMkrC0791yambGVga2k0zGyfhY6Np7uJtDDRpV70WPkHCr5qyCldg 4Hbj/4RKJklgV+5o+aXPwH9neD80cUlYVJbLWb+bgQ30uTvNUySoFK/jidjL wPkv2Re/TpOgl5cqyTVhYJTA9RtRMyQcfr9c5agZA//rKE6xnCXhUmbs+nwL Brp8mRmZovhhypLdDCsGlri9s/abIyE5IdjqljUD3YXvtfLMk/A5Vuh853EG +koq3jhFcWvk7Su7TjEQ7LW3vaaYGzodkHSGgSdMH8i0Uiwd6P5UyJ6BeucW i0xRrHWPne7gwMCejOeK8xTvvWVXWOnMwBV++kbU/204c/VvtfZFBj7Jmn1U QvF1j8OdIR4MDDYe4t6hOMqlenTEk4Gzvwy9tSjOPm8keMibgWfU9RWKqPNW nf685OMNBsb7rq7fRjHr+LYV8j4MzAq0z3xO1Sdw+N0Wb18GXtpsm0RSejDM tUxa/Rho2h6Vp0bxFpOk4xjIwAGMG9pN6XdwN9MlPpiB2U60nWaUvi4YdVMg jNKj6cUHfUp//y20ULsIBv62ljahU/4krg94URbNwKqlG/g6KT8LtQTeacYy UDG290/wGAljy8brOQnU932NE0n/SJBQcuk3S2Zgjo+WniSVF80lfVNZrxlo EaL3+jRJwkmRP0qXsxg4dd157vMACVcWWGg3vmOg5Pjy3up+qv/nK3BrDgMz RMqmy3pJ+DaSf3qeYKBIX3r6xU4SNjQ/f3WvkoHXtlaey26k9tOyn8z4yUDe 882juvUkDN7QMmyoofQpHU55XUtCgsqnLo0/DMzXjo60qKL8c6hlfOth4Dmm z6t3n0kYyot2/tfPQPu8xekxhVS/Lz6VL89m4JbvH/pd8km4k8W2dhhhIN38 8OKm9yRMefMbC/Iw8Z04GJYkk9CXaHcxmpeJIou/64klklD3ozxKi5+JK7Li r+s/JyGLGdx3SJCJqampbW5RJFwoUfRNEmVin+ipS133SWgU2Vi0R46J9/Zx VBXOkPB1fVRfszwTB0Jc6lecoPJwfEbcTZGJ5je/zqpaU/Mto9gmisFE0Saf rDZTEvZYmM70L6fW214h67yF6rdo+00BOky04a//2ilCwsvPFTaM9UwMeem/ 59MCEsIGtXzfbmBirV2lx415DrjqjdU2bWbi55Twjsp/HNDouHVxFTDRL67P 3aeBA5ErnqRVGjOxujlN2fIZB0ZHdif67GOiT2HoD2vq98S04F+MrikTW96c XWoUzAFBC6MHieZM7NS7Fd19gwMe3twL144wUW3qF6+kDQfMq6w0Vp9j4oF5 92Y9WQ6kR/EzO+0pFpHBdlEOCJ3JWBJ1gYlH7r556cTPgcJJwYULnJmokBEV uXOIDdrL3nc3uzNxQpyoHShjg6iH9IsHt5mYcTxUXO0/NthDYbT+HarexA5+ Rwc2fBZyDJ28y8SX0RZuMSfZcCX2860zfkwsjjKuzzNiQ3/pxZPbgpnoXktm ziix4at8tdLgUyZeUFzXeKNwEG4UBkYa5zDRuqNAoXpsAFgthWHcXCby/A5Y +65vAA5NjwQn5THxceXOyLtNA7Bqo7UffyETHxV+ujZQOAANr9WufPrKxJpv /pnf/QdAJ4I4urKWieuO6iiz5Qegx2mYLkIyUXVBucXgKhbsU7BM/KaqjEf5 djxE+T6I7NIqPKqmjHJhD7M/CfVBV9qipkF1ZXzarSqvwe0Fb8yTEF+pjEJX eipq63sh9ayKt8VaZewN+dl3PawXhDOHzZu3K6NugwqdIdILpTsf8nAOKePy EIt/jexu2OFcelTSTxmnrqb+Tg7ohDhFW9XvY8qY5l3+YNf9JlD9nPR4/9Fl 2CnjdIp/SxUEry7W/FWwDD2YPGInk1LAbkSgulxTBfVeBOt6bC8hbN/8Y/oF qWBvkPx8uGEd4XjbYsYoRAX3Clcd2+ldR3hYZNeLPlLBsonPNu2ZdcT9MY+H wY9VcJ3egrQeuXoiXZc7HfFUBXdIzwjID9QT3KL5upevVfCp2bVTPP6NRPAv 0cCSUhXMfqxQkf2xmYh+4WTvW059/1lClsVqJl64/zDYU6GC94plFwsptRDv ZYOnKn6o4Pjh5xf4b7UQzdZS9jW/VVClLSK227CVUO+UM+joVEH7gfJLetV/ iTXZV+gJ3Sp45fCnvXM8bcSWu03cs70q2KRcNGcm10aYqD/N7GepYEj1HrHq vW3ERQcmfXhYBS09ezJ1U9oI7223uFkjKmjs86h1XXEbcVeso9Z9TAV3n96j qPynjYjMeBEwOamCs2bD12qF2on4mwLncqdUsOoE7Vq0cjvx+sBZfe8ZFUTb cH0L3XYie9lXJb05FXytr5k3Z9ZOFIyocefnVVBHO6A+xr6d+B+3p/iE "]]}}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesOrigin->{0, 0}, ImageSize->{1061., Automatic}, PlotRange->{{0., 0.0047}, {0., 7.925371344592115*^251}}, PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Scaled[0.02]}]], "Output", CellChangeTimes->{{3.6691042993455715`*^9, 3.6691043168925753`*^9}, 3.669104575075342*^9, {3.6691067869358535`*^9, 3.6691068488863974`*^9}, 3.6691068856284986`*^9, {3.671364330300709*^9, 3.671364397524393*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"plotAppr", "=", RowBox[{"ListPlot", "[", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", "+", "1"}], "]"}], "]"}], ",", " ", RowBox[{"y", "[", RowBox[{"[", RowBox[{"i", "+", "1"}], "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "0", ",", "n"}], "}"}]}], "]"}], "]"}]}]], "Input"], Cell[BoxData[ GraphicsBox[{{}, {Hue[0.67, 0.6, 0.6], PointBox[CompressedData[" 1:eJxTTMoPSmViYGAwAGIQjQ50nWVePzKTslcwN90hxyrgAOFr2Tuz5zKsmJ7n mGQUaqC10tg+ec1kTQazs04QeSv7zxyvVwU80Xb5s/LjJd8kB3vdPMNnDxZ2 uULUO9tvv544mSHvpdvxXTt62T642c8XDEpg9/DwgOj3sk9asGZKiP5yz8mn 6z321/raR1QJ7fJTZvOGmBdg79ph4pGsnOqjn7eYcQ9rkP0cvu0hbPpHfCHm B9v/NdrstslV2X/qN7b4GT6h9vdPr+l9ntYUALEvzP5i95VonUkPA3/VZe0p mRxhH55ekMxz3CEYYn+UvdZZ1+kSrAtCEtjPSQbeirYXu3309tnp/0Ih7om1 L9zn7a9iFht+tN+wXFcx3r7+wWG94ie7IyDuS7DfsY/zNuMiqSixHq9XLCaJ 9usM7LgV8iqjIe5NspfekZaj6nEjxiPgj0Tx9SR78e9proX65nEQ9yfbb/ij KJ+lPC2+WmSd+8OqFPvtsZMLapW/JkD8k2r/yGqTn7x+SNK6G/Fl/nJp9jPT P8mfct2cDPFfmr3Lfcvgf2lCqQ/nCC7dezDd3m6WkpTTpMI0iH8z7F0XJqkp Hr+QLpJ4+LJ2aqY933fTLB1Wg0yI/7Psf2n4b/88vS/LTbWUaRZHtv2qboNS D7O32ZDwyLaXzTu+aeoT79zKl2qGHGty7Cdc57qgtGhVHiR8cu2rqhyb7fM4 C9asBTkwz96td0GFs0dGISS88u2XnHP51qt/vOh+YVffk0/59vmpwS2tymol kPArsF89UWjmdOXWUk7peXGapwvsS/3PeFjoPymDhGehvam8/8+Hrs4VSn+/ lT6YU2hfp/BYQjR9USUkfIvsd03ovhE+ibHa5r5/7/S8InujnryXpscTaiDh XWT/eYWdrQPrgdrQgyuW+DkU22/xZ58sOUOuHgBRlUi8 "]]}, {}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesOrigin->{0, 0}, PlotRange->{{0., 0.0047}, {0., 7.925371344592115*^251}}, PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Scaled[0.02]}]], "Output", CellChangeTimes->{3.6713643454772434`*^9}] }, Open ]] }, WindowSize->{1350, 631}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, FrontEndVersion->"8.0 for Microsoft Windows (64-bit) (November 7, 2010)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[557, 20, 3288, 95, 347, "Input"], Cell[CellGroupData[{ Cell[3870, 119, 2498, 75, 230, "Input"], Cell[6371, 196, 13216, 223, 659, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[19624, 424, 454, 14, 31, "Input"], Cell[20081, 440, 1360, 28, 230, "Output"] }, Open ]] } ] *) (* End of internal cache information *)