(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 8.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 157, 7] NotebookDataLength[ 37777, 782] NotebookOptionsPosition[ 36625, 740] NotebookOutlinePosition[ 36963, 755] CellTagsIndexPosition[ 36920, 752] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ RowBox[{"Plot", "[", RowBox[{ RowBox[{"Cos", "[", "t", "]"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", RowBox[{"10", "\[Pi]"}]}], "}"}]}], "]"}]], "Input", CellChangeTimes->{{3.508470183098377*^9, 3.508470257467126*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"$RecursionLimit", " ", "=", "50000"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "[", "t_", "]"}], ":=", RowBox[{ RowBox[{"y", "[", RowBox[{"t", "-", "1"}], "]"}], "-", RowBox[{"0.1", " ", RowBox[{"Sin", "[", RowBox[{"t", "-", "1"}], "]"}]}]}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"y", "[", "0", "]"}], "=", "1"}], ";"}], "\[IndentingNewLine]", RowBox[{"ListLinePlot", "[", RowBox[{"Table", "[", RowBox[{ RowBox[{"y", "[", RowBox[{"1", ",", "t"}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", RowBox[{"10", "\[Pi]"}]}], "}"}]}], "]"}], "]"}]}], "Input", CellChangeTimes->{{3.5084737092076063`*^9, 3.5084737097170258`*^9}, { 3.5117245943766537`*^9, 3.511724598787417*^9}}], Cell[BoxData[""], "Input", CellChangeTimes->{{3.5084709397632236`*^9, 3.5084709397894278`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{"$RecursionLimit", " ", "=", "50000"}], ";"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"y", "[", RowBox[{"a_", ",", "t_"}], "]"}], ":=", RowBox[{ RowBox[{"y", "[", RowBox[{"a", ",", RowBox[{"t", "-", "1"}]}], "]"}], "-", RowBox[{"a", " ", RowBox[{"Sin", "[", RowBox[{"t", "-", "1"}], "]"}]}]}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"y", "[", RowBox[{"a_", ",", "0"}], "]"}], ":=", "1"}], ";"}]}]}]], "Input", CellChangeTimes->{{3.5084702833515253`*^9, 3.508470371260749*^9}, { 3.5084704049457808`*^9, 3.508470585633668*^9}, {3.508470630221963*^9, 3.508470685659617*^9}, {3.508471934868383*^9, 3.5084719502552233`*^9}, { 3.508471987707246*^9, 3.508471993591114*^9}, {3.508472071185714*^9, 3.50847208369586*^9}, {3.508472114869843*^9, 3.50847222598955*^9}, { 3.508472274639151*^9, 3.5084722924350443`*^9}, {3.508472818627358*^9, 3.5084728208405523`*^9}, {3.50847307799758*^9, 3.5084730788812113`*^9}, { 3.508473124957897*^9, 3.508473310091757*^9}, {3.508473348585227*^9, 3.508473379803138*^9}, {3.508473410322609*^9, 3.5084734234021263`*^9}, { 3.511724621571731*^9, 3.511724728636238*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"y", "[", RowBox[{".1", ",", "10"}], "]"}]], "Input", CellChangeTimes->{{3.5117247694504633`*^9, 3.511724786961331*^9}}], Cell[BoxData["0.8044790517892619`"], "Output", CellChangeTimes->{3.511724787377471*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Table", "[", RowBox[{ RowBox[{"y", "[", RowBox[{".1", ",", "i"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "50"}], "}"}]}], "]"}]], "Input", CellChangeTimes->{{3.511724794269846*^9, 3.5117248569672422`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ "1.`", ",", "0.9158529015192103`", ",", "0.8249231588366421`", ",", "0.8108111580306554`", ",", "0.8864914075614483`", ",", "0.9823838350277622`", ",", "1.0103253848476548`", ",", "0.9446267249757758`", ",", "0.8456909003134376`", ",", "0.8044790517892619`", ",", "0.8588811628781989`", ",", "0.9588801835332693`", ",", "1.0125374753333127`", ",", "0.9705207716506487`", ",", "0.8714600360811616`", ",", "0.80643125206545`", ",", "0.8352215837319565`", ",", "0.9313613329199122`", ",", "1.0064600575970797`", ",", "0.9914723366307845`", ",", "0.9001778115580217`", ",", "0.8165122477044161`", ",", "0.8173973786334565`", ",", "0.9020194190509736`", ",", "0.992577255251636`", ",", "1.0058124302614133`", ",", "0.9295565852134531`", ",", "0.8339189923730028`", ",", "0.8068284135422159`", ",", "0.8731918019635126`", ",", "0.9719949643727989`", ",", "1.0123987289051053`", ",", "0.9572560607809363`", ",", "0.8572648747702095`", ",", "0.8043566061582071`", ",", "0.8471748731078222`", ",", "0.9463527584521338`", ",", "1.0107065717878339`", ",", "0.9810697139168953`", ",", "0.8846901752884866`", ",", "0.8101788592405517`", ",", "0.8260411261210225`", ",", "0.917693280912586`", ",", "1.0008707551754459`", ",", "0.9991005626649045`", ",", "0.9140102102114926`", ",", "0.8238313754466118`", ",", "0.8114740631720894`", ",", "0.888299529304456`", ",", "0.9836747945804032`"}], "}"}]], "Output", CellChangeTimes->{{3.5117248282080812`*^9, 3.5117248579924517`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ListPlot", "[", RowBox[{"%6", ",", RowBox[{"Joined", "\[Rule]", " ", "True"}]}], "]"}]], "Input", CellChangeTimes->{{3.5117248823715754`*^9, 3.511724890714814*^9}}], Cell[BoxData[ GraphicsBox[ {Hue[0.67, 0.6, 0.6], LineBox[CompressedData[" 1:eJxTTMoPSmViYGAwAmIQDQEf7FFpBgc/ifp9q/zeQvkcDiLBcx2Opr6C8gUc ItTv22p9egnlizhkaZVaKce9gfIlHDrZiiw35L+H8mUclK8kf/bSgpmv4FAR FXY50eodlK/kEHT/dcwz/tdQvoqDioDFLa99MPPVHNxOm0V+qYbJazgI3lpY p7oepl/L4VXdi7ooY5j5Og7HG05NbuSD2a/n0LsqXIPhMUy/gUPZ30MvPc7B zDd02BaROFdpF8x/Rg67zph4bDsL87+xQ8XBxjulUjDzTRwsfk4zVNkFM9/U 4Z398k2O52D+N3PYO0HU+54CzDxzhzxlpR45DRjfwqHjpOqV8Lsw9ZYOyk8d EwwPw8yzcvilEHfsnDjMPmsHk30rit/ug7nHxmGa9o+HZeth5tk6RGw3udR5 EeYfO4f3hkLv9T/C/Gvv0F/8W2+qFMx8B4cnPUVMJ42g5jc4OHCXxZVdXwgL T0cH1p0q+7flQfU3ODqUrtkh4LUXZr6Tw/Jt+6v4ZGDyTg55DI1TWj1g+p0d Kp5WLb+lDTPf2YHpAnfAmxSY/S4OcnPsjib6Q/3f4OLAoslt/ucNzHxXhzvP d515mwf1X4OrQ8nLFT77YmH+d3OYyjqBdyIzzHw3h6zPjvsW/oCZ7+6wzUAp fJI9VH2Du8ON4/fmXo6BhZeHgwZ3TMmM71D7Gjwc/NxylT7nwOLD08HrlfAE 56r39gCxcOg7 "]]}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesOrigin->{0, Automatic}, PlotRange->Automatic, PlotRangeClipping->True]], "Output", CellChangeTimes->{3.511724891023458*^9}, ImageCache->GraphicsData["CompressedBitmap", "\<\ eJztfQfUXlW17U8qpIdUEgIRQ5MrKnIBL3qxo6BSVHoLiHQpgoh0BUORjrRA qNJ7EQSkSOigBEgghN4CgRDSSEj73z/PmWeeNY/7Oyf33jHec7z3xrjh+p1/ nr32XnvvtVfd54ejD95z9/1GH7zXbqNHbnzQ6AP23Gu3X4787v4HdTzqvExb 2zJ7trW1vT6yDf+7veN/8j/tHf/Xvhf+gyfp//3T7P/16tS+ePZ3Xrj70qcP X2fgFSttf/p/3PDwh20L29u7tC98//d3HLbtwfuefM7ZJ+2z7fcOv3Manv97 9lrvjv/OOffpPT87YdDaGz6w0UYbdfz72iPrf2V0n0dmtny1bV29++lP3772 wCu23P7w0WeNG3vmiftttfN+l7zcNq812bb1snd7dm5fNGOnJy49Yuq3hg49 8pu/vbXrewuWht6StWdPOOPmXTbbe/cTL7jgrBP33XHTI26f2raghl6nnEtb 5WQ7/jv78id2WG30+uPeyv7c8XS9WY/v9tia64/t8ubC9pUXvH7uxPU+t0ev p2YXXe01eOH7t5938B7jlp0yL8P3XjTrpR5vzV3c6r2yt/O2fu3qs4fc+d4C PFnl01dOe+OLK23S774ZrUmun73bo1v7p688fup6X17tuMnz67pYklqw45uX 7bvhb8bPaFvSnkGmrDtq674PfrykJalGxsw/YvLv1h636uin2ua0t/dqn/3E lB1GfWHUCVPms5c9V57/+vWnHnHgGYMenLYwe9L8Tg/M4Vc/fujnE1f73EE9 JsxtiefIlmlfMvPoJ47Z9aiDjh1z+JjjDnlz0xHDT/ruIaetcMu7s1qS2iB7 d7n+7R/d+/bXh39p5XNfn90S28iFuWc/s8/qp6xx0LNtn7S392ifO+HF/VZf 4/MHPDuXZHpgeTy031e3/OIf/jGrg/kt8fmG64VNs8XbN56z53nnHnPHNpvs P+LGNz+tI8PXPjv/5Yt3++WFU7rNs5XR9FqfxTMnrnnPpJlti9sXbzTjgd1f GrXq3h2z3/K1pVoUn69yctUxLxVcxwQv3vz9m77x9IrrnNH11QUt8eV6mP/r yWPWOeULBz3TNrd18yX35u/w2uWHfe7iiRMfP0VsaHptxYVvXXPMdkefuc/Y s4//40F7H7nWFc/P7mBJy9f+iQ2zrnxs21WrQuPRNTe4ENtqpYVvnKdtVcrV Jd+afu+uv/rRiZMhD9eb/fh2z63yb8d2f3F+63fL7bz4G9Pv22/r3a5+p4Pe 4q2nXr3BAyM2un2ZD5a0fpdytWv7grd+d8MRRzy7x+dNetT0tfPy7dMfOGS3 kyahn22ds2fdWr3T6Z8Onc0n3nTmRYftuv3G+901rW1Rzh1I4+P+fPgOBxxw 2tjTT91/h42PvCuTxpSoyy+a/sDFx++y+8+OHzt23O9Hb7bl/je8Wwjx1Guc xUELp911+ZE77NMh+8+68dAttt7guEc/7pjFlq+trw16Mya2YEcDld7tsx6b tO2aXxi/wYYbjt9ww6+e/v39H15mxpKWr3WqP57j/x6W87stP9vb83/47w+W 7u8D8P8yzWCrDJN8ocf/Y896/gv15f/Us66t18+/Ujf//3L4v3k5dFq8eMm/ HCv+dyyH/MjsfOYZpz+sZ13+hXjwX1wi/1k+W3HatDnVZ8JtlHg2In92//6/ uE3vrlRP47/QpTe2+smfvvrshKnZr1H5s0OuuvIZrbwt82ff+Pvf39G7xL2+ 1U//9O2nnnw7+7VHa1Jfy5/1njt3wZKv/+f5O9z9l5faVs2fXXb8cfet/fLL 0yOpmw7/zV/07mfzZ10XLly86BsbXbDfDdc/17ZnI6n1J018v32jr51/xGWX Pl2Q+uuB+9+27403PJ/92iJ/1mXRosUzN/nexQM//nheMarPvvP2TLx7zqmn PLQUpHa94/YXF3zrG2MvPOmEBwtSk3fY7uprjjnqnkhqg4nPZ136wstTPixI fefJJ9+a/+1vXoi5Ld7t9cknC0Tj663JD8mfjZw6dVaVWaf88exHH/jFfrdi xEWzczb+zkXvbrn55bFLv/7TFf9Al77/2KNvFl3a45ZbJt39ywPveG/zzS4r 3v3VlX8q10PsUtYjduSK4377V7AzduTOQw7+8y+vvnrCK9tufRUJdOk3e/an s7/33XFv/vTHf1r1rbc+Zl/+cvBBGc3dbr/9hQI65vzzHscMzvr+xuP6z5o1 n1288YjD/7LZ+Idey36tlT/TEsa/5VtzbXD+bLU33/wYy2no9OlzC2rYCWu9 9upHmJDOHaPNx9AFj17cYftr/vS7Y/+Kud48f9xtwYLF6NfZp586/uiLxz3F Rq4+9uh7t7vnnimP7/nzmzZ87tmp7PKje+15UzbJ+PW5/FmfOXM+fXb0ztdl vwY0dvmSMcffj8n60ksvfUD+dkUTn3z32xehu2/9ZMsrVnrvvdn8W5fvPvHE W/ceeMDte95y80S8y36DUU/t/rMb9r7pxonn/+Hkv7Hf6O9/PPfce0D+7Pbb XmC/MXGYLDCh6PdxF459YnEH70C0rt+D8pZXeeedme9v/qPLHttzj5vQJ3b9 K88//94/dtv1evx6aN99bvnPZ555t+j66D/f8eKlvz/+/s+99tpHL2+7zVXs +pGXXvL0ieed+xjm/rZfH3onu/7hDze9FPN4aMdq/sM5Zz/KroMvx4+94IkL Tj7pwbY1c+iI99+f/cGPfnAptjxeae792JNOfPCoSy5+GnJqlzv/PJm9xyrF gsevy4/73X073XXn5GJS0M3fX3D+48ssWdIOYsM++GBu22Z5c1gD2G3rvvjC tKd33+0GvgLJOLdjg+LXDx8e//odhx5yZ5xjbF6IqOk/3PSSTFStkTeHXh3T sfiwjL445aUPa8YzMG9t5femzp622Q8vxS6EmICoLyidevZZj0Ak4NfvLhz7 JMZd/O3cU/7w0D433pgJT8jnrf/615fbfpT/bdlPP12ERYJRrPDhh3Mz0cH3 IO2e33mna/ELu/21rX96ZfE3bMIp221zFX5hljJhvXr+t3UmT/7g7R9veUXP efMW3nPQgbdDUtYtNg4PMhuTjl8QPBhSQQ4S5kfjx7+OX5i/i8f8/v7ib7ce 9uu7Nn/ob5lIOfDaa57942mnPsS/LXvyOX98FN2GNL3oxBMeOPyyS/+OLY4V POSj6Z/ksG44IB/eZ6+bM7nYsTVwJuDUWiX/Mxbt7Yf+6s7iz1iePToGxz9j +Y/fZ+9b8AvL/p0fb3E5RMxqef/uO+AXt2WSseMX1l223mrYkf+p2/APps3B 9sjWTMcTvIbXOWqIDJx2+PXNp59+J5NQ/BtW57+/8MI0/MJqzSQVOYIVj92G 7Yzz4oQOOY3j7ok9dr8RMpEcwRoqJgP/MJ9YfxzyQddcPSHqdxN23eU6TDv/ /JMH7n/l+iOPuLv4810H//KObe+9Z0p+JnXDTKJPnXkuYaYx481MOePMMx4+ 6dxzHiueY2lBUJJu3479hr1YHHg4YDH9HBREGRY5fhXaw/IzZ87P3+6MoxNs yfdn/gwHSvcFCxZlvz6TNwORl61pdgKsw4JiJ846/bTxWIfFn6869ph7t7/7 7in88y+uv+45HDzFnyEzcMRQ9KFfk3ba4ZqsI7kvJdPpIB2zX/GAzHiUP+gM mQgZkx2LfBE7+Llddrou9L47TgoskII+CM77zrcuzM7/HNMZqxd7IJsgthW3 YIFDO1+e/OIHkWNY/dJ1R+bPIAmwuJabP39h0R40gR8/+MCrkbPYZNIOOv7h cMYqzFnXHadXJtj4Z0wzhBDUs0JtwdmI5Zb9IikIfxyO/8S/yiIjJ8FqLLTY BIQFZG8c/e633foCpHjEYbXh1IyjxyqExhRxkJcQTrE97DGccpGbf//Zrjdg Y8X2Nn7i8bce/MW+t8b2sKL2v+665yI3J+6847XSzzv+bXXffa/ccOThd3O5 Qx5CXYjNQEWX8stu4RCD9vn5V16ZXth233v8sTexXZo42j9/NmjGjHkf/WCT S7J1wSawthZ+8+tj47o7/awzHi5OsaJLEGwQcJED0M6gk0TcRs/8413so9h1 7F6ccJGjOpFCe1BWjx130ZOxvd9cftnfobFEjkLw46QsxgeWYLOSo1i+2iJs BiYBNIlMvWO3oGFkJkfOsOy/2EaFTrUUHEWPca7HkeHsLsRbMTJM0g8eefiN 2CUcnsWBVOB+ev99pcQmrtANcNQVOCwGPMt+5RZqZ2jXEHCxPcxOphZWtiGO n4iDggX2RBzkLTiNE7hg2ZM/3/3G4kgrWIb9Az5kv1bOn4HNaC/jA/mC/z11 i80ub2Jtv/zZgJkzs8X6mXffnVUM8ZldR19f3YbYrhkmdB3SCUdnxOEskDhh 13GEfOvpp96OXR/19tszoTFH1mJDHHDdtc9GXDEF2bHE9mDJZxZjoAtBD9Mg 0oVxu+brr8+IOBkfAYdx4bTMBDbp4twVu8laTBE2cXbgNrMWEqpQvYshYoVC nMQhQsOdU5zj7BL0n+rqgaSGgRu7js0FvTq2B5eCWEG6UPU1VcRt+ugjb4iN bA/sAtsiXax2HOsRd/NvDvvLlg8++GrEffrtb46V9tClxOHkyH7RFyS1MPQP q3jwRx99EnZ9HWvBioOvvmpCbAJ2GDZdHCJsSGym2HVsVKzI2HUYw9rUQepD m43tQd5BA490oVdmplvAoX8QqbG9vlGYkC4ORq2yzqWw07tsrzCpIl3oTI/s vefNkbVQW7S8iMPB+G+vvjp9KVk77oQxD8B7EJvAcA4tTip26ee33jJJakbn Uh5hg8Uh/m2/fW75+j/okgusgIYEzatoDzRBO9LFCoPmFOlCR8Bmj6sMuweC IuLQllZe2FXakcTBhnlhx+2viXSh+0uHIWtNJhOH3SPB1szaW35T2m5FE7At MtsvdAnyMzN3Q9ehQEBxzVhG1sLnAVM1sgLsgXoeu46pkx3RNbELiIPnCtOS HUXEQRis8cYbM2L/oKTIm8r+gTVStIkDa7SD2B6UYxx3ke6r22x1ZWHdFTjb kc2sxUaAnhSb2PEvd7105W+PtdmG+QS5F1kGpRrHUcRBl8qOg4DDtEj+suvw TMCYi3TBLvjhIg7T/udfHfLniANrqsciXHvSqMlaKFUQHHDTFDiY9WBR9qtb os8jEjuDdE8768xHZLJF1mac7ZeYd76IOQcDY4ehK8uzS0bhkJPSSQZAKa5K MBhEcI9FHESKpHjX8l2pTcTBHtdkEIetvfOdPFiIw1rGwRTXNJiMLR5xYAgY ExkK4w6OiezXivkzHHBoL+IOu+Lyf5xQSPHmtYqhYEix67BcIKhjl2Csa96J w0FTPU/hT5PdSxzOUqh1mU5BHAzFzIEecGgb7uXMLiEOzor1Jk16P+IgPjDM SBd7JvOSBZzRWKl8V3YOWQYXSuYfDqw1zZU4nAqZ+6uBtX3zZxiuPL7dSlbI mGWXoBdUVSScuXICEbfJo4++Ifua7eGwkKHJbWaHA3HQKGFIxyHiHMe0xP5B 44UXJdLF2S4vfGrbht0CH3Ska4fSiuXy0hIhDoYSwh5LyVp0HR662ATYjdWT bS52CV58zVjXUlrBbopdh/4Aozy2h9WtwNeIcrNKwBAHKzRz9Ichwt6XakYc PHPwAUS6UNdkjqZUuJVKISZViu1dd9SR98B2jHRhreuQIw5uSIQulpK1cC+/ tP22udLdPX+GjZD54wMrsEqqWi9UocKZm3e9E+wBtoIjUB4OtlL4c2OHMQC5 HTgwuAhkWhAHzUL7hoyC5IOnOfYK+6uqXZhVyfbOO+Xkv8GHFemaO5U4OIpg gS4lQ9FNeUzICrgb5L0bUW6VqkQE0yH4Y9exTrFeY5ewpYp4QdF1KDdFdKSg m1ozmEh4RmN7OJCkO5Mupg57J/YPBrv2IXHwhWSe0NAerD1YfZEu9iZCaRGH Aw3ez6VkLSwgndkcIqwOsXFEuX5lZ4fDBjpsZiSy65BAUsjYJTsI2HUoRTps SBcxlCLUU+Cg8ElvZHuFgzcehoUbO+IgsWWEsj1z3JEupDB8FREHVsNjG/sH SV8ENHj417EWglwWFZswAb9iuQKgpkccDqXV36S6xikoYm3V2Zb6wvZw1oMd sT1s/qOLzb9iOc2SiN1LlmnlkW4R5ox04cGVh4ztAYP3Y3tQ3+RlJs4OPuKg X0s9bWYtVHSo6rEJU0tICl2sCg6o7XBCxyFCk9VKIQ4bv2chV4fnz7CypfYT B6+zznvShf4glYs4WHI4SCKuCLFGHByj2NgRB9UKKlbEQZfReU8cWC2NmTjk Y+gsiKzNOEuGjolCny/ClJUKTQJ73XzTRJ3yYdvixI04eF9hv2S/ls2fwQKX mk6GwjSQdsH2cPJfezQPObYH16JsOuLMYU1cjPoWOEhmuTCJwwJQCIA4+Nm0 b4iDsghGRxxOgCwKvXRrFYF8qfhswlxwJIUgCFgecbCI5E8iDh6CLK0lsNYC LWQtNAvF3dge9Ei5CNnejE2/n8URIw7KU3UNwuSVQkYcrEKIkYjDlMJ3F3Hm HSYOrK7uEUyTAkHNrIUuqCgcmzBnMklB9a5KOqyy6hTArwX/VmStqd9kLXz3 MlbZHvQHJdUQ15nmS2bxsT1I4d9edOGTka4FpNkeok3yWREHx4/ibsRh6j7e 5HsXRxyWR3W34ECTOdTMWkgNZd6xCXgdoapHUhD4CvsRh1lVMGvFcsXLwbFs uZGkSpFlEAY6INkeDpYstSPgzIfA9hAPlHAKu0orijhEzhUIY3sQdNp9xEFP yEyhgMNBDQd/7J/RbWatOXDZBA4fDScMUTO7bPmsqkrBtajDkDhohdrAbA/y EvH6SLcIgUQc9A7FDNgeZGB1o1u8mzgctjKc2V6RQhBxZkyTV4VaF3Gmnjaz 1oKIHCI2H4aYeVrZJYRb5WtfLn+Gc1juNuKgFUpVYZewiXRUDk/IMuKgLUK1 i0M0Q4Y4c/eQtdhlWcQ64KA6KswyPCGciIPqqCgZ6So5KuBMTSRr85S/VuWX q82ffNwr66yy2YC7PlyC8p2rjthm572OPueiY885aqdjtz3q7s4fLGpd0tl6 0jDn2eSEzmGNy1XHQZguyUnDWa/BkikwNiHdY3umMA4vpaTc1MTBbIf5HunC sZLlxVX2odyfw8uFphg7+weJrQN8eDm58nawPRjXWJSRLhaAfIjEIbQJD8o/ T1qr0tCu7Qs+/PJzT73aae7i+YdPPu4L562+x99RIdezfc5Tk3db7YujTpgy p1W5aWrS+uTPLPrFwVpCDQcBb5hsK+KwoWVfkik4L+Tw52BxDkg5JA5auRKV iYO/RGm4NYvFEiqIg3tA3kTiICQR6Ix0YQTJyUu6SEfc4m8Pvhbbw5kpTx9x kCpShotJg3RqXWDbSRJs4Q7vXP6tW1f74W2dpi1uH7R42s0vbbLqt1e48l1F z0Nl7IJkpQJnDEtdeiNHauotR2CmInFQTuTXIUegnMiRx5FCRxTnhpWjh5cg tpfinFnhxCGSJEWGOEufIc5C5KnjhDjLGiLOnIrEQUPQKd83zFjrAuJ8Nnp1 bLd3Lrl253XHbPzrP9x16Mbr/niX69+F7GxRu5uUi5wwnKMKdrNjph5xADA+ s+qEgIMqI02RuE7MDc1kLXFfmzBhapH9WeAsWkocPKxVxkGLVUIbcWhLfgHi UHUhTZs4c5ESBzVSnjriYNQrUk2ceXeJM7kdJ6x1iTMnbI15Lxw14UtrH7bc s5+0rz13wq+e/vw6v0Vdb4sq43zCumK6euWkYVhINWJ3kHql8Aq7DYd5VoAS cNgD0IIjDme13CjE4YhTUJs4c0oTh/QR6SbEIaNYTnjiYF5v89d7X4445C7g /YgzLXtYuSyrdKGCVpcb2pPBEZaHjrk4Ta3LmvNpWm7d2U/u9Pjq6xXH2x9f WHeNXXs+Nafdyp27t95M8JwrVWC5UuLIwTGslIgyFYmD90UjDtxXfJw4s1OI g6OhujpSdGEbIck44swmGVYKAHkTiCvMzIhTMn/AIW1OKRnDEjSIgwIFqfjP s9SqbLrYTN3b57961qX7brPvr84ad/KZB2+9476XvdY2v729ZRl26wmD4iCv DTtmUXIOAGIf4j/ikPggnzpxOJpkSRBnR0RQHOS+IM6yt4izfGvi4HJR8gdx MKSlorN4D1Yd+pj9WqGcHBmRy5V9rh6n6JtcfMsltjgnLGUlkbUpwW2eoWEJ oRq6LtuZXYdmoHRO4oqyjIiz3JjlEvuDdLHX5CmtmXpLriBdaAx6Rrrm+ggH n+JsbM9iF8SZu6aZtfCFVoUtdP2q8I61EUXXEVCT0ceuo9uK+hCH4kelbhKH oJ0OHNKFKi4RNaxkj2we4szrtUK5M6CyZmGSnolpIQ66h9wrbM/yZUgXhmZ1 l+Ich3N0KVmLM0tmXTifqhtd+QKBZdhEcn2w6xaKIc4sGOKwKZVwTLqWGEe6 pjMSB1+t/J5sz9IsyFpL6yYOUljeMbZnuQ+ka67xoFeomLWZtUgCUgIYm4A1 K22RXYLMkkeKXccBIu1zhXI4ko09Wg/REuhI18wF4ooa3NgeQt4yqlYohYli 9z0Ty4Y4s7ZJFytdvtphpbCTnCbOXDjNrEWYWp55dh05VioyGZo/M+uMXQdb 5bRg1yFI5Dzo0XqI6GJWWRJwFtEamlhRpGtxTuKsCoU4RPvlZCTd1JkBR4k8 YUPLKZXjsUe5MyCglpK1VnHSIyGjSMoUZXYdfkVtEHbdXLtsD8eaDF/ikF+o NE3izJdEuiYHSRf+KiUbEWfqIHHmZCRdHKl9i9y/Hom+sD1TQ4mzNM3I2oyz ZCi2nTJJ+KKdiSSAbYztHDuMtaZtQVwfFrRFHGLqOtAC46sDM48q24PlL68t 27O0HuJMvyHOomFDSxEiUUMc1EntG+Kw3xT36FnuJfmh+rReq71LUuI/m7DM bJJCXFEJPsSZPz6sLT0jztQr4iw/gDgEsKtTan72mv6Zr4E4y/QemtjKxMEc UTUPcdA9VA1GHA40relm1uIsVi4Rmzgh5tyTFELXkFgRZ9KPOAsgEmfSjzhI cLkriIN1pmzXGslptYDEWSEScXaOE4cMMe004qxWgjgERpVsT5x5WZtZi4ML oeXYhDl0SQqmvpJliDMBT5w5ZYMkVlh5aGu6ViIyNLHiiYO01jQTZ0Z9z3Ll wTaLOKtPIA7ty0AmDga4NFfi8J4qaJpZa7nXbAL2vzY/SVlCOj1OZlvxEgys duUssz2Tl2wPapkSH4hDKFvem6HlCpWcJs4yAIgz85s4y04dmpiWnollQxwO vTGVcSCYid27lKyFDFX+CJuwrE6SMvlG1kIxl21F1sK6rso8O9vZHpx7yoqt kY2mUxBnTjXiTJUiDjExCR3iLLGTOPORDi1XspZSz8S7zayFaZvdUBCaQMFN lRWQWVLsyVrLUSFrreysRiZD6FRVJDuyiMO0KzWHOER4pA4RZyn7xNkRPbRc edJwibN8VuKsyiEIInlFmllrjgY2YWViJIXTUrYLWWthYLLWUthrZDI0zez+ i4CzcApxZr8RF2+bKHCqPgj9g5KlwNuQcokoBFSzS6E1yHQnDkebQjbNrLWs InYJ3VZxI7sE5VwOGOLMF0qcBaeJM+ufOMv8YtctE4pDNK8D2zPjhhfxoH2Z 3xwbbEnllZOu6bCka5cpDG2Ng19MZ0YzayHv4Fpv2fUhiW72SrBncLlCdXcD aZhLr2b1FFWIcYjmKyNd80WRLtqXvkq6iL3L1Cbd7D6Ows4LAkYXXJGu2WXE xVtnjLUZZ0kUvlBlp7LDOBi0lYeUW0rZKsQh6qQMWA7MMu17l+3JSGZ7FkNn exYvJw66hwrQiLPLMEjXnI+kCyejRFdqHGQUaMpoIUORpSgfAnFWzt+8Vk0i 9iqli2wSdt08UHwXKpJKLomzfFHiLCrAIdo5Trp2zwjbs2gE20P6mioLiDPV kTgcpDpwg7jQfiBuTjRbw56Td613Ylp6N7LW3MocItQe+U1JytQXvmvKNHHI DlFuM3Hw/itmSdZajlVKz2B75hJhe6YxEwd/q+KxxFnInTiY5Apf906IGuKQ 4qZYbsCJV82sTW04C/OSFA4uZT7zXVNViLNAKHEWP6pZjaYDEGehaeJM6yUO 9TaqVSHObCviLBqREljEWZ03cWaM1LCWbLRSGTZhYUeSgsySqUicxcuIw0ZV Gi9x5osizuphwkGq7ADiLDmDOPMIDC5Xsmwm4sy2Ig4+4urqNoFVMwXmqm9m rSUcpuyywYnhEIdwuy6TIQ7asewU4ixhiDgkU6pQgji7hoQ4DAVDijhMs3xR xMEsVuk+cShRyq7oDDiol9WpL25kjDjzMBCH3BpVozWz1tIC2QQUZ63QwYlN SBzytZQ+SJxpx8RZELzm2LFrSIizrL3eCeE0uBy2UorCElGUlzgIOxkoxIHV 1WMR54gKQYnTZYJLx1ok71U3uln6JGXZpMTBISNlhDj4upSyS5xlkxKHY03D Jg4rQhZ8DV3TPomzugjizPFDnB3RGS4rNk25lNgKmK5ip2aGmqHGJnBeQ5uL HbFwB6MRln2fX4rprmu2h4Ccqg3CdpQU753YemFtiS5xpv4RB6NU1RDEWXre 4LJ/UicDyyQaiDOvHnF2iVVkbeeQbIZVIWWJjLJ7esgouNoqLrPOcKqoyI04 ME5lk+yIxRvYYTvGiLMrR4gzh0kIK0mXJF3kOqiWke2Zacz27OQnDm7SqpAw Ja13uZsrx1jdWjXPI7uOGL9K1th1i7YSZ+WLxIGFKuFklyyfLrBM5XM1Qseq ZUjXbrQYVPZZlbpsD0VC0JVjexaFrhES8LGpzC61RJpZa55Hdt2sMXbd9Ms+ iRU6KLGS2aUifh+7npKD8KlVp8D8Z6RrziDSNeuJ7SFfAHeLZ8XIg8pnKjlh e5YwSLoWpA9TJfHYzFrLYiEpK4wclNhIfUppJWlKnEld4qxCkDiz5QMr5P7j EM0aY3s4YnWFVH5/rpsbLAyxvAHirFCK7cHqglUVWWtGdxCU8mf0amStsZFN WDIQWWH+KeKQLSi/GLtu/ikO0ZJ3BiWECduzeinirJ6WOPOLka4Zv6Rr/ini TLFke1jtKlAeXE6L0kNSgrKZtVj42tRswlLdOERLwCHOdL+B5QpQTlXf1gLG /P5sz8pxiEtE9z1cRLpm/JKuZS4QpzuOA868awNdiGVJSKlxNLM2tfBNHSKp zrzrpexSpvuRAea+Y0fMNRk6vKjSYbt8hDjLQyAOEV2dCsQhuC33cw1dUyyJ Mx9zEBcKJhFnEYVmhqaWO8w6lTSQlJ3tfVsz3uwo4qxEkjjzHRFnOgBx5k0g DrqMMl6Jw5CLe8AjXZmoxJk/jjgUNmefj6isfekofcs9IvFYw9oQP8Jyz85O NmEJwSRl1SLEwa5SLgFx5r4jzjJWgqSTOUqcXT5CHKp0JRGJszSSgYm+EAfP rXxqxME5r1t5a3ZLvJM93tCgDN9m1to5ySbMYiEpC/X1TawU4qChyn1HnFWK BUmnEF4NDrqDlHPi7OIH4izdhDi7Coc4c9QQZ7E74hKZsZ4bFlmbcTZE/KXy 8kW4cWWxkIC5rvsmpBpx5rSrERfGFOKwnnULe83ALFFvYGIyiLMMROJM0yHO roUgzorfiTPXU/NatTwVNmFqMElZgkXfclokhYizVKEgLlTCMTDBxpo1aKUZ xFlSKXGmmRBn2V/EWaoQcfDVyssVxIWCTsRZSUgza83PyfsKLUKff5rAnRQk ZQdQzaqwi2iIgypRXWXmMh+Y2HrEWSyWONOnicPdAxobcWZaEGdZYjVLBO52 iaRm1toiT10FSdaCDbqenaTM2qlZFSadibODijjLuCLOSkKIswR+4swKDKzQ AUQc2pdBTJz5iYmDK0Z6MnEI1iiDtpm1tirCBaZyp5C15tknzmz0Aa1XhekU 7LqpV8SZ9UQc1CN57Gt2ldn3KVawPdNciTM/MXFIvtONG6n+kbVLf43DqvNe Ovv2HTfdZ/eTLzj74iN22W/UuCmf1L2bKtjjpNn9W2SKOW/IFKwXOcWJM58V cdA5lSqfurSSTIFdJd8scZboRJxpAzX71RYa2zMdm+1ZUlPNfgUPdMF1av/b pC3FNQ7tn/3k5au/fNFT2edyvzLzkZ3vW+sbl3Z5e2HLd2smzYIwZIrlRJAp dncMcYjR6c5X4swhycFCu5eHjEwxo5vtWTCf7Zl2TJxZuMSZR5B0zeYkzpL8 2B78IMrLGVjS0CIgDtEdSZaenLT6axzy2/KWG7p46jWvfWXkd/vf+1FrbM0s Qa/T0dCv3EaqA+XooERpqfVLcJA4uzKUOHO2pbga5HCV++bN75eQfcSZc5E4 K0riN2zgllbSUviwhXIhOUtWRUkcNAelYcRZary6oc/imU+ddt4B25228wHH /PyM889+46srrTf8srdntXqvZsLMO8qO2WFFhiSuPXSnNRkCOaOYFD+iYi68 5cuttqiofyTOMrtI11L0SBcqnG5lYXsW1Oxfvqu4AXGYfIXs2Z55g0nXDH62 Z9mP1W2Vvrohvx0FXud5v5942JpXj9z6zk4fLm4fNX/K6VPXHvqVwTe+/0nL d8v7G1jQB/8dDuXYd4S6lIzLvpsLtl9CfJIXFpDlGC0xhThzTfUvN7rctzWb FQk8UofZnnnB2R5SIbQxiTMLvX9io4f+yYNAHIxJFYrGuWp9f0OYq+MmHb72 0V885kXcCNCBOf+pr3z5qGVfmNf63ZodlpoNuFtkz5F7drkUceYr4GjN1cPR 2uXnxJm92b/1KrAoMnFmHxJnpS7EmUZDHCSCosP9E9Kjpn+oKVHE0mat6T6H 5RdNv3/s8T87aPQp48adesqBu33t+L991Lao9Xs1E5aS/Zb3UCNy7EJ+4ixb p3/rbQHRp0J/4hD2V3lYzbY1XwlxZsUSZ6YbcXCjKn+DOJhUMK0izi4iIA75 eJUKkaS5xFpJuwKSTVgkl6Qs+4c4SzMhznwMYe0qS6gGZ8Y2cWadEmeuLuIs XZE4czgQhwuBKt9n8ATYmv7Bd1gpd6hjreUksAlz/dRsY/NGE2dOgwyXRV0G tPprW+rzf56nRJz5oMMalCu/f2vhh2fyhtQw3hxhxMF3WKnBq2OoJU2zCat9 r9lmFjkhzlKiiDNnFnF2QxFxMHLk5ybOsoQCTqZ4WIPyzRFnLl3icDaIjTVT YDV4xEG9r1SO1rHW/HBsAjngMwvRULPNLN5HnLkkiYMLVt/gIc7u6CLOysGJ wxGpBD3iLHkrtKfbdWvEFKZJylL/xLSE1S12hylVKWMza61wlc2aw6d/6w1n 7ky+a+lZxKFAVd8mWT6xCYmDzi/HVVg91dVomgBxVnJTI6YgXGS/9E9sfuIs 3hemvlKlX8daiz2wWVNH2azFzvjMMir4LpQdlQASZ1edEmeqNnFmY9asbvRN nnHizL1UI7BwEFY+3+axG+LsQCMOHonKjRIhUEWGWnSMjZmm2D+xzWrWtCVl hTWoY5HvIhWwqtyk2jNtmc/MX8d3Ex9jS4opuJx05y9xlppIHNak7GDi4Pur 3H5St1atEJnNmouzxiAwVyNxlpRVswbtSl8+s49b1OihVixOnCVl8Zml0PBd fdwj4CzlhTj7uBhxxYc2lpK15npns5Z5yGZT3hFzCBJnF7LwGTwSuh0qxcbA 7upwUjaZVfwQZ583SokpPusRPy3D9izgllo2fIblUbmopo619l06NmEFy/1L nDyFxJnjgjhLyiLOrsftn9j8KSFBHCSiHP3Lt16htpKJ06dc8ItWJfzPur90 +YSAYXv2cTFq4BBqWjbNrE35VxKf43MHKXHmb+tXbholsQ5IrG7iLJ+XuNSB Zu74lFztl2BPjfPX/HzEWXYocXYvG3E4gCXsmlkLW7zadSs34RAtdEScFVOF LsmdlYrPpaIBxKWUeMtqXz4hxNieuflSISbiIFxk8ITVqOxQ4iwPhTiojZX7 1upYa5mH7LoVNXGIcxKueUunYpcslDqgXN0KIhFn0THioNJJCyHOgjTEQfmv hlwtC71mSnFeSMMdkBBsxJnWQBxcSrqXtYa1dCbbB33YhN3pRFKp2JTVWASW yY4izqqKiDMTlThY9DIpibPICnGWKECc2WXEWaIAcdByZJeFqdel4sSZWU0c 3D9aSs2sxUaQl4BNWIITSVlAnDgzMoizRK2ahALLcSAOGfeVL3546SFx5nUg zr4hRZyltxAHo0P3MYfYuYJ+xJnNSRx8SlpKzaxNbUzzcpGU6X4DEuwhDrJX NcYhDUYhFeLMWUgc8hvktCPOQhYDypUnbwJxYKGSvIgzpztxMH3k7yLOSrKI s4tjQ6aUEhIjazuHmBK0VnmW+KL5AUjAnNMDynflmuhXEpVHizj77hhx5rIh zr4oSZy5gAYktmhLupm7Mpg5yvRhK2amh8inNCLikO0sj0rzCrXvxrKJPlGt 7pdg3oAEk/uVq1FJWcRZXRlxkG0qziLOrPCaTWifgqiha65w4jARsv5T3mHi 7P4V4pBBLQHTzFrzPDKwbsGaVCYjSVkyE7tk9cShS4pGEGcOE+IQzJeyRJwF f4gzWVtD15zixGFK9ImCmj1nlfHEwZGk47iZtZApkjNkrQUuU/m3JGUpR+yS nfzEWVJ3Tc5HJ9rjWaVGKhBKXCr3IkUXx5g0DuJgySmUTpydLcRZyJ84HO9S IptZawYsWWtuErLWlCqSMuWrJunJ4hJhiDJWSRdmhVJnSNc8t2zPPifL9iyx ijhTvoiDJaePcAaWyQtHHKaumkkH60zWWDNrTc/jEOE808WRobKiWuFg2TzE pZIALWTJrqfqdyzbiO2latzMpOlb9lnZUKRrTnHSNS8XcfAzzKlMFWjI4iNd pBdkO6rkXx1r7bugbAIpzNWkzlRdjuXdEGd+XeLM7UKclYiHKa0mp9rlcsRB CMm/2jex08iyVKkV2CU2DkpMVZ+S3XMrOLSlXdXMWvguVYU5MCEk2CW7C5Q4 K2/qW5KX+2Ng4lmqQpJdh/4qxyXbs6/+hSxIyUvi7Hwga809Qxx8Iwrmp0q/ U599CqXpSh1rZq1tanbdru6pqXxMpV6bdhE2erXo1b4Rza5DAapOqa1Qtgc7 Xqd833KzKheIOHMqki6WgoL5NfXTFlYKFyqIB5G1GWf5mSy7eJov2oU97HCq Sjdxz6xfhJa62SHUoMumIw46gOyeGnFhYqpm7ZtTnDhowoqDEWdVncSZaUwc ggCVD/Ik1ypZazkubMJu4uQQLXU0VdwSasurazBVGG75NsTZx/+Is9tPiLNb o4mzL0wQl6qlh29L9UrE2VegibOPkROH3zIemlmbWj1WLkVSqRsRrN6GOMuZ CayoXmcAU1E5B8Rhmqr3oJhndFBCTBEH4wEO+YhL3QABD4NKHmvaw7iUghJw 8jo0sxa+KBW9sgkEj+XbIinz+xNntb59ynerG86+SUac2ffEQbWqXK7qCUfE mTOeOHPUEAfPq+w84mBzStgRZ3Zen1I4Kfge6B7j3yesY635p1LV7SRlyjRx VutLnKVLEGc3kRFnuYPEIfIlNhIH73C1PTjjdc1y2OiaFuLMKU4cfFuS52G3 yF2Zus4gCEoth2bWmvuOTVgWcupDhjWCw5J8iDOnDHEWgxhUbv7q5RV2KRtx lnIaNrpciWG3qJqJOGAknIjDkSjPAXFWphXGKyHWzFrbrGzCMmNrjh1zihNn qWnEQd5psxKHOIUu2SQOGq8c6jWXdVjlPnHmAB+UeEYcVnaVZZZDRpyFd4Kg VOJUM2vtlOeVJqYZpu6iIyk7AoOMkhOFOLsTnTjLoyEOdprqwGqumLEdRNyY +I2LQYmVTBz8XQq5E2fCjjj77EYQlNJwm1lrl0SStVg5clykblAkKbNTalhh d6ITZ/6kcEuNkiRI1xJKalhhcpW4xJ3jmd0oM5g4S/Ygzr4sFASl6Daz1m7M 4RDNCk/d+0lSpjWwS2ZGhF0gtZ/t2SVD4W6lyg227tkgXUtlIl0LkRNnWkO4 2UzLhnTNcxDucJUvjziYJRpvKxvMbCu+aPn5qTtqw/1eyhwJ20yGLtszZxzb M2M63AMm5qUYH+hKsyZd+3YbcabhEoeAukTc4JLJKgMjQ60OgjdBgsH6pN6y jWt1ZkwQGlyKBgn4mpuVzS7rXW4zuWcGl+utepm31dwMKbe8zMfeiXUUxIW8 Tb0TYoWsNbss3Bipg3lI4l32D8WympbwLYHKd8WTrO2eWO5DEmxkl1L3gduK 712uAGXZEJe4H8kv7iZdmI+KCvQqRY38ukNa07U708ha7Ap5E3qXdGV+sz0L zJCuGTzEwbumLJtm1qYWfuqeTvMskWXm+epdbrgelY1uB1rNqgBd5SEQZ9ku Q1rTtS9RDC43sCQ720MUQ0YQ20PwvDql8EMrfYU41LdKDWtmLSz96hBhtuqg 6l3OmKrC2fXUBxGwefXtR7ZnGilxVkQSbrWV2Rru/dXNH0Na0y1uO4s4029C e9LAibO0QOIs8hu+bKFAVDNrTdlnE7iip3rVOnxOkqtkrX1zukb8mx1FnJnB xCEILj9bDc6iG8RZeh5xpt8QB51FdmOYUhk3xFncb0i5++QXI2v/x1dlLFl7 9oQzbt5ls713P/GCC846cd8dNz3i9qltC9K1kN3LLaXr44YkGMVB2DdOibN7 xYjDnpHxS5xd9k6cOUKHlAyt4jDhSi2poWuJf0E0SA8iDmZJla7dG0WclcUQ B8+hRKFN2v/gqowFO7552b4b/mb8jLYl7dm7U9YdtXXfBz9eUjNp9mUAds6O HQ4CR50MGeKsJoM4uwGNOLt5nzjzswWcBBtx9k3BIYldQJwF5IeUkqC6I81H R1zi9qnkF5LgulCyVzFpS3NVRudwVUY7q5A792//6N63vz78Syuf+/qCZE14 9wSX2RNzMteIe/Nl9UqMljj78gdxljVIHDwAylAhzspCiLPUD+KsEow49E1e BuLg6pNpQpx9npg4RHVVuxBktVyMcYb+u9dkLOp4PnHNeybN7Nhqizea8cDu L41adW/dPpNduNA1J205VeyOhW3Ybfj95YAhzj4u0CvxjDireyTOsliIw8kt vZY4mAzSnYmz5Cfi7Dq6MJ0SEuE7VwrvDCmFhLLMibOLKIiD10E8qG6k/+7l GItXXPjWNcdsd/SZ+4w9+/g/HrT3kWtd8fxsSMi6zWUfeWLvLP7HUWCRKhpS IzascJcfebIMC94bYJcwsD0Ej3Uuhq96SRUmztK6iLO0yXBWKh2SOPNI1qwq ywkYUgoYeWXjrDVdk4HknJbXZCzfPv2BQ3Y7aVLbvPb27q0nyuwTMtb0/56l 1JJCyo5bBnw4V5RTGj6rJvOtZ2ICiIMpKGcHcdjm1U9sWoom6dodLMQhBakq 3eB4UWpSaryserC7cYiDo0T+XJuo/+7NGL3bZz02ads1vzB+gw03HL/hhl89 /fv7P7zMjCV12qAl8bBjZlWQcVihqmwlQ0zzI0Ps+6I1n2m1Wqeh5TOVr/Us d2X1M3BQyysX+rmFQxw8WHKXEIdjsTpe+xgBJ8xMAuLgD1ASdbN1ZGVpbMKS bjhEqw1l180SYtetuJjtwbNV/RSh2eBhCuQR6pEQKKnPt4db/qtHSHGHXLae w2fllDzE9ixrgHRTn6+F4JXjuoa1VNTs5m9+yB1KlURP2O7Vrzeao45dN68O u4STU5oCu158ITL7Rboo46jUP7milvoaOunifFBMIOw+ecU4DrhzlLQ0NNFn 0oVqLhG6QuLd7o2stTpQNmEOPZIC6+WzJy71PUXLXQofDVZ4ukdi5bE9hM4r 9U9uuRBntWKkazePDS37Ir9tj8TKY3tWwkechWBWSKz4ZtbaNw7YhEXCe5RT oGTGQKq6KqzMmEOElazoFtuzRHm2B3mpHBHirDCDOMvqI10rUQ6CSCGdHqUQ k7wMdCWnibPvRRMHOS33SzNrLYrDJizth6QgA+UMJM4+tsIhWl1p+Ka9VkDY cEqqYHvQV46prG6wWreDEWd1uaQLT750B9K1T4CzPfhPpTsEuqrZIs7sd+LA J+kszaw13xCbgP+0Sgpee+XkEwc3gcLOxFmZMXHYqLq2gzgrCVqhXClyBoYr J+RmCENUDUI40WVthXEo7SdMqZyaxFlomzgzvInD7pHeHlmbcZYMtZu/+aJl n5AAPKlKMyMOnageSvZRIeLMxOtRDqLaHra8TLweJQ0FoImzmw+Ig74tvZk4 eDhVfBcmUj584uzSEeKQoSKpSxx4ovBS01rNysOGJVYtCUCh17ETNs2Yygow zyJxliROnCWNrVAKIq0e4iwzIUyk3DHEwRKUTU8cnHwSbGEiJeOJQ62ByjaI s6QP4tBfuW2aN799nJVNWGCDpHCMaTjEYS7laA43cMn/R7pQUHS1bSheUuyL 7VnMiO3Z57SJs7tliMNv5SkShyOsUvSTORCq8tI+nU2c6bXE4ayRz7KZtZaM wCbs66MkhViQKprDbFdKwbykmHSthp+stYgt20MCuzxgbM/MbOKswJI4rESJ JOIwnSqi71FOi05+4uxKeOLs69TEQUNSUKmZtSb+2QSMQV3MQVKop1Jsnzgo LVK+epQrWRuTdLHZdIqStZZ7x/ZwmlcFh+nJxFkRJ3FIR1cYkHRhsUkzIV28 p9gScZZPQ5wZ06QLwaSYVjNrrZKETVhpGbuOFBgZv8TBZEAXMlU9KEY60UnX vH3LlexWChJxUJbERqZ/wKhVOgRxpteSrsXnh5VjU3kd6eLgk4VGnNWKka7V bBEHU0UBmmbWWvHNsMSMkRSGLAfk8HI4mhZ23cp92R7cM/p+A3F2IQhxUNtl bpCuHYvEpc4CLAWtPOJwFFUKKbJdWp1SCD/JVdK1uxeIg0mjyEoNa7uWK0+l aWSZFe6QlIXPhpebUBkp7LrpfsMS7CbOEoqIQ/KQbq4nXfvmEnF2FrA9ixiw f3b0EmderjAF0gZIN3Ucg4bOgmbWmqNseEJIkJRFF4nD0SHLhl23QAnpQhOW dsz2OsXbvNkejkoVUiyb2PxszwI+pGvaLNszhZE4pNXJ5c72THMlXfuOFXHw mklI/JN5QIaae4wdsdJIEkDIVKo7cVBfFEhkh63Ilx2B+q2YyLKtJxKqvHQU 4qAOKnjH9iwWQ7rQZqX+DS8nUocccdATVDI6LLG/SNcibWzPzOVujWvVIgxs wgp6l02wcXi5jsQydh3ap7yqpIvMJG0ptmfZU8PLIcpRvmxivRFn2eHEQZuV 0UIc0nak/oWTQh5e4uxaEeLskzfEQf2TeGxmrTlH2IQlnJGU+X+IM6cMWYtu qzSSdOHrUj4e2zO7J2uvE5wBQTGSKUtqkJ/V09t0WOIgfxVEJw7vVc8HU9yI AzNlghAH00X5Zc0MNafY8LJZebRICmaEogLEQa3WvBMHJ55k3vByfcgVTpx9 2Io4aAwqMSIOe0H7hji4pKtHjOmwxEF+KjebOHjg5NEaXoopeQmIM/OFOBxt Eo/NrDWnGJswW56k4JFSRvaK+TNYK0rVZaAHw9PWI86y7ru3lmCQkro1iHQh PpTgSJwlZRJnOizpmsJIujh6FRFme+ZxY3t2UhAHx42muZm1dg9K2CDV09uM S3YdIkDSlF2HqFCgnji7k7B7awkGTUCZLKQL8SGXOXEQUXLuEWc6LOlCJ1Yq MenCt1U9hE1PDlMqP0XQfiQem1lriYvhUJJywy5B11UIekT+DJtI5zhxkFnS GzlEHI7aSMTZDW1BNqo0h0OEIFHKMXF2gxzbg36pu2XYPyiGShLqXk6L3KRs zwqYSdd8HMRhd0tQNrPWVGM2YREFdgnWmVSkEeXGlKgnDqq8Lv8gazFN1Skw JyDpwurSEUOcldizPRyOkPFZCJQ4CA0dhiPKZaOv+xIH/UbZDWzP9IzupRBD SVHsH3akBGUza7GpZe2wCasQSWVps+swDeTnJA6zL7nKriOjRrGskDGkeBTp Qi4q1Ne93Kxy7o0ohZPyEGicoy5R2ixxEDgqEGZ78I5JTrN/Zml2K/uiaDBx cOFIk2hmrbkhOEQU1OgCJXbJAjPsOrRFXXbALqEtGZfEQZeq+DQzY0RJtuw6 5KVWD9uzj4uwPRwxut6HuF/HmwqIAwvlFgoCS16ucGZoebE9OA81jhVLHkhQ RtZGy8sSd8hQuw6UHYFqpnlih6HZyWpjRyxBlzjLPGZ7lhTCDkNVUw4O2zMd m+1hK6pUmDg4eaTNEgfRo3e7lyJJ2f6ki/5KFHYr96acLcTBcSPxWLNWmXpk jl52yVQVkoJUUuyUuJ5MRM8kXYihaX0QZ/cecoiWysSuQxGXd4LtwahVkjjb w+0Uur+NONNmR5Ss0GFIunZTF+naAc72rIKCOKh0MjyaWWuhPnYJB4g0Q5LC qtNhs1L+DC5HiZCu5epRHjzbsxJgDtEu0mHXoWXqHCddS/YYUa4o+ba6lbtA PiviYLRInBFXVHllqUfEQe2U1CXO6hzYP/g4Ku3VsdbymUkKOSnaXCQFLVMh GrIW57DyT8haS1Fie5akxvYsMyuc43KKE2efxCauuPUqYw9xEGBylBNXpDJl rhfiIJlVRE+cGUvEWekQWQt6EpTNrLWsXJJClEdhlm7lRlKSxEolK+RuI2st 2Y7tWcYV27OKa+IgSPRhEuLsq5eki12mXcVxmDa7cv4MvxXBYv9wwClLgXTh maic915GQJyCAmV7day1CzfYBGwtnYkkBQtcLOMQkdShfFu2Z3ercYgI78ix zSHajVSkC1mOJItI16rJSdfUNdI1bZZ0IYRkR4VdpWgf6VrySNgt2rmrNrLR PALsJjaqbA3irB6Q3cRqkqVPHFQmfTOOOGgMynLicBDHldkaNGGZwRwOzFH1 hf2zMhzSNc115XJ1S5gQB/UNRkYWnWN7psKRrqWCNLPRStjZLDa5jqJwjClZ jN3E6ahLbImDz1RGAXHAaAWQjXb76oiSFSqkCLJRZutKJU4ZK6RrWirpwjyV 4UEc7D7tArZn/l/SNS26mY12+fZK5SqT74g4WO9ycbOb5jclzpLKiMOQlW5J NuJUletvRDktUrSIs/wYtmf5jqSLKjsdWSuXO00fkCIOJrVKkzheCBxdrUA2 WplmZGPGRTaGM03RxzAn8ocSZ/f9sHM4J3GIZKoOcfZRE+IsZEGmpIQ5VB/p FF3LtS8VLtCFSyWzqkjXylFWLte5Dk3i4L9VVRxxUNdkoZGulY80r0FL5Q1n nRx5LJmD50Xhk5H5M8vEJg7KubRPtodyWkndoJEqVBLUJjnZiEM6hnyapIv1 p+RO0jWpS7qYTiVyEgeDA3pFZkCMLMcrW5586czijiwUuVojG8EKXa5B8lh1 yrEneSx/Sc6RpdSQpMvLdtwpRpxdWR4EPGhn6hR3ASQf3s1sqxDK0dk5shQN sl3YPyuY4jjMnUIcWK1nbA+3akitS+kZzWy0JBqSh90iAU/ySG5RgITkwRpp fGSjlQERZ1817FJuOKU8kC60VCX+hINKJiXbM3cK6ZoGSVzPeEcnx4EDU2Mj zu5f6VKyW99JqGFj55KUBDebhTNDBwZxYI2CIWE4iiVyOOYcIQ5CQ5l7YXOB bZk1TRwCDVgZ2UbicGA+yiMwslyhilGRrmmLn8mf2b0DxIGFmvqwunXwpdS1 ZjbaBxZGlkPEgs7y54mDb1EBiJHlipIIZzdhXilOShxcw1opbA/HAXwz2QYO dBFwyDYsh2NfjyMOR6A2Jumatkg2IkFeaXLE4T35rEaWfVbwp0s5VfI7kY15 BX2n9sWzv/PC3Zc+ffg6A69YafvT/+OGhz+sXnvQ6mqEls9bT5BZJ/nAusMB plyTjn8WOvlMKRZ0qHDwls3EwRfXJGcTQZrQSyU12R6UJxxcmd1NnFk/bA9/ h26aqeOka+mXHIOlh3b8gzjBoZX14zPlPpUNTx5BCdP1kpyYnAwupJ19+RM7 rDZ6/XFvZWWHHU/jvQat7j5o9bxmVuzS3by33aCWKh7X8a/4Fvwq+SuwrLUL yBi764qMwStSbjr+YQHg+MsY/9lygmSIsC0oZ0rsZJfgq1JosuNfcVdE3qVu sLYlUzr+QfdCLugqpUhUALRLOSUKwRj7u1TZn7MKNdWtLjCY2fJigxrWQ92Q isBxYpfLO9LxD8oqXJH5QLpDLmN/FMducd1cttbIC1O4Ov7h+MWMjsqpwvDH LGRbjBy3CB57AnkohajjX6H6koqdCh3/oCvBZFillDpSYTheUIS1lEmo1ZeS 460uJPig5UUFNRy3IDXHCZbCEVP4i4rATL5Au0JVLZJGZFnmf+sWD3+wD85c JGLlinw3WGgyTjv+FWn25CBkiS7C7PhXKMX8M3xTciJ1/ANh7Ir8z8sVmVRQ fbGxEAaD8xerQf6dpWVyq+sEZra8ZqC81iHnLYRx5ncm2+BZLASvvGz8G1ZS EYKSTsW/xS/aQHfLhGa+O7tib0DuQCArv51TgdBDoScWWc7YJGQm5HPh0sHr RdoatxW83mgNti1UIpiSCFhlodc1q0ycdeVj265aMDFfbb1aVfa3ep5apGvk w4TFUQgr7TKyB2KvyB6GOyoLfvJvhZ8HHkUs6EzvyldAF1g3WG+F85WvQGcq UvJhB2ZSnX+DcQBe4BcMnsynQamNcDj0lEy6fz7nDyuAOnSKzSfedOZFh+26 /cb73TUNRfFkUaua+pbPW7MIUrGI9Un0s98QQsUe12cI8791Kf6GVQX/S95c F9hM2LtYJdBg2RIkAxz5WAXY+llSIf8G5hYmH06SzGL+cc6JlJLKjmO3InQC A1YfymOTSEQqRKaMdXYcghh/Qy8xh/mazNIXcDpkr+XyvUshJ4qL2/gYGwK2 JmQEWFLT17zxLhBfmOTCUmQ3IWyg78LfhD2TmQpcF6iGRkpw5rr7XP4MhjaI Z9Uj7AxOPUiFQp3iY6xMSMCiwqK5j8ijgC8S+hYWI9lVGORFX3n0IU7dvtHX zs8UibXyZ1CZdY7xdRwTSGkolAy+Du8ESqKgjbf9pKlvWeoNDJLMNZ0fCVk+ B4yKzNhgs9AXZc6uVT7T4cteQcQCVxzEfB2pz5jzjL/NvYLyhV5kphJ7hRnG usjuE2az0KYVZ/23/Jm8V/iXi84uUPpgYhWnF1+HNxF3lWayvKZXXCQYk9Lv 2QTsS+g6WfLbYfkz7B04ejJ/MXHw+WFWs/0SSWWUSAAeTjmU+SKkPF7MYp4k AENDLhDioMvoKojmseDklfONTUDlwVgy1ZekMGHKUx5VrhkVTTaTQnNKLWIT OLcwodmvw0ry1dHDlpVzq4bUFvkzaJTyzLJZHAuyhtgsJJwylpqbxUmiPDE2 iwNM/qLDWzeR/fcH/xp/b1vmfwFhN8Os\ \>"]] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ListPlot", "[", "%", "]"}]], "Input", CellChangeTimes->{{3.511724863088964*^9, 3.511724869495454*^9}}], Cell[BoxData[ GraphicsBox[ {Hue[0.67, 0.6, 0.6], PointBox[CompressedData[" 1:eJxTTMoPSmViYGAwAmIQDQEf7FFpBgc/ifp9q/zeQvkcDiLBcx2Opr6C8gUc ItTv22p9egnlizhkaZVaKce9gfIlHDrZiiw35L+H8mUclK8kf/bSgpmv4FAR FXY50eodlK/kEHT/dcwz/tdQvoqDioDFLa99MPPVHNxOm0V+qYbJazgI3lpY p7oepl/L4VXdi7ooY5j5Og7HG05NbuSD2a/n0LsqXIPhMUy/gUPZ30MvPc7B zDd02BaROFdpF8x/Rg67zph4bDsL87+xQ8XBxjulUjDzTRwsfk4zVNkFM9/U 4Z398k2O52D+N3PYO0HU+54CzDxzhzxlpR45DRjfwqHjpOqV8Lsw9ZYOyk8d EwwPw8yzcvilEHfsnDjMPmsHk30rit/ug7nHxmGa9o+HZeth5tk6RGw3udR5 EeYfO4f3hkLv9T/C/Gvv0F/8W2+qFMx8B4cnPUVMJ42g5jc4OHCXxZVdXwgL T0cH1p0q+7flQfU3ODqUrtkh4LUXZr6Tw/Jt+6v4ZGDyTg55DI1TWj1g+p0d Kp5WLb+lDTPf2YHpAnfAmxSY/S4OcnPsjib6Q/3f4OLAoslt/ucNzHxXhzvP d515mwf1X4OrQ8nLFT77YmH+d3OYyjqBdyIzzHw3h6zPjvsW/oCZ7+6wzUAp fJI9VH2Du8ON4/fmXo6BhZeHgwZ3TMmM71D7Gjwc/NxylT7nwOLD08HrlfAE 56r39gCxcOg7 "]]}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesOrigin->{0, Automatic}, PlotRange->Automatic, PlotRangeClipping->True]], "Output", CellChangeTimes->{3.511724869953926*^9}] }, Open ]], Cell[BoxData[ RowBox[{"Manipulate", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{"ListLinePlot", "[", RowBox[{"Table", "[", RowBox[{ RowBox[{"y", "[", RowBox[{"a", ",", "t"}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", RowBox[{"10", "\[Pi]"}]}], "}"}]}], "]"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"a", ",", "0.1", ",", "2", ",", "0.01"}], "}"}]}], "\[IndentingNewLine]", "]"}]], "Input", CellChangeTimes->{{3.5084702833515253`*^9, 3.508470371260749*^9}, { 3.5084704049457808`*^9, 3.508470585633668*^9}, {3.508470630221963*^9, 3.508470685659617*^9}, {3.508471934868383*^9, 3.5084719502552233`*^9}, { 3.508471987707246*^9, 3.508471993591114*^9}, {3.508472071185714*^9, 3.50847208369586*^9}, {3.508472114869843*^9, 3.50847222598955*^9}, { 3.508472274639151*^9, 3.5084722924350443`*^9}, {3.508472818627358*^9, 3.5084728208405523`*^9}, {3.50847307799758*^9, 3.5084730788812113`*^9}, { 3.508473124957897*^9, 3.508473310091757*^9}, {3.508473348585227*^9, 3.508473379803138*^9}, {3.508473410322609*^9, 3.5084734234021263`*^9}, { 3.511724621571731*^9, 3.511724728636238*^9}}], Cell[BoxData[""], "Input", CellChangeTimes->{3.508473051146953*^9}], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{"Manipulate", "[", RowBox[{ RowBox[{"ListLinePlot", "[", RowBox[{"{", "\[IndentingNewLine]", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{"t", ",", RowBox[{"Cos", "[", "t", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", RowBox[{"10", "\[Pi]"}], ",", "0.1"}], "}"}]}], "]"}], ",", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{"t", ",", RowBox[{"y", "[", RowBox[{"a", ",", "t"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", RowBox[{"10", "\[Pi]"}]}], "}"}]}], "]"}]}], "\[IndentingNewLine]", "}"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"a", ",", "0.001", ",", "2", ",", "0.01"}], "}"}]}], "\[IndentingNewLine]", "]"}]}]], "Input", CellChangeTimes->{{3.508470711826589*^9, 3.5084707915768433`*^9}, { 3.508470826600807*^9, 3.508470929022626*^9}, {3.508470986360839*^9, 3.508470987944158*^9}, {3.508471284848448*^9, 3.508471295849*^9}, { 3.508471965053668*^9, 3.508471967654603*^9}, {3.508472095609482*^9, 3.508472108575961*^9}, {3.508472254938147*^9, 3.508472285074898*^9}, { 3.508472320290619*^9, 3.508472320500173*^9}, {3.508472881119774*^9, 3.508472898927884*^9}, {3.5084729624861727`*^9, 3.508472975288372*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ListLinePlot", "[", "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{"t", ",", RowBox[{"y", "[", RowBox[{"0.1", ",", "t"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", RowBox[{"10", "\[Pi]"}]}], "}"}]}], "]"}], "\[IndentingNewLine]", "]"}]], "Input", CellChangeTimes->{{3.508472924788168*^9, 3.508472952651958*^9}, { 3.5084729948444633`*^9, 3.508473028354723*^9}}], Cell[BoxData[ GraphicsBox[{}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesOrigin->{0, 0}, PlotRange->{{0., 0.}, {0., 0.}}, PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Scaled[0.02]}]], "Output", CellChangeTimes->{{3.508472926166854*^9, 3.5084730287928*^9}}] }, Open ]], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"RUNGE", "-", RowBox[{"KUTTA", ":"}]}]}]], "Input", CellChangeTimes->{{3.508471214733235*^9, 3.5084712294976597`*^9}, 3.5084712749310827`*^9}], Cell[BoxData[{ RowBox[{ RowBox[{"k1", "[", RowBox[{"a_", ",", "t_"}], "]"}], " ", ":=", " ", RowBox[{ RowBox[{"-", " ", "a"}], " ", RowBox[{"Sin", "[", "t", "]"}]}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"k2", "[", RowBox[{"a_", ",", "t_"}], "]"}], " ", ":=", " ", RowBox[{ RowBox[{"-", " ", "a"}], " ", RowBox[{"Sin", "[", RowBox[{"t", "-", FractionBox["a", "2"]}], "]"}]}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"k3", "[", RowBox[{"a_", ",", "t_"}], "]"}], " ", ":=", " ", RowBox[{ RowBox[{"-", " ", "a"}], " ", RowBox[{"Sin", "[", RowBox[{"t", "-", FractionBox["a", "2"]}], "]"}]}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"k4", "[", RowBox[{"a_", ",", "t_"}], "]"}], " ", ":=", " ", RowBox[{ RowBox[{"-", " ", "a"}], " ", RowBox[{"Sin", "[", RowBox[{"t", "-", "a"}], "]"}]}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "[", RowBox[{"a_", ",", "t_"}], "]"}], " ", ":=", " ", RowBox[{"y", "[", RowBox[{"a", ",", RowBox[{"t", "-", "1"}]}], "]"}]}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.5084712609342213`*^9, 3.508471267860257*^9}, 3.5084713119703217`*^9, {3.5084725616463614`*^9, 3.508472781036109*^9}}] }, WindowSize->{1280, 747}, WindowMargins->{{1, Automatic}, {-24, Automatic}}, FrontEndVersion->"8.0 for Linux x86 (64-bit) (November 14, 2010)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[557, 20, 262, 7, 30, "Input"], Cell[822, 29, 820, 24, 88, "Input"], Cell[1645, 55, 96, 1, 30, "Input"], Cell[1744, 58, 1291, 30, 69, "Input"], Cell[CellGroupData[{ Cell[3060, 92, 147, 3, 30, "Input"], Cell[3210, 97, 88, 1, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[3335, 103, 264, 7, 30, "Input"], Cell[3602, 112, 1601, 28, 88, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[5240, 145, 196, 4, 30, "Input"], Cell[5439, 151, 24863, 414, 239, 1001, 21, "CachedBoxData", "BoxData", \ "Output"] }, Open ]], Cell[CellGroupData[{ Cell[30339, 570, 127, 2, 30, "Input"], Cell[30469, 574, 1044, 22, 239, "Output"] }, Open ]], Cell[31528, 599, 1192, 24, 88, "Input"], Cell[32723, 625, 68, 1, 30, "Input"], Cell[32794, 628, 1478, 34, 126, "Input"], Cell[CellGroupData[{ Cell[34297, 666, 488, 13, 69, "Input"], Cell[34788, 681, 325, 10, 237, "Output"] }, Open ]], Cell[35128, 694, 231, 5, 69, "Input"], Cell[35362, 701, 1259, 37, 173, "Input"] } ] *) (* End of internal cache information *)