(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 8.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 157, 7] NotebookDataLength[ 35197, 824] NotebookOptionsPosition[ 34631, 802] NotebookOutlinePosition[ 34975, 817] CellTagsIndexPosition[ 34932, 814] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"h", "=", "0"}], ";"}], "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"R", "=", "5"}], ",", RowBox[{"R", "\[LessEqual]", "5"}], ",", RowBox[{"R", "=", RowBox[{"R", "+", "1"}]}], ",", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"S", "=", "6"}], ",", RowBox[{"S", "\[LessEqual]", "6"}], ",", RowBox[{"S", "=", RowBox[{"s", "+", "1"}]}], ",", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"a", "=", "2.8"}], ",", RowBox[{"a", "\[LessEqual]", "3"}], ",", RowBox[{"a", "=", RowBox[{"a", "+", "0.1"}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"T", "=", "3000"}], ",", RowBox[{"T", "\[LessEqual]", "4000"}], ",", RowBox[{"T", "=", RowBox[{"T", "+", "100"}]}], ",", "\[IndentingNewLine]", RowBox[{"m", "=", RowBox[{"NDSolve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ SuperscriptBox["y1", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", RowBox[{"y7", "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ SuperscriptBox["y2", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", RowBox[{"y8", "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ SuperscriptBox["y3", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", RowBox[{"y9", "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ RowBox[{"y4", "'"}], "[", "t", "]"}], "==", RowBox[{"y10", "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ RowBox[{"y5", "'"}], "[", "t", "]"}], "\[Equal]", RowBox[{"y11", "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ RowBox[{"y6", "'"}], "[", "t", "]"}], "\[Equal]", RowBox[{"y12", "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ RowBox[{"y7", "'"}], "[", "t", "]"}], "\[Equal]", RowBox[{"y13", "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ RowBox[{"y8", "'"}], "[", "t", "]"}], "\[Equal]", RowBox[{"y14", "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ RowBox[{"y9", "'"}], "[", "t", "]"}], "\[Equal]", RowBox[{"y15", "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ RowBox[{"y10", "'"}], "[", "t", "]"}], "\[Equal]", RowBox[{"y16", "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ RowBox[{"y11", "'"}], "[", "t", "]"}], "\[Equal]", RowBox[{"y17", "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ RowBox[{"y12", "'"}], "[", "t", "]"}], "\[Equal]", RowBox[{"y18", "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ RowBox[{"y13", "'"}], "[", "t", "]"}], "\[Equal]", RowBox[{"y19", "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ RowBox[{"y14", "'"}], "[", "t", "]"}], "\[Equal]", RowBox[{"y20", "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ RowBox[{"y15", "'"}], "[", "t", "]"}], "\[Equal]", RowBox[{"y21", "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ RowBox[{"y16", "'"}], "[", "t", "]"}], "\[Equal]", RowBox[{"y22", "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ RowBox[{"y17", "'"}], "[", "t", "]"}], "\[Equal]", RowBox[{"y23", "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ RowBox[{"y18", "'"}], "[", "t", "]"}], "\[Equal]", RowBox[{"y24", "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ SuperscriptBox["y19", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{"a", "^", "2"}], ")"}], "*", "T", "*", RowBox[{"(", RowBox[{"1", "-", "t"}], ")"}]}], ")"}]}], " ", RowBox[{"y25", "[", "t", "]"}]}], "+", RowBox[{"2", "*", RowBox[{"(", RowBox[{"a", "^", "2"}], ")"}], " ", "*", RowBox[{"y13", "[", "t", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"a", "^", "4"}], ")"}], "*", RowBox[{"y1", "[", "t", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"S", "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"h", "*", RowBox[{"(", RowBox[{"T", "^", "0.5"}], ")"}]}], "+", RowBox[{"a", " ", "*", "R"}]}], ")"}], RowBox[{"(", RowBox[{"1", "-", "t"}], ")"}]}]}], ")"}], RowBox[{"y16", "[", "t", "]"}]}], "+", RowBox[{ RowBox[{"a", "^", "2"}], RowBox[{"(", RowBox[{"S", "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"h", "*", RowBox[{"(", RowBox[{"T", "^", "0.5"}], ")"}]}], "+", RowBox[{"a", " ", "*", "R"}]}], ")"}], RowBox[{"(", RowBox[{"1", "-", "t"}], ")"}]}]}], ")"}], RowBox[{"y4", "[", "t", "]"}]}]}]}], ",", RowBox[{ RowBox[{ SuperscriptBox["y20", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{"a", "^", "2"}], ")"}], "*", "T", "*", RowBox[{"(", RowBox[{"1", "-", "t"}], ")"}]}], ")"}]}], " ", RowBox[{"y26", "[", "t", "]"}]}], "+", RowBox[{"2", "*", RowBox[{"(", RowBox[{"a", "^", "2"}], ")"}], " ", "*", RowBox[{"y14", "[", "t", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"a", "^", "4"}], ")"}], "*", RowBox[{"y2", "[", "t", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"S", "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"h", "*", RowBox[{"(", RowBox[{"T", "^", "0.5"}], ")"}]}], "+", RowBox[{"a", " ", "*", "R"}]}], ")"}], RowBox[{"(", RowBox[{"1", "-", "t"}], ")"}]}]}], ")"}], RowBox[{"y17", "[", "t", "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"a", "^", "2"}], ")"}], "*", RowBox[{"(", RowBox[{"S", "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"h", "*", RowBox[{ RowBox[{"(", "T", ")"}], "^", "0.5"}]}], "+", RowBox[{"a", " ", "*", "R"}]}], ")"}], RowBox[{"(", RowBox[{"1", "-", "t"}], ")"}]}]}], ")"}], RowBox[{"y5", "[", "t", "]"}]}]}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ SuperscriptBox["y21", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{"a", "^", "2"}], ")"}], "*", "T", "*", RowBox[{"(", RowBox[{"1", "-", "t"}], ")"}]}], ")"}]}], " ", RowBox[{"y27", "[", "t", "]"}]}], "+", RowBox[{"2", "*", RowBox[{"(", RowBox[{"a", "^", "2"}], ")"}], " ", "*", RowBox[{"y15", "[", "t", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"a", "^", "4"}], ")"}], "*", RowBox[{"y3", "[", "t", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"S", "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"h", "*", RowBox[{"(", RowBox[{"T", "^", "0.5"}], ")"}]}], "+", RowBox[{"a", " ", "*", "R"}]}], ")"}], RowBox[{"(", RowBox[{"1", "-", "t"}], ")"}]}]}], ")"}], RowBox[{"y18", "[", "t", "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"a", "^", "2"}], ")"}], "*", RowBox[{"(", RowBox[{"S", "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"h", "*", RowBox[{ RowBox[{"(", "T", ")"}], "^", "0.5"}]}], "+", RowBox[{"a", " ", "*", "R"}]}], ")"}], RowBox[{"(", RowBox[{"1", "-", "t"}], ")"}]}]}], ")"}], RowBox[{"y6", "[", "t", "]"}]}]}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ SuperscriptBox["y22", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{"a", "^", "2"}], ")"}], "*", "T", "*", RowBox[{"(", RowBox[{"1", "-", "t"}], ")"}]}], ")"}]}], " ", RowBox[{"y28", "[", "t", "]"}]}], "+", RowBox[{"2", "*", RowBox[{"(", RowBox[{"a", "^", "2"}], " ", ")"}], "*", RowBox[{"y16", "[", "t", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"a", "^", "4"}], ")"}], "*", RowBox[{"y4", "[", "t", "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"S", "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"h", "*", RowBox[{"(", RowBox[{"T", "^", "0.5"}], ")"}]}], "+", RowBox[{"a", " ", "*", "R"}]}], ")"}], RowBox[{"(", RowBox[{"1", "-", "t"}], ")"}]}]}], ")"}], RowBox[{"y13", "[", "t", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"a", "^", "2"}], ")"}], "*", RowBox[{"(", RowBox[{"S", "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"h", "*", RowBox[{ RowBox[{"(", "T", ")"}], "^", "0.5"}]}], "+", RowBox[{"a", " ", "*", "R"}]}], ")"}], RowBox[{"(", RowBox[{"1", "-", "t"}], ")"}]}]}], ")"}], RowBox[{"y1", "[", "t", "]"}]}]}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ SuperscriptBox["y23", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{"a", "^", "2"}], ")"}], "*", "T", "*", RowBox[{"(", RowBox[{"1", "-", "t"}], ")"}]}], ")"}]}], " ", RowBox[{"y29", "[", "t", "]"}]}], "+", RowBox[{"2", "*", RowBox[{"(", RowBox[{"a", "^", "2"}], " ", ")"}], "*", RowBox[{"y17", "[", "t", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"a", "^", "4"}], ")"}], "*", RowBox[{"y5", "[", "t", "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"S", "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"h", "*", RowBox[{"(", RowBox[{"T", "^", "0.5"}], ")"}]}], "+", RowBox[{"a", " ", "*", "R"}]}], ")"}], RowBox[{"(", RowBox[{"1", "-", "t"}], ")"}]}]}], ")"}], RowBox[{"y14", "[", "t", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"a", "^", "2"}], ")"}], "*", RowBox[{"(", RowBox[{"S", "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"h", "*", RowBox[{ RowBox[{"(", "T", ")"}], "^", "0.5"}]}], "+", RowBox[{"a", " ", "*", "R"}]}], ")"}], RowBox[{"(", RowBox[{"1", "-", "t"}], ")"}]}]}], ")"}], RowBox[{"y2", "[", "t", "]"}]}]}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ SuperscriptBox["y24", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{"a", "^", "2"}], ")"}], "*", "T", "*", RowBox[{"(", RowBox[{"1", "-", "t"}], ")"}]}], ")"}]}], " ", RowBox[{"y30", "[", "t", "]"}]}], "+", RowBox[{"2", "*", RowBox[{"(", RowBox[{"a", "^", "2"}], ")"}], " ", "*", RowBox[{"y18", "[", "t", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"a", "^", "4"}], ")"}], "*", RowBox[{"y6", "[", "t", "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"S", "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"h", "*", RowBox[{"(", RowBox[{"T", "^", "0.5"}], ")"}]}], "+", RowBox[{"a", " ", "*", "R"}]}], ")"}], RowBox[{"(", RowBox[{"1", "-", "t"}], ")"}]}]}], ")"}], RowBox[{"y15", "[", "t", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"a", "^", "2"}], ")"}], "*", RowBox[{"(", RowBox[{"S", "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"h", "*", RowBox[{ RowBox[{"(", "T", ")"}], "^", "0.5"}]}], "+", RowBox[{"a", " ", "*", "R"}]}], ")"}], RowBox[{"(", RowBox[{"1", "-", "t"}], ")"}]}]}], ")"}], RowBox[{"y3", "[", "t", "]"}]}]}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"y25", "'"}], "[", "t", "]"}], "\[Equal]", RowBox[{"y31", "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ RowBox[{"y26", "'"}], "[", "t", "]"}], "\[Equal]", RowBox[{"y32", "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ RowBox[{"y27", "'"}], "[", "t", "]"}], "\[Equal]", RowBox[{"y33", "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ RowBox[{"y28", "'"}], "[", "t", "]"}], "\[Equal]", RowBox[{"y34", "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ RowBox[{"y29", "'"}], "[", "t", "]"}], "\[Equal]", RowBox[{"y35", "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ RowBox[{"y30", "'"}], "[", "t", "]"}], "\[Equal]", RowBox[{"y36", "[", "t", "]"}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"y31", "'"}], "[", "t", "]"}], "\[Equal]", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"a", "^", "2"}], ")"}], "*", RowBox[{"y25", "[", "t", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"S", "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"h", "*", RowBox[{ RowBox[{"(", "T", ")"}], "^", "0.5"}]}], "+", RowBox[{"a", " ", "*", "R"}]}], ")"}], RowBox[{"(", RowBox[{"1", "-", "t"}], ")"}]}]}], ")"}], RowBox[{"y28", "[", "t", "]"}]}], "+", RowBox[{"y1", "[", "t", "]"}]}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"y32", "[", "t", "]"}], "\[Equal]", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"a", "^", "2"}], ")"}], "*", RowBox[{"y26", "[", "t", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"S", "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"h", "*", RowBox[{ RowBox[{"(", "T", ")"}], "^", "0.5"}]}], "+", RowBox[{"a", " ", "*", "R"}]}], ")"}], RowBox[{"(", RowBox[{"1", "-", "t"}], ")"}]}]}], ")"}], RowBox[{"y29", "[", "t", "]"}]}], "+", RowBox[{"y2", "[", "t", "]"}]}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"y33", "'"}], "[", "t", "]"}], "\[Equal]", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"a", "^", "2"}], ")"}], "*", RowBox[{"y27", "[", "t", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"S", "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"h", "*", RowBox[{ RowBox[{"(", "T", ")"}], "^", "0.5"}]}], "+", RowBox[{"a", " ", "*", "R"}]}], ")"}], RowBox[{"(", RowBox[{"1", "-", "t"}], ")"}]}]}], ")"}], RowBox[{"y30", "[", "t", "]"}]}], "+", RowBox[{"y3", "[", "t", "]"}]}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"y34", "'"}], "[", "t", "]"}], "\[Equal]", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"a", "^", "2"}], ")"}], "*", RowBox[{"y28", "[", "t", "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"S", "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"h", "*", RowBox[{ RowBox[{"(", "T", ")"}], "^", "0.5"}]}], "+", RowBox[{"a", " ", "*", "R"}]}], ")"}], RowBox[{"(", RowBox[{"1", "-", "t"}], ")"}]}]}], ")"}], RowBox[{"y25", "[", "t", "]"}]}], "+", RowBox[{"y4", "[", "t", "]"}]}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"y35", "'"}], "[", "t", "]"}], "\[Equal]", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"a", "^", "2"}], ")"}], "*", RowBox[{"y29", "[", "t", "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"S", "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"h", "*", RowBox[{ RowBox[{"(", "T", ")"}], "^", "0.5"}]}], "+", RowBox[{"a", " ", "*", "R"}]}], ")"}], RowBox[{"(", RowBox[{"1", "-", "t"}], ")"}]}]}], ")"}], RowBox[{"y26", "[", "t", "]"}]}], "+", RowBox[{"y5", "[", "t", "]"}]}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"y36", "'"}], "[", "t", "]"}], "\[Equal]", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"a", "^", "2"}], ")"}], "*", RowBox[{"y30", "[", "t", "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"S", "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"h", "*", RowBox[{ RowBox[{"(", "T", ")"}], "^", "0.5"}]}], "+", RowBox[{"a", " ", "*", "R"}]}], ")"}], RowBox[{"(", RowBox[{"1", "-", "t"}], ")"}]}]}], ")"}], RowBox[{"y27", "[", "t", "]"}]}], "+", RowBox[{"y6", "[", "t", "]"}]}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"y13", "[", "0", "]"}], "\[Equal]", RowBox[{"y16", "[", "0", "]"}], "\[Equal]", RowBox[{"y20", "[", "0", "]"}], "==", RowBox[{"y23", "[", "0", "]"}], "\[Equal]", RowBox[{"y33", "[", "0", "]"}], "==", RowBox[{"y36", "[", "0", "]"}], "\[Equal]", "1"}], ",", RowBox[{ RowBox[{"y1", "[", "0", "]"}], "\[Equal]", RowBox[{"y2", "[", "0", "]"}], "\[Equal]", RowBox[{"y3", "[", "0", "]"}], "\[Equal]", RowBox[{"y4", "[", "0", "]"}], "\[Equal]", RowBox[{"y5", "[", "0", "]"}], "\[Equal]", RowBox[{"y6", "[", "0", "]"}], "\[Equal]", RowBox[{"y7", "[", "0", "]"}], "\[Equal]", RowBox[{"y8", "[", "0", "]"}], "\[Equal]", RowBox[{"y9", "[", "0", "]"}], "\[Equal]", RowBox[{"y10", "[", "0", "]"}], "\[Equal]", RowBox[{"y11", "[", "0", "]"}], "\[Equal]", RowBox[{"y12", "[", "0", "]"}], "\[Equal]", RowBox[{"y14", "[", "0", "]"}], "\[Equal]", RowBox[{"y15", "[", "0", "]"}], "\[Equal]", RowBox[{"y17", "[", "0", "]"}], "\[Equal]", RowBox[{"y18", "[", "0", "]"}], "\[Equal]", RowBox[{"y19", "[", "0", "]"}], "\[Equal]", RowBox[{"y21", "[", "0", "]"}], "\[Equal]", RowBox[{"y22", "[", "0", "]"}], "\[Equal]", RowBox[{"y24", "[", "0", "]"}], "\[Equal]", RowBox[{"y25", "[", "0", "]"}], "\[Equal]", RowBox[{"y26", "[", "0", "]"}], "\[Equal]", RowBox[{"y27", "[", "0", "]"}], "\[Equal]", RowBox[{"y28", "[", "0", "]"}], "\[Equal]", RowBox[{"y29", "[", "0", "]"}], "\[Equal]", RowBox[{"y30", "[", "0", "]"}], "\[Equal]", RowBox[{"y31", "[", "0", "]"}], "\[Equal]", RowBox[{"y32", "[", "0", "]"}], "\[Equal]", RowBox[{"y34", "[", "0", "]"}], "\[Equal]", RowBox[{"y35", "[", "0", "]"}], "==", "0"}]}], "}"}], ",", RowBox[{"{", RowBox[{ "y1", ",", "y2", ",", "y3", ",", "y4", ",", "y5", ",", "y6", ",", "y7", ",", "y8", ",", "y9", ",", "y10", ",", "y11", ",", "y12", ",", "y13", ",", "y14", ",", "y15", ",", "y16", ",", "y17", ",", "y18", ",", "y19", ",", "y20", ",", "y21", ",", "y22", ",", "y23", ",", "y24", ",", "y25", ",", "y26", ",", "y27", ",", "y28", ",", "y29", ",", "y30", ",", "y31", ",", "y32", ",", "y33", ",", "y34", ",", "y35", ",", "y36"}], "}"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "1"}], "}"}]}], "]"}]}]}], "]"}], "\[IndentingNewLine]", "c"}], "=", RowBox[{"Det", "[", RowBox[{"First", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"y1", "[", "1", "]"}], "+", RowBox[{"I", " ", RowBox[{"y4", "[", "1", "]"}]}]}], " ", ",", RowBox[{ RowBox[{"y2", "[", "1", "]"}], "+", RowBox[{"I", " ", RowBox[{"y5", "[", "1", "]"}]}]}], ",", RowBox[{ RowBox[{"y3", "[", "1", "]"}], "+", RowBox[{"I", " ", RowBox[{"y6", "[", "1", "]"}]}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"y7", "[", "1", "]"}], "+", RowBox[{"I", " ", RowBox[{"y10", "[", "1", "]"}]}]}], ",", RowBox[{ RowBox[{"y8", "[", "1", "]"}], "+", RowBox[{"I", " ", RowBox[{"y11", "[", "1", "]"}]}]}], ",", RowBox[{ RowBox[{"y9", "[", "1", "]"}], "+", RowBox[{"I", " ", RowBox[{"y12", "[", "1", "]"}]}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"y25", "[", "1", "]"}], "+", " ", RowBox[{"I", " ", RowBox[{"y28", "[", "1", "]"}]}]}], ",", RowBox[{ RowBox[{"y26", "[", "1", "]"}], "+", RowBox[{"I", " ", RowBox[{"y29", "[", "1", "]"}]}]}], ",", RowBox[{ RowBox[{"y27", "[", "1", "]"}], "+", " ", RowBox[{"I", " ", RowBox[{"y30", "[", "1", "]"}]}]}]}], "}"}]}], "}"}], "/.", "m"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"d", "=", RowBox[{"Abs", "[", "c", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"g", "=", RowBox[{"d", "^", "2"}]}], ";", "\[IndentingNewLine]", RowBox[{"Print", "[", RowBox[{ "R", ",", "\"\< \>\"", ",", "S", ",", "\"\< \>\"", ",", "T", ",", "\"\< \>\"", ",", "g"}], "]"}]}]}], "]"}]}], "]"}]}], "]"}]}], "Input", CellChangeTimes->{{3.5948140884701176`*^9, 3.5948142447043924`*^9}, { 3.5948145457069206`*^9, 3.594814580214182*^9}, {3.5948146200410514`*^9, 3.5948146246898594`*^9}}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"Set", "::", "write"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Tag \[NoBreak]\\!\\(Times\\)\[NoBreak] in \[NoBreak]\\!\\(\ \\(\\((\\(\\(45.18695441796399` \[VeryThinSpace]\\)\\) - \ \\(\\(43.83242278204331`\\\\ \[ImaginaryI]\\)\\))\\)\\)\\\\ Null\\)\\n\ \[NoBreak] is Protected. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/General/write\\\", ButtonNote -> \ \\\"Set::write\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.594814633301075*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "5", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "6", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "4100", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "3963.1421365149445`"}], SequenceForm[ 5, " ", 6, " ", 4100, " ", 3963.1421365149445`], Editable->False]], "Print", CellChangeTimes->{3.5948146333166747`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"Set", "::", "write"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Tag \[NoBreak]\\!\\(Times\\)\[NoBreak] in \[NoBreak]\\!\\(\ \\(\\((\\(\\(45.18695441796399` \[VeryThinSpace]\\)\\) - \ \\(\\(43.83242278204331`\\\\ \[ImaginaryI]\\)\\))\\)\\)\\\\ Null\\)\\n\ \[NoBreak] is Protected. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/General/write\\\", ButtonNote -> \ \\\"Set::write\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5948146335506754`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "5", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "6", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "4100", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "3963.1421365149445`"}], SequenceForm[ 5, " ", 6, " ", 4100, " ", 3963.1421365149445`], Editable->False]], "Print", CellChangeTimes->{3.5948146335506754`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"Set", "::", "write"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Tag \[NoBreak]\\!\\(Times\\)\[NoBreak] in \[NoBreak]\\!\\(\ \\(\\((\\(\\(45.18695441796399` \[VeryThinSpace]\\)\\) - \ \\(\\(43.83242278204331`\\\\ \[ImaginaryI]\\)\\))\\)\\)\\\\ Null\\)\\n\ \[NoBreak] is Protected. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/General/write\\\", ButtonNote -> \ \\\"Set::write\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5948146338158755`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"General", "::", "stop"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Further output of \[NoBreak]\\!\\(\\*StyleBox[\\(Set :: \ write\\), \\\"MessageName\\\"]\\)\[NoBreak] will be suppressed during this \ calculation. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/General/stop\\\", ButtonNote -> \ \\\"General::stop\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5948146338158755`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "5", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "6", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "4100", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "3963.1421365149445`"}], SequenceForm[ 5, " ", 6, " ", 4100, " ", 3963.1421365149445`], Editable->False]], "Print", CellChangeTimes->{3.5948146338158755`*^9}] }, Open ]] }, WindowSize->{1350, 611}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, FrontEndVersion->"8.0 for Microsoft Windows (64-bit) (February 23, 2011)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[579, 22, 30359, 686, 612, "Input"], Cell[30941, 710, 590, 12, 39, "Message"], Cell[31534, 724, 451, 10, 23, "Print"], Cell[31988, 736, 592, 12, 39, "Message"], Cell[32583, 750, 451, 10, 23, "Print"], Cell[33037, 762, 592, 12, 39, "Message"], Cell[33632, 776, 529, 11, 21, "Message"], Cell[34164, 789, 451, 10, 23, "Print"] }, Open ]] } ] *) (* End of internal cache information *)