(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 10.2' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 158, 7] NotebookDataLength[ 80047, 2316] NotebookOptionsPosition[ 78013, 2276] NotebookOutlinePosition[ 78377, 2292] CellTagsIndexPosition[ 78334, 2289] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell["DOUBLE PENDULUM MODEL - CONVENTIONAL LAGRANGE", "Section", CellChangeTimes->{{3.738713096946988*^9, 3.738713097958741*^9}, { 3.738713210889555*^9, 3.738713211560526*^9}, {3.7736608585146313`*^9, 3.773660896170555*^9}, {3.773660926642432*^9, 3.773660927370582*^9}},ExpressionUUID->"ab0118e6-f1e9-40e8-94a5-\ 172a16b58603"], Cell[BoxData[ RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["r", "1"], "[", "t_", "]"}], " ", "=", " ", RowBox[{"{", RowBox[{ RowBox[{ SubscriptBox["a", "1"], " ", RowBox[{"Sin", "[", RowBox[{ SubscriptBox["\[Phi]", "1"], "[", "t", "]"}], "]"}]}], ",", " ", RowBox[{ SubscriptBox["a", RowBox[{"1", " "}]], RowBox[{"Cos", "[", RowBox[{ SubscriptBox["\[Phi]", "1"], "[", "t", "]"}], "]"}]}]}], "}"}]}], ";", "\n", RowBox[{ RowBox[{ SubscriptBox["r", "2"], "[", "t_", "]"}], " ", "=", " ", RowBox[{ RowBox[{ SubscriptBox["r", "1"], "[", "t", "]"}], " ", "+", " ", RowBox[{"{", RowBox[{ RowBox[{ SubscriptBox["a", "2"], " ", RowBox[{"Sin", "[", RowBox[{ RowBox[{ SubscriptBox["\[Phi]", "1"], "[", "t", "]"}], "+", RowBox[{ SubscriptBox["\[Phi]", "2"], "[", "t", "]"}]}], "]"}]}], ",", " ", RowBox[{ SubscriptBox["a", RowBox[{"2", " "}]], RowBox[{"Cos", "[", RowBox[{ RowBox[{ SubscriptBox["\[Phi]", "1"], "[", "t", "]"}], "+", RowBox[{ SubscriptBox["\[Phi]", "2"], "[", "t", "]"}]}], "]"}]}]}], "}"}]}]}], ";", "\n", RowBox[{"L", " ", "=", " ", RowBox[{ RowBox[{ FractionBox["1", "2"], SubscriptBox["m", "1"], RowBox[{ RowBox[{ RowBox[{ SubscriptBox["r", "1"], "'"}], "[", "t", "]"}], ".", RowBox[{ RowBox[{ SubscriptBox["r", "1"], "'"}], "[", "t", "]"}]}]}], " ", "+", " ", RowBox[{ FractionBox["1", "2"], SubscriptBox["m", "2"], RowBox[{ RowBox[{ RowBox[{ SubscriptBox["r", "2"], "'"}], "[", "t", "]"}], ".", RowBox[{ RowBox[{ SubscriptBox["r", "2"], "'"}], "[", "t", "]"}]}]}], " ", "+", " ", RowBox[{ SubscriptBox["m", "1"], " ", RowBox[{ RowBox[{"{", RowBox[{"0", ",", " ", "g"}], "}"}], ".", RowBox[{ SubscriptBox["r", "1"], "[", "t", "]"}]}]}], " ", "+", " ", RowBox[{ SubscriptBox["m", "2"], " ", RowBox[{ RowBox[{"{", RowBox[{"0", ",", " ", "g"}], "}"}], ".", RowBox[{ SubscriptBox["r", "2"], "[", "t", "]"}]}]}]}]}], ";", "\n", RowBox[{"rq", "=", " ", RowBox[{ RowBox[{"Flatten", "@", RowBox[{"Solve", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"#", "\[Equal]", "0"}], ")"}], "&"}], "/@", " ", RowBox[{"(", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"D", "[", RowBox[{"L", ",", RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["\[Phi]", "1"], "'"}], "[", "t", "]"}], ",", " ", RowBox[{ RowBox[{ SubscriptBox["\[Phi]", "2"], "'"}], "[", "t", "]"}]}], "}"}], "}"}]}], "]"}], ",", " ", "t"}], "]"}], " ", "-", " ", RowBox[{"D", "[", RowBox[{"L", ",", RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{ SubscriptBox["\[Phi]", "1"], "[", "t", "]"}], ",", " ", RowBox[{ SubscriptBox["\[Phi]", "2"], "[", "t", "]"}]}], "}"}], "}"}]}], "]"}]}], ")"}]}], ",", " ", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["\[Phi]", "1"], "''"}], "[", "t", "]"}], ",", " ", RowBox[{ RowBox[{ SubscriptBox["\[Phi]", "2"], "''"}], "[", "t", "]"}]}], "}"}]}], "]"}]}], "//", "Simplify"}]}]}], "*)"}]], "Code", CellChangeTimes->{{3.738713168571104*^9, 3.738713354470234*^9}, { 3.738713411872229*^9, 3.7387134304739733`*^9}, {3.7387134910857687`*^9, 3.738713734698469*^9}, {3.738713789938566*^9, 3.738713803556654*^9}, { 3.738714303820204*^9, 3.73871435213754*^9}, {3.738716427220037*^9, 3.738716536214409*^9}, {3.7387165662598886`*^9, 3.7387166045869493`*^9}, { 3.738716840187546*^9, 3.738716842172291*^9}, 3.773660964116867*^9}, CellLabel->"In[1]:=",ExpressionUUID->"8a0291bb-a54c-4f63-97d8-90f5c57b242d"] }, Open ]], Cell[CellGroupData[{ Cell["DOUBLE PENDULUM MODEL - MODULAR MODELING METHODOLOGY", "Section", CellChangeTimes->{{3.7387147230374804`*^9, 3.738714726180401*^9}, { 3.7736609019673147`*^9, 3.77366093016317*^9}},ExpressionUUID->"b23ed7d1-72d4-433a-8ada-\ b443518b34b0"], Cell[BoxData[ RowBox[{"(*", RowBox[{ RowBox[{"\[DoubleStruckCapitalM]", " ", "=", " ", RowBox[{"(", GridBox[{ {"0", "0", "0", "0", "0", "0"}, {"0", SubscriptBox["m", "1"], "0", "0", "0", "0"}, {"0", "0", SubscriptBox["m", "1"], "0", "0", "0"}, {"0", "0", "0", "0", "0", "0"}, {"0", "0", "0", "0", SubscriptBox["m", "2"], "0"}, {"0", "0", "0", "0", "0", SubscriptBox["m", "2"]} }], ")"}]}], ";", " ", RowBox[{"\[DoubleStruckF]", " ", "=", " ", RowBox[{"(", GridBox[{ {"0"}, {"0"}, { RowBox[{ SubscriptBox["m", "1"], "g"}]}, {"0"}, {"0"}, { RowBox[{ SubscriptBox["m", "2"], "g"}]} }], ")"}]}], ";", "\n", RowBox[{ RowBox[{ SubscriptBox["p", "0"], "[", "t_", "]"}], "=", RowBox[{"{", RowBox[{"0", ",", " ", RowBox[{"-", "1"}], ",", " ", "0"}], "}"}]}], ";", " ", RowBox[{ RowBox[{ SubscriptBox["p", "1"], "[", "t_", "]"}], "=", RowBox[{"{", RowBox[{"0", ",", " ", "0", ",", " ", "0"}], "}"}]}], ";", " ", RowBox[{ RowBox[{ SubscriptBox["p", "2"], "[", "t_", "]"}], "=", RowBox[{"{", RowBox[{ RowBox[{ SubscriptBox["x", "2"], "[", "t", "]"}], ",", RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}], ",", " ", "0"}], "}"}]}], ";", " ", RowBox[{ RowBox[{ SubscriptBox["p", "3"], "[", "t_", "]"}], "=", RowBox[{"{", RowBox[{ RowBox[{ SubscriptBox["x", "3"], "[", "t", "]"}], ",", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}], ",", " ", "0"}], "}"}]}], ";", "\n", RowBox[{ SubscriptBox["l", "0"], "=", "1"}], ";", " ", "\n", RowBox[{ RowBox[{"\[DoubleStruckQ]", "[", "t_", "]"}], " ", "=", " ", RowBox[{"{", RowBox[{ RowBox[{ SubscriptBox["\[Phi]", "1"], "[", "t", "]"}], ",", " ", RowBox[{ SubscriptBox["x", "2"], "[", "t", "]"}], ",", " ", RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}], ",", " ", RowBox[{ SubscriptBox["\[Phi]", "2"], "[", "t", "]"}], ",", " ", RowBox[{ SubscriptBox["x", "3"], "[", "t", "]"}], ",", " ", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}]}], "}"}]}], ";", "\n", RowBox[{ RowBox[{"\[DoubleStruckH]", "[", "t_", "]"}], " ", "=", " ", RowBox[{"{", "\n", "\t", RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ SubscriptBox["l", "0"], SubscriptBox["l", "1"]}], ")"}], RowBox[{"Cos", "[", RowBox[{ SubscriptBox["\[Phi]", "1"], "[", "t", "]"}], "]"}]}], "-", RowBox[{"Dot", "[", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ SubscriptBox["p", "2"], "[", "t", "]"}], "-", RowBox[{ SubscriptBox["p", "1"], "[", "t", "]"}]}], ")"}], ",", RowBox[{"(", RowBox[{ RowBox[{ SubscriptBox["p", "1"], "[", "t", "]"}], "-", RowBox[{ SubscriptBox["p", "0"], "[", "t", "]"}]}], ")"}]}], "]"}]}], ",", "\n", "\t", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ SubscriptBox["l", "0"], SubscriptBox["l", "1"]}], ")"}], RowBox[{"Sin", "[", RowBox[{ SubscriptBox["\[Phi]", "1"], "[", "t", "]"}], "]"}]}], "-", RowBox[{"(", RowBox[{ RowBox[{"Cross", "[", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ SubscriptBox["p", "2"], "[", "t", "]"}], "-", RowBox[{ SubscriptBox["p", "1"], "[", "t", "]"}]}], ")"}], ",", RowBox[{"(", RowBox[{ RowBox[{ SubscriptBox["p", "1"], "[", "t", "]"}], "-", RowBox[{ SubscriptBox["p", "0"], "[", "t", "]"}]}], ")"}]}], "]"}], "\[LeftDoubleBracket]", "3", "\[RightDoubleBracket]"}], ")"}]}], ",", "\n", "\t", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ SubscriptBox["l", "1"], SubscriptBox["l", "2"]}], ")"}], RowBox[{"Cos", "[", RowBox[{ SubscriptBox["\[Phi]", "2"], "[", "t", "]"}], "]"}]}], "-", RowBox[{"Dot", "[", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ SubscriptBox["p", "3"], "[", "t", "]"}], "-", RowBox[{ SubscriptBox["p", "2"], "[", "t", "]"}]}], ")"}], ",", RowBox[{"(", RowBox[{ RowBox[{ SubscriptBox["p", "2"], "[", "t", "]"}], "-", RowBox[{ SubscriptBox["p", "1"], "[", "t", "]"}]}], ")"}]}], "]"}]}], ",", "\n", "\t", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ SubscriptBox["l", "1"], SubscriptBox["l", "2"]}], ")"}], RowBox[{"Sin", "[", RowBox[{ SubscriptBox["\[Phi]", "2"], "[", "t", "]"}], "]"}]}], "-", RowBox[{"(", RowBox[{ RowBox[{"Cross", "[", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ SubscriptBox["p", "3"], "[", "t", "]"}], "-", RowBox[{ SubscriptBox["p", "2"], "[", "t", "]"}]}], ")"}], ",", RowBox[{"(", RowBox[{ RowBox[{ SubscriptBox["p", "2"], "[", "t", "]"}], "-", RowBox[{ SubscriptBox["p", "1"], "[", "t", "]"}]}], ")"}]}], "]"}], "\[LeftDoubleBracket]", "3", "\[RightDoubleBracket]"}], ")"}]}]}], "\n", "\t", "}"}]}], ";", "\n", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"\[DoubleStruckCapitalA]", "[", "t_", "]"}], ",", " ", RowBox[{"\[DoubleStruckB]", "[", "t_", "]"}]}], "}"}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"Normal", "@", RowBox[{"{", RowBox[{ RowBox[{"Last", "[", "#", "]"}], ",", RowBox[{"-", RowBox[{"First", "[", "#", "]"}]}]}], "}"}]}], "&"}], " ", "@", " ", RowBox[{"CoefficientArrays", "[", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"\[DoubleStruckH]", "[", "t", "]"}], ",", " ", RowBox[{"{", RowBox[{"t", ",", "2"}], "}"}]}], "]"}], ",", " ", RowBox[{ RowBox[{"\[DoubleStruckQ]", "''"}], "[", "t", "]"}]}], "]"}]}]}], ";", "\n", RowBox[{ RowBox[{"\[DoubleStruckCapitalS]", "[", "t_", "]"}], " ", "=", " ", RowBox[{ RowBox[{"Transpose", " ", "@", " ", RowBox[{"NullSpace", "[", RowBox[{"\[DoubleStruckCapitalA]", "[", "t", "]"}], "]"}]}], "//", "Simplify"}]}], ";", "\n", RowBox[{"rl", " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"Join", "[", RowBox[{"#", ",", RowBox[{"D", "[", RowBox[{"#", ",", "t"}], "]"}], ",", RowBox[{"D", "[", RowBox[{"#", ",", RowBox[{"{", RowBox[{"t", ",", "2"}], "}"}]}], "]"}]}], "]"}], "&"}], " ", "@", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["x", "2"], "[", "t", "]"}], "\[Rule]", RowBox[{ SubscriptBox["a", "1"], " ", RowBox[{"Sin", "[", RowBox[{ SubscriptBox["\[Phi]", "1"], "[", "t", "]"}], "]"}]}]}], ",", RowBox[{ RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}], "\[Rule]", RowBox[{ SubscriptBox["a", "1"], " ", RowBox[{"Cos", "[", RowBox[{ SubscriptBox["\[Phi]", "1"], "[", "t", "]"}], "]"}]}]}], ",", RowBox[{ RowBox[{ SubscriptBox["x", "3"], "[", "t", "]"}], "\[Rule]", RowBox[{ RowBox[{ SubscriptBox["a", "1"], " ", RowBox[{"Sin", "[", RowBox[{ SubscriptBox["\[Phi]", "1"], "[", "t", "]"}], "]"}]}], "+", RowBox[{ SubscriptBox["a", "2"], " ", RowBox[{"Sin", "[", RowBox[{ RowBox[{ SubscriptBox["\[Phi]", "1"], "[", "t", "]"}], "+", RowBox[{ SubscriptBox["\[Phi]", "2"], "[", "t", "]"}]}], "]"}]}]}]}], ",", RowBox[{ RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}], "\[Rule]", RowBox[{ RowBox[{ SubscriptBox["a", "1"], " ", RowBox[{"Cos", "[", RowBox[{ SubscriptBox["\[Phi]", "1"], "[", "t", "]"}], "]"}]}], "+", RowBox[{ SubscriptBox["a", "2"], " ", RowBox[{"Cos", "[", RowBox[{ RowBox[{ SubscriptBox["\[Phi]", "1"], "[", "t", "]"}], "+", RowBox[{ SubscriptBox["\[Phi]", "2"], "[", "t", "]"}]}], "]"}]}]}]}]}], "}"}]}]}], ";", "\n", RowBox[{ RowBox[{ RowBox[{ RowBox[{"Transpose", "[", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"\[DoubleStruckCapitalS]", "[", "t", "]"}], "/.", "rl"}], "//", "Simplify"}], ")"}], "]"}], ".", RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"\[DoubleStruckCapitalM]", ".", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{ RowBox[{"\[DoubleStruckQ]", "''"}], "[", "t", "]"}], "/.", "rl"}], "/.", "rq"}], ")"}]}], "-", "\[DoubleStruckF]"}], ")"}], "//", "Simplify"}], ")"}]}], "//", "Simplify"}], "//", "Norm"}]}], "*)"}]], "Code", CellChangeTimes->{{3.738714733679265*^9, 3.738714757240458*^9}, { 3.738714793713504*^9, 3.73871508080786*^9}, {3.7387151164959917`*^9, 3.738715213637285*^9}, {3.738715319716114*^9, 3.738715378915904*^9}, { 3.7387154190004997`*^9, 3.738715498203505*^9}, {3.73871561748785*^9, 3.7387156233904467`*^9}, {3.7387157519143887`*^9, 3.738715862871901*^9}, { 3.738715894283777*^9, 3.7387159815447063`*^9}, {3.738716076815089*^9, 3.73871609334184*^9}, {3.738716131226635*^9, 3.738716203971139*^9}, { 3.738716234896759*^9, 3.7387162870986433`*^9}, {3.738716329796945*^9, 3.738716336869817*^9}, {3.7387163679095*^9, 3.738716393133065*^9}, { 3.738716613387546*^9, 3.738716619531884*^9}, {3.738716909208609*^9, 3.738716909985396*^9}, {3.738717236880988*^9, 3.738717269266539*^9}, { 3.738717322552346*^9, 3.738717323597872*^9}, 3.773660969768978*^9}, CellLabel->"In[2]:=",ExpressionUUID->"95515b98-4828-43a1-a754-b9ffaafce4ae"], Cell[BoxData[ RowBox[{"(*", RowBox[{ RowBox[{"\[DoubleStruckCapitalA]", "[", "t", "]"}], "//", "MatrixForm"}], "*)"}]], "Input", CellChangeTimes->{{3.738717207931938*^9, 3.7387172135684853`*^9}, 3.7736609749717913`*^9}, CellLabel->"In[3]:=",ExpressionUUID->"9aaccd05-a673-44f1-8df0-6ba7b7a41952"], Cell[BoxData[ RowBox[{"(*", RowBox[{ RowBox[{"\[DoubleStruckB]", "[", "t", "]"}], "//", "MatrixForm"}], "*)"}]], "Input", CellChangeTimes->{{3.7387172161157103`*^9, 3.738717222180008*^9}, 3.773660977726322*^9}, CellLabel->"In[4]:=",ExpressionUUID->"9c5f0a77-452d-411b-a09d-0e94b018cce1"], Cell[BoxData[ RowBox[{"(*", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"\[DoubleStruckCapitalS]", "[", "t", "]"}], "/.", "rl"}], "//", "Simplify"}], ")"}], "//", "MatrixForm"}], "*)"}]], "Input", CellChangeTimes->{{3.738716738941844*^9, 3.738716739871696*^9}, { 3.738716941528304*^9, 3.738716943718101*^9}, 3.773660981382303*^9}, CellLabel->"In[5]:=",ExpressionUUID->"a53523e0-5ac1-462d-8152-e417939d03b8"] }, Open ]], Cell[CellGroupData[{ Cell["n-PENDULUM MODEL - MODULAR MODELING METHODOLOGY", "Section", CellChangeTimes->{{3.7389576531280823`*^9, 3.738957690425064*^9}, 3.773660935790267*^9},ExpressionUUID->"89ecdc42-0a7d-421c-9292-\ 790d76cf9a3f"], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"n", " ", "=", " ", "4"}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"\[DoubleStruckQ]", "[", "t_", "]"}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ SubscriptBox["x", "#"], "[", "t", "]"}], ",", " ", RowBox[{ SubscriptBox["y", "#"], "[", "t", "]"}], ",", " ", RowBox[{ SubscriptBox["\[Phi]", "#"], "[", "t", "]"}]}], "}"}], " ", "&"}], " ", "/@", " ", RowBox[{"Range", "[", "n", "]"}]}], ")"}], " ", "//", " ", "Join"}], " ", "//", " ", "Flatten"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["\[DoubleStruckR]", "0"], "[", "t_", "]"}], " ", "=", " ", RowBox[{"{", RowBox[{"0", ",", " ", "0", ",", " ", "0"}], "}"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ SubscriptBox["\[DoubleStruckR]", "#"], "[", "t_", "]"}], " ", "=", " ", RowBox[{"{", RowBox[{ RowBox[{ SubscriptBox["x", "#"], "[", "t", "]"}], ",", " ", RowBox[{ SubscriptBox["y", "#"], "[", "t", "]"}], ",", " ", "0"}], "}"}]}], ")"}], " ", "&"}], " ", "/@", " ", RowBox[{"Range", "[", "n", "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"\[DoubleStruckH]", "[", "t_", "]"}], " ", "=", " ", RowBox[{"Join", "[", "\n", "\t", RowBox[{ RowBox[{"{", "\n", "\t\t", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["l", "1"], " ", RowBox[{"Cos", "[", RowBox[{ SubscriptBox["\[Phi]", "1"], "[", "t", "]"}], "]"}]}], "-", " ", RowBox[{ RowBox[{ SubscriptBox["\[DoubleStruckR]", "1"], "[", "t", "]"}], ".", RowBox[{"{", RowBox[{"0", ",", " ", "1", ",", " ", "0"}], "}"}]}]}], ",", " ", "\n", "\t\t", RowBox[{ RowBox[{ SubscriptBox["l", "1"], " ", RowBox[{"Sin", "[", RowBox[{ SubscriptBox["\[Phi]", "1"], "[", "t", "]"}], "]"}]}], "-", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ SubscriptBox["\[DoubleStruckR]", "1"], "[", "t", "]"}], "\[Cross]", RowBox[{"{", RowBox[{"0", ",", " ", "1", ",", " ", "0"}], "}"}]}], ")"}], ".", RowBox[{"{", RowBox[{"0", ",", " ", "0", ",", " ", "1"}], "}"}]}]}]}], "\n", "\t", "}"}], ",", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"{", "\n", "\t\t", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["l", "#"], SubscriptBox["l", RowBox[{"#", "-", "1"}]], " ", RowBox[{"Cos", "[", RowBox[{ SubscriptBox["\[Phi]", "#"], "[", "t", "]"}], "]"}]}], "-", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ SubscriptBox["\[DoubleStruckR]", "#"], "[", "t", "]"}], "-", RowBox[{ SubscriptBox["\[DoubleStruckR]", RowBox[{"#", "-", "1"}]], "[", "t", "]"}]}], ")"}], ".", RowBox[{"(", RowBox[{ RowBox[{ SubscriptBox["\[DoubleStruckR]", RowBox[{"#", "-", "1"}]], "[", "t", "]"}], "-", RowBox[{ SubscriptBox["\[DoubleStruckR]", RowBox[{"#", "-", "2"}]], "[", "t", "]"}]}], ")"}]}]}], ",", " ", "\n", "\t\t", RowBox[{ RowBox[{ SubscriptBox["l", "#"], SubscriptBox["l", RowBox[{"#", "-", "1"}]], " ", RowBox[{"Sin", "[", RowBox[{ SubscriptBox["\[Phi]", "#"], "[", "t", "]"}], "]"}]}], "-", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ SubscriptBox["\[DoubleStruckR]", "#"], "[", "t", "]"}], "-", RowBox[{ SubscriptBox["\[DoubleStruckR]", RowBox[{"#", "-", "1"}]], "[", "t", "]"}]}], ")"}], "\[Cross]", RowBox[{"(", RowBox[{ RowBox[{ SubscriptBox["\[DoubleStruckR]", RowBox[{"#", "-", "1"}]], "[", "t", "]"}], "-", RowBox[{ SubscriptBox["\[DoubleStruckR]", RowBox[{"#", "-", "2"}]], "[", "t", "]"}]}], ")"}]}], ")"}], ".", RowBox[{"{", RowBox[{"0", ",", " ", "0", ",", " ", "1"}], "}"}]}]}]}], "\n", "\t", "}"}], "&"}], " ", "/@", " ", RowBox[{"Range", "[", RowBox[{"2", ",", "n"}], "]"}]}], ")"}], " ", "//", " ", "Flatten"}]}], "\n", "\t", "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"\[DoubleStruckH]", "[", "t", "]"}], " ", "//", " ", "TableForm"}], "\n"}], "\n", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"\[DoubleStruckCapitalA]", "[", "t_", "]"}], ",", " ", RowBox[{"\[DoubleStruckB]", "[", "t_", "]"}]}], "}"}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"Normal", "@", RowBox[{"{", RowBox[{ RowBox[{"Last", "[", "#", "]"}], ",", RowBox[{"-", RowBox[{"First", "[", "#", "]"}]}]}], "}"}]}], "&"}], " ", "@", " ", RowBox[{"CoefficientArrays", "[", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"\[DoubleStruckH]", "[", "t", "]"}], ",", " ", RowBox[{"{", RowBox[{"t", ",", "2"}], "}"}]}], "]"}], ",", " ", RowBox[{ RowBox[{"\[DoubleStruckQ]", "''"}], "[", "t", "]"}]}], "]"}]}]}], ";"}], "\n", RowBox[{ RowBox[{"\[DoubleStruckCapitalA]", "[", "t", "]"}], " ", "//", " ", "MatrixForm"}], "\n", RowBox[{ RowBox[{ RowBox[{"\[DoubleStruckB]", "[", "t", "]"}], " ", "//", " ", "MatrixForm"}], "\n"}], "\n", RowBox[{ RowBox[{"\[ScriptCapitalL]", " ", "=", " ", RowBox[{"Sum", "[", RowBox[{ RowBox[{ FractionBox["1", "2"], " ", SubscriptBox["m", "j"], " ", RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"#", ".", "#"}], ")"}], "&"}], "@", RowBox[{"D", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"1", "-", SubscriptBox["\[Chi]", "j"]}], ")"}], RowBox[{ SubscriptBox["\[DoubleStruckR]", RowBox[{"j", "-", "1"}]], "[", "t", "]"}]}], "+", RowBox[{ SubscriptBox["\[Chi]", "j"], RowBox[{ SubscriptBox["\[DoubleStruckR]", "j"], "[", "t", "]"}]}]}], ",", "t"}], "]"}]}], ")"}], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", FractionBox[ SubscriptBox["\[Rho]", "j"], SubscriptBox["l", "j"]], ")"}], "2"], RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["\[DoubleStruckR]", "j"], "'"}], "[", "t", "]"}], "-", RowBox[{ RowBox[{ SubscriptBox["\[DoubleStruckR]", RowBox[{"j", "-", "1"}]], "'"}], "[", "t", "]"}]}], ")"}], ".", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["\[DoubleStruckR]", "j"], "'"}], "[", "t", "]"}], "-", RowBox[{ RowBox[{ SubscriptBox["\[DoubleStruckR]", RowBox[{"j", "-", "1"}]], "'"}], "[", "t", "]"}]}], ")"}]}]}], "+", RowBox[{"2", RowBox[{ RowBox[{"{", RowBox[{"0", ",", " ", "g", ",", " ", "0"}], "}"}], ".", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"1", "-", SubscriptBox["\[Chi]", "j"]}], ")"}], RowBox[{ SubscriptBox["\[DoubleStruckR]", RowBox[{"j", "-", "1"}]], "[", "t", "]"}]}], "+", RowBox[{ SubscriptBox["\[Chi]", "j"], RowBox[{ SubscriptBox["\[DoubleStruckR]", "j"], "[", "t", "]"}]}]}], ")"}]}]}]}], ")"}]}], ",", " ", RowBox[{"{", RowBox[{"j", ",", "n"}], "}"}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{ RowBox[{"\[DoubleStruckE]", "[", "t_", "]"}], " ", "=", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"D", "[", RowBox[{"\[ScriptCapitalL]", ",", RowBox[{"{", RowBox[{ RowBox[{"\[DoubleStruckQ]", "'"}], "[", "t", "]"}], "}"}]}], "]"}], ",", " ", "t"}], "]"}], " ", "-", " ", RowBox[{"D", "[", RowBox[{"\[ScriptCapitalL]", ",", RowBox[{"{", RowBox[{"\[DoubleStruckQ]", "[", "t", "]"}], "}"}]}], "]"}]}], ")"}], " ", "//", " ", "FullSimplify"}]}], ";"}], "\n", RowBox[{"(*", RowBox[{ RowBox[{"\[DoubleStruckE]", "[", "t", "]"}], " ", "//", " ", "TableForm"}], "*)"}], "\n"}], "\n", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"\[DoubleStruckCapitalM]", "[", "t_", "]"}], ",", " ", RowBox[{"\[DoubleStruckF]", "[", "t_", "]"}]}], "}"}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"FullSimplify", "@", RowBox[{"Normal", "@", RowBox[{"{", RowBox[{ RowBox[{"Last", "[", "#", "]"}], ",", RowBox[{"-", RowBox[{"First", "[", "#", "]"}]}]}], "}"}]}]}], "&"}], " ", "@", " ", RowBox[{"CoefficientArrays", "[", RowBox[{ RowBox[{"\[DoubleStruckE]", "[", "t", "]"}], ",", " ", RowBox[{ RowBox[{"\[DoubleStruckQ]", "''"}], "[", "t", "]"}]}], "]"}]}]}], ";"}], "\n", RowBox[{ RowBox[{"\[DoubleStruckCapitalM]", "[", "t", "]"}], " ", "//", " ", "MatrixForm"}], "\n", RowBox[{ RowBox[{"\[DoubleStruckF]", "[", "t", "]"}], " ", "//", " ", "MatrixForm"}]}], "Code", CellChangeTimes->{{3.7389576980532637`*^9, 3.7389577860664988`*^9}, { 3.7389578167710876`*^9, 3.738958107337634*^9}, {3.738958164309383*^9, 3.738958208194663*^9}, {3.738958279533186*^9, 3.73895841129618*^9}, { 3.738958486732791*^9, 3.738958493117483*^9}, {3.738958576325405*^9, 3.738958628494872*^9}, {3.738958731995323*^9, 3.738958753018186*^9}, { 3.7389587945497017`*^9, 3.738958801957364*^9}, {3.738958842040492*^9, 3.738958895247943*^9}, {3.738958929876376*^9, 3.7389590101741343`*^9}, { 3.738959051844101*^9, 3.7389590541969337`*^9}, {3.7389591746992693`*^9, 3.7389592281961403`*^9}, {3.7389593426994343`*^9, 3.738959439192668*^9}, { 3.738959526481279*^9, 3.7389595582775*^9}, {3.738959673895013*^9, 3.7389596741162043`*^9}, {3.7389610774890614`*^9, 3.7389610776456842`*^9}, {3.7389613273106413`*^9, 3.738961327492675*^9}, 3.77241231964491*^9, {3.772412456270812*^9, 3.7724124587156487`*^9}}, CellLabel->"In[6]:=",ExpressionUUID->"2c63b02f-3cfe-40f5-8f2c-f810b5ce01ec"], Cell[BoxData[ TagBox[ TagBox[GridBox[{ { RowBox[{ RowBox[{ RowBox[{"Cos", "[", RowBox[{ SubscriptBox["\[Phi]", "1"], "[", "t", "]"}], "]"}], " ", SubscriptBox["l", "1"]}], "-", RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}]}]}, { RowBox[{ RowBox[{ RowBox[{"Sin", "[", RowBox[{ SubscriptBox["\[Phi]", "1"], "[", "t", "]"}], "]"}], " ", SubscriptBox["l", "1"]}], "-", RowBox[{ SubscriptBox["x", "1"], "[", "t", "]"}]}]}, { RowBox[{ RowBox[{ RowBox[{"Cos", "[", RowBox[{ SubscriptBox["\[Phi]", "2"], "[", "t", "]"}], "]"}], " ", SubscriptBox["l", "1"], " ", SubscriptBox["l", "2"]}], "-", RowBox[{ RowBox[{ SubscriptBox["x", "1"], "[", "t", "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{ SubscriptBox["x", "1"], "[", "t", "]"}]}], "+", RowBox[{ SubscriptBox["x", "2"], "[", "t", "]"}]}], ")"}]}], "-", RowBox[{ RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}]}], "+", RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}]}], ")"}]}]}]}, { RowBox[{ RowBox[{ RowBox[{"Sin", "[", RowBox[{ SubscriptBox["\[Phi]", "2"], "[", "t", "]"}], "]"}], " ", SubscriptBox["l", "1"], " ", SubscriptBox["l", "2"]}], "-", RowBox[{ RowBox[{ SubscriptBox["x", "2"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}]}], "+", RowBox[{ RowBox[{ SubscriptBox["x", "1"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}]}]}]}, { RowBox[{ RowBox[{ RowBox[{"Cos", "[", RowBox[{ SubscriptBox["\[Phi]", "3"], "[", "t", "]"}], "]"}], " ", SubscriptBox["l", "2"], " ", SubscriptBox["l", "3"]}], "-", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{ SubscriptBox["x", "1"], "[", "t", "]"}]}], "+", RowBox[{ SubscriptBox["x", "2"], "[", "t", "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{ SubscriptBox["x", "2"], "[", "t", "]"}]}], "+", RowBox[{ SubscriptBox["x", "3"], "[", "t", "]"}]}], ")"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}]}], "+", RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}]}], "+", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}]}], ")"}]}]}]}, { RowBox[{ RowBox[{ RowBox[{"Sin", "[", RowBox[{ SubscriptBox["\[Phi]", "3"], "[", "t", "]"}], "]"}], " ", SubscriptBox["l", "2"], " ", SubscriptBox["l", "3"]}], "-", RowBox[{ RowBox[{ SubscriptBox["x", "2"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}]}], "+", RowBox[{ RowBox[{ SubscriptBox["x", "3"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}]}], "+", RowBox[{ RowBox[{ SubscriptBox["x", "1"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}]}], "-", RowBox[{ RowBox[{ SubscriptBox["x", "3"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}]}], "-", RowBox[{ RowBox[{ SubscriptBox["x", "1"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}]}], "+", RowBox[{ RowBox[{ SubscriptBox["x", "2"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}]}]}]}, { RowBox[{ RowBox[{ RowBox[{"Cos", "[", RowBox[{ SubscriptBox["\[Phi]", "4"], "[", "t", "]"}], "]"}], " ", SubscriptBox["l", "3"], " ", SubscriptBox["l", "4"]}], "-", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{ SubscriptBox["x", "2"], "[", "t", "]"}]}], "+", RowBox[{ SubscriptBox["x", "3"], "[", "t", "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{ SubscriptBox["x", "3"], "[", "t", "]"}]}], "+", RowBox[{ SubscriptBox["x", "4"], "[", "t", "]"}]}], ")"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}]}], "+", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}]}], "+", RowBox[{ SubscriptBox["y", "4"], "[", "t", "]"}]}], ")"}]}]}]}, { RowBox[{ RowBox[{ RowBox[{"Sin", "[", RowBox[{ SubscriptBox["\[Phi]", "4"], "[", "t", "]"}], "]"}], " ", SubscriptBox["l", "3"], " ", SubscriptBox["l", "4"]}], "-", RowBox[{ RowBox[{ SubscriptBox["x", "3"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}]}], "+", RowBox[{ RowBox[{ SubscriptBox["x", "4"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}]}], "+", RowBox[{ RowBox[{ SubscriptBox["x", "2"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}]}], "-", RowBox[{ RowBox[{ SubscriptBox["x", "4"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}]}], "-", RowBox[{ RowBox[{ SubscriptBox["x", "2"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "4"], "[", "t", "]"}]}], "+", RowBox[{ RowBox[{ SubscriptBox["x", "3"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "4"], "[", "t", "]"}]}]}]} }, GridBoxAlignment->{ "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.5599999999999999]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], Column], Function[BoxForm`e$, TableForm[BoxForm`e$]]]], "Output", CellChangeTimes->{{3.7389582034011583`*^9, 3.7389582101668053`*^9}, { 3.73895829112352*^9, 3.7389583275862513`*^9}, 3.738958900185498*^9, { 3.738958969477621*^9, 3.738958996939023*^9}, 3.7389590545384197`*^9, { 3.738959182986147*^9, 3.738959228544485*^9}, 3.7389593776810217`*^9, { 3.738959423764462*^9, 3.738959439755392*^9}, {3.738959529005062*^9, 3.738959558555397*^9}, 3.738959674381528*^9, 3.738961077906*^9, 3.738961327676906*^9, 3.772412335816132*^9, 3.772412401275339*^9, 3.7724124596135263`*^9, 3.773660947000496*^9, 3.77366099616361*^9}, CellLabel-> "Out[11]//TableForm=",ExpressionUUID->"db5deadd-6775-4bd5-9e17-\ 182852fbd247"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"0", RowBox[{"-", "1"}], RowBox[{ RowBox[{"-", RowBox[{"Sin", "[", RowBox[{ SubscriptBox["\[Phi]", "1"], "[", "t", "]"}], "]"}]}], " ", SubscriptBox["l", "1"]}], "0", "0", "0", "0", "0", "0", "0", "0", "0"}, { RowBox[{"-", "1"}], "0", RowBox[{ RowBox[{"Cos", "[", RowBox[{ SubscriptBox["\[Phi]", "1"], "[", "t", "]"}], "]"}], " ", SubscriptBox["l", "1"]}], "0", "0", "0", "0", "0", "0", "0", "0", "0"}, { RowBox[{ RowBox[{"2", " ", RowBox[{ SubscriptBox["x", "1"], "[", "t", "]"}]}], "-", RowBox[{ SubscriptBox["x", "2"], "[", "t", "]"}]}], RowBox[{ RowBox[{"2", " ", RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}]}], "-", RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}]}], "0", RowBox[{"-", RowBox[{ SubscriptBox["x", "1"], "[", "t", "]"}]}], RowBox[{"-", RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}]}], RowBox[{ RowBox[{"-", RowBox[{"Sin", "[", RowBox[{ SubscriptBox["\[Phi]", "2"], "[", "t", "]"}], "]"}]}], " ", SubscriptBox["l", "1"], " ", SubscriptBox["l", "2"]}], "0", "0", "0", "0", "0", "0"}, { RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}], RowBox[{"-", RowBox[{ SubscriptBox["x", "2"], "[", "t", "]"}]}], "0", RowBox[{"-", RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}]}], RowBox[{ SubscriptBox["x", "1"], "[", "t", "]"}], RowBox[{ RowBox[{"Cos", "[", RowBox[{ SubscriptBox["\[Phi]", "2"], "[", "t", "]"}], "]"}], " ", SubscriptBox["l", "1"], " ", SubscriptBox["l", "2"]}], "0", "0", "0", "0", "0", "0"}, { RowBox[{ RowBox[{"-", RowBox[{ SubscriptBox["x", "2"], "[", "t", "]"}]}], "+", RowBox[{ SubscriptBox["x", "3"], "[", "t", "]"}]}], RowBox[{ RowBox[{"-", RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}]}], "+", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}]}], "0", RowBox[{ RowBox[{"-", RowBox[{ SubscriptBox["x", "1"], "[", "t", "]"}]}], "+", RowBox[{"2", " ", RowBox[{ SubscriptBox["x", "2"], "[", "t", "]"}]}], "-", RowBox[{ SubscriptBox["x", "3"], "[", "t", "]"}]}], RowBox[{ RowBox[{"-", RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}]}], "+", RowBox[{"2", " ", RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}]}], "-", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}]}], "0", RowBox[{ RowBox[{ SubscriptBox["x", "1"], "[", "t", "]"}], "-", RowBox[{ SubscriptBox["x", "2"], "[", "t", "]"}]}], RowBox[{ RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}], "-", RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}]}], RowBox[{ RowBox[{"-", RowBox[{"Sin", "[", RowBox[{ SubscriptBox["\[Phi]", "3"], "[", "t", "]"}], "]"}]}], " ", SubscriptBox["l", "2"], " ", SubscriptBox["l", "3"]}], "0", "0", "0"}, { RowBox[{ RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}], "-", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}]}], RowBox[{ RowBox[{"-", RowBox[{ SubscriptBox["x", "2"], "[", "t", "]"}]}], "+", RowBox[{ SubscriptBox["x", "3"], "[", "t", "]"}]}], "0", RowBox[{ RowBox[{"-", RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}]}], "+", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}]}], RowBox[{ RowBox[{ SubscriptBox["x", "1"], "[", "t", "]"}], "-", RowBox[{ SubscriptBox["x", "3"], "[", "t", "]"}]}], "0", RowBox[{ RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}], "-", RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}]}], RowBox[{ RowBox[{"-", RowBox[{ SubscriptBox["x", "1"], "[", "t", "]"}]}], "+", RowBox[{ SubscriptBox["x", "2"], "[", "t", "]"}]}], RowBox[{ RowBox[{"Cos", "[", RowBox[{ SubscriptBox["\[Phi]", "3"], "[", "t", "]"}], "]"}], " ", SubscriptBox["l", "2"], " ", SubscriptBox["l", "3"]}], "0", "0", "0"}, {"0", "0", "0", RowBox[{ RowBox[{"-", RowBox[{ SubscriptBox["x", "3"], "[", "t", "]"}]}], "+", RowBox[{ SubscriptBox["x", "4"], "[", "t", "]"}]}], RowBox[{ RowBox[{"-", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}]}], "+", RowBox[{ SubscriptBox["y", "4"], "[", "t", "]"}]}], "0", RowBox[{ RowBox[{"-", RowBox[{ SubscriptBox["x", "2"], "[", "t", "]"}]}], "+", RowBox[{"2", " ", RowBox[{ SubscriptBox["x", "3"], "[", "t", "]"}]}], "-", RowBox[{ SubscriptBox["x", "4"], "[", "t", "]"}]}], RowBox[{ RowBox[{"-", RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}]}], "+", RowBox[{"2", " ", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}]}], "-", RowBox[{ SubscriptBox["y", "4"], "[", "t", "]"}]}], "0", RowBox[{ RowBox[{ SubscriptBox["x", "2"], "[", "t", "]"}], "-", RowBox[{ SubscriptBox["x", "3"], "[", "t", "]"}]}], RowBox[{ RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}], "-", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}]}], RowBox[{ RowBox[{"-", RowBox[{"Sin", "[", RowBox[{ SubscriptBox["\[Phi]", "4"], "[", "t", "]"}], "]"}]}], " ", SubscriptBox["l", "3"], " ", SubscriptBox["l", "4"]}]}, {"0", "0", "0", RowBox[{ RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}], "-", RowBox[{ SubscriptBox["y", "4"], "[", "t", "]"}]}], RowBox[{ RowBox[{"-", RowBox[{ SubscriptBox["x", "3"], "[", "t", "]"}]}], "+", RowBox[{ SubscriptBox["x", "4"], "[", "t", "]"}]}], "0", RowBox[{ RowBox[{"-", RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}]}], "+", RowBox[{ SubscriptBox["y", "4"], "[", "t", "]"}]}], RowBox[{ RowBox[{ SubscriptBox["x", "2"], "[", "t", "]"}], "-", RowBox[{ SubscriptBox["x", "4"], "[", "t", "]"}]}], "0", RowBox[{ RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}], "-", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}]}], RowBox[{ RowBox[{"-", RowBox[{ SubscriptBox["x", "2"], "[", "t", "]"}]}], "+", RowBox[{ SubscriptBox["x", "3"], "[", "t", "]"}]}], RowBox[{ RowBox[{"Cos", "[", RowBox[{ SubscriptBox["\[Phi]", "4"], "[", "t", "]"}], "]"}], " ", SubscriptBox["l", "3"], " ", SubscriptBox["l", "4"]}]} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output", CellChangeTimes->{{3.7389582034011583`*^9, 3.7389582101668053`*^9}, { 3.73895829112352*^9, 3.7389583275862513`*^9}, 3.738958900185498*^9, { 3.738958969477621*^9, 3.738958996939023*^9}, 3.7389590545384197`*^9, { 3.738959182986147*^9, 3.738959228544485*^9}, 3.7389593776810217`*^9, { 3.738959423764462*^9, 3.738959439755392*^9}, {3.738959529005062*^9, 3.738959558555397*^9}, 3.738959674381528*^9, 3.738961077906*^9, 3.738961327676906*^9, 3.772412335816132*^9, 3.772412401275339*^9, 3.7724124596135263`*^9, 3.773660947000496*^9, 3.7736609962073402`*^9}, CellLabel-> "Out[13]//MatrixForm=",ExpressionUUID->"a7a6a88e-f116-4a57-802b-\ 66ab5694c9db"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", TagBox[GridBox[{ { RowBox[{ RowBox[{"Cos", "[", RowBox[{ SubscriptBox["\[Phi]", "1"], "[", "t", "]"}], "]"}], " ", SubscriptBox["l", "1"], " ", SuperscriptBox[ RowBox[{ SuperscriptBox[ SubscriptBox["\[Phi]", "1"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "2"]}]}, { RowBox[{ RowBox[{"Sin", "[", RowBox[{ SubscriptBox["\[Phi]", "1"], "[", "t", "]"}], "]"}], " ", SubscriptBox["l", "1"], " ", SuperscriptBox[ RowBox[{ SuperscriptBox[ SubscriptBox["\[Phi]", "1"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "2"]}]}, { RowBox[{ RowBox[{"2", " ", RowBox[{ SuperscriptBox[ SubscriptBox["x", "1"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{ SuperscriptBox[ SubscriptBox["x", "1"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "+", RowBox[{ SuperscriptBox[ SubscriptBox["x", "2"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], ")"}]}], "+", RowBox[{"2", " ", RowBox[{ SuperscriptBox[ SubscriptBox["y", "1"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{ SuperscriptBox[ SubscriptBox["y", "1"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "+", RowBox[{ SuperscriptBox[ SubscriptBox["y", "2"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], ")"}]}], "+", RowBox[{ RowBox[{"Cos", "[", RowBox[{ SubscriptBox["\[Phi]", "2"], "[", "t", "]"}], "]"}], " ", SubscriptBox["l", "1"], " ", SubscriptBox["l", "2"], " ", SuperscriptBox[ RowBox[{ SuperscriptBox[ SubscriptBox["\[Phi]", "2"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "2"]}]}]}, { RowBox[{ RowBox[{"2", " ", RowBox[{ SuperscriptBox[ SubscriptBox["x", "2"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], " ", RowBox[{ SuperscriptBox[ SubscriptBox["y", "1"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "-", RowBox[{"2", " ", RowBox[{ SuperscriptBox[ SubscriptBox["x", "1"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], " ", RowBox[{ SuperscriptBox[ SubscriptBox["y", "2"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "+", RowBox[{ RowBox[{"Sin", "[", RowBox[{ SubscriptBox["\[Phi]", "2"], "[", "t", "]"}], "]"}], " ", SubscriptBox["l", "1"], " ", SubscriptBox["l", "2"], " ", SuperscriptBox[ RowBox[{ SuperscriptBox[ SubscriptBox["\[Phi]", "2"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "2"]}]}]}, { RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox[ SubscriptBox["x", "1"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "-", RowBox[{ SuperscriptBox[ SubscriptBox["x", "2"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{ SuperscriptBox[ SubscriptBox["x", "2"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "+", RowBox[{ SuperscriptBox[ SubscriptBox["x", "3"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], ")"}]}], "-", RowBox[{"2", " ", RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox[ SubscriptBox["y", "1"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "-", RowBox[{ SuperscriptBox[ SubscriptBox["y", "2"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{ SuperscriptBox[ SubscriptBox["y", "2"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "+", RowBox[{ SuperscriptBox[ SubscriptBox["y", "3"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], ")"}]}], "+", RowBox[{ RowBox[{"Cos", "[", RowBox[{ SubscriptBox["\[Phi]", "3"], "[", "t", "]"}], "]"}], " ", SubscriptBox["l", "2"], " ", SubscriptBox["l", "3"], " ", SuperscriptBox[ RowBox[{ SuperscriptBox[ SubscriptBox["\[Phi]", "3"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "2"]}]}]}, { RowBox[{ RowBox[{"2", " ", RowBox[{ SuperscriptBox[ SubscriptBox["x", "2"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], " ", RowBox[{ SuperscriptBox[ SubscriptBox["y", "1"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "-", RowBox[{"2", " ", RowBox[{ SuperscriptBox[ SubscriptBox["x", "3"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], " ", RowBox[{ SuperscriptBox[ SubscriptBox["y", "1"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "-", RowBox[{"2", " ", RowBox[{ SuperscriptBox[ SubscriptBox["x", "1"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], " ", RowBox[{ SuperscriptBox[ SubscriptBox["y", "2"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "+", RowBox[{"2", " ", RowBox[{ SuperscriptBox[ SubscriptBox["x", "3"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], " ", RowBox[{ SuperscriptBox[ SubscriptBox["y", "2"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "+", RowBox[{"2", " ", RowBox[{ SuperscriptBox[ SubscriptBox["x", "1"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], " ", RowBox[{ SuperscriptBox[ SubscriptBox["y", "3"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "-", RowBox[{"2", " ", RowBox[{ SuperscriptBox[ SubscriptBox["x", "2"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], " ", RowBox[{ SuperscriptBox[ SubscriptBox["y", "3"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "+", RowBox[{ RowBox[{"Sin", "[", RowBox[{ SubscriptBox["\[Phi]", "3"], "[", "t", "]"}], "]"}], " ", SubscriptBox["l", "2"], " ", SubscriptBox["l", "3"], " ", SuperscriptBox[ RowBox[{ SuperscriptBox[ SubscriptBox["\[Phi]", "3"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "2"]}]}]}, { RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox[ SubscriptBox["x", "2"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "-", RowBox[{ SuperscriptBox[ SubscriptBox["x", "3"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{ SuperscriptBox[ SubscriptBox["x", "3"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "+", RowBox[{ SuperscriptBox[ SubscriptBox["x", "4"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], ")"}]}], "-", RowBox[{"2", " ", RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox[ SubscriptBox["y", "2"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "-", RowBox[{ SuperscriptBox[ SubscriptBox["y", "3"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{ SuperscriptBox[ SubscriptBox["y", "3"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "+", RowBox[{ SuperscriptBox[ SubscriptBox["y", "4"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], ")"}]}], "+", RowBox[{ RowBox[{"Cos", "[", RowBox[{ SubscriptBox["\[Phi]", "4"], "[", "t", "]"}], "]"}], " ", SubscriptBox["l", "3"], " ", SubscriptBox["l", "4"], " ", SuperscriptBox[ RowBox[{ SuperscriptBox[ SubscriptBox["\[Phi]", "4"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "2"]}]}]}, { RowBox[{ RowBox[{"2", " ", RowBox[{ SuperscriptBox[ SubscriptBox["x", "3"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], " ", RowBox[{ SuperscriptBox[ SubscriptBox["y", "2"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "-", RowBox[{"2", " ", RowBox[{ SuperscriptBox[ SubscriptBox["x", "4"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], " ", RowBox[{ SuperscriptBox[ SubscriptBox["y", "2"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "-", RowBox[{"2", " ", RowBox[{ SuperscriptBox[ SubscriptBox["x", "2"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], " ", RowBox[{ SuperscriptBox[ SubscriptBox["y", "3"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "+", RowBox[{"2", " ", RowBox[{ SuperscriptBox[ SubscriptBox["x", "4"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], " ", RowBox[{ SuperscriptBox[ SubscriptBox["y", "3"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "+", RowBox[{"2", " ", RowBox[{ SuperscriptBox[ SubscriptBox["x", "2"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], " ", RowBox[{ SuperscriptBox[ SubscriptBox["y", "4"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "-", RowBox[{"2", " ", RowBox[{ SuperscriptBox[ SubscriptBox["x", "3"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], " ", RowBox[{ SuperscriptBox[ SubscriptBox["y", "4"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "+", RowBox[{ RowBox[{"Sin", "[", RowBox[{ SubscriptBox["\[Phi]", "4"], "[", "t", "]"}], "]"}], " ", SubscriptBox["l", "3"], " ", SubscriptBox["l", "4"], " ", SuperscriptBox[ RowBox[{ SuperscriptBox[ SubscriptBox["\[Phi]", "4"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "2"]}]}]} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.5599999999999999]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], Column], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output", CellChangeTimes->{{3.7389582034011583`*^9, 3.7389582101668053`*^9}, { 3.73895829112352*^9, 3.7389583275862513`*^9}, 3.738958900185498*^9, { 3.738958969477621*^9, 3.738958996939023*^9}, 3.7389590545384197`*^9, { 3.738959182986147*^9, 3.738959228544485*^9}, 3.7389593776810217`*^9, { 3.738959423764462*^9, 3.738959439755392*^9}, {3.738959529005062*^9, 3.738959558555397*^9}, 3.738959674381528*^9, 3.738961077906*^9, 3.738961327676906*^9, 3.772412335816132*^9, 3.772412401275339*^9, 3.7724124596135263`*^9, 3.773660947000496*^9, 3.7736609962117243`*^9}, CellLabel-> "Out[14]//MatrixForm=",ExpressionUUID->"e7c8faff-5277-42bc-821b-\ bd2c9dd42172"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{ RowBox[{ SubscriptBox["m", "1"], " ", RowBox[{"(", RowBox[{ FractionBox[ SubsuperscriptBox["\[Rho]", "1", "2"], SubsuperscriptBox["l", "1", "2"]], "+", SubsuperscriptBox["\[Chi]", "1", "2"]}], ")"}]}], "+", RowBox[{ SubscriptBox["m", "2"], " ", RowBox[{"(", RowBox[{ FractionBox[ SubsuperscriptBox["\[Rho]", "2", "2"], SubsuperscriptBox["l", "2", "2"]], "+", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", SubscriptBox["\[Chi]", "2"]}], ")"}], "2"]}], ")"}]}]}], "0", "0", RowBox[{ SubscriptBox["m", "2"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", FractionBox[ SubsuperscriptBox["\[Rho]", "2", "2"], SubsuperscriptBox["l", "2", "2"]]}], "+", SubscriptBox["\[Chi]", "2"], "-", SubsuperscriptBox["\[Chi]", "2", "2"]}], ")"}]}], "0", "0", "0", "0", "0", "0", "0", "0"}, {"0", RowBox[{ RowBox[{ SubscriptBox["m", "1"], " ", RowBox[{"(", RowBox[{ FractionBox[ SubsuperscriptBox["\[Rho]", "1", "2"], SubsuperscriptBox["l", "1", "2"]], "+", SubsuperscriptBox["\[Chi]", "1", "2"]}], ")"}]}], "+", RowBox[{ SubscriptBox["m", "2"], " ", RowBox[{"(", RowBox[{ FractionBox[ SubsuperscriptBox["\[Rho]", "2", "2"], SubsuperscriptBox["l", "2", "2"]], "+", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", SubscriptBox["\[Chi]", "2"]}], ")"}], "2"]}], ")"}]}]}], "0", "0", RowBox[{ SubscriptBox["m", "2"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", FractionBox[ SubsuperscriptBox["\[Rho]", "2", "2"], SubsuperscriptBox["l", "2", "2"]]}], "+", SubscriptBox["\[Chi]", "2"], "-", SubsuperscriptBox["\[Chi]", "2", "2"]}], ")"}]}], "0", "0", "0", "0", "0", "0", "0"}, {"0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0"}, { RowBox[{ SubscriptBox["m", "2"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", FractionBox[ SubsuperscriptBox["\[Rho]", "2", "2"], SubsuperscriptBox["l", "2", "2"]]}], "+", SubscriptBox["\[Chi]", "2"], "-", SubsuperscriptBox["\[Chi]", "2", "2"]}], ")"}]}], "0", "0", RowBox[{ RowBox[{ SubscriptBox["m", "2"], " ", RowBox[{"(", RowBox[{ FractionBox[ SubsuperscriptBox["\[Rho]", "2", "2"], SubsuperscriptBox["l", "2", "2"]], "+", SubsuperscriptBox["\[Chi]", "2", "2"]}], ")"}]}], "+", RowBox[{ SubscriptBox["m", "3"], " ", RowBox[{"(", RowBox[{ FractionBox[ SubsuperscriptBox["\[Rho]", "3", "2"], SubsuperscriptBox["l", "3", "2"]], "+", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", SubscriptBox["\[Chi]", "3"]}], ")"}], "2"]}], ")"}]}]}], "0", "0", RowBox[{ SubscriptBox["m", "3"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", FractionBox[ SubsuperscriptBox["\[Rho]", "3", "2"], SubsuperscriptBox["l", "3", "2"]]}], "+", SubscriptBox["\[Chi]", "3"], "-", SubsuperscriptBox["\[Chi]", "3", "2"]}], ")"}]}], "0", "0", "0", "0", "0"}, {"0", RowBox[{ SubscriptBox["m", "2"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", FractionBox[ SubsuperscriptBox["\[Rho]", "2", "2"], SubsuperscriptBox["l", "2", "2"]]}], "+", SubscriptBox["\[Chi]", "2"], "-", SubsuperscriptBox["\[Chi]", "2", "2"]}], ")"}]}], "0", "0", RowBox[{ RowBox[{ SubscriptBox["m", "2"], " ", RowBox[{"(", RowBox[{ FractionBox[ SubsuperscriptBox["\[Rho]", "2", "2"], SubsuperscriptBox["l", "2", "2"]], "+", SubsuperscriptBox["\[Chi]", "2", "2"]}], ")"}]}], "+", RowBox[{ SubscriptBox["m", "3"], " ", RowBox[{"(", RowBox[{ FractionBox[ SubsuperscriptBox["\[Rho]", "3", "2"], SubsuperscriptBox["l", "3", "2"]], "+", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", SubscriptBox["\[Chi]", "3"]}], ")"}], "2"]}], ")"}]}]}], "0", "0", RowBox[{ SubscriptBox["m", "3"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", FractionBox[ SubsuperscriptBox["\[Rho]", "3", "2"], SubsuperscriptBox["l", "3", "2"]]}], "+", SubscriptBox["\[Chi]", "3"], "-", SubsuperscriptBox["\[Chi]", "3", "2"]}], ")"}]}], "0", "0", "0", "0"}, {"0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0"}, {"0", "0", "0", RowBox[{ SubscriptBox["m", "3"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", FractionBox[ SubsuperscriptBox["\[Rho]", "3", "2"], SubsuperscriptBox["l", "3", "2"]]}], "+", SubscriptBox["\[Chi]", "3"], "-", SubsuperscriptBox["\[Chi]", "3", "2"]}], ")"}]}], "0", "0", RowBox[{ RowBox[{ SubscriptBox["m", "3"], " ", RowBox[{"(", RowBox[{ FractionBox[ SubsuperscriptBox["\[Rho]", "3", "2"], SubsuperscriptBox["l", "3", "2"]], "+", SubsuperscriptBox["\[Chi]", "3", "2"]}], ")"}]}], "+", RowBox[{ SubscriptBox["m", "4"], " ", RowBox[{"(", RowBox[{ FractionBox[ SubsuperscriptBox["\[Rho]", "4", "2"], SubsuperscriptBox["l", "4", "2"]], "+", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", SubscriptBox["\[Chi]", "4"]}], ")"}], "2"]}], ")"}]}]}], "0", "0", RowBox[{ SubscriptBox["m", "4"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", FractionBox[ SubsuperscriptBox["\[Rho]", "4", "2"], SubsuperscriptBox["l", "4", "2"]]}], "+", SubscriptBox["\[Chi]", "4"], "-", SubsuperscriptBox["\[Chi]", "4", "2"]}], ")"}]}], "0", "0"}, {"0", "0", "0", "0", RowBox[{ SubscriptBox["m", "3"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", FractionBox[ SubsuperscriptBox["\[Rho]", "3", "2"], SubsuperscriptBox["l", "3", "2"]]}], "+", SubscriptBox["\[Chi]", "3"], "-", SubsuperscriptBox["\[Chi]", "3", "2"]}], ")"}]}], "0", "0", RowBox[{ RowBox[{ SubscriptBox["m", "3"], " ", RowBox[{"(", RowBox[{ FractionBox[ SubsuperscriptBox["\[Rho]", "3", "2"], SubsuperscriptBox["l", "3", "2"]], "+", SubsuperscriptBox["\[Chi]", "3", "2"]}], ")"}]}], "+", RowBox[{ SubscriptBox["m", "4"], " ", RowBox[{"(", RowBox[{ FractionBox[ SubsuperscriptBox["\[Rho]", "4", "2"], SubsuperscriptBox["l", "4", "2"]], "+", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", SubscriptBox["\[Chi]", "4"]}], ")"}], "2"]}], ")"}]}]}], "0", "0", RowBox[{ SubscriptBox["m", "4"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", FractionBox[ SubsuperscriptBox["\[Rho]", "4", "2"], SubsuperscriptBox["l", "4", "2"]]}], "+", SubscriptBox["\[Chi]", "4"], "-", SubsuperscriptBox["\[Chi]", "4", "2"]}], ")"}]}], "0"}, {"0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0"}, {"0", "0", "0", "0", "0", "0", RowBox[{ SubscriptBox["m", "4"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", FractionBox[ SubsuperscriptBox["\[Rho]", "4", "2"], SubsuperscriptBox["l", "4", "2"]]}], "+", SubscriptBox["\[Chi]", "4"], "-", SubsuperscriptBox["\[Chi]", "4", "2"]}], ")"}]}], "0", "0", RowBox[{ SubscriptBox["m", "4"], " ", RowBox[{"(", RowBox[{ FractionBox[ SubsuperscriptBox["\[Rho]", "4", "2"], SubsuperscriptBox["l", "4", "2"]], "+", SubsuperscriptBox["\[Chi]", "4", "2"]}], ")"}]}], "0", "0"}, {"0", "0", "0", "0", "0", "0", "0", RowBox[{ SubscriptBox["m", "4"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", FractionBox[ SubsuperscriptBox["\[Rho]", "4", "2"], SubsuperscriptBox["l", "4", "2"]]}], "+", SubscriptBox["\[Chi]", "4"], "-", SubsuperscriptBox["\[Chi]", "4", "2"]}], ")"}]}], "0", "0", RowBox[{ SubscriptBox["m", "4"], " ", RowBox[{"(", RowBox[{ FractionBox[ SubsuperscriptBox["\[Rho]", "4", "2"], SubsuperscriptBox["l", "4", "2"]], "+", SubsuperscriptBox["\[Chi]", "4", "2"]}], ")"}]}], "0"}, {"0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0"} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output", CellChangeTimes->{{3.7389582034011583`*^9, 3.7389582101668053`*^9}, { 3.73895829112352*^9, 3.7389583275862513`*^9}, 3.738958900185498*^9, { 3.738958969477621*^9, 3.738958996939023*^9}, 3.7389590545384197`*^9, { 3.738959182986147*^9, 3.738959228544485*^9}, 3.7389593776810217`*^9, { 3.738959423764462*^9, 3.738959439755392*^9}, {3.738959529005062*^9, 3.738959558555397*^9}, 3.738959674381528*^9, 3.738961077906*^9, 3.738961327676906*^9, 3.772412335816132*^9, 3.772412401275339*^9, 3.7724124596135263`*^9, 3.773660947000496*^9, 3.773661000650528*^9}, CellLabel-> "Out[18]//MatrixForm=",ExpressionUUID->"da68d85e-4e97-45b7-873f-\ 6800e69abb81"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", TagBox[GridBox[{ {"0"}, { RowBox[{"g", " ", RowBox[{"(", RowBox[{ RowBox[{ SubscriptBox["m", "1"], " ", SubscriptBox["\[Chi]", "1"]}], "-", RowBox[{ SubscriptBox["m", "2"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", SubscriptBox["\[Chi]", "2"]}], ")"}]}]}], ")"}]}]}, {"0"}, {"0"}, { RowBox[{"g", " ", RowBox[{"(", RowBox[{ RowBox[{ SubscriptBox["m", "2"], " ", SubscriptBox["\[Chi]", "2"]}], "-", RowBox[{ SubscriptBox["m", "3"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", SubscriptBox["\[Chi]", "3"]}], ")"}]}]}], ")"}]}]}, {"0"}, {"0"}, { RowBox[{"g", " ", RowBox[{"(", RowBox[{ RowBox[{ SubscriptBox["m", "3"], " ", SubscriptBox["\[Chi]", "3"]}], "-", RowBox[{ SubscriptBox["m", "4"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", SubscriptBox["\[Chi]", "4"]}], ")"}]}]}], ")"}]}]}, {"0"}, {"0"}, { RowBox[{"g", " ", SubscriptBox["m", "4"], " ", SubscriptBox["\[Chi]", "4"]}]}, {"0"} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.5599999999999999]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], Column], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output", CellChangeTimes->{{3.7389582034011583`*^9, 3.7389582101668053`*^9}, { 3.73895829112352*^9, 3.7389583275862513`*^9}, 3.738958900185498*^9, { 3.738958969477621*^9, 3.738958996939023*^9}, 3.7389590545384197`*^9, { 3.738959182986147*^9, 3.738959228544485*^9}, 3.7389593776810217`*^9, { 3.738959423764462*^9, 3.738959439755392*^9}, {3.738959529005062*^9, 3.738959558555397*^9}, 3.738959674381528*^9, 3.738961077906*^9, 3.738961327676906*^9, 3.772412335816132*^9, 3.772412401275339*^9, 3.7724124596135263`*^9, 3.773660947000496*^9, 3.773661000651754*^9}, CellLabel-> "Out[19]//MatrixForm=",ExpressionUUID->"85f255a8-1d60-4a06-bf75-\ c26e32ecca97"] }, Open ]], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"\[DoubleStruckM]", "[", "j_", "]"}], "=", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ SubscriptBox["m", "j"], " ", RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", FractionBox[ SubscriptBox["\[Rho]", "j"], SubscriptBox["l", "j"]], ")"}], "2"], "+", SubsuperscriptBox["\[Chi]", "j", "2"]}], ")"}]}], "+", RowBox[{ SubscriptBox["m", RowBox[{"j", "+", "1"}]], " ", RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", FractionBox[ SubscriptBox["\[Rho]", RowBox[{"j", "+", "1"}]], SubscriptBox["l", RowBox[{"j", "+", "1"}]]], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{ SubscriptBox["\[Chi]", RowBox[{"j", "+", "1"}]], "-", "1"}], ")"}], "2"]}], ")"}]}]}], ")"}], RowBox[{"(", GridBox[{ {"1", "0", "0"}, {"0", "1", "0"}, {"0", "0", "0"} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], ")"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"\[DoubleStruckM]\[DoubleStruckN]", "[", "j_", "]"}], "=", RowBox[{ SubscriptBox["m", "j"], " ", RowBox[{"(", RowBox[{ RowBox[{ SubscriptBox["\[Chi]", "j"], RowBox[{"(", RowBox[{"1", "-", SubscriptBox["\[Chi]", "j"]}], ")"}]}], "-", SuperscriptBox[ RowBox[{"(", FractionBox[ SubscriptBox["\[Rho]", "j"], SubscriptBox["l", "j"]], ")"}], "2"]}], ")"}], RowBox[{"(", GridBox[{ {"1", "0", "0"}, {"0", "1", "0"}, {"0", "0", "0"} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], ")"}]}]}], ";"}]}], "Code", CellChangeTimes->{{3.738960490961378*^9, 3.738960723343835*^9}, { 3.7389608776653547`*^9, 3.738960888944178*^9}, {3.738963104849278*^9, 3.738963106767695*^9}}, CellLabel->"In[20]:=",ExpressionUUID->"46e414bf-6509-4aa9-af6a-1e9a43438b46"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"ArrayFlatten", "[", RowBox[{"(", GridBox[{ { RowBox[{"\[DoubleStruckM]", "[", "1", "]"}], RowBox[{"\[DoubleStruckM]\[DoubleStruckN]", "[", "2", "]"}], "0", "0"}, { RowBox[{"\[DoubleStruckM]\[DoubleStruckN]", "[", "2", "]"}], RowBox[{"\[DoubleStruckM]", "[", "2", "]"}], RowBox[{"\[DoubleStruckM]\[DoubleStruckN]", "[", "3", "]"}], "0"}, {"0", RowBox[{"\[DoubleStruckM]\[DoubleStruckN]", "[", "3", "]"}], RowBox[{"\[DoubleStruckM]", "[", "3", "]"}], RowBox[{"\[DoubleStruckM]\[DoubleStruckN]", "[", "4", "]"}]}, {"0", "0", RowBox[{"\[DoubleStruckM]\[DoubleStruckN]", "[", "4", "]"}], RowBox[{"\[DoubleStruckM]", "[", "4", "]"}]} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], ")"}], "]"}], "/.", RowBox[{ SubscriptBox["m", "5"], "\[Rule]", "0"}]}], ")"}], "-", RowBox[{"\[DoubleStruckCapitalM]", "[", "t", "]"}]}], "//", "Simplify"}], "//", "Norm"}]], "Code", CellChangeTimes->{{3.738960730287775*^9, 3.738960868899097*^9}, { 3.773661005342558*^9, 3.773661006262004*^9}}, CellLabel->"In[23]:=",ExpressionUUID->"da989899-3633-4276-98ff-a9971236ae0b"], Cell[BoxData["0"], "Output", CellChangeTimes->{{3.738960835325244*^9, 3.738960890448263*^9}, 3.738963108484775*^9, 3.7724123650183153`*^9, 3.772412436162724*^9, 3.773660951003384*^9, {3.773661000723195*^9, 3.773661006739204*^9}}, CellLabel->"Out[23]=",ExpressionUUID->"711a6ac8-c78b-4d23-a808-7a0d6e7eaf15"] }, Open ]] }, Open ]] }, WindowSize->{1680, 998}, WindowMargins->{{0, Automatic}, {1, Automatic}}, Magnification:>1.25 Inherited, FrontEndVersion->"11.3 for Linux x86 (64-bit) (March 6, 2018)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[580, 22, 336, 5, 84, "Section",ExpressionUUID->"ab0118e6-f1e9-40e8-94a5-172a16b58603"], Cell[919, 29, 4533, 134, 160, "Code",ExpressionUUID->"8a0291bb-a54c-4f63-97d8-90f5c57b242d"] }, Open ]], Cell[CellGroupData[{ Cell[5489, 168, 249, 4, 84, "Section",ExpressionUUID->"b23ed7d1-72d4-433a-8ada-b443518b34b0"], Cell[5741, 174, 10707, 313, 504, "Code",ExpressionUUID->"95515b98-4828-43a1-a754-b9ffaafce4ae"], Cell[16451, 489, 312, 7, 39, "Input",ExpressionUUID->"9aaccd05-a673-44f1-8df0-6ba7b7a41952"], Cell[16766, 498, 303, 7, 39, "Input",ExpressionUUID->"9c5f0a77-452d-411b-a09d-0e94b018cce1"], Cell[17072, 507, 445, 10, 39, "Input",ExpressionUUID->"a53523e0-5ac1-462d-8152-e417939d03b8"] }, Open ]], Cell[CellGroupData[{ Cell[17554, 522, 217, 3, 84, "Section",ExpressionUUID->"89ecdc42-0a7d-421c-9292-790d76cf9a3f"], Cell[CellGroupData[{ Cell[17796, 529, 11534, 325, 747, "Code",ExpressionUUID->"2c63b02f-3cfe-40f5-8f2c-f810b5ce01ec"], Cell[29333, 856, 7960, 243, 228, "Output",ExpressionUUID->"db5deadd-6775-4bd5-9e17-182852fbd247"], Cell[37296, 1101, 8537, 262, 228, "Output",ExpressionUUID->"a7a6a88e-f116-4a57-802b-66ab5694c9db"], Cell[45836, 1365, 13599, 376, 260, "Output",ExpressionUUID->"e7c8faff-5277-42bc-821b-bd2c9dd42172"], Cell[59438, 1743, 10892, 309, 572, "Output",ExpressionUUID->"da68d85e-4e97-45b7-873f-6800e69abb81"], Cell[70333, 2054, 2663, 79, 324, "Output",ExpressionUUID->"85f255a8-1d60-4a06-bf75-c26e32ecca97"] }, Open ]], Cell[73011, 2136, 2869, 86, 190, "Code",ExpressionUUID->"46e414bf-6509-4aa9-af6a-1e9a43438b46"], Cell[CellGroupData[{ Cell[75905, 2226, 1760, 40, 126, "Code",ExpressionUUID->"da989899-3633-4276-98ff-a9971236ae0b"], Cell[77668, 2268, 317, 4, 44, "Output",ExpressionUUID->"711a6ac8-c78b-4d23-a808-7a0d6e7eaf15"] }, Open ]] }, Open ]] } ] *)