(* Content-type: application/mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 6.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 145, 7] NotebookDataLength[ 83369, 2194] NotebookOptionsPosition[ 80138, 2082] NotebookOutlinePosition[ 80479, 2097] CellTagsIndexPosition[ 80436, 2094] WindowFrame->Normal ContainsDynamic->False*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[TextData[{ StyleBox["An Improved Method for Lyapunov Exponents Computation \n", "Subsubtitle", FontSize->18, FontVariations->{"CompatibilityType"->0}], StyleBox["Authors' Data", FontSize->14, FontWeight->"Bold"], StyleBox[": ", FontSize->14], StyleBox["Housam Binous & Nasri Zakia\nDepartment of Chemical Engineering\n\ National Institute of Applied Sciences and Technology\nTunis, TUNISIA\nEmail: \ binoushousam@yahoo.com ", FontSize->14, FontWeight->"Plain"] }], "Title", CellChangeTimes->{{3.4000354481602945`*^9, 3.400035458715472*^9}, 3.400035572909675*^9, {3.407590174329053*^9, 3.4075902150375886`*^9}, { 3.4130982131726866`*^9, 3.413098218971024*^9}, {3.4132653420031567`*^9, 3.413265386106574*^9}, 3.413288850406581*^9, {3.4140652190176115`*^9, 3.414065233588563*^9}}, TextAlignment->Center, Background->RGBColor[0.605478, 0.996109, 0.605478]], Cell["\<\ Acknowlegment : 1- Sandri, M. \"Numerical Calculation of Lyapunov Exponents.\" Mathematica J. \ 6, 78-84, 1996. http://library.wolfram.com/infocenter/Articles/2902/. \ \>", "Subsubtitle", CellChangeTimes->{{3.4132872054011793`*^9, 3.413287212130856*^9}, { 3.4132874243259773`*^9, 3.413287452055851*^9}, 3.4132877204017134`*^9, { 3.4132877553519697`*^9, 3.4132878168904576`*^9}, {3.4140652451752243`*^9, 3.4140652516745696`*^9}}, Background->RGBColor[0.7176470588235294, 1., 0.7176470588235294]], Cell["\<\ Test with Lorentz system. We get a chaotic behavior for this choice of \ parametres. Max lyapunov exponent is positive = 1.48\ \>", "Subsubtitle", CellGroupingRules->{GroupTogetherGrouping, 10000.}, CellChangeTimes->{{3.4130981153820705`*^9, 3.413098206422981*^9}, { 3.413098785766035*^9, 3.413098785836136*^9}, {3.413256895657909*^9, 3.413256900094288*^9}, 3.4132867526000834`*^9, {3.413288451913576*^9, 3.41328845285493*^9}, 3.4132889032725983`*^9, {3.414065274076782*^9, 3.4140652748078337`*^9}}, Background->RGBColor[0, 1, 1]], Cell[CellGroupData[{ Cell["\<\ Compute estimate of maximum Lyapunov exponent using first method\ \>", "Subsubtitle", CellChangeTimes->{{3.4140653243490705`*^9, 3.4140653474623055`*^9}}, Background->RGBColor[1., 0.7568627450980392, 1.]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"F", "[", RowBox[{"{", RowBox[{"x_", ",", "y_", ",", "z_"}], "}"}], "]"}], ":=", RowBox[{"{", RowBox[{ RowBox[{"16", " ", RowBox[{"(", RowBox[{"y", "-", "x"}], ")"}]}], ",", RowBox[{ RowBox[{"x", " ", RowBox[{"(", RowBox[{"45.92", "-", "z"}], ")"}]}], "-", "y"}], ",", RowBox[{ RowBox[{"x", " ", "y"}], " ", "-", RowBox[{"4", "z"}]}]}], "}"}]}]], "Input", CellGroupingRules->{GroupTogetherGrouping, 10000.}, CellChangeTimes->{{3.4130930611143904`*^9, 3.4130931183667154`*^9}, 3.4132867526000834`*^9}], Cell[BoxData[ RowBox[{ RowBox[{"JacobianMatrix", "[", RowBox[{"funs_List", ",", " ", "vars_List"}], "]"}], " ", ":=", " ", RowBox[{"Outer", "[", RowBox[{"D", ",", " ", "funs", ",", " ", "vars"}], "]"}]}]], "Input", CellGroupingRules->{GroupTogetherGrouping, 10000.}, CellChangeTimes->{3.4132867526000834`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"J", "=", RowBox[{"JacobianMatrix", "[", RowBox[{ RowBox[{"F", "[", RowBox[{"{", RowBox[{ RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}], ",", RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}], ",", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}]}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}], ",", RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}], ",", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}]}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.414042202040824*^9, 3.4140422167720065`*^9}, { 3.4140424828445997`*^9, 3.414042483265205*^9}, {3.4140426098372064`*^9, 3.4140426323696065`*^9}, {3.414042831806382*^9, 3.414042872304616*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "16"}], ",", "16", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"45.92`", "\[InvisibleSpace]", "-", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}]}], ",", RowBox[{"-", "1"}], ",", RowBox[{"-", RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}], ",", RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}], ",", RowBox[{"-", "4"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.414042217513072*^9, 3.41404248384604*^9, 3.4140426334912195`*^9, 3.4140428733160706`*^9, 3.414045037788429*^9, 3.4140475680968337`*^9, 3.4140540724796767`*^9, 3.4140543322732415`*^9, 3.4140652852027807`*^9, 3.4140658047198095`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"J", "//", "MatrixForm"}]], "Input", CellChangeTimes->{{3.414042518666109*^9, 3.414042522451552*^9}}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{"-", "16"}], "16", "0"}, { RowBox[{"45.92`", "\[InvisibleSpace]", "-", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}]}], RowBox[{"-", "1"}], RowBox[{"-", RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}]}]}, { RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}], RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}], RowBox[{"-", "4"}]} }, GridBoxAlignment->{ "Columns" -> {{Left}}, "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.4140425230424013`*^9, 3.4140426395299025`*^9, 3.414042882469232*^9, 3.414045038058818*^9, 3.41404756899813*^9, 3.4140540772465315`*^9, 3.414054333084408*^9, 3.414065286334408*^9, 3.414065805010227*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Y", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ SubscriptBox["y", "i"], "[", "t", "]"}], ",", RowBox[{ SubscriptBox["y", RowBox[{"i", "+", "3"}]], "[", "t", "]"}], ",", RowBox[{ SubscriptBox["y", RowBox[{"i", "+", "6"}]], "[", "t", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "4", ",", "6"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.4140422737739716`*^9, 3.414042428135933*^9}, 3.414042886645237*^9, {3.41404335084272*^9, 3.414043375277856*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ SubscriptBox["y", "4"], "[", "t", "]"}], ",", RowBox[{ SubscriptBox["y", "7"], "[", "t", "]"}], ",", RowBox[{ SubscriptBox["y", "10"], "[", "t", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ SubscriptBox["y", "5"], "[", "t", "]"}], ",", RowBox[{ SubscriptBox["y", "8"], "[", "t", "]"}], ",", RowBox[{ SubscriptBox["y", "11"], "[", "t", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ SubscriptBox["y", "6"], "[", "t", "]"}], ",", RowBox[{ SubscriptBox["y", "9"], "[", "t", "]"}], ",", RowBox[{ SubscriptBox["y", "12"], "[", "t", "]"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.414042429007186*^9, 3.414043375898749*^9, 3.41404504193439*^9, 3.414047569739195*^9, 3.4140540788388205`*^9, 3.414054333915603*^9, 3.4140652874059486`*^9, 3.4140658052405586`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"%", "//", "MatrixForm"}]], "Input", CellChangeTimes->{{3.41404243693859*^9, 3.4140424411947107`*^9}}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{ SubscriptBox["y", "4"], "[", "t", "]"}], RowBox[{ SubscriptBox["y", "7"], "[", "t", "]"}], RowBox[{ SubscriptBox["y", "10"], "[", "t", "]"}]}, { RowBox[{ SubscriptBox["y", "5"], "[", "t", "]"}], RowBox[{ SubscriptBox["y", "8"], "[", "t", "]"}], RowBox[{ SubscriptBox["y", "11"], "[", "t", "]"}]}, { RowBox[{ SubscriptBox["y", "6"], "[", "t", "]"}], RowBox[{ SubscriptBox["y", "9"], "[", "t", "]"}], RowBox[{ SubscriptBox["y", "12"], "[", "t", "]"}]} }, GridBoxAlignment->{ "Columns" -> {{Left}}, "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.4140424419658194`*^9, 3.414043378262147*^9, 3.414045042244837*^9, 3.414047570540347*^9, 3.4140540798903327`*^9, 3.41405433463664*^9, 3.4140652883773456`*^9, 3.4140658078142595`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"J", ".", "Y"}], "//", "MatrixForm"}]], "Input", CellChangeTimes->{{3.4140436943065968`*^9, 3.414043700144992*^9}}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{ RowBox[{ RowBox[{"-", "16"}], " ", RowBox[{ SubscriptBox["y", "4"], "[", "t", "]"}]}], "+", RowBox[{"16", " ", RowBox[{ SubscriptBox["y", "5"], "[", "t", "]"}]}]}], RowBox[{ RowBox[{ RowBox[{"-", "16"}], " ", RowBox[{ SubscriptBox["y", "7"], "[", "t", "]"}]}], "+", RowBox[{"16", " ", RowBox[{ SubscriptBox["y", "8"], "[", "t", "]"}]}]}], RowBox[{ RowBox[{ RowBox[{"-", "16"}], " ", RowBox[{ SubscriptBox["y", "10"], "[", "t", "]"}]}], "+", RowBox[{"16", " ", RowBox[{ SubscriptBox["y", "11"], "[", "t", "]"}]}]}]}, { RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"45.92`", "\[InvisibleSpace]", "-", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}]}], ")"}], " ", RowBox[{ SubscriptBox["y", "4"], "[", "t", "]"}]}], "-", RowBox[{ SubscriptBox["y", "5"], "[", "t", "]"}], "-", RowBox[{ RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "6"], "[", "t", "]"}]}]}], RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"45.92`", "\[InvisibleSpace]", "-", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}]}], ")"}], " ", RowBox[{ SubscriptBox["y", "7"], "[", "t", "]"}]}], "-", RowBox[{ SubscriptBox["y", "8"], "[", "t", "]"}], "-", RowBox[{ RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "9"], "[", "t", "]"}]}]}], RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"45.92`", "\[InvisibleSpace]", "-", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}]}], ")"}], " ", RowBox[{ SubscriptBox["y", "10"], "[", "t", "]"}]}], "-", RowBox[{ SubscriptBox["y", "11"], "[", "t", "]"}], "-", RowBox[{ RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "12"], "[", "t", "]"}]}]}]}, { RowBox[{ RowBox[{ RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "4"], "[", "t", "]"}]}], "+", RowBox[{ RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "5"], "[", "t", "]"}]}], "-", RowBox[{"4", " ", RowBox[{ SubscriptBox["y", "6"], "[", "t", "]"}]}]}], RowBox[{ RowBox[{ RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "7"], "[", "t", "]"}]}], "+", RowBox[{ RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "8"], "[", "t", "]"}]}], "-", RowBox[{"4", " ", RowBox[{ SubscriptBox["y", "9"], "[", "t", "]"}]}]}], RowBox[{ RowBox[{ RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "10"], "[", "t", "]"}]}], "+", RowBox[{ RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "11"], "[", "t", "]"}]}], "-", RowBox[{"4", " ", RowBox[{ SubscriptBox["y", "12"], "[", "t", "]"}]}]}]} }, GridBoxAlignment->{ "Columns" -> {{Left}}, "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.414043700815957*^9, 3.4140450424851823`*^9, 3.4140475713414993`*^9, 3.4140540808417006`*^9, 3.414054335487864*^9, 3.414065289488944*^9, 3.4140658090760736`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"E1", "=", RowBox[{"Flatten", "[", RowBox[{"Transpose", "[", RowBox[{"J", ".", "Y"}], "]"}], "]"}]}]], "Input", CellChangeTimes->{{3.4140424861393375`*^9, 3.4140424969348607`*^9}, { 3.41404254802833*^9, 3.4140425486091647`*^9}, {3.4140435488073792`*^9, 3.4140435706688147`*^9}, {3.4140446072293153`*^9, 3.414044620348179*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{"-", "16"}], " ", RowBox[{ SubscriptBox["y", "4"], "[", "t", "]"}]}], "+", RowBox[{"16", " ", RowBox[{ SubscriptBox["y", "5"], "[", "t", "]"}]}]}], ",", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"45.92`", "\[InvisibleSpace]", "-", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}]}], ")"}], " ", RowBox[{ SubscriptBox["y", "4"], "[", "t", "]"}]}], "-", RowBox[{ SubscriptBox["y", "5"], "[", "t", "]"}], "-", RowBox[{ RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "6"], "[", "t", "]"}]}]}], ",", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "4"], "[", "t", "]"}]}], "+", RowBox[{ RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "5"], "[", "t", "]"}]}], "-", RowBox[{"4", " ", RowBox[{ SubscriptBox["y", "6"], "[", "t", "]"}]}]}], ",", RowBox[{ RowBox[{ RowBox[{"-", "16"}], " ", RowBox[{ SubscriptBox["y", "7"], "[", "t", "]"}]}], "+", RowBox[{"16", " ", RowBox[{ SubscriptBox["y", "8"], "[", "t", "]"}]}]}], ",", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"45.92`", "\[InvisibleSpace]", "-", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}]}], ")"}], " ", RowBox[{ SubscriptBox["y", "7"], "[", "t", "]"}]}], "-", RowBox[{ SubscriptBox["y", "8"], "[", "t", "]"}], "-", RowBox[{ RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "9"], "[", "t", "]"}]}]}], ",", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "7"], "[", "t", "]"}]}], "+", RowBox[{ RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "8"], "[", "t", "]"}]}], "-", RowBox[{"4", " ", RowBox[{ SubscriptBox["y", "9"], "[", "t", "]"}]}]}], ",", RowBox[{ RowBox[{ RowBox[{"-", "16"}], " ", RowBox[{ SubscriptBox["y", "10"], "[", "t", "]"}]}], "+", RowBox[{"16", " ", RowBox[{ SubscriptBox["y", "11"], "[", "t", "]"}]}]}], ",", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"45.92`", "\[InvisibleSpace]", "-", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}]}], ")"}], " ", RowBox[{ SubscriptBox["y", "10"], "[", "t", "]"}]}], "-", RowBox[{ SubscriptBox["y", "11"], "[", "t", "]"}], "-", RowBox[{ RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "12"], "[", "t", "]"}]}]}], ",", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "10"], "[", "t", "]"}]}], "+", RowBox[{ RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "11"], "[", "t", "]"}]}], "-", RowBox[{"4", " ", RowBox[{ SubscriptBox["y", "12"], "[", "t", "]"}]}]}]}], "}"}]], "Output", CellChangeTimes->{ 3.414042498146603*^9, 3.41404254944036*^9, 3.4140426438360944`*^9, 3.4140433800847683`*^9, {3.4140435590621247`*^9, 3.4140435715500813`*^9}, 3.414044621539893*^9, 3.4140450426053553`*^9, 3.4140475721526656`*^9, 3.414054081923256*^9, 3.414054336289016*^9, 3.414065290590528*^9, 3.414065810337888*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"EQ3", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{ SubscriptBox["y", "i"], "[", "t", "]"}], ",", RowBox[{"{", RowBox[{"t", ",", "1"}], "}"}]}], "]"}], "\[Equal]", " ", RowBox[{ RowBox[{"F", "[", RowBox[{"{", RowBox[{ RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}], ",", RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}], ",", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}]}], "}"}], "]"}], "[", RowBox[{"[", "i", "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.4140426749808784`*^9, 3.414042751681168*^9}, { 3.4140433882364893`*^9, 3.4140434096572914`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ SuperscriptBox[ SubscriptBox["y", "1"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", RowBox[{"16", " ", RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}]}], "+", RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}]}], ")"}]}]}], ",", RowBox[{ RowBox[{ SuperscriptBox[ SubscriptBox["y", "2"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", RowBox[{ RowBox[{"-", RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}]}], "+", RowBox[{ RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}], " ", RowBox[{"(", RowBox[{"45.92`", "\[InvisibleSpace]", "-", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}]}], ")"}]}]}]}], ",", RowBox[{ RowBox[{ SuperscriptBox[ SubscriptBox["y", "3"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}]}], "-", RowBox[{"4", " ", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}]}]}]}]}], "}"}]], "Output", CellChangeTimes->{{3.4140427431889567`*^9, 3.414042752632536*^9}, { 3.4140433824381523`*^9, 3.414043410628688*^9}, 3.4140446239333344`*^9, 3.414045042725528*^9, 3.414047572963832*^9, 3.414054084206539*^9, 3.414054337090168*^9, 3.414065291692112*^9, 3.4140658114795294`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"EQ9", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{ SubscriptBox["y", "i"], "[", "t", "]"}], ",", RowBox[{"{", RowBox[{"t", ",", "1"}], "}"}]}], "]"}], "\[Equal]", RowBox[{"E1", "[", RowBox[{"[", RowBox[{"i", "-", "3"}], "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "4", ",", "12"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.4140435126553955`*^9, 3.414043529649832*^9}, { 3.4140435769979153`*^9, 3.4140436056991854`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ SuperscriptBox[ SubscriptBox["y", "4"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", RowBox[{ RowBox[{ RowBox[{"-", "16"}], " ", RowBox[{ SubscriptBox["y", "4"], "[", "t", "]"}]}], "+", RowBox[{"16", " ", RowBox[{ SubscriptBox["y", "5"], "[", "t", "]"}]}]}]}], ",", RowBox[{ RowBox[{ SuperscriptBox[ SubscriptBox["y", "5"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"45.92`", "\[InvisibleSpace]", "-", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}]}], ")"}], " ", RowBox[{ SubscriptBox["y", "4"], "[", "t", "]"}]}], "-", RowBox[{ SubscriptBox["y", "5"], "[", "t", "]"}], "-", RowBox[{ RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "6"], "[", "t", "]"}]}]}]}], ",", RowBox[{ RowBox[{ SuperscriptBox[ SubscriptBox["y", "6"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "4"], "[", "t", "]"}]}], "+", RowBox[{ RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "5"], "[", "t", "]"}]}], "-", RowBox[{"4", " ", RowBox[{ SubscriptBox["y", "6"], "[", "t", "]"}]}]}]}], ",", RowBox[{ RowBox[{ SuperscriptBox[ SubscriptBox["y", "7"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", RowBox[{ RowBox[{ RowBox[{"-", "16"}], " ", RowBox[{ SubscriptBox["y", "7"], "[", "t", "]"}]}], "+", RowBox[{"16", " ", RowBox[{ SubscriptBox["y", "8"], "[", "t", "]"}]}]}]}], ",", RowBox[{ RowBox[{ SuperscriptBox[ SubscriptBox["y", "8"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"45.92`", "\[InvisibleSpace]", "-", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}]}], ")"}], " ", RowBox[{ SubscriptBox["y", "7"], "[", "t", "]"}]}], "-", RowBox[{ SubscriptBox["y", "8"], "[", "t", "]"}], "-", RowBox[{ RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "9"], "[", "t", "]"}]}]}]}], ",", RowBox[{ RowBox[{ SuperscriptBox[ SubscriptBox["y", "9"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "7"], "[", "t", "]"}]}], "+", RowBox[{ RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "8"], "[", "t", "]"}]}], "-", RowBox[{"4", " ", RowBox[{ SubscriptBox["y", "9"], "[", "t", "]"}]}]}]}], ",", RowBox[{ RowBox[{ SuperscriptBox[ SubscriptBox["y", "10"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", RowBox[{ RowBox[{ RowBox[{"-", "16"}], " ", RowBox[{ SubscriptBox["y", "10"], "[", "t", "]"}]}], "+", RowBox[{"16", " ", RowBox[{ SubscriptBox["y", "11"], "[", "t", "]"}]}]}]}], ",", RowBox[{ RowBox[{ SuperscriptBox[ SubscriptBox["y", "11"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"45.92`", "\[InvisibleSpace]", "-", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}]}], ")"}], " ", RowBox[{ SubscriptBox["y", "10"], "[", "t", "]"}]}], "-", RowBox[{ SubscriptBox["y", "11"], "[", "t", "]"}], "-", RowBox[{ RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "12"], "[", "t", "]"}]}]}]}], ",", RowBox[{ RowBox[{ SuperscriptBox[ SubscriptBox["y", "12"], "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "10"], "[", "t", "]"}]}], "+", RowBox[{ RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}], " ", RowBox[{ SubscriptBox["y", "11"], "[", "t", "]"}]}], "-", RowBox[{"4", " ", RowBox[{ SubscriptBox["y", "12"], "[", "t", "]"}]}]}]}]}], "}"}]], "Output", CellChangeTimes->{3.4140436064202223`*^9, 3.4140446249548035`*^9, 3.4140450428156576`*^9, 3.4140475737950273`*^9, 3.414054085348181*^9, 3.4140543379313774`*^9, 3.414065292693552*^9, 3.4140658125410557`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"YI9", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["y", "i"], "[", "0", "]"}], "\[Equal]", " ", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Mod", "[", RowBox[{"i", ",", "4"}], "]"}], "\[Equal]", " ", "0"}], ",", "1", ",", "0"}], "]"}]}], ",", RowBox[{ RowBox[{ SubscriptBox["y", RowBox[{"i", "+", "3"}]], "[", "0", "]"}], "\[Equal]", " ", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Mod", "[", RowBox[{ RowBox[{"i", "+", "3"}], ",", "4"}], "]"}], "\[Equal]", " ", "0"}], ",", "1", ",", "0"}], "]"}]}], ",", RowBox[{ RowBox[{ SubscriptBox["y", RowBox[{"i", "+", "6"}]], "[", "0", "]"}], "\[Equal]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Mod", "[", RowBox[{ RowBox[{"i", "+", "6"}], ",", "4"}], "]"}], "\[Equal]", " ", "0"}], ",", "1", ",", "0"}], "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "4", ",", "6"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.4140437616233935`*^9, 3.414043788251683*^9}, { 3.4140438379531507`*^9, 3.414043962121696*^9}, {3.4140440405244336`*^9, 3.4140440466432323`*^9}, {3.414044142110507*^9, 3.4140441429917746`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["y", "4"], "[", "0", "]"}], "\[Equal]", "1"}], ",", RowBox[{ RowBox[{ SubscriptBox["y", "7"], "[", "0", "]"}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{ SubscriptBox["y", "10"], "[", "0", "]"}], "\[Equal]", "0"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["y", "5"], "[", "0", "]"}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{ SubscriptBox["y", "8"], "[", "0", "]"}], "\[Equal]", "1"}], ",", RowBox[{ RowBox[{ SubscriptBox["y", "11"], "[", "0", "]"}], "\[Equal]", "0"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["y", "6"], "[", "0", "]"}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{ SubscriptBox["y", "9"], "[", "0", "]"}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{ SubscriptBox["y", "12"], "[", "0", "]"}], "\[Equal]", "1"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.4140439627926607`*^9, 3.414044047724787*^9, 3.41404414357261*^9, 3.4140446265070353`*^9, 3.414045042955859*^9, 3.414047574596179*^9, 3.4140495880013137`*^9, 3.414049882905365*^9, 3.414049948008979*^9, 3.4140502243863907`*^9, 3.414050296029408*^9, 3.414053123274787*^9, 3.414053198382787*^9, 3.414054086309563*^9, 3.414054338692472*^9, 3.414065293745064*^9, 3.4140658136326256`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"YI3", "=", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["y", "1"], "[", "0", "]"}], "\[Equal]", " ", "19"}], ",", RowBox[{ RowBox[{ SubscriptBox["y", "2"], "[", "0", "]"}], "\[Equal]", " ", "20"}], ",", RowBox[{ RowBox[{ SubscriptBox["y", "3"], "[", "0", "]"}], "\[Equal]", " ", "50"}]}], "}"}]}]], "Input", CellChangeTimes->{{3.41404414548536*^9, 3.41404418303936*^9}, { 3.4140449056684494`*^9, 3.414044911917435*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["y", "1"], "[", "0", "]"}], "\[Equal]", "19"}], ",", RowBox[{ RowBox[{ SubscriptBox["y", "2"], "[", "0", "]"}], "\[Equal]", "20"}], ",", RowBox[{ RowBox[{ SubscriptBox["y", "3"], "[", "0", "]"}], "\[Equal]", "50"}]}], "}"}]], "Output", CellChangeTimes->{3.4140441836502385`*^9, 3.414044627558547*^9, 3.4140449123580685`*^9, 3.4140450433664494`*^9, 3.4140475754774466`*^9, 3.414049589092883*^9, 3.414049884758029*^9, 3.4140499490204334`*^9, 3.4140502252676573`*^9, 3.414050296910675*^9, 3.4140531240258675`*^9, 3.4140531990938096`*^9, 3.4140540874011326`*^9, 3.414054339573739*^9, 3.414065294746504*^9, 3.4140658147241955`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"sol", "=", RowBox[{"NDSolve", "[", RowBox[{ RowBox[{"Join", "[", RowBox[{"EQ9", ",", "EQ3", ",", "YI3", ",", "YI9"}], "]"}], ",", RowBox[{"Table", "[", RowBox[{ RowBox[{ SubscriptBox["y", "i"], "[", "t", "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "12"}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "40"}], "}"}], ",", RowBox[{"MaxSteps", "\[Rule]", " ", "Infinity"}], ",", RowBox[{"AccuracyGoal", "\[Rule]", "2"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.4140440598722544`*^9, 3.414044125396474*^9}, { 3.4140441907905054`*^9, 3.414044242154363*^9}, {3.4140443092408285`*^9, 3.414044328428419*^9}, {3.4140445281255693`*^9, 3.4140445324317617`*^9}, 3.414044632545718*^9, {3.414044666464491*^9, 3.414044683909576*^9}, { 3.414044721363432*^9, 3.414044757996107*^9}, {3.4140448004671774`*^9, 3.414044801458603*^9}, {3.4140449189074864`*^9, 3.4140449193481197`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}], "\[Rule]", RowBox[{ TagBox[ RowBox[{"InterpolatingFunction", "[", RowBox[{ RowBox[{"{", RowBox[{"{", RowBox[{"0.`", ",", "40.`"}], "}"}], "}"}], ",", "\<\"<>\"\>"}], "]"}], False, Editable->False], "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}], "\[Rule]", RowBox[{ TagBox[ RowBox[{"InterpolatingFunction", "[", RowBox[{ RowBox[{"{", RowBox[{"{", RowBox[{"0.`", ",", "40.`"}], "}"}], "}"}], ",", "\<\"<>\"\>"}], "]"}], False, Editable->False], "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}], "\[Rule]", RowBox[{ TagBox[ RowBox[{"InterpolatingFunction", "[", RowBox[{ RowBox[{"{", RowBox[{"{", RowBox[{"0.`", ",", "40.`"}], "}"}], "}"}], ",", "\<\"<>\"\>"}], "]"}], False, Editable->False], "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ SubscriptBox["y", "4"], "[", "t", "]"}], "\[Rule]", RowBox[{ TagBox[ RowBox[{"InterpolatingFunction", "[", RowBox[{ RowBox[{"{", RowBox[{"{", RowBox[{"0.`", ",", "40.`"}], "}"}], "}"}], ",", "\<\"<>\"\>"}], "]"}], False, Editable->False], "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ SubscriptBox["y", "5"], "[", "t", "]"}], "\[Rule]", RowBox[{ TagBox[ RowBox[{"InterpolatingFunction", "[", RowBox[{ RowBox[{"{", RowBox[{"{", RowBox[{"0.`", ",", "40.`"}], "}"}], "}"}], ",", "\<\"<>\"\>"}], "]"}], False, Editable->False], "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ SubscriptBox["y", "6"], "[", "t", "]"}], "\[Rule]", RowBox[{ TagBox[ RowBox[{"InterpolatingFunction", "[", RowBox[{ RowBox[{"{", RowBox[{"{", RowBox[{"0.`", ",", "40.`"}], "}"}], "}"}], ",", "\<\"<>\"\>"}], "]"}], False, Editable->False], "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ SubscriptBox["y", "7"], "[", "t", "]"}], "\[Rule]", RowBox[{ TagBox[ RowBox[{"InterpolatingFunction", "[", RowBox[{ RowBox[{"{", RowBox[{"{", RowBox[{"0.`", ",", "40.`"}], "}"}], "}"}], ",", "\<\"<>\"\>"}], "]"}], False, Editable->False], "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ SubscriptBox["y", "8"], "[", "t", "]"}], "\[Rule]", RowBox[{ TagBox[ RowBox[{"InterpolatingFunction", "[", RowBox[{ RowBox[{"{", RowBox[{"{", RowBox[{"0.`", ",", "40.`"}], "}"}], "}"}], ",", "\<\"<>\"\>"}], "]"}], False, Editable->False], "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ SubscriptBox["y", "9"], "[", "t", "]"}], "\[Rule]", RowBox[{ TagBox[ RowBox[{"InterpolatingFunction", "[", RowBox[{ RowBox[{"{", RowBox[{"{", RowBox[{"0.`", ",", "40.`"}], "}"}], "}"}], ",", "\<\"<>\"\>"}], "]"}], False, Editable->False], "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ SubscriptBox["y", "10"], "[", "t", "]"}], "\[Rule]", RowBox[{ TagBox[ RowBox[{"InterpolatingFunction", "[", RowBox[{ RowBox[{"{", RowBox[{"{", RowBox[{"0.`", ",", "40.`"}], "}"}], "}"}], ",", "\<\"<>\"\>"}], "]"}], False, Editable->False], "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ SubscriptBox["y", "11"], "[", "t", "]"}], "\[Rule]", RowBox[{ TagBox[ RowBox[{"InterpolatingFunction", "[", RowBox[{ RowBox[{"{", RowBox[{"{", RowBox[{"0.`", ",", "40.`"}], "}"}], "}"}], ",", "\<\"<>\"\>"}], "]"}], False, Editable->False], "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ SubscriptBox["y", "12"], "[", "t", "]"}], "\[Rule]", RowBox[{ TagBox[ RowBox[{"InterpolatingFunction", "[", RowBox[{ RowBox[{"{", RowBox[{"{", RowBox[{"0.`", ",", "40.`"}], "}"}], "}"}], ",", "\<\"<>\"\>"}], "]"}], False, Editable->False], "[", "t", "]"}]}]}], "}"}], "}"}]], "Output", CellChangeTimes->{ 3.414044126257712*^9, {3.4140441851824417`*^9, 3.414044194936467*^9}, { 3.4140443070777187`*^9, 3.414044329169485*^9}, 3.414044529767931*^9, 3.4140445986870317`*^9, 3.414044633897662*^9, {3.4140446687177315`*^9, 3.4140446861327724`*^9}, {3.414044725289077*^9, 3.4140447593981233`*^9}, 3.414044803080936*^9, 3.414044859061432*^9, {3.414044914991856*^9, 3.4140449208803234`*^9}, 3.4140450456697617`*^9, 3.4140475781813345`*^9, 3.414049591636541*^9, 3.414054089614315*^9, 3.414054341626691*^9, 3.4140652978709965`*^9, 3.414065817017493*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"u", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Random", "[", "]"}], ",", RowBox[{"{", "3", "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.414044411247507*^9, 3.414044411728198*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ "0.6944940072974286`", ",", "0.2313875695644118`", ",", "0.557651959887384`"}], "}"}]], "Output", CellChangeTimes->{ 3.4140444122088895`*^9, {3.414044751216358*^9, 3.414044761551219*^9}, 3.4140448054743776`*^9, 3.414044860493491*^9, {3.4140449227329874`*^9, 3.414044936152283*^9}, 3.414045045880064*^9, 3.4140475793129616`*^9, 3.414047961492509*^9, 3.414049592928398*^9, 3.414054091156533*^9, 3.4140543426882176`*^9, 3.414065299353128*^9, 3.414065818519653*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"T", "=", "40"}], ";"}]], "Input", CellChangeTimes->{{3.414044429513773*^9, 3.41404443139648*^9}, { 3.4140447633437967`*^9, 3.4140447640347905`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"PhiT", "=", RowBox[{"Transpose", "[", RowBox[{"First", "[", RowBox[{ RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ SubscriptBox["y", "i"], "[", "t", "]"}], ",", RowBox[{ SubscriptBox["y", RowBox[{"i", "+", "3"}]], "[", "t", "]"}], ",", RowBox[{ SubscriptBox["y", RowBox[{"i", "+", "6"}]], "[", "t", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "4", ",", "6"}], "}"}]}], "]"}], "/.", "sol"}], "/.", RowBox[{"t", "\[Rule]", " ", "40"}]}], "]"}], "]"}]}]], "Input", CellChangeTimes->{{3.414044438576805*^9, 3.4140444624711633`*^9}, { 3.4140447679003487`*^9, 3.4140447684210978`*^9}, {3.4140448337951007`*^9, 3.414044846072755*^9}, {3.4140479199027057`*^9, 3.4140479451790514`*^9}, { 3.4140480199766045`*^9, 3.414048089636771*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "2.976914726906416`*^25"}], ",", RowBox[{"-", "4.558462650419199`*^25"}], ",", "3.594217768283253`*^23"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "5.1311318977585605`*^25"}], ",", RowBox[{"-", "7.857152540816572`*^25"}], ",", "6.195140650707902`*^23"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "2.470356748947348`*^25"}], ",", RowBox[{"-", "3.782785200901558`*^25"}], ",", "2.9826182257818362`*^23"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{{3.414048069768202*^9, 3.4140480902176065`*^9}, 3.414049594741005*^9, 3.4140540930192113`*^9, 3.4140543456424656`*^9, 3.4140653074047055`*^9, 3.4140658205025043`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"Log", "[", RowBox[{"Norm", "[", RowBox[{"PhiT", ".", "u"}], "]"}], "]"}], "/", "T"}]], "Input", CellChangeTimes->{{3.414044427070259*^9, 3.4140444352520237`*^9}}], Cell[BoxData["1.4866548292088604`"], "Output", CellChangeTimes->{ 3.414044466086362*^9, 3.414044770243718*^9, 3.414044809019475*^9, { 3.414044848466197*^9, 3.414044863317552*^9}, {3.414044925557048*^9, 3.4140449388862143`*^9}, 3.4140450460202656`*^9, 3.414047582807987*^9, { 3.4140479592993555`*^9, 3.41404796372572*^9}, 3.414048103376528*^9, 3.4140495958826466`*^9, 3.4140540954226675`*^9, 3.4140543485165987`*^9, 3.4140653086565056`*^9, 3.414065821724261*^9}] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["Compute all three Lyapunov exponents using second method", "Subsubtitle", CellChangeTimes->{{3.4140653243490705`*^9, 3.4140653474623055`*^9}, { 3.4140658923357954`*^9, 3.414065902400267*^9}}, Background->RGBColor[1., 0.7568627450980392, 1.]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"<<", "\"\\""}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"General", "::", "\<\"obspkg\"\>"}], RowBox[{ ":", " "}], "\<\"\\!\\(\\\"LinearAlgebra`Orthogonalization`\\\"\\) is now \ obsolete. The legacy version being loaded may conflict with current \ Mathematica functionality. See the Compatibility Guide for updating \ information. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:Compatibility/Tutorials/LinearAlgebra/\ Orthogonalization\\\", ButtonNote -> \\\"General::obspkg\\\"]\\)\"\>"}]], \ "Message", "MSG", CellChangeTimes->{3.414047603257392*^9, 3.4140481208015842`*^9, 3.4140540982867856`*^9, 3.4140543641891346`*^9, 3.414065779643752*^9, 3.4140658251591997`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"YI9", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["y", "i"], "[", "0", "]"}], "\[Equal]", " ", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Mod", "[", RowBox[{"i", ",", "4"}], "]"}], "\[Equal]", " ", "0"}], ",", "1", ",", "0"}], "]"}]}], ",", RowBox[{ RowBox[{ SubscriptBox["y", RowBox[{"i", "+", "3"}]], "[", "0", "]"}], "\[Equal]", " ", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Mod", "[", RowBox[{ RowBox[{"i", "+", "3"}], ",", "4"}], "]"}], "\[Equal]", " ", "0"}], ",", "1", ",", "0"}], "]"}]}], ",", RowBox[{ RowBox[{ SubscriptBox["y", RowBox[{"i", "+", "6"}]], "[", "0", "]"}], "\[Equal]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Mod", "[", RowBox[{ RowBox[{"i", "+", "6"}], ",", "4"}], "]"}], "\[Equal]", " ", "0"}], ",", "1", ",", "0"}], "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "4", ",", "6"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.4140437616233935`*^9, 3.414043788251683*^9}, { 3.4140438379531507`*^9, 3.414043962121696*^9}, {3.4140440405244336`*^9, 3.4140440466432323`*^9}, {3.414044142110507*^9, 3.4140441429917746`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["y", "4"], "[", "0", "]"}], "\[Equal]", "1"}], ",", RowBox[{ RowBox[{ SubscriptBox["y", "7"], "[", "0", "]"}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{ SubscriptBox["y", "10"], "[", "0", "]"}], "\[Equal]", "0"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["y", "5"], "[", "0", "]"}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{ SubscriptBox["y", "8"], "[", "0", "]"}], "\[Equal]", "1"}], ",", RowBox[{ RowBox[{ SubscriptBox["y", "11"], "[", "0", "]"}], "\[Equal]", "0"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["y", "6"], "[", "0", "]"}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{ SubscriptBox["y", "9"], "[", "0", "]"}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{ SubscriptBox["y", "12"], "[", "0", "]"}], "\[Equal]", "1"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{ 3.4140439627926607`*^9, 3.414044047724787*^9, 3.41404414357261*^9, 3.4140446265070353`*^9, 3.414045042955859*^9, 3.414047574596179*^9, 3.4140495880013137`*^9, 3.414049882905365*^9, 3.414049948008979*^9, 3.4140502243863907`*^9, 3.414050296029408*^9, 3.414053123274787*^9, 3.414053198382787*^9, {3.414053280450795*^9, 3.414053297785722*^9}, 3.4140534190100336`*^9, 3.414053461871666*^9, {3.414053655119542*^9, 3.414053676590416*^9}, 3.414053844101285*^9, 3.414054108381301*^9, 3.4140543661118994`*^9, 3.414054490761136*^9, 3.414054574151045*^9, 3.4140546756169453`*^9, 3.4140554369617057`*^9, 3.4140658264510574`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"YI3", "=", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["y", "1"], "[", "0", "]"}], "\[Equal]", " ", "19"}], ",", RowBox[{ RowBox[{ SubscriptBox["y", "2"], "[", "0", "]"}], "\[Equal]", " ", "20"}], ",", RowBox[{ RowBox[{ SubscriptBox["y", "3"], "[", "0", "]"}], "\[Equal]", " ", "50"}]}], "}"}]}]], "Input", CellChangeTimes->{{3.41404414548536*^9, 3.41404418303936*^9}, { 3.4140449056684494`*^9, 3.414044911917435*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["y", "1"], "[", "0", "]"}], "\[Equal]", "19"}], ",", RowBox[{ RowBox[{ SubscriptBox["y", "2"], "[", "0", "]"}], "\[Equal]", "20"}], ",", RowBox[{ RowBox[{ SubscriptBox["y", "3"], "[", "0", "]"}], "\[Equal]", "50"}]}], "}"}]], "Output", CellChangeTimes->{ 3.4140441836502385`*^9, 3.414044627558547*^9, 3.4140449123580685`*^9, 3.4140450433664494`*^9, 3.4140475754774466`*^9, 3.414049589092883*^9, 3.414049884758029*^9, 3.4140499490204334`*^9, 3.4140502252676573`*^9, 3.414050296910675*^9, 3.4140531240258675`*^9, 3.4140531990938096`*^9, { 3.4140532811618176`*^9, 3.4140532985368013`*^9}, 3.4140534198011713`*^9, 3.4140534625826883`*^9, {3.4140536558305645`*^9, 3.4140536774716835`*^9}, 3.4140538449024367`*^9, 3.414054112216816*^9, 3.4140543671233535`*^9, 3.414054491812648*^9, 3.414054575022298*^9, 3.4140546763780403`*^9, 3.414055437802915*^9, 3.4140658275025697`*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"s", "=", RowBox[{"{", "}"}]}], ";"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"K", "=", "500"}], ";"}]], "Input", CellChangeTimes->{{3.4140476346124783`*^9, 3.4140476369858913`*^9}, { 3.4140481645344687`*^9, 3.4140481659765425`*^9}, {3.4140495769053583`*^9, 3.4140495784676046`*^9}, {3.414049724127053*^9, 3.4140497246478014`*^9}, { 3.414049890816741*^9, 3.4140498910270433`*^9}, {3.414049953587*^9, 3.414049954868843*^9}, {3.414053134040267*^9, 3.414053134611088*^9}, { 3.4140532046417875`*^9, 3.4140532048120317`*^9}, {3.414053261463493*^9, 3.4140532618039827`*^9}, {3.414053293739904*^9, 3.4140532942205954`*^9}, { 3.4140536811770115`*^9, 3.414053681297184*^9}, 3.4140538490183554`*^9, { 3.414054584325675*^9, 3.414054585226971*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"Do", "[", "\n", " ", RowBox[{ RowBox[{ RowBox[{"sol", "=", RowBox[{"NDSolve", "[", RowBox[{ RowBox[{"Join", "[", RowBox[{"EQ9", ",", "EQ3", ",", "YI3", ",", "YI9"}], "]"}], ",", RowBox[{"Table", "[", RowBox[{ RowBox[{ SubscriptBox["y", "i"], "[", "t", "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "12"}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "0.1"}], "}"}], ",", RowBox[{"MaxSteps", "\[Rule]", " ", "Infinity"}], ",", RowBox[{"AccuracyGoal", "\[Rule]", "2"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"PhiT", "=", RowBox[{"Transpose", "[", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"First", "[", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["y", "i"], "[", "t", "]"}], "/.", "sol"}], "/.", RowBox[{"t", "\[Rule]", " ", "0.1"}]}], "]"}], ",", RowBox[{"First", "[", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["y", RowBox[{"i", "+", "3"}]], "[", "t", "]"}], "/.", "sol"}], "/.", RowBox[{"t", "\[Rule]", " ", "0.1"}]}], "]"}], ",", RowBox[{"First", "[", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["y", RowBox[{"i", "+", "6"}]], "[", "t", "]"}], "/.", "sol"}], "/.", RowBox[{"t", "\[Rule]", " ", "0.1"}]}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "4", ",", "6"}], "}"}]}], "]"}], "]"}]}], ";", "\n", " ", RowBox[{"W", " ", "=", " ", RowBox[{"GramSchmidt", "[", RowBox[{"PhiT", ",", " ", RowBox[{"Normalized", " ", "->", " ", "False"}]}], "]"}]}], ";", "\n", " ", RowBox[{"norms", " ", "=", " ", RowBox[{"Map", "[", RowBox[{"Norm", ",", " ", "W"}], "]"}]}], ";", "\n", " ", RowBox[{"s", " ", "=", " ", RowBox[{"Append", "[", RowBox[{"s", ",", " ", "norms"}], "]"}]}], ";", "\n", " ", RowBox[{"PhiT", " ", "=", " ", RowBox[{"W", "/", "norms"}]}], ";", RowBox[{"YI9", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["y", RowBox[{"i", "+", "3"}]], "[", "0", "]"}], "\[Equal]", " ", RowBox[{"PhiT", "[", RowBox[{"[", RowBox[{"1", ",", "i"}], "]"}], "]"}]}], ",", RowBox[{ RowBox[{ SubscriptBox["y", RowBox[{"i", "+", "6"}]], "[", "0", "]"}], "\[Equal]", " ", RowBox[{"PhiT", "[", RowBox[{"[", RowBox[{"2", ",", "i"}], "]"}], "]"}]}], ",", RowBox[{ RowBox[{ SubscriptBox["y", RowBox[{"i", "+", "9"}]], "[", "0", "]"}], "\[Equal]", RowBox[{"PhiT", "[", RowBox[{"[", RowBox[{"3", ",", "i"}], "]"}], "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}]}], ";", RowBox[{"YI3", "=", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{ SubscriptBox["y", "1"], "[", "0", "]"}], "\[Equal]", " ", RowBox[{ SubscriptBox["y", "1"], "[", "t", "]"}]}], "/.", "sol"}], "/.", RowBox[{"t", "\[Rule]", " ", "0.1"}]}], ",", RowBox[{ RowBox[{ RowBox[{ RowBox[{ SubscriptBox["y", "2"], "[", "0", "]"}], "\[Equal]", " ", RowBox[{ SubscriptBox["y", "2"], "[", "t", "]"}]}], "/.", "sol"}], "/.", RowBox[{"t", "\[Rule]", " ", "0.1"}]}], ",", RowBox[{ RowBox[{ RowBox[{ RowBox[{ SubscriptBox["y", "3"], "[", "0", "]"}], "\[Equal]", " ", RowBox[{ SubscriptBox["y", "3"], "[", "t", "]"}]}], "/.", "sol"}], "/.", RowBox[{"t", "\[Rule]", " ", "0.1"}]}]}], " ", "}"}]}]}], ",", "\n", " ", RowBox[{"{", "K", "}"}]}], "\n", " ", "]"}], ";"}], " "}]], "Input", CellChangeTimes->{{3.414047644646907*^9, 3.4140478827192383`*^9}, { 3.4140479729990544`*^9, 3.4140479845556717`*^9}, {3.41404815408945*^9, 3.414048246322074*^9}, {3.414048281732992*^9, 3.414048285378234*^9}, { 3.4140486548895645`*^9, 3.4140486645835037`*^9}, {3.41404972034161*^9, 3.4140497808586287`*^9}, {3.4140498200149326`*^9, 3.4140498747937007`*^9}, {3.414050176607688*^9, 3.41405024484581*^9}, { 3.4140502860650797`*^9, 3.4140502892696877`*^9}, 3.414052329954048*^9, { 3.414052366526637*^9, 3.414052371623966*^9}, {3.414053188528618*^9, 3.414053191613053*^9}, {3.414053408705216*^9, 3.4140534140529056`*^9}, { 3.414053453920232*^9, 3.4140534577657614`*^9}, {3.414053585829909*^9, 3.4140535891146317`*^9}, {3.414053635741678*^9, 3.414053699232974*^9}, { 3.4140538140580845`*^9, 3.4140538165817137`*^9}, {3.414054294198493*^9, 3.4140542974932303`*^9}, {3.4140544999443407`*^9, 3.4140545610722384`*^9}, {3.4140546613464255`*^9, 3.414054671851531*^9}, { 3.4140553374886703`*^9, 3.41405535218981*^9}, {3.4140553907953215`*^9, 3.414055428309264*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"lce", "=", RowBox[{ RowBox[{"Rest", "[", RowBox[{"FoldList", "[", RowBox[{"Plus", ",", "0", ",", RowBox[{"Log", "[", "s", "]"}]}], "]"}], "]"}], "/", RowBox[{"(", RowBox[{"0.1", " ", RowBox[{"Range", "[", "K", "]"}]}], ")"}]}]}], ";"}]], "Input", CellChangeTimes->{{3.414053825334299*^9, 3.4140538402557554`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Last", "[", "lce", "]"}]], "Input", CellChangeTimes->{{3.414053328459829*^9, 3.4140533291007504`*^9}, { 3.4140538676952114`*^9, 3.414053871580798*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"1.4838109104991495`", ",", RowBox[{"-", "0.0006945567630707417`"}], ",", RowBox[{"-", "22.491562791417365`"}]}], "}"}]], "Output", CellChangeTimes->{3.4140533297216434`*^9, 3.4140538830773296`*^9, 3.4140544020936384`*^9, 3.4140546143488464`*^9, 3.4140546889060545`*^9, 3.4140554500004544`*^9, 3.414065843285264*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ListPlot", "[", RowBox[{"Transpose", "[", "lce", "]"}], "]"}]], "Input"], Cell[BoxData[ GraphicsBox[{{}, {Hue[0.67, 0.6, 0.6], PointBox[CompressedData[" 1:eJw9mXk8lOv7x+k4JUpDnLQoWSqhI3FOSnwipLKTJcvYd2bsM2PGpE2LpZW0 cErfihaUdtFGq6SE6pQ1SVS0SKqf3+s8l/7p9eme67nu97Xd9/M00y/aMXCU hITEJ0kJif//+78/H0xF4pu9s3uVKpl/QOz7kJdznkowWhplJWeml7n1Vvyn WTDCJW3fuX2MVsSnzSyt4zXvGa2Mg5YeVSb6rYyehqbOMY9usppN/9OqyPKu rNn4vz5Gq+F2tNHTv5e9ZbQGdD7Fqrnm1TF6FlLmpY+/ntLB6DkojxtYfuT+ AKPnYlWN5rP+7F+M1kE7p8DAWkjr86B8LrxCU/Ce0Xr4WSHYodhFej7WLA2v +2fJV0br48rOlWH3NAYZvQCCaS/PPXvcz2gDOLdFfgraTPaG+Gar9ehoCPH8 Bb0FwfN3nf/M6L9RtzI0YPpcsl+IW2dsJlof72W0ETK/S29u/PaR0YugPXX0 2nVZtP/FmPBj4q5tJhL4TxtDn9XzVs9BitFLcHp0ZU9v3WhGm+CETueCXg1a N8WK46nbwx9SfAD7XuX4eRU//tNiYPZxvtGFTkZLLMXmGvk0r9YhZn0pVsX/ ePx15ndm3Qzy/u+cIoqY+IjN4H1Ye/LHQYqXOW5rq2tdlqN1c0ilBFr3biKe ZRiTNmHDvSom3uJlCG8YW2UVQesWyDYbevFyE61bQNKwUbsu6AuzbonFY7Om 1PR9YtYtYfGGu+9CC8XbCot1wdokzWixFXpVNgTtO8f8XmI5Dj9Lk5HLJfvl uHiC5ZpTQvmxxvmSpC7FQSafYmvIWPxI9NhA6yvAX9zk82stPX8F3B0/2qhX 0f5WQufo2D+TDtD6Suztf7PFv53sV6HZYIH6fBlaX4Xz9o/BkiJ+G1g5junQ Z1H8bLBsnvvKRSP1YIuTvSmCo2MYf2JbvIhSSL6kS/x2w7pGLfQ98dlBpY33 UGUx8dvDenze3As7mf3AHvus9x7Zsp9+b48UJ9GC9RJM/Cvtwdpz7LzCBMqv A6YLw25a9jH1AQeUCTI6rm6menGAQCbl4KWEb4y9AzTlqg/syqb+ckSUsWAw +2+GB474UrBxMreW8u2IiZPP7zhvRP4d8dojSMx3o/g6wThaeVBrLKPhBMtm lYDn2hQPJ2hseeixh/qv0gmZly3Ur5yj+DjDTDmJ90WD7J1Rse7+3vK9lA9n +Cenl2hWM/GodIZXUYLEZHOyd8HXUzZHhb1k74InK15vfL6B4RG74J6Sdrqe BKMrXcASbtPuDaD9r8aES8aympJkvxrrv2yrv61G/Kuxnfd8+Vt/sl+NodlO 5k0j8XNF5JQZ3y1mkL0r7j6JulqkR/yuWP9Xo3PEbxS/4d9XfTbwvkr+3dCb 8lnyURnDAzdcVErsc6shfjeUyGcqBBqTvRtCL1y+fseSyaeEO0rDZv4qy2Lm AdzhcL0k9a0bzRN3aLb7L+mwZeqh0h1nXnXM7nhD9h54el4oG/Q7o+GBt0Gx f9z2YLTYA3cSGnRVh4jfA79ynP+3ciPV/xr4z9CKGFSg+lmDtxnFm64IKP5r oNvk6mI5lezXIENaK3rtKYqfJzoGQ/r6FBgNT3THF9rqr6P4e+L3RxPjjmoQ vycOmOxSbP6X4ucFxcLIdsEGir8XbrhqTSgzpfh7Yfc7paxnsoyu9EL2/Y7G 45PJ3hsDrZIf8o6RvTee6zldmzWa/HvjYppSuloG2XujZ1LDku1dVH8+uNWQ tdRtHdn74K5l437DSrL3waGC+V/rzIjfB3U+PWNyEomfjWQp2YUhOxh7VTZ6 6nVTF96j57GxKUfo8uod83s2GzuTS2VKMii+bFxf87jVqIhZz2ejvzzutEY5 7ZcNDRuV0KrjjG5m4x+28vHdqcTvC/OWB0bQJP++sC03nR48QPXoi2+l2Rsz TjPrbF+Yejb9NOggPl+UdpxOjXzN7CffFxOOfzowKox4faFuPe3LWEvm982+ iFy8yOHRE/Lvh3/1nhc4i8i/H2KVPNcHexK/HzTkbx2KVCX/fpB6V2JnkEb5 9cOyz982bMoifj+c2ZR+Zkk9+fdDf+WYM3Y2TD03++HznMDg3e5Uv/54MbRK xzyIsVf1R86ruuPPn5N/f2i9bEz4Q4fi74/S+UYsyQ7y749taWoOChSffH+M 3h1Sclie6tUfypa6N6QoP83+WKPSan7Pn/IfgIVLNzgWFRN/AO5e9e6+RfFB ACJE/luv3CD/AXi/adr7uDHUnwHQUFwzeYo+cx7kB+Cgq87um1007wNg8G1T 8YT9DG9zAM5JtWrtekX+A7Gbf3/gsC6zrhqIMqPCRYe+Uj8GDvvbxz5QRv4D sf1sdke5A+U/EN2ZjfrXQfEPhLd989tpi4k/EFINmqF7ZYk/EOKhMRV+FB+J IGhLHrCtNqD4B8Gse0XUvyPzPAj1RXob5LIo/0EYuOc/weEVxT8IoZsHnh45 Rv6DENP2KPcEKP9B6HHel3DlX/IfhFOyAyu0ZpH/YOTVW9o7a5H/YEz2CKoZ /Yb4g7FtB6oyTjPPYwfj7sVnv5K8qP+C0T9oP6lDSP6DMVVBcmyOIfEHw9H8 9ZEbP6j/glEV5DKvkupLIgQrvLdGRpyk/IdALnvIcMMo8h+Cghk+E62uU/xD 0GD26Y0hl/wPa+PPPz/Nof4LweXeftnAMPIfgsiznl/YpsQfAsn8x4Hc/cQf ik4Nra+/zhJ/KOzy/zwaS/lAKJbUOq/s/kXxD0Wl8Qt2+FrKfyieLui6lu1I /kOx8/rEW3sLKf6hmGpifjifRfUXio8zrv8dlUD+w2D7znrfKVfyH4bchb+0 5eaS/zDYPHyp4vuF/IehMFl/Ys9ryn8YvGOtygzZFP8wvPMt/VamR/7DUBvi opwylvo/DFKlQw/aflL/hyMh29tUZaT+wzE+t6QpiuoJw7pq1LQ94RT/cLTK vjnC3kf84XDvnfgtyIf8h0Ndal1tpQ3FPxy53D3ylUco/uEISA7w+lFP/BHQ 3Fx88C8e8Udg1Q87pWpN8h+BP+cuXfJ2C/mPgENOwedWV8p/BAIW23695cDw 5Udg+qTzLu2+zDyojMD6vM+BJc7M/aA5AgVHwz3aFOn9IRJPGw/62nUw9qqR ONO7L/KnMd0HImG4TDJeypXR7Egcn5Epn6JL8ycSEZKJAR9lyH8kUk6nNga9 ovhHIjgsy9OV7g/NkZA+XTsWTyn+UYhI0LpTWETxj8LaY/bc4vV0n4hCunLY kZZ86r8o3DtrYmOfTPxRGJeqn7PfhuovCil1lo231Mh/FC4dmXZ0tDn5j4JC 2bgPW1XIfzRYig+jVD9S/KNhUajGy6B5h2hcjn3n+w/ND3Y04idNtztJ6+Jo 3NTp/Oc13Qfzo/Hp9gLntWso/9GYEfawfB2dL83R0JeOfcGnfpTgYEu7kvUh Oh9ZHCw6OGWwNYD2w0EcXyATa8JoPQ7WcY1+vpej/XFQVP95firNV3sOpvaf vHsohvbLgWOIf+VyOj84HGj5NP/PdRvtn4Pfj2lUn6Z5mMWBM49be3UF8XBg fiFeJbuT0cUcrPZu6/Two/hywDVpqSpXZnQtB3bvrX9T3Ui8HIS5WH1ptGX0 Bw7eDH2ZJjInfi7K5qS1OSoRPxdGbesWcai/VbkQf82/2TeZ+LlIvjnjQ3oB 8XNx2GryYIwn49+eC9792C3SX6leuLBKW5yqXspoDhdyts07P30nfi5m+vuF Z1I/Z3FRWyG2D7pL/Fzw0z+Kb50g/mH/Q9dMo1Mov1xIm2XsSfVidC0XE33v HnqjS/xc/Phe7nNbnfi56P3Ur3BHhvhjII7kGxxYRvwxuL3hY2bVRcp/DNS7 D2jZU33qxWDZpYCq3nzij8HAa79oQyvKfwxksnk5z2iesGMwNqBbT1eO+GMg WD/1Vvph6p8YWPOuXEgXMzorBr9f5vEUW4g/BkFv0mriqF6KY2Ax0wyTuMQf AzOX7vGH6fyrjYHStulh7e7EH4Nt35bny3sTfwy+S2//02Yv8cfil77PqZKt xB+LbW0K1auoPlRj0bT9d9WoRcQfiyNT3i7u+R/xx2KBzlzTc3S/to+FWovY aLCA8h+Ld1UPfh5ZQfyx8B3smn+KzlNxLE6qj5sdKqD8xyInf+i5J83b/Fj4 j5ZrcOUTfyxitlx8aedC/LG43qBgWhBP/LFYPs5NLPmE+If349g9e+Al8cfi vlWycQ31h0QcTuQ/btLwJP441N9qaThC9aYaB5hc+sRdyOxfLw7h++ctUxuZ j3HoKzHbNkaD+ONgoDTUnVhM+Y8D55nt1kzKD2f4eQ+atK/S/VMch2qJNFMh 5SsrDulGzx03U7/mx0Fzz46shJH8xyGgaK7o+hfij0PXxTnn0qWp/+PgdXBt 4RLy3xyHWXeuZWyl950PcVCSTelpOUP88WhcP1Svo83Ys+JRwDXza1tG50E8 JvueW2z0gvIfj07fF4OPd1H+42E6xzx6I81f+3g4tbU+eLmB+OPxhHvDS8GN +OPRMjfEcLc98Q+vd/dviaJ5khWPpwUqShZz6TyJx65WmSf3TjC6OB451d6/ vd1F8y8em2V8R6lQf9XGQ5fz5MqxCuKPxw7dao7oJPHHo3BBy5aWK8SfADX5 k3MP0v5ZCXC1qxo7PpPyn4Akxa1HN14j/gT82aw3y42+1yEBl9QkTsukEX8C xOFX/zGLIv4E7GCdqFlH76OcBIwxiewbv5z6PwEff0/IbH5J/Z+AYIOCO5Ys Os8TwL9g8H0L8RYnQP7E9p6CKcSfgBuRi1ob6H5Xm4B241spfnS+NifgpHhJ f+ZP4k9A2eserWuPiD8Rs/bYBVw6RPyJyFxTv1/3JvEnIqYizms15UcvEaKF cpw/RvKfiGl3b2dJUb/YJ8JN1t3Pi+YDOxErds7JzqH+5STiwEBe7kU6j8WJ WLbkjIlECeU/EeMPh3eWz6D8J0L1TeT9n1bEn4gh+/xwp0Gq/0Tkm3V8lCf/ tYkIPH7p3e6R+k/EIj2U1NE8/ZCItQ6vpMfS/JFIwoeX0d2RNH9YSfC3MpQf t5v4k3BDyrDnYx/xJ+HRdL/cG6up/5Mg2eZbX2xC/Z+EcYVFs8IaiD8JzhpP W2JoP5zh509zctGm+744CVOmG/pEkb+sJAylt3TqRFL/JyFgpfuEqdRPxUkI /vfxVE96f6xMwhy50W0H6H2oNglHlhcfONVE+U9Ce9j8n9fuMfpDEiJkz5q/ pHkswUPzXznLC74RPw+b4lK3G46n/uch/NxEF9ZvNP940B7v4PK4nfLPw2pl jb8m0v3Gnge9rif3p9P5wuah62qXceZh4ufhR12Pqv/I/WfYn0iwSoHO0ywe BH0CQ/VW4ufhbP3VcZepfot54EQ/Srek+0ElDxZTx83ctYfyz0NpX7dzOa03 8yD/WZibZkz8PORF9brY3yB+PpyvN4QXT2f6jcXHbtaLTodHxM/HX4tk/rlo S/x8bPryh2K4JuWfj5n58bOuLKX881G1+89xq2bT+cfHMZ1crjp9z+bwobNI XdTVRPx82G81fzhpiPj5MDUpsIgYOf/5cGHZmWrdJn4+rs75cTTtHPHzoRAr PXUr1U8tHysnVWiMvK8089Gzum5yw2aqfz4uOT1Q6L1M9S/ARcUFdUY/KP8C bJ2VZWfWQ/UvwN1716MUaX7qCbC2rea3gZH7nwBmOrOqPtN5ZC/Aj6861eaR xC/A+tDv1awa4hfA4uq19+EHaf4J4DO2dIZIi+afAA/yXlRsJf/5AmSc/aFV Jk/9L4Cab1uWA90fKwVQydBMK31O/AKId2qO31dO/AJUz3onF5FL/AKUVb33 Ksoj/mR0iFRW59H3K1YylqwY0L5J30dVk8Fudzd1G+n/ZBzmS8/ePYryn4zl XuY8k3vEnwzjMKe9s6ne2cPP67990bGR6j8Z1zwW5urTfUmcjDWdnA0hV4g/ GXkhykYVgzT/khGtbVteeIj4k8Fd1qf82IDmfzL2KjxO20PzsjYZzdYppUXz qP+TcTYkXNGE3g8+JMPQ9WbvQ8qXhBAz4rLdNel9mCWE2v0FvokfiF+IfTde p1bSPNMT4pRMxY4B6ncIoRBY9eb8BeIXImOhXZLrAeIX4nbhH12J1O8cIRz3 tPf9TeeHWIhH30qOFdyh+hdi50Ov4mL6/4J8Id45Sp/vGnn/EcLBzML+cyHl Xzh8n91WHEXnR+3w82X/7kikfmsWYpee8OtV+l7zQYgz+uvi156h/hchXrqm tSma7j8ijJ4gm/tqZP6J4Jo92yaO5p2eCPm2Qk9Dum9BBM2L+5/l9RK/CK+a SmsmPSR+ESr2PG/IPU/8ItxVFW6Oo/koFmFMWYhl/cj8F+FXcB67ne4f+SLI 5VRHbaT3kWIR5sezFYdG/j9FhLcrHj6Woe//tSLEZL+dcIbmS7MI572njpI6 TfUvQut4Pcneke8fKeg37FR4P3L/S8GeeUE8Sfo+q5oC90m9VgF039BLQUDB hT4O3a+RgtVF7NTi4X77P82KTvA= "]]}, {Hue[0.9060679774997897, 0.6, 0.6], PointBox[CompressedData[" 1:eJw1WXlcjPsXHnuWyBIiGSolJFSIeCyFFk0rqZj2bappZpppaqZGhGRJWbpt d2xdtOhaUkiTimwJ1y4mJMmStbjkd3+fzusfn6fTec/3OcvzPe/bBL8ol8Ce LBbrSw8W6///d/9rX8iRjVv8p8HHSvoBUk1GlKTN+E5YA101VVVGs2sJayH2 sGmKT03Hwm48AmZ3PXqNtnhLeDTczRfp/VovJ6yLkHFtNc0PH5M/G5s99+c4 +d0hPBFjP91827z5T8IGqCno6VD1vYzwJHTIItzcs5jfN8ajT64XQk4/J2yC +tK9Jp431YSnolr1dvSt1KeETaFdtT8vL/E+YTN8G37R6aSqmvAMNA3K7ME9 fJHwTJyxcNo8TTOJ8CyYhj/rIbiSSdgcy8szZDoaZwhbYONDR22tZ38TtoT7 pVEljm5pxH823ttUpNxeVkl4DjaEpDbrLTpAeC4G3klpL+linmeF9n0zRr95 eoPwPDQJcjeW1DHnnw9bB7aFZ8kzwtb4/PfEA2w5w38BvD/43TC0ZPBC7Lj3 tHDKxieEgTa7DcMX337QjRXAiM+uliEXCLMWoeBovFea8yOyL0JykUlQ2fm7 ZF+MN+ZZWTvCGsi+GLqNvboEZwmzloAvGLz79t3bZF8CT+t9J/cFMHyW4mMe L/HHkQtkX4ob5eNTHnVUkd0GQ3gWh181XCK7Db4unnD/YzJjt8Xoue+bXxZm dudPYQv/5rl1GUEeZF+GN0vvsmSf8sl/GaZZvTKSFB4n+3IYuAXWZjjvJv/l cE+uCjriVkX1WIHfBYssUgtLyb4C1w37zko8mED+dtD1/mJ/o4Gep7BDxdGK 6I2Cw2S3x1xz4eZHqTvIbo/qO3cs+C5M/zhge9/dluedS8jugIu/O5fYCVRk d8SVkJmhZjnEV+GIzLK4o/WLGftK5F7IHXPizTmyr8Rl0/K7W1T7yO6E5sBv HQK7Mjq/E9ZaN9w+ZVBL/Djoa2m9rKbuZDcGB54eYz/OPMKclwMVR8YamX6y G6s4mGi1b9LM0Uz+nWHTFDE4qZX0AM648kXHLimY5knhjNH6eRcdOIRVzvia dXDJwW+l5O8CReJRhe0o6ne4wNPs/UqREdkVLihI8njx6BzxU7lgov+repuU Q+TvimLJUHbnAiWd3xVVyTwZ22878XWFzgeJJO5cNvm7QqP9y4aY9r3k7wb1 rS8/EvNSyN8NTpIcnaj4E+TvBuemkl4m9fndWOWGvZHBtV8dtpC/O55qPR9m 6lVE53eH94lXWxZMoHwp3DFR0Ht5TuUxiu+OP4t1ZihnFJK/B6xfjdGfnqEk fw80Gf8a+aOF9E/hgZx3tiUeGTnk74HGsJlaKZnJ5L8KwvKZLXsEe8h/Fa5O ObHJ5usu8l+FxrkVSdWd0eS/Cqebq4IKAw5R/VcjaYf77bf8fcR/NSK1Z+tv 000h/9WQJ3bYx587Qv6rccVt67rBPKZ+nngSpNg9eRD1Izzh9Whu7hU90k+F J14OLDM2O0j9ovJE67e14uTlTP+uwd89fSOPniwn/zXgy8o2dWacJ/81UHe+ 2vSmsIL812BhuthlwSrqB5YX3qO6tGY7zR+88OfVzEu7Ppwgfy9E3ArQ/Hj2 FPl7oaXQb2qNEdWH5Y0cV00ri01p5O+NzkOxqWOMs8jfG9eDeoodX+8nf28U CXd92+HMzLcP2jJQGqwVRvnzwcWLvxvv6R2h/vGBRM/04b7SXOofH1Te2HLG 6zGjH2sx8d4O8fCXuyn+Wlyauac8bGswxV8LR+71/RqLqT9Va5E22UvJP8Ch +q3Ddem9cvce1B9Yh0Tf+gFjfxaQ/zq0KFJjz2/7i86/DtYp/9gnRVM9WVxY vi66ONLqYDdmczHF6uS+kSzSJ3CRerLs/qMbxJ/LxcIy46QuOZMfLkoLHDgz SwgruYjZuGKl/UBm3rhQBgYM231pdTdWc1F/Zuuejk9/0fl9EX3NPz/f4Vg3 Zvvizo++e06tT6V8+iIkdZ799T8pP1xfOC0J1dt+m/gofDHSS/vsawuaJ6Uv Qt+axWb1YPrVF6vcj35xaCL9V/vC53fkWLOXqcTfDzzfwT9UjXQeth/YrubG RX1LKL4fPvzMOZp07Y9uzPXDzYrWmphtGyi+H6q7htVdn3aA4vvBbvzq4cWB pAcqPyQdW8TOvnSU4vuhTZm/ZuhTqg/LH+m6P39v70X9xPZHVmlYzumPWym+ P9x/afRo6UP7A9cfGofHW8x0WEv95Y9kHc1hk+M2UXx/DGmcLj9zayvF90dx PWuTYirphdoflUkrbvpuXk/5DwBrwai6651fFnTHD8CkrZGpjYO2Uf0DMEq/ 7+4Af6b+ARgmerWixYT4KQIgHmtj0ziP5k0ZgHkdV+2HZdF9pgrAyJpblzVt CasD8P3NDk/UUL5YgXi2Ibh5aybpOTsQvDXrI3uuIIxAOD++MWavF9WTG4jB XX+vlp/ZSfEDMc5kU0jEpuRuPspAPHIysXc6zKN5CcTmRujsSaf7TB2IOaOf DjjisZniB2H8CB0dUXIi1T8IsDZdpTP+T8p/EGQNLn89jWfyH4SoGH3zQe83 UvwglCiT4q1a8oh/EPrHpT78I430XBWE2oG2H3xyaT7UQSh69Y+OmfsfFD8Y 5h5TrF9eI71mB+PFxg2tjpn0+wjGbMPKydZj6HncYAROM/vjRhfd74pgNH/Y ecvrCT1PGYyB+rGxD+V0PlUw2qoHNd+ZQHzUwdCryV/SYnyY6h8C5C7b4K5P /NghsIy65yD4RfsuQpA+q6Tfu2SaZ24IojJNHrVdovMoQrDCRvDVSYOJH4Lm CTfErU2EVSFoVQtHRF2hflKHoMbheymrUkz8QzHRTDf+XRTNOzsUhz/0u3ui bhPlPxS5lpflQ7YLKX4o0jsGe/1+mE7xQ2Hh+6JdHErnUYaideC4KP4Nyp8q FMXnCnrcP0/6qA7FxUVV/5jeZ/ajMBiJD2W0c9MudMcPw60Bxf9+ek77OMKQ Y1k0JSGBMDcMJ8SlczbZ0PuMIgzJuXk3nG2pv5RheLt+bJVuPVP/MAwd8yTA fSTplToMNx/pqgs1SG9Z4VgvG5hdPJXp/3Doe3lNkfvQfoJwGP/W+VW6hvSf G44l+lfiZ7XFUPxwLB/Y96/Yc+spfjj+nRT0ar4DzbsqHD0HK80zCkiv1OGY HT0u945LENWfh426bY9n2Uso/zy4Rq2/79YrjOLzcETtZTp4CPUnl4eo75/m bTEjPVDwMPlT/ulpo2i/UfKw043/xSaY3n9UPPhMbZz2i0P3r5qH5RPchqYX MPdvBE5b7+6y/k76wY6A12uvWaIQ8kcEIrMGRT2uJzs3ApWRdueKfZn7PQIN jQ1axxcy+hMBtyfzOX7rGP2NgPJYQfGoo7QfqCNg1NbBu9Cf9IgViT5LZo+4 kFlM8SOhbrka3tOI6oVI2M24eub+fNJvbiTMl2l28UuIryIS1+qKLMylpE/K SAzgzatbPJf2N1UkJp2p3Ni7k/KljoTk0vgvAabM/hqFo7ecl/sZkh6yo3Cx bCW/pYn2S0ShatABeesnJv9RaJEYrj2Rysx/FOTTpB8C/yKsjIKb/Ovrn5l+ pH9RYPW3yi6vjqP5j4L9vg1ZmduXUHw+in8O3KZhQPG0+Li6/3sTy1ZO5+Fj 3s/WZUun0T5txod2dD4LzRtoPvg4Z+Tv+jON9ksOH59aXvTfH0/zyeXD916E YI4tYT4fWtLs3GdmdF8p+NBnTYvM4lD8ND7C2wd21oYSXyUfT9aPnZJURPks 4eP41EfuRVeZfYyPncmzKnI4NG8NfDyt/rTzTVYu5ZsP63EN2nWgeWjn47lg pV7x0Z3U/9H41dTecopL860VjY+tF+eMatpF8xANf4fbboF+pFdm0Tjb46lz +WVmX4vG1Yyc7EOBdB5ONG429O9dl0z15UbDumd92cJBNO/8aDTpVuSpe1G9 FdFICbjPS9xG+pEWjUaTTbuDCrcT/2iMmSs2vMCn/bQkGiHeiMMU0gdVNGxq J9ev0kgk/tFocb2p9SNyD9U7Gscra59Uxu7txu3RSDuyXrp7OukHS4AE08QR ekLSey0BZhkZF68KonqxBejpPN3a+C3tE2YCrDj6q62mhN4PIIDWAWOD2mmE OQKMXv2mz5BhdH6uAANXDC9b+oz2Sb4AjccWnTKsZt5vBOjvxj018TLpdZoA W95dtjMMp/tGKcBHCzv2lt4Uv0QAl9FRth/S6b5WCdD3QXv7tOmkZw0CnOid M25QBOmhWoB/L2VE8C4Sv3YBXjoe98m7KKX+FyK0/FG0a6eC6i/EobhMl72H BFR/IdwEkz+mTYol/kK4fMs2NjpL+YEQMysn3XzLIj4cIbpmvL5W0Zu5r4Wo NZj+R/NVsvOFaGgVaWXFEz+FECmnLhn02RlB/IX4NdyKty1X1h1fKYTNS71+ NxO8unGJELNH6HQ5XKN6q4TY0zJj/Y/7pO8NQowaasArfkz1UAvx7MTD0KQ5 ocRfCLavt0FEF3P/ivC+dP7l5FCm/iJUzffgW7hRP7JFyOxccsdeSPUzE6Gz Z6CF6DxzP4ng9D3INnIB6Q9HhDH3jgecM2T2dRH2pRX0nN9k0n1+vghjzSb0 YUcw31tEeHZzTsy/tlu6cZoIyS32/r9rAqj+IuSbGC5O76D5LREht+3J9tLp pA8qETy+593njmL4i/DDpj7z3Xna79QifNE8uHN4E+lNuwgpQu0BFo2UH1YM un4Oe8z7SPXUikFDncI09wfpETsGVwtmmDnepn42i8EpjzY921Km/2OgN+Ls h+ghbsQ/Bvb5P0Qe/dLpvv7veRHvbI6k0n7Fj4E+52Cc0yHKhyIGU+P0b6fe puelxeDc3+uTzHwo38oYvPvalmi0gtG/GKTo1hR+GsC8X8XAV2p6l6sg3BCD BQW6l28uYfatGDhml+tXjw6m+Y9BgUGfrlPfOcRfjK3pTVXmLNoXtcQwVGnX 1H6m+WGL4fP5+NbUgzQvZmJUPzCM3NQqIv0X48qJH9umy6K6MUeMe7N0ix58 pvuAK4b8Z6A1R0j544tx5JfLiNAu2g8VYsy0ZXes/kj6mSZG1eS72S3FzD4n hrRZL/Z7JdlLxBCmT/L50MbwF6Ng7d7SM7syiL8YEzZaSUv+pf5Ti7HJfJ6z 8VbSl3YxNmg883p+hurNkuCC18JoPVEIzb8EqWaeKS3XRtL8SxA2ViDwO8rw l6DuZZHOKl1mP5ZgfuQO95btzP0nQYWTF79fZRLVXwJ74/U5f6QS5ktwaM7g /pIKR+IvwcsszbP+Zcz9J8F7zUy+60z6nqOUYPH45a/LTlO/l0ig99G7+UIH nV8lgfdWpyWzpTSvDRL85ua0RzrQedUSvL5UeDVwqTnxl8DKoPw3T07nYcWC PXpIfZ1FIPGPheBR+LI7jsz9H4tfQf+e1/Gg55vF4kP/NovxOyjfiIVF7TF9 19FUX04sNF1OGO3eS+9n3FhYjm97+rCW8sWP/U/vDN9wx1K/KGKxXHBw+m0d 0oO0WGR/Xdy+/+Rm0r9YTIsPs5ypGUn8YzH2fP6/6EXxVLHQuT9/+NZdNC8N sfhacMDWhkvzoo6F3E7DLuEL9Ut7LPSP6Tas8KP3AZYUkpiPgzR1Kb9aUsx6 1vF1ch59H2D/h3NMzgxi5slMiugBrbcmSxj+UnCzQ3iHdpH+cKSQ6fRZOfUz nZ8rhdGY6dtaufQ+w5civd/j2ztuuBB/KRaIOMW971G90qSwV0vvaXfSPCql uBCTWSkMExF/KcIy9Kx6iZjveVKkuv/zXHqd/BukmPbMNvJTGPWTWooZJ1tZ 1aepn9qlSB4bUZg0ivn7ShzYmj69uPtp3rTikHpkkJh3jjA7DkHXDK6YfqJ9 yiwOtiOKNW46MvtPHD4rvfO/ZzP6H4eEkMXL+wyheePGYUuPofN7zyc7Pw7l JeayuU10XkUcRnUqii4Mpv5Ii4O+y6tVo5Sk38o4DHhTMG9ZG3P/x6HXwYof FUJ/4h+HGy/kLo9u76D7Pw4dC9YZVC2j7wHqOFjY+l079XQV8Y/Dnuq3f7NM mL/nxKNy7vdPO/pTPK14VI8wPvU0m0/846HbJ7yH06Q44h+PMW+1eXZhpI+I B1znbnM0pP2ME48Fni+0ouaQvnHjUbJgtU/9SMofPx4JiU1WNbnM/R+P+QvD RBMa6PfT4nFsbkW920+qnzIeTx31LbM3Uf+VxGPL9En97GfTPKriYXipb9lu GfVvQzyWv/5Vmv+K7nd1PNIWPtmXzqX36/Z4HDbaPb25kfZ3lgwFsx9ulb5d RvxlOBfou1FZJyP+MgweUHfAagrdn2Yy7NPIT/3xjvoRMlxpW6XQLqf9iCPD jq0TN9tPoO+bXNl//ZaRJh5JfPgyzL1ZudnSgfpHIUMq/8JewWqaxzQZ+COX 77LSIayUobomJPe3A+WvRIaThl8+z9Zmvi/IUMc7VtV1i+Evw4wBdr3WKWk+ 1TIMPzm1Lncx9V+7DF75xTXKFEb/5Vh6qOgux4BH/OUosLqmWfOB+pMtR+HJ Yt8e72gezeR4e+NAUbiU9hHI4ZUw5tvSs4uIvxypL7dteZNH+ebK0XT4uufX mfS9jS/HX88361QNCiT+cvRetFJjxjem/+U4bi6adHM0PV8px8+dMfoBjsz8 /3eer5LJOg+pH1VybI93GpD5jPqxQY4XD1n9sh2oH9Vy7L0+ipPnQfrYLsdY y7PfL2sy358S8LayzDpdQVgrASu6zs5f00p6wU5A7YMHGw7z6flmCdjc0/9u ViPNCxIwd+h1Sw8J9ScnActMPTI8a0k/uQnI9zPuzFu0luqfgG+NhxelnqF5 VSSg/1PPC/FFNK9pCRDWhB7MUM4i/glYpD4dNMmZ7qOSBMw4b512cQj1kyoB DRO9decFk70hAVXRwwa+MWb2/wSIVL7av5YriH8CnrmuzlyQy8x/Ing/Tvzz oIPOr5WI4949s09vpn2MnYilE2Iap65k3n8S8UJ93Vp7B8M/EcF16XY+F20r /wcXcU++ "]]}, {Hue[0.14213595499957954`, 0.6, 0.6], PointBox[CompressedData[" 1:eJxFWXlczOv3T3aiseVKNFlKKFkqUXqHKEubSmmb9r1ma6ZZG5G9a2zZmeyy hSxlG+ulLKGILqWufauQb5br1+/1dZ6vf7yO85znLO/3Oc/5DIuodL9YQwMD g89tDAz+/+///mlw9Txco75aOVL/+x/AmzBJ5j5r9G+5E5xKTqx6XjP2t8yB 8OSVxf7z7H/LvbFon/OB53Xjfst/YEfZ4wJ5OslmGG5nP3gfh85zEdnH9dLh JyQPwpzI+d7nS0keAs33z6kPR5NsCY7/T+O7YSQPgyrGxWDnV5KH4/2n2Hmp Yoff8kg4r5FxPnQn2Ra9x5c+fr6NztthkanL15tVJI/GiMrKUde96fwYuPbN CzqZT/JYDC4usVo6nuRxeGeTYJ3/guzt0WIrWtnDhPQOGN3AaxsaSbIj+vYc 0KHCj+TxSGjbWBTwi+ydEOvV6VVmJ9JPgLOq79Y8L5In4pv1rRdfCkh2xuKq +oLCbo6/ZReMmLrE2DqY5EmQfBgde6A9ya5I9BaUbcole2Dyvl490xS/ZQ0w b9InS4ma9G6o+ilwLI4nvRs2rtnY6WUY6SdjoNZr4KNk0k/GdEFc7aV00k/B 04brqrQg0k9BjfTKXqEP6afirGXJi9U80k+Fu+nJfQ9jSe+Os1aLdp6eTXp3 9J3WdxN3Gumn4VPdz2Qrdv80XI8pP3iE7jOYjqsv+lZVTif9dJh/K3KdMIn0 Hlgc4fk+xpf0HnjitMlwP7P3hE3VH3mdPUnvCZ1JpZ2hE+ln4MyUXc9uziD9 DKT1XqGLjyH9TDzsZyvewuKbiV8VdunHJpB+FmyvGNddnUz6Wcj6MfzLfcaX 2Yi0fhCSriH9bLSc0ub9QfU28EKBzcnZQ1l+Xlja9nORnPHXG9tnvLz/PJH0 3ujR5NcSyie9DzL/9TrgRefhgxGjK2Q2HnTeB+Y/j3avCP8t633g+67Jz2U+ 2fti/9DtI9/sJntf+LtVR3fdRfa+6J0sCvwmJHtfzO9kX1jI6uOHR2uueu5O I3s/jDIUzxrE4vWDxCl5c6d5ZO+H2Gb+CiHj3xxEJCV1/JpE9nMw7Pu6EYUJ ZD8H708GJZ8JJPs5eNxm15eWuWTvj5cOJV6Pid/wxxffa20TKR6NPy6oFsb1 Yf79URc33dt/DtkHYI9oV0tsFNkHoGT53PYaKdkHIEAf+eQBzSN9AL40bf8e SvcZBOKeYIHmZTDZB0Ix4LnOL4XsA/HCdm2PQKq3PhAd6hrLSjPJfi68Cy7H byZ8MBem9j24uayf5kLXUpA+WED2c4G+3jUjmX0Qru13kbhFk30QRtwcsPBC CNkHIWXFxXAV4aEPQuZud2974qNBMFKO3MyatJnsg3F6Sq7xqeNkH4x/L+9L fbGT7INhNdQRLsx+HkpG+Zq/pHgwD10thjTEsXk0D/V7n+0MZvWbh95RHs/d 2XwJwZ7jRga7RWTfKsc8zOvF6h+CrvNOnB9P5/UhMPi8a2sT5WMQivk9eCVT SI9Q9OxypuUT3acJRbht0Ohc6jd9KOROgweFsfkUhk5hETv8SI8w9A1p7j+D +K4Jw/e6IRahxE99GAqedIq4y/o7HKd3Bp0VM/6FI/PEtc+PCS9NONyL+oeY Eh/14aj/Rxn1nM2nCEQ3Rf/rxPCLgEmM+5Hu1N+aCBw1OvE1W072ERgYuH6U fyrZ85BnPLlqLcXD5SEksVp2lvoHPLz02aP/RPXk8bDm2MWW+IV0Pw+3mj/b mEt+yzoeyty+2A+ifPQ8TH9YM8yL+FzLg83mQ/940HmDSGSfv7NtFeXLjcSg ZYuWdaD6IhJGxim5f9N9vEgkhFzcJiE+aCIxfqS2c3I2+Y+E4px9x41K8h+J 8FMXmt5S/WsjMTy88/pSNj+jsF1Xv6sqg/xHocu74uLblC+i0Hi32kRA9rwo vNOtPNGGzZcofBvptLSW4tdFYc4+G+d1xG99FATgz36xlPxHwT35xK+tMvIf jYXrVuvvER7caHRbqzLjMD5Hoz5kjRvrd140br/uOvwEwzcaeyZLurSj+HTR OLPMMi6f8T0aL7RPlF/JX200Sr/eTfAjPhjE4EHWko/vCB9uDHAqnbeH8TkG lpd8b2VQfXgxsBd2ca9eQP5jcHx5Qcj71eQ/Bl75i/u115L/GIROv+Qso36u jUHzmP7GeQz/WJysfCrYQPFwY/HjZLd2Yaz+sZg700lqTfXlxWJ46KG5Fqw/ YzFWNyn1MtVHF4ufG4/UaOl+fSw4W7plV1H/1MbC2kp1YSR7/+Iw2ux4Py/y x43DA9GIP/RsHsVh5/GE+1yqJy8OV+P5Rn6ElyYOwugs3CQ+6uIQGZy29KaK /MfBNrfqsCOrfxx66SqLOBSPQTz+nbWzYRDFw42Hperj/Go6j3h4/bnvTQv1 Gy8eZQ0XltTRe6CJh3bjv8MzqB66eCzZ6mX1gL238Yi+VT+5mPmPx5V6E10U wz8Brr92Sd4y/wk4sk06+w6bhwn4nH71rSfr/wRYrbLceIntIwmYMmpJYTDt i7oELL9setKB+U+Ae0FV52iSaxMwNabXvBjmPxELdQt33aP7uIkoeDiiwyhW /0RUPzHad4fi4yXifKFNhJj4qElE3ZrUlq5Z5L9Vvz4gqc8i8p+Ino979vzC +JcIqf+uqiVUL4MklAdVuK5l+CdhS1YTHCkfJOGGzGTVRYqXl4RmpXNWGJvv rXrbsIF+DP8kjLrwaMcZNn+SMH9A4cORi8l/EjbvdZ0oWEL+k1GxyW7qOcb/ ZDxI6df4gfiLZDxy1Ux8Q/HzknG8OvyoGfFLk4yCm6L+F6geumTcnrFVacPe y2SYmO3/WMfwT8aqJeMGmbH3NwUmDhWLutP93BQEiFZztWz+pKA0WrQlheLh pcCC8yL0M9VHk4KEV3ddnJeR/xQY+rZ7+yyP/KdAsOWJZBS9/7UpMDvSzrJy O/lPhfvM/H2r6HuEmwpjD3nJNuI7UtFSml1azfifisMXPUSHad5rUpGSY56x k/LRpWKAcOBdBduXUnEp56CbNZ2vTYXbFr1qJfHFIA29nc7O2k715KZhk8dH izcM/zT4tcuu2U3neWnYuv6MYz7hq0nDpamB024Rfro0cIvzp2gZ/mnIK7DO S6L7a9PwxuremM9s/qcjpFR7czTxnZuO6YPXrxbQfUjH0ezv9/MY/ungBP2T dpDtN+lIdBObXGH4p0M2//JlNd2nT4f0X3m3WXRfbTqUr4OTDdn85SMo4XHC KsKbw0f4xmLzvaTn8rHjk2PsFYrXjo/CbTe8NpM/8OGYP7riOs1nHz6yLo45 W0T84/ExxfJZegXZ8/n40B0X+lN8Gj6Sp2bVy8ley0fRuDvvEqi/dHwc2P23 yx46X8hHsWLIzGjCQ8/Hqz5G97OpPuV89OvXeZ8p5VPLR68FlaujKZ4GPmw9 tfu92fwR4MoGm+4MT44Aq94eT2+m81wBCp5VC2ZRPHYCcPnhtWPY+ySAxchB hXYUj48Ar9sp5/lTv/MESMquP3qM3ku+AI+OdKyvYPkLcDS9W8Qvuk8rwHbf dinGlI9OgPaTbu1g73+hAJ3vLr47n/W3ABLu8eu/CI9yAaL6TE4OYngL8LXt j0PNVP8GAdp8nDulL/t9RIh8P354D5I5QhzUrbU3pfNcIXL2lfepIv7aCdGU dGp5M5vPQnjJzRvjqN4+QvAHfAvoxeaFEBGGnoaXqb58Idbb58f8b38U4k12 WZKI6qMVoqTU6bk3+dcJIT0w13c71adQCEfRt8fDqD56IYYaxnTeT/mWC9He Lcd1KdWjVog9LlLzvSQ3CFExcvPKfgx/EbhG86a1p/w4Ihj3U1zPYe9Bq3xT 3TSN6mMngtsml3PPGP4inPyjsCSX5oWPCFMv/TSVsnklwoODpVdrSc9vPW+a nihi/StCXPfrgiaqp1aE8T6xVtupXjoRynaZRDtTvIUiyE9bu6goH70I1f4b 9hrSfeUi/Hlnw+HhZF8rwjKzYSl+xM8GEfqld/a1YvNHDBtVbe4Wuo8jRkap fYwlw1+M3pV9PFaTvZ0YRt9a2l1k81mMTf0Nzbay/MXY2FQvHsb2VTHOXbv2 qQfFzxej5TO3ail7v8TI2xiRy+anVgy/BuFa9t7oxEj38D7myPgvxpnAoi5+ bL8S47yyU24Yzf9yMTbbF1e7Mf6LgbE+qjt0X4MYJW4P88cx/DMw6NmBqBFk z8nAkCFXCiOJj9wMhLWboZhA9bTLwJRZpjx3xv8MTOxvtiiI+OiTAc9us5/3 Ze9FBgYe+HisivT8DLxese6rB8M/AyleaokL1VubgY7XTSvPsX0qA4uGZwwe nUP5Z2DSrZbpF6j++gycLDAcM4biK2+Nx2Hi0nyWfwaOVf2Ze4Xq1dBqb7l1 8TmSDSSoW7lXY03+ORIMfbiBV8PeIwnMX1gUHWT5S7DjrFPyG8oPEjzpUL/k MMtfAmXd0V7/+16U4EeDZ59chr8EZeHLBo1j+6MEkTkd1wyneamVYPEow5c1 bJ+QoPOhxj6b6P0ulMAkv6ZfDXtfJVgq3F4bQvUql+BQiO7P2yTXSmB8cVUk +/5okMDe9fKAbQx/Ka79SBZasP6XwvlIm2Q12XOluPyRm1NN5+2ksJ9bFDyc 7edS3DFZvqGC5S/FhkZ+zyay50kR1cngL9aPfCm6HmvmhbH5L0X3CJ/rwxn/ paiKaCyvY/hLcc5h+dMZLH8p1nkcKDdh/Jdi4uLbXHOKp1yKtoGB4fls35Wi xa9s0ASqd4MUs98M/rGW7X+ZMHd9pFxI/jmZKDu3Rcm+T7mZiHU+/yyHzf/W 81v6+B+heJCJlhLBsiSSfTJbv2fG5Oxk+0omNLwc3U2GfyamqeN+Mbw0mcji NjYvoHi0mTixf8SXd+z9y8SSj78csimewkxYmYbsUjD8W/ULuh2vpn2vPBN7 Gst/vGf7diacbzib/U1yQyZsp5rqj7D9T4Zzs7c7pDL8ZTBKDs9QkZ4rw6sB 6QGXKV47GaZOsE7Iofggw4Piv8KrGf4yrCgc37WQ4S9D2ardm4bSfXwZPjSt rrBhv2/JYFMZt0TF8Jchejb/760Uj06GwdZNVdeoPoUyvKhb0PYQw1/Wyu+9 7d8w/GUI3CN+XcLwl6FZ16NIzfCX4bL/64A19H1kIEd28LSCy8spfzkMjc0j npGeK8eW+/9Zv4HNfzmOb66xTmP9L4cwaGljFN3vI0fGYyOFGctfjqD+y0w2 M/7L0f3GpVBrhr8cvm3ahR2m+7RyXBSH7pjB8m+Nx+LDCE/Kp1AOp5m68Y0M fzmiB5bOCSX7cjmWW58eq2T7thwFFzYFRFB9G+QoHjAygfHTQAFFaYjAnOw5 CiR5OR5tZvgrYBy8QTGXztspUNvWbrYH288VOKqduXQ5+fNRwLLx7osD1A88 BQRl9xb+RfzjKxD0S7yRS/Nco0D14Ldid6qHVoEPQ7PGVbP8FWhnNTrjKMVT qMCNXhynnnRer0CHfF6dku2/ChQGmw75k73/CgRYPL3rxN7/1vssf9zbzP5/ R4mGDqKHP1j/K5HV690bN/Y9qITm9KITCxj/lZh++/T2ttRvUCKs4kl/T5a/ EvG1MRvS2PeyEsuNa/Vsn+Er4ebN9R3Hfr9QIvfBs2OVxDetEs7ql2W32fxX 4qmHYZdf9N4VKjHrZdoTW7b/K2HkOLbgEMO/1d+ZDxvsiI+1SvwzdCHnKXv/ lFhjdr6uieIxUGFM0IOfmyg/jgoXMx0W5jH8VSh2uVLTic0/FXI8HVacYv2v guV7W2Mty1+FihVuQ6Xs/VdBvTPyDJvXfBUmhy9o0FN9NCp01R9o5lE8WhV+ DGrrsIr9nqGC6eEdA3kMfxWGbPSfbcn2XxVeRcgO/WT7rwpzSl7rAtn7p0Lk rLgQa9pnGlQ4pK/imxEfDdT4kGn69ArxlaPGTJPHtsvY97AaW75POg+GvxoX r7VdXcT6Xw2r4uLB81n+anTZ6752M5v/aiyT5yoGs/zVOPKSs2Yd2WvUyDNx ++TA+K/G/vowYyM2/9WYlF2xp571vxp/eXXYOorhr0b4Hn00n81/NW50fHR+ E+t/Neb0n1LUmb3/apy+FmM7hu0/Wfi1bP6r6WTPyULpG6OTy4hv3Cw4rtyg W0R8ssvCs789hqoJL2QhpYeF8cPW+/4P7N+yEQ== "]]}, {}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesOrigin->{0, 0}, PlotRange->{{0., 500.}, {-23.3482827861482, 2.852994782850906}}, PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Scaled[0.02]}]], "Output", CellChangeTimes->{3.414053928152144*^9, 3.414054402584344*^9, 3.414054615640704*^9, 3.4140546936128225`*^9, 3.4140554532351055`*^9, 3.414065847170851*^9}, ImageCache->GraphicsData["CompressedBitmap", "\<\ eJztXGtsFNcVHnZmdr22wYQWNX2kQVXpj/5rEkhaab02BmNc7BgCtBKgmiYt oCok9aPQBhJVrQqt7YSkBOJd29jYuKIOgpbEPBoBFRIyUhBCKoJfIFkKCkQE YSC84t3OufeemXNn78zueg0YAsruzJx7nt8598ydu+NULWtY8auXlzWsfHHZ lNl1y15dsfLF+inlr9RZJH2cpo1boWna6SkanCetU/GVtP5LroQvTfs6O+hi IMk/8D3/7o0XcdqaE8d6kq3NkS1H/tPMZUrSypzrjS8EmRttG8sylDE6Du9d DzJoiw1G04pxUy0RNJcixqQEc8Xg2YGTfZ21p3d0LKk987/+dDbyFBDoCtrj XD3SkD5JTeaY6/v6328A0qm+rUu1/WmdkDDVafCcBuoCcKHpPFVCvMAzBpPT ADsHPx/+EPfY4R+OtRTbPBM85VIACCgACFKXOKutJiSrARqc82gNGq2LjBID u7Yv1zk2cIokjdPCg9tbF6Aj1Ag6CPFqBalq4ZxnJ0T9RhPIqom0sDgxLbpf TCGJVcTl8lqORMjJXjuuEB6d8mgGpyEC/IqPUBrI56nJAu6rHW+XU/IErupC 9+ZqW3Ehp1FWLazgK/KmSaVteE8ak31PoPDC5D+0Z8cqmPhe9DvxN0sg2xyW Qq/CoOlA2uedmyo5FkEkw+dm28YZydamSCLWHNHVUHEp0w0sqzqTj/EZ12RP CzY2EbPVYtN57gJOCPJoEL4ZKawYdcTYWcDN0iJYuBGkcw7uKFWH/EIbsDVF sHtoeUoLzWwAJQy3SkYOeXluZuivE5z/qOkebcFYDUqD8zw1OSUMTp6gUFyo oIUVtCJ/Gr/h2vXvNzXw3sFLHu/vySTUMVxLzd2irz1+tOtax8bytOl2pUDH PuPQ4CPICRsXX8hUmTLRETa5ivd9+H49WwDA2ERFyemU1qygCb6gqIAma+Bm 21tl9kDIMUh9Y5wGVnVzJHXc4VMWNiFTSSSP5xE2Rc53b6lpOXpo09ne2KLb 8TdLIXNDVkbASTzyrDUpKh8nHZpplsEVASeYvOA3FDpUfOmiUpFDGJUKJo0A AjwYlS9ZpxqFn6Y/Nin8RkLwIJ2Tq6g4ZgdYqHp6TIjJ48XvlT5HVo4wF1t+ R6qf6vGSd/vjvr7U+fdKXFyLJZ1PDyrAW+9nXZvn0l70RfvbM+GYaCULU7Fo xT40bFegnFE637FxTMYmkWil91F+bimKwrSCB4OLXZuq3CXg5neVzJzBnvde ADmKE0xVhoHlMj50qPiyOV5vf3vWxzu3vQTBu+nQ/NB/rC3aLvx0th3ev4E2 +Ex9wFzv3L/7NSrnpsO1FoDoNR3oiLW16gRaAaXhOfhMbdl4hlkRTMLCoTYy xZfixd2aDGIIBaj2goya8eNPE/ITFB4MDeFHGyoeDJvC5UcHnZctfSzFfP5M zoSXLoklHaKcUR55xFQP0udrOh52hpFEzzOSDrmlqfOioUg81Fl2nZ+2H+G6 /oO9OxtpX4KJfmXrOxXYCg7u+ecqITPZ7bW7dCGft+NvldI2SXMO2GI52uHw XBWq8uDmpflhwXzTvfmhTaSPkM7zC++t2HMtfYuRZrVA+34MvVc828voAu1H FFJT3l5KJmGJKdaEhUMd78zGbouwcjNyB5a2ZmZQ9UVu9eAOrvhuWU0PXOs7 8O91oPKTnvfmSRXKJMq9KyCgLqDUnaw5VEch6pB3yeRrqjNVRRhVyKZUpucq 5dwu40YaK1IfOffeAs+Qlz2OVpAyIuZ2XJKdkDsuys71VyurJ7XKaihfkPpP YWaE5z1ZKZyUlTUfAxhZCQRYmKniKUkU4iHRP5CVHsU+l0YB1xW1R3diBQP+ Y+6IvSJJzwjVh10yvnpYSmgGspK21ukWKQX2eQR269sVU/rYmNaUbMwj2XAw I8m8t3vnWY/rNtRjyauHOWqnW3uPO60wdVxuUWMlqq9mLkc0znoF7zBiI9n5 l65niN9kVK03q3bs3e0z1J3SokfLM65HvuckM+yoXNZ0y7LrUCZZCdhJYWem I3P3Zu1DV99fzVn9KOqHJyqnQ7uX5tn1Le3+LcU9+vOo+Dc2VY3ZJfgDV/6P on4U9T2LKhmNtG5paT7iN35jVlmb58LqlVV7avf2n5HG2btEeiYdQOwwZtVp nJ0bw0s6pf2OxE4hHDLqY4UK7WYamfG29nScD3yNjd3xBznqkc/cXKJOtkZ+ u+VIc5/f+I22skbvB7FV79ae6T+mPcc8AYvOIs9ZOUyHM8NrdJotqytGn2XW 1I+wLgueOnTquL1FvWZ7zwkbdpUZLXXtpKLB5xmuWuLz4vfSMS0DPn/7I3Z/ ejamx8HZOH7mxhUKZuj12evtcs0R12mjIT4ySJ+yi1Zsd7h8uHeAJjdH6gDQ g6+u2pMroFk55pJ92p5QfDSrff+nPI1kpWYUgOWdUXc6qq32vqCa49TNKSnP ZCXsV6cVxwYGT/6ytu/00iU72Aoyxzp9OrtgvNTk0gByR9Zj9zTbRpBHb+XS XQuuT/R8FIAuAa9Vw508SVesxr7+hlbWQ8SdfqodlOfN1nN0NG7Uz8mg+t2o pbDP/Wxhb3Jl5C8Yx7nehX/SFWHbf9jR8xGGDa8apvCx759Qc9JP6lDOdzaX NILkcKy4ThzrNQE1XxQJXTyfBrWZiEXqQl5mn6Vmx4NZFAXeT7dVv85UCFW3 4qUNJ/tqm661l6+G61N9S/+Kf8sBKHDXonWwDlt7vOvAtY7yNeCyprLPn7FC tCqcNSAvJo3DaoAd6hd/6DJplPC5apnTVAXIvqW3ddhP8ejynXhJPRyvbK34 vRA/aK0mgRGG4A1doJ3ZsXg9nLExsW9XMThwGhC52FW17qaVCBHowK7lG+2c C/Ev49F69FjLV/AFOQ2isGmFCr6wgibqb3D7gj/atG9z2oVunkRPu4c8k8Of JVlbWJzEfqAuUgsHWIG3HN20G/A4vWPJerYaDzg4ne1d9GeQgYoQMGOxiKZm ZwZVUmh5iMHbVg1aWWk4313zBogCG4hwjQVQ7WgcjV5vn7X6i/ay3+FsAG28 WzKZyXA9ZKHuPibEbAOZS52Vr4EM6sXHDSxzt7zICANJpVvBq6NfhPaddLLI n45vmMTi5VPCTkw+wkSHVWLcy/+O1DNMZrbRpHNLFXU2crmi7WfHHXO2MY22 3Ghj7pXnhy0mtGO3Op/X03EpQBdMBlmxilZnq8rz1gEN53Jn5R/aDm/4B66s gM+6N7Mlwu1Yab398qrTpzTRx21R8feRKh7xYqtBTQGdexWWPECRJLzvTeiq CELOvQBNwrgp30oRBEOhIeCJgcyXnyFN2JBAUNnNVJ+KVjDG9Ullm/oypyjZ oTmzO+zH1CkK2MIZmnt8lPlUtKn3gO8bWUPJnhssCG/NLI3bv278WKGmREGL 5kDzdwueL8CdvjfWHZC3fqXXb+11GV2IQ3UEQMh6DMH9C/aB70oqno/iuOiG xxZEZLgkGgM0HHXBgZXLd6E3tG0yzWVEMw81BHXY17muD/g+mT+vh/90xEHX uE4jURKJwRnWs/AcH18dz0uo5yGvPLKIcX1L9hyOvPyb3TpX/WVpNG6z4h/C i2gYh3CBRouPFvsaG/qRD+h5sriTpid8vcX0MkDEjQEjwM9waXHM5hETmZq3 tNpIIV3TFXkPKGRNDv7gogXb4ezUL5aysjG4OFza4iEFQlzcRDGkA2JiTKc1 IhDNQ71Id5tX5Z6jY1IzNtKmHBlTn6/IyA8UtHaaIFM1jwRwhe4I4fjB6sa9 kJwrlRVb4dqqqRiqDioAE8BCsd+cWdZG8mcmopEYNnIETfsaH5NnAQdOm+qU sQ3wk5x2oaa621Ufah2rOf/VivIOm3+6An1DiRJtFzjZaPlCmKJlgDas9s+q 53bhxDMUIPHiDQMgl6oqO+2tQotz7bau49csZ4EOWIuKCgLv+fk13cD3+dzK Tt64imOGMzw0p7wDMAfRwYUv9NyZEY0j9JafCYtbqt+wEBWpgmmIMQFvUM4k DGN2b5eVxrGdIZpYMOzqX4rZqUZYuovjX62Qdtby7qajoMbV8nSpseiKiS+m GLiOuOES8NOa6m080mLW+LnKPCTbyY3yyFGcB1Dg9u7szxf1Xi+f1b6vsf5D 9IjtmpWVxDUO4vdBDU0PzTiNBM6pPrADqXTzA/3i81VdwNf2tw2HkX7LSgzm EOXcPgKKvACfpG653QMxNKPiA1RQtWpcdQQXPv71SzuhcEWCQlhfbn5Wb3mM Z5cXTy5H9AUm22jqHekR5hSUzGjrptP6fseIR2xFMMPuRsz3+yhNcZ8H8zxF L1PdVEUvo7MdF0qKPmgCDXsmfKzSimmP8d6a0uSssrD1fsvTH4N2JuD9burK E9hNRwOQE473JhT4rZkz2mif1KoVFr+noBV7g2gq1jzCaUqDm6bAMREtjtk4 TlIH8lgqlPicoPHXP3S4VV/+aWUnU0X/j1e6s6+gS8v1oLwABpr7XqT9UBE9 C3RiWggUeQu67y54zgR9UGXfPu/RFMkISw9I0bRidA3F4M1MzJ2kzMQUDy0l aWWkvKXIMJEir0z5aE8Hazbj2rj/A0iiwME=\ \>"]] }, Open ]] }, Closed]] }, Open ]] }, WindowSize->{1016, 647}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, FrontEndVersion->"6.0 for Microsoft Windows (32-bit) (April 20, 2007)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[590, 23, 896, 22, 212, "Title"], Cell[1489, 47, 519, 11, 120, "Subsubtitle"], Cell[2011, 60, 559, 10, 44, "Subsubtitle", CellGroupingRules->{GroupTogetherGrouping, 10000.}], Cell[CellGroupData[{ Cell[2595, 74, 216, 4, 44, "Subsubtitle"], Cell[CellGroupData[{ Cell[2836, 82, 614, 19, 31, "Input", CellGroupingRules->{GroupTogetherGrouping, 10000.}], Cell[3453, 103, 325, 7, 31, "Input", CellGroupingRules->{GroupTogetherGrouping, 10000.}] }, Open ]], Cell[CellGroupData[{ Cell[3815, 115, 846, 23, 31, "Input"], Cell[4664, 140, 892, 25, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[5593, 170, 125, 2, 31, "Input"], Cell[5721, 174, 1230, 35, 72, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[6988, 214, 601, 17, 31, "Input"], Cell[7592, 233, 983, 29, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[8612, 267, 126, 2, 31, "Input"], Cell[8741, 271, 1358, 39, 72, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[10136, 315, 150, 3, 31, "Input"], Cell[10289, 320, 4309, 132, 72, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[14635, 457, 367, 7, 31, "Input"], Cell[15005, 466, 3654, 119, 126, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[18696, 590, 854, 25, 31, "Input"], Cell[19553, 617, 1639, 49, 50, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[21229, 671, 585, 17, 31, "Input"], Cell[21817, 690, 5012, 161, 145, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[26866, 856, 1457, 41, 92, "Input"], Cell[28326, 899, 1520, 44, 50, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[29883, 948, 512, 15, 31, "Input"], Cell[30398, 965, 754, 18, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[31189, 988, 1011, 21, 52, "Input"], Cell[32203, 1011, 5227, 168, 240, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[37467, 1184, 226, 6, 31, "Input"], Cell[37696, 1192, 526, 10, 30, "Output"] }, Open ]], Cell[38237, 1205, 186, 4, 31, "Input"], Cell[CellGroupData[{ Cell[38448, 1213, 962, 25, 72, "Input"], Cell[39413, 1240, 787, 20, 77, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[40237, 1265, 207, 5, 31, "Input"], Cell[40447, 1272, 482, 7, 30, "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[40978, 1285, 251, 3, 44, "Subsubtitle"], Cell[CellGroupData[{ Cell[41254, 1292, 188, 4, 92, "Input"], Cell[41445, 1298, 731, 14, 56, "Message"] }, Open ]], Cell[CellGroupData[{ Cell[42213, 1317, 1457, 41, 92, "Input"], Cell[43673, 1360, 1772, 48, 50, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[45482, 1413, 512, 15, 31, "Input"], Cell[45997, 1430, 1014, 22, 30, "Output"] }, Open ]], Cell[47026, 1455, 86, 3, 31, "Input"], Cell[47115, 1460, 720, 11, 31, "Input"], Cell[47838, 1473, 5703, 139, 352, "Input"], Cell[53544, 1614, 396, 11, 31, "Input"], Cell[CellGroupData[{ Cell[53965, 1629, 178, 3, 31, "Input"], Cell[54146, 1634, 370, 7, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[54553, 1646, 99, 2, 31, "Input"], Cell[54655, 1650, 25443, 427, 230, 20550, 345, "CachedBoxData", "BoxData", \ "Output"] }, Open ]] }, Closed]] }, Open ]] } ] *) (* End of internal cache information *)