(* Content-type: application/mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 7.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 145, 7] NotebookDataLength[ 98190, 2772] NotebookOptionsPosition[ 91790, 2595] NotebookOutlinePosition[ 92134, 2610] CellTagsIndexPosition[ 92091, 2607] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell["\<\ \:041a\:0438\:043d\:0435\:043c\:0430\:0442\:0438\:0447\:0435\:0441\:043a\:0438\ \:0435 \:0443\:0440\:0430\:0432\:043d\:0435\:043d\:0438\:044f \:043f\:043e\ \:0441\:0442\:0443\:043f\:0430\:0442\:0435\:043b\:044c\:043d\:043e\:0433\:043e\ \:0434\:0432\:0438\:0436\:0435\:043d\:0438\:044f\ \>", "Title"], Cell[BoxData[ FormBox[ RowBox[{"vg", "=", RowBox[{"{", RowBox[{ RowBox[{"vgx", "[", "t", "]"}], ",", " ", RowBox[{"vgy", "[", "t", "]"}], ",", " ", RowBox[{"vgz", "[", "t", "]"}]}], "}"}]}], TraditionalForm]], "Input"], Cell[BoxData[ FormBox[ RowBox[{"r", "=", RowBox[{"{", RowBox[{ RowBox[{"l", "[", "t", "]"}], ",", " ", RowBox[{"h", "[", "t", "]"}], ",", " ", RowBox[{"z", "[", "t", "]"}]}], "}"}]}], TraditionalForm]], "Input"], Cell[BoxData[ FormBox[ RowBox[{"kur1", "=", RowBox[{"TableForm", "[", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"r", "[", RowBox[{"[", "1", "]"}], "]"}], ",", "t"}], "]"}], "==", RowBox[{"vg", "[", RowBox[{"[", "1", "]"}], "]"}]}], "]"}]}], TraditionalForm]], "Input"], Cell[BoxData[ RowBox[{"kur2", "=", RowBox[{"TableForm", "[", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"r", "[", RowBox[{"[", "2", "]"}], "]"}], ",", "t"}], "]"}], "==", RowBox[{"vg", "[", RowBox[{"[", "2", "]"}], "]"}]}], "]"}]}]], "Input"], Cell[BoxData[ RowBox[{"kur3", "=", RowBox[{"TableForm", "[", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"r", "[", RowBox[{"[", "3", "]"}], "]"}], ",", "t"}], "]"}], "==", RowBox[{"vg", "[", RowBox[{"[", "3", "]"}], "]"}]}], "]"}]}]], "Input"], Cell[BoxData[ RowBox[{"pr", "=", RowBox[{"{", RowBox[{"kur1", ",", " ", "kur2", ",", " ", "kur3"}], "}"}]}]], "Input"], Cell[BoxData[ RowBox[{"TableForm", "[", RowBox[{"m1", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Cos", "[", RowBox[{"\[Theta]", "[", "t", "]"}], "]"}], ",", RowBox[{"-", RowBox[{"Sin", "[", RowBox[{"\[Theta]", "[", "t", "]"}], "]"}]}], ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Sin", "[", RowBox[{"\[Theta]", "[", "t", "]"}], "]"}], ",", RowBox[{"Cos", "[", RowBox[{"\[Theta]", "[", "t", "]"}], "]"}], ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "0", ",", "1"}], "}"}]}], "}"}]}], "]"}]], "Input"], Cell[BoxData[ FormBox[ RowBox[{"\[IndentingNewLine]", RowBox[{"TableForm", "[", RowBox[{"m2", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Cos", "[", RowBox[{"\[Psi]", "[", "t", "]"}], "]"}], ",", " ", "0", ",", " ", RowBox[{"Sin", "[", RowBox[{"\[Psi]", "[", "t", "]"}], "]"}]}], "}"}], ",", " ", RowBox[{"{", RowBox[{"0", ",", " ", "1", ",", " ", "0"}], "}"}], ",", " ", RowBox[{"{", RowBox[{ RowBox[{"-", RowBox[{"Sin", "[", RowBox[{"\[Psi]", "[", "t", "]"}], "]"}]}], ",", " ", "0", ",", " ", RowBox[{"Cos", "[", RowBox[{"\[Psi]", "[", "t", "]"}], "]"}]}], "}"}]}], "}"}]}], "]"}]}], TraditionalForm]], "Input"], Cell[BoxData[ FormBox[ RowBox[{"TableForm", "[", RowBox[{"m3", "=", RowBox[{"m2", ".", "m1"}]}], "]"}], TraditionalForm]], "Input"], Cell[BoxData[ FormBox[ RowBox[{"TableForm", "[", RowBox[{"ReplaceAll", "[", RowBox[{"m3", ",", " ", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"\[Theta]", "[", "t", "]"}], "\[Rule]", "0"}], ",", " ", RowBox[{ RowBox[{"\[Psi]", "[", "t", "]"}], "\[Rule]", "0"}]}], "}"}]}], "]"}], "]"}], TraditionalForm]], "Input"], Cell[BoxData[ FormBox[ RowBox[{"rg", "=", RowBox[{"m3", ".", "vg"}]}], TraditionalForm]], "Input"], Cell[BoxData[ FormBox[ RowBox[{"\[IndentingNewLine]", RowBox[{"q1", "=", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"r", "[", RowBox[{"[", "1", "]"}], "]"}], ",", " ", "t"}], "]"}], "\[Equal]", RowBox[{"rg", "[", RowBox[{"[", "1", "]"}], "]"}]}]}]}], TraditionalForm]], "Input"], Cell[BoxData[ FormBox[ RowBox[{"\[IndentingNewLine]", RowBox[{"q2", "=", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"r", "[", RowBox[{"[", "2", "]"}], "]"}], ",", " ", "t"}], "]"}], "\[Equal]", RowBox[{"rg", "[", RowBox[{"[", "2", "]"}], "]"}]}]}]}], TraditionalForm]], "Input"], Cell[BoxData[ FormBox[ RowBox[{"q3", "=", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"r", "[", RowBox[{"[", "3", "]"}], "]"}], ",", " ", "t"}], "]"}], "\[Equal]", RowBox[{"rg", "[", RowBox[{"[", "3", "]"}], "]"}]}]}], TraditionalForm]], "Input"], Cell[BoxData[ FormBox[ RowBox[{"\[IndentingNewLine]", RowBox[{"TableForm", "[", RowBox[{"{", RowBox[{"q1", ",", " ", "q2", ",", " ", "q3"}], "}"}], "]"}]}], TraditionalForm]], "Input"] }, Open ]], Cell[CellGroupData[{ Cell["\<\ \:0414\:0438\:043d\:0430\:043c\:0438\:0447\:0435\:0441\:043a\:0438\:0435 \ \:0443\:0440\:0430\:0432\:043d\:0435\:043d\:0438\:044f \:043f\:043e\:0441\ \:0442\:0443\:043f\:0430\:0442\:0435\:043b\:044c\:043d\:043e\:0433\:043e \ \:0434\:0432\:0438\:0436\:0435\:043d\:0438\:044f \ \>", "Title"], Cell[BoxData[ FormBox[ RowBox[{"v", "=", RowBox[{"{", RowBox[{ RowBox[{"vx", "[", "t", "]"}], ",", " ", RowBox[{"vy", "[", "t", "]"}], ",", " ", RowBox[{"vz", "[", "t", "]"}]}], "}"}]}], TraditionalForm]], "Input"], Cell[BoxData[ FormBox[ RowBox[{"w", "=", RowBox[{"{", RowBox[{ RowBox[{"wx", "[", "t", "]"}], ",", " ", RowBox[{"wy", "[", "t", "]"}], ",", " ", RowBox[{"wz", "[", "t", "]"}]}], "}"}]}], TraditionalForm]], "Input"], Cell["\<\ (*\:041c\:0430\:0442\:0440\:0438\:0446\:0430 \:043f\:0435\:0440\:0435\:0445\ \:043e\:0434\:0430 \:043e\:0442 \:0441\:0432\:044f\:0437\:0430\:043d\:043d\ \:043e\:0439 \:0441\:0438\:0441\:0442\:0435\:043c\:044b \:043a\:043e\:043e\ \:0440\:0434\:0438\:043d\:0430\:0442 \:043a \:043d\:043e\:0440\:043c\:0430\ \:043b\:044c\:043d\:043e\:0439 \:0441\:0438\:0441\:0442\:0435\:043c\:0435 \ \:043a\:043e\:043e\:0440\:0434\:0438\:043d\:0430\:0442*)\ \>", "Text"], Cell[BoxData[ RowBox[{"TableForm", "[", RowBox[{"mm1", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Cos", "[", RowBox[{"psi", "[", "t", "]"}], "]"}], ",", " ", "0", ",", " ", RowBox[{"Sin", "[", RowBox[{"psi", "[", "t", "]"}], "]"}]}], "}"}], ",", " ", RowBox[{"{", RowBox[{"0", ",", " ", "1", ",", " ", "0"}], "}"}], ",", " ", RowBox[{"{", RowBox[{ RowBox[{"-", RowBox[{"Sin", "[", RowBox[{"psi", "[", "t", "]"}], "]"}]}], ",", " ", "0", ",", " ", RowBox[{"Cos", "[", RowBox[{"psi", "[", "t", "]"}], "]"}]}], "}"}]}], "}"}]}], "]"}]], "Input"], Cell[BoxData[ RowBox[{"TableForm", "[", RowBox[{"mm2", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Cos", "[", RowBox[{"teta", "[", "t", "]"}], "]"}], ",", " ", RowBox[{"-", RowBox[{"Sin", "[", RowBox[{"teta", "[", "t", "]"}], "]"}]}], ",", " ", "0"}], "}"}], ",", " ", RowBox[{"{", RowBox[{ RowBox[{"Sin", "[", RowBox[{"teta", "[", "t", "]"}], "]"}], ",", " ", RowBox[{"Cos", "[", RowBox[{"teta", "[", "t", "]"}], "]"}], ",", " ", "0"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"0", ",", " ", "0", ",", " ", "1"}], "}"}]}], "}"}]}], "]"}]], "Input"], Cell[BoxData[ RowBox[{"TableForm", "[", RowBox[{"mm3", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", RowBox[{"Cos", "[", RowBox[{"gam", "[", "t", "]"}], "]"}], ",", RowBox[{"-", RowBox[{"Sin", "[", RowBox[{"gam", "[", "t", "]"}], "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", RowBox[{"Sin", "[", RowBox[{"gam", "[", "t", "]"}], "]"}], ",", RowBox[{"Cos", "[", RowBox[{"gam", "[", "t", "]"}], "]"}]}], "}"}]}], "}"}]}], "]"}]], "Input"], Cell[BoxData[ RowBox[{"TableForm", "[", RowBox[{"mm4", "=", RowBox[{"Transpose", "[", RowBox[{"mm3", ".", "mm2", ".", "mm1"}], "]"}]}], "]"}]], "Input"], Cell["\<\ (*\:041c\:0430\:0442\:0440\:0438\:0446\:0430 \:043f\:0435\:0440\:0435\:0445\ \:043e\:0434\:0430 \:043e\:0442 \:0441\:043a\:043e\:0440\:043e\:0441\:0442\ \:043d\:043e\:0439 \:0441\:0438\:0441\:0442\:0435\:043c\:044b \:043a\:043e\ \:043e\:0440\:0434\:0438\:043d\:0430\:0442 \:043a \:0441\:0432\:044f\:0437\ \:0430\:043d\:043d\:043e\:0439 \:0441\:0438\:0441\:0442\:0435\:043c\:0435 \ \:043a\:043e\:043e\:0440\:0434\:0438\:043d\:0430\:0442*)\ \>", "Text"], Cell[BoxData[ RowBox[{ RowBox[{"TableForm", "[", RowBox[{"sm1", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Cos", "[", RowBox[{"alfa", "[", "t", "]"}], "]"}], ",", RowBox[{"-", RowBox[{"Sin", "[", RowBox[{"alfa", "[", "t", "]"}], "]"}]}], ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Sin", "[", RowBox[{"alfa", "[", "t", "]"}], "]"}], ",", RowBox[{"Cos", "[", RowBox[{"alfa", "[", "t", "]"}], "]"}], ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "0", ",", "1"}], "}"}]}], "}"}]}], "]"}], " "}]], "Input"], Cell[BoxData[ RowBox[{"TableForm", "[", RowBox[{"sm2", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Cos", "[", RowBox[{"beta", "[", "t", "]"}], "]"}], ",", "0", ",", RowBox[{"Sin", "[", RowBox[{"beta", "[", "t", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", RowBox[{"Sin", "[", RowBox[{"beta", "[", "t", "]"}], "]"}]}], ",", "0", ",", RowBox[{"Cos", "[", RowBox[{"beta", "[", "t", "]"}], "]"}]}], "}"}]}], "}"}]}], "]"}]], "Input"], Cell[BoxData[ RowBox[{"TableForm", "[", RowBox[{"sm3", "=", RowBox[{"Transpose", "[", RowBox[{"sm2", ".", "sm1"}], "]"}]}], "]"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"Cy", "[", RowBox[{ RowBox[{"alfa", "[", "t", "]"}], ",", RowBox[{"mah", "[", "t", "]"}], ",", "ugstrel", ",", RowBox[{"dzak", "[", "t", "]"}], ",", RowBox[{"dh", "[", "t", "]"}]}], "]"}], "=", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"Cy", "[", RowBox[{ RowBox[{"alfa", "[", "t", "]"}], ",", RowBox[{"mah", "[", "t", "]"}], ",", "ugstrel", ",", RowBox[{"dzak", "[", "t", "]"}], ",", RowBox[{"dv", "[", "t", "]"}]}], "]"}], ",", RowBox[{"alfa", "[", "t", "]"}]}], "]"}], "*", RowBox[{"(", RowBox[{ RowBox[{"alfa", "[", "t", "]"}], "-", "alfa0"}], ")"}], " ", RowBox[{"(*", RowBox[{ "\:041a\:043e\:044d\:0444\:0444\:0438\:0446\:0438\:0435\:043d\:0442", " ", "\:043f\:043e\:0434\:044a\:0451\:043c\:043d\:043e\:0439", " ", "\:0441\:0438\:043b\:044b"}], "*)"}]}]}]], "Input"], Cell[BoxData[ RowBox[{"lyam", "=", RowBox[{ RowBox[{"(", RowBox[{"l", "^", "2"}], ")"}], "/", "S"}], " ", RowBox[{"(*", RowBox[{ "\:041e\:0442\:043d\:043e\:0441\:0438\:0442\:0435\:043b\:044c\:043d\:043e\ \:0435", " ", "\:0443\:0434\:043b\:0438\:043d\:0435\:043d\:0438\:0435", " ", "\:043a\:0440\:044b\:043b\:0430"}], "*)"}]}]], "Input"], Cell[BoxData[ RowBox[{"Cx", "=", RowBox[{"Cx0", "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"Cy", "[", RowBox[{ RowBox[{"alfa", "[", "t", "]"}], ",", RowBox[{"mah", "[", "t", "]"}], ",", "ugstrel", ",", RowBox[{"dzak", "[", "t", "]"}], ",", RowBox[{"dv", "[", "t", "]"}]}], "]"}], "^", "2"}], ")"}], "/", RowBox[{"(", RowBox[{"lyam", "*", "Pi"}], ")"}]}], " ", RowBox[{"(*", RowBox[{ "\:041a\:043e\:044d\:0444\:0444\:0438\:0446\:0438\:0435\:043d\:0442", " ", "\:0441\:0438\:043b\:044b", " ", "\:043b\:043e\:0431\:043e\:0432\:043e\:0433\:043e", " ", "\:0441\:043e\:043f\:0440\:043e\:0442\:0438\:0432\:043b\:0435\:043d\:0438\ \:044f"}], "*)"}]}]}]], "Input"], Cell[BoxData[ RowBox[{"q", "=", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"ro", "[", RowBox[{"h", "[", "t", "]"}], "]"}], "*", RowBox[{"vk", "^", "2"}]}], ")"}], "/", "2"}], " ", RowBox[{"(*", RowBox[{ "\:0421\:043a\:043e\:0440\:043e\:0441\:0442\:043d\:043e\:0439", " ", "\:043d\:0430\:043f\:043e\:0440"}], "*)"}]}]], "Input"], Cell[BoxData[ RowBox[{"faer", "=", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "Cx"}], "*", "q", "*", "s"}], ",", RowBox[{ RowBox[{"Cy", "[", RowBox[{ RowBox[{"alfa", "[", "t", "]"}], ",", RowBox[{"mah", "[", "t", "]"}], ",", "ugstrel", ",", RowBox[{"dzak", "[", "t", "]"}], ",", RowBox[{"dv", "[", "t", "]"}]}], "]"}], "*", "q", "*", "s"}], ",", RowBox[{ RowBox[{"Cz", "[", RowBox[{ RowBox[{"beta", "[", "t", "]"}], ",", RowBox[{"dh", "[", "t", "]"}]}], "]"}], "*", "q", "*", "s"}]}], "}"}], " ", RowBox[{"(*", RowBox[{ "\:0410\:044d\:0440\:043e\:0434\:0438\:043d\:0430\:043c\:0438\:0447\:0435\ \:0441\:043a\:0430\:044f", " ", "\:0441\:0438\:043b\:0430", " ", "\:0432", " ", "\:0441\:043a\:043e\:0440\:043e\:0441\:0442\:043d\:043e\:0439", " ", "\:0441\:0438\:0441\:0442\:0435\:043c\:0435", " ", "\:043a\:043e\:043e\:0440\:0434\:0438\:043d\:0430\:0442"}], "*)"}]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"TableForm", "[", RowBox[{"faersv", "=", RowBox[{"sm3", ".", "faer"}]}], "]"}], RowBox[{"(*", RowBox[{ "\:0410\:044d\:0440\:043e\:0434\:0438\:043d\:0430\:043c\:0438\:0447\:0435\ \:0441\:043a\:0430\:044f", " ", "\:0441\:0438\:043b\:0430", " ", "\:0432", " ", "\:0441\:0432\:044f\:0437\:0430\:043d\:043d\:043e\:0439", " ", "\:0441\:0438\:0441\:0442\:0435\:043c\:0435", " ", "\:043a\:043e\:043e\:0440\:0434\:0438\:043d\:0430\:0442"}], "*)"}]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"TableForm", "[", RowBox[{"gs", "=", RowBox[{ RowBox[{"Transpose", "[", "mm4", "]"}], ".", RowBox[{"{", RowBox[{"0", ",", RowBox[{"-", "g"}], ",", "0"}], "}"}]}]}], "]"}], " ", RowBox[{"(*", RowBox[{ "\:0421\:0438\:043b\:0430", " ", "\:0442\:044f\:0436\:0435\:0441\:0442\:0438", " ", "\:0432", " ", "\:0441\:0432\:044f\:0437\:0430\:043d\:043d\:043e\:0439", " ", "\:0441\:0438\:0441\:0442\:0435\:043c\:0435", " ", "\:043a\:043e\:043e\:0440\:0434\:0438\:043d\:0430\:0442"}], "*)"}]}]], "Input"], Cell[BoxData[ RowBox[{"fx", "=", RowBox[{ RowBox[{"faersv", "[", RowBox[{"[", "1", "]"}], "]"}], "+", RowBox[{"gs", "[", RowBox[{"[", "1", "]"}], "]"}], "+", RowBox[{ RowBox[{"p", "[", RowBox[{ RowBox[{"dsg", "[", "t", "]"}], ",", "v", ",", RowBox[{"h", "[", "t", "]"}]}], "]"}], "*", RowBox[{"Cos", "[", "phi", "]"}]}]}]}]], "Input"], Cell[BoxData[ RowBox[{"fy", "=", RowBox[{ RowBox[{"faersv", "[", RowBox[{"[", "2", "]"}], "]"}], "+", RowBox[{"gs", "[", RowBox[{"[", "2", "]"}], "]"}], "+", RowBox[{ RowBox[{"p", "[", RowBox[{ RowBox[{"dsg", "[", "t", "]"}], ",", RowBox[{"vk", "[", "t", "]"}], ",", RowBox[{"h", "[", "t", "]"}]}], "]"}], "*", RowBox[{"Sin", "[", "phi", "]"}]}]}]}]], "Input"], Cell[BoxData[ RowBox[{"fz", "=", RowBox[{ RowBox[{"faersv", "[", RowBox[{"[", "3", "]"}], "]"}], "+", RowBox[{"gs", "[", RowBox[{"[", "3", "]"}], "]"}]}]}]], "Input"], Cell[BoxData[ RowBox[{"f", "=", RowBox[{"{", RowBox[{"fx", ",", "fy", ",", "fz"}], "}"}]}]], "Input"], Cell[BoxData[ FormBox[ RowBox[{"f1", "=", RowBox[{"{", RowBox[{ RowBox[{"fx1", "[", RowBox[{ RowBox[{"dv", "[", "t", "]"}], ",", RowBox[{"dh", "[", "t", "]"}], ",", RowBox[{"de", "[", "t", "]"}], ",", RowBox[{"dsg", "[", "t", "]"}]}], "]"}], ",", " ", RowBox[{"fy1", "[", RowBox[{ RowBox[{"dv", "[", "t", "]"}], ",", RowBox[{"dh", "[", "t", "]"}], ",", RowBox[{"de", "[", "t", "]"}], ",", RowBox[{"dsg", "[", "t", "]"}]}], "]"}], ",", " ", RowBox[{"fz1", "[", RowBox[{ RowBox[{"dv", "[", "t", "]"}], ",", RowBox[{"dh", "[", "t", "]"}], ",", RowBox[{"de", "[", "t", "]"}], ",", RowBox[{"dsg", "[", "t", "]"}]}], "]"}]}], "}"}]}], TraditionalForm]], "Input"], Cell[BoxData[ FormBox[ RowBox[{"\[IndentingNewLine]", RowBox[{"d", "=", RowBox[{"m", "*", RowBox[{"(", RowBox[{ RowBox[{"D", "[", RowBox[{"v", ",", " ", "t"}], "]"}], "+", RowBox[{"Cross", "[", RowBox[{"w", ",", " ", "v"}], "]"}]}], ")"}]}]}]}], TraditionalForm]], "Input"], Cell[BoxData[ FormBox[ RowBox[{"q4", "=", RowBox[{ RowBox[{"d", "[", RowBox[{"[", "1", "]"}], "]"}], "==", RowBox[{"f", "[", RowBox[{"[", "1", "]"}], "]"}]}]}], TraditionalForm]], "Input"], Cell[BoxData[ FormBox[ RowBox[{"\[IndentingNewLine]", RowBox[{"q5", "=", RowBox[{ RowBox[{"d", "[", RowBox[{"[", "2", "]"}], "]"}], "==", RowBox[{"f", "[", RowBox[{"[", "2", "]"}], "]"}]}]}]}], TraditionalForm]], "Input"], Cell[BoxData[ FormBox[ RowBox[{"\[IndentingNewLine]", RowBox[{"q6", "=", RowBox[{ RowBox[{"d", "[", RowBox[{"[", "3", "]"}], "]"}], "==", RowBox[{"f", "[", RowBox[{"[", "3", "]"}], "]"}]}]}]}], TraditionalForm]], "Input"], Cell[BoxData[ FormBox[ RowBox[{"\[IndentingNewLine]", RowBox[{"TableForm", "[", RowBox[{"{", RowBox[{"q4", ",", " ", "q5", ",", " ", "q6"}], "}"}], "]"}]}], TraditionalForm]], "Input"] }, Open ]], Cell[CellGroupData[{ Cell["\<\ \:041a\:0438\:043d\:0435\:043c\:0430\:0442\:0438\:0447\:0435\:0441\:043a\:0438\ \:0435 \:0443\:0440\:0430\:0432\:043d\:0435\:043d\:0438\:044f \:0432\:0440\ \:0430\:0449\:0430\:0442\:0435\:043b\:044c\:043d\:043e\:0433\:043e \:0434\ \:0432\:0438\:0436\:0435\:043d\:0438\:044f \ \>", "Title"], Cell[BoxData[ FormBox[ RowBox[{"\[IndentingNewLine]", RowBox[{"ug", "=", RowBox[{"{", RowBox[{ RowBox[{"psi", "[", "t", "]"}], ",", " ", RowBox[{"teta", "[", "t", "]"}], ",", " ", RowBox[{"gam", "[", "t", "]"}]}], "}"}]}]}], TraditionalForm]], "Input"], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{"w", "=", RowBox[{"{", RowBox[{ RowBox[{"wx", "[", "t", "]"}], ",", " ", RowBox[{"wy", "[", "t", "]"}], ",", " ", RowBox[{"wz", "[", "t", "]"}]}], "}"}]}]}]], "Input"], Cell[BoxData[ FormBox[ RowBox[{"\[IndentingNewLine]", RowBox[{"(*", RowBox[{"TableForm", "[", RowBox[{"mm4", "=", RowBox[{"mm3", ".", "mm2", ".", "mm1"}]}], "]"}], "*)"}]}], TraditionalForm]], "Input"], Cell[CellGroupData[{ Cell["\<\ \:0421\:043f\:0440\:043e\:0435\:0446\:0438\:0440\:0443\:0435\:043c \:0443\ \:0433\:043b\:043e\:0432\:044b\:0435 \:0441\:043a\:043e\:0440\:043e\:0441\ \:0442\:0438 \:0440\:044b\:0441\:043a\:0430\:043d\:044c\:044f, \:0442\:0430\ \:043d\:0433\:0430\:0436\:0430 \:0438 \:043a\:0440\:0435\:043d\:0430 \:043d\ \:0430 \:043e\:0441\:0438 \:0441\:0432\:044f\:0437\:0430\:043d\:043d\:043e\ \:0439 \:0441.\:043a.\ \>", "Subsubtitle", CellChangeTimes->{{3.4970722216083984`*^9, 3.497072265201172*^9}}], Cell[BoxData[ RowBox[{"prpsi", "=", RowBox[{ RowBox[{"Transpose", "[", RowBox[{"mm2", ".", "mm3"}], "]"}], ".", RowBox[{"{", RowBox[{"0", ",", RowBox[{"D", "[", RowBox[{ RowBox[{"psi", "[", "t", "]"}], ",", "t"}], "]"}], ",", "0"}], "}"}]}]}]], "Input", CellChangeTimes->{{3.4970723837089844`*^9, 3.4970724387285156`*^9}, { 3.4970725976191406`*^9, 3.4970726448808594`*^9}}], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{"prteta", "=", RowBox[{ RowBox[{"Transpose", "[", "mm3", "]"}], ".", RowBox[{"{", RowBox[{"0", ",", "0", ",", RowBox[{"D", "[", RowBox[{ RowBox[{"teta", "[", "t", "]"}], ",", "t"}], "]"}]}], "}"}]}]}]}]], "Input", CellChangeTimes->{{3.497072758939453*^9, 3.49707277871875*^9}}], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{"prgam", "=", RowBox[{"{", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"gam", "[", "t", "]"}], ",", "t"}], "]"}], ",", "0", ",", "0"}], "}"}]}]}]], "Input", CellChangeTimes->{{3.4970727853554688`*^9, 3.4970727986240234`*^9}}], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{"pr", "=", RowBox[{"prgam", "+", "prpsi", "+", "prteta"}]}]}]], "Input", CellChangeTimes->{{3.497072909881836*^9, 3.4970729462128906`*^9}}], Cell[BoxData[ RowBox[{"q7", "=", RowBox[{ RowBox[{"w", "[", RowBox[{"[", "1", "]"}], "]"}], "==", RowBox[{"pr", "[", RowBox[{"[", "1", "]"}], "]"}]}]}]], "Input", CellChangeTimes->{{3.497072515258789*^9, 3.4970725713554688`*^9}, 3.497072659428711*^9, {3.4970729607060547`*^9, 3.4970729921777344`*^9}}], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{"q8", "=", RowBox[{ RowBox[{"w", "[", RowBox[{"[", "2", "]"}], "]"}], "==", RowBox[{"pr", "[", RowBox[{"[", "2", "]"}], "]"}]}]}]}]], "Input", CellChangeTimes->{{3.4970728045791016`*^9, 3.4970728351689453`*^9}, { 3.4970729731845703`*^9, 3.4970729759833984`*^9}}], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{"q9", "=", RowBox[{ RowBox[{"w", "[", RowBox[{"[", "3", "]"}], "]"}], "==", RowBox[{"pr", "[", RowBox[{"[", "3", "]"}], "]"}]}]}]}]], "Input", CellChangeTimes->{{3.497072839439453*^9, 3.4970728727626953`*^9}, { 3.49707297828125*^9, 3.4970729808779297`*^9}}], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{"TableForm", "[", RowBox[{"{", RowBox[{"q7", ",", " ", "q8", ",", " ", "q9"}], "}"}], "]"}]}]], "Input"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["\<\ \:0414\:0438\:043d\:0430\:043c\:0438\:0447\:0435\:0441\:043a\:0438\:0435 \ \:0443\:0440\:0430\:0432\:043d\:0435\:043d\:0438\:044f \:0432\:0440\:0430\ \:0449\:0430\:0442\:0435\:043b\:044c\:043d\:043e\:0433\:043e \:0434\:0432\ \:0438\:0436\:0435\:043d\:0438\:044f \ \>", "Title", CellChangeTimes->{{3.4970730647041016`*^9, 3.497073078645508*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"TableForm", "[", RowBox[{"TI", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"Ixx", ",", RowBox[{"-", "Ixy"}], ",", RowBox[{"-", "Ixz"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "Iyx"}], ",", "Iyy", ",", RowBox[{"-", "Iyz"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "Izx"}], ",", RowBox[{"-", "Izy"}], ",", "Izz"}], "}"}]}], "}"}]}], "]"}], " ", RowBox[{"(*", RowBox[{ "\:0422\:0435\:043d\:0437\:043e\:0440", " ", "\:0438\:043d\:0435\:0440\:0446\:0438\:0438"}], "*)"}], "\[IndentingNewLine]"}]], "Input", CellChangeTimes->{{3.497073102783203*^9, 3.4970731236416016`*^9}, { 3.497073161123047*^9, 3.497073162623047*^9}}], Cell["\<\ \:0422.\:043a. OZ \:043f\:0435\:0440\:043f\:0435\:043d\:0434\:0438\:043a\ \:0443\:043b\:044f\:0440\:043d\:0430 \:043f\:043b\:043e\:0441\:043a\:043e\ \:0441\:0442\:0438 \:0441\:0438\:043c\:043c\:0435\:0442\:0440\:0438\:0438 \ XOY, \:0442\:043e\ \>", "Text", CellChangeTimes->{{3.4970732330478516`*^9, 3.4970732528828125`*^9}}], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{"Ixz", "=", RowBox[{"Izx", "=", RowBox[{"Iyz", "=", RowBox[{"Izy", "=", "0", " "}]}]}]}]}]], "Input", CellChangeTimes->{{3.4970731367734375`*^9, 3.4970731370976562`*^9}, 3.4970732316464844`*^9}], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{"Iyx", "=", "Ixy"}]}]], "Input", CellChangeTimes->{{3.497073147826172*^9, 3.497073148583008*^9}}], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{"TableForm", "[", "TI", "]"}]}]], "Input", CellChangeTimes->{{3.4970731754121094`*^9, 3.4970731954052734`*^9}}], Cell[CellGroupData[{ Cell["\<\ \:0413\:043b\:0430\:0432\:043d\:044b\:0439 \:043c\:043e\:043c\:0435\:043d\ \:0442 \:043a\:043e\:043b\:0438\:0447\:0435\:0441\:0442\:0432 \:0434\:0432\ \:0438\:0436\:0435\:043d\:0438\:044f \:0432 \:0441\:0432\:044f\:0437\:0430\ \:043d\:043d\:043e\:0439 \:0441\:0438\:0441\:0442\:0435\:043c\:0435 \:043a\ \:043e\:043e\:0440\:0434\:0438\:043d\:0430\:0442:\ \>", "Subsubtitle", CellChangeTimes->{{3.4970732178271484`*^9, 3.4970732256503906`*^9}}], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{"TableForm", "[", RowBox[{"K", "=", RowBox[{"TI", ".", "w"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.497073185026367*^9, 3.4970732160976562`*^9}}], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{"M", "=", RowBox[{"{", RowBox[{ RowBox[{"Mx", "[", RowBox[{ RowBox[{"alfa", "[", "t", "]"}], ",", RowBox[{"mah", "[", "t", "]"}], ",", " ", RowBox[{"dv", "[", "t", "]"}], ",", RowBox[{"dh", "[", "t", "]"}], ",", RowBox[{"de", "[", "t", "]"}], ",", RowBox[{"dsg", "[", "t", "]"}]}], "]"}], ",", RowBox[{"My", "[", RowBox[{ RowBox[{"alfa", "[", "t", "]"}], ",", RowBox[{"mah", "[", "t", "]"}], ",", " ", RowBox[{"dv", "[", "t", "]"}], ",", RowBox[{"dh", "[", "t", "]"}], ",", RowBox[{"de", "[", "t", "]"}], ",", RowBox[{"dsg", "[", "t", "]"}]}], "]"}], ",", RowBox[{"Mz", "[", RowBox[{ RowBox[{"alfa", "[", "t", "]"}], ",", " ", RowBox[{"mah", "[", "t", "]"}], ",", " ", RowBox[{"dv", "[", "t", "]"}], ",", RowBox[{"dh", "[", "t", "]"}], ",", RowBox[{"de", "[", "t", "]"}], ",", RowBox[{"dsg", "[", "t", "]"}]}], "]"}]}], "}"}]}]}]], "Input", CellChangeTimes->{{3.4970732691376953`*^9, 3.497073269541992*^9}}] }, Open ]], Cell[CellGroupData[{ Cell["\<\ \:041f\:043e\:043b\:043d\:0430\:044f \:043f\:0440\:043e\:0438\:0437\:0432\ \:043e\:0434\:043d\:0430\:044f \:043f\:043e \:0432\:0440\:0435\:043c\:0435\ \:043d\:0438 \:0433\:043b\:0430\:0432\:043d\:043e\:0433\:043e \:0432\:0435\ \:043a\:0442\:043e\:0440\:0430 \:043c\:043e\:043c\:0435\:043d\:0442\:0430 \ \:043a\:043e\:043b\:0438\:0447\:0435\:0441\:0442\:0432 \:0434\:0432\:0438\ \:0436\:0435\:043d\:0438\:044f\ \>", "Subsubtitle", CellChangeTimes->{{3.4970733012558594`*^9, 3.497073301758789*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"TableForm", "[", RowBox[{"pdtgv", "=", RowBox[{ RowBox[{"D", "[", RowBox[{"K", ",", "t"}], "]"}], "+", RowBox[{"Cross", "[", RowBox[{"w", ",", "K"}], "]"}]}]}], "]"}], "\[IndentingNewLine]"}]], "Input", CellChangeTimes->{{3.497073305866211*^9, 3.497073356705078*^9}}], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{"q10", "=", RowBox[{ RowBox[{"pdtgv", "[", RowBox[{"[", "1", "]"}], "]"}], "==", RowBox[{"M", "[", RowBox[{"[", "1", "]"}], "]"}]}]}]}]], "Input", CellChangeTimes->{{3.4970733664179688`*^9, 3.49707341978125*^9}}], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{"q11", "=", RowBox[{ RowBox[{"pdtgv", "[", RowBox[{"[", "2", "]"}], "]"}], "==", RowBox[{"M", "[", RowBox[{"[", "2", "]"}], "]"}]}]}]}]], "Input", CellChangeTimes->{{3.497073430461914*^9, 3.4970734459541016`*^9}}], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{"q12", "=", RowBox[{ RowBox[{"pdtgv", "[", RowBox[{"[", "3", "]"}], "]"}], "==", RowBox[{"M", "[", RowBox[{"[", "3", "]"}], "]"}]}]}]}]], "Input", CellChangeTimes->{{3.497073442741211*^9, 3.497073452970703*^9}}], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{"TableForm", "[", RowBox[{"{", RowBox[{"q10", ",", " ", "q11", ",", " ", "q12"}], "}"}], "]"}]}]], "Input"], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"MatrixForm", "[", RowBox[{"urav", "=", RowBox[{"{", RowBox[{ "q1", ",", " ", "q2", ",", " ", "q3", ",", " ", "q4", ",", " ", "q5", ",", " ", "q6", ",", " ", "q7", ",", " ", "q8", ",", " ", "q9", ",", " ", "q10", ",", " ", "q11", ",", " ", "q12"}], "}"}]}], "]"}]}]], "Input"], Cell[BoxData[ RowBox[{"vxq", "=", RowBox[{"Flatten", "[", RowBox[{"{", RowBox[{ RowBox[{"D", "[", RowBox[{"r", ",", "t"}], "]"}], ",", " ", RowBox[{"D", "[", RowBox[{"v", ",", "t"}], "]"}], ",", " ", RowBox[{"D", "[", RowBox[{"w", ",", "t"}], "]"}], ",", " ", RowBox[{"D", "[", RowBox[{"ug", ",", "t"}], "]"}]}], "}"}], "]"}]}]], "Input"], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{"vq", "=", RowBox[{"Flatten", "[", RowBox[{"{", RowBox[{"r", ",", "v", ",", " ", "w", ",", " ", "ug"}], "}"}], "]"}]}]}]], "Input"], Cell[BoxData[ RowBox[{"MatrixForm", "[", RowBox[{"fk", "=", RowBox[{"Flatten", "[", RowBox[{"Solve", "[", RowBox[{"urav", ",", "vxq", ",", " ", RowBox[{"Sort", "\[Rule]", "False"}]}], "]"}], "]"}]}], "]"}]], "Input"], Cell[BoxData[ RowBox[{"TableForm", "[", RowBox[{"uravn", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"fk", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "\[Equal]", RowBox[{"fk", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}]}], ",", " ", RowBox[{"{", RowBox[{"i", ",", "1", ",", "12"}], "}"}]}], "]"}]}], "]"}]], "Input"], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{"MatrixForm", "[", RowBox[{"matrA", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"uravn", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], ",", RowBox[{"vq", "[", RowBox[{"[", "j", "]"}], "]"}]}], "]"}], ",", " ", RowBox[{"{", RowBox[{"i", ",", "1", ",", "12"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"j", ",", "1", ",", "12"}], "}"}]}], "]"}]}], "]"}]}]], "Input"], Cell[BoxData[ RowBox[{"upr", "=", RowBox[{"{", RowBox[{ RowBox[{"dv", "[", "t", "]"}], ",", RowBox[{"dh", "[", "t", "]"}], ",", " ", RowBox[{"de", "[", "t", "]"}], ",", " ", RowBox[{"dsg", "[", "t", "]"}]}], "}"}]}]], "Input"], Cell[BoxData[ RowBox[{"MatrixForm", "[", RowBox[{"matrB", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"uravn", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], ",", " ", RowBox[{"upr", "[", RowBox[{"[", "j", "]"}], "]"}]}], "]"}], ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "1", ",", " ", "12"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"j", ",", " ", "1", ",", " ", "4"}], "}"}]}], "]"}]}], "]"}]], "Input"], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{"MatrixForm", "[", RowBox[{"matper", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1"}], "}"}]}], "}"}]}], "]"}]}]], "Input"], Cell[BoxData[""], "Input"], Cell[BoxData[ RowBox[{"MatrixForm", "[", RowBox[{"prodbok", "=", RowBox[{"matper", ".", "fk"}]}], "]"}]], "Input"], Cell[BoxData[ RowBox[{"MatrixForm", "[", RowBox[{"mrazdA", "=", RowBox[{"matper", ".", "matrA"}]}], "]"}]], "Input"], Cell[BoxData[ RowBox[{"MatrixForm", "[", RowBox[{"mrazdB", "=", RowBox[{"matper", ".", "matrB"}]}], "]"}]], "Input"], Cell[BoxData[ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "<", "7"}], ",", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", "7"}], ",", " ", RowBox[{"j", "<", "13"}], ",", RowBox[{ RowBox[{ RowBox[{"mrazdA", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], "=", "0"}], ";", RowBox[{"j", "++"}]}]}], "]"}], ";", RowBox[{"i", "++"}]}]}], "]"}]], "Input"], Cell[BoxData[ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "7"}], ",", RowBox[{"i", "<", "13"}], ",", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", "1"}], ",", " ", RowBox[{"j", "<", "7"}], ",", RowBox[{ RowBox[{ RowBox[{"mrazdA", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], "=", "0"}], ";", RowBox[{"j", "++"}]}]}], "]"}], ";", RowBox[{"i", "++"}]}]}], "]"}]], "Input"], Cell["\<\ \:041c\:0430\:0442\:0440\:0438\:0446\:0430 \:0410 \:0434\:043b\:044f \:043f\ \:0440\:043e\:0434\:043e\:043b\:044c\:043d\:043e\:0433\:043e \:0434\:0432\ \:0438\:0436\:0435\:043d\:0438\:044f:\ \>", "Text"], Cell[BoxData[ RowBox[{"MatrixForm", "[", RowBox[{"mApr", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"mrazdA", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "1", ",", " ", "6"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"j", ",", " ", "1", ",", " ", "6"}], "}"}]}], "]"}]}], "]"}]], "Input"], Cell["\<\ \:041c\:0430\:0442\:0440\:0438\:0446\:0430 \:0410 \:0434\:043b\:044f \:0431\ \:043e\:043a\:043e\:0432\:043e\:0433\:043e \:0434\:0432\:0438\:0436\:0435\ \:043d\:0438\:044f:\ \>", "Text"], Cell[BoxData[ RowBox[{"MatrixForm", "[", RowBox[{"mAb", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"mrazdA", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "7", ",", " ", "12"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"j", ",", " ", "7", ",", " ", "12"}], "}"}]}], "]"}]}], "]"}]], "Input"], Cell[BoxData[ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "<", "7"}], ",", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", "3"}], ",", " ", RowBox[{"j", "<", "5"}], ",", RowBox[{ RowBox[{ RowBox[{"mrazdB", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], "=", "0"}], ";", RowBox[{"j", "++"}]}]}], "]"}], ";", RowBox[{"i", "++"}]}]}], "]"}]], "Input"], Cell[BoxData[ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "7"}], ",", RowBox[{"i", "<", "13"}], ",", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", "1"}], ",", " ", RowBox[{"j", "<", "3"}], ",", RowBox[{ RowBox[{ RowBox[{"mrazdB", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], "=", "0"}], ";", RowBox[{"j", "++"}]}]}], "]"}], ";", RowBox[{"i", "++"}]}]}], "]"}]], "Input"], Cell["\<\ \:041c\:0430\:0442\:0440\:0438\:0446\:0430 B \:0434\:043b\:044f \:043f\:0440\ \:043e\:0434\:043e\:043b\:044c\:043d\:043e\:0433\:043e \:0434\:0432\:0438\ \:0436\:0435\:043d\:0438\:044f:\ \>", "Text"], Cell[BoxData[ RowBox[{"MatrixForm", "[", RowBox[{"mBpr", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"mrazdB", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "1", ",", " ", "6"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"j", ",", " ", "1", ",", " ", "2"}], "}"}]}], "]"}]}], "]"}]], "Input"], Cell["\<\ \:041c\:0430\:0442\:0440\:0438\:0446\:0430 B \:0434\:043b\:044f \:0431\:043e\ \:043a\:043e\:0432\:043e\:0433\:043e \:0434\:0432\:0438\:0436\:0435\:043d\ \:0438\:044f:\ \>", "Text"], Cell[BoxData[ RowBox[{"MatrixForm", "[", RowBox[{"mBb", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"mrazdB", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "6", ",", " ", "12"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"j", ",", " ", "3", ",", " ", "4"}], "}"}]}], "]"}]}], "]"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"MatrixForm", "[", RowBox[{"Xb", "=", RowBox[{"{", RowBox[{ RowBox[{"wx", "[", "t", "]"}], ",", RowBox[{"wy", "[", "t", "]"}], ",", RowBox[{"beta", "[", "t", "]"}], ",", RowBox[{"gam", "[", "t", "]"}], ",", RowBox[{"psi", "[", "t", "]"}]}], "}"}]}], "]"}], " ", RowBox[{"(*", RowBox[{ "\:0412\:0435\:043a\:0442\:043e\:0440", " ", "\:0441\:043e\:0441\:0442\:043e\:044f\:043d\:0438\:044f", " ", "\:0431\:043e\:043a\:043e\:0432\:043e\:0433\:043e", " ", "\:0434\:0432\:0438\:0436\:0435\:043d\:0438\:044f"}], "*)"}]}]], "Input"], Cell[BoxData[ RowBox[{"MatrixForm", "[", RowBox[{"Ub", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"dh", "[", "t", "]"}], "}"}], ",", RowBox[{"{", RowBox[{"de", "[", "t", "]"}], "}"}], ",", RowBox[{"{", "0", "}"}], ",", RowBox[{"{", "0", "}"}], ",", RowBox[{"{", "0", "}"}]}], "}"}]}], "]"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"a1", "=", "0.635"}], ";", RowBox[{"a2", "=", "5.47"}], ";", RowBox[{"a3", "=", "2.72"}], ";", RowBox[{"a4", "=", "0.269"}], ";", RowBox[{"a5", "=", "3.26"}], ";", RowBox[{"a6", "=", "0.709"}], ";", RowBox[{"a7", "=", "0.043"}], ";"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"b1", "=", "3.1"}], ";", RowBox[{"b2", "=", "20.2"}], ";", RowBox[{"b3", "=", "17.6"}], ";", RowBox[{"b4", "=", "0.072"}], ";", RowBox[{"b5", "=", RowBox[{"-", "0.518"}]}], ";", RowBox[{"b6", "=", "0.057"}], ";", RowBox[{"b7", "=", "0.065"}], ";"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"c1", "=", "0.82"}], ";", RowBox[{"c2", "=", "11.0"}], ";", RowBox[{"c3", "=", "8.8"}], ";", RowBox[{"c4", "=", "1.2"}], ";", RowBox[{"c5", "=", "0.3"}], ";", RowBox[{"c6", "=", "7.3"}], ";"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"MatrixForm", "[", RowBox[{"Ab", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "b1"}], ",", RowBox[{"-", "a6"}], ",", RowBox[{"-", "b2"}], ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "b6"}], ",", RowBox[{"-", "a1"}], ",", RowBox[{"-", "a2"}], ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"b7", ",", "a8", ",", RowBox[{"-", "a4"}], ",", "b4", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", RowBox[{"-", "a9"}], ",", "0", ",", "b8", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "a10", ",", "0", ",", RowBox[{"-", "b9"}], ",", "0"}], "}"}]}], "}"}]}], "]"}], RowBox[{"(*", RowBox[{ "\:041c\:0430\:0442\:0440\:0438\:0446\:0430", " ", "\:0410", " ", "\:0434\:043b\:044f", " ", "\:0431\:043e\:043a\:043e\:0432\:043e\:0433\:043e", " ", "\:0434\:0432\:0438\:0436\:0435\:043d\:0438\:044f"}], "*)"}]}]], "Input", CellChangeTimes->{{3.501904463591797*^9, 3.5019044773984375`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"MatrixForm", "[", RowBox[{"Bb", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "a5"}], ",", RowBox[{"-", "a3"}], ",", RowBox[{"-", "a7"}], ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "b3"}], ",", RowBox[{"-", "b5"}], ",", "0", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}]}], "}"}]}], "]"}], RowBox[{"(*", RowBox[{ "\:041c\:0430\:0442\:0440\:0438\:0446\:0430", " ", "\:0412", " ", "\:0434\:043b\:044f", " ", "\:0431\:043e\:043a\:043e\:0432\:043e\:0433\:043e", " ", "\:0434\:0432\:0438\:0436\:0435\:043d\:0438"}], "*)"}]}]], "Input", CellChangeTimes->{{3.5019044812265625`*^9, 3.5019044936904297`*^9}}], Cell[BoxData[ RowBox[{"dxb", "=", RowBox[{"D", "[", RowBox[{"Xb", ",", "t"}], "]"}]}]], "Input"], Cell[BoxData[ RowBox[{"lxb", "=", RowBox[{"LaplaceTransform", "[", RowBox[{"dxb", ",", "t", ",", "s"}], "]"}]}]], "Input"], Cell[BoxData[ RowBox[{"lxb", "=", RowBox[{"Xb", ".", RowBox[{"DiagonalMatrix", "[", RowBox[{"{", RowBox[{"s", ",", "s", ",", "s", ",", "s", ",", "s"}], "}"}], "]"}]}]}]], "Input"], Cell[BoxData[ RowBox[{"urB", "=", RowBox[{"MatrixForm", "[", RowBox[{"lxb", "-", RowBox[{"Ab", ".", "Xb"}]}], "]"}]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"MatrixForm", "[", RowBox[{"murB", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"urB", "[", RowBox[{"[", RowBox[{"1", ",", "i"}], "]"}], "]"}], ",", RowBox[{"Xb", "[", RowBox[{"[", "j", "]"}], "]"}]}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "5"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "5"}], "}"}]}], "]"}]}], "]"}], " ", RowBox[{"(*", RowBox[{ RowBox[{ "\:043c\:0430\:0442\:0440\:0438\:0446\:0430", " ", "\:043a\:043e\:044d\:0444\:0444"}], "-", RowBox[{"\:0432", " ", "\:0441\:0438\:0441\:0442\:0435\:043c\:044b", " ", RowBox[{"\:043b\:0438\:043d", ".", " ", "\:0443\:0440"}]}], "-", RowBox[{"\:0439", " ", "\:0432", " ", "\:043e\:043f\:0435\:0440"}], "-", RowBox[{ "\:0439", " ", "\:0444\:043e\:0440\:043c\:0435", " ", "\:0431\:043e\:043a\:043e\:0432\:043e\:0433\:043e", " ", "\:0434\:0432\:0438\:0436\:0435\:043d\:0438\:044f"}]}], "*)"}]}]], "Input"], Cell[BoxData[{ RowBox[{ RowBox[{"MurKr", "=", RowBox[{"Delete", "[", RowBox[{"murB", ",", RowBox[{"{", RowBox[{ RowBox[{"{", "2", "}"}], ",", RowBox[{"{", "3", "}"}], ",", RowBox[{"{", "5", "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "2"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "3"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "5"}], "}"}], ",", RowBox[{"{", RowBox[{"4", ",", "2"}], "}"}], ",", RowBox[{"{", RowBox[{"4", ",", "3"}], "}"}], ",", RowBox[{"{", RowBox[{"4", ",", "5"}], "}"}]}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"MatrixForm", "[", RowBox[{"MurKr", "=", RowBox[{"ReplaceAll", "[", RowBox[{"MurKr", ",", RowBox[{"{", RowBox[{"b8", "\[Rule]", "0"}], "}"}]}], "]"}]}], "]"}]}], "Input"], Cell[BoxData[ RowBox[{"Wkrenmx", "=", RowBox[{ RowBox[{"Det", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"MurKr", "[", RowBox[{"[", RowBox[{"1", ",", "1"}], "]"}], "]"}], ",", "57.3"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"MurKr", "[", RowBox[{"[", RowBox[{"2", ",", "1"}], "]"}], "]"}], ",", "0"}], "}"}]}], "}"}], "]"}], "/", RowBox[{"Det", "[", "MurKr", "]"}]}], " ", RowBox[{"(*", RowBox[{ "\:041f\:0435\:0440\:0435\:0434\:0430\:0442\:043e\:0447\:043d\:0430\:044f", " ", "\:0444\:0443\:043d\:043a\:0446\:0438\:044f", " ", RowBox[{"\:043a\:0440\:0435\:043d", "/", "Mx"}]}], "*)"}]}]], "Input"], Cell[BoxData[ RowBox[{"hkrenmx", "=", RowBox[{"InverseLaplaceTransform", "[", RowBox[{ RowBox[{"Wkrenmx", "/", "s"}], ",", "s", ",", "t"}], "]"}], RowBox[{"(*", RowBox[{ "\:041f\:0435\:0440\:0435\:0445\:043e\:0434\:043d\:0430\:044f", " ", "\:0444\:0443\:043d\:043a\:0446\:0438\:044f", " ", "\:0441\:0432\:043e\:0431\:043e\:0434\:043d\:043e\:0433\:043e", " ", "\:0441\:0430\:043c\:043e\:043b\:0435\:0442\:0430"}], "*)"}]}]], "Input", CellChangeTimes->{{3.5018751071152344`*^9, 3.5018751163251953`*^9}, { 3.5018752838828125`*^9, 3.501875287477539*^9}}], Cell[BoxData[ RowBox[{"(*", RowBox[{"\:0410\:041f", "-", "\:043a\:0440\:0435\:043d"}], "*)"}]], "Input"], Cell[BoxData[ RowBox[{"tr", "=", "4", " ", RowBox[{"(*", RowBox[{ "\:0412\:0440\:0435\:043c\:044f", " ", "\:0440\:0435\:0433\:0443\:043b\:0438\:0440\:043e\:0432\:0430\:043d\:0438\ \:044f"}], "*)"}]}]], "Input"], Cell[BoxData[ RowBox[{"(*", RowBox[{ "\:0421\:043e\:043e\:0442\:043d\:043e\:0448\:0435\:043d\:0438\:044f", " ", "\:0434\:043b\:044f", " ", "\:043f\:0435\:0440\:0435\:0434\:0430\:0442\:043e\:0447\:043d\:044b\:0445", " ", "\:0447\:0438\:0441\:0435\:043b", " ", "\:0410\:041f"}], "*)"}]], "Input"], Cell[BoxData[ RowBox[{"myue", "=", RowBox[{ RowBox[{"(", RowBox[{"18", "-", RowBox[{"b1", "*", "tr"}]}], ")"}], "/", RowBox[{"(", RowBox[{"b3", "*", "tr"}], ")"}]}]}]], "Input"], Cell[BoxData[ RowBox[{"ie", "=", RowBox[{"108", "/", RowBox[{"(", RowBox[{"b3", "*", "tr", "*", "tr"}], ")"}]}]}]], "Input"], Cell[BoxData[ RowBox[{"nyue", "=", RowBox[{"216", "/", RowBox[{"(", RowBox[{"b3", "*", "tr", "*", "tr", "*", "tr"}], ")"}]}]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"MatrixForm", "[", RowBox[{"MurKrap", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"(", RowBox[{"s", "+", "b1"}], ")"}], ",", "0", ",", "b3"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", "s", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"myue", ",", RowBox[{"(", RowBox[{"ie", "+", RowBox[{"nyue", "/", "s"}]}], ")"}], ",", RowBox[{"-", "1"}]}], "}"}]}], "}"}]}], "]"}], RowBox[{"(*", RowBox[{ "\:041c\:0430\:0442\:0440\:0438\:0446\:0430", " ", "\:0441\:043e\:0441\:0442\:043e\:044f\:043d\:0438\:044f", " ", "\:0441\:0432\:043e\:0431\:043e\:0434\:043d\:043e\:0433\:043e", " ", "\:0441\:0430\:043c\:043e\:043b\:0435\:0442\:0430"}], "*)"}]}]], "Input", CellChangeTimes->{{3.5018751341308594`*^9, 3.5018751482158203`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"MatrixForm", "[", RowBox[{"MurKrap0", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"(", RowBox[{"s", "+", "b1"}], ")"}], ",", "0", ",", "b3"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"myue", ",", RowBox[{"nyue", "/", "s"}], ",", RowBox[{"-", "1"}]}], "}"}]}], "}"}]}], "]"}], RowBox[{"(*", RowBox[{ "\:041c\:0430\:0442\:0440\:0438\:0446\:0430", " ", "\:0441\:043e\:0441\:0442\:043e\:044f\:043d\:0438\:044f", " ", "\:0432\:043e\:0437\:043c\:0443\:0449\:0435\:043d\:043d\:043e\:0433\:043e\ ", " ", "\:0441\:0430\:043c\:043e\:043b\:0435\:0442\:0430"}], "*)"}]}]], "Input", CellChangeTimes->{{3.501875151958008*^9, 3.5018751758320312`*^9}}], Cell[BoxData[ RowBox[{"Wgam", "=", RowBox[{ RowBox[{ RowBox[{"Det", "[", "MurKrap0", "]"}], "/", RowBox[{"Det", "[", "MurKrap", "]"}]}], RowBox[{"(*", RowBox[{ "\:0420\:0430\:0441\:0447\:0435\:0442", " ", "\:043f\:0435\:0440\:0435\:0434\:0430\:0442\:043e\:0447\:043d\:043e\:0439\ ", " ", "\:0444\:0443\:043d\:043a\:0446\:0438\:0438"}], "*)"}]}]}]], "Input", CellChangeTimes->{{3.501875185088867*^9, 3.5018751961152344`*^9}}], Cell[BoxData[ RowBox[{"hkrap", "=", RowBox[{"InverseLaplaceTransform", "[", RowBox[{ RowBox[{"Wgam", "/", "s"}], ",", "s", ",", "t"}], "]"}], RowBox[{"(*", RowBox[{ RowBox[{ "\:041f\:0435\:0440\:0435\:0445\:043e\:0434\:043d\:0430\:044f", " ", "\:0444\:0443\:043d\:043a\:0446\:0438\:044f", " ", "\:0441\:0438\:0441\:0442\:0435\:043c\:044b", " ", "\:041b\:0410"}], "-", "\:0410\:041f"}], "*)"}]}]], "Input", CellChangeTimes->{{3.5018752559208984`*^9, 3.501875280783203*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"Plot", "[", RowBox[{"hkrenmx", ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "10"}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{"0", ",", "100"}], "}"}]}]}], "]"}], " ", RowBox[{"(*", RowBox[{ "\:0421\:0432\:043e\:0431\:043e\:0434\:043d\:044b\:0439", " ", "\:0441\:0430\:043c\:043e\:043b\:0435\:0442"}], "*)"}]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"Plot", "[", RowBox[{"hkrap", ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "10"}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{"0", ",", "2"}], "}"}]}]}], "]"}], RowBox[{"(*", RowBox[{ RowBox[{ "\:0421\:0438\:0441\:0442\:0435\:043c\:0430", " ", "\:041b\:0410"}], "-", "\:0410\:041f\:043a\:0440\:0435\:043d\:0430"}], "*)"}]}]], "Input", CellChangeTimes->{3.501873076553711*^9}], Cell[BoxData[ RowBox[{ RowBox[{"ht", "[", "t_", "]"}], "=", RowBox[{"Abs", "[", "hkrap", "]"}]}]], "Input"], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"ksh", "=", "1000"}], ";"}], RowBox[{"(*", RowBox[{ "\:041a\:043e\:043b\:0438\:0447\:0435\:0441\:0442\:0432\:043e", " ", "\:0448\:0430\:043e\:0432"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{"dt", "=", RowBox[{ RowBox[{"10", "/", "ksh"}], RowBox[{"(*", RowBox[{ "\:0432\:0435\:043b\:0438\:0447\:0438\:043d\:0430", " ", "\:0448\:0430\:0433\:0430"}], "*)"}]}]}]}], "Input", CellChangeTimes->{{3.5018753292890625`*^9, 3.501875347734375*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"us", "=", RowBox[{"ht", "[", "1000", "]"}]}], ";", RowBox[{"(*", RowBox[{ "\:0423\:0441\:0442\:0430\:043d\:043e\:0432\:0438\:0432\:0448\:0435\:0435\ \:0441\:044f", " ", "\:0437\:043d\:0430\:0447\:0435\:043d\:0438\:0435", " ", "\:043f\:0435\:0440\:0435\:0445\:043e\:0434\:043d\:043e\:0433\:043e", " ", "\:043f\:0440\:043e\:0446\:0435\:0441\:0441\:0430"}], "*)"}], "\[IndentingNewLine]", RowBox[{"hmax", "=", RowBox[{"ht", "[", "0", "]"}]}], ";"}], RowBox[{"(*", RowBox[{ "\:041d\:0430\:0447\:0430\:043b\:044c\:043d\:043e\:0435", " ", "\:043c\:0430\:043a\:0441\:0438\:043c\:0430\:043b\:044c\:043d\:043e\:0435\ ", " ", "\:0437\:043d\:0430\:0447\:0435\:043d\:0438\:0435", " ", "\:0444\:0443\:043d\:043a\:0446\:0438\:0438"}], "*)"}]}]], "Input", CellChangeTimes->{{3.5018753532109375`*^9, 3.501875413692383*^9}}], Cell[BoxData[ RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{ "\:0427\:0438\:0441\:043b\:0435\:043d\:043d\:043e\:0435", " ", "\:0440\:0435\:0448\:0435\:043d\:0438\:0435", " ", "\:0443\:0440\:0430\:0432\:043d\:0435\:043d\:0438\:0439", " ", "h", RowBox[{"(", "t", ")"}]}], "=", RowBox[{ RowBox[{"0.95", " ", "\:0438", " ", "h", RowBox[{"(", "t", ")"}]}], "=", "1.05"}]}], ",", " ", RowBox[{ RowBox[{ RowBox[{ "\:043e\:0431\:044a\:0435\:0434\:0438\:043d\:0435\:043d\:043d\:044b\:0445\ ", " ", "\:0432", " ", "\:043e\:0434\:043d\:043e"}], " ", "|", RowBox[{ RowBox[{"h", RowBox[{"(", "t", ")"}]}], "-", "1"}], "|", RowBox[{"-", "0.05"}]}], "=", "0"}], ",", " ", RowBox[{ RowBox[{ "\:0442", ".", "\:0435", ".", " ", "\:043d\:0430\:0445\:043e\:0436\:0434\:0435\:043d\:0438\:0435"}], " ", "\:0432\:0440\:0435\:043c\:0435\:043d\:0438", " ", "\:0440\:0435\:0433\:0443\:043b\:0438\:0440\:043e\:0432\:0430\:043d\:0438\ \:044f"}], ",", " ", RowBox[{ "\:0438", " ", "\:0442\:0430\:043a\:0436\:0435", " ", "\:0440\:0430\:0441\:0447\:0435\:0442", " ", "\:043f\:0435\:0440\:0435\:0440\:0435\:0433\:0443\:043b\:0438\:0440\:043e\ \:0432\:0430\:043d\:0438\:044f"}]}], "*)"}]], "Input", CellChangeTimes->{{3.5018754763935547`*^9, 3.5018755326845703`*^9}}], Cell[BoxData[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "<", RowBox[{"ksh", "+", "1"}]}], ",", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Sign", "[", RowBox[{ RowBox[{"Abs", "[", RowBox[{ RowBox[{"ht", "[", RowBox[{"dt", "*", RowBox[{"(", RowBox[{"i", "-", "1"}], ")"}]}], "]"}], "-", "us"}], "]"}], "-", "0.05"}], "]"}], "\[NotEqual]", RowBox[{"Sign", "[", RowBox[{ RowBox[{"Abs", "[", RowBox[{ RowBox[{"ht", "[", RowBox[{"dt", "*", RowBox[{"(", "i", ")"}]}], "]"}], "-", "us"}], "]"}], "-", "0.05"}], "]"}]}], ",", RowBox[{"tr", "=", RowBox[{"dt", "*", RowBox[{"(", RowBox[{"i", "-", "0.5"}], ")"}]}]}]}], "]"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{"hmax", "<", RowBox[{"ht", "[", RowBox[{"dt", "*", "i"}], "]"}]}], ",", RowBox[{"hmax", "=", RowBox[{"ht", "[", RowBox[{"dt", "*", "i"}], "]"}]}]}], "]"}]}], "}"}], ";", RowBox[{"i", "++"}]}]}], "]"}], "\[IndentingNewLine]", RowBox[{"sig", "=", RowBox[{"Abs", "[", RowBox[{"100", "*", RowBox[{ RowBox[{"(", RowBox[{"hmax", "-", "us"}], ")"}], "/", "us"}]}], "]"}]}], "\[IndentingNewLine]", "tr"}], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"kk", "=", "4"}], ";"}]], "Input"], Cell[BoxData[ RowBox[{"step", "=", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"1.9", "-", "0.1"}], ")"}], "/", "2"}], RowBox[{"(*", RowBox[{ "\:0437\:0430\:0434\:0430\:0435\:043c", " ", "\:0448\:0430\:0433", " ", "\:0438\:0437\:043c\:0435\:043d\:0435\:043d\:0438\:044f", " ", "\:043f\:0435\:0440\:0435\:0434\:0430\:0442\:043e\:0447\:043d\:044b\:0445\ ", " ", "\:0447\:0438\:0441\:0435\:043b"}], "*)"}]}]}]], "Input", CellChangeTimes->{{3.501875575317383*^9, 3.5018755966621094`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"myun", "=", "1.2"}], ";", RowBox[{"nyun", "=", "1.8"}], ";", RowBox[{"in", "=", "3"}], ";"}], RowBox[{"(*", RowBox[{ "\:041d\:0430\:0447\:0430\:043b\:044c\:043d\:044b\:0435", " ", "\:0437\:043d\:0430\:0447\:0435\:043d\:0438\:044f", " ", "\:043f\:0435\:0440\:0435\:0434\:0430\:0442\:043e\:0447\:043d\:044b\:0445\ ", " ", "\:0447\:0438\:0441\:0435\:043b"}], "*)"}]}]], "Input", CellChangeTimes->{{3.5018755469248047`*^9, 3.5018755607529297`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"minTreg", "=", "100"}], ";"}]], "Input", CellChangeTimes->{{3.5018713925371094`*^9, 3.5018714023408203`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell["\<\ \:041c\:0435\:0442\:043e\:0434 \:043f\:043e\:043b\:043d\:043e\:0433\:043e \ \:043f\:043b\:0430\:043d\:0438\:0440\:043e\:0432\:0430\:043d\:0438\:044f \ \>", "Subsubtitle", CellChangeTimes->{{3.5018748883916016`*^9, 3.5018749090126953`*^9}}], Cell[BoxData[ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "<", "kk"}], ",", RowBox[{"(*", RowBox[{"\:0438\:0437\:043c\:0435\:043d\:0435\:043d\:0438\:0435", " ", RowBox[{ "\:043f\:0435\:0440\:0435\:0434", ".", " ", "\:0447\:0438\:0441\:043b\:0430"}], " ", "myu"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"{", RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", "1"}], ",", RowBox[{"j", "<", "kk"}], ",", RowBox[{"(*", RowBox[{"\:0438\:0437\:043c\:0435\:043d\:0435\:043d\:0438\:0435", " ", RowBox[{ "\:043f\:0435\:0440\:0435\:0434", ".", " ", "\:0447\:0438\:0441\:043b\:0430"}], " ", "nyu"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"{", RowBox[{"For", "[", RowBox[{ RowBox[{"k", "=", "1"}], ",", RowBox[{"k", "<", "kk"}], ",", RowBox[{"(*", RowBox[{ "\:0438\:0437\:043c\:0435\:043d\:0435\:043d\:0438\:0435", " ", RowBox[{ "\:043f\:0435\:0440\:0435\:0434", ".", " ", "\:0447\:0438\:0441\:043b\:0430"}], " ", "i"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"myu", "=", RowBox[{"myun", "*", RowBox[{"(", RowBox[{"0.1", "+", RowBox[{"step", "*", RowBox[{"(", RowBox[{"i", "-", "1"}], ")"}]}]}], ")"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"nyu", "=", RowBox[{"nyun", "*", RowBox[{"(", RowBox[{"0.1", "+", RowBox[{"step", "*", RowBox[{"(", RowBox[{"j", "-", "1"}], ")"}]}]}], ")"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"ii", "=", RowBox[{"in", "*", RowBox[{"(", RowBox[{"0.1", "+", RowBox[{"step", "*", RowBox[{"(", RowBox[{"k", "-", "1"}], ")"}]}]}], ")"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"M1", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"(", RowBox[{"s", "+", "b1"}], ")"}], ",", "0", ",", "b3"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", "s", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"myu", ",", RowBox[{"(", RowBox[{"ii", "+", RowBox[{"nyu", "/", "s"}]}], ")"}], ",", RowBox[{"-", "1"}]}], "}"}]}], "}"}]}], ";", RowBox[{"(*", RowBox[{ "\:041c\:0430\:0442\:0440\:0438\:0446\:0430", " ", "\:0441\:043e\:0441\:0442\:043e\:044f\:043d\:0438\:044f", " ", "\:0441\:0432\:043e\:0431\:043e\:0434\:043d\:043e\:0433\:043e\ ", " ", "\:0441\:0430\:043c\:043e\:043b\:0435\:0442\:0430"}], "*)"}], "\[IndentingNewLine]", RowBox[{"M0", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"(", RowBox[{"s", "+", "b1"}], ")"}], ",", "0", ",", "b3"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"myu", ",", RowBox[{"nyu", "/", "s"}], ",", RowBox[{"-", "1"}]}], "}"}]}], "}"}]}], ";", RowBox[{"(*", RowBox[{ "\:041c\:0430\:0442\:0440\:0438\:0446\:0430", " ", "\:0441\:043e\:0441\:0442\:043e\:044f\:043d\:0438\:044f", " ", "\:0432\:043e\:0437\:043c\:0443\:0449\:0435\:043d\:043d\:043e\ \:0433\:043e", " ", "\:0441\:0430\:043c\:043e\:043b\:0435\:0442\:0430"}], "*)"}], "\[IndentingNewLine]", FormBox[ RowBox[{"wn", "=", RowBox[{ RowBox[{"Det", "[", "M0", "]"}], "/", RowBox[{"Det", "[", "M1", "]"}]}]}], TraditionalForm], ";", RowBox[{"(*", RowBox[{ "\:041f\:0435\:0440\:0435\:0434\:0430\:0442\:043e\:0447\:043d\ \:0430\:044f", " ", "\:0444\:0443\:043d\:043a\:0446\:0438\:044f"}], "*)"}], "\[IndentingNewLine]", RowBox[{"hkrap", "=", RowBox[{"InverseLaplaceTransform", "[", RowBox[{ RowBox[{"wn", "/", "s"}], ",", "s", ",", "t"}], "]"}]}], ";", RowBox[{"(*", RowBox[{ "\:041f\:0435\:0440\:0435\:0445\:043e\:0434\:043d\:0430\:044f", " ", "\:0444\:0443\:043d\:043a\:0446\:0438\:044f"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"ht", "[", "t_", "]"}], "=", RowBox[{"Abs", "[", "hkrap", "]"}]}], ";", "\[IndentingNewLine]", FormBox[ RowBox[{ RowBox[{ RowBox[{"us", "=", RowBox[{"ht", "[", "10000", "]"}]}], ";"}], RowBox[{"(*", RowBox[{ "\:0437\:0430\:0434\:0430\:0435\:043c", " ", "\:0443\:0441\:0442\:0430\:043d\:043e\:0432\:0438\:0432\ \:0448\:0435\:0435\:0441\:044f", " ", "\:0437\:043d\:0430\:0447\:0435\:043d\:0438\:0435"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"hmax", "=", "1"}], ";"}], RowBox[{"(*", RowBox[{ "\:0437\:0430\:0434\:0430\:0435\:043c", " ", "\:043c\:0430\:043a\:0441\:0438\:043c\:0430\:043b\:044c\ \:043d\:043e\:0435", " ", "\:0437\:043d\:0430\:0447\:0435\:043d\:0438\:0435"}], "*)"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{ "\:0427\:0438\:0441\:043b\:0435\:043d\:043d\:043e\:0435", " ", "\:0440\:0435\:0448\:0435\:043d\:0438\:0435", " ", "\:0443\:0440\:0430\:0432\:043d\:0435\:043d\:0438\:0439", " ", RowBox[{"h", "(", "t", ")"}]}], "=", RowBox[{ RowBox[{"0.95", " ", "\:0438", " ", RowBox[{"h", "(", "t", ")"}]}], "=", "1.05"}]}], ",", " ", RowBox[{ RowBox[{ RowBox[{ "\:043e\:0431\:044a\:0435\:0434\:0438\:043d\:0435\:043d\ \:043d\:044b\:0445", " ", "\:0432", " ", "\:043e\:0434\:043d\:043e"}], " ", "|", RowBox[{ RowBox[{"h", "(", "t", ")"}], "-", "1"}], "|", RowBox[{"-", "0.05"}]}], "=", "0"}], ",", " ", RowBox[{ RowBox[{ "\:0442", ".", "\:0435", ".", " ", "\:043d\:0430\:0445\:043e\:0436\:0434\:0435\:043d\:0438\ \:0435"}], " ", "\:0432\:0440\:0435\:043c\:0435\:043d\:0438", " ", "\:0440\:0435\:0433\:0443\:043b\:0438\:0440\:043e\:0432\ \:0430\:043d\:0438\:044f"}], ",", " ", RowBox[{ "\:0438", " ", "\:0442\:0430\:043a\:0436\:0435", " ", "\:0440\:0430\:0441\:0447\:0435\:0442", " ", "\:043f\:0435\:0440\:0435\:0440\:0435\:0433\:0443\:043b\ \:0438\:0440\:043e\:0432\:0430\:043d\:0438\:044f"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"l", "=", "1"}], ",", RowBox[{"l", "<", RowBox[{"ksh", "+", "1"}]}], ",", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Sign", "[", RowBox[{ RowBox[{"Abs", "[", RowBox[{ RowBox[{"ht", "[", RowBox[{"dt", "*", RowBox[{"(", RowBox[{"l", "-", "1"}], ")"}]}], "]"}], "-", "us"}], "]"}], "-", "0.05"}], "]"}], "\[NotEqual]", RowBox[{"Sign", "[", RowBox[{ RowBox[{"Abs", "[", RowBox[{ RowBox[{"ht", "[", RowBox[{"dt", "*", RowBox[{"(", "l", ")"}]}], "]"}], "-", "us"}], "]"}], "-", "0.05"}], "]"}]}], ",", RowBox[{"trs", "=", RowBox[{"dt", "*", RowBox[{"(", RowBox[{"l", "-", "0.5"}], ")"}]}]}]}], "]"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{"hmax", "<", RowBox[{"ht", "[", RowBox[{"dt", "*", "i"}], "]"}]}], ",", RowBox[{"hmax", "=", RowBox[{"ht", "[", RowBox[{"dt", "*", "i"}], "]"}]}]}], "]"}]}], "}"}], ";", RowBox[{"l", "++"}]}]}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"trs", "<", "minTreg"}], ",", RowBox[{"minTreg", "=", "trs"}], ",", RowBox[{"sigm", "=", RowBox[{"Round", "[", RowBox[{"Abs", "[", RowBox[{"100", "*", RowBox[{ RowBox[{"(", RowBox[{"hmax", "-", "us"}], ")"}], "/", "us"}]}], "]"}], "]"}]}]}], "]"}], ";"}]}], TraditionalForm]}], RowBox[{"(*", RowBox[{ "\:041d\:0430\:0445\:043e\:0434\:0438\:043c", " ", "\:043c\:0438\:043d\:0438\:043c\:0430\:043b\:044c\:043d\:043e\ \:0435", " ", "\:0432\:0440\:0435\:043c\:044f", " ", "\:0440\:0435\:0433\:0443\:043b\:0438\:0440\:043e\:0432\:0430\ \:043d\:0438\:044f", " ", "\:0438", " ", "\:043f\:0435\:0440\:0435\:0440\:0435\:0433\:0443\:043b\:0438\ \:0440\:043e\:0432\:0430\:043d\:0438\:0435"}], "*)"}], "\[IndentingNewLine]", "}"}], ";", "\[IndentingNewLine]", RowBox[{"k", "++"}]}]}], "]"}], "}"}], ";", "\[IndentingNewLine]", RowBox[{"j", "++"}]}]}], "]"}], "}"}], ";", "\[IndentingNewLine]", RowBox[{"i", "++"}]}]}], "]"}]], "Input", CellChangeTimes->{{3.501863327841797*^9, 3.5018634624384766`*^9}, 3.5018636246533203`*^9, 3.501863740973633*^9, {3.5018639826933594`*^9, 3.501863986158203*^9}, {3.501864557310547*^9, 3.5018645668847656`*^9}, { 3.5018659076279297`*^9, 3.5018659617373047`*^9}, 3.50186605253125*^9, { 3.5018665430234375`*^9, 3.501866589272461*^9}, {3.501866790216797*^9, 3.5018667995634766`*^9}, {3.5018668390664062`*^9, 3.501866876404297*^9}, 3.501866939185547*^9, {3.5018670946748047`*^9, 3.501867114328125*^9}, { 3.5018674136679688`*^9, 3.5018674658828125`*^9}, {3.5018674964628906`*^9, 3.5018675167910156`*^9}, {3.50186813640625*^9, 3.501868154366211*^9}, { 3.501868458345703*^9, 3.501868459166992*^9}, 3.501868531026367*^9, { 3.501868673665039*^9, 3.5018686741279297`*^9}, {3.501868835801758*^9, 3.501868881044922*^9}, {3.501869006323242*^9, 3.5018690099345703`*^9}, { 3.501869161291992*^9, 3.5018691672177734`*^9}, {3.501869225665039*^9, 3.5018692279765625`*^9}, {3.501869281817383*^9, 3.5018692829160156`*^9}, { 3.501869419678711*^9, 3.501869421953125*^9}, {3.5018694595615234`*^9, 3.501869461260742*^9}, {3.5018695609560547`*^9, 3.501869582373047*^9}, 3.5018696608945312`*^9, 3.501869712821289*^9, {3.501869840354492*^9, 3.5018698715117188`*^9}, {3.5018699082666016`*^9, 3.5018699303564453`*^9}, {3.501869988213867*^9, 3.5018699896367188`*^9}, { 3.5018702714316406`*^9, 3.5018702759472656`*^9}, 3.5018703147890625`*^9, { 3.5018705865146484`*^9, 3.5018706219248047`*^9}, {3.5018707085751953`*^9, 3.5018707096728516`*^9}, 3.501870754705078*^9, {3.5018714090410156`*^9, 3.5018716765585938`*^9}, {3.5018717242851562`*^9, 3.5018717434257812`*^9}, {3.501871832088867*^9, 3.501871833428711*^9}, { 3.5018756124765625`*^9, 3.5018756977148438`*^9}, {3.50187573321875*^9, 3.501875807538086*^9}, {3.5018758380341797`*^9, 3.501875901276367*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"myuoptim", "=", RowBox[{ RowBox[{"(", RowBox[{"18", "-", RowBox[{"b1", "*", "tr"}]}], ")"}], "/", RowBox[{"(", RowBox[{"b3", "*", "minTreg"}], ")"}]}]}], ";"}], RowBox[{"(*", RowBox[{ "\:043e\:043f\:0442\:0438\:043c\:0430\:043b\:044c\:043d\:044b\:0435", " ", "\:0437\:043d\:0430\:0447\:0435\:043d\:0438\:044f", " ", "\:043f\:0435\:0440\:0435\:0434\:0430\:0442\:043e\:0447\:043d\:044b\:0445\ ", " ", "\:0447\:0438\:0441\:0435\:043b"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"nyuoptim", "=", RowBox[{"216", "/", RowBox[{"(", RowBox[{"b3", "*", RowBox[{"minTreg", "^", "3"}]}], ")"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"ioptim", "=", RowBox[{"108", "/", RowBox[{"(", RowBox[{"b3", "*", RowBox[{"minTreg", "^", "2"}]}], ")"}]}]}], ";"}]}], "Input", CellChangeTimes->{{3.5018726404785156`*^9, 3.5018727579140625`*^9}, { 3.501875912083008*^9, 3.5018759410810547`*^9}}], Cell[BoxData[ RowBox[{"Print", "[", RowBox[{ "\"\<\:0423\:0447\:0438\:0442\:044b\:044b\:0432\:0430\:044f \:0437\:043d\ \:0430\:0447\:0435\:043d\:0438\:044f \:0432\:0440\:0435\:043c\:0435\:043d\ \:0438 \:0440\:0435\:0433\:0443\:043b\:0438\:0440\:043e\:0432\:0430\:043d\ \:0438\:044f \:0438 \:043f\:0435\:0440\:0435\:0440\:0435\:0433\:0443\:043b\ \:0438\:0440\:043e\:0432\:0430\:043d\:0438\:044f, \:043f\:043e\:043b\:0443\ \:0447\:0430\:0435\:043c, \:0447\:0442\:043e \:043d\:0430\:0438\:0431\:043e\ \:043b\:0435\:0435 \:043e\:043f\:0442\:0438\:043c\:0430\:043b\:044c\:043d\ \:044b\:0435 \:0437\:043d\:0430\:0447\:0435\:043d\:0438\:044f \:043f\:0435\ \:0440\:0435\:0434\:0430\:0442\:043e\:0447\:043d\:044b\:0445 \:0447\:0438\ \:0441\:0435\:043b \:0430\:0432\:0442\:043e\:043f\:0438\:043b\:043e\:0442\ \:0430 \:0440\:0430\:0432\:043d\:044b: myu=\>\"", ",", " ", "myuoptim", ",", "\"\<, nyu=\>\"", ",", "nyuoptim", ",", "\"\<, i=\>\"", ",", "ioptim", ",", "\"\<. \:041f\:0440\:0438 \:044d\:0442\:043e\:043c \:0432\:0440\:0435\ \:043c\:044f \:0440\:0435\:0433\:0443\:043b\:0438\:0440\:043e\:0432\:0430\ \:043d\:0438\:044f \:0441\:043e\:0441\:0442\:0430\:0432\:043b\:044f\:0435\ \:0442 \>\"", ",", " ", "minTreg", ",", "\"\\"", ",", "sigm", ",", "\"\<%.\>\""}], "]"}]], "Input", CellChangeTimes->{{3.5018711045273438`*^9, 3.5018711048095703`*^9}, { 3.501871873276367*^9, 3.501871917333008*^9}, {3.501872559294922*^9, 3.501872635073242*^9}, {3.5018727677158203`*^9, 3.501872834397461*^9}, { 3.50187287478125*^9, 3.5018729733447266`*^9}, 3.501873115131836*^9}] }, Open ]], Cell[CellGroupData[{ Cell["\<\ \:041b\:0430\:0442\:0438\:043d\:0441\:043a\:0438\:0439 \:043a\:0432\:0430\ \:0434\:0440\:0430\:0442 \ \>", "Subsubtitle", CellChangeTimes->{{3.5018748476123047`*^9, 3.501874876573242*^9}, { 3.501876390213867*^9, 3.501876391435547*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"ksh", "=", "3"}], ";"}], " ", RowBox[{"(*", RowBox[{ "\:043a\:043e\:043b\:0438\:0447\:0435\:0441\:0442\:0432\:043e", " ", "\:0443\:0440\:043e\:0432\:043d\:0435\:0439", " ", "\:0434\:043b\:044f", " ", "\:043a\:0430\:0436\:0434\:043e\:0433\:043e", " ", "\:0444\:0430\:043a\:0442\:043e\:0440\:0430"}], "*)"}]}]], "Input", CellChangeTimes->{{3.5018749162734375`*^9, 3.5018749166660156`*^9}, { 3.501874987245117*^9, 3.501875035779297*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"step", "=", RowBox[{ RowBox[{"(", RowBox[{"1.9", "-", "0.1"}], ")"}], "/", "2"}]}], " ", ";"}], RowBox[{"(*", RowBox[{ "\:0432\:0435\:043b\:0438\:0447\:0438\:043d\:0430", " ", "\:0448\:0430\:0433\:0430"}], "*)"}]}]], "Input", CellChangeTimes->{{3.5018750422353516`*^9, 3.5018750425820312`*^9}, { 3.501875974326172*^9, 3.5018759888515625`*^9}, 3.5018760821328125`*^9}], Cell[BoxData[ StyleBox[ RowBox[{"(*", RowBox[{ "\:0437\:043d\:0430\:0447\:0435\:043d\:0438\:0435", " ", "\:0432\:0441\:0435\:0445", " ", "\:0444\:0430\:043a\:0442\:043e\:0440\:043e\:0432", " ", "\:043d\:0430", " ", "\:043a\:0430\:0436\:0434\:043e\:043c", " ", RowBox[{"\:0443\:0440\:043e\:0432\:043d\:0435", ":"}]}], "*)"}], FontWeight->"Bold"]], "Input", CellChangeTimes->{{3.5018763574052734`*^9, 3.5018763665722656`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"myuf", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"myun", "*", RowBox[{"(", RowBox[{"0.1", "+", RowBox[{"step", "*", RowBox[{"(", "i", ")"}]}]}], ")"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "0", ",", "ksh"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"nyuf", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"nyun", "*", RowBox[{"(", RowBox[{"0.1", "+", RowBox[{"step", "*", RowBox[{"(", "i", ")"}]}]}], ")"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "0", ",", "ksh"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"iif", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"in", "*", RowBox[{"(", RowBox[{"0.1", "+", RowBox[{"step", "*", RowBox[{"(", "i", ")"}]}]}], ")"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "0", ",", "ksh"}], "}"}]}], "]"}]}], ";"}]}], "Input", CellChangeTimes->{{3.5018761392460938`*^9, 3.5018761487382812`*^9}, { 3.501876187915039*^9, 3.5018763474189453`*^9}}], Cell[BoxData[ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "<", RowBox[{"ksh", "+", "1"}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", "1"}], ",", RowBox[{"j", "<", RowBox[{"ksh", "+", "1"}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"myu", "=", RowBox[{"myuf", "[", RowBox[{"[", "j", "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"nyu", "=", RowBox[{"nyuf", "[", RowBox[{"[", "j", "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"i", "=", RowBox[{"iif", "[", RowBox[{"[", "j", "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"M1", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"(", RowBox[{"s", "+", "b1"}], ")"}], ",", "0", ",", "b3"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", "s", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"myu", ",", RowBox[{"(", RowBox[{"ii", "+", RowBox[{"nyu", "/", "s"}]}], ")"}], ",", RowBox[{"-", "1"}]}], "}"}]}], "}"}]}], ";", RowBox[{"(*", RowBox[{ "\:041c\:0430\:0442\:0440\:0438\:0446\:0430", " ", "\:0441\:043e\:0441\:0442\:043e\:044f\:043d\:0438\:044f", " ", "\:0441\:0432\:043e\:0431\:043e\:0434\:043d\:043e\:0433\:043e", " ", "\:0441\:0430\:043c\:043e\:043b\:0435\:0442\:0430"}], "*)"}], "\[IndentingNewLine]", RowBox[{"M0", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"(", RowBox[{"s", "+", "b1"}], ")"}], ",", "0", ",", "b3"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"myu", ",", RowBox[{"nyu", "/", "s"}], ",", RowBox[{"-", "1"}]}], "}"}]}], "}"}]}], ";", RowBox[{"(*", RowBox[{ "\:041c\:0430\:0442\:0440\:0438\:0446\:0430", " ", "\:0441\:043e\:0441\:0442\:043e\:044f\:043d\:0438\:044f", " ", "\:0432\:043e\:0437\:043c\:0443\:0449\:0435\:043d\:043d\:043e\ \:0433\:043e", " ", "\:0441\:0430\:043c\:043e\:043b\:0435\:0442\:0430"}], "*)"}], "\[IndentingNewLine]", FormBox[ RowBox[{"wn", "=", RowBox[{ RowBox[{"Det", "[", "M0", "]"}], "/", RowBox[{"Det", "[", "M1", "]"}]}]}], TraditionalForm], ";", RowBox[{"(*", RowBox[{ "\:041f\:0435\:0440\:0435\:0434\:0430\:0442\:043e\:0447\:043d\ \:0430\:044f", " ", "\:0444\:0443\:043d\:043a\:0446\:0438\:044f"}], "*)"}], "\[IndentingNewLine]", RowBox[{"hkrap", "=", RowBox[{"InverseLaplaceTransform", "[", RowBox[{ RowBox[{"wn", "/", "s"}], ",", "s", ",", "t"}], "]"}]}], ";", RowBox[{"(*", RowBox[{ "\:041f\:0435\:0440\:0435\:0445\:043e\:0434\:043d\:0430\:044f", " ", "\:0444\:0443\:043d\:043a\:0446\:0438\:044f"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"ht", "[", "t_", "]"}], "=", RowBox[{"Abs", "[", "hkrap", "]"}]}], ";", "\[IndentingNewLine]", FormBox[ RowBox[{ RowBox[{ RowBox[{"us", "=", RowBox[{"ht", "[", "10000", "]"}]}], ";"}], RowBox[{"(*", RowBox[{ "\:0437\:0430\:0434\:0430\:0435\:043c", " ", "\:0443\:0441\:0442\:0430\:043d\:043e\:0432\:0438\:0432\:0448\ \:0435\:0435\:0441\:044f", " ", "\:0437\:043d\:0430\:0447\:0435\:043d\:0438\:0435"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"hmax", "=", "1"}], ";"}], RowBox[{"(*", RowBox[{ "\:0437\:0430\:0434\:0430\:0435\:043c", " ", "\:043c\:0430\:043a\:0441\:0438\:043c\:0430\:043b\:044c\:043d\ \:043e\:0435", " ", "\:0437\:043d\:0430\:0447\:0435\:043d\:0438\:0435"}], "*)"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{ "\:0427\:0438\:0441\:043b\:0435\:043d\:043d\:043e\:0435", " ", "\:0440\:0435\:0448\:0435\:043d\:0438\:0435", " ", "\:0443\:0440\:0430\:0432\:043d\:0435\:043d\:0438\:0439", " ", RowBox[{"h", "(", "t", ")"}]}], "=", RowBox[{ RowBox[{"0.95", " ", "\:0438", " ", RowBox[{"h", "(", "t", ")"}]}], "=", "1.05"}]}], ",", " ", RowBox[{ RowBox[{ RowBox[{ "\:043e\:0431\:044a\:0435\:0434\:0438\:043d\:0435\:043d\:043d\ \:044b\:0445", " ", "\:0432", " ", "\:043e\:0434\:043d\:043e"}], " ", "|", RowBox[{ RowBox[{"h", "(", "t", ")"}], "-", "1"}], "|", RowBox[{"-", "0.05"}]}], "=", "0"}], ",", " ", RowBox[{ RowBox[{ "\:0442", ".", "\:0435", ".", " ", "\:043d\:0430\:0445\:043e\:0436\:0434\:0435\:043d\:0438\:0435\ "}], " ", "\:0432\:0440\:0435\:043c\:0435\:043d\:0438", " ", "\:0440\:0435\:0433\:0443\:043b\:0438\:0440\:043e\:0432\:0430\ \:043d\:0438\:044f"}], ",", " ", RowBox[{ "\:0438", " ", "\:0442\:0430\:043a\:0436\:0435", " ", "\:0440\:0430\:0441\:0447\:0435\:0442", " ", "\:043f\:0435\:0440\:0435\:0440\:0435\:0433\:0443\:043b\:0438\ \:0440\:043e\:0432\:0430\:043d\:0438\:044f"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"l", "=", "1"}], ",", RowBox[{"l", "<", RowBox[{"ksh", "+", "1"}]}], ",", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Sign", "[", RowBox[{ RowBox[{"Abs", "[", RowBox[{ RowBox[{"ht", "[", RowBox[{"dt", "*", RowBox[{"(", RowBox[{"l", "-", "1"}], ")"}]}], "]"}], "-", "us"}], "]"}], "-", "0.05"}], "]"}], "\[NotEqual]", RowBox[{"Sign", "[", RowBox[{ RowBox[{"Abs", "[", RowBox[{ RowBox[{"ht", "[", RowBox[{"dt", "*", RowBox[{"(", "l", ")"}]}], "]"}], "-", "us"}], "]"}], "-", "0.05"}], "]"}]}], ",", RowBox[{"trs", "=", RowBox[{"dt", "*", RowBox[{"(", RowBox[{"l", "-", "0.5"}], ")"}]}]}]}], "]"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{"hmax", "<", RowBox[{"ht", "[", RowBox[{"dt", "*", "i"}], "]"}]}], ",", RowBox[{"hmax", "=", RowBox[{"ht", "[", RowBox[{"dt", "*", "i"}], "]"}]}]}], "]"}]}], "}"}], ";", RowBox[{"l", "++"}]}]}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"trs", "<", "minTreg"}], ",", RowBox[{"minTreg", "=", "trs"}], ",", RowBox[{"sigm", "=", RowBox[{"Round", "[", RowBox[{"Abs", "[", RowBox[{"100", "*", RowBox[{ RowBox[{"(", RowBox[{"hmax", "-", "us"}], ")"}], "/", "us"}]}], "]"}], "]"}]}]}], "]"}]}], TraditionalForm]}], RowBox[{"(*", RowBox[{ "\:041d\:0430\:0445\:043e\:0434\:0438\:043c", " ", "\:043c\:0438\:043d\:0438\:043c\:0430\:043b\:044c\:043d\:043e\ \:0435", " ", "\:0432\:0440\:0435\:043c\:044f", " ", "\:0440\:0435\:0433\:0443\:043b\:0438\:0440\:043e\:0432\:0430\ \:043d\:0438\:044f", " ", "\:0438", " ", "\:043f\:0435\:0440\:0435\:0440\:0435\:0433\:0443\:043b\:0438\ \:0440\:043e\:0432\:0430\:043d\:0438\:0435"}], "*)"}], "}"}], ";", "\[IndentingNewLine]", RowBox[{"j", "++"}]}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"iif", "=", RowBox[{"RotateLeft", "[", RowBox[{"iif", ",", "1"}], "]"}]}], ";"}], RowBox[{"(*", RowBox[{ "\:0421\:0434\:0432\:0438\:0433", " ", "\:043c\:0430\:0441\:0441\:0438\:0432\:0430", " ", "\:043f\:0435\:0440\:0435\:0434\:0430\:0442\:043e\:0447\:043d\:043e\ \:0433\:043e", " ", "\:0447\:0438\:0441\:043b\:0430", " ", "i", " ", "\:0432\:043b\:0435\:0432\:043e", " ", "\:043d\:0430", " ", "1", " ", "\:043f\:043e\:0437\:0438\:0446\:0438\:044e"}], "*)"}], "}"}], ";", "\[IndentingNewLine]", RowBox[{"i", "++"}]}]}], "]"}]], "Input", CellChangeTimes->{{3.5018763793291016`*^9, 3.5018765761816406`*^9}, { 3.501876606413086*^9, 3.501876661216797*^9}, {3.501876704364258*^9, 3.501876866826172*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"myuoptim", "=", RowBox[{ RowBox[{"(", RowBox[{"18", "-", RowBox[{"b1", "*", "tr"}]}], ")"}], "/", RowBox[{"(", RowBox[{"b3", "*", "minTreg"}], ")"}]}]}], ";"}], RowBox[{"(*", RowBox[{ "\:043e\:043f\:0442\:0438\:043c\:0430\:043b\:044c\:043d\:044b\:0435", " ", "\:0437\:043d\:0430\:0447\:0435\:043d\:0438\:044f", " ", "\:043f\:0435\:0440\:0435\:0434\:0430\:0442\:043e\:0447\:043d\:044b\:0445\ ", " ", "\:0447\:0438\:0441\:0435\:043b"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"nyuoptim", "=", RowBox[{"216", "/", RowBox[{"(", RowBox[{"b3", "*", RowBox[{"minTreg", "^", "3"}]}], ")"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"ioptim", "=", RowBox[{"108", "/", RowBox[{"(", RowBox[{"b3", "*", RowBox[{"minTreg", "^", "2"}]}], ")"}]}]}], ";"}]}], "Input", CellChangeTimes->{{3.5018768910908203`*^9, 3.5018768921992188`*^9}}], Cell[BoxData[ RowBox[{"Print", "[", RowBox[{ "\"\<\:0423\:0447\:0438\:0442\:044b\:044b\:0432\:0430\:044f \:0437\:043d\ \:0430\:0447\:0435\:043d\:0438\:044f \:0432\:0440\:0435\:043c\:0435\:043d\ \:0438 \:0440\:0435\:0433\:0443\:043b\:0438\:0440\:043e\:0432\:0430\:043d\ \:0438\:044f \:0438 \:043f\:0435\:0440\:0435\:0440\:0435\:0433\:0443\:043b\ \:0438\:0440\:043e\:0432\:0430\:043d\:0438\:044f, \:043f\:043e\:043b\:0443\ \:0447\:0430\:0435\:043c, \:0447\:0442\:043e \:043d\:0430\:0438\:0431\:043e\ \:043b\:0435\:0435 \:043e\:043f\:0442\:0438\:043c\:0430\:043b\:044c\:043d\ \:044b\:0435 \:0437\:043d\:0430\:0447\:0435\:043d\:0438\:044f \:043f\:0435\ \:0440\:0435\:0434\:0430\:0442\:043e\:0447\:043d\:044b\:0445 \:0447\:0438\ \:0441\:0435\:043b \:0430\:0432\:0442\:043e\:043f\:0438\:043b\:043e\:0442\ \:0430 \:0440\:0430\:0432\:043d\:044b: myu=\>\"", ",", " ", "myuoptim", ",", "\"\<, nyu=\>\"", ",", "nyuoptim", ",", "\"\<, i=\>\"", ",", "ioptim", ",", "\"\<. \:041f\:0440\:0438 \:044d\:0442\:043e\:043c \:0432\:0440\:0435\ \:043c\:044f \:0440\:0435\:0433\:0443\:043b\:0438\:0440\:043e\:0432\:0430\ \:043d\:0438\:044f \:0441\:043e\:0441\:0442\:0430\:0432\:043b\:044f\:0435\ \:0442 \>\"", ",", " ", "minTreg", ",", "\"\\"", ",", "sigm", ",", "\"\<%.\>\""}], "]"}]], "Input", CellChangeTimes->{{3.5018769065859375`*^9, 3.501876907665039*^9}}] }, Open ]], Cell["\<\ \:0412 \:0440\:0435\:0437\:0443\:043b\:044c\:0442\:0430\:0442\:0435 \:043f\ \:043e\:043b\:0443\:0447\:0438\:043b\:0438, \:0447\:0442\:043e \:043e\:043f\ \:0442\:0438\:043c\:0430\:043b\:044c\:043d\:044b\:0435 \:0437\:043d\:0430\ \:0447\:0435\:043d\:0438\:044f \:043f\:0435\:0440\:0435\:0434\:0430\:0442\ \:043e\:0447\:043d\:044b\:0445 \:0447\:0438\:0441\:0435\:043b \:0430\:0432\ \:0442\:043e\:043f\:0438\:043b\:043e\:0442\:0430 \:043f\:0440\:0438 \:0438\ \:0441\:043f\:043e\:043b\:044c\:0437\:043e\:0432\:0430\:043d\:0438\:0438 \ \:043f\:043e\:043b\:043d\:043e\:0433\:043e \:0444\:0430\:043a\:0442\:043e\ \:0440\:043d\:043e\:0433\:043e \:044d\:043a\:0441\:043f\:0435\:0440\:0438\ \:043c\:0435\:043d\:0442\:0430 \:0438 \:043b\:0430\:0442\:0438\:043d\:0441\ \:043a\:043e\:0433\:043e \:043a\:0432\:0430\:0434\:0440\:0430\:0442\:0430 \ \:0441\:043e\:0432\:043f\:0430\:0434\:0430\:044e\:0442. \:041f\:0440\:0438 \ \:044d\:0442\:043e\:043c \:0432 \:043f\:0435\:0440\:0432\:043e\:043c \:043c\ \:0435\:0442\:043e\:0434\:0435 \:0431\:044b\:043b\:043e N^3 \ \:043f\:0440\:043e\:0433\:043e\:043d\:043e\:0432 (\:0447\:0438\:0441\:043b\ \:043e \:0443\:0440\:043e\:0432\:043d\:0435\:0439 N = 3), \:0430 \:0432\:043e \ \:0432\:0442\:043e\:0440\:043e\:043c - N^2. \:0421\:043b\:0435\:0434\:043e\ \:0432\:0430\:0442\:0435\:043b\:044c\:043d\:043e, \:0438\:0441\:043f\:043e\ \:043b\:044c\:0437\:0443\:044f \:043c\:0435\:0442\:043e\:0434 \"\:041b\:0430\ \:0442\:0438\:043d\:0441\:043a\:0438\:0439 \:043a\:0432\:0430\:0434\:0440\ \:0430\:0442\" \:043c\:044b \:044d\:043a\:043e\:043d\:043e\:043c\:0438\:043c \ \:0432\:0440\:0435\:043c\:044f \:0438 \:0440\:0435\:0441\:0443\:0440\:0441\ \:044b.\ \>", "Subsubtitle", CellChangeTimes->{{3.5018770627109375`*^9, 3.501877180864258*^9}, { 3.5018772217333984`*^9, 3.5018772265742188`*^9}, {3.5019039257177734`*^9, 3.5019039310273438`*^9}, {3.5019040306289062`*^9, 3.5019040473115234`*^9}, { 3.5019040830195312`*^9, 3.501904206899414*^9}}] }, Open ]] }, WindowSize->{1350, 633}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, FrontEndVersion->"7.0 for Microsoft Windows (32-bit) (November 10, 2008)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[567, 22, 310, 5, 83, "Title"], Cell[880, 29, 247, 7, 30, "Input"], Cell[1130, 38, 240, 7, 30, "Input"], Cell[1373, 47, 323, 10, 30, "Input"], Cell[1699, 59, 286, 9, 31, "Input"], Cell[1988, 70, 286, 9, 31, "Input"], Cell[2277, 81, 126, 3, 31, "Input"], Cell[2406, 86, 663, 19, 31, "Input"], Cell[3072, 107, 812, 23, 50, "Input"], Cell[3887, 132, 143, 4, 30, "Input"], Cell[4033, 138, 370, 11, 30, "Input"], Cell[4406, 151, 105, 3, 30, "Input"], Cell[4514, 156, 332, 10, 50, "Input"], Cell[4849, 168, 332, 10, 50, "Input"], Cell[5184, 180, 288, 9, 30, "Input"], Cell[5475, 191, 205, 6, 50, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[5717, 202, 300, 7, 191, "Title"], Cell[6020, 211, 243, 7, 30, "Input"], Cell[6266, 220, 243, 7, 30, "Input"], Cell[6512, 229, 461, 7, 29, "Text"], Cell[6976, 238, 687, 20, 31, "Input"], Cell[7666, 260, 703, 22, 31, "Input"], Cell[8372, 284, 647, 20, 31, "Input"], Cell[9022, 306, 165, 4, 31, "Input"], Cell[9190, 312, 461, 7, 29, "Text"], Cell[9654, 321, 687, 21, 31, "Input"], Cell[10344, 344, 651, 20, 31, "Input"], Cell[10998, 366, 153, 4, 31, "Input"], Cell[11154, 372, 930, 25, 31, "Input"], Cell[12087, 399, 363, 9, 31, "Input"], Cell[12453, 410, 765, 20, 31, "Input"], Cell[13221, 432, 368, 11, 31, "Input"], Cell[13592, 445, 1004, 26, 31, "Input"], Cell[14599, 473, 523, 12, 31, "Input"], Cell[15125, 487, 595, 16, 31, "Input"], Cell[15723, 505, 387, 12, 31, "Input"], Cell[16113, 519, 420, 13, 31, "Input"], Cell[16536, 534, 187, 6, 31, "Input"], Cell[16726, 542, 109, 3, 31, "Input"], Cell[16838, 547, 802, 23, 30, "Input"], Cell[17643, 572, 334, 11, 50, "Input"], Cell[17980, 585, 216, 7, 30, "Input"], Cell[18199, 594, 258, 8, 50, "Input"], Cell[18460, 604, 258, 8, 50, "Input"], Cell[18721, 614, 205, 6, 50, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[18963, 625, 300, 7, 191, "Title"], Cell[19266, 634, 290, 8, 50, "Input"], Cell[19559, 644, 250, 7, 52, "Input"], Cell[19812, 653, 229, 7, 50, "Input"], Cell[CellGroupData[{ Cell[20066, 664, 500, 8, 28, "Subsubtitle"], Cell[20569, 674, 422, 12, 31, "Input"], Cell[20994, 688, 383, 11, 52, "Input"], Cell[21380, 701, 317, 9, 52, "Input"], Cell[21700, 712, 201, 4, 52, "Input"], Cell[21904, 718, 327, 8, 31, "Input"], Cell[22234, 728, 347, 9, 52, "Input"], Cell[22584, 739, 342, 9, 52, "Input"], Cell[22929, 750, 170, 4, 52, "Input"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[23148, 760, 356, 8, 191, "Title"], Cell[23507, 770, 805, 24, 52, "Input"], Cell[24315, 796, 338, 7, 47, "Text"], Cell[24656, 805, 271, 7, 52, "Input"], Cell[24930, 814, 155, 3, 52, "Input"], Cell[25088, 819, 169, 3, 52, "Input"], Cell[CellGroupData[{ Cell[25282, 826, 454, 8, 47, "Subsubtitle"], Cell[25739, 836, 216, 5, 52, "Input"], Cell[25958, 843, 1153, 29, 52, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[27148, 877, 508, 9, 47, "Subsubtitle"], Cell[27659, 888, 340, 10, 52, "Input"], Cell[28002, 900, 295, 8, 52, "Input"], Cell[28300, 910, 296, 8, 52, "Input"], Cell[28599, 920, 294, 8, 52, "Input"], Cell[28896, 930, 177, 5, 52, "Input"], Cell[29076, 937, 394, 9, 72, "Input"], Cell[29473, 948, 402, 12, 31, "Input"], Cell[29878, 962, 208, 6, 52, "Input"], Cell[30089, 970, 242, 6, 31, "Input"], Cell[30334, 978, 429, 13, 31, "Input"], Cell[30766, 993, 563, 16, 52, "Input"], Cell[31332, 1011, 253, 7, 31, "Input"], Cell[31588, 1020, 541, 16, 31, "Input"], Cell[32132, 1038, 2347, 54, 92, "Input"], Cell[34482, 1094, 26, 0, 31, "Input"], Cell[34511, 1096, 122, 3, 31, "Input"], Cell[34636, 1101, 124, 3, 31, "Input"], Cell[34763, 1106, 124, 3, 31, "Input"], Cell[34890, 1111, 489, 16, 31, "Input"], Cell[35382, 1129, 489, 16, 31, "Input"], Cell[35874, 1147, 213, 4, 29, "Text"], Cell[36090, 1153, 407, 12, 31, "Input"], Cell[36500, 1167, 195, 4, 29, "Text"], Cell[36698, 1173, 408, 12, 31, "Input"], Cell[37109, 1187, 488, 16, 31, "Input"], Cell[37600, 1205, 489, 16, 31, "Input"], Cell[38092, 1223, 208, 4, 29, "Text"], Cell[38303, 1229, 407, 12, 31, "Input"], Cell[38713, 1243, 190, 4, 29, "Text"], Cell[38906, 1249, 407, 12, 31, "Input"], Cell[39316, 1263, 631, 16, 31, "Input"], Cell[39950, 1281, 362, 11, 31, "Input"], Cell[40315, 1294, 297, 8, 31, "Input"], Cell[40615, 1304, 315, 9, 31, "Input"], Cell[40933, 1315, 252, 7, 31, "Input"], Cell[41188, 1324, 1166, 31, 52, "Input"], Cell[42357, 1357, 1062, 28, 31, "Input"], Cell[43422, 1387, 104, 3, 31, "Input"], Cell[43529, 1392, 130, 3, 31, "Input"], Cell[43662, 1397, 203, 6, 31, "Input"], Cell[43868, 1405, 141, 4, 31, "Input"], Cell[44012, 1411, 1112, 29, 31, "Input"], Cell[45127, 1442, 905, 28, 52, "Input"], Cell[46035, 1472, 750, 22, 31, "Input"], Cell[46788, 1496, 590, 12, 31, "Input"], Cell[47381, 1510, 108, 2, 31, "Input"], Cell[47492, 1514, 226, 6, 31, "Input"], Cell[47721, 1522, 315, 7, 31, "Input"], Cell[48039, 1531, 204, 7, 31, "Input"], Cell[48246, 1540, 137, 4, 31, "Input"], Cell[48386, 1546, 150, 4, 31, "Input"], Cell[48539, 1552, 931, 25, 31, "Input"], Cell[49473, 1579, 881, 24, 31, "Input"], Cell[50357, 1605, 460, 11, 31, "Input"], Cell[50820, 1618, 517, 12, 31, "Input"], Cell[51340, 1632, 428, 12, 31, "Input"], Cell[51771, 1646, 486, 14, 31, "Input"], Cell[52260, 1662, 113, 3, 31, "Input"], Cell[52376, 1667, 529, 15, 52, "Input"], Cell[52908, 1684, 918, 20, 52, "Input"], Cell[53829, 1706, 1368, 33, 52, "Input"], Cell[55200, 1741, 1548, 48, 72, "Input"], Cell[56751, 1791, 68, 2, 31, "Input"], Cell[56822, 1795, 522, 12, 31, "Input"], Cell[57347, 1809, 520, 12, 31, "Input"], Cell[57870, 1823, 145, 3, 31, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[58052, 1831, 251, 6, 66, "Subsubtitle"], Cell[58306, 1839, 13287, 283, 412, "Input"], Cell[71596, 2124, 1055, 30, 72, "Input"], Cell[72654, 2156, 1706, 26, 52, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[74397, 2187, 249, 7, 66, "Subsubtitle"], Cell[74649, 2196, 509, 11, 31, "Input"], Cell[75161, 2209, 450, 12, 31, "Input"], Cell[75614, 2223, 457, 10, 31, "Input"], Cell[76074, 2235, 1151, 37, 72, "Input"], Cell[77228, 2274, 10059, 229, 372, "Input"], Cell[87290, 2505, 1004, 29, 72, "Input"], Cell[88297, 2536, 1483, 23, 52, "Input"] }, Open ]], Cell[89795, 2562, 1979, 30, 85, "Subsubtitle"] }, Open ]] } ] *) (* End of internal cache information *)