(************** Content-type: application/mathematica ************** CreatedBy='Mathematica 5.2' Mathematica-Compatible Notebook This notebook can be used with any Mathematica-compatible application, such as Mathematica, MathReader or Publicon. The data for the notebook starts with the line containing stars above. To get the notebook into a Mathematica-compatible application, do one of the following: * Save the data starting with the line of stars above into a file with a name ending in .nb, then open the file inside the application; * Copy the data starting with the line of stars above to the clipboard, then use the Paste menu command inside the application. Data for notebooks contains only printable 7-bit ASCII and can be sent directly in email or through ftp in text mode. Newlines can be CR, LF or CRLF (Unix, Macintosh or MS-DOS style). NOTE: If you modify the data for this notebook not in a Mathematica- compatible application, you must delete the line below containing the word CacheID, otherwise Mathematica-compatible applications may try to use invalid cache data. For more information on notebooks and Mathematica-compatible applications, contact Wolfram Research: web: http://www.wolfram.com email: info@wolfram.com phone: +1-217-398-0700 (U.S.) Notebook reader applications are available free of charge from Wolfram Research. *******************************************************************) (*CacheID: 232*) (*NotebookFileLineBreakTest NotebookFileLineBreakTest*) (*NotebookOptionsPosition[ 83758, 2447]*) (*NotebookOutlinePosition[ 84439, 2470]*) (* CellTagsIndexPosition[ 84395, 2466]*) (*WindowFrame->Normal*) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[ FormBox[ RowBox[{"\[IndentingNewLine]", StyleBox[\(\: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\), "Subtitle"]}], TraditionalForm]], "Input"], Cell[BoxData[ \(\:0434\:0432\:0438\:0436\:0435\:043d\:0438\:044f\ \:0443\:0440\:0430\ \:0432\:043d\:0435\:043d\:0438\:044f\ \:041a\:0438\:043d\:0435\:043c\:0430\ \:0442\:0438\:0447\:0435\:0441\:043a\:0438\:0435\ \:043f\:043e\:0441\:0442\ \:0443\:043f\:0430\:0442\:0435\:043b\:044c\:043d\:043e\:0433\:043e\)], \ "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(TraditionalForm\`vg = {vgx[t], \ vgy[t], \ vgz[t]}\)], "Input"], Cell[BoxData[ \({vgx[t], vgy[t], vgz[t]}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(TraditionalForm\`r = {l[t], \ h[t], \ z[t]}\)], "Input"], Cell[BoxData[ \({l[t], h[t], z[t]}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(TraditionalForm\`kur1 = TableForm[D[r[\([1]\)], t] == vg[\([1]\)]]\)], "Input"], Cell[BoxData[ TagBox[ RowBox[{ RowBox[{ SuperscriptBox["l", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", \(vgx[t]\)}], Function[ BoxForm`e$, TableForm[ BoxForm`e$]]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(kur2 = TableForm[D[r[\([2]\)], t] == vg[\([2]\)]]\)], "Input"], Cell[BoxData[ TagBox[ RowBox[{ RowBox[{ SuperscriptBox["h", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", \(vgy[t]\)}], Function[ BoxForm`e$, TableForm[ BoxForm`e$]]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(kur3 = TableForm[D[r[\([3]\)], t] == vg[\([3]\)]]\)], "Input"], Cell[BoxData[ TagBox[ RowBox[{ RowBox[{ SuperscriptBox["z", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", \(vgz[t]\)}], Function[ BoxForm`e$, TableForm[ BoxForm`e$]]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(pr = {kur1, \ kur2, \ kur3}\)], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{ TagBox[ RowBox[{ RowBox[{ SuperscriptBox["l", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", \(vgx[t]\)}], Function[ BoxForm`e$, TableForm[ BoxForm`e$]]], ",", TagBox[ RowBox[{ RowBox[{ SuperscriptBox["h", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", \(vgy[t]\)}], Function[ BoxForm`e$, TableForm[ BoxForm`e$]]], ",", TagBox[ RowBox[{ RowBox[{ SuperscriptBox["z", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", \(vgz[t]\)}], Function[ BoxForm`e$, TableForm[ BoxForm`e$]]]}], "}"}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(TableForm[ m1 = {{Cos[\[Theta][t]], \(-Sin[\[Theta][t]]\), 0}, {Sin[\[Theta][t]], Cos[\[Theta][t]], 0}, {0, 0, 1}}]\)], "Input"], Cell[BoxData[ TagBox[GridBox[{ {\(Cos[\[Theta][t]]\), \(-Sin[\[Theta][t]]\), "0"}, {\(Sin[\[Theta][t]]\), \(Cos[\[Theta][t]]\), "0"}, {"0", "0", "1"} }, RowSpacings->1, ColumnSpacings->3, RowAlignments->Baseline, ColumnAlignments->{Left}], Function[ BoxForm`e$, TableForm[ BoxForm`e$]]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(TraditionalForm\`\(\(\[IndentingNewLine]\)\(TableForm[ m2 = {{Cos[\[Psi][t]], \ 0, \ Sin[\[Psi][t]]}, \ {0, \ 1, \ 0}, \ {\(-Sin[\[Psi][t]]\), \ 0, \ Cos[\[Psi][t]]}}]\)\)\)], "Input"], Cell[BoxData[ TagBox[GridBox[{ {\(Cos[\[Psi][t]]\), "0", \(Sin[\[Psi][t]]\)}, {"0", "1", "0"}, {\(-Sin[\[Psi][t]]\), "0", \(Cos[\[Psi][t]]\)} }, RowSpacings->1, ColumnSpacings->3, RowAlignments->Baseline, ColumnAlignments->{Left}], Function[ BoxForm`e$, TableForm[ BoxForm`e$]]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(TraditionalForm\`TableForm[m3 = m2 . m1]\)], "Input"], Cell[BoxData[ TagBox[GridBox[{ {\(Cos[\[Theta][t]]\ Cos[\[Psi][ t]]\), \(\(-Cos[\[Psi][t]]\)\ Sin[\[Theta][ t]]\), \(Sin[\[Psi][t]]\)}, {\(Sin[\[Theta][t]]\), \(Cos[\[Theta][t]]\), "0"}, {\(\(-Cos[\[Theta][t]]\)\ Sin[\[Psi][t]]\), \(Sin[\[Theta][ t]]\ Sin[\[Psi][t]]\), \(Cos[\[Psi][t]]\)} }, RowSpacings->1, ColumnSpacings->3, RowAlignments->Baseline, ColumnAlignments->{Left}], Function[ BoxForm`e$, TableForm[ BoxForm`e$]]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(TraditionalForm\`TableForm[ ReplaceAll[ m3, \ {\[Theta][t] \[Rule] 0, \ \[Psi][t] \[Rule] 0}]]\)], "Input"], Cell[BoxData[ TagBox[GridBox[{ {"1", "0", "0"}, {"0", "1", "0"}, {"0", "0", "1"} }, RowSpacings->1, ColumnSpacings->3, RowAlignments->Baseline, ColumnAlignments->{Left}], Function[ BoxForm`e$, TableForm[ BoxForm`e$]]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(TraditionalForm\`rg = m3 . vg\)], "Input"], Cell[BoxData[ \({Cos[\[Theta][t]]\ Cos[\[Psi][t]]\ vgx[t] - Cos[\[Psi][t]]\ Sin[\[Theta][t]]\ vgy[t] + Sin[\[Psi][t]]\ vgz[t], Sin[\[Theta][t]]\ vgx[t] + Cos[\[Theta][t]]\ vgy[t], \(-Cos[\[Theta][t]]\)\ Sin[\[Psi][t]]\ vgx[ t] + Sin[\[Theta][t]]\ Sin[\[Psi][t]]\ vgy[t] + Cos[\[Psi][t]]\ vgz[t]}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(TraditionalForm\`\(\(\[IndentingNewLine]\)\(q1 = D[r[\([1]\)], \ t] \[Equal] rg[\([1]\)]\)\)\)], "Input"], Cell[BoxData[ RowBox[{ RowBox[{ SuperscriptBox["l", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", \(Cos[\[Theta][t]]\ Cos[\[Psi][t]]\ vgx[t] - Cos[\[Psi][t]]\ Sin[\[Theta][t]]\ vgy[t] + Sin[\[Psi][t]]\ vgz[t]\)}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(TraditionalForm\`\(\(\[IndentingNewLine]\)\(q2 = D[r[\([2]\)], \ t] \[Equal] rg[\([2]\)]\)\)\)], "Input"], Cell[BoxData[ RowBox[{ RowBox[{ SuperscriptBox["h", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", \(Sin[\[Theta][t]]\ vgx[t] + Cos[\[Theta][t]]\ vgy[t]\)}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(TraditionalForm\`q3 = D[r[\([3]\)], \ t] \[Equal] rg[\([3]\)]\)], "Input"], Cell[BoxData[ RowBox[{ RowBox[{ SuperscriptBox["z", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", \(\(-Cos[\[Theta][t]]\)\ Sin[\[Psi][t]]\ vgx[t] + Sin[\[Theta][t]]\ Sin[\[Psi][t]]\ vgy[t] + Cos[\[Psi][t]]\ vgz[t]\)}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(TraditionalForm\`\(\(\[IndentingNewLine]\)\(TableForm[{q1, \ q2, \ q3}]\)\)\)], "Input"], Cell[BoxData[ InterpretationBox[GridBox[{ { RowBox[{ RowBox[{ SuperscriptBox["l", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", \(Cos[\[Theta][t]]\ Cos[\[Psi][t]]\ vgx[t] - Cos[\[Psi][t]]\ Sin[\[Theta][t]]\ vgy[t] + Sin[\[Psi][t]]\ vgz[t]\)}]}, { RowBox[{ RowBox[{ SuperscriptBox["h", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", \(Sin[\[Theta][t]]\ vgx[t] + Cos[\[Theta][t]]\ vgy[t]\)}]}, { RowBox[{ RowBox[{ SuperscriptBox["z", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", \(\(-Cos[\[Theta][t]]\)\ Sin[\[Psi][t]]\ vgx[t] + Sin[\[Theta][t]]\ Sin[\[Psi][t]]\ vgy[t] + Cos[\[Psi][t]]\ vgz[t]\)}]} }, RowSpacings->1, ColumnSpacings->3, RowAlignments->Baseline, ColumnAlignments->{Left}], TableForm[ { Equal[ Derivative[ 1][ l][ t], Plus[ Times[ Cos[ \[Theta][ t]], Cos[ \[Psi][ t]], vgx[ t]], Times[ -1, Cos[ \[Psi][ t]], Sin[ \[Theta][ t]], vgy[ t]], Times[ Sin[ \[Psi][ t]], vgz[ t]]]], Equal[ Derivative[ 1][ h][ t], Plus[ Times[ Sin[ \[Theta][ t]], vgx[ t]], Times[ Cos[ \[Theta][ t]], vgy[ t]]]], Equal[ Derivative[ 1][ z][ t], Plus[ Times[ -1, Cos[ \[Theta][ t]], Sin[ \[Psi][ t]], vgx[ t]], Times[ Sin[ \[Theta][ t]], Sin[ \[Psi][ t]], vgy[ t]], Times[ Cos[ \[Psi][ t]], vgz[ t]]]]}]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ FormBox[ RowBox[{ "\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ StyleBox[\(\: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\), "Subtitle"], "\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]"}]}], TraditionalForm]], "Input"], Cell[BoxData[ \(\:0434\:0432\:0438\:0436\:0435\:043d\:0438\:044f\ \:0443\:0440\:0430\ \:0432\:043d\:0435\:043d\:0438\:044f\ \:0414\:0438\:043d\:0430\:043c\:0438\ \:0447\:0435\:0441\:043a\:0438\:0435\ \:043f\:043e\:0441\:0442\:0443\:043f\ \:0430\:0442\:0435\:043b\:044c\:043d\:043e\:0433\:043e\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(TraditionalForm\`v = {vx[t], \ vy[t], \ vz[t]}\)], "Input"], Cell[BoxData[ \({vx[t], vy[t], vz[t]}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(TraditionalForm\`w = {wx[t], \ wy[t], \ wz[t]}\)], "Input"], Cell[BoxData[ \({wx[t], wy[t], wz[t]}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(TraditionalForm\`f = {fx, \ fy, \ fz}\)], "Input"], Cell[BoxData[ \({fx, fy, fz}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(TraditionalForm\`\(\(\[IndentingNewLine]\)\(d = m*\((D[v, \ t] + Cross[w, \ v])\)\)\)\)], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"m", " ", RowBox[{"(", RowBox[{\(vz[t]\ wy[t]\), "-", \(vy[t]\ wz[t]\), "+", RowBox[{ SuperscriptBox["vx", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], ")"}]}], ",", RowBox[{"m", " ", RowBox[{"(", RowBox[{\(\(-vz[t]\)\ wx[t]\), "+", \(vx[t]\ wz[t]\), "+", RowBox[{ SuperscriptBox["vy", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], ")"}]}], ",", RowBox[{"m", " ", RowBox[{"(", RowBox[{\(vy[t]\ wx[t]\), "-", \(vx[t]\ wy[t]\), "+", RowBox[{ SuperscriptBox["vz", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], ")"}]}]}], "}"}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(TraditionalForm\`q4 = d[\([1]\)] == f[\([1]\)]\)], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"m", " ", RowBox[{"(", RowBox[{\(vz[t]\ wy[t]\), "-", \(vy[t]\ wz[t]\), "+", RowBox[{ SuperscriptBox["vx", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], ")"}]}], "\[Equal]", "fx"}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(TraditionalForm\`\(\(\[IndentingNewLine]\)\(q5 = d[\([2]\)] == f[\([2]\)]\)\)\)], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"m", " ", RowBox[{"(", RowBox[{\(\(-vz[t]\)\ wx[t]\), "+", \(vx[t]\ wz[t]\), "+", RowBox[{ SuperscriptBox["vy", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], ")"}]}], "\[Equal]", "fy"}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(TraditionalForm\`\(\(\[IndentingNewLine]\)\(q6 = d[\([3]\)] == f[\([3]\)]\)\)\)], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"m", " ", RowBox[{"(", RowBox[{\(vy[t]\ wx[t]\), "-", \(vx[t]\ wy[t]\), "+", RowBox[{ SuperscriptBox["vz", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], ")"}]}], "\[Equal]", "fz"}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(TraditionalForm\`\(\(\[IndentingNewLine]\)\(TableForm[{q4, \ q5, \ q6}]\)\)\)], "Input"], Cell[BoxData[ InterpretationBox[GridBox[{ { RowBox[{ RowBox[{"m", " ", RowBox[{"(", RowBox[{\(vz[t]\ wy[t]\), "-", \(vy[t]\ wz[t]\), "+", RowBox[{ SuperscriptBox["vx", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], ")"}]}], "\[Equal]", "fx"}]}, { RowBox[{ RowBox[{"m", " ", RowBox[{"(", RowBox[{\(\(-vz[t]\)\ wx[t]\), "+", \(vx[t]\ wz[t]\), "+", RowBox[{ SuperscriptBox["vy", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], ")"}]}], "\[Equal]", "fy"}]}, { RowBox[{ RowBox[{"m", " ", RowBox[{"(", RowBox[{\(vy[t]\ wx[t]\), "-", \(vx[t]\ wy[t]\), "+", RowBox[{ SuperscriptBox["vz", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], ")"}]}], "\[Equal]", "fz"}]} }, RowSpacings->1, ColumnSpacings->3, RowAlignments->Baseline, ColumnAlignments->{Left}], TableForm[ { Equal[ Times[ m, Plus[ Times[ vz[ t], wy[ t]], Times[ -1, vy[ t], wz[ t]], Derivative[ 1][ vx][ t]]], fx], Equal[ Times[ m, Plus[ Times[ -1, vz[ t], wx[ t]], Times[ vx[ t], wz[ t]], Derivative[ 1][ vy][ t]]], fy], Equal[ Times[ m, Plus[ Times[ vy[ t], wx[ t]], Times[ -1, vx[ t], wy[ t]], Derivative[ 1][ vz][ t]]], fz]}]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ FormBox[ RowBox[{ "\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]", StyleBox[\(\(\: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\)\(\ \[IndentingNewLine]\)\(\[IndentingNewLine]\)\(\[IndentingNewLine]\) \), "Subtitle"]}], TraditionalForm]], "Input"], Cell[BoxData[ \(\:0434\:0432\:0438\:0436\:0435\:043d\:0438\:044f\ \:0443\:0440\:0430\ \:0432\:043d\:0435\:043d\:0438\:044f\ \:0432\:0440\:0430\:0449\:0430\:0442\ \:0435\:043b\:044c\:043d\:043e\:0433\:043e\ \:041a\:0438\:043d\:0435\:043c\ \:0430\:0442\:0438\:0447\:0435\:0441\:043a\:0438\:0435\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(TraditionalForm\`\(\(\[IndentingNewLine]\)\(ug = {psi[t], \ teta[t], \ gam[t]}\)\)\)], "Input"], Cell[BoxData[ \({psi[t], teta[t], gam[t]}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(w = {wx[t], \ wy[t], \ wz[t]}\)\)\)], "Input"], Cell[BoxData[ \({wx[t], wy[t], wz[t]}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(TraditionalForm\`\(\(TableForm[ mm1 = {{Cos[psi[t]], \ 0, \ Sin[psi[t]]}, \ {0, \ 1, \ 0}, \ {\(-Sin[psi[t]]\), \ 0, \ Cos[psi[t]]}}]\)\(\[IndentingNewLine]\) \)\)], "Input"], Cell[BoxData[ TagBox[GridBox[{ {\(Cos[psi[t]]\), "0", \(Sin[psi[t]]\)}, {"0", "1", "0"}, {\(-Sin[psi[t]]\), "0", \(Cos[psi[t]]\)} }, RowSpacings->1, ColumnSpacings->3, RowAlignments->Baseline, ColumnAlignments->{Left}], Function[ BoxForm`e$, TableForm[ BoxForm`e$]]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(TraditionalForm\`\(\(\[IndentingNewLine]\)\(TableForm[ mm2 = {{Cos[teta[t]], \ \(-Sin[teta[t]]\), \ 0}, \ {Sin[teta[t]], \ Cos[teta[t]], \ 0}, \ {0, \ 0, \ 1}}]\)\)\)], "Input", CellChangeTimes -> {{3.4970717489335938`*^9, 3.4970717556064453`*^9}, \ {3.497071793182617*^9, 3.497071809453125*^9}}], Cell[BoxData[ TagBox[GridBox[{ {\(Cos[teta[t]]\), \(-Sin[teta[t]]\), "0"}, {\(Sin[teta[t]]\), \(Cos[teta[t]]\), "0"}, {"0", "0", "1"} }, RowSpacings->1, ColumnSpacings->3, RowAlignments->Baseline, ColumnAlignments->{Left}], Function[ BoxForm`e$, TableForm[ BoxForm`e$]]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(TraditionalForm\`\(\(\[IndentingNewLine]\)\(TableForm[ mm3 = {{1, 0, 0}, {0, Cos[gam[t]], \(-Sin[gam[t]]\)}, {0, Sin[gam[t]], Cos[gam[t]]}}]\)\)\)], "Input", CellChangeTimes -> {{3.497071746642578*^9, 3.497071790451172*^9}}], Cell[BoxData[ TagBox[GridBox[{ {"1", "0", "0"}, {"0", \(Cos[gam[t]]\), \(-Sin[gam[t]]\)}, {"0", \(Sin[gam[t]]\), \(Cos[gam[t]]\)} }, RowSpacings->1, ColumnSpacings->3, RowAlignments->Baseline, ColumnAlignments->{Left}], Function[ BoxForm`e$, TableForm[ BoxForm`e$]]]], "Output"] }, Open ]], Cell[BoxData[ \(TraditionalForm\`\(\(\[IndentingNewLine]\)\( (*TableForm[ mm4 = mm3 . mm2 . mm1]*) \)\)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ StyleBox[\(\: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 . \), FontSize->16]], "Input", CellChangeTimes -> {{3.4970722216083984`*^9, 3.497072265201172*^9}}], Cell[BoxData[ RowBox[{\(Syntax::"tsntxi"\), \(\(:\)\(\ \)\), \ "\<\"\\\"\\!\\(\\(\\(\: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 .\\)\\)\\)\\)\\)\\\" is incomplete; \ more input is needed.\\!\\(\\\"\\\"\\) \\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\ \\\", ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \ ButtonData:>\\\"General::tsntxi\\\"]\\)\"\>"}]], "Message"], Cell[BoxData[ StyleBox[ ErrorBox[\(\: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 . \)], ShowStringCharacters->True]], "Message"] }, Open ]], Cell[BoxData[ \(\[IndentingNewLine]\)], "Input", CellChangeTimes -> {{3.4970723490908203`*^9, 3.4970723627089844`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ \(prpsi = Transpose[mm2 . mm3] . {0, D[psi[t], t], 0}\)], "Input", CellChangeTimes -> {{3.4970723837089844`*^9, 3.4970724387285156`*^9}, \ {3.4970725976191406`*^9, 3.4970726448808594`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{\(Sin[teta[t]]\), " ", RowBox[{ SuperscriptBox["psi", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], ",", RowBox[{\(Cos[gam[t]]\), " ", \(Cos[teta[t]]\), " ", RowBox[{ SuperscriptBox["psi", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], ",", RowBox[{\(-Cos[teta[t]]\), " ", \(Sin[gam[t]]\), " ", RowBox[{ SuperscriptBox["psi", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}], "}"}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(prteta = Transpose[mm3] . {0, 0, D[teta[t], t]}\)\)\)], "Input", CellChangeTimes -> {{3.497072758939453*^9, 3.49707277871875*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"0", ",", RowBox[{\(Sin[gam[t]]\), " ", RowBox[{ SuperscriptBox["teta", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], ",", RowBox[{\(Cos[gam[t]]\), " ", RowBox[{ SuperscriptBox["teta", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}], "}"}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(prgam = {D[gam[t], t], 0, 0}\)\)\)], "Input",\ CellChangeTimes -> {{3.4970727853554688`*^9, 3.4970727986240234`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{ SuperscriptBox["gam", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], ",", "0", ",", "0"}], "}"}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(pr = prgam + prpsi + prteta\)\)\)], "Input", CellChangeTimes -> {{3.497072909881836*^9, 3.4970729462128906`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ SuperscriptBox["gam", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "+", RowBox[{\(Sin[teta[t]]\), " ", RowBox[{ SuperscriptBox["psi", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}], ",", RowBox[{ RowBox[{\(Cos[gam[t]]\), " ", \(Cos[teta[t]]\), " ", RowBox[{ SuperscriptBox["psi", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "+", RowBox[{\(Sin[gam[t]]\), " ", RowBox[{ SuperscriptBox["teta", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}], ",", RowBox[{ RowBox[{\(-Cos[teta[t]]\), " ", \(Sin[gam[t]]\), " ", RowBox[{ SuperscriptBox["psi", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "+", RowBox[{\(Cos[gam[t]]\), " ", RowBox[{ SuperscriptBox["teta", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}]}], "}"}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(q7 = w[\([1]\)] == pr[\([1]\)]\)], "Input", CellChangeTimes -> {{3.497072515258789*^9, 3.4970725713554688`*^9}, \ 3.497072659428711*^9, {3.4970729607060547`*^9, 3.4970729921777344`*^9}}], Cell[BoxData[ RowBox[{\(wx[t]\), "\[Equal]", RowBox[{ RowBox[{ SuperscriptBox["gam", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "+", RowBox[{\(Sin[teta[t]]\), " ", RowBox[{ SuperscriptBox["psi", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}]}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(q8 = w[\([2]\)] == pr[\([2]\)]\)\)\)], "Input", CellChangeTimes -> {{3.4970728045791016`*^9, 3.4970728351689453`*^9}, \ {3.4970729731845703`*^9, 3.4970729759833984`*^9}}], Cell[BoxData[ RowBox[{\(wy[t]\), "\[Equal]", RowBox[{ RowBox[{\(Cos[gam[t]]\), " ", \(Cos[teta[t]]\), " ", RowBox[{ SuperscriptBox["psi", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "+", RowBox[{\(Sin[gam[t]]\), " ", RowBox[{ SuperscriptBox["teta", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}]}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(q9 = w[\([3]\)] == pr[\([3]\)]\)\)\)], "Input", CellChangeTimes -> {{3.497072839439453*^9, 3.4970728727626953`*^9}, \ {3.49707297828125*^9, 3.4970729808779297`*^9}}], Cell[BoxData[ RowBox[{\(wz[t]\), "\[Equal]", RowBox[{ RowBox[{\(-Cos[teta[t]]\), " ", \(Sin[gam[t]]\), " ", RowBox[{ SuperscriptBox["psi", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "+", RowBox[{\(Cos[gam[t]]\), " ", RowBox[{ SuperscriptBox["teta", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}]}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(TableForm[{q7, \ q8, \ q9}]\)\)\)], "Input"], Cell[BoxData[ InterpretationBox[GridBox[{ { RowBox[{\(wx[t]\), "\[Equal]", RowBox[{ RowBox[{ SuperscriptBox["gam", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "+", RowBox[{\(Sin[teta[t]]\), " ", RowBox[{ SuperscriptBox["psi", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}]}]}, { RowBox[{\(wy[t]\), "\[Equal]", RowBox[{ RowBox[{\(Cos[gam[t]]\), " ", \(Cos[teta[t]]\), " ", RowBox[{ SuperscriptBox["psi", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "+", RowBox[{\(Sin[gam[t]]\), " ", RowBox[{ SuperscriptBox["teta", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}]}]}, { RowBox[{\(wz[t]\), "\[Equal]", RowBox[{ RowBox[{\(-Cos[teta[t]]\), " ", \(Sin[gam[t]]\), " ", RowBox[{ SuperscriptBox["psi", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "+", RowBox[{\(Cos[gam[t]]\), " ", RowBox[{ SuperscriptBox["teta", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}]}]} }, RowSpacings->1, ColumnSpacings->3, RowAlignments->Baseline, ColumnAlignments->{Left}], TableForm[ { Equal[ wx[ t], Plus[ Derivative[ 1][ gam][ t], Times[ Sin[ teta[ t]], Derivative[ 1][ psi][ t]]]], Equal[ wy[ t], Plus[ Times[ Cos[ gam[ t]], Cos[ teta[ t]], Derivative[ 1][ psi][ t]], Times[ Sin[ gam[ t]], Derivative[ 1][ teta][ t]]]], Equal[ wz[ t], Plus[ Times[ -1, Cos[ teta[ t]], Sin[ gam[ t]], Derivative[ 1][ psi][ t]], Times[ Cos[ gam[ t]], Derivative[ 1][ teta][ t]]]]}]]], "Output"] }, Open ]], Cell[BoxData[ \(\[IndentingNewLine]\)], "Input", CellChangeTimes -> {3.4970727545322266`*^9}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"\[IndentingNewLine]", StyleBox[\(\(\: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\)\(\[IndentingNewLine]\)\(\ \[IndentingNewLine]\)\(\[IndentingNewLine]\)\(\[IndentingNewLine]\)\(\ \[IndentingNewLine]\)\(\[IndentingNewLine]\) \), "Title", FontSize->24]}]], "Input", CellChangeTimes -> {{3.4970730647041016`*^9, 3.497073078645508*^9}}], Cell[BoxData[ \(\:0434\:0432\:0438\:0436\:0435\:043d\:0438\:044f\ \:0443\:0440\:0430\ \:0432\:043d\:0435\:043d\:0438\:044f\ \:0414\:0438\:043d\:0430\:043c\:0438\ \:0447\:0435\:0441\:043a\:0438\:0435\ \:0432\:0440\:0430\:0449\:0430\:0442\ \:0435\:043b\:044c\:043d\:043e\:0433\:043e\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(TableForm[ TI = {{Ixx, \(-Ixy\), \(-Ixz\)}, {\(-Iyx\), Iyy, \(-Iyz\)}, {\(-Izx\), \(-Izy\), Izz}}]\)\(\ \)\( (*\: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[BoxData[ TagBox[GridBox[{ {"Ixx", \(-Ixy\), \(-Ixz\)}, {\(-Iyx\), "Iyy", \(-Iyz\)}, {\(-Izx\), \(-Izy\), "Izz"} }, RowSpacings->1, ColumnSpacings->3, RowAlignments->Baseline, ColumnAlignments->{Left}], Function[ BoxForm`e$, TableForm[ BoxForm`e$]]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(\: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\)\)\)], "Input", CellChangeTimes -> {{3.4970732330478516`*^9, 3.4970732528828125`*^9}}], Cell[BoxData[ RowBox[{\(Syntax::"tsntxi"\), \(\(:\)\(\ \)\), \ "\<\"\\\"\\!\\(\\(\\(\\(\\(\: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\:0438XOY\\)\\), \:0442\:043e\\)\\\" is incomplete; more input is \ needed.\\!\\(\\\"\\\"\\) \\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\", \ ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \ ButtonData:>\\\"General::tsntxi\\\"]\\)\"\>"}]], "Message"], Cell[BoxData[ StyleBox[ ErrorBox[\(\: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\)], ShowStringCharacters->True]], "Message"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(Ixz = \(Izx = \(Iyz = \(\(Izy\)\(=\)\(0\)\(\ \ \)\)\)\)\)\)\)], "Input", CellChangeTimes -> {{3.4970731367734375`*^9, 3.4970731370976562`*^9}, \ 3.4970732316464844`*^9}], Cell[BoxData[ \(0\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(Iyx = Ixy\)\)\)], "Input", CellChangeTimes -> {{3.497073147826172*^9, 3.497073148583008*^9}}], Cell[BoxData[ \(Ixy\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(TableForm[TI]\)\)\)], "Input", CellChangeTimes -> {{3.4970731754121094`*^9, 3.4970731954052734`*^9}}], Cell[BoxData[ TagBox[GridBox[{ {"Ixx", \(-Ixy\), "0"}, {\(-Ixy\), "Iyy", "0"}, {"0", "0", "Izz"} }, RowSpacings->1, ColumnSpacings->3, RowAlignments->Baseline, ColumnAlignments->{Left}], Function[ BoxForm`e$, TableForm[ BoxForm`e$]]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(\: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\)\(:\)\)\)\)\)], \ "Input", CellChangeTimes -> {{3.4970732178271484`*^9, 3.4970732256503906`*^9}}], Cell[BoxData[ RowBox[{\(Syntax::"tsntxi"\), \(\(:\)\(\ \)\), \ "\<\"\\\"\\!\\(\:043a\:043e\:043e\:0440\:0434\:0438\:043d\:0430\:0442 \ :\\)\\\" is incomplete; more input is needed.\\!\\(\\\"\\\"\\) \ \\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\", ButtonStyle->\\\"RefGuideLinkText\ \\\", ButtonFrame->None, ButtonData:>\\\"General::tsntxi\\\"]\\)\"\>"}]], \ "Message"], Cell[BoxData[ StyleBox[ RowBox[{ "\: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", ErrorBox[\(\(\:043a\:043e\:043e\:0440\:0434\:0438\:043d\:0430\:0442\)\ \(:\)\)]}], ShowStringCharacters->True]], "Message"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(TableForm[K = TI . w]\)\)\)], "Input", CellChangeTimes -> {{3.497073185026367*^9, 3.4970732160976562`*^9}}], Cell[BoxData[ InterpretationBox[GridBox[{ {\(Ixx\ wx[t] - Ixy\ wy[t]\)}, {\(\(-Ixy\)\ wx[t] + Iyy\ wy[t]\)}, {\(Izz\ wz[t]\)} }, RowSpacings->1, ColumnSpacings->3, RowAlignments->Baseline, ColumnAlignments->{Left}], TableForm[ { Plus[ Times[ Ixx, wx[ t]], Times[ -1, Ixy, wy[ t]]], Plus[ Times[ -1, Ixy, wx[ t]], Times[ Iyy, wy[ t]]], Times[ Izz, wz[ t]]}]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(M = {Mx[alfa[t], mah[t]], My[alfa[t], mah[t]], Mz[alfa[t], \ mah[t]]}\)\)\)], "Input", CellChangeTimes -> {{3.4970732691376953`*^9, 3.497073269541992*^9}}], Cell[BoxData[ \({Mx[alfa[t], mah[t]], My[alfa[t], mah[t]], Mz[alfa[t], mah[t]]}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \({Mx[alfa[t], mah[t], de[t], dh[t]], My[alfa[t], mah[t]], Mz[alfa[t], mah[t]]}\)], "Input"], Cell[BoxData[ \({Mx[alfa[t], mah[t], de[t], dh[t]], My[alfa[t], mah[t]], Mz[alfa[t], mah[t]]}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(\: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\)\)\)], "Input",\ CellChangeTimes -> {{3.4970733012558594`*^9, 3.497073301758789*^9}}], Cell[BoxData[ \(\:043f\:043e\ \:041f\:043e\:043b\:043d\:0430\:044f\ \:0432\:0435\:043a\ \:0442\:043e\:0440\:0430\ \:0432\:0440\:0435\:043c\:0435\:043d\:0438\ \:043c\ \:043e\:043c\:0435\:043d\:0442\:0430\ \:0433\:043b\:0430\:0432\:043d\:043e\ \:0433\:043e\ \:0434\:0432\:0438\:0436\:0435\:043d\:0438\:044f\ \:043a\:043e\ \:043b\:0438\:0447\:0435\:0441\:0442\:0432\ \:043f\:0440\:043e\:0438\:0437\ \:0432\:043e\:0434\:043d\:0430\:044f\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(TableForm[pdtgv = D[K, t] + Cross[w, K]]\)\(\[IndentingNewLine]\) \)\)], "Input", CellChangeTimes -> {{3.497073305866211*^9, 3.497073356705078*^9}}], Cell[BoxData[ InterpretationBox[GridBox[{ { RowBox[{\(Ixy\ wx[t]\ wz[t]\), "-", \(Iyy\ wy[t]\ wz[t]\), "+", \(Izz\ wy[t]\ wz[t]\), "+", RowBox[{"Ixx", " ", RowBox[{ SuperscriptBox["wx", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "-", RowBox[{"Ixy", " ", RowBox[{ SuperscriptBox["wy", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}]}, { RowBox[{\(Ixx\ wx[t]\ wz[t]\), "-", \(Izz\ wx[t]\ wz[t]\), "-", \(Ixy\ wy[t]\ wz[t]\), "-", RowBox[{"Ixy", " ", RowBox[{ SuperscriptBox["wx", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "+", RowBox[{"Iyy", " ", RowBox[{ SuperscriptBox["wy", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}]}, { RowBox[{\(\(-Ixy\)\ wx[t]\^2\), "-", \(Ixx\ wx[t]\ wy[t]\), "+", \(Iyy\ wx[t]\ wy[t]\), "+", \(Ixy\ wy[t]\^2\), "+", RowBox[{"Izz", " ", RowBox[{ SuperscriptBox["wz", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}]} }, RowSpacings->1, ColumnSpacings->3, RowAlignments->Baseline, ColumnAlignments->{Left}], TableForm[ { Plus[ Times[ Ixy, wx[ t], wz[ t]], Times[ -1, Iyy, wy[ t], wz[ t]], Times[ Izz, wy[ t], wz[ t]], Times[ Ixx, Derivative[ 1][ wx][ t]], Times[ -1, Ixy, Derivative[ 1][ wy][ t]]], Plus[ Times[ Ixx, wx[ t], wz[ t]], Times[ -1, Izz, wx[ t], wz[ t]], Times[ -1, Ixy, wy[ t], wz[ t]], Times[ -1, Ixy, Derivative[ 1][ wx][ t]], Times[ Iyy, Derivative[ 1][ wy][ t]]], Plus[ Times[ -1, Ixy, Power[ wx[ t], 2]], Times[ -1, Ixx, wx[ t], wy[ t]], Times[ Iyy, wx[ t], wy[ t]], Times[ Ixy, Power[ wy[ t], 2]], Times[ Izz, Derivative[ 1][ wz][ t]]]}]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(q10 = pdtgv[\([1]\)] == M[\([1]\)]\)\)\)], "Input", CellChangeTimes -> {{3.4970733664179688`*^9, 3.49707341978125*^9}}], Cell[BoxData[ RowBox[{ RowBox[{\(Ixy\ wx[t]\ wz[t]\), "-", \(Iyy\ wy[t]\ wz[t]\), "+", \(Izz\ wy[t]\ wz[t]\), "+", RowBox[{"Ixx", " ", RowBox[{ SuperscriptBox["wx", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "-", RowBox[{"Ixy", " ", RowBox[{ SuperscriptBox["wy", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}], "\[Equal]", \(Mx[alfa[t], mah[t]]\)}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(q11 = pdtgv[\([2]\)] == M[\([2]\)]\)\)\)], "Input", CellChangeTimes -> {{3.497073430461914*^9, 3.4970734459541016`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{\(Ixx\ wx[t]\ wz[t]\), "-", \(Izz\ wx[t]\ wz[t]\), "-", \(Ixy\ wy[t]\ wz[t]\), "-", RowBox[{"Ixy", " ", RowBox[{ SuperscriptBox["wx", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "+", RowBox[{"Iyy", " ", RowBox[{ SuperscriptBox["wy", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}], "\[Equal]", \(My[alfa[t], mah[t]]\)}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(q12 = pdtgv[\([3]\)] == M[\([3]\)]\)\)\)], "Input", CellChangeTimes -> {{3.497073442741211*^9, 3.497073452970703*^9}}], Cell[BoxData[ RowBox[{ RowBox[{\(\(-Ixy\)\ wx[t]\^2\), "-", \(Ixx\ wx[t]\ wy[t]\), "+", \(Iyy\ wx[t]\ wy[t]\), "+", \(Ixy\ wy[t]\^2\), "+", RowBox[{"Izz", " ", RowBox[{ SuperscriptBox["wz", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}], "\[Equal]", \(Mz[alfa[t], mah[t]]\)}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(TableForm[{q10, \ q11, \ q12}]\)\)\)], "Input"], Cell[BoxData[ InterpretationBox[GridBox[{ { RowBox[{ RowBox[{\(Ixy\ wx[t]\ wz[t]\), "-", \(Iyy\ wy[t]\ wz[t]\), "+", \(Izz\ wy[t]\ wz[t]\), "+", RowBox[{"Ixx", " ", RowBox[{ SuperscriptBox["wx", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "-", RowBox[{"Ixy", " ", RowBox[{ SuperscriptBox["wy", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}], "\[Equal]", \(Mx[alfa[t], mah[t]]\)}]}, { RowBox[{ RowBox[{\(Ixx\ wx[t]\ wz[t]\), "-", \(Izz\ wx[t]\ wz[t]\), "-", \(Ixy\ wy[t]\ wz[t]\), "-", RowBox[{"Ixy", " ", RowBox[{ SuperscriptBox["wx", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "+", RowBox[{"Iyy", " ", RowBox[{ SuperscriptBox["wy", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}], "\[Equal]", \(My[alfa[t], mah[t]]\)}]}, { RowBox[{ RowBox[{\(\(-Ixy\)\ wx[t]\^2\), "-", \(Ixx\ wx[t]\ wy[t]\), "+", \(Iyy\ wx[t]\ wy[t]\), "+", \(Ixy\ wy[t]\^2\), "+", RowBox[{"Izz", " ", RowBox[{ SuperscriptBox["wz", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}], "\[Equal]", \(Mz[alfa[t], mah[t]]\)}]} }, RowSpacings->1, ColumnSpacings->3, RowAlignments->Baseline, ColumnAlignments->{Left}], TableForm[ { Equal[ Plus[ Times[ Ixy, wx[ t], wz[ t]], Times[ -1, Iyy, wy[ t], wz[ t]], Times[ Izz, wy[ t], wz[ t]], Times[ Ixx, Derivative[ 1][ wx][ t]], Times[ -1, Ixy, Derivative[ 1][ wy][ t]]], Mx[ alfa[ t], mah[ t]]], Equal[ Plus[ Times[ Ixx, wx[ t], wz[ t]], Times[ -1, Izz, wx[ t], wz[ t]], Times[ -1, Ixy, wy[ t], wz[ t]], Times[ -1, Ixy, Derivative[ 1][ wx][ t]], Times[ Iyy, Derivative[ 1][ wy][ t]]], My[ alfa[ t], mah[ t]]], Equal[ Plus[ Times[ -1, Ixy, Power[ wx[ t], 2]], Times[ -1, Ixx, wx[ t], wy[ t]], Times[ Iyy, wx[ t], wy[ t]], Times[ Ixy, Power[ wy[ t], 2]], Times[ Izz, Derivative[ 1][ wz][ t]]], Mz[ alfa[ t], mah[ t]]]}]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(\[IndentingNewLine]\)\(MatrixForm[ urav = {q1, \ q2, \ q3, \ q4, \ q5, \ q6, \ q7, \ q8, \ q9, \ q10, \ q11, \ q12}]\)\)\)], "Input"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", TagBox[GridBox[{ { RowBox[{ RowBox[{ SuperscriptBox["l", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", \(Cos[\[Theta][t]]\ Cos[\[Psi][t]]\ vgx[t] - Cos[\[Psi][t]]\ Sin[\[Theta][t]]\ vgy[t] + Sin[\[Psi][t]]\ vgz[t]\)}]}, { RowBox[{ RowBox[{ SuperscriptBox["h", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", \(Sin[\[Theta][t]]\ vgx[t] + Cos[\[Theta][t]]\ vgy[t]\)}]}, { RowBox[{ RowBox[{ SuperscriptBox["z", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", \(\(-Cos[\[Theta][t]]\)\ Sin[\[Psi][t]]\ vgx[ t] + Sin[\[Theta][t]]\ Sin[\[Psi][t]]\ vgy[t] + Cos[\[Psi][t]]\ vgz[t]\)}]}, { RowBox[{ RowBox[{"m", " ", RowBox[{"(", RowBox[{\(vz[t]\ wy[t]\), "-", \(vy[t]\ wz[t]\), "+", RowBox[{ SuperscriptBox["vx", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], ")"}]}], "\[Equal]", "fx"}]}, { RowBox[{ RowBox[{"m", " ", RowBox[{"(", RowBox[{\(\(-vz[t]\)\ wx[t]\), "+", \(vx[t]\ wz[t]\), "+", RowBox[{ SuperscriptBox["vy", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], ")"}]}], "\[Equal]", "fy"}]}, { RowBox[{ RowBox[{"m", " ", RowBox[{"(", RowBox[{\(vy[t]\ wx[t]\), "-", \(vx[t]\ wy[t]\), "+", RowBox[{ SuperscriptBox["vz", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], ")"}]}], "\[Equal]", "fz"}]}, { RowBox[{\(wx[t]\), "\[Equal]", RowBox[{ RowBox[{ SuperscriptBox["gam", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "+", RowBox[{\(Sin[teta[t]]\), " ", RowBox[{ SuperscriptBox["psi", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}]}]}, { RowBox[{\(wy[t]\), "\[Equal]", RowBox[{ RowBox[{\(Cos[gam[t]]\), " ", \(Cos[teta[t]]\), " ", RowBox[{ SuperscriptBox["psi", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "+", RowBox[{\(Sin[gam[t]]\), " ", RowBox[{ SuperscriptBox["teta", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}]}]}, { RowBox[{\(wz[t]\), "\[Equal]", RowBox[{ RowBox[{\(-Cos[teta[t]]\), " ", \(Sin[gam[t]]\), " ", RowBox[{ SuperscriptBox["psi", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "+", RowBox[{\(Cos[gam[t]]\), " ", RowBox[{ SuperscriptBox["teta", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}]}]}, { RowBox[{ RowBox[{\(Ixy\ wx[t]\ wz[t]\), "-", \(Iyy\ wy[t]\ wz[t]\), "+", \(Izz\ wy[t]\ wz[t]\), "+", RowBox[{"Ixx", " ", RowBox[{ SuperscriptBox["wx", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "-", RowBox[{"Ixy", " ", RowBox[{ SuperscriptBox["wy", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}], "\[Equal]", \(Mx[alfa[t], mah[t]]\)}]}, { RowBox[{ RowBox[{\(Ixx\ wx[t]\ wz[t]\), "-", \(Izz\ wx[t]\ wz[t]\), "-", \(Ixy\ wy[t]\ wz[t]\), "-", RowBox[{"Ixy", " ", RowBox[{ SuperscriptBox["wx", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "+", RowBox[{"Iyy", " ", RowBox[{ SuperscriptBox["wy", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}], "\[Equal]", \(My[alfa[t], mah[t]]\)}]}, { RowBox[{ RowBox[{\(\(-Ixy\)\ wx[t]\^2\), "-", \(Ixx\ wx[t]\ wy[t]\), "+", \(Iyy\ wx[t]\ wy[t]\), "+", \(Ixy\ wy[t]\^2\), "+", RowBox[{"Izz", " ", RowBox[{ SuperscriptBox["wz", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}], "\[Equal]", \(Mz[alfa[t], mah[t]]\)}]} }, RowSpacings->1, ColumnAlignments->{Left}], Column], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(vxq = Flatten[{D[r, t], \ D[v, t], \ D[w, t], \ D[ug, t]}]\)], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{ SuperscriptBox["l", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], ",", RowBox[{ SuperscriptBox["h", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], ",", RowBox[{ SuperscriptBox["z", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], ",", RowBox[{ SuperscriptBox["vx", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], ",", RowBox[{ SuperscriptBox["vy", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], ",", RowBox[{ SuperscriptBox["vz", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], ",", RowBox[{ SuperscriptBox["wx", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], ",", RowBox[{ SuperscriptBox["wy", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], ",", RowBox[{ SuperscriptBox["wz", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], ",", RowBox[{ SuperscriptBox["psi", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], ",", RowBox[{ SuperscriptBox["teta", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], ",", RowBox[{ SuperscriptBox["gam", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "}"}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(vq = Flatten[{r, v, \ w, \ ug}]\)\)\)], "Input"], Cell[BoxData[ \({l[t], h[t], z[t], vx[t], vy[t], vz[t], wx[t], wy[t], wz[t], psi[t], teta[t], gam[t]}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(MatrixForm[ fk = Flatten[Solve[urav, vxq, \ Sort \[Rule] False]]]\)], "Input"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", TagBox[GridBox[{ { RowBox[{ RowBox[{ SuperscriptBox["l", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Rule]", \(Cos[\[Theta][t]]\ Cos[\[Psi][t]]\ vgx[t] - Cos[\[Psi][t]]\ Sin[\[Theta][t]]\ vgy[t] + Sin[\[Psi][t]]\ vgz[t]\)}]}, { RowBox[{ RowBox[{ SuperscriptBox["h", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Rule]", \(Sin[\[Theta][t]]\ vgx[t] + Cos[\[Theta][t]]\ vgy[t]\)}]}, { RowBox[{ RowBox[{ SuperscriptBox["z", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Rule]", \(\(-Cos[\[Theta][t]]\)\ Sin[\[Psi][t]]\ vgx[ t] + Sin[\[Theta][t]]\ Sin[\[Psi][t]]\ vgy[t] + Cos[\[Psi][t]]\ vgz[t]\)}]}, { RowBox[{ RowBox[{ SuperscriptBox["vx", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Rule]", \(-\(\(\(-fx\) + m\ vz[t]\ wy[t] - m\ vy[t]\ wz[t]\)\/m\)\)}]}, { RowBox[{ RowBox[{ SuperscriptBox["vy", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Rule]", \(-\(\(\(-fy\) - m\ vz[t]\ wx[t] + m\ vx[t]\ wz[t]\)\/m\)\)}]}, { RowBox[{ RowBox[{ SuperscriptBox["vz", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Rule]", \(-\(\(\(-fz\) + m\ vy[t]\ wx[t] - m\ vx[t]\ wy[t]\)\/m\)\)}]}, { RowBox[{ RowBox[{ SuperscriptBox["wx", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Rule]", \(-\(\(Iyy\ Mx[alfa[t], mah[t]] + Ixy\ My[alfa[t], mah[t]] - Ixx\ Ixy\ wx[t]\ wz[t] - Ixy\ Iyy\ wx[t]\ wz[t] + Ixy\ Izz\ wx[t]\ wz[t] + Ixy\^2\ wy[t]\ wz[t] + Iyy\^2\ wy[t]\ wz[t] - Iyy\ Izz\ wy[t]\ wz[t]\)\/\(Ixy\^2 - Ixx\ Iyy\)\)\)}]}, { RowBox[{ RowBox[{ SuperscriptBox["wy", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Rule]", \(-\(\(Ixy\ Mx[alfa[t], mah[t]] + Ixx\ My[alfa[t], mah[t]] - Ixx\^2\ wx[t]\ wz[t] - Ixy\^2\ wx[t]\ wz[t] + Ixx\ Izz\ wx[t]\ wz[t] + Ixx\ Ixy\ wy[t]\ wz[t] + Ixy\ Iyy\ wy[t]\ wz[t] - Ixy\ Izz\ wy[t]\ wz[t]\)\/\(Ixy\^2 - Ixx\ Iyy\)\)\)}]}, { RowBox[{ RowBox[{ SuperscriptBox["wz", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Rule]", \(-\(\(\(-Mz[alfa[t], mah[t]]\) - Ixy\ wx[t]\^2 - Ixx\ wx[t]\ wy[t] + Iyy\ wx[t]\ wy[t] + Ixy\ wy[t]\^2\)\/Izz\)\)}]}, { RowBox[{ RowBox[{ SuperscriptBox["psi", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Rule]", \(Sec[ teta[t]]\ \((Cos[gam[t]]\ wy[t] - Sin[gam[t]]\ wz[t])\)\)}]}, { RowBox[{ RowBox[{ SuperscriptBox["teta", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Rule]", \(Sin[gam[t]]\ wy[t] + Cos[gam[t]]\ wz[t]\)}]}, { RowBox[{ RowBox[{ SuperscriptBox["gam", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Rule]", \(\(-\(1\/2\)\)\ Sec[ teta[t]]\ \((\(-2\)\ Cos[teta[t]]\ wx[t] - Sin[gam[t] - teta[t]]\ wy[t] + Sin[gam[t] + teta[t]]\ wy[t] - Cos[gam[t] - teta[t]]\ wz[t] + Cos[gam[t] + teta[t]]\ wz[t])\)\)}]} }, RowSpacings->1, ColumnAlignments->{Left}], Column], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(TableForm[ uravn = Table[ fk[\([i, 1]\)] \[Equal] fk[\([i, 2]\)], \ {i, 1, 12}]]\)], "Input"], Cell[BoxData[ InterpretationBox[GridBox[{ { RowBox[{ RowBox[{ SuperscriptBox["l", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", \(Cos[\[Theta][t]]\ Cos[\[Psi][t]]\ vgx[t] - Cos[\[Psi][t]]\ Sin[\[Theta][t]]\ vgy[t] + Sin[\[Psi][t]]\ vgz[t]\)}]}, { RowBox[{ RowBox[{ SuperscriptBox["h", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", \(Sin[\[Theta][t]]\ vgx[t] + Cos[\[Theta][t]]\ vgy[t]\)}]}, { RowBox[{ RowBox[{ SuperscriptBox["z", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", \(\(-Cos[\[Theta][t]]\)\ Sin[\[Psi][t]]\ vgx[t] + Sin[\[Theta][t]]\ Sin[\[Psi][t]]\ vgy[t] + Cos[\[Psi][t]]\ vgz[t]\)}]}, { RowBox[{ RowBox[{ SuperscriptBox["vx", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", \(-\(\(\(-fx\) + m\ vz[t]\ wy[t] - m\ vy[t]\ wz[t]\)\/m\)\)}]}, { RowBox[{ RowBox[{ SuperscriptBox["vy", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", \(-\(\(\(-fy\) - m\ vz[t]\ wx[t] + m\ vx[t]\ wz[t]\)\/m\)\)}]}, { RowBox[{ RowBox[{ SuperscriptBox["vz", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", \(-\(\(\(-fz\) + m\ vy[t]\ wx[t] - m\ vx[t]\ wy[t]\)\/m\)\)}]}, { RowBox[{ RowBox[{ SuperscriptBox["wx", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", \(-\(\(Iyy\ Mx[alfa[t], mah[t]] + Ixy\ My[alfa[t], mah[t]] - Ixx\ Ixy\ wx[t]\ wz[t] - Ixy\ Iyy\ wx[t]\ wz[t] + Ixy\ Izz\ wx[t]\ wz[t] + Ixy\^2\ wy[t]\ wz[t] + Iyy\^2\ wy[t]\ wz[t] - Iyy\ Izz\ wy[t]\ wz[t]\)\/\(Ixy\^2 - Ixx\ Iyy\)\)\)}]}, { RowBox[{ RowBox[{ SuperscriptBox["wy", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", \(-\(\(Ixy\ Mx[alfa[t], mah[t]] + Ixx\ My[alfa[t], mah[t]] - Ixx\^2\ wx[t]\ wz[t] - Ixy\^2\ wx[t]\ wz[t] + Ixx\ Izz\ wx[t]\ wz[t] + Ixx\ Ixy\ wy[t]\ wz[t] + Ixy\ Iyy\ wy[t]\ wz[t] - Ixy\ Izz\ wy[t]\ wz[t]\)\/\(Ixy\^2 - Ixx\ Iyy\)\)\)}]}, { RowBox[{ RowBox[{ SuperscriptBox["wz", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", \(-\(\(\(-Mz[alfa[t], mah[t]]\) - Ixy\ wx[t]\^2 - Ixx\ wx[t]\ wy[t] + Iyy\ wx[t]\ wy[t] + Ixy\ wy[t]\^2\)\/Izz\)\)}]}, { RowBox[{ RowBox[{ SuperscriptBox["psi", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", \(Sec[ teta[t]]\ \((Cos[gam[t]]\ wy[t] - Sin[gam[t]]\ wz[t])\)\)}]}, { RowBox[{ RowBox[{ SuperscriptBox["teta", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", \(Sin[gam[t]]\ wy[t] + Cos[gam[t]]\ wz[t]\)}]}, { RowBox[{ RowBox[{ SuperscriptBox["gam", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", \(\(-\(1\/2\)\)\ Sec[ teta[t]]\ \((\(-2\)\ Cos[teta[t]]\ wx[t] - Sin[gam[t] - teta[t]]\ wy[t] + Sin[gam[t] + teta[t]]\ wy[t] - Cos[gam[t] - teta[t]]\ wz[t] + Cos[gam[t] + teta[t]]\ wz[t])\)\)}]} }, RowSpacings->1, ColumnSpacings->3, RowAlignments->Baseline, ColumnAlignments->{Left}], TableForm[ { Equal[ Derivative[ 1][ l][ t], Plus[ Times[ Cos[ \[Theta][ t]], Cos[ \[Psi][ t]], vgx[ t]], Times[ -1, Cos[ \[Psi][ t]], Sin[ \[Theta][ t]], vgy[ t]], Times[ Sin[ \[Psi][ t]], vgz[ t]]]], Equal[ Derivative[ 1][ h][ t], Plus[ Times[ Sin[ \[Theta][ t]], vgx[ t]], Times[ Cos[ \[Theta][ t]], vgy[ t]]]], Equal[ Derivative[ 1][ z][ t], Plus[ Times[ -1, Cos[ \[Theta][ t]], Sin[ \[Psi][ t]], vgx[ t]], Times[ Sin[ \[Theta][ t]], Sin[ \[Psi][ t]], vgy[ t]], Times[ Cos[ \[Psi][ t]], vgz[ t]]]], Equal[ Derivative[ 1][ vx][ t], Times[ -1, Power[ m, -1], Plus[ Times[ -1, fx], Times[ m, vz[ t], wy[ t]], Times[ -1, m, vy[ t], wz[ t]]]]], Equal[ Derivative[ 1][ vy][ t], Times[ -1, Power[ m, -1], Plus[ Times[ -1, fy], Times[ -1, m, vz[ t], wx[ t]], Times[ m, vx[ t], wz[ t]]]]], Equal[ Derivative[ 1][ vz][ t], Times[ -1, Power[ m, -1], Plus[ Times[ -1, fz], Times[ m, vy[ t], wx[ t]], Times[ -1, m, vx[ t], wy[ t]]]]], Equal[ Derivative[ 1][ wx][ t], Times[ -1, Power[ Plus[ Power[ Ixy, 2], Times[ -1, Ixx, Iyy]], -1], Plus[ Times[ Iyy, Mx[ alfa[ t], mah[ t]]], Times[ Ixy, My[ alfa[ t], mah[ t]]], Times[ -1, Ixx, Ixy, wx[ t], wz[ t]], Times[ -1, Ixy, Iyy, wx[ t], wz[ t]], Times[ Ixy, Izz, wx[ t], wz[ t]], Times[ Power[ Ixy, 2], wy[ t], wz[ t]], Times[ Power[ Iyy, 2], wy[ t], wz[ t]], Times[ -1, Iyy, Izz, wy[ t], wz[ t]]]]], Equal[ Derivative[ 1][ wy][ t], Times[ -1, Power[ Plus[ Power[ Ixy, 2], Times[ -1, Ixx, Iyy]], -1], Plus[ Times[ Ixy, Mx[ alfa[ t], mah[ t]]], Times[ Ixx, My[ alfa[ t], mah[ t]]], Times[ -1, Power[ Ixx, 2], wx[ t], wz[ t]], Times[ -1, Power[ Ixy, 2], wx[ t], wz[ t]], Times[ Ixx, Izz, wx[ t], wz[ t]], Times[ Ixx, Ixy, wy[ t], wz[ t]], Times[ Ixy, Iyy, wy[ t], wz[ t]], Times[ -1, Ixy, Izz, wy[ t], wz[ t]]]]], Equal[ Derivative[ 1][ wz][ t], Times[ -1, Power[ Izz, -1], Plus[ Times[ -1, Mz[ alfa[ t], mah[ t]]], Times[ -1, Ixy, Power[ wx[ t], 2]], Times[ -1, Ixx, wx[ t], wy[ t]], Times[ Iyy, wx[ t], wy[ t]], Times[ Ixy, Power[ wy[ t], 2]]]]], Equal[ Derivative[ 1][ psi][ t], Times[ Sec[ teta[ t]], Plus[ Times[ Cos[ gam[ t]], wy[ t]], Times[ -1, Sin[ gam[ t]], wz[ t]]]]], Equal[ Derivative[ 1][ teta][ t], Plus[ Times[ Sin[ gam[ t]], wy[ t]], Times[ Cos[ gam[ t]], wz[ t]]]], Equal[ Derivative[ 1][ gam][ t], Times[ Rational[ -1, 2], Sec[ teta[ t]], Plus[ Times[ -2, Cos[ teta[ t]], wx[ t]], Times[ -1, Sin[ Plus[ gam[ t], Times[ -1, teta[ t]]]], wy[ t]], Times[ Sin[ Plus[ gam[ t], teta[ t]]], wy[ t]], Times[ -1, Cos[ Plus[ gam[ t], Times[ -1, teta[ t]]]], wz[ t]], Times[ Cos[ Plus[ gam[ t], teta[ t]]], wz[ t]]]]]}]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(ReplaceAll[ uravn, \ {vgy[t] \[Rule] 0, vgz[t] \[Rule] 0, \ vx[t] \[Rule] 0, \ vy[t] \[Rule] 0, \ vz[t] \[Rule] 0, \ wx[t] \[Rule] 0, \ wy[t] \[Rule] 0, \ wz[t] \[Rule] o, \ psi[t] \[Rule] 0, \ teta[t] \[Rule] 0, \ gam[t] \[Rule] 0, \ fx \[Rule] 0, \ fy \[Rule] 0, \ fz \[Rule] 0, \ Mx[alfa[t], mah[t]] \[Rule] 0, \ My[alfa[t], mah[t]] \[Rule] 0, \ Mz[alfa[t], mah[t]] \[Rule] 0}]\)], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ SuperscriptBox["l", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", \(Cos[\[Theta][t]]\ Cos[\[Psi][t]]\ vgx[t]\)}], ",", RowBox[{ RowBox[{ SuperscriptBox["h", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", \(Sin[\[Theta][t]]\ vgx[t]\)}], ",", RowBox[{ RowBox[{ SuperscriptBox["z", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", \(\(-Cos[\[Theta][t]]\)\ Sin[\[Psi][t]]\ vgx[t]\)}], ",", RowBox[{ RowBox[{ SuperscriptBox["vx", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{ SuperscriptBox["vy", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{ SuperscriptBox["vz", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{ SuperscriptBox["wx", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{ SuperscriptBox["wy", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{ SuperscriptBox["wz", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{ SuperscriptBox["psi", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{ SuperscriptBox["teta", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", "o"}], ",", RowBox[{ RowBox[{ SuperscriptBox["gam", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", "0"}]}], "}"}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(MatrixForm[ matrA1 = Table[ D[uravn[\([i, 2]\)], vq[\([j]\)]], \ {i, 1, 12}, \ {j, 1, 12}]]\)\)\)], "Input"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0"}, {"0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0"}, {"0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0"}, {"0", "0", "0", "0", \(wz[t]\), \(-wy[t]\), "0", \(-vz[t]\), \(vy[t]\), "0", "0", "0"}, {"0", "0", "0", \(-wz[t]\), "0", \(wx[t]\), \(vz[t]\), "0", \(-vx[t]\), "0", "0", "0"}, {"0", "0", "0", \(wy[t]\), \(-wx[t]\), "0", \(-vy[t]\), \(vx[t]\), "0", "0", "0", "0"}, {"0", "0", "0", "0", "0", "0", \(-\(\(\(-Ixx\)\ Ixy\ wz[t] - Ixy\ Iyy\ wz[t] + Ixy\ Izz\ wz[t]\)\/\(Ixy\^2 - Ixx\ Iyy\)\)\), \(-\(\(Ixy\^2\ wz[t] + Iyy\^2\ wz[t] - Iyy\ Izz\ wz[t]\)\/\(Ixy\^2 - Ixx\ Iyy\)\)\), \(-\(\(\(-Ixx\)\ Ixy\ wx[t] - Ixy\ Iyy\ wx[t] + Ixy\ Izz\ wx[t] + Ixy\^2\ wy[t] + Iyy\^2\ wy[t] - Iyy\ Izz\ wy[t]\)\/\(Ixy\^2 - Ixx\ Iyy\)\)\), "0", "0", "0"}, {"0", "0", "0", "0", "0", "0", \(-\(\(\(-Ixx\^2\)\ wz[t] - Ixy\^2\ wz[t] + Ixx\ Izz\ wz[t]\)\/\(Ixy\^2 - Ixx\ Iyy\)\)\), \(-\(\(Ixx\ Ixy\ wz[t] + Ixy\ Iyy\ wz[t] - Ixy\ Izz\ wz[t]\)\/\(Ixy\^2 - Ixx\ Iyy\)\)\), \(-\(\(\(-Ixx\^2\)\ wx[t] - Ixy\^2\ wx[t] + Ixx\ Izz\ wx[t] + Ixx\ Ixy\ wy[t] + Ixy\ Iyy\ wy[t] - Ixy\ Izz\ wy[t]\)\/\(Ixy\^2 - Ixx\ Iyy\)\)\), "0", "0", "0"}, {"0", "0", "0", "0", "0", "0", \(-\(\(\(-2\)\ Ixy\ wx[t] - Ixx\ wy[t] + Iyy\ wy[t]\)\/Izz\)\), \(-\(\(\(-Ixx\)\ wx[t] + Iyy\ wx[t] + 2\ Ixy\ wy[t]\)\/Izz\)\), "0", "0", "0", "0"}, {"0", "0", "0", "0", "0", "0", "0", \(Cos[gam[t]]\ Sec[teta[t]]\), \(\(-Sec[teta[t]]\)\ Sin[ gam[t]]\), "0", \(Sec[teta[t]]\ Tan[ teta[t]]\ \((Cos[gam[t]]\ wy[t] - Sin[gam[t]]\ wz[t])\)\), \(Sec[ teta[t]]\ \((\(-Sin[gam[t]]\)\ wy[t] - Cos[gam[t]]\ wz[t])\)\)}, {"0", "0", "0", "0", "0", "0", "0", \(Sin[gam[t]]\), \(Cos[gam[t]]\), "0", "0", \(Cos[gam[t]]\ wy[t] - Sin[gam[t]]\ wz[t]\)}, {"0", "0", "0", "0", "0", "0", "1", \(\(-\(1\/2\)\)\ Sec[ teta[t]]\ \((\(-Sin[gam[t] - teta[t]]\) + Sin[gam[t] + teta[t]])\)\), \(\(-\(1\/2\)\)\ \((\(-Cos[ gam[t] - teta[t]]\) + Cos[gam[t] + teta[t]])\)\ Sec[ teta[t]]\), "0", \(\(-\(1\/2\)\)\ Sec[teta[t]]\ Tan[ teta[t]]\ \((\(-2\)\ Cos[teta[t]]\ wx[t] - Sin[gam[t] - teta[t]]\ wy[t] + Sin[gam[t] + teta[t]]\ wy[t] - Cos[gam[t] - teta[t]]\ wz[t] + Cos[gam[t] + teta[t]]\ wz[t])\) - 1\/2\ Sec[ teta[t]]\ \((2\ Sin[teta[t]]\ wx[t] + Cos[gam[t] - teta[t]]\ wy[t] + Cos[gam[t] + teta[t]]\ wy[t] - Sin[gam[t] - teta[t]]\ wz[t] - Sin[gam[t] + teta[t]]\ wz[t])\)\), \(\(-\(1\/2\)\)\ Sec[ teta[t]]\ \((\(-Cos[gam[t] - teta[t]]\)\ wy[t] + Cos[gam[t] + teta[t]]\ wy[t] + Sin[gam[t] - teta[t]]\ wz[t] - Sin[gam[t] + teta[t]]\ wz[t])\)\)} }, RowSpacings->1, ColumnSpacings->1, ColumnAlignments->{Left}], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(MatrixForm[ matrA2 = Table[ D[uravn[\([i, 2]\)], \ vq[\([j]\)]], \ {i, \ 4, 12}, \ {j, \ 4, \ 12}]]\)\)\)], "Input"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"0", \(wz[t]\), \(-wy[t]\), "0", \(-vz[t]\), \(vy[t]\), "0", "0", "0"}, {\(-wz[t]\), "0", \(wx[t]\), \(vz[t]\), "0", \(-vx[t]\), "0", "0", "0"}, {\(wy[t]\), \(-wx[t]\), "0", \(-vy[t]\), \(vx[t]\), "0", "0", "0", "0"}, {"0", "0", "0", \(-\(\(\(-Ixx\)\ Ixy\ wz[t] - Ixy\ Iyy\ wz[t] + Ixy\ Izz\ wz[t]\)\/\(Ixy\^2 - Ixx\ Iyy\)\)\), \(-\(\(Ixy\^2\ wz[t] + Iyy\^2\ wz[t] - Iyy\ Izz\ wz[t]\)\/\(Ixy\^2 - Ixx\ Iyy\)\)\), \(-\(\(\(-Ixx\)\ Ixy\ wx[t] - Ixy\ Iyy\ wx[t] + Ixy\ Izz\ wx[t] + Ixy\^2\ wy[t] + Iyy\^2\ wy[t] - Iyy\ Izz\ wy[t]\)\/\(Ixy\^2 - Ixx\ Iyy\)\)\), "0", "0", "0"}, {"0", "0", "0", \(-\(\(\(-Ixx\^2\)\ wz[t] - Ixy\^2\ wz[t] + Ixx\ Izz\ wz[t]\)\/\(Ixy\^2 - Ixx\ Iyy\)\)\), \(-\(\(Ixx\ Ixy\ wz[t] + Ixy\ Iyy\ wz[t] - Ixy\ Izz\ wz[t]\)\/\(Ixy\^2 - Ixx\ Iyy\)\)\), \(-\(\(\(-Ixx\^2\)\ wx[t] - Ixy\^2\ wx[t] + Ixx\ Izz\ wx[t] + Ixx\ Ixy\ wy[t] + Ixy\ Iyy\ wy[t] - Ixy\ Izz\ wy[t]\)\/\(Ixy\^2 - Ixx\ Iyy\)\)\), "0", "0", "0"}, {"0", "0", "0", \(-\(\(\(-2\)\ Ixy\ wx[t] - Ixx\ wy[t] + Iyy\ wy[t]\)\/Izz\)\), \(-\(\(\(-Ixx\)\ wx[t] + Iyy\ wx[t] + 2\ Ixy\ wy[t]\)\/Izz\)\), "0", "0", "0", "0"}, {"0", "0", "0", "0", \(Cos[gam[t]]\ Sec[teta[t]]\), \(\(-Sec[teta[t]]\)\ Sin[ gam[t]]\), "0", \(Sec[teta[t]]\ Tan[ teta[t]]\ \((Cos[gam[t]]\ wy[t] - Sin[gam[t]]\ wz[t])\)\), \(Sec[ teta[t]]\ \((\(-Sin[gam[t]]\)\ wy[t] - Cos[gam[t]]\ wz[t])\)\)}, {"0", "0", "0", "0", \(Sin[gam[t]]\), \(Cos[gam[t]]\), "0", "0", \(Cos[gam[t]]\ wy[t] - Sin[gam[t]]\ wz[t]\)}, {"0", "0", "0", "1", \(\(-\(1\/2\)\)\ Sec[ teta[t]]\ \((\(-Sin[gam[t] - teta[t]]\) + Sin[gam[t] + teta[t]])\)\), \(\(-\(1\/2\)\)\ \((\(-Cos[ gam[t] - teta[t]]\) + Cos[gam[t] + teta[t]])\)\ Sec[ teta[t]]\), "0", \(\(-\(1\/2\)\)\ Sec[teta[t]]\ Tan[ teta[t]]\ \((\(-2\)\ Cos[teta[t]]\ wx[t] - Sin[gam[t] - teta[t]]\ wy[t] + Sin[gam[t] + teta[t]]\ wy[t] - Cos[gam[t] - teta[t]]\ wz[t] + Cos[gam[t] + teta[t]]\ wz[t])\) - 1\/2\ Sec[ teta[t]]\ \((2\ Sin[teta[t]]\ wx[t] + Cos[gam[t] - teta[t]]\ wy[t] + Cos[gam[t] + teta[t]]\ wy[t] - Sin[gam[t] - teta[t]]\ wz[t] - Sin[gam[t] + teta[t]]\ wz[t])\)\), \(\(-\(1\/2\)\)\ Sec[ teta[t]]\ \((\(-Cos[gam[t] - teta[t]]\)\ wy[t] + Cos[gam[t] + teta[t]]\ wy[t] + Sin[gam[t] - teta[t]]\ wz[t] - Sin[gam[t] + teta[t]]\ wz[t])\)\)} }, RowSpacings->1, ColumnSpacings->1, ColumnAlignments->{Left}], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(upr = {dv[t], dh[t], \ de[t], \ dsg[t]}\)], "Input"], Cell[BoxData[ \({dv[t], dh[t], de[t], dsg[t]}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(MatrixForm[ matrB1 = Table[ D[uravn[\([i, 2]\)], \ upr[\([j]\)]], \ {i, \ 1, \ 12}, \ {j, \ 1, \ 4}]]\)], "Input"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"0", "0", "0", "0"}, {"0", "0", "0", "0"}, {"0", "0", "0", "0"}, {"0", "0", "0", "0"}, {"0", "0", "0", "0"}, {"0", "0", "0", "0"}, {"0", "0", "0", "0"}, {"0", "0", "0", "0"}, {"0", "0", "0", "0"}, {"0", "0", "0", "0"}, {"0", "0", "0", "0"}, {"0", "0", "0", "0"} }, RowSpacings->1, ColumnSpacings->1, ColumnAlignments->{Left}], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]], "Output"] }, Open ]] }, FrontEndVersion->"5.2 for Microsoft Windows", ScreenRectangle->{{0, 1280}, {0, 926}}, WindowSize->{1100, 873}, WindowMargins->{{58, Automatic}, {-19, Automatic}}, StyleDefinitions -> "Default.nb" ] (******************************************************************* Cached data follows. If you edit this Notebook file directly, not using Mathematica, you must remove the line containing CacheID at the top of the file. The cache data will then be recreated when you save this file from within Mathematica. *******************************************************************) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[1776, 53, 430, 8, 56, "Input"], Cell[2209, 63, 324, 5, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[2570, 73, 83, 1, 29, "Input"], Cell[2656, 76, 58, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[2751, 82, 76, 1, 29, "Input"], Cell[2830, 85, 52, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[2919, 91, 106, 2, 29, "Input"], Cell[3028, 95, 263, 8, 43, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[3328, 108, 82, 1, 30, "Input"], Cell[3413, 111, 263, 8, 43, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[3713, 124, 82, 1, 30, "Input"], Cell[3798, 127, 263, 8, 43, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[4098, 140, 60, 1, 30, "Input"], Cell[4161, 143, 891, 26, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[5089, 174, 167, 3, 30, "Input"], Cell[5259, 179, 389, 11, 71, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[5685, 195, 236, 4, 48, "Input"], Cell[5924, 201, 381, 11, 71, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[6342, 217, 73, 1, 29, "Input"], Cell[6418, 220, 583, 14, 71, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[7038, 239, 141, 3, 29, "Input"], Cell[7182, 244, 320, 11, 71, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[7539, 260, 62, 1, 29, "Input"], Cell[7604, 263, 359, 6, 48, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[8000, 274, 132, 2, 48, "Input"], Cell[8135, 278, 298, 7, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[8470, 290, 132, 2, 48, "Input"], Cell[8605, 294, 232, 6, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[8874, 305, 101, 2, 29, "Input"], Cell[8978, 309, 303, 7, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[9318, 321, 117, 2, 48, "Input"], Cell[9438, 325, 2358, 79, 71, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[11833, 409, 583, 11, 152, "Input"], Cell[12419, 422, 310, 4, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[12766, 431, 79, 1, 29, "Input"], Cell[12848, 434, 55, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[12940, 440, 79, 1, 29, "Input"], Cell[13022, 443, 55, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[13114, 449, 70, 1, 29, "Input"], Cell[13187, 452, 46, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[13270, 458, 125, 2, 48, "Input"], Cell[13398, 462, 876, 21, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[14311, 488, 79, 1, 29, "Input"], Cell[14393, 491, 314, 8, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[14744, 504, 117, 2, 48, "Input"], Cell[14864, 508, 319, 8, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[15220, 521, 117, 2, 48, "Input"], Cell[15340, 525, 314, 8, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[15691, 538, 117, 2, 48, "Input"], Cell[15811, 542, 2073, 64, 71, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[17921, 611, 555, 10, 152, "Input"], Cell[18479, 623, 310, 4, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[18826, 632, 124, 2, 48, "Input"], Cell[18953, 636, 59, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[19049, 642, 93, 1, 50, "Input"], Cell[19145, 645, 55, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[19237, 651, 230, 5, 48, "Input"], Cell[19470, 658, 369, 11, 71, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[19876, 674, 337, 5, 48, "Input"], Cell[20216, 681, 373, 11, 71, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[20626, 697, 265, 4, 48, "Input"], Cell[20894, 703, 369, 11, 71, "Output"] }, Open ]], Cell[21278, 717, 128, 2, 48, "Input"], Cell[CellGroupData[{ Cell[21431, 723, 548, 8, 32, "Input"], Cell[21982, 733, 742, 10, 22, "Message"], Cell[22727, 745, 517, 9, 22, "Message"] }, Open ]], Cell[23259, 757, 125, 2, 50, "Input"], Cell[CellGroupData[{ Cell[23409, 763, 209, 3, 30, "Input"], Cell[23621, 768, 615, 14, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[24273, 787, 186, 3, 50, "Input"], Cell[24462, 792, 410, 10, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[24909, 807, 167, 3, 50, "Input"], Cell[25079, 812, 209, 6, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[25325, 823, 162, 2, 50, "Input"], Cell[25490, 827, 1187, 29, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[26714, 861, 208, 3, 30, "Input"], Cell[26925, 866, 365, 9, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[27327, 880, 226, 4, 50, "Input"], Cell[27556, 886, 435, 10, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[28028, 901, 221, 4, 50, "Input"], Cell[28252, 907, 436, 10, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[28725, 922, 91, 1, 50, "Input"], Cell[28819, 925, 2492, 73, 71, "Output"] }, Open ]], Cell[31326, 1001, 99, 2, 50, "Input"], Cell[CellGroupData[{ Cell[31450, 1007, 612, 11, 178, "Input"], Cell[32065, 1020, 298, 4, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[32400, 1029, 416, 8, 50, "Input"], Cell[32819, 1039, 356, 11, 71, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[33212, 1055, 392, 5, 50, "Input"], Cell[33607, 1062, 576, 8, 22, "Message"], Cell[34186, 1072, 367, 7, 22, "Message"] }, Open ]], Cell[CellGroupData[{ Cell[34590, 1084, 218, 4, 50, "Input"], Cell[34811, 1090, 35, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[34883, 1096, 142, 2, 50, "Input"], Cell[35028, 1100, 37, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[35102, 1106, 150, 2, 50, "Input"], Cell[35255, 1110, 336, 11, 71, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[35628, 1126, 512, 8, 50, "Input"], Cell[36143, 1136, 369, 6, 22, "Message"], Cell[36515, 1144, 542, 11, 25, "Message"] }, Open ]], Cell[CellGroupData[{ Cell[37094, 1160, 156, 2, 50, "Input"], Cell[37253, 1164, 584, 22, 71, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[37874, 1191, 213, 3, 50, "Input"], Cell[38090, 1196, 104, 2, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[38231, 1203, 117, 2, 30, "Input"], Cell[38351, 1207, 118, 2, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[38506, 1214, 553, 8, 50, "Input"], Cell[39062, 1224, 450, 6, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[39549, 1235, 178, 3, 50, "Input"], Cell[39730, 1240, 2625, 82, 73, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[42392, 1327, 175, 3, 50, "Input"], Cell[42570, 1332, 503, 12, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[43110, 1349, 176, 3, 50, "Input"], Cell[43289, 1354, 503, 12, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[43829, 1371, 174, 3, 50, "Input"], Cell[44006, 1376, 370, 8, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[44413, 1389, 103, 2, 50, "Input"], Cell[44519, 1393, 3228, 100, 73, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[47784, 1498, 197, 3, 70, "Input"], Cell[47984, 1503, 5971, 131, 217, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[53992, 1639, 91, 1, 30, "Input"], Cell[54086, 1642, 1513, 38, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[55636, 1685, 102, 2, 50, "Input"], Cell[55741, 1689, 126, 2, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[55904, 1696, 104, 2, 30, "Input"], Cell[56011, 1700, 4957, 108, 297, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[61005, 1813, 129, 3, 30, "Input"], Cell[61137, 1818, 10680, 344, 297, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[71854, 2167, 473, 8, 50, "Input"], Cell[72330, 2177, 2250, 63, 48, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[74617, 2245, 184, 4, 50, "Input"], Cell[74804, 2251, 4011, 73, 263, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[78852, 2329, 192, 4, 50, "Input"], Cell[79047, 2335, 3639, 69, 215, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[82723, 2409, 72, 1, 30, "Input"], Cell[82798, 2412, 63, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[82898, 2418, 163, 4, 30, "Input"], Cell[83064, 2424, 678, 20, 215, "Output"] }, Open ]] } ] *) (******************************************************************* End of Mathematica Notebook file. *******************************************************************)