(* 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[ 127088, 2824] NotebookOptionsPosition[ 123297, 2752] NotebookOutlinePosition[ 123645, 2767] CellTagsIndexPosition[ 123602, 2764] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell["\:041c\:0435\:0442\:043e\:0434\:0438", "Title", CellChangeTimes->{{3.747918045893719*^9, 3.7479180792016163`*^9}, { 3.748019567260964*^9, 3.7480195984310884`*^9}, {3.7491799306613116`*^9, 3.749179945330095*^9}, {3.7491820390416107`*^9, 3.749182140816432*^9}, { 3.7491822953940344`*^9, 3.7491823155841894`*^9}, {3.749182378833807*^9, 3.749182424987447*^9}, {3.749182543546228*^9, 3.74918254777547*^9}, { 3.7491855618392*^9, 3.749185791014152*^9}, 3.7494816355928173`*^9, { 3.7494826582226067`*^9, 3.749482678435937*^9}, {3.749483068682683*^9, 3.7494830721446886`*^9}, {3.7494884082644987`*^9, 3.7494884424507885`*^9}, {3.749488494590744*^9, 3.7494885492422333`*^9}, { 3.749488662765763*^9, 3.7494886848150673`*^9}, {3.7494887168485804`*^9, 3.749488740273077*^9}, 3.749490389264352*^9, {3.749789440131727*^9, 3.749789499828309*^9}, {3.74979010532883*^9, 3.7497901250714474`*^9}, 3.7499847604043837`*^9, {3.7499852544136286`*^9, 3.749985272769679*^9}, { 3.7515612145185156`*^9, 3.751561319842022*^9}, {3.751561379568722*^9, 3.751561400725582*^9}, {3.7521741297192936`*^9, 3.752174150598387*^9}}, FontFamily-> "Times New Roman",ExpressionUUID->"467d1ec9-6d84-4b2a-90de-bdb8eba1d886"], Cell[BoxData[ RowBox[{ RowBox[{"explicitEuler", "[", RowBox[{"f_", ",", "h_", ",", "t0_", ",", "T_", ",", "u0_"}], "]"}], ":=", " ", RowBox[{"(", "\[IndentingNewLine]", RowBox[{ RowBox[{"n", "=", RowBox[{"Ceiling", "[", RowBox[{ RowBox[{"(", RowBox[{"T", "-", "t0"}], ")"}], "/", "h"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"t", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"t0", "+", RowBox[{"i", "*", "h"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "0", ",", "n"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"y", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"n", "+", "1"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "1", "]"}], "]"}], " ", "=", "u0"}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", "<", RowBox[{"n", "+", "1"}]}], ",", " ", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "[", RowBox[{"[", RowBox[{"i", "+", "1"}], "]"}], "]"}], " ", "=", " ", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}], "+", RowBox[{"h", "*", RowBox[{"f", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}]}]}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"Transpose", "[", RowBox[{"{", RowBox[{"t", ",", "y"}], "}"}], "]"}]}], "\[IndentingNewLine]", ")"}]}]], "Input", CellChangeTimes->{ 3.7515563087600636`*^9, {3.751556360006592*^9, 3.7515564799865913`*^9}, { 3.751556639141062*^9, 3.751556639357938*^9}, {3.75155671085882*^9, 3.751556714680622*^9}, {3.751556782540594*^9, 3.7515567830503016`*^9}, { 3.751557223977725*^9, 3.751557232127038*^9}, {3.751557273543219*^9, 3.75155727616871*^9}, {3.751557439493783*^9, 3.751557439852577*^9}, { 3.751558318906334*^9, 3.751558319058245*^9}, {3.7515586051607094`*^9, 3.75155861412757*^9}, {3.7521570890106964`*^9, 3.75215709038891*^9}, { 3.752176074470169*^9, 3.7521760830142946`*^9}, {3.784026399199913*^9, 3.784026404691227*^9}},ExpressionUUID->"d0a40c0c-c125-4a61-9539-\ c7544c932fc8"], Cell[BoxData[ RowBox[{ RowBox[{"implicitEuler", "[", RowBox[{"f_", ",", "h_", ",", "t0_", ",", "T_", ",", "u0_"}], "]"}], ":=", " ", RowBox[{"(", "\[IndentingNewLine]", RowBox[{ RowBox[{"n", "=", RowBox[{"Ceiling", "[", RowBox[{ RowBox[{"(", RowBox[{"T", "-", "t0"}], ")"}], "/", "h"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"t", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"t0", "+", RowBox[{"i", "*", "h"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "0", ",", "n"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"y", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"n", "+", "1"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "1", "]"}], "]"}], " ", "=", "u0"}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", "<", RowBox[{"n", "+", "1"}]}], ",", " ", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"initial", " ", "=", " ", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}], "+", RowBox[{"h", "*", RowBox[{"f", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}]}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "[", RowBox[{"[", RowBox[{"i", "+", "1"}], "]"}], "]"}], "=", RowBox[{"yNew", "/.", RowBox[{"FindRoot", "[", RowBox[{ RowBox[{"yNew", "==", " ", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}], "+", RowBox[{"h", "*", RowBox[{"f", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", "+", "1"}], "]"}], "]"}], ",", "yNew"}], "]"}]}]}]}], ",", RowBox[{"{", RowBox[{"yNew", ",", "initial"}], "}"}]}], "]"}]}]}], ";"}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", RowBox[{"Transpose", "[", RowBox[{"{", RowBox[{"t", ",", "y"}], "}"}], "]"}]}], "\[IndentingNewLine]", ")"}]}]], "Input", CellChangeTimes->{{3.7515574834674945`*^9, 3.7515575998356133`*^9}, { 3.751557648158707*^9, 3.7515576636607933`*^9}, {3.751557738529467*^9, 3.751557763031377*^9}, 3.7515580336403885`*^9, {3.7515580964512835`*^9, 3.751558123797539*^9}, {3.7515581689555693`*^9, 3.751558169385322*^9}, { 3.751558322920025*^9, 3.7515583230319777`*^9}, {3.751558458330151*^9, 3.7515585241552973`*^9}, {3.751570772243475*^9, 3.751570817505479*^9}, { 3.7521766927406673`*^9, 3.752176696682413*^9}, {3.784026407497388*^9, 3.784026425861438*^9}},ExpressionUUID->"4a6b5a71-ee89-4044-925a-\ 9e0eca7f52ca"], Cell[BoxData[ RowBox[{ RowBox[{"RK4", "[", RowBox[{"f_", ",", "h_", ",", "t0_", ",", "T_", ",", "u0_"}], "]"}], ":=", " ", RowBox[{"(", "\[IndentingNewLine]", RowBox[{ RowBox[{"n", "=", RowBox[{"Ceiling", "[", RowBox[{ RowBox[{"(", RowBox[{"T", "-", "t0"}], ")"}], "/", "h"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"t", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"t0", "+", RowBox[{"i", "*", "h"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "0", ",", "n"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"y", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"i", ",", " ", "0", ",", " ", "n"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "1", "]"}], "]"}], " ", "=", " ", "u0"}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", "<=", "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[{ RowBox[{"1", "/", "6"}], "k1"}], "+", RowBox[{ RowBox[{"1", "/", "3"}], "k2"}], "+", RowBox[{ RowBox[{"1", "/", "3"}], "k3"}], "+", RowBox[{ RowBox[{"1", "/", "6"}], "k4"}]}]}]}]}], "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"Transpose", "[", RowBox[{"{", RowBox[{"t", ",", "y"}], "}"}], "]"}]}], "\[IndentingNewLine]", ")"}]}]], "Input", CellChangeTimes->{{3.749790138229862*^9, 3.7497901566012683`*^9}, 3.7499849756016817`*^9, {3.750398226091392*^9, 3.750398389332283*^9}, { 3.7503984331110306`*^9, 3.750398439266488*^9}, {3.750398574497512*^9, 3.7503985813046026`*^9}, {3.7503986122967186`*^9, 3.75039861385282*^9}, { 3.7503987198976784`*^9, 3.7503987590611153`*^9}, {3.750861148761247*^9, 3.7508611531587133`*^9}, {3.750861592248006*^9, 3.7508615945696664`*^9}, 3.7508617880070868`*^9, 3.750861941006491*^9, {3.7508691323068924`*^9, 3.7508691387071886`*^9}, {3.750869283890499*^9, 3.750869292626458*^9}, 3.751556304862305*^9, {3.751558296542212*^9, 3.7515583069642015`*^9}, { 3.751560941093563*^9, 3.7515609415253134`*^9}, {3.7515692698583946`*^9, 3.7515692734823074`*^9}, {3.7521767011338625`*^9, 3.752176704922719*^9}, { 3.783328759467596*^9, 3.783328760144196*^9}, {3.7840264309797306`*^9, 3.7840264421173677`*^9}},ExpressionUUID->"218c4638-48db-45f8-9652-\ b892485b01bc"], Cell[BoxData[ RowBox[{ RowBox[{"AB4", "[", RowBox[{"f_", ",", "h_", ",", "t0_", ",", "T_", ",", "u0_"}], "]"}], ":=", RowBox[{"(", "\[IndentingNewLine]", RowBox[{ RowBox[{"n", "=", RowBox[{"Ceiling", "[", RowBox[{ RowBox[{"(", RowBox[{"T", "-", "t0"}], ")"}], "/", "h"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"t", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"t0", "+", RowBox[{"i", " ", "h"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "0", ",", "n"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"y", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"n", "+", "1"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"fi", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"n", "+", "1"}], "}"}]}], "]"}]}], ";", " ", RowBox[{"(*", RowBox[{ RowBox[{ "list", " ", "of", " ", "values", " ", "of", " ", "the", " ", "right"}], "-", RowBox[{"hand", " ", "side", " ", "f"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "1", "]"}], "]"}], "=", "u0"}], ";", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{ "Compute", " ", "the", " ", "first", " ", "approximate", " ", "values"}], ",", " ", RowBox[{"using", " ", "the", " ", "RK4", " ", "method"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", "\[LessEqual]", "3"}], ",", " ", 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]", RowBox[{ RowBox[{"fi", "[", RowBox[{"[", RowBox[{"i", "+", "1"}], "]"}], "]"}], "=", RowBox[{"f", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", "+", "1"}], "]"}], "]"}], ",", RowBox[{"y", "[", RowBox[{"[", RowBox[{"i", "+", "1"}], "]"}], "]"}]}], "]"}]}], ";"}]}], "\[IndentingNewLine]", "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"Compute", " ", "the", " ", "remaining", " ", "values"}], ",", " ", RowBox[{"using", " ", "the", " ", "AB4", " ", "formula"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "4"}], ",", RowBox[{"i", "<", RowBox[{"n", "+", "1"}]}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"y", "[", RowBox[{"[", RowBox[{"i", "+", "1"}], "]"}], "]"}], "=", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}], "+", RowBox[{ FractionBox[ RowBox[{"h", " "}], "24"], RowBox[{"(", RowBox[{ RowBox[{"55", " ", RowBox[{"fi", "[", RowBox[{"[", "i", "]"}], "]"}]}], "-", " ", RowBox[{"59", " ", RowBox[{"fi", "[", RowBox[{"[", RowBox[{"i", "-", "1"}], "]"}], "]"}]}], "+", " ", RowBox[{"37", " ", RowBox[{"fi", "[", RowBox[{"[", RowBox[{"i", "-", "2"}], "]"}], "]"}]}], "-", " ", RowBox[{"9", RowBox[{"fi", "[", RowBox[{"[", RowBox[{"i", "-", "3"}], "]"}], "]"}]}]}], ")"}]}]}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"fi", "[", RowBox[{"[", RowBox[{"i", "+", "1"}], "]"}], "]"}], "=", RowBox[{"f", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", "+", "1"}], "]"}], "]"}], ",", RowBox[{"y", "[", RowBox[{"[", RowBox[{"i", "+", "1"}], "]"}], "]"}]}], "]"}]}], ";"}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", RowBox[{"Transpose", "[", RowBox[{"{", RowBox[{"t", ",", "y"}], "}"}], "]"}]}], "\[IndentingNewLine]", ")"}]}]], "Input", CellChangeTimes->{{3.7515586981543894`*^9, 3.7515587517855453`*^9}, { 3.75156929377765*^9, 3.751569300563752*^9}, {3.7521767135187993`*^9, 3.752176717614464*^9}, {3.7521771841507807`*^9, 3.752177189484725*^9}, { 3.7521773558228617`*^9, 3.752177356058737*^9}, {3.7833287667353816`*^9, 3.783328767463978*^9}, {3.784026447369668*^9, 3.784026466158743*^9}},ExpressionUUID->"36837f17-bac1-40b1-b9f2-\ 4ea41252841b"], Cell[BoxData[""], "Input", CellChangeTimes->{{3.7840231080506697`*^9, 3.784023109024726*^9}},ExpressionUUID->"5625b212-3c6c-4992-af0f-\ 9cca96159dfb"] }, Open ]], Cell[CellGroupData[{ Cell[TextData[{ "\:041d\:0415 \:0415 \:0420\:0415\:0410\:041b\:0418\:0417\:0418\:0420\:0410\ \:041d\:041e \:041d\:0410 \:0423\:041f\:0420\:0410\:0416\:041d\:0415\:041d\ \:0418\:0415!\n\:0422\:044f\:043b\:043e, \:043f\:0430\:0434\:0430\:0449\:043e \ \:0441\:0432\:043e\:0431\:043e\:0434\:043d\:043e \:043a\:044a\:043c \:0437\ \:0435\:043c\:043d\:0430\:0442\:0430 \:043f\:043e\:0432\:044a\:0440\:0445\ \:043d\:043e\:0441\:0442, \:043c\:043e\:0436\:0435 \:0434\:0430 \:0431\:044a\ \:0434\:0435 \:043c\:043e\:0434\:0435\:043b\:0438\:0440\:0430\:043d\:043e \ \:0447\:0440\:0435\:0437 \:0441\:0438\:0441\:0442\:0435\:043c\:0430\:0442\ \:0430 \:041e\:0414\:0423 \:043e\:0442 \:0432\:0442\:043e\:0440\:0438 \:0440\ \:0435\:0434:\n", Cell[BoxData[ FormBox[ RowBox[{ RowBox[{ RowBox[{"x", "''"}], RowBox[{"(", "t", ")"}]}], "=", "9.8"}], TraditionalForm]], ExpressionUUID->"9d74bb27-a696-48f0-9ff4-83ceefcc3801"], ", 0 "e9d01ce1-39a4-41c7-8d7b-51a4b32a04bc"], "(t)=v. \:0422\:043e\:0433\:0430\:0432\:0430:\n", Cell[BoxData[{ FormBox[ RowBox[{ RowBox[{"x", "'"}], "=", "v"}], TraditionalForm], "\[IndentingNewLine]", FormBox[ RowBox[{ RowBox[{"v", "'"}], "=", "9.8"}], TraditionalForm], "\[IndentingNewLine]", FormBox[ RowBox[{ RowBox[{"x", "(", "0", ")"}], "=", "0"}], TraditionalForm], "\[IndentingNewLine]", FormBox[ RowBox[{ RowBox[{"v", "(", "0", ")"}], "=", "0"}], TraditionalForm]}], ExpressionUUID->"1fc68b38-1eb6-4c00-8bbf-b7cb9d7b846a"], "\n\:0438\:043b\:0438 \n\:0421\:0432\:0435\:0436\:0434\:0430\:043c\:0435 \ \:0437\:0430\:0434\:0430\:0447\:0430\:0442\:0430 \:043a\:044a\:043c \:0437\ \:0430\:0434\:0430\:0447\:0430\:0442\:0430 \:043d\:0430 \:041a\:043e\:0448\ \:0438 \:0441 \:043f\:043e\:043b\:0430\:0433\:0430\:043d\:0435\:0442\:043e ", Cell[BoxData[ FormBox[ SuperscriptBox["x", "'"], TraditionalForm]],ExpressionUUID-> "4abbdd28-3c06-47e4-aaf8-15a977473994"], "(t)=v. \:0422\:043e\:0433\:0430\:0432\:0430:\n", Cell[BoxData[{ FormBox[ RowBox[{ RowBox[{"u", "'"}], "=", RowBox[{"f", "(", RowBox[{"t", ",", "u"}], ")"}]}], TraditionalForm], "\[IndentingNewLine]", FormBox[ RowBox[{ RowBox[{"u", "(", "0", ")"}], "=", SubscriptBox["u", "0"]}], TraditionalForm]}],ExpressionUUID-> "04b71b96-ad85-4927-9544-e7d95b94b8a3"], "\n\:043a\:044a\:0434\:0435\:0442\:043e u=(x,v)" }], "Title", CellChangeTimes->{{3.747918045893719*^9, 3.7479180792016163`*^9}, { 3.748019567260964*^9, 3.7480195984310884`*^9}, {3.7491799306613116`*^9, 3.749179945330095*^9}, {3.7491820390416107`*^9, 3.749182140816432*^9}, { 3.7491822953940344`*^9, 3.7491823155841894`*^9}, {3.749182378833807*^9, 3.749182424987447*^9}, {3.749182543546228*^9, 3.74918254777547*^9}, { 3.7491855618392*^9, 3.749185791014152*^9}, 3.7494816355928173`*^9, { 3.7494826582226067`*^9, 3.749482678435937*^9}, {3.749483068682683*^9, 3.7494830721446886`*^9}, {3.7494884082644987`*^9, 3.7494884424507885`*^9}, {3.749488494590744*^9, 3.7494885492422333`*^9}, { 3.749488662765763*^9, 3.7494886848150673`*^9}, {3.7494887168485804`*^9, 3.749488740273077*^9}, 3.749490389264352*^9, {3.749789440131727*^9, 3.749789499828309*^9}, {3.74979010532883*^9, 3.7497901250714474`*^9}, 3.7499847604043837`*^9, {3.7499852544136286`*^9, 3.749985272769679*^9}, { 3.7515612145185156`*^9, 3.751561319842022*^9}, {3.751561379568722*^9, 3.751561400725582*^9}, {3.7521741959964952`*^9, 3.7521742399644175`*^9}, { 3.7521774008681836`*^9, 3.7521775625919323`*^9}, {3.7840265607431526`*^9, 3.784026599398364*^9}, {3.7840270730204535`*^9, 3.7840270754825945`*^9}, { 3.784040591724004*^9, 3.7840406078804197`*^9}, {3.7840410291265182`*^9, 3.7840410350541124`*^9}}, FontFamily-> "Times New Roman",ExpressionUUID->"33bc2259-f836-4b21-81b1-d28abbbb91b9"], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"Clear", "[", RowBox[{"x", ",", "t"}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"exact1", "[", "t_", "]"}], "=", RowBox[{ RowBox[{"x", "[", "t", "]"}], "/.", RowBox[{ RowBox[{"DSolve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{"x", "''"}], "[", "t", "]"}], "\[Equal]", "9.8"}], ",", RowBox[{ RowBox[{"x", "[", "0", "]"}], "\[Equal]", "0"}], ",", " ", RowBox[{ RowBox[{ RowBox[{"x", "'"}], "[", "0", "]"}], "\[Equal]", "0"}]}], "}"}], ",", RowBox[{"x", "[", "t", "]"}], ",", "t"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"h", "=", "0.01"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"t", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"0", "+", RowBox[{"i", "*", "h"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "0", ",", RowBox[{"10", "/", "h"}]}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ "\:0429\:0435", " ", "\:043c\:043e\:0434\:0438\:0444\:0438\:0446\:0438\:0440\:0430\:043c\:0435\ ", " ", "\:043a\:043e\:0434\:0430", " ", "\:0437\:0430", " ", "\:043c\:0435\:0442\:043e\:0434\:0430", " ", "\:043d\:0430", " ", "\:041e\:0439\:043b\:0435\:0440", " ", "\:0434\:0430", " ", "\:0440\:0430\:0431\:043e\:0442\:0438", " ", "\:0437\:0430", " ", "2", " ", "\:0443\:0440\:0430\:0432\:043d\:0435\:043d\:0438\:044f"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"explicitEuler", "[", RowBox[{ StyleBox["f1_", Background->RGBColor[1, 1, 0]], StyleBox[",", Background->RGBColor[1, 1, 0]], StyleBox["f2_", Background->RGBColor[1, 1, 0]], ",", "h_", ",", "t0_", ",", "T_", ",", StyleBox["x0_", Background->RGBColor[1, 1, 0]], StyleBox[",", Background->RGBColor[1, 1, 0]], StyleBox["v0_", Background->RGBColor[1, 1, 0]]}], "]"}], ":=", " ", RowBox[{"(", "\[IndentingNewLine]", RowBox[{ RowBox[{"n", "=", RowBox[{"Ceiling", "[", RowBox[{ RowBox[{"(", RowBox[{"T", "-", "t0"}], ")"}], "/", "h"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"t", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"t0", "+", RowBox[{"i", "*", "h"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "0", ",", "n"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", StyleBox[ RowBox[{"x", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"n", "+", "1"}], "}"}]}], "]"}]}], Background->RGBColor[1, 1, 0]], StyleBox[";", Background->RGBColor[1, 1, 0]], StyleBox["\[IndentingNewLine]", Background->RGBColor[1, 1, 0]], StyleBox[ RowBox[{"v", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"n", "+", "1"}], "}"}]}], "]"}]}], Background->RGBColor[1, 1, 0]], StyleBox[";", Background->RGBColor[1, 1, 0]], "\[IndentingNewLine]", "\[IndentingNewLine]", StyleBox[ RowBox[{ RowBox[{"x", "[", RowBox[{"[", "1", "]"}], "]"}], " ", "=", "x0"}], Background->RGBColor[1, 1, 0]], StyleBox[";", Background->RGBColor[1, 1, 0]], StyleBox["\[IndentingNewLine]", Background->RGBColor[1, 1, 0]], StyleBox[ RowBox[{ RowBox[{"v", "[", RowBox[{"[", "1", "]"}], "]"}], "=", "v0"}], Background->RGBColor[1, 1, 0]], StyleBox[";", Background->RGBColor[1, 1, 0]], "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", "<", RowBox[{"n", "+", "1"}]}], ",", " ", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", StyleBox[ RowBox[{ RowBox[{ RowBox[{"x", "[", RowBox[{"[", RowBox[{"i", "+", "1"}], "]"}], "]"}], " ", "=", " ", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], "+", RowBox[{"h", "*", RowBox[{"f1", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"v", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}]}]}], ";", " ", "\[IndentingNewLine]", RowBox[{ RowBox[{"v", "[", RowBox[{"[", RowBox[{"i", "+", "1"}], "]"}], "]"}], " ", "=", " ", RowBox[{ RowBox[{"v", "[", RowBox[{"[", "i", "]"}], "]"}], "+", RowBox[{"h", "*", RowBox[{"f2", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"v", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}]}]}]}], Background->RGBColor[1, 1, 0]]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"Transpose", "[", RowBox[{"{", RowBox[{"t", ",", "x", ",", "v"}], "}"}], "]"}]}], "\[IndentingNewLine]", ")"}]}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"f1", "[", RowBox[{"t_", ",", "x_", ",", "v_"}], "]"}], ":=", "v"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"f2", "[", RowBox[{"t_", ",", "x_", ",", "v_"}], "]"}], ":=", "9.8"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"apprPos1", " ", "=", " ", RowBox[{ RowBox[{"explicitEuler", "[", RowBox[{ "f1", ",", "f2", ",", "h", ",", "0", ",", "10", ",", "0", ",", "0"}], "]"}], "[", RowBox[{"[", RowBox[{"All", ",", " ", "2"}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Animate", "[", RowBox[{ RowBox[{"Graphics", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", RowBox[{"-", RowBox[{"exact1", "[", RowBox[{"t", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}]}], "}"}], ",", "5"}], "]"}], ",", " ", RowBox[{"{", RowBox[{"Red", ",", "Dashed", ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", RowBox[{"-", RowBox[{"apprPos1", "[", RowBox[{"[", "i", "]"}], "]"}]}]}], "}"}], ",", "5"}], "]"}]}], "}"}]}], "}"}], ",", "\[IndentingNewLine]", "\t\t\t\t\t\t\t ", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "10"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "550"}], ",", "0"}], "}"}]}], "}"}]}], ",", " ", RowBox[{"Axes", "\[Rule]", RowBox[{"{", RowBox[{"False", ",", "True"}], "}"}]}]}], "]"}], ",", "\[IndentingNewLine]", "\t\t ", RowBox[{"{", RowBox[{"i", ",", "1", ",", RowBox[{"Length", "[", "t", "]"}], ",", "1"}], "}"}]}], "]"}]}], "Input", CellChangeTimes->{{3.78403846399135*^9, 3.7840388448801155`*^9}, { 3.7840389393036213`*^9, 3.7840389414753685`*^9}, {3.784038998827242*^9, 3.7840390022022915`*^9}, {3.7840397835273604`*^9, 3.78403982749765*^9}, { 3.784040432608474*^9, 3.7840404428458624`*^9}, {3.784040662930119*^9, 3.784040668872994*^9}},ExpressionUUID->"49b1b934-b068-4cfd-b851-\ 5abc5cdb112c"], Cell[BoxData[ TagBox[ StyleBox[ DynamicModuleBox[{$CellContext`i$$ = 327, Typeset`show$$ = True, Typeset`bookmarkList$$ = {}, Typeset`bookmarkMode$$ = "Menu", Typeset`animator$$, Typeset`animvar$$ = 1, Typeset`name$$ = "\"untitled\"", Typeset`specs$$ = {{ Hold[$CellContext`i$$], 1, 1001, 1}}, Typeset`size$$ = { 36., {213., 219.}}, Typeset`update$$ = 0, Typeset`initDone$$, Typeset`skipInitDone$$ = True, $CellContext`i$4048$$ = 0}, DynamicBox[Manipulate`ManipulateBoxes[ 1, StandardForm, "Variables" :> {$CellContext`i$$ = 1}, "ControllerVariables" :> { Hold[$CellContext`i$$, $CellContext`i$4048$$, 0]}, "OtherVariables" :> { Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$, Typeset`animator$$, Typeset`animvar$$, Typeset`name$$, Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$, Typeset`skipInitDone$$}, "Body" :> Graphics[{ Circle[{0, -$CellContext`exact1[ Part[$CellContext`t, $CellContext`i$$]]}, 5], {Red, Dashed, Circle[{0, -Part[$CellContext`apprPos1, $CellContext`i$$]}, 5]}}, PlotRange -> {{0, 10}, {-550, 0}}, Axes -> {False, True}], "Specifications" :> {{$CellContext`i$$, 1, 1001, 1, AppearanceElements -> { "ProgressSlider", "PlayPauseButton", "FasterSlowerButtons", "DirectionButton"}}}, "Options" :> { ControlType -> Animator, AppearanceElements -> None, DefaultBaseStyle -> "Animate", DefaultLabelStyle -> "AnimateLabel", SynchronousUpdating -> True, ShrinkingDelay -> 10.}, "DefaultOptions" :> {}], ImageSizeCache->{332., {252., 258.}}, SingleEvaluation->True], Deinitialization:>None, DynamicModuleValues:>{}, SynchronousInitialization->True, UndoTrackedVariables:>{Typeset`show$$, Typeset`bookmarkMode$$}, UnsavedVariables:>{Typeset`initDone$$}, UntrackedVariables:>{Typeset`size$$}], "Animate", Deployed->True, StripOnInput->False], Manipulate`InterpretManipulate[1]]], "Output", CellChangeTimes->{{3.7840411793008204`*^9, 3.7840411931528234`*^9}},ExpressionUUID->"da393554-3426-4f13-bc00-\ 8005af2c961f"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"apprPos1", "[", RowBox[{"[", "90", "]"}], "]"}]], "Input", CellChangeTimes->{{3.784041099337987*^9, 3.784041149686917*^9}},ExpressionUUID->"19b20860-11d6-4cda-ba18-\ 208c62396fcc"], Cell[BoxData["3.8376799999999975`"], "Output", CellChangeTimes->{{3.7840387911821175`*^9, 3.7840388198075905`*^9}, { 3.784040433499961*^9, 3.7840404431472464`*^9}, 3.7840405399933453`*^9, { 3.7840406461920223`*^9, 3.7840406706886873`*^9}, 3.784040960608075*^9, 3.784041003821128*^9, {3.7840410918453107`*^9, 3.7840411503305454`*^9}, 3.784041193189803*^9},ExpressionUUID->"a20753d2-d829-4324-94a5-\ 699220be10c2"] }, Open ]], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"explicitEuler", "[", RowBox[{ StyleBox["f1_", Background->RGBColor[1, 1, 0]], StyleBox[",", Background->RGBColor[1, 1, 0]], StyleBox["f2_", Background->RGBColor[1, 1, 0]], ",", "h_", ",", "t0_", ",", "T_", ",", StyleBox["x0_", Background->RGBColor[1, 1, 0]], StyleBox[",", Background->RGBColor[1, 1, 0]], StyleBox["v0_", Background->RGBColor[1, 1, 0]]}], "]"}], ":=", " ", RowBox[{"(", "\[IndentingNewLine]", RowBox[{ RowBox[{"n", "=", RowBox[{"Ceiling", "[", RowBox[{ RowBox[{"(", RowBox[{"T", "-", "t0"}], ")"}], "/", "h"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"t", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"t0", "+", RowBox[{"i", "*", "h"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "0", ",", "n"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", StyleBox[ RowBox[{"x", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"n", "+", "1"}], "}"}]}], "]"}]}], Background->RGBColor[1, 1, 0]], StyleBox[";", Background->RGBColor[1, 1, 0]], StyleBox["\[IndentingNewLine]", Background->RGBColor[1, 1, 0]], StyleBox[ RowBox[{"v", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"n", "+", "1"}], "}"}]}], "]"}]}], Background->RGBColor[1, 1, 0]], StyleBox[";", Background->RGBColor[1, 1, 0]], "\[IndentingNewLine]", "\[IndentingNewLine]", StyleBox[ RowBox[{ RowBox[{"x", "[", RowBox[{"[", "1", "]"}], "]"}], " ", "=", "x0"}], Background->RGBColor[1, 1, 0]], StyleBox[";", Background->RGBColor[1, 1, 0]], StyleBox["\[IndentingNewLine]", Background->RGBColor[1, 1, 0]], StyleBox[ RowBox[{ RowBox[{"v", "[", RowBox[{"[", "1", "]"}], "]"}], "=", "v0"}], Background->RGBColor[1, 1, 0]], StyleBox[";", Background->RGBColor[1, 1, 0]], "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", "<", RowBox[{"n", "+", "1"}]}], ",", " ", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", StyleBox[ RowBox[{ RowBox[{ RowBox[{"x", "[", RowBox[{"[", RowBox[{"i", "+", "1"}], "]"}], "]"}], " ", "=", " ", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], "+", RowBox[{"h", "*", RowBox[{"f1", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"v", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}]}]}], ";", " ", "\[IndentingNewLine]", RowBox[{ RowBox[{"v", "[", RowBox[{"[", RowBox[{"i", "+", "1"}], "]"}], "]"}], " ", "=", " ", RowBox[{ RowBox[{"v", "[", RowBox[{"[", "i", "]"}], "]"}], "+", RowBox[{"h", "*", RowBox[{"f2", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"v", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}]}]}]}], Background->RGBColor[1, 1, 0]]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"Transpose", "[", RowBox[{"{", RowBox[{"t", ",", "x", ",", "v"}], "}"}], "]"}]}], "\[IndentingNewLine]", ")"}]}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{ "\:0422\:043e\:0432\:0430", " ", "\:043c\:043e\:0434\:0438\:0444\:0438\:0446\:0438\:0440\:0430\:043d\ \:0435", " ", "\:043d\:0430", " ", "\:043a\:043e\:0434\:0430", " ", "\:043f\:043e"}], "-", RowBox[{ "\:0433\:043e\:0440\:0435", " ", "\:043d\:0435", " ", "\:0435", " ", "\:0438\:0437\:0432\:044a\:043d\:0440\:0435\:0434\:043d\:043e", " ", "\:0443\:0434\:043e\:0431\:043d\:043e"}]}], ",", " ", RowBox[{ RowBox[{ "\:0437\:0430\:0449\:043e\:0442\:043e", " ", "\:043d\:0430\:043b\:0430\:0433\:0430", " ", "\:0434\:0430", " ", "\:043f\:0440\:0435\:043d\:0430\:043f\:0438\:0441\:0432\:0430\:043c\ \:0435", " ", "\:043a\:043e\:0434\:0430", " ", "\:0437\:0430", " ", "\:0440\:0430\:0437\:043b\:0438\:0447\:0435\:043d", " ", "\:0431\:0440\:043e\:0439", " ", RowBox[{ "\:0443\:0440\:0430\:0432\:043d\:0435\:043d\:0438\:044f", ".", " ", "\:041e\:0441\:0432\:0435\:043d"}], " ", "\:0442\:043e\:0432\:0430", " ", "\:0440\:0435\:0430\:043b\:043d\:043e", " ", "\:043a\:043e\:0434\:044a\:0442", " ", "\:043c\:0430\:0440\:043a\:0438\:0440\:0430\:043d", " ", "\:0432", " ", "\:0436\:044a\:043b\:0442\:043e", " ", "\:043f\:0440\:0430\:0432\:0438", " ", "\:0435\:0434\:043d\:043e", " ", "\:0438", " ", "\:0441\:044a\:0449\:043e", " ", "\:043d\:0435\:0449\:043e", " ", "\:0437\:0430", " ", "x", " ", "\:0438", " ", RowBox[{"v", ".", " ", "\:0422\:043e\:0433\:0430\:0432\:0430"}], " ", "\:0440\:0435\:0430\:043b\:043d\:043e", " ", "y", " ", "\:043c\:043e\:0436\:0435", " ", "\:0434\:0430", " ", "\:0435", " ", "\:0441\:043f\:0438\:0441\:044a\:043a", " ", "\:043e\:0442", " ", "\:0432\:0438\:0434\:0430", " ", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"_", ",", "_"}], "}"}], ",", RowBox[{"{", RowBox[{"_", ",", "_"}], "}"}], ",", "...", ",", RowBox[{"{", RowBox[{"_", ",", "_"}], "}"}]}], "}"}], ".", " ", "\:0422\:043e\:0433\:0430\:0432\:0430"}], " ", "f", RowBox[{"(", RowBox[{"t", ",", " ", "u"}], ")"}]}], "=", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"f1", RowBox[{"(", RowBox[{"t", ",", "u"}], ")"}]}], ",", " ", RowBox[{"f2", RowBox[{"(", RowBox[{"t", ",", "u"}], ")"}]}]}], ")"}], "T"]}]}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"f", "[", RowBox[{"x_", ",", "u_"}], "]"}], ":=", RowBox[{"{", RowBox[{ RowBox[{"u", "[", RowBox[{"[", "2", "]"}], "]"}], ",", "9.8"}], "}"}]}], " "}], "\n", RowBox[{ RowBox[{"explicitEuler", "[", RowBox[{"f_", ",", "h_", ",", "t0_", ",", "T_", ",", "u0_"}], "]"}], ":=", " ", RowBox[{"(", "\[IndentingNewLine]", RowBox[{ RowBox[{"n", "=", RowBox[{"Ceiling", "[", RowBox[{ RowBox[{"(", RowBox[{"T", "-", "t0"}], ")"}], "/", "h"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"t", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"t0", "+", RowBox[{"i", "*", "h"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "0", ",", "n"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"y", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"n", "+", "1"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "1", "]"}], "]"}], " ", "=", "u0"}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", "<", RowBox[{"n", "+", "1"}]}], ",", " ", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "[", RowBox[{"[", RowBox[{"i", "+", "1"}], "]"}], "]"}], " ", "=", " ", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}], "+", RowBox[{"h", "*", RowBox[{"f", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}]}]}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"Transpose", "[", RowBox[{"{", RowBox[{"t", ",", "y"}], "}"}], "]"}]}], "\[IndentingNewLine]", ")"}]}]}], "Input",ExpressionUUID->"1dc66ca4-ec0e-4a3f-88c4-d6723b1da711"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"apprPos2", "=", " ", RowBox[{ RowBox[{"explicitEuler", "[", RowBox[{"f", ",", "h", ",", "0", ",", "10", ",", RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}]}], "]"}], "[", RowBox[{"[", RowBox[{"All", ",", " ", "2", ",", " ", "1"}], "]"}], "]"}]}], ";"}], "\n", RowBox[{"apprPos2", "-", "apprPos1", " ", RowBox[{"(*", RowBox[{ "\:043f\:043e\:043b\:0443\:0447\:0435\:043d\:0438\:0442\:0435", " ", "\:0440\:0435\:0437\:0443\:043b\:0442\:0430\:0442\:0438", " ", "\:0441", " ", "\:0434\:0432\:0430\:0442\:0430", " ", "\:043c\:0435\:0442\:043e\:0434\:0430", " ", "\:0441\:0430", " ", "\:0435\:0434\:043d\:0438", " ", "\:0438", " ", "\:0441\:044a\:0449\:0438"}], "*)"}]}]}]}]], "Input", CellChangeTimes->{{3.752177918609593*^9, 3.7521781792219567`*^9}, { 3.752178260797431*^9, 3.752178316932417*^9}, {3.7521783577131577`*^9, 3.7521784131095643`*^9}, {3.752178834848409*^9, 3.7521788437113714`*^9}, 3.7521794777117834`*^9, {3.7521796243461375`*^9, 3.7521796578090715`*^9}, { 3.7521803712161956`*^9, 3.752180420421115*^9}, 3.7521806271312213`*^9, { 3.7523356272674203`*^9, 3.752335629724021*^9}, {3.752336368471693*^9, 3.7523363936373568`*^9}, {3.7840266139031935`*^9, 3.784026639667667*^9}, { 3.7840266812010427`*^9, 3.7840267314429164`*^9}, {3.7840267732203054`*^9, 3.784026816056756*^9}, {3.784026863300458*^9, 3.784027041368643*^9}, { 3.784038437305148*^9, 3.7840384448280315`*^9}, {3.7840389733229647`*^9, 3.784039062167674*^9}, {3.7840390937434616`*^9, 3.78403911226377*^9}, { 3.7840392036559906`*^9, 3.7840394424901104`*^9}, {3.7840394864267445`*^9, 3.784039490196568*^9}, {3.7840395438525925`*^9, 3.784039687572624*^9}, { 3.784039857678109*^9, 3.7840400361324606`*^9}, {3.7840400670740833`*^9, 3.7840400678723507`*^9}, {3.7840401040679193`*^9, 3.7840401987253733`*^9}, {3.784040255103859*^9, 3.784040329260044*^9}, { 3.784040361455741*^9, 3.784040417028134*^9}, {3.7840404540127726`*^9, 3.7840405279705815`*^9}},ExpressionUUID->"43927774-48b4-4b09-b3a5-\ b0e0d403cb3b"], Cell[BoxData[ RowBox[{"{", RowBox[{ "0", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`", ",", "0.`"}], "}"}]], "Output", CellChangeTimes->{{3.7840268958803215`*^9, 3.7840270421136856`*^9}, 3.784038339030588*^9, 3.7840389827985096`*^9, 3.784040336833477*^9, { 3.7840403985853567`*^9, 3.7840404190064144`*^9}, {3.784040456079209*^9, 3.7840404914026527`*^9}, 3.784040540097307*^9, 3.784040960694043*^9, 3.784041003923068*^9, 3.7840410919392567`*^9, 3.7840411932877464`*^9},ExpressionUUID->"5e613eb3-9ff3-474f-8465-\ 8992faaac46a"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Animate", "[", RowBox[{ RowBox[{"Graphics", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", RowBox[{"-", RowBox[{"exact1", "[", RowBox[{"t", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}]}], "}"}], ",", "5"}], "]"}], ",", " ", RowBox[{"{", RowBox[{"Red", ",", "Dashed", ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", RowBox[{"-", RowBox[{"apprPos2", "[", RowBox[{"[", "i", "]"}], "]"}]}]}], "}"}], ",", "5"}], "]"}]}], "}"}]}], "}"}], ",", "\[IndentingNewLine]", "\t\t\t\t\t\t\t ", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "10"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "550"}], ",", "0"}], "}"}]}], "}"}]}], ",", " ", RowBox[{"Axes", "\[Rule]", RowBox[{"{", RowBox[{"False", ",", "True"}], "}"}]}]}], "]"}], ",", "\[IndentingNewLine]", "\t\t ", RowBox[{"{", RowBox[{"i", ",", "1", ",", RowBox[{"Length", "[", "t", "]"}], ",", "1"}], "}"}]}], "]"}]], "Input", CellChangeTimes->{3.7840405003121386`*^9, 3.7840406869122343`*^9},ExpressionUUID->"a27a273f-311e-4fc4-94f2-\ 2a24d322b436"], Cell[BoxData[ TagBox[ StyleBox[ DynamicModuleBox[{$CellContext`i$$ = 692, Typeset`show$$ = True, Typeset`bookmarkList$$ = {}, Typeset`bookmarkMode$$ = "Menu", Typeset`animator$$, Typeset`animvar$$ = 1, Typeset`name$$ = "\"untitled\"", Typeset`specs$$ = {{ Hold[$CellContext`i$$], 1, 1001, 1}}, Typeset`size$$ = { 36., {213., 219.}}, Typeset`update$$ = 0, Typeset`initDone$$, Typeset`skipInitDone$$ = True, $CellContext`i$4068$$ = 0}, DynamicBox[Manipulate`ManipulateBoxes[ 1, StandardForm, "Variables" :> {$CellContext`i$$ = 1}, "ControllerVariables" :> { Hold[$CellContext`i$$, $CellContext`i$4068$$, 0]}, "OtherVariables" :> { Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$, Typeset`animator$$, Typeset`animvar$$, Typeset`name$$, Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$, Typeset`skipInitDone$$}, "Body" :> Graphics[{ Circle[{0, -$CellContext`exact1[ Part[$CellContext`t, $CellContext`i$$]]}, 5], {Red, Dashed, Circle[{0, -Part[$CellContext`apprPos2, $CellContext`i$$]}, 5]}}, PlotRange -> {{0, 10}, {-550, 0}}, Axes -> {False, True}], "Specifications" :> {{$CellContext`i$$, 1, 1001, 1, AppearanceElements -> { "ProgressSlider", "PlayPauseButton", "FasterSlowerButtons", "DirectionButton"}}}, "Options" :> { ControlType -> Animator, AppearanceElements -> None, DefaultBaseStyle -> "Animate", DefaultLabelStyle -> "AnimateLabel", SynchronousUpdating -> True, ShrinkingDelay -> 10.}, "DefaultOptions" :> {}], ImageSizeCache->{332., {252., 258.}}, SingleEvaluation->True], Deinitialization:>None, DynamicModuleValues:>{}, SynchronousInitialization->True, UndoTrackedVariables:>{Typeset`show$$, Typeset`bookmarkMode$$}, UnsavedVariables:>{Typeset`initDone$$}, UntrackedVariables:>{Typeset`size$$}], "Animate", Deployed->True, StripOnInput->False], Manipulate`InterpretManipulate[1]]], "Output", CellChangeTimes->{3.7840405030418286`*^9, 3.784040540180274*^9, 3.784040687978239*^9, 3.7840409607350187`*^9, 3.7840410039670525`*^9, 3.7840410919832306`*^9, 3.78404119331673*^9},ExpressionUUID->"1b958677-d76f-4868-96fb-90eef5c203a6"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[TextData[{ "\:0420\:0435\:0448\:0435\:0442\:0435 \:0437\:0430\:0434\:0430\:0447\:0430\ \:0442\:0430\n", Cell[BoxData[ FormBox[ RowBox[{ RowBox[{ RowBox[{"r", "'"}], RowBox[{"(", "s", ")"}]}], "=", RowBox[{"cos", "(", "phi"}]}], TraditionalForm]],ExpressionUUID-> "580175f4-fb8f-44ca-a579-cfd83ee703b1"], "), \n", Cell[BoxData[ FormBox[ RowBox[{ RowBox[{ RowBox[{"z", "'"}], RowBox[{"(", "s", ")"}]}], "=", RowBox[{"sin", "(", "phi", ")"}]}], TraditionalForm]],ExpressionUUID-> "580175f4-fb8f-44ca-a579-cfd83ee703b1"], ", \n", Cell[BoxData[ FormBox[ RowBox[{ RowBox[{ RowBox[{"phi", "'"}], RowBox[{"(", "s", ")"}]}], "=", RowBox[{ RowBox[{"2", "b"}], "+", " ", RowBox[{"c", " ", "z"}], " ", "-", FractionBox[ RowBox[{"sin", "(", "phi", ")"}], "r"]}]}], TraditionalForm]], ExpressionUUID->"580175f4-fb8f-44ca-a579-cfd83ee703b1"], ", \nr(0) =z(0)=phi(0)= 0,\n\:0438\:0437\:043f\:043e\:043b\:0437\:0432\:0430\ \:0439\:043a\:0438 \:0410\:0434\:0430\:043c\:0441-\:0411\:0430\:0448\:0444\ \:043e\:0440\:0442. \:0420\:0435\:0448\:0430\:0432\:0430\:0439\:0442\:0435, \ \:0434\:043e\:043a\:0430\:0442\:043e r<1" }], "Title", CellChangeTimes->{{3.747918045893719*^9, 3.7479180792016163`*^9}, { 3.748019567260964*^9, 3.7480195984310884`*^9}, {3.7491799306613116`*^9, 3.749179945330095*^9}, {3.7491820390416107`*^9, 3.749182140816432*^9}, { 3.7491822953940344`*^9, 3.7491823155841894`*^9}, {3.749182378833807*^9, 3.749182424987447*^9}, {3.749182543546228*^9, 3.74918254777547*^9}, { 3.7491855618392*^9, 3.749185791014152*^9}, 3.7494816355928173`*^9, { 3.7494826582226067`*^9, 3.749482678435937*^9}, {3.749483068682683*^9, 3.7494830721446886`*^9}, {3.7494884082644987`*^9, 3.7494884424507885`*^9}, {3.749488494590744*^9, 3.7494885492422333`*^9}, { 3.749488662765763*^9, 3.7494886848150673`*^9}, {3.7494887168485804`*^9, 3.749488740273077*^9}, 3.749490389264352*^9, {3.749789440131727*^9, 3.749789499828309*^9}, {3.74979010532883*^9, 3.7497901250714474`*^9}, 3.7499847604043837`*^9, {3.7499852544136286`*^9, 3.749985272769679*^9}, { 3.7515612145185156`*^9, 3.751561319842022*^9}, {3.751561379568722*^9, 3.751561400725582*^9}, {3.784022901846876*^9, 3.784022907924224*^9}, { 3.7840231556633935`*^9, 3.7840232781353984`*^9}, {3.78402332851328*^9, 3.7840233557128353`*^9}}, FontFamily-> "Times New Roman",ExpressionUUID->"4ed04c9b-35d8-4b0c-ad3b-76a639522d18"], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"b", "=", RowBox[{"1", "/", "0.5424"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"c", "=", RowBox[{"-", "2.9"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"f", "[", RowBox[{"s_", ",", "u_"}], "]"}], ":=", RowBox[{"{", "\[IndentingNewLine]", RowBox[{ RowBox[{"Cos", "[", RowBox[{"u", "[", RowBox[{"[", "3", "]"}], "]"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sin", "[", RowBox[{"u", "[", RowBox[{"[", "3", "]"}], "]"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "1", "]"}], "]"}], "\[NotEqual]", "0"}], ",", RowBox[{ RowBox[{"2", "b"}], "+", RowBox[{"c", "*", RowBox[{"u", "[", RowBox[{"[", "2", "]"}], "]"}]}], "-", RowBox[{ RowBox[{"Sin", "[", RowBox[{"u", "[", RowBox[{"[", "3", "]"}], "]"}], "]"}], "/", RowBox[{"u", "[", RowBox[{"[", "1", "]"}], "]"}]}]}], ",", "b"}], "]"}]}], "\[IndentingNewLine]", "}"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"res", "=", RowBox[{"AB4", "[", RowBox[{"f", ",", "0.01", ",", "0", ",", "5", ",", RowBox[{"{", RowBox[{"0", ",", "0", ",", "0"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"res", "[", RowBox[{"[", RowBox[{"i", ",", "2", ",", "1"}], "]"}], "]"}], ",", RowBox[{"res", "[", RowBox[{"[", RowBox[{"i", ",", "2", ",", "2"}], "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", RowBox[{"Length", "[", "res", "]"}]}], "}"}]}], "]"}], ",", " ", RowBox[{"AspectRatio", "\[Rule]", "Automatic"}]}], "]"}]}], "Input", CellChangeTimes->{{3.784023375609974*^9, 3.7840233874306498`*^9}, { 3.784024825852923*^9, 3.784024887992477*^9}, {3.7840249907113523`*^9, 3.7840251951860476`*^9}, {3.7840252731995096`*^9, 3.784025286265257*^9}, { 3.7840254005237923`*^9, 3.784025416689717*^9}, 3.784025596036975*^9, { 3.784026487053938*^9, 3.784026489222062*^9}},ExpressionUUID->"52f186ac-f2de-44ba-a065-\ 046acfe79cb1"], Cell[BoxData[ GraphicsBox[{{}, {{}, {RGBColor[0.368417, 0.506779, 0.709798], PointSize[ 0.0055000000000000005`], AbsoluteThickness[1.6], PointBox[CompressedData[" 1:eJxdV3dczf/3b+8dEhqiSEKiEnm+tIhKRWbSENIgipasQkRDZYQWZSSVVEKJ tD+0l/be997q1r239b2/f3/3n/t4Ps7rPM/7cV7nec7rrHQ4b+XExcHBQefk 4Pi/////S78/QKJ8AyF09PFYv+pScOg3aqX7RGLkmfyx5ao66OCtX9if+RKZ Y9Itr9fvw0Ro3mSSfRJMORouf1Q+gpNef0yI73vk5Gr3flY5BaHpuZ6ZPWko fpdnJnDYFWebKxQXqjOgaEFnZHhcwvee4ACf45/BVyDw7ayND26lqz5R6MrG 2jhZ3lTO65ggGvc/n85FgraE7dLMWziS/+m00OA35GnUhQ3+uoOjmdHxvWfz cXNRl5qeRgjsrnpFCU38gPxy27C7V8NwOEF8+eZNP6FTPlvmLfEIMfQF10qX X9hsIKPMXxWFkMPyNypeF8JTTO5TU94TrGLckhBr+w3jd5bRd5pikPVr99ar i4oxetnuS5JiLCp3+CRwmZTA4sXCfJdJPN6pZT2641eKsy7mTG3ZRKhbhbbQ 35ex43XPLZp/BRetuUHt5nJYw8FLl5GEaw1SpTp8/+G0pGzQaf63mLu7Ta1n 4x8s83qy0VrtPSKFh86pHP4Lb55y1x+OH1Cb1Lt7Sq8SOSnPjezSPuKgiKmq 1N9K/Ji/c8VLIB0sWWuJNtsq7Fafv3hPPwNayc9CLUerIPMjIGjP7U+Y15V5 fMq3Gr1aFifv1mTibOIFL3G+GvgsaZ5fuS4Lj4Jkm/aG1SD6R0PeSHA2qMZ7 q8Rla5G/8sGp/PEcrPhV9edwXC3GF9PVL5zKhYpt3qsVKnXovN/ZXtv2FS3L O+4eeVcH8w+8tT/tvyPNf/UiPvV6KEd28YmP5mHldIW+Ymo92lKKpXaTH/Bg XvdMUG9AS/k7MyGxAtRf3OF6930D2tof7P/aXgAe+nKLP2saMXdtn6Lh559Q C+V/eyGhEetPxuY+fPgLnBF6Dq7LmzDIEyAV5FIIe5+8kryIJmS37F4uafob hwxi7jkJNOPF6rr/Vmwqwsrba7Js/ZvRdf7IhgiZYnS8lPjyZqwZLnGLrY5z lUDs6PGT+if/4erp4VMO1BKs8jSpXvfnH17e+pAa3lkKhQtrH9tvb4GyrvmT mroyDAsLP+xLaoHWtbo0qf/KYdzpTv0i0YozaaszthVXYEfw3/0t3q2YzwkK 0iz8D2cz1kaat7fi7DeePcO//kAs9CxL0rANSwIfmZsW/QXHXI3U2uQ2rBGe WRTvVInGXyfyIwTakcuxw9pmsBL/MK5x8Gw7OPk+08zcqxBYkjPkXtSOrWJb ZFwnq9AiovLvomQHGpM0ZT74VaMygF7XurMDBUrGzvw8NajYWrXY16UD7gun Ujwe1EB91TZ3k8cdOJey+k2XTC0OjmjVmf3swNKhCjPLxFpsjnB1CB7pwC/O Q18yN9ZBik/JfH5xJ6JW79vM870OUapW1NSdnRDaG7Fee1897qVn57843YnL VFkTo+Z6DNk0W/z3oBNBam6HVM414OfQCmuS2Qlerw2W1awGbNeseMps6gT1 q8+pXSGNcCjYyD8734kXWzYMucg34Znh89i9q7qw5pKiy970JoiuuxswYNyF W7ZbY+sMm1GRnNDQ6NyF2NDTbnxNzcjzH6XLh3Rh8lPacJ3bP8zFN+nmfehC +MMRWz3uFlSsOiOX+6cLkVsf6Wo9bYGg1R7VJZQu+K53V8zd2IpcD4H5OrFu JM9pdvwsaoWRtJb5nHo36FOSV/batmF6XXhbqGk3dqSuEt831YbX3HpB4ee6 4TTROJ/3sB2JclNr+e52Y/ac8qlK6Q6U8Zn+oL3qhuRmEuJl1wFJzXaLIwXd 8O6xdDT90AG5xoQxvdZu+GlHOluxOnCoVOlrEqMbAqc9NIN2d8LO7L+eaOke zOwQdO6K7IRtJ0mS2tCDbbEff5/p6oRW8wsj+T09CKg+JSO9qQv0v/+k8+x7 oNa7e2fX1S7kfVVynfLtQaU/9/qa8i5IfFjzqfxRD4JttnZ3y3ajdzrkkH5K D37SKi9Lne2GdK5RhW1hDyy3XhU8kdWNoEez4cotPTDdMFZRwNMDcu7q+piJ Htxz1OXEgR40SyitKRTqhcbd9ZPV8T3oSdwl/XplLw59DOfwp/aAeiHhjr5O LyytAh/ooBf6OmMK8ea9WOKVKcr/sBcjT9b6/z7Vi7eaH+Z7W3ox6LWGZPr2 Ipp76UCNWh/4lUSlL4b1YvPk3LG/vn2oVih8yZ3UC3sdflp9aR/O6fuHun3t RbLhl32DS/shMrv2UWZlL/ZVfJLmPduP48kcR1t7e/FgYGR6bXY/2p/Mho2w esH1eSTlAN8ABouODveJ98E7RWDklvUAmPWyS6tW94EuGXsx59UAhOwDiz5u 60NuuhUndWIA93TfvLpj3gdP493HVQ0GEf68xfaEYx+W/LU1cowYRD5rV+xW 7z6sVfIOet45iCsO+qKLH/QhTPhjZ82mIQTskDTjSOiDmv5LBYHrQ8g0W7cw k9WHdzcTuXX+DsGj/m++SEUfrORULBzkhyH36Y6KdicbG3YUBrkNo9bx8e9r U30wPasmH/9tGCd7AlSHhfvxkLp/0SfhEZjsuth4dWU/Nj9vPJZ7bARpzVWW O7T7sdA4GPv57Qh2mMSLqpn1AzYv7icyRpDZOutn7tgPGROH/Ju7R5GT/YDn vU8/eDY61FlGj+KA5yyvUVg/5lSHT4r3juK0j0OdYnI//FzbB75ojmFV/5fv +nn9kPJdR9t/k431TBdS6/rx6vAO6t/KMdAKFf46j/ZDYfLeAR0FCtQbHQ9d 4R1AgtCygDtuFDxa0ZL+T24Ap74Jlnz9SoFb7V7aM60B8IgvyR/koYLu1Lzk y/4BRDMqL/SpUaFAbi/RdR5A/8l86RorKoR236pbeWsAXUa631J9qGDSk1U9 XwzANc7A2zuOimVnaxu35AxAQ+jU/o3FVDz1P/bkTM0AQvxD9KtHqQgyOy0i QRlAyerZfXaLaEjXPFOsITyIMsd1Jxt1aTCI9DGpXDOIEZ7PZ7fZ01CdVbUw ZjgI5ywr85t3aPjds8Yw3GEQDuWlY+kfaLBYpGKTe30QTNPhlcU1NDx1aPB2 ix1E8ImInB9MGqLiZbM+5A3CfyzZ/4XCODJ2Si4LaBtEYtE6DRujcRQXx0Z2 zrPttdsjZs6N49mB1vlmhSHo/N5g7Rs2Dj7pqGVuu4bQUP9ybcPncUjt/x3+ zHEImk13MiX+jWOPzKopp9tDOLXf8ZkyxwSOdMX/qHo7hLtxJl6LlCdw0X7B ouO/IWQYyjT/M5nAVOGB0UfjQ7jRGHXY152NU5MKxmTYdbbw5vJYxATeiQa3 j+sNY2QuIFcrewIeizP1Ek4N415/YZzlvwlsCv2QPX9/GO8t2+L1FibASvEQ FcocxsAJRQ2m0iT+PlLq/t0yDBa3xPebxpNYUvT9hybfCFwuJ35sdJ6E3OMG SetNIzC+dtudGTKJsAobZbXjI3BOsz/e/XESdqHH1LNvj0Bdf3VtZPUkZMoq drIyRrDa10JThD4JBfFLlvT2EUhd2lu5W4aOWbkUvBUdRaj82wWjbXSILj9S tGTHKM4eCZbmPc62BxcE7nQZxanuruO3/OmIk3kwqxgzihGVdpX8F3R07kgL zysfBSNQoyArj46xH+73ZOZGEavgcflsOx3C8ifF128cw4/eS+HVC3Rol3qp zdqPISDvnc2EwhRWHpOpuBM1Bp/3AUqlmILEQ42CP6VjUMxNVbc6OYVeuePf K+fH0Po8tTQkYAo/bVjWD7ZQcOA3r8GVF1P4EHGRm8eFgvtNqf3i36awkefI yi0JFJy2vcJv3TyFyah3YyuaKTh38eSkIWMK6f7CxbaCVDzzZHK1LZ7GRxWF W0IqVPy+VHpPXnMaFmJcKyUNqHA3eNXHYzGNLxYeAhftqNiwSMwr1HUad54l PFofQIWG+5Mn3+5Oo+iMpzh5TsVmN+TcezWNVr34r+9zqVgpXwRG/jRiKull 7k1U2LZLxnH/m8YfrnuOdxlU4HXg/vf0aZzMmH05LUPDId5DZaPiDAAZ975q 0/B5eCiiVJUB5tfFa2sO06DRprZZz4ABjawWVz1vGk5lNGnvtmGgJdDMgPmE BkpMh26fJwPnBO/e482lwcznbdPSBwzkVlYusf/H9jd2zm19xYDE9bZ0oTka dsW439f8xkCvj+pyTrausla5ccvWMFD/Wk/AYNc4NMaLOkMGGbhjLKb813Ec zS9rxh4sMFBevG/Rq9vjSJZVnl22mIk296/nf7wdx/kVZT2b1zEhXB9WovTf OAxUHz6o3clE8IrS5FLqOA5xaffOWzExnHHn4adFE3gYnN+TdpoJynvPqTad CZRMhYYO+jDB67PCau+JCcQf/DSRHMKEHdNCdPLGBKK3CK0bfskEreD+g+ak CWwJ8T2UkcaEWqDl5tnyCSx8WPFmtoCJae9o68O0CRRbOh0trWaiJXOV6dCS SSyNN02T6WaiyWTO6dOOSdzUuzjfN87EyuWlY2kOk4jS3hmvxcWCr7HWhs67 k7gyNs/JJ8lCHJ/TWX22rlSruKMPKrBQTzNvq66bxLeZzLaV6iycurLxT9js JBrfJzu76bKQGDByyXcVHcbBHyu1drMwJnhN8OFeOix+jHy5eoCFnPzjxaUe dKwcrW/YdZKF1u3/etSe0hHSaDh48xwLhsIJTzJ+0DHI7/GFeLHQoSkgazPA 1tGqGQW/ayyoql1/riIxhadb+jg2B7NgIjd3VEJnCufKHpk6RbDgeOZIsIzd FBpmvy4simEhek2Wo+7dKeg6OazalcjC898Dq66kTWGDaEnF4DsWjh6u4ato nMJp+2Ny4hkszBSW7dPmnIbullj99Bx2/DD6lmzVaZi1yl2uzGNhdM07/r1W 05gMNJpzLWThdWy5yJjvNFa/t1C9U8rCdPHq8PjEacyzlA8r/GHh0pWCXoeK aUhd1aJsrWbBkmlkt4ld92+2H3Uqr2PhRZDhcQF5Bk66vdLpaGThbUqY1qAx A77bDSc8/7HQHxt2uOo8u+5196+718rCVoMQ1fwnDGx79+SqYjsLBeNMufQC BhJirj/W6WAhYagk5fUQA31CP9fXsLH3zrcHY6SZqJLELJWNL1tYxoTtYOKh cfmfCDY+m/NDNNCJiTzHyctf2Hy+IdeJ50MmVliXVjq1sVA03/H3ZDb7fGFs wZMWFmoMttgZdTDBY9iofaCZhZBsvperBFmo3RS/LrqBBVvtgzuZGizccst5 6lDLtr/NnS88xoKM+PoHXyrZ+Xv39nXQLRaMly9b/LyChdQToVy6KSwUep43 ES5h8x/kGutg+3FLJ5kI/GIh4umAic8cO08JN3dEfmehoiSGk0tlBpffV+1P y2ZBcY+QkK/5DE4fXZ1zIp2Fw4/q7Dsuz4Cb50HaC/b9cr5yk9wSO4Ps8OO+ Huz7P3hZVc2zeAajVU02zez6cP1oUx5LmcG7Qf7XTY9YUDl3bmmWzCwOjx55 feE+C72b6vSzMAt6wfnPcTdZuKlmGhR7ZhaWJ1s0XHxYeBf8ebVH6CwWhzT5 VZ1noXOz1oV12bMYUBNm1jqx67VIq6CsbRazWnrSV46zsGvSyc+Sbw4JHqt0 vliwoERyRPLU5/CCKdv03IiF1ZudRsWt5yCZ13duFVtPw0Jdscb+c5ih+zmZ b2Dnh7s5zT5xDtpvvhooKLGwWW11tl0Z+/zhHbceL2Zh/3DWAQPaHK51iL3K EWDh+L6nGoJL58GZGLQ2cIYJPVYVb8bOeWybmo2cGWXCUvNgqJ7TPDR5X59X ZN/riwdD0e/vz8NWxoxMVjGhVaRAmUmfh+zBTWHev5joiFb23NA4jyzhwsHU TCYE330RNpifx9/QO41PXjMR8F3g9bbVC3ifWDegHc3Edt4wHam9C/ALyU5+ cJuJP+2zxRXnF6Bc8iUh5jIThaNCts5RC9hz4ecZR3Y/DH1hJzaQu4Ct5oOB fdbs79U7R9/TsYCbiZZ3VYyYSKqSXVHExUFob3k5lLew6xqKXEcXcxCPihqf biUmGjPSd7eu4SDTT9VvOUoyoVFZKnxAl4MoXBR4m8Tu31uodh65phwkeCj/ SuooA4FD1rHSJzmI+DDT1/8fA19PHv5h68FBRANC9KRKGZC7W8H1/BYH0bB0 MvHIYoA/9PuT8igO4tx1cU1MIgMPb37NHk3mIKUHQk+GhzEQ2Vb8hCuXzd/h 6m19lQGf3O3eQhUc5MLIrekBZwb2xqVG8bax/ZO47M0OMdhztVVzksJBoJWn HajPQEj1wyf1nJxkzFV57uEGBgzK/UU+SHMSq7EtGueXseN3aP27osxJ3Mtc rVT4GKh2ETygrc1JGH7FrZ9p01j6tKtweA8nOborw1yhdRo6IdqBUcc4SVUh n8Kpkmn4nxns2+LKSZj8fMVBn9jzXX2dUdlVTuJ1++D3oJfT8Lprz7QO5SS7 bk1fdgqehvllSbvGOE6SfOnAGRXPaSzvy22wzOAk4b5qHGW201hj8PPHz1+c JKEvPt3SZBofpg5cVavjJLqvtw3lsd8TD7v2eIT0cZIdzxuXSMtPQ7lJiKtv mpOgpSxhv8A0Ghf7ROsIcpHY50XTXuNTmP9SER24jIukjlOv3mmZwhXLq65l alxk3dbumKCiKex1UL4ppMdFGLNpWR7sPu4kdMbAyJyL8GWOaZg9m8KL1MoV fie5iFZiznnZwCnUfi/wS7nARcQv5FY2uE1h6/39f5pucBGeGa3M4MNT0Ndf d4v7ERcRDJXw3rRrCgHTNwTXveIi5geqvP9bNwWrxuFKs89cZJDySdJh0RTq /XebnC/iIk6HUp5T5uiQsn3V9bCBi2yKvXDxSj977l1aI/phgItwj94YYlXS UcKg6pQxuUjo7q3H/XLp4HDkqu8T4iZ1tt9kWYl08F9ffZB7BTf5sXaDr/cD Om6W/aekqM5NaveG/p2+TAehCqTr7eQm47uaXXzs6NAJPe5ms5+bXNPlj5s3 Yc9NtQ+//e24idd96W93Ndl2+w1rYj24iVqxgLisHB2bD5L5Xze5Cbc1d38a Hx2f7W2yhh9xE/c7OuGW1Em8XDz7ZclrbrJUafbobBN7rp++cdEoi5s4G/4O yfjFfjfoKezxLuYmUymbPb0+TELP6Gf8x0ZuUrNX7Jjx40l8lGb1DA9yE20p izsqNyZRGtPlpT7DTSyFFA1lXCaxf04021OEh0Tda2iQtZ6Ecoq74E85HuL0 kjN4IyaRQN9VsXgjDzl4+9iLo6qTkO1n2VwgPKTP/M/Rp9KTMHqqLlNtyUM0 vC5zUecmcKjH9+B2Rx6yKyK62XZgAv7X/V1TPHnIX74nq/urJ9AanZ6pfJuH 2MfzKN37PgFakcv15Mc85Oa269ImbyZg9+qb4ua3PMR1e4GO8qMJBKU79vzO 5SHSincrFQImUNr8Tcaxgoeo/4laqus8AbkrH3mE2niIX2ecmdfBCXRXvRjP pfAQ3XqT99WYQIy+tcplTl5ySkbM9YDaBH7Wz3dul+Ylc2PbShlLJrA9QsJJ RJmXXDzSxPrFNYG8vRjp1+Il+Z+TtbPGxpGk5JH+Zw8v+XZ2PrOyeRyPk9OG fxzjJalWP5MWFY/D/Kni3zxXXjK9wwK3Po1jwYYaUhLAS85tOPtRLm4cwSst D7SF8ZLb6TIqvSHjyM1xMeNI5CVPJRIrm33G8Uf1cOKGz7zE9YFb88Lpcehy Svg6F/OSwWcdrtYHx/HF/sJEWhMv+U/B6mML+93r2bx2A+8IL1EuGyyK3jiO Cd5co9PzvGSxf1B9oNw4LtHdrWol+AjPv+fzycLjYLgccd6/io+MCq0+wcmi 4ejInucNW/mIYZybaMQADaYvqpiue/jIpRrhtccbaNDiIw/Ej/ORh5zpP08U 0eCQ/Nw2342PuN8Wpz77TAPt/LuL/tf5CGs2J1/6NQ2PpPvLjB/xkdbxWb3S SBqGd4U4ySfxkYtyxZe+BdKwoiBcj/sLH1GXv3p6zJOGIte0vfRyPqLU+0/G zom9Jyw7fWOyjY/U/KgNkjpEwwdHVjXHOB9JX3fuHd9uGowOB6ov4+Unq+1C /PR02HuHOM2PLOUnWxvpQ9mq7D1ik3DSJTV+giT+vsvLadi0+1bYp538RKv5 tbWfKA2e0swVHFb85DBFWr5ogYruD7zrjjrxE9VjVfJHxqkIjDB+lOfNTxqM L2/e3EtF9l1njU0h/GT+SdgG60YqDpzZ1vwhlp8s6N3oLiinojD9vZ32J34y 0XV0hU8+FU1ClMiKIn6yTeHrS79PVFzz8TV0a+YnN1Z/WleWTIWFn5yu7Bg/ keR94evE3sMud0jsrOQUIN/jOAzMw6kQdBVYErFYgHTvubQ56DYV7zSNgmxV BchBNR+akD8V+TlvTLfqCZCGR2eUOjyocP1csm6JpQAJkN5jz3uWCko0dwuH kwB58/ybgZ8tFSUvo5bRvQWIcELPEX1rKr4npr4cDxEgCcYGqsdNqTjvclee GSdA2t8naRWx98iZyKvmAp/Z/OckZm9tp6JCZGWnQqkAiTCMHY3UpOI256E7 aGXzLfb1ZqpR8dbyLPU0TYAkS3/jylhNxQln5bIoXkFyxMVhOE+OirA79B/l soLkt+b2e4oyVLjZ068JbBAkGrV1GxskqGAsk0jfpy9Iphs8lbuFqPC5LzsX eUiQfNnO+cGQl4odhQPCPecEycpL85yMBQrEKi8G6lwTJDaagm+lpij4Hd8r EvFIkMgstU9qGqFAYu2vvdRkNp/cXn2zHgry/Ge7rL4Jkpky7iy3FgoSlEV9 v1QKkuadP9W211Gw88SSzFW9gsTP+V3e5z8UWB1lbQlnChLDs1W+rSUUhAdd ieMUEyKagjiR8YuClDcGzzyVhIisx8SBLfkUSA2pFA1qCZGwDoqe41cK/qQu KbHbJ0SSwvWHd+ZQsMjQwrDppBCJSlosV/yZAhXzS+1WnkLEWLEgaOYTBe7j MTv/uytE/lS1l/zLoKDpkz//nhdC5LHhGb8zbJzAuNLxM12I2F/7ueolG3Nr fLmgVyRERuMo266x/X2izeyym4VIhl7cGiE2/zW/HUc2UYRIkbvNnp3ZFEys iuV6wy1MSgLeDC7LpSBSPkdIfqkwuaGovi/2OwVHlwrJRawXJmWyXS/LCih4 q+Tez71LmGiZbhd8XkRBrZOq7CVrYcLLcyReqoICt21cju3OwkTAa4uHejUF Uc8PXtkTIEy2m7Z6DTey4wuoCKVGCBO5Ha9j93VQMPxYp0o8WZisdbapNh2g IDRALtDtqzBRiwnuHaVS0F1u3FL0V5gUyrx+pc6igM9radTyHmHSHfjuDC8X FdXxDBcXhjC5nlInos9PxattJxWzRERI8bnNSypEqIhL6bk0oyhC/lZEVidI UfFmav+67VtFyB79M/d/L6UipClR0MtEhEg1lQdqKlDhYKrW/OaECDEa+D43 rEyFvnzQ6XoPEbJqxZZNY+up4CxP9ZsLEiGutfqW27ZQ0XpdbFDumQipKdoZ XsHWw+lD/te1U0XI7ob25clsveil7V2y96cIWbP01NKSfVRE1ZvdsK4XIXd6 Mt5uOEhFXm7Xs6NDImSB7/pYkw37/OjkzkPzIqSw3Z2z1IndT35zm+6TEiXS V7cyWe5UiH5VSdRRESX7nuYPuHtTYUqXWCuvK0ouX7RsX3WTChfRRQUzZqLk 8Pl3A0tDqLggpHS4yl6UBNxeI2oWTQUr6k/dSy9RcsSm1zI/jordbvvVTgWL knDRrhz392y9SfjvWflClLyIkd91LIuK41kVivVpoiSdlt93vYAd7+qRdzcL Rcmc98bEzgp2/kerW1UaRYmd5nY3P3b/a/1xN+vnsChxPrh1m1kPFSnTxzcf WhAlFAtppjWVihtG+ZadUmJE0jbkeeQsFd5/02ROqYiR+SPakgKC7P5unBnc tk2MOEi+0f+4mIYy24qPFmZiZIB6f1mwEg1i7y9F5tqJkf9Yv69HbqRhUUuk 4QpPMTInf+945Q4aXrne/Ot1R4yUtlwJ3b6XhpGvo8bFz8TIH43MhdrDNITE tOZJpIqRldVcD2PY8+Sj9SszqwIx4prctuzhJRpSNtYJh9SKEYWo8KD3N2jI d9u2KK9fjLhc35cyHkqDS0vhrQGWGFGJ0XF1eEnD8p5gH2ExcWJDi45nptDw RO3DkjUrxUnStvI12V9p+B8aGWU5 "]]}, {}}, {}, {}, {}, {}}, Axes->{True, True}, AxesLabel->{None, None}, AxesOrigin->{0, 0}, DisplayFunction->Identity, Frame->{{False, False}, {False, False}}, FrameLabel->{{None, None}, {None, None}}, FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}}, GridLines->{None, None}, GridLinesStyle->Directive[ GrayLevel[0.5, 0.4]], ImagePadding->All, Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({ (Identity[#]& )[ Part[#, 1]], (Identity[#]& )[ Part[#, 2]]}& ), "CopiedValueFunction" -> ({ (Identity[#]& )[ Part[#, 1]], (Identity[#]& )[ Part[#, 2]]}& )}}, PlotRange->{{0, 3.8911878367993715`}, {0, 1.6799173215268761`}}, PlotRangeClipping->True, PlotRangePadding->{{ Scaled[0.02], Scaled[0.02]}, { Scaled[0.02], Scaled[0.05]}}, Ticks->{Automatic, Automatic}]], "Output", CellChangeTimes->{{3.784025411620427*^9, 3.784025416954732*^9}, 3.7840264925852547`*^9, 3.784038251938941*^9, 3.78404054052127*^9, 3.7840409616034994`*^9, 3.7840410040939703`*^9, 3.78404109210917*^9, 3.7840411934616456`*^9},ExpressionUUID->"87f05e1b-9fc3-4a81-b896-\ 949e6f1f10b9"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{ "\:041c\:0435\:0440\:043d\:0438\:0442\:0435", " ", "\:0435\:0434\:0438\:043d\:0438\:0446\:0438", " ", "\:0432", " ", "\:0437\:0430\:0434\:0430\:0447\:0430\:0442\:0430", " ", "\:0441\:0430", " ", "\:0438\:0437\:0431\:0440\:0430\:043d\:0438", " ", "\:043e\:0431\:0438\:043a\:043d\:043e\:0432\:0435\:043d\:043e", " ", "\:0442\:0430\:043a\:0430"}], ",", " ", RowBox[{ "\:0447\:0435", " ", "\:0440\:0430\:0434\:0438\:0443\:0441\:044a\:0442", " ", "\:043d\:0430", " ", "\:043a\:0430\:043f\:0438\:043b\:044f\:0440\:0430", " ", "\:0435", " ", "1.", " ", "\:0422\:043e\:0432\:0430", " ", "\:043e\:0437\:043d\:0430\:0447\:0430\:0432\:0430"}], ",", " ", RowBox[{ "\:0447\:0435", " ", "\:0442\:0440\:044f\:0431\:0432\:0430", " ", "\:0434\:0430", " ", "\:0440\:0435\:0448\:0430\:0432\:0430\:043c\:0435"}], ",", " ", RowBox[{ RowBox[{"\:0434\:043e\:043a\:0430\:0442\:043e", " ", "r", RowBox[{"(", "s", ")"}]}], "\[LessEqual]", RowBox[{ "1.", " ", "\:0421\:043b\:0435\:0434", " ", "\:0442\:043e\:0432\:0430", " ", "\:0440\:0435\:0437\:0443\:043b\:0442\:0430\:0442\:0438\:0442\:0435\ "}]}], ",", " ", RowBox[{ "\:043a\:043e\:0438\:0442\:043e", " ", "\:0441\:0435", " ", "\:043f\:043e\:043b\:0443\:0447\:0430\:0432\:0430\:0442"}], ",", " ", RowBox[{ "\:043d\:0435", " ", "\:0441\:0430", " ", "\:0444\:0438\:0437\:0438\:0447\:0435\:0441\:043a\:0438", " ", RowBox[{ "\:0441\:043c\:0438\:0441\:043b\:0435\:043d\:0438", ".", " ", "\:0417\:0430"}], " ", "\:0446\:0435\:043b\:0442\:0430", " ", "\:043c\:043e\:0434\:0438\:0444\:0438\:0446\:0438\:0440\:0430\:043c\:0435\ ", " ", "\:0444\:0443\:043d\:043a\:0446\:0438\:044f\:0442\:0430", " ", "AB4"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"AB4", "[", RowBox[{"f_", ",", "h_", ",", "t0_", ",", "u0_"}], "]"}], ":=", RowBox[{"(", "\[IndentingNewLine]", RowBox[{ RowBox[{"t", "=", RowBox[{"{", "t0", "}"}]}], ";", "\[IndentingNewLine]", RowBox[{"y", "=", RowBox[{"{", "u0", "}"}]}], ";", "\[IndentingNewLine]", RowBox[{"fi", "=", RowBox[{"{", RowBox[{"f", "[", RowBox[{"t0", ",", "u0"}], "]"}], "}"}]}], ";", " ", RowBox[{"(*", RowBox[{ RowBox[{ "list", " ", "of", " ", "values", " ", "of", " ", "the", " ", "right"}], "-", RowBox[{"hand", " ", "side", " ", "f"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{ "Compute", " ", "the", " ", "first", " ", "approximate", " ", "values"}], ",", " ", RowBox[{"using", " ", "the", " ", "RK4", " ", "method"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", "\[LessEqual]", "3"}], ",", " ", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"AppendTo", "[", RowBox[{"t", ",", RowBox[{"t0", "+", RowBox[{"i", "*", "h"}]}]}], "]"}], ";", "\[IndentingNewLine]", 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[{"AppendTo", "[", RowBox[{"y", ",", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}], "+", RowBox[{ FractionBox["1", "6"], RowBox[{"(", RowBox[{"k1", "+", RowBox[{"2", "k2"}], "+", RowBox[{"2", "k3"}], "+", "k4"}], ")"}]}]}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"AppendTo", "[", RowBox[{"fi", ",", " ", RowBox[{"f", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", "+", "1"}], "]"}], "]"}], ",", RowBox[{"y", "[", RowBox[{"[", RowBox[{"i", "+", "1"}], "]"}], "]"}]}], "]"}]}], "]"}], ";"}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"Compute", " ", "the", " ", "remaining", " ", "values"}], ",", " ", RowBox[{"using", " ", "the", " ", "AB4", " ", "formula"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{"i", "=", "4"}], ";", "\[IndentingNewLine]", RowBox[{"While", "[", RowBox[{ RowBox[{ RowBox[{"y", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "\[LessEqual]", "1"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"AppendTo", "[", RowBox[{"t", ",", RowBox[{"t0", "+", RowBox[{"i", "*", "h"}]}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"AppendTo", "[", RowBox[{"y", ",", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}], "+", RowBox[{ FractionBox[ RowBox[{"h", " "}], "24"], RowBox[{"(", RowBox[{ RowBox[{"55", " ", RowBox[{"fi", "[", RowBox[{"[", "i", "]"}], "]"}]}], "-", " ", RowBox[{"59", " ", RowBox[{"fi", "[", RowBox[{"[", RowBox[{"i", "-", "1"}], "]"}], "]"}]}], "+", " ", RowBox[{"37", " ", RowBox[{"fi", "[", RowBox[{"[", RowBox[{"i", "-", "2"}], "]"}], "]"}]}], "-", " ", RowBox[{"9", RowBox[{"fi", "[", RowBox[{"[", RowBox[{"i", "-", "3"}], "]"}], "]"}]}]}], ")"}]}]}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"AppendTo", "[", RowBox[{"fi", ",", RowBox[{"f", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", "+", "1"}], "]"}], "]"}], ",", RowBox[{"y", "[", RowBox[{"[", RowBox[{"i", "+", "1"}], "]"}], "]"}]}], "]"}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"i", "++"}], ";"}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", RowBox[{"Transpose", "[", RowBox[{"{", RowBox[{"t", ",", "y"}], "}"}], "]"}]}], "\[IndentingNewLine]", ")"}]}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"b", "=", RowBox[{"1", "/", "0.5424"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"c", "=", RowBox[{"-", "2.9"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"f", "[", RowBox[{"s_", ",", "u_"}], "]"}], ":=", RowBox[{"{", "\[IndentingNewLine]", RowBox[{ RowBox[{"Cos", "[", RowBox[{"u", "[", RowBox[{"[", "3", "]"}], "]"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sin", "[", RowBox[{"u", "[", RowBox[{"[", "3", "]"}], "]"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "1", "]"}], "]"}], "\[NotEqual]", "0"}], ",", RowBox[{ RowBox[{"2", "b"}], "+", RowBox[{"c", "*", RowBox[{"u", "[", RowBox[{"[", "2", "]"}], "]"}]}], "-", RowBox[{ RowBox[{"Sin", "[", RowBox[{"u", "[", RowBox[{"[", "3", "]"}], "]"}], "]"}], "/", RowBox[{"u", "[", RowBox[{"[", "1", "]"}], "]"}]}]}], ",", "b"}], "]"}]}], "\[IndentingNewLine]", "}"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"res", "=", RowBox[{"AB4", "[", RowBox[{"f", ",", "0.01", ",", "0", ",", RowBox[{"{", RowBox[{"0", ",", "0", ",", "0"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"res", "[", RowBox[{"[", RowBox[{"i", ",", "2", ",", "1"}], "]"}], "]"}], ",", RowBox[{"res", "[", RowBox[{"[", RowBox[{"i", ",", "2", ",", "2"}], "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", RowBox[{"Length", "[", "res", "]"}]}], "}"}]}], "]"}], ",", " ", RowBox[{"AspectRatio", "\[Rule]", "Automatic"}]}], "]"}], "\n"}]}]], "Input", CellChangeTimes->{{3.7840256045524616`*^9, 3.784026021380303*^9}, { 3.7840260731422634`*^9, 3.784026202603668*^9}, {3.7840262430129795`*^9, 3.7840262481792755`*^9}, {3.7840262913157425`*^9, 3.784026379054761*^9}},ExpressionUUID->"faac927f-4c4f-4fa7-9b67-\ 07cf15a280a5"], Cell[BoxData[ GraphicsBox[{{}, {{}, {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.009166666666666668], AbsoluteThickness[1.6], PointBox[CompressedData[" 1:eJxdUmk0FAoDtUS2LMnYXrx4KVFJi8rrXbQnisqjkBKZj7RI2SstKo9QIsS8 h8JTlsaWJckS2bdsYcw+Y4xh0KT4+v5+95x77rk/7rn3x1119qK9h4SYmFjF T/5P/x8FkSyL+KA7kHNKmGAaakDMqm9bQeATjCfpnNQ23I5Rqd7FI+RUkCdU hzKNrSHzh3TvqmsZOCz2+VreakdYLv/w9XhMNhriV7++MHgWGt5SH45mvcaS ZuNwob43MmUOpXUOF+CXBE5Fmv8VqG/cki44T8b2w74Gs+oBKElWOm7LL8aY nPTMpoAwyGbl6xVeLcO1qzbVxU3hEDVrazHnyuHjMVN2vOQepi9G8UsCqlC0 7HuCp3gk5kULMgHB1TDyPJ95LOsRIsxeDa0vew+n3ou3y0/EwbQ67Zv7TA3s BtrElKXiQWFnuIib1OI4IUevkZKAMWm/IJ5XHVTIiStvc5MQb2K3T5lUj/Ez Z3K7V6XCtCl4l31PAwgXyM/X1JDga7ibmCLTCELW3EM8+AfpgRSPvh1NMDsb 3jTjkYF3XCmvSa9PmCuVzpJ3fAHCFlWrtvhmiK6eWmNxOgsdFUpEt+oWPF69 UudocA6Ia/eJHrFacTeIZjP7by6W3DnksnppOyrWit3gTr6GnP+25e8S2/H2 RLvNdEY+uKpVJgqGHbBSK/SwHSwAX+GO4ZLSDgRfSmxw034DixZRY+reTujg ha+kBxkqYWW7ujt++qdtqYqlRSjhFCs9d+6CV6zMtNOKEqiQ9htN0bugvRjg XRFYioo8iS/NF7ox9yY3WIFVhpiUgC4DYTcM75Rq/Hq6HElKtbrfAnrQteJw Q/NwBbSWCbdafu/BVlFWNtOzCmIBR3mzob0QdlyP2jv3DnUNgTcJC724uO1t yJeGagSzzkUlh3xGjbr/SsXE99h5ZItauOgz1m4OLdPxqUGfCz+4/moflpl2 2Yl2f8AsNO97TPThWb/1cIRuLSQD15519uxH5fX2Y+SFWuR5jj7OGuoH85fA HC9KHUhS0UsP2A1gzf4wanRDPdi2hkt31A6goH1RUaWgAYxIZ/OgrYPo+jXN rv/5R7RI+uXIZQ7i8ScOuyOqERb/MSlhLx+ChOjZVvrNJkwHxOmq3xhCnubT 4h/XPiFTN2w+jj0E52cd6tKXmpFc5FTuYv8F+j6cTWzvFpx+Gm8QUPYFOiEJ SnH/aQU/5reTVJ1h1BmMCrjebSj+i15PCh/GhRjFdQKzdtTVLG4spA3DSJdy 3aS8Hb1U02zNvSMgun0dyUYHakJ+ie9PH8ERkVf4sboOpAvK7ILERnF+3jld /3AnTGdTmds2jOKO94Nc5e5O6HeGbPjt1Ci85yM2qLh0oSbJfXx3xCheaaiN aTO7EKeyVy6xcBTfb5pmr7nSjexZycd6X0ZBeLYhxmihG+mE8EK6NAUP22Ja dSJ7cPbJKzJ1IwUfT3Mbvmr0IpIVtUTXkQKb7jiFgpe9uCJnZ5R0g4Jj+eYn YfYZTzY6tf35ggKdwrCU5IbPOH7f+s3JZgok1uzhVTr24ZvTr0WZAgrcOwjx qZw+nGvb0mJOGMNZpqv8ttB+mHa5Lvt15xheFN+dDVEeQCQkSSdcxlAWJmHs kzEAzewNDSM3xtCnbpf/ffsgLJYUDL7/ewydi53Fa1oHoePy++HFmjE0E5P/ ZrgPIaX7XkoydQyukwbim78NIUjx0rkkSSqSIhcfS8d+gcm2dUcW9ajY7Wuz 79TaYdx666nVaklFtu+JTUbVwyjuvxcv50aFP4medMlxBP3L20mloVSsrcsX rhaMgPwww3goiYq0lGae3t5RaPi+2x5YQkW/+LVIeuwoyCf8ux53UZHrUG/d NzwKas5XDRM+FUSrB+VfjShYJvOPqb0cDU9F3sF7Ail4//B3h4XfaNC3CkNZ PQV7ZB3qtoOGaL1XrcdXjMFa83GzlBMNYX+ZSaidHUNJh/QH4hUaiuxX1s3l jcFvfbQsMZKG8rL7qvM/xnBwaoOcTAYNEwLTb78cpsJsPtDMuoIGk8YP8ad+ 7p6vmZQ366ahkmK8vIhFRYOMr2w7lwaLcrGcNWY0vPaMLdWSpKNioC6z4C4N mnnrs7S16PBy2XfK7mcuh1bm89mEjt1WQ8ZS+nTsSa1Ns99Px8eit/c/XabD 5yVxJNaFDlrOFtY/1XSszm1seeZHxzoF8bYHSgz8WXVb+eIDOpa2hzSFuTKw P9zkKCHtZ99O7qWbrxjQahs2e0Km41Mi2S3qOwP2MbFXeY10pFekI8OaCbOi KrLBCB0eD/XefkhiIny+9eUeIR3tQ/7hHDYTcE/9YS3LwH1F0RH1HSy80vrh YanDwHY/p+ED91lQKTSJMdzMwEGFzvHQzyx88PJfJ3uAgV1153cVGbDh7HmD QnNm4Pttgyc8fzYayPOn319moNLpdPHqOjYIncq26fcY8F3FIjqv4GDvm+Gj j5IZaOhzvhrjzgFJdUH2YT4DPuLuSdWFHNjyFwiJdQyQ/uSncsW5CJZjKFcN MGBkHG2lbMdFATE7UmySAfm811brSVx4p/+rdk6KCb8AgaMlnwu72a1mPC0m XJ4H7bT+Yxy2qoLLJBMmSCn7Ew5GjSM82tkifB8T05ceaJoPjUPE26OW4MxE sqqkh44RD+vstp+iXmHi3VsfwlQgD875Ij/vB0xsXTlUQ/7Ig33JVMlmEhO3 lpjOnVOfQJ9Vhr9FCRNWV6PMxD0noB6lJohvZSItt3YhgjwBnoaV+3YGE77W tMo5CT7KL90TGS0wkd9CHrS148OrQ77nKoEFUv4+cmQaHzeHTeXUN7Kw4Hs3 MZfHh/vNH2mEAyyYN7EYn1dOIlwl4NH1MyzwFTvqG/dMgiip92FXMAvs3KcS /3pPQm2tA+FCPAutFX/lBMVNIqKt8LhsPgs2N5I+mpVNomf81g7CJxa2PNV2 o4xMQiFw1cUEBgvTdrUXrkkLIGHjEx0jwQZR/yFzxliAllf6OlK6bGQipsL1 mAC0Rqe4r+Zs9N3U6MkLFGDA84u1rxMbvlEBy1lpAniJv0wgXmfjrbXqEel6 AXr/2f+aF8/GphcOrrLjAozkOb35QWaD9HelwqTKFMZnpskp3Wzsr+3VLTOb Asdc998WIRs7m2wd3F2moJ1ie/2xGgdNpMXzvPApTHAHqJPbODhKkpR2yJpC yLVb6SxHDjREzQPJLVPYZ8M4GhrMQc9ey4eVU1N4c1nLtyCVg/DN/QVl6tNw 1TA4F1HDwVnjWlb079Mg7k0IF2dyECpQqjY/Mw39xRWLqxS4aJPvrqm6Ow1p +zkhZxMXbt/drhFyphHtSjvn6shFoLLfScvWaXS8uLIl7AYXDpVe7F1T0+hq bV9/+CUXMh4hiXIEIZ4sXllsbOPi6TnxvswdQuRr1puLvnKxcDBDVdFFiNUU luOg3jg8lV4XWt4Uwu+Ad/5lm3EkmB3SQLoQtB+qNyoDxtH9oi19Sb0QDidz ae8yxrFjuUJBDEuIky2d1MCOnz8ejk2kyM1g2UjXmfGFcew7JJM/YzyD2KEH BhrreZAvDD30yXYGyRa7h36c4sGtt7zU7dIMtjLylFIjebDiPj9Ijp3BXeWe jaJyHhbtvA7VFs5g7RGZaCUeDxUdq5Y+6ppBym6BC0VnAhEGE03LhTMwORhX 5Gs3gadVxFGbFbOQXPgzquTOBGJDpeP+2DKLL1kbq8tLJ6DpSlxHPTYLMf8M +WDeBPSbPgp3+M0iMkhSeUaPD8vb1rss42bxMuLeUUMnPojC83Zz+bO4nHJ3 s3oMHz4Nu2+ebJuFt2yeRFUDH9xEVW0ibxYGSun+imKT2PkiyH+V/Bz+CxS6 99M= "]]}, {}}, {}, {}, {}, {}}, Axes->{True, True}, AxesLabel->{None, None}, AxesOrigin->{0, 0}, DisplayFunction->Identity, Frame->{{False, False}, {False, False}}, FrameLabel->{{None, None}, {None, None}}, FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}}, GridLines->{None, None}, GridLinesStyle->Directive[ GrayLevel[0.5, 0.4]], ImagePadding->All, Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({ (Identity[#]& )[ Part[#, 1]], (Identity[#]& )[ Part[#, 2]]}& ), "CopiedValueFunction" -> ({ (Identity[#]& )[ Part[#, 1]], (Identity[#]& )[ Part[#, 2]]}& )}}, PlotRange->{{0, 1.0000137896722179`}, {0, 1.4402174557711889`}}, PlotRangeClipping->True, PlotRangePadding->{{ Scaled[0.02], Scaled[0.02]}, { Scaled[0.02], Scaled[0.05]}}, Ticks->{Automatic, Automatic}]], "Output", CellChangeTimes->{ 3.784026281221165*^9, {3.7840263504211235`*^9, 3.7840263803948374`*^9}, 3.7840264985415955`*^9, 3.7840382521258388`*^9, 3.784040540574254*^9, 3.7840409618233747`*^9, 3.784041004159931*^9, 3.7840410921691227`*^9, 3.7840411935196123`*^9},ExpressionUUID->"b7957176-ff6d-41cf-a911-\ 23b6125aeca0"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[TextData[{ "\:0420\:0435\:0448\:0435\:0442\:0435 \:0437\:0430\:0434\:0430\:0447\:0430\ \:0442\:0430\n", Cell[BoxData[ FormBox[ RowBox[{ RowBox[{ RowBox[{"u", "'"}], RowBox[{"(", "x", ")"}]}], "=", SuperscriptBox["10", "6"]}], TraditionalForm]],ExpressionUUID-> "580175f4-fb8f-44ca-a579-cfd83ee703b1"], "(-u+sin t ), 0{{3.747918045893719*^9, 3.7479180792016163`*^9}, { 3.748019567260964*^9, 3.7480195984310884`*^9}, {3.7491799306613116`*^9, 3.749179945330095*^9}, {3.7491820390416107`*^9, 3.749182140816432*^9}, { 3.7491822953940344`*^9, 3.7491823155841894`*^9}, {3.749182378833807*^9, 3.749182424987447*^9}, {3.749182543546228*^9, 3.74918254777547*^9}, { 3.7491855618392*^9, 3.749185791014152*^9}, 3.7494816355928173`*^9, { 3.7494826582226067`*^9, 3.749482678435937*^9}, {3.749483068682683*^9, 3.7494830721446886`*^9}, {3.7494884082644987`*^9, 3.7494884424507885`*^9}, {3.749488494590744*^9, 3.7494885492422333`*^9}, { 3.749488662765763*^9, 3.7494886848150673`*^9}, {3.7494887168485804`*^9, 3.749488740273077*^9}, 3.749490389264352*^9, {3.749789440131727*^9, 3.749789499828309*^9}, {3.74979010532883*^9, 3.7497901250714474`*^9}, 3.7499847604043837`*^9, {3.7499852544136286`*^9, 3.749985272769679*^9}, { 3.7515612145185156`*^9, 3.751561319842022*^9}, {3.751561379568722*^9, 3.751561400725582*^9}, {3.784022901846876*^9, 3.784022907924224*^9}}, FontFamily-> "Times New Roman",ExpressionUUID->"4ed04c9b-35d8-4b0c-ad3b-76a639522d18"], Cell[BoxData[{ RowBox[{ RowBox[{"Clear", "[", "x", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"f", "[", RowBox[{"x_", ",", "u_"}], "]"}], ":=", RowBox[{ RowBox[{"10", "^", "6"}], RowBox[{"(", RowBox[{ RowBox[{"-", "u"}], "+", RowBox[{"Sin", "[", "x", "]"}]}], ")"}]}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"exact3", "[", RowBox[{"x_", ",", "l_"}], "]"}], ":=", RowBox[{ RowBox[{"l", "*", RowBox[{ RowBox[{"Exp", "[", RowBox[{ RowBox[{"-", "l"}], "*", "x"}], "]"}], "/", RowBox[{"(", RowBox[{"1", "+", RowBox[{"l", "^", "2"}]}], ")"}]}]}], "+", RowBox[{ RowBox[{ RowBox[{"l", "^", "2"}], "/", RowBox[{"(", RowBox[{"1", "+", RowBox[{"l", "^", "2"}]}], ")"}]}], RowBox[{"Sin", "[", "x", "]"}]}], "-", RowBox[{ RowBox[{"l", "/", RowBox[{"(", RowBox[{"1", "+", RowBox[{"l", "^", "2"}]}], ")"}]}], RowBox[{"Cos", "[", "x", "]"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"plots", "[", "appr_", "]"}], ":=", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Max", "[", RowBox[{"Abs", "[", "appr", "]"}], "]"}], "<", RowBox[{"10", "^", "6"}]}], ",", "\[IndentingNewLine]", "\t\t\t\t", RowBox[{"Show", "[", "\[IndentingNewLine]", "\t\t\t\t\t", RowBox[{ RowBox[{"ListPlot", "[", RowBox[{"appr", ",", RowBox[{"PlotRange", "\[Rule]", "Automatic"}]}], "]"}], ",", "\[IndentingNewLine]", "\t\t\t\t\t", RowBox[{"Plot", "[", RowBox[{ RowBox[{"exact3", "[", RowBox[{"x", ",", RowBox[{"10.", "^", "6"}]}], "]"}], ",", RowBox[{"{", RowBox[{"x", ",", "0", ",", "10"}], "}"}], ",", " ", RowBox[{"PlotStyle", "\[Rule]", "Red"}]}], "]"}]}], "\[IndentingNewLine]", "\t\t\t\t", "]"}], ",", "\[IndentingNewLine]", "\t\t\t\t", "\"\\""}], "]"}]}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.751556858181094*^9, 3.7515569723424397`*^9}, { 3.7515570580181684`*^9, 3.7515571041856165`*^9}, {3.751557152390895*^9, 3.7515571560038157`*^9}, {3.7515588752945166`*^9, 3.751559009715211*^9}, { 3.7515592708620443`*^9, 3.7515592965272665`*^9}, {3.7515599556701965`*^9, 3.751559970390748*^9}, {3.7521742748175573`*^9, 3.7521744452983103`*^9}, { 3.752174477639867*^9, 3.7521745458329954`*^9}, {3.7521746109708395`*^9, 3.752174769825224*^9}, 3.752174814186925*^9, {3.7521748995772405`*^9, 3.752174908051389*^9}, {3.752174995267666*^9, 3.7521750136281757`*^9}, { 3.752175099765067*^9, 3.752175158597514*^9}, 3.752175259941711*^9, { 3.7521753368668237`*^9, 3.7521753827886343`*^9}, {3.7521754325722404`*^9, 3.752175603081993*^9}, {3.7521756885552454`*^9, 3.7521757608530293`*^9}, { 3.7521758256240716`*^9, 3.752175837822115*^9}, {3.752175877769332*^9, 3.752175907667282*^9}, {3.7521764971382256`*^9, 3.7521765346098366`*^9}, { 3.7521765752366657`*^9, 3.752176575890293*^9}, {3.784040929132264*^9, 3.7840409453159204`*^9}, {3.7840409895223813`*^9, 3.784040997660703*^9}},ExpressionUUID->"07d23880-93e2-4b4d-a3fc-\ b8906f99386d"], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"h", "=", "0.1"}], ";", " ", RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{ "\:0417\:0430", " ", "\:043d\:0435\:044f\:0432\:043d\:0438\:044f", " ", "\:043c\:0435\:0442\:043e\:0434", " ", "\:043d\:0430", " ", "\:041e\:0439\:043b\:0435\:0440", " ", "\:043f\:043e\:043b\:0443\:0447\:0430\:0432\:0430\:043c\:0435", " ", "\:0434\:043e\:0431\:0440\:0438", " ", "\:0440\:0435\:0437\:0443\:043b\:0442\:0430\:0442\:0438", " ", "\:0434\:043e\:0440\:0438", " ", "\:0437\:0430", " ", "\:0433\:043e\:043b\:044f\:043c\:0430", " ", "\:0441\:0442\:044a\:043f\:043a\:0430", " ", "h"}], "=", RowBox[{ RowBox[{ "0.1", ".", " ", "\:041f\:0440\:0438\:0447\:0438\:043d\:0430\:0442\:0430"}], " ", "\:0437\:0430", " ", "\:0442\:043e\:0432\:0430"}]}], ",", " ", RowBox[{ "\:0447\:0435", " ", "\:043c\:0435\:0442\:043e\:0434\:044a\:0442", " ", "\:0435", " ", "\:0430\:0431\:0441\:043e\:043b\:044e\:0442\:043d\:043e", " ", "\:0443\:0441\:0442\:043e\:0439\:0447\:0438\:0432"}], ",", " ", RowBox[{ RowBox[{ "\:0442", ".", "\:0435", ".", " ", "\:043c\:043e\:0436\:0435\:043c"}], " ", "\:0434\:0430", " ", "\:0440\:0430\:0437\:0447\:0438\:0442\:0430\:043c\:0435"}], ",", " ", RowBox[{ "\:0447\:0435", " ", "\:0434\:043e\:0440\:0438", " ", "\:043f\:0440\:0438", " ", "\:0433\:043e\:043b\:0435\:043c\:0438", " ", "\:0441\:0442\:044a\:043f\:043a\:0438", " ", "\:043d\:044f\:043c\:0430", " ", "\:0434\:0430", " ", "\:0438\:043c\:0430", " ", "\:043b\:043e\:0448\:043e", " ", "\:043f\:043e\:0432\:0435\:0434\:0435\:043d\:0438\:0435", " ", "\:043d\:0430", " ", RowBox[{"\:043c\:0435\:0442\:043e\:0434\:0430", "."}]}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"yImplicit", " ", "=", " ", RowBox[{"implicitEuler", "[", RowBox[{"f", ",", "h", ",", "0.", ",", "10.", ",", "0."}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"plots", "[", "yImplicit", "]"}]}], "Input", CellChangeTimes->{{3.7521763186839857`*^9, 3.7521763295148087`*^9}, { 3.752176361808408*^9, 3.75217642428076*^9}, {3.7521764942798386`*^9, 3.7521765018974943`*^9}, {3.7521773289252157`*^9, 3.7521773314847393`*^9}, { 3.7840407296097493`*^9, 3.784040852884266*^9}},ExpressionUUID->"d6e139a0-90b9-4a15-a972-\ eed403abef6d"], Cell[BoxData[ GraphicsBox[{{{}, {{}, {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.009166666666666668], AbsoluteThickness[1.6], PointBox[CompressedData[" 1:eJxd1HswXHcUB3DvEKp2gxVZy1oSSWsIW0TxO0nqMWmqwoTWH5qReEYaHRpN Os1gkrYqMhIkCN1STFqNV+IRwYqKkKLxJmEt67XsKyxBhFb37D89M3d27tzf +Zw7c79nmaHn/cPUVFRUwv69tn//X7k52/WQLDW16f58/SFR3HeQ1KzHwyvn Ogjbcbv6yAH1rqsdwl58PkKe0hZXnUTDRKHwyUuNkXtTZTw8LyCGYpZvKmWS RIRv1wxJLdNvismfwn4hoQLNNUM+Q7o6t0tEHsjtnNQoQvRkpPAkPzBAPI/n XxFfhs3h1isi9JcI5Yqg7UifGPvl5Pg3aX4JAxKct0KK3Bqjun6UordKNqOT gi8tSNFbI+kf1zYkbUjRe0OyXIpV6x5J0dsgA5VnJlr2S9HbJGyrrSqajwQ9 FRAUXax8baKcrwLZnsHO0ZwF9FWBXhyxxR4XYr8aOBfEXe4dncV56rC3YPzL 8Pxp9DSg7It9pRV7BehpgCyHenHzMh89TWjXDBpkB79ETwsmPoxWvxA7gN4O qGw8zGvmdKOnDcWfa/0GBi3oaUMRfXnO/lgxejpwbehFuZHgHlfh7YRWVn38 un0rV+HpwpqPeKU1+DlX4elBeDM10TlukKvw9KAmdCHa3X+Uq/DegZajkbQo OR89ffAavul9MlOA3rvQbJyd9b7aDHoGYP7BobN1znN43gBu5Q3b+9vNo28A zYO/P34lWOD+1+5IgfcO5B9zPSHGeRQIXdeMuxonQY8KJYKMKpG3FD0q0GNz G/w6pehRYc29RRyyJcX32QWdS4fKMvlS9HbB6chmLsRK0TOEX9MPxieXS9Az BNtvc7q9fhGjZwh/HLS+ddRThJ4R9DBSI2QF8+gZwT7rksHv6ufQM4YLrYs1 zzJm0DOG1KGgkiDHKfSMoT/JZN2FM4EeDSp4vNP6DWPo0eC+KtshTzyEngnY RC0nLk/0oGcC1J7KLeaJp+iZAHN85lKUZzV6u8GGl1Ib7nwDv/9uqObfiTMd bMS8mMJHjI03Gil/Yb5M4auF21besj7Mjym4DfnV6wpHMG97IDEpfjSFM47e HgipYgTdTJ5Ejw4e6o/oO7Sm0aNDTKgspDpgFj06XD/+07O6MCF6ZvBkObD0 iZNyf8yA4zaW4NQiQo8BW+ydD2gayv1mgO335/Uc1iToMWCt3SWdylHusznk WXbJDeTKfTcHFTdaX8Br5f+BBXhk3unrvSslinxZQFJ795+BOsp9t4Ab7l9n JphKsN8COuS2Xtk8Ec6zAEpGVjf/swX0mDDswNF0ShOix4QQ145yu8RZ9Jgg 9rjbpOc8jR4TuIX908W1k+gxgRGaVhKizUfPEpK1CmT7T71AzxLOJZ6NNY/p R88SHFen7Cp+6ETPEkZuZxz5dJSLniWYRm5cexuSgx4Lcu37n+u438c8sCCj 5pMaX782zA8LJAlJ2clnejBvLFg0K509FT+EeWLB27Kov8FnDPNnBcW8Quv8 yAnuP33D2gg= "]]}, {}}, {}, {}, {}, {}}, {{{}, {}, TagBox[ {RGBColor[1, 0, 0], AbsoluteThickness[1.6], Opacity[1.], LineBox[CompressedData[" 1:eJwVlnk8VO8Xxy2RrTAjZZsZEVJRWr5JOScqFUqkpFIpUnZRlFZRSJbIklQi EpJoES6SJZVClizJbpi52Q3id3//zH29X3fu8zzncz6fc6+yrau5nQAfH996 6uf/14dXbIXXnArW51PNamX379cvrWWZHGSdAbZiS6Me4wzMlRX6sVg+oC06 UyzK9IGHRSFMIVYglD3LjLdkBoLTFi/nCWYs/DvTrSxJXWunHSL/Mp/Dpixi zyal5zC/uUOdw3wPCrrddkGL3sOTz4/cB5iVYGdlZZg2XgGZLwObnbR+wW1W Mins1QSn4UfsoEcnyGz0d5vb3gF8ieN5zcGdUM//kO1u3QEx8xRbPid1Qkes PDa5dEBFpT3jeX0nOP0uMg6I7gCN/TOJdpu6oPVu8B2rvg7oc1BPb+Pvhkp3 a/7aW51wNuJy4fewHghPd75fVdAFLt3LO3My+yFjfN53A8FeaLdpEimu7Afh y/2/hWR7waLxttbXrn54sSLS/Z1GL2ys6vXplmdDvm1T14BpLwi+SpaSvcWG 4XaRmKDoXojxVda/cGwAnnp6Eg/U+6CYLherK8WBl02vhHM29cO6kIrC7Ss4 sNht9uAe4354JuzdtW8HBx4LbIyuP9wPd3j12md9OXA2E7OyfPvhwO/Isrg+ DriWqgoFFfQDO01qdLqICyc/z19erscGMeWb6r83/oWu/te3lTQHoGJzuPaV 7X8BJEzeOv03AAFWCf8pmf8Fxim38bRtAyAQ9s7osONfUJl19x+1GYCp2cHT DQ//wiu9/wQ+RAzAQItlarXAEFxfmuK9cGIAvsSoLy/6MgRxjK742OxBCMpZ t+ZY0xBY2fDnfigYhJ3ft+rOdg9B5pMv7KqKQSidf2TX5rkh+DKPvPembRDy z4edebdmGJSETjByRDmQsZ/3POv+MPi0vN65+wgH7kpVaT45PgLffYrZb0Y4 8P4/pdw9LiPgffWvuOk/DnTZuMLMpRGQUQw1+CrMBb0M+v6D0SOw1dSuLliO C327jlxe8G0ErGT4vvDpc8HQj1PtozcKbm+WmZTe4IJLGlir7xoF470Hnk0G cSHuR3hX3YFRqNV0yZS9x4W/zA08bY9ReHH1eZbkUy4k5F9Z2pM6Cto1m5pH i7kwOSbpZS47Bg2PJzgB01xQUbKdm1MZg7Wu0TcsBEjYsy0nMGPNGLSLdjuL ipKQHGGVIGI6BrX2M+vWyJJgrv2kvNBvDJJE2udctUm4cmBkn1P4GEyl65cG rCfh+eXtLXKPxsDRqv19gB4JfF/6/3rmjcHTwHPqBkYkZDjoyK8YGgPOVe2o Z0dJaAy9mdQ4OwYR5C0n4ZMkzHtbrxUgMQ4Rpp0Wex1IsBa6ZPhHfRyWf0XD eA8SRBI/OkXbjMPkvcjNSwNIWFu5aGKb0zh8Y99Z8TmIBJu/p68P+4xDUb33 gsOhJOTqS9w3jRoHjzOd8f/FkGDbvJ8Q/DIOp0s+hDWnkiD24I+ze9M4vLxy IvJdOgnZ1i6Kv3uo9R9a77yaRYLgr1s+efwTkC74+k3lWxKSGvPWuv83AQJb T9ibfyLBJMaoo23bBISOxCT7VpAwerAuzMR8AmzWR10PqyJhewOHo+48AYVW MN/rBwk9P1mpbYkTYOAbKLWmhYS7URkHTLIm4ED6s5mUNhI2WG4SyiuYAL9K tpfIHxJu1VnYRjVOwNb120UDuklYHdkuJdgzAWTpGtbTXko/C2fCbWQCVpcN FTzvJ0GjNkDRZOEk6L3fcfsch4TvEfSq9wqTEBy6OnsDSYK3+WMf9eWTcCwx +mjnXxJYtFUaURsmoXDWNNBnmITKH+/rBbZRHLR+HW+EBI/wHf5u+ybBVBiP 2o6RIL+vdm2bzST0i7sK5o6TcPb7YNh7n0k47BD2UZZHAi3MB9RvTcIK+Nm8 bIqEvL3C3MjISVj1pNOZMU3pL3kvXiBxEs5ckPDhmyFBvJpp7PZyEoQCs2aq KM65m85rzZ+E3Es4ee0fCUf26KYaf56EZXWbHZmzlD8Wlh143zAJZuPrLVMo Tv9qLqTePQnVy77kLpkjYX/I79eRw5OwJODlnXMUz5g42Qrw8aBA/nHTW4qT JHhSbgt4EF3uHN9NsekXf6JVngc9ZpN1MxSPB9NcjDV44HBY2+8fxQnGjxTf r+eBnRMtpZfiHeIrq9QMeZDdn4UfKOZ+fucTacaDJDlNM2+K7wdt1xCw4cHY 0oyfDIr1d9fUuzryIF3C50cmdd4e0WP+rd484BS0GGpQHFo5sNY4gAcf/m1Y FUzVuzHQu+PdPR4Ue4yHNFB6tO8UCld7woOZxFj7BRTfFomAyEwe2I5kvl1J 6bm6gsHlz+fB1ZS66+spvRtvvYh3reSBxJ93n9Spflwz2mjcWk+d/5HcTcFJ yh/zP/F2d/GATpQWVFD9+162L/XdEA/W8aLPXaD66x3QdkBtjgdtHdeeSY5S /tjhKBQpMQU+g74Hwyl/uH+6aeuqPgU37Eb+20f5idG1Q3h87RR4VpkcyqT8 ViUgmuaLU9CyXLF/YoCEZRgyFHhoCp6nHLYy6yOhxmZvlJT9FJiuFN10rIeE q5eldaM9puCaR8pd6y4SmvKiriUFT8HYn5wo2XYSApoOqq6ImYIQk+LdDa1U 3iflKl4lTUE5b/LCzWYSQtYnSBIFU2BglopJ9SRszXqW0EROgevOxlwelT/u NweD4zNT8GBf/gFGJQkPOJo9PSLTID5yxVe7jIQxzZerRpWnQTCz5ejiImp+ Jb/NX7h/Gm781c/4l03CgVKf41HHp4EpLpoW8JKaB5168xSdp6F35Qm9mRck HGUVGS8PmIYXkV6qGUkkSD+o+GX4bhoKd/7eqxZFQuH7oCufS6fhdNT51TLh JDg2mizd92MaNohtf8e9Q8KnRT/O2LCnYf7NgvIzN0nwCWua9Facge+vtFZ7 U/Pujz9bNvPaDFQtd0mY3E3Ni6T0vHUhM/BJx2RL9nYS9D662HyInYEljZGv rJCEKL6R5IrsGdDcEKruRM3jXZem13V2zcDmzqaUDAYJr9wkLJbs+gdFvR8W JXC48CWO7pNx4B9UfYjSjOvhQk+p/CODU/8gmZc4dfU3F+Tllg84Xf0Hu1ra 20V+cOFm8Xa/4tx/EOS64pDUay4coF19fXbpLKyrdyqf8uCC2+aApjntWbBM ebYk6iwXguxD5iK3zMKW4XNX5Gy5QOQ9MCasZmHgkcnjdjMuaJx810kPnQVW hsCgwyouTL0eohdMz4LyAqe9pzs4INM2qWsuOgdsqfr+5EYOaInwHe+VnQPD sFn+6m8cOHlkYbq0zhy89psob8/jUO/pFdvsHebgdHiHp2AEBxIsT3lK/pwD uvpAg/FmDjy55MhefZAPIyQcF7KuD8LqItHJW0f5UMtGKDbMaxCK5qUK/T7J h+d53wzJM4PQHtLFCnHjw3b5tYsvmA8C4/HRg/2BfDgndevmq6WDEFe6t/RJ Ph+6xqguaSQGIGLBugTaUn7kGWUdVuGw4cajmX2jg/zYJjHPKXRRP5itfue+ YoQfh0MUZJYI9wOj+Fy4LY8fRdyuPAkf74O8Dvb370IC6Ln3RfjBhj4YXta0 N4MhgC2W+F07pg9s03NN7c0E0IHZoCG/uA+2vnfZ1fhaABmiNSsVxHqBr7Yd 8y8K4uX0iY/S5V2g5mnvvV9TCGWnaAfAsw26xqIPiBcIY59SbMXTthrY3X32 coS7CK5o9FRbI/wG9nSFrWAwxDAiWcdr0JUg2B9tr1VkiOPE7LP2auE6Qku6 Xsnw2AK0FuHz/LSynegJ6bkRZbcAFx9jyEcYtBMPxSZ6ex0X4EDsr5j9Vu2E xLwl2Xe8F2DWp8b+PL92gj1+aEd9+AJsE5x/TPxXO5Hc0upypnQBeracqHt5 8w+hmNpNhC5fiDlmRlvSfnQQYjB2omVkIWZJFHPuG3UT3ls1RN5MLUQ7695/ Qke6iR7Dw5mh/JJo8er3KTu3bqJkZ/GUgaQk1q6KmB6N6SYumofce64piaRM TLplfzcxYLfs0/kTkqgidmSn1K0e4ssdSw3paklsTT62TvhNL7Ep9PY39k9J nO1bQD9V0Uukhn/wLG2RxM98r69n/eolbt5XLr7AlsQmtmGTzFwvwfg+vpIU kMIDUXz/RI36iPqmkk+WDCnMjb3yMK62jzDiWE8stZRC9XHxVTNt/USCjF/t VSsp7GXVb6jg9hOjei9ethyWQr+cjMHrs/3E46Dp0/dtpTDhrNLqQiU2MaUe 3yjqJoUrtQQ6rA6ziXTb1vd/g6TwmcKc74JaNiHVeMy3sEgKRaTjEzVzBgj7 uVtWCqVSqCljlhxdPEDkq2Wt8y6XwrZUNbGhbwPEGS8+zppvUijmfafOoX+A KKE9OZrcLIXtd31POyoNEp6mHfp3xqXwdehASM71QaKx5BS/9Upp9MqrFZTX 4xB/t/6Sv6wtjcc69p802s4h5hfvXfdYRxo9d30ptt3LITYQm073bJTGKakq 25MnOUTkB6mvHtul8clkdmBwEIfYm5MfE2wjjXfkZNd7/OQQ9mt1sjNPSKP/ 8MFs3zYOcSU7perHKWnkOlp89+zlEOlZEbOLHaWxZ7d74yoehxDLcDiVdEEa f1nEXi9R4BJlyTKrC8KkMTXE6ZH7US7Rphq0q/2eNN7zkXthYcclxp7O2QpG S+NMSeg8FWcuoZrIjtr1UBrP1qaHBPpyiRsJRdM/n0uj+uiMp1ccl9CPdirn lkhjXMiNx6U1XOKAbEe7dJk0rr8ZccKniUs4Rx2cWlcpjVbDUxoK7Vwi/p7B qkvVVD2FoXnLOFxiKnTJvfkt0hh9293ngDBJ5N4utWGNSeNzyWpplfUkkbq7 asmtSWkUenTvuKIeScRJ1NRwpqUxIDlYfg5J4lrY7x0fBGh4KY9FXDMhCZOY Ka0DUjQsUSw46GlLEmDN359Pp+Fmo8A7706ThI6iyFOVxTRc8MxIoseJJBY/ XrR4SImG062vWmfOk0RHyuq54BU0rGv2PCkXRBI/z/z3fliLhrruMUEpd0mi YoX+uUM6NMzLuNXOuEcSGS+Ne9V0aSjfwxWvfEAS3m/tq4t30PB8vr+YeDpJ OPo4B2rspqGe/8Z9HS9JwkbP0zDUlNp/XVZm4muSMCSuvz2yn4b7TZc9bcsj iYXl8Y8mTtAwa/1OTcMKkuAPfHrIxo6GWyW+t9lWkcTo7jT6Jwdq/xWb7jt+ I4mmb29vRbjS8PD1uXkb60jiS1jhVp4HDXkOy99N15NEofmn6WPnaeg2V2OX 2kQSSfU1rqsu03DnC15qfhtJRMc0LY+8RsPqB583qf4hiSDr9s4pPxqqjpYQ np0k4drGOVgRRMN3rbfjv/eShO3jUWntuzQ0TdAeau4nif2201VR4TRMVM/S +T5AEpt6RPBkDA3b3FPPnyNJYlWq5FTlAxoSv66cXzpEEqyzsjmrH1H9sBk8 mjdMErSVSi7RiTSsWfB5hd4oSQhxVTRmk2n4KZP/d9IYSUy+1Ow49ZyGv/NC zk+MkwTbfU18VToN15zwGdKZJInWtRsP6GRR/w9o33OQRxLfx/SlYl/TEMvJ kFNTJPHx7fbPc29oeNWrK8N6miTe+JjctM+jITNeIlN3hvKXnoX+1wIa8s+W hfD/I4kH/w5Nri2mYZGsn0k2xXeJ49lxpTSUVijuN56l/Hb9tBN/BQ1X8KTt vlF8ztBFzaGKun9ocaHuHEnYC3m1f/tGw9p4v9EQiuWZlQ3TP2i4vGRfdjnF 3zYqVWv8pOHzuPZtHIpvmLuXWTbSsPLCkvRZijc4fSq40UxD1xPc3mmK2f5y uS/baDi+6NRkF8UJj5zTW/5Q64ndbn1Psfn74qei3VQ/z3hFX6RYuHbRgw19 NOyS2a+mTnHe4JmIkwM0rE+yCiKo87oIFwaGcWlIP1360ZDipSza9YIhGgYf 6655TdVbr2vvzR6lIcdI9p0kxUEWea6LJ2l4ObTw3EFKL33nhae3TdOwzNpc OJjSczjA1sZ9loaWfuZuaZTeyY/fWCbw0/FMtEl2LtUPqzwx06p5dIxp+vUl neqXRJ3Ntsn5dHTFpKLQCZLwnD9/rflCOj5q+LhBjuq3hvJhzavSdJzvlZJH jJBEy6aXyukydNSX+iNvTvllm8tBKWEFan1i5twmyl+Tt17MX8ugo8DXHy5h HJJIfzI3e0yZjp6fO3bVUH6U+ZnCea9Ox410zWS5PpLo0pv87LyOjgpmWrS3 v0kixtKkJO4/OupeO7vaqZWaH66P35dvomPeckndBc0kkZO4M1V5Kx0r62Mk 6VSeborG+v80pePg/cTZVVQeN6pwfAX20ZFkF8YbUHkd3LzVU3s/HX/axekY fqLy4dZvG2hNx4ebNE1FCJJQbdDFLafpeKHmpVfiK5Jo/Bvy39mzdBw/6Sq6 OJMk7oh1aEU707Go1in5QhqV9y1BSkPn6LhCwGxuJpEkSpN+TSVdo+O/jMLL S6j5c9LjUq5EHB0bu1N6Il0p/+i05S5/SMfu+iCXzrMkITaMb3Y8puoxWcll 2JNElrvw22vP6PirwLD9yBGSmHELfzeaTcf9B707R41I4r5ryoeWKjpuVjVe HqNAEme0xfJ53+hYu/HH47eLSGIz6ZQvW0PHbbqZhqWSJNHpolNg1kjHtNxO uWRBkljtUlBY2kVHkcs8V382l6h0qi1O/0fHwvlzy0pzqPm/ckNJJZ8MFsTZ K77L4BJugzElPYIymLR+2eXYZ1xC1snmI0tMBl10R+6qxnAJW8f+0khZGfzx 8ZDrkYtcYubMXLmvtgyauO9RPajHJbRPr/hmckIG9/Zavox7ySHazjBdM07J oH1OxF2fZA5xx4kutdBBBvet7t5t+IBD9LtP76t2kcFP8bwL0QEc4qlv1c99 vjJ4+p3LOcMjHGJxhGPrgWgZ3Lpx0cMsQQ4xl/9i8PhXGXx6bZHOEsNBIpN4 FFL8XQYDHarIExsGiSMl97SW1skgu3dPWPTyQeJ9+SW3zl8y+KpxR0KZ5CBx rsZ41K5PBsV6a/+uax4gensHps8KLsJ3UlW51i4DRDVtpbiX7iIk84KrJAPY VL6c+x9sXoS9sbTDgxfYRGZ5ZnkJLMKordGX3p1hE3d919yU2rEIz3/rj2Ca sgmTng3/XpgvQj0pRRlNOpv4/H4r2eG4CHP56KYb4/uJT8cP1polLMKmtYTw pcQ+4kPmjQcrBWWxr3jDh9LrPYSwiaV8tLAsiivMc93r0kPs61eP5ReTxcfe bWWV1j1Er8rX+/VSsng15v2DkLU9BD12ccQ1hizuPxnvaNHVTTj6pd+u05XF LcbMg10G3YSCVb2Xr5ssmpnL7PUb7SQu8WuafWmVxRMuW1fKrfpDiLVc3OmS sxgNk+xqTx2tJ7qfKUc+8l2CewMWPrsdWUpEasjNe28ph2FTdxd0Nz+HsDj7 8hMsedz8OXRml081WEce/SbRK4/LDbILW+2bYXJbMu9RoQKeKd6xWMy9E3w/ PPIrLVZAq6L6/TM3O4FPJ25Bf6kCCh0v9++I6YT5zFBlnSoFlLrwMiyoqBNk ed67ShsUUC/8+qiRZBeszTCJ7furgLK2u484p3WBy6LRjWtUFFHV3sRfsqEb hoO5JZZqiviFzTfh0N8NXgL9JheXK+Jcj2hHznQ3+JKtxz9qK+InrdFxFVYP BFeW37bcrIhHJmb+HHXogdTLDxp9LBUxfJ7F/P/Ge6Cz28C75JYi7voq9lxn fh8YPPqWqBmsiKP/YsU15frgsZX114i7iui/rH5YckUf2FS5Lz0ZpYj/ojwU Uvb0wa+sx18Enyrid26NE39UH9RcmmVtL1REyWnVCWtGP5RI51WWjyriB++J MGMVNihXbR/VnlTENNcHR6Z12HDt5g9GzLQitkR/UosxYIP+RN85BwElFHnh FBhzgg15LYsZolJKOFezZ8/rBDZkp3h57F6hhC2fPTYYLRqAxC1rFL6eUMLV MWtPJ4wMQMkeMb/pU0rYHqvnzZsbgI5jnezlDkoY9sK01kBiEFRuROUFuCih kaz/zseqg5BUxjuEl5Rw8/KzVv77ByF5b0nM60glfLLU0XcoaxDKjj/g64hW wvVvNBNX5Q9Cj7ung9QDJQz62xd3sHwQ1CLVNjo/obhpVb9n6yA8awxqUMtU wm1F4+vrRDiQesJCNq5cCVs3nnvEPMKBCo+VVyo+K2F9tSix3p4DfX5CPeNf ldB3HW2HrhsHNJ69zbWoU0LmMi13kZsceM5WsFzwRwn1djFDytI4kHau6961 KSWMOqIemj/KgQx/L+nTqxh4G8wqTl/kwuDmozf8VzPwzuGHA5PXubBydNvI 07UMrEl8duNCIBfSbGV+tusy8OMKgrc5hgup8DrGegcD96jphfblcCGJ95e5 5xgD1dYtu1/I5kJXVmOYky0Dz4brDKgOcUHFoYg/2I6BfMvs3LwnuJBYH9pZ 7shAjd8d0z2CJDx6rZWy1ZuB0pnHM/8okBDn7KS1PpyBGxgyx3fvJOGX6v5H FpEMfHFDUPCGKQlyLXpSHtEMTLn69+EzcxJidosPZz5k4K1/Ji/fHCEhSj0t VyONgcK9RwaGXEmoawtXM8pg4L3pdbwwTxJk7vtE22UxcHOwca+iDwkR83b5 JL6h6tNac2j6Oglhf3o3K35kYEK75a6aCBK+x1RnbCpj4OJqo4WN90mQNHvL OFTJwOdrniR8jiPhbmEA3/1qBko4ziw+n0jCnQfLSiVbGHilkX9VTxYJX8wX rNf6zcBzbUES5jkkiIuNJZt0MPBLRGhy2lsSgrxLbwX2UXpGj04rF5Jwy/Kk seAYA3ur5eYt/0xC4t6WI3aTDOx8v/jVhq8k5O+ydCmfZuCfgvMJq76TMLzF KDxYgIlffg6uafhJgsTGokSOEBPtZsdrAxtJUNfRzdkrysRtzjszlzWTcERt RQNdiomEr/U8uXYSLrCS+rzoTHS8G+br1kHpIa801SDLRP6qyWWvukgoXyip FK/ExIeLZoW4fSR0iNzWmmUxMV7uni6bTcKMAD8eV2Xigw27Y6oHSdCZGLZV XcFEZaOSH5Z/STAdcvQM0GLif43/no8PkeAw0OXft4aJQJilXhsh4Ub30ejd 65lYrvakamSUhIe/61PTNzKxLKKbZjZOwrumvXkLNzPR+I/q5agJEmpqK6rc gImbjc+Llk6SwPm6tbXGgIlDp/jfNvNIEKnI467bwcQGxuzNlikSVErW8kXv YuLspxyPsmkS9PPTpXkmTIyZd+Vq7AwJVm+WqRw2Y6JXR0ma5T8SzmUlrCuw YKLLp5mRKYrvpi3ewTzIxI/uudaBsySkJoUdvG7NxHNit9r45kj4mCB6pvMo E/2lFa+eoLgt5sbF7SeYqO+gvfkFxbyI6eCUU0y02XtnURvFMiGeD0UdmNQ3 1ULxKYq1bnEyHR2Z2P3yieIcxbuu2xd9dWGiRqDBbup7G05d+v1D24OJcdn/ 7pVSfMXLqjPci4mhxn95NymOdf0xOuLNxJkaI18tinPO7BY+4MvEk+oGCsXU eatPflz87ioTnwSsathMMfvo5uXyfkx8FeKQ9YSqb55V7ibfACZ6bLRNISk9 mOZaJm2BTNzzJzZfjeJNJilHMYSJbHz0dwel3/4dLNfEMCa+dqZtM6P0dcXY a/MiKT2an741oPQP2kSLsI9mYp21jAmD6k/yuuCnFXFM/LZkvUAn1c8irXm5 mglMVEjo/RU2RsLY0vEGbhL1fF/TRMowCVJKrv1mqUx856elL035RXNx31T2 Cyaah+u/OEmScFz8l9KFbCbyrrjMfhwg4aKQhXZTLhOlxpf1/uin8j9XhXrv mZgpXjFd0UvC55GCk3MEE8X6MjM9OklY3/Lk+a2vTPTV+3o6p4naT8th8uV3 JvKdaRnVbSBh8KqWUWMtpU/5UNqLOhKSVD50afxiYr52XIxFNdU/xzrm5x4m 2rP8nud+JOFvfpzLcD8TT+cvzIwvovK+8ESBPIeJul/e9rsWkHAzm2PtOMJE JXOrhc1vSJjyFTQW5mNhjiQYlaaS0Jds7xHHz0KxhV/0FySTUP+tMlZLkIUa 2YlXDJ6QkM0K6zsgzMIX6ent7rEknC1VDEiRYGGvxIlzXbdJaBLfULxTjoUB plxVhVMklK2L7WuRZyE73LVh+THKD0dnJN0VWbjvWtk/VWtqvr0ssYllslC8 2S+7fS8JOy32zvQvo9bbUiXrsonKW5zDf8E6LDwq2FDWKU7Cs49VNsx1LAx7 FrTzgxAJkYNaAa/Xs7DW/qvX1TkuuOmP1TVvZGFJWlTH12EuaHRc91gJLLz9 qM/Tr5ELMcsfZHw1ZuH3lgxly8dcGB3ZkexnysIbRRHfrKn3yd7C4XjdvSxs eWW3ZHcYF4Qtdt9JNmdhh/71uO6rXPDy5Z29fIiFy6Zq+KVtuGBefVBj1Wmq njnPFn1ZLmTGCrI6HVhoJr4I/0hwQfTUy8WxZ1lo5f/qmbMgF4omhecLubBQ /mVszLa/HNBe+qa7xZOFE5JE3UAFByS8ZJ7eucHCzKMRkmrnOeAARXEGN6l6 kzsEnRw58FHUKWLSn4XJcRbu8cc5cDHh4/VTgSwsjjVuyN/Ngf5yj+Obw1jo WUdmzShxoEz+h9LgQxaeUVzbdLVoEK4WhcQYv2fhoY5ChR9jA8BuLYrk5bGQ 72fwmty+ATgwPRKWks/CqK/bYvybB2DlButAwSIWRhR9uDxQNACNL9Qufihj Yc3noKwvQQOgE00cXlHHQp3DOsoc+QHocR5iiJMsVBWqtBhcyQZTBcvkz6rK aC2w9S7K90FMl1bRYTVlXBJ5N+eDaB90ZYg0D6or48NuVXkNXi/4Yr6U5Apl FLnYU1XX0Avpdiq+FmuUsSf8e9+VyF4Qyxoyb9mijBsbVRhM8V4o33aXj3tA GVXDLYabON2w1aX8sHSgMvIupf9MDe6ER4q2ql/GlDHdt/LO9tvNoPox5f6e w0uxY5HzCcFN1RC2qkSzpnAperL4FhxPSQP7kXk/KjVVcMvTMF2vLaWE7ath VmCoCvaEys9FGdUTTjcsZnaHq+BOseoj23zrCS+LnAaJeypYPvHR5k9WPXF7 zOtu2H0V1NEXyuiRayAydXnT0Q9VEGVm5skPNBC84rn6Zy9UMN7s8gm+oCYi rEYipLRcBV/fV6jKeddCxD11dgiopJ7/KCXLZrcQTz2/Ge6sUsFbJbILRZVa iTeyYVNV31RwzOrJWcHrrUSLNd2h9qcKLm2PTug2aiPUO+UMOzpV8PRA5Tn9 H7+J1TkXGUndKnjR6sOuWb52YpN/M8+uVwV/KRfPmsm1EybqD7P62SoY/mPn gh+72gkPRxZjaEgFLb17snTT2gnfzdd52SMquNvvXtvaknbCf0FHneeYCm4/ uVNR+Vc7EfPyafDkpAr+Mxu6XCf6h0i8Nu903pQKVh+jXY5T/kO82Gdn4Duj gmAbZWCh+4fIWVqmpD+rgmkGmvmzZn+IwhE13tycCq7RDm6Id/hD/A9J0PiD "]]}, Annotation[#, "Charting`Private`Tag$4237#1"]& ]}, {}, {}}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->{True, True}, AxesLabel->{None, None}, AxesOrigin->{0, 0}, DisplayFunction->Identity, Frame->{{False, False}, {False, False}}, FrameLabel->{{None, None}, {None, None}}, FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}}, GridLines->{None, None}, GridLinesStyle->Directive[ GrayLevel[0.5, 0.4]], ImagePadding->All, Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({ (Identity[#]& )[ Part[#, 1]], (Identity[#]& )[ Part[#, 2]]}& ), "CopiedValueFunction" -> ({ (Identity[#]& )[ Part[#, 1]], (Identity[#]& )[ Part[#, 2]]}& )}}, PlotRange->{{0, 10.}, {-0.9999231952405686, 0.9995735822543441}}, PlotRangeClipping->True, PlotRangePadding->{{ Scaled[0.02], Scaled[0.02]}, { Scaled[0.05], Scaled[0.05]}}, Ticks->{Automatic, Automatic}]], "Output", CellChangeTimes->{ 3.7521763419737196`*^9, {3.7521763737745824`*^9, 3.752176414147539*^9}, 3.75217654338085*^9, 3.7521770510237017`*^9, {3.7521773015038404`*^9, 3.752177331926486*^9}, 3.7521828107931056`*^9, 3.752321459515777*^9, 3.7523271966668377`*^9, 3.752402057198332*^9, 3.7840227821650305`*^9, 3.7840265393709307`*^9, 3.784038265057227*^9, 3.784040554937485*^9, 3.7840409619303284`*^9, 3.784041004263873*^9, 3.784041092285056*^9, 3.7840411936145573`*^9},ExpressionUUID->"8bba002a-dd0c-44a5-82b4-\ 59d4ab8323b0"] }, Open ]], Cell[BoxData[ RowBox[{"(*", RowBox[{ RowBox[{ "\:0412\:0441\:0438\:0447\:043a\:0438", " ", "\:044f\:0432\:043d\:0438", " ", "\:043c\:0435\:0442\:043e\:0434\:0438", " ", "\:0441\:0435", " ", "\:043f\:0440\:043e\:0432\:0430\:043b\:044f\:0442"}], ",", " ", RowBox[{ "\:0437\:0430\:0449\:043e\:0442\:043e", " ", "\:0442\:0435", " ", "\:0441\:0430", " ", "\:0443\:0441\:043b\:043e\:0432\:043d\:043e", " ", RowBox[{ "\:0443\:0441\:0442\:043e\:0439\:0447\:0438\:0432\:0438", ".", " ", "\:0417\:0430"}], " ", "\:0434\:0430", " ", "\:043f\:043e\:043b\:0443\:0447\:0438\:043c", " ", "\:0434\:043e\:0431\:0440\:0438", " ", "\:0440\:0435\:0437\:0443\:043b\:0442\:0430\:0442\:0438\:0442\:0435", " ", "\:0441\:0442\:044a\:043f\:043a\:0430\:0442\:0430", " ", "\:0442\:0440\:044f\:0431\:0432\:0430", " ", "\:0434\:0430", " ", "\:0435", " ", "\:043c\:043d\:043e\:0433\:043e", " ", "\:043c\:0430\:043b\:043a\:0430"}], ",", " ", RowBox[{ "\:043a\:043e\:0435\:0442\:043e", " ", "\:0437\:0430\:0431\:0430\:0432\:044f", " ", "\:0437\:043d\:0430\:0447\:0438\:0442\:0435\:043b\:043d\:043e", " ", RowBox[{ "\:0438\:0437\:0447\:0438\:0441\:043b\:0435\:043d\:0438\:044f\:0442\:0430\ ", "."}]}]}], "*)"}]], "Input", CellChangeTimes->{{3.7840408647853937`*^9, 3.784040918665291*^9}},ExpressionUUID->"c388219c-3291-4e96-aeac-\ 83a37b1099d9"], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"h", "=", "0.00001"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"yExplicit", " ", "=", " ", RowBox[{"explicitEuler", "[", RowBox[{"f", ",", "h", ",", "0.", ",", "10.", ",", "0."}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"plots", "[", "yExplicit", "]"}]}], "Input", CellChangeTimes->{{3.7521758334236255`*^9, 3.7521758477734404`*^9}, { 3.75217592655151*^9, 3.75217598810041*^9}, {3.7521761047648745`*^9, 3.7521761121856413`*^9}, 3.752176157989519*^9, {3.7521761971483183`*^9, 3.752176301133994*^9}, {3.7521763957360415`*^9, 3.752176398958202*^9}, { 3.752176548134141*^9, 3.752176550599721*^9}},ExpressionUUID->"c38522ce-fa28-4bec-8893-\ f45de4ab0043"], Cell[BoxData["\<\"The step-size is too big!\"\>"], "Output", CellChangeTimes->{{3.7521762902901783`*^9, 3.7521763086617184`*^9}, 3.7521764075133233`*^9, {3.7521765681666994`*^9, 3.75217658845013*^9}, 3.7521770506758924`*^9, 3.752182810628197*^9, 3.752321458849157*^9, 3.7523271959162664`*^9, 3.75240205646929*^9, 3.784022780947961*^9, 3.7840265359997377`*^9, 3.784038264875729*^9, 3.7840405546906075`*^9, 3.784040975822291*^9, 3.784041015384451*^9, 3.784041102428203*^9, 3.7840412047841096`*^9},ExpressionUUID->"43f8f882-d275-401f-8abe-\ 5a509e92323e"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"h", "=", "0.001"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"yRK4", " ", "=", " ", RowBox[{"RK4", "[", RowBox[{"f", ",", "h", ",", "0.", ",", "10.", ",", "0."}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"plots", "[", "yRK4", "]"}]}], "Input", CellChangeTimes->{{3.7521764438895884`*^9, 3.7521764771875863`*^9}, { 3.752176600489264*^9, 3.752176605684318*^9}, {3.75217673527639*^9, 3.752176746837804*^9}, {3.7521770333337827`*^9, 3.752177035617481*^9}, { 3.7521772477844915`*^9, 3.7521772483471546`*^9}, {3.75217729145457*^9, 3.7521772944988337`*^9}},ExpressionUUID->"2e8a9dcf-e092-4f44-a644-\ a3d91f71cae9"], Cell[BoxData["\<\"The step-size is too big!\"\>"], "Output", CellChangeTimes->{{3.752177283566084*^9, 3.752177295276388*^9}, 3.7521773367747216`*^9, 3.7521828111369247`*^9, 3.7523214598785686`*^9, 3.752327197686269*^9, 3.75240205769736*^9, 3.7840227827280626`*^9, 3.7840265422590957`*^9, 3.784038265513874*^9, 3.784040555322753*^9, 3.7840409762340527`*^9, 3.7840410158032265`*^9, 3.7840411028919325`*^9, 3.784041205380767*^9},ExpressionUUID->"a3171f7f-647c-457e-b92a-\ 62e24a65188e"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"h", "=", "0.00001"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"yAB4", " ", "=", " ", RowBox[{"AB4", "[", RowBox[{"f", ",", "h", ",", "0.", ",", "10.", ",", "0."}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"plots", "[", "yAB4", "]"}]}], "Input", CellChangeTimes->{{3.7521771048100185`*^9, 3.75217712253391*^9}, { 3.752177219693497*^9, 3.7521772202421846`*^9}, {3.752177341714903*^9, 3.7521773730430365`*^9}},ExpressionUUID->"3b5f1f64-8219-4c2f-9630-\ 3382ece604c8"], Cell[BoxData["\<\"The step-size is too big!\"\>"], "Output", CellChangeTimes->{{3.7521771292740645`*^9, 3.752177133940403*^9}, { 3.752177207252591*^9, 3.752177225826999*^9}, 3.7521772715649133`*^9, { 3.752177344293442*^9, 3.7521773897555046`*^9}, 3.75218283514522*^9, 3.752321477646229*^9, 3.7523272447864046`*^9, 3.7524020790445814`*^9, 3.784022801353128*^9, 3.784026562292241*^9, 3.7840383377656765`*^9, 3.7840405730001783`*^9, 3.7840409947373734`*^9, 3.7840410348342237`*^9, 3.78404111441928*^9, 3.7840412238351116`*^9},ExpressionUUID->"885fe4fc-9d30-46e5-9f98-\ e49c5db20bed"] }, Open ]] }, Open ]] }, WindowSize->{1366, 685}, WindowMargins->{{-8, Automatic}, {Automatic, -8}}, FrontEndVersion->"11.2 for Microsoft Windows (64-bit) (September 10, 2017)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[579, 22, 1239, 17, 92, "Title",ExpressionUUID->"467d1ec9-6d84-4b2a-90de-bdb8eba1d886"], Cell[1821, 41, 2560, 67, 238, "Input",ExpressionUUID->"d0a40c0c-c125-4a61-9539-c7544c932fc8"], Cell[4384, 110, 3130, 83, 238, "Input",ExpressionUUID->"4a6b5a71-ee89-4044-925a-9e0eca7f52ca"], Cell[7517, 195, 4434, 116, 276, "Input",ExpressionUUID->"218c4638-48db-45f8-9652-b892485b01bc"], Cell[11954, 313, 6709, 191, 590, "Input",ExpressionUUID->"36837f17-bac1-40b1-b9f2-4ea41252841b"], Cell[18666, 506, 154, 3, 28, "Input",ExpressionUUID->"5625b212-3c6c-4992-af0f-9cca96159dfb"] }, Open ]], Cell[CellGroupData[{ Cell[18857, 514, 4650, 91, 1130, "Title",ExpressionUUID->"33bc2259-f836-4b21-81b1-d28abbbb91b9"], Cell[CellGroupData[{ Cell[23532, 609, 7966, 230, 542, "Input",ExpressionUUID->"49b1b934-b068-4cfd-b851-5abc5cdb112c"], Cell[31501, 841, 2225, 44, 529, "Output",ExpressionUUID->"da393554-3426-4f13-bc00-8005af2c961f"] }, Open ]], Cell[CellGroupData[{ Cell[33763, 890, 211, 5, 28, "Input",ExpressionUUID->"19b20860-11d6-4cda-ba18-208c62396fcc"], Cell[33977, 897, 427, 6, 32, "Output",ExpressionUUID->"a20753d2-d829-4324-94a5-699220be10c2"] }, Open ]], Cell[34419, 906, 8678, 238, 618, "Input",ExpressionUUID->"1dc66ca4-ec0e-4a3f-88c4-d6723b1da711"], Cell[CellGroupData[{ Cell[43122, 1148, 2221, 40, 67, "Input",ExpressionUUID->"43927774-48b4-4b09-b3a5-b0e0d403cb3b"], Cell[45346, 1190, 13157, 176, 451, "Output",ExpressionUUID->"5e613eb3-9ff3-474f-8465-8992faaac46a"] }, Open ]], Cell[CellGroupData[{ Cell[58540, 1371, 1479, 43, 67, "Input",ExpressionUUID->"a27a273f-311e-4fc4-94f2-2a24d322b436"], Cell[60022, 1416, 2340, 45, 529, "Output",ExpressionUUID->"1b958677-d76f-4868-96fb-90eef5c203a6"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[62411, 1467, 2519, 56, 384, "Title",ExpressionUUID->"4ed04c9b-35d8-4b0c-ad3b-76a639522d18"], Cell[CellGroupData[{ Cell[64955, 1527, 2368, 66, 200, "Input",ExpressionUUID->"52f186ac-f2de-44ba-a065-046acfe79cb1"], Cell[67326, 1595, 11794, 208, 186, "Output",ExpressionUUID->"87f05e1b-9fc3-4a81-b896-949e6f1f10b9"] }, Open ]], Cell[CellGroupData[{ Cell[79157, 1808, 10600, 273, 899, "Input",ExpressionUUID->"faac927f-4c4f-4fa7-9b67-07cf15a280a5"], Cell[89760, 2083, 5335, 103, 449, "Output",ExpressionUUID->"b7957176-ff6d-41cf-a911-23b6125aeca0"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[95144, 2192, 2050, 36, 316, "Title",ExpressionUUID->"4ed04c9b-35d8-4b0c-ad3b-76a639522d18"], Cell[97197, 2230, 3267, 79, 200, "Input",ExpressionUUID->"07d23880-93e2-4b4d-a3fc-b8906f99386d"], Cell[CellGroupData[{ Cell[100489, 2313, 2424, 50, 86, "Input",ExpressionUUID->"d6e139a0-90b9-4a15-a972-eed403abef6d"], Cell[102916, 2365, 15162, 266, 230, "Output",ExpressionUUID->"8bba002a-dd0c-44a5-82b4-59d4ab8323b0"] }, Open ]], Cell[118093, 2634, 1423, 29, 48, "Input",ExpressionUUID->"c388219c-3291-4e96-aeac-83a37b1099d9"], Cell[CellGroupData[{ Cell[119541, 2667, 736, 15, 67, "Input",ExpressionUUID->"c38522ce-fa28-4bec-8893-f45de4ab0043"], Cell[120280, 2684, 577, 8, 32, "Output",ExpressionUUID->"43f8f882-d275-401f-8abe-5a509e92323e"] }, Open ]], Cell[CellGroupData[{ Cell[120894, 2697, 684, 14, 67, "Input",ExpressionUUID->"2e8a9dcf-e092-4f44-a644-a3d91f71cae9"], Cell[121581, 2713, 504, 7, 32, "Output",ExpressionUUID->"a3171f7f-647c-457e-b92a-62e24a65188e"] }, Open ]], Cell[CellGroupData[{ Cell[122122, 2725, 537, 12, 67, "Input",ExpressionUUID->"3b5f1f64-8219-4c2f-9630-3382ece604c8"], Cell[122662, 2739, 607, 9, 32, "Output",ExpressionUUID->"885fe4fc-9d30-46e5-9f98-e49c5db20bed"] }, Open ]] }, Open ]] } ] *)