(* 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[ 65836, 1488] NotebookOptionsPosition[ 64015, 1423] NotebookOutlinePosition[ 64361, 1438] CellTagsIndexPosition[ 64318, 1435] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"Raq", "=", RowBox[{"6.3782", " ", RowBox[{"10", "^", "6"}]}]}], ";", RowBox[{"Rpol", "=", RowBox[{"6.3568", " ", RowBox[{"10", "^", "6"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Raq", "-", "Rpol"}]}], "Input", CellChangeTimes->{{3.5589800338894997`*^9, 3.5589801006420164`*^9}}], Cell[BoxData["21400.`"], "Output", CellChangeTimes->{3.55898010259202*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"Clear", "[", "\"\\"", "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"G", "=", RowBox[{"6.67", " ", RowBox[{"10", "^", RowBox[{"-", "11"}]}]}]}], ";", RowBox[{"w", "=", RowBox[{"2", RowBox[{"Pi", "/", RowBox[{"(", RowBox[{"3600", " ", "24"}], ")"}]}]}]}], ";", RowBox[{"M", "=", RowBox[{"5.976", " ", RowBox[{"10", "^", "24"}]}]}], ";", RowBox[{"r", "=", RowBox[{"6.371", " ", RowBox[{"10", "^", "6"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"NSolve", "[", RowBox[{ RowBox[{ RowBox[{"G", " ", "M", RowBox[{"(", RowBox[{ RowBox[{"1", "/", RowBox[{"r", "^", "2"}]}], "-", RowBox[{"1", "/", RowBox[{ RowBox[{"(", RowBox[{"r", "+", "h"}], ")"}], "^", "2"}]}]}], ")"}]}], "\[Equal]", RowBox[{ RowBox[{"(", RowBox[{"r", "+", "h"}], ")"}], RowBox[{"w", "^", "2"}]}]}], ",", "h"}], "]"}], "//", "Simplify"}]}], "Input", CellChangeTimes->{{3.5591146106276064`*^9, 3.559114758859065*^9}, { 3.559114992875073*^9, 3.5591150019386888`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"h", "\[Rule]", "1.8505090415576959`*^9"}], "}"}], ",", RowBox[{"{", RowBox[{"h", "\[Rule]", RowBox[{"-", "1.2731117219190326`*^7"}]}], "}"}], ",", RowBox[{"{", RowBox[{"h", "\[Rule]", "10976.53046269707`"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.5591147121993837`*^9, 3.5591150775832205`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"Clear", "[", "\"\\"", "]"}], "\[IndentingNewLine]", RowBox[{"DSolve", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"x", "^", "2"}], "+", RowBox[{ RowBox[{"y", "[", "x", "]"}], "^", "2"}]}], ")"}], "^", RowBox[{"(", RowBox[{"3", "/", "2"}], ")"}]}], "\[Equal]", " ", FractionBox[ RowBox[{"A", RowBox[{"(", RowBox[{"x", "+", RowBox[{ RowBox[{"y", "[", "x", "]"}], " ", RowBox[{ SuperscriptBox["y", "\[Prime]", MultilineFunction->None], "[", "x", "]"}]}]}], ")"}]}], "x"]}], ",", RowBox[{"y", "[", "x", "]"}], ",", "x"}], "]"}]}], "Input", CellChangeTimes->{{3.559059839320614*^9, 3.55905991749235*^9}, { 3.559059972716446*^9, 3.5590599790032573`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"y", "[", "x", "]"}], "\[Rule]", RowBox[{"-", FractionBox[ RowBox[{"\[ImaginaryI]", " ", SqrtBox[ RowBox[{ RowBox[{ RowBox[{"-", "4"}], " ", SuperscriptBox["A", "2"]}], "+", SuperscriptBox["x", "6"], "+", RowBox[{"4", " ", SuperscriptBox["x", "4"], " ", RowBox[{"C", "[", "1", "]"}]}], "+", RowBox[{"4", " ", SuperscriptBox["x", "2"], " ", SuperscriptBox[ RowBox[{"C", "[", "1", "]"}], "2"]}]}]]}], RowBox[{ SuperscriptBox["x", "2"], "+", RowBox[{"2", " ", RowBox[{"C", "[", "1", "]"}]}]}]]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"y", "[", "x", "]"}], "\[Rule]", FractionBox[ RowBox[{"\[ImaginaryI]", " ", SqrtBox[ RowBox[{ RowBox[{ RowBox[{"-", "4"}], " ", SuperscriptBox["A", "2"]}], "+", SuperscriptBox["x", "6"], "+", RowBox[{"4", " ", SuperscriptBox["x", "4"], " ", RowBox[{"C", "[", "1", "]"}]}], "+", RowBox[{"4", " ", SuperscriptBox["x", "2"], " ", SuperscriptBox[ RowBox[{"C", "[", "1", "]"}], "2"]}]}]]}], RowBox[{ SuperscriptBox["x", "2"], "+", RowBox[{"2", " ", RowBox[{"C", "[", "1", "]"}]}]}]]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{{3.5590599654936333`*^9, 3.5590599837768655`*^9}, 3.5590620331936398`*^9, 3.5590623276441526`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"Clear", "[", "\"\\"", "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"y", "[", "x_", "]"}], ":=", FractionBox[ SqrtBox[ RowBox[{ RowBox[{"4", " ", SuperscriptBox["A", "2"]}], "-", SuperscriptBox["x", "6"], "-", RowBox[{"4", " ", SuperscriptBox["x", "4"], " ", "c"}], "-", RowBox[{"4", " ", SuperscriptBox["x", "2"], " ", SuperscriptBox["c", "2"]}]}]], RowBox[{ SuperscriptBox["x", "2"], "+", RowBox[{"2", " ", "c"}]}]]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"x", "^", "2"}], "+", RowBox[{ RowBox[{"y", "[", "x", "]"}], "^", "2"}]}], ")"}], "^", RowBox[{"(", RowBox[{"3", "/", "2"}], ")"}]}], "//", "Simplify"}], "\[IndentingNewLine]", RowBox[{ FractionBox[ RowBox[{"A", RowBox[{"(", RowBox[{"x", "+", RowBox[{ RowBox[{"y", "[", "x", "]"}], " ", RowBox[{ SuperscriptBox["y", "\[Prime]", MultilineFunction->None], "[", "x", "]"}]}]}], ")"}]}], "x"], "//", "Simplify"}]}], "Input", CellChangeTimes->{{3.559060041793367*^9, 3.559060090184651*^9}, 3.5590601935504313`*^9, {3.5590605172821956`*^9, 3.559060542133039*^9}, { 3.559061997999978*^9, 3.559062011618802*^9}, {3.5590626668511443`*^9, 3.5590626676779456`*^9}}], Cell[BoxData[ RowBox[{"8", " ", SuperscriptBox[ RowBox[{"(", FractionBox[ SuperscriptBox["A", "2"], SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"2", " ", "c"}], "+", SuperscriptBox["x", "2"]}], ")"}], "2"]], ")"}], RowBox[{"3", "/", "2"}]]}]], "Output", CellChangeTimes->{ 3.559060094646259*^9, 3.559060195016834*^9, 3.559060524536208*^9, { 3.55906202224242*^9, 3.559062036438445*^9}, 3.5590623294069557`*^9, 3.559062668988348*^9, 3.5590637880718994`*^9, 3.5591084427356524`*^9}], Cell[BoxData[ RowBox[{"-", FractionBox[ RowBox[{"8", " ", SuperscriptBox["A", "3"]}], SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"2", " ", "c"}], "+", SuperscriptBox["x", "2"]}], ")"}], "3"]]}]], "Output", CellChangeTimes->{ 3.559060094646259*^9, 3.559060195016834*^9, 3.559060524536208*^9, { 3.55906202224242*^9, 3.559062036438445*^9}, 3.5590623294069557`*^9, 3.559062668988348*^9, 3.5590637880718994`*^9, 3.559108442751252*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"Clear", "[", "c", "]"}], "\[IndentingNewLine]", RowBox[{"Solve", "[", RowBox[{ RowBox[{ RowBox[{"y", "[", "a", "]"}], "\[Equal]", "0"}], ",", "c"}], "]"}]}], "Input", CellChangeTimes->{{3.559060508936181*^9, 3.559060510074983*^9}, { 3.5590605769679*^9, 3.559060592583527*^9}, {3.559060749582201*^9, 3.559060759535018*^9}, {3.5590609227113028`*^9, 3.5590609443017406`*^9}, { 3.5590612419034595`*^9, 3.5590612501402736`*^9}, {3.559062313853729*^9, 3.5590623190485377`*^9}, {3.5590626515319176`*^9, 3.559062654636323*^9}, { 3.5590637741722746`*^9, 3.559063780443486*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"c", "\[Rule]", FractionBox[ RowBox[{ RowBox[{"-", SuperscriptBox["a", "3"]}], "-", RowBox[{"2", " ", "A"}]}], RowBox[{"2", " ", "a"}]]}], "}"}], ",", RowBox[{"{", RowBox[{"c", "\[Rule]", FractionBox[ RowBox[{ RowBox[{"-", SuperscriptBox["a", "3"]}], "+", RowBox[{"2", " ", "A"}]}], RowBox[{"2", " ", "a"}]]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{ 3.559060593862729*^9, 3.5590607620622225`*^9, {3.559060928467713*^9, 3.559060945674543*^9}, 3.5590612518562765`*^9, 3.559062047342864*^9, { 3.5590623205305405`*^9, 3.55906233185616*^9}, 3.5590626727635546`*^9, { 3.5590637817226877`*^9, 3.5590637903651032`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"c", "=", FractionBox[ RowBox[{ RowBox[{"-", SuperscriptBox["a", "3"]}], "-", RowBox[{"2", " ", "A"}]}], RowBox[{"2", " ", "a"}]]}], ";"}], "\[IndentingNewLine]", RowBox[{"FullSimplify", "[", RowBox[{ FractionBox[ SqrtBox[ RowBox[{ RowBox[{"4", " ", SuperscriptBox["A", "2"]}], "-", SuperscriptBox["x", "6"], "-", RowBox[{"4", " ", SuperscriptBox["x", "4"], " ", "c"}], "-", RowBox[{"4", " ", SuperscriptBox["x", "2"], " ", SuperscriptBox["c", "2"]}]}]], RowBox[{ SuperscriptBox["x", "2"], "+", RowBox[{"2", " ", "c"}]}]], ",", RowBox[{"Assumptions", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"a", ">", "0"}], ",", RowBox[{"A", ">", "0"}], ",", RowBox[{"Element", "[", RowBox[{"x", ",", "Reals"}], "]"}]}], "}"}]}]}], "]"}]}], "Input", CellChangeTimes->{{3.559060793168677*^9, 3.559060831014343*^9}, 3.559060972303789*^9, 3.5590612617154937`*^9, {3.5590613003099613`*^9, 3.5590613033363667`*^9}, 3.5590620626776905`*^9, 3.5590623550690002`*^9, 3.559062427687127*^9, {3.559062478137615*^9, 3.5590624991508517`*^9}, { 3.559062679611966*^9, 3.5590626849627757`*^9}, 3.559063807649933*^9, 3.5590640903848267`*^9, {3.559064166403759*^9, 3.559064213469041*^9}, { 3.559064248085501*^9, 3.559064257164717*^9}, 3.5590646693174353`*^9, 3.5590647656632037`*^9}], Cell[BoxData[ RowBox[{"-", FractionBox[ SqrtBox[ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"a", "-", "x"}], ")"}]}], " ", RowBox[{"(", RowBox[{"a", "+", "x"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox["a", "4"], " ", SuperscriptBox["x", "2"]}], "-", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "A"}], "+", RowBox[{"a", " ", SuperscriptBox["x", "2"]}]}], ")"}], "2"]}], ")"}]}]], RowBox[{ SuperscriptBox["a", "3"], "+", RowBox[{"2", " ", "A"}], "-", RowBox[{"a", " ", SuperscriptBox["x", "2"]}]}]]}]], "Output", CellChangeTimes->{{3.559060813651512*^9, 3.559060832979946*^9}, 3.5590609745813932`*^9, 3.559061317610391*^9, 3.559062065657296*^9, 3.5590623571438046`*^9, 3.559062429964731*^9, {3.5590624806648197`*^9, 3.559062501116455*^9}, 3.559062695118394*^9, 3.5590638106451387`*^9, 3.5590640039918756`*^9, 3.5590640930680313`*^9, 3.559064267008334*^9, 3.559064671361039*^9, 3.559064766958006*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"y", "[", "x_", "]"}], ":=", FractionBox[ SqrtBox[ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"a", "-", "x"}], ")"}]}], " ", RowBox[{"(", RowBox[{"a", "+", "x"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox["a", "4"], " ", SuperscriptBox["x", "2"]}], "-", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "A"}], "+", RowBox[{"a", " ", SuperscriptBox["x", "2"]}]}], ")"}], "2"]}], ")"}]}]], RowBox[{ SuperscriptBox["a", "3"], "+", RowBox[{"2", " ", "A"}], "-", RowBox[{"a", " ", SuperscriptBox["x", "2"]}]}]]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "[", "a", "]"}], "//", "Simplify"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "[", "0", "]"}], "//", "Simplify"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"y", "'"}], "[", "0", "]"}], "//", "Simplify"}]}], "Input", CellChangeTimes->{{3.5590613558148575`*^9, 3.5590613799168997`*^9}, 3.559062079057719*^9, {3.559062246134011*^9, 3.559062250268018*^9}, 3.5590625131596766`*^9, 3.5590627038076086`*^9, {3.5590640309019227`*^9, 3.5590640403087387`*^9}, 3.5590641000412436`*^9, 3.559064283825164*^9, { 3.559064608290129*^9, 3.559064618851348*^9}, 3.559064678630652*^9, 3.5590647736348176`*^9, 3.559115798663278*^9}], Cell[BoxData["0"], "Output", CellChangeTimes->{{3.5590613676084785`*^9, 3.559061381102502*^9}, 3.5590620819593244`*^9, 3.55906225148482*^9, 3.559062518011285*^9, 3.559062709891619*^9, 3.559064042102742*^9, 3.5590641018508463`*^9, 3.5590642859935675`*^9, 3.559064621347352*^9, 3.559064680471455*^9, 3.55906477510122*^9, 3.55911580009848*^9}], Cell[BoxData[ FractionBox[ RowBox[{"2", " ", SqrtBox[ RowBox[{ SuperscriptBox["a", "2"], " ", SuperscriptBox["A", "2"]}]]}], RowBox[{ SuperscriptBox["a", "3"], "+", RowBox[{"2", " ", "A"}]}]]], "Output", CellChangeTimes->{{3.5590613676084785`*^9, 3.559061381102502*^9}, 3.5590620819593244`*^9, 3.55906225148482*^9, 3.559062518011285*^9, 3.559062709891619*^9, 3.559064042102742*^9, 3.5590641018508463`*^9, 3.5590642859935675`*^9, 3.559064621347352*^9, 3.559064680471455*^9, 3.55906477510122*^9, 3.55911580009848*^9}], Cell[BoxData["0"], "Output", CellChangeTimes->{{3.5590613676084785`*^9, 3.559061381102502*^9}, 3.5590620819593244`*^9, 3.55906225148482*^9, 3.559062518011285*^9, 3.559062709891619*^9, 3.559064042102742*^9, 3.5590641018508463`*^9, 3.5590642859935675`*^9, 3.559064621347352*^9, 3.559064680471455*^9, 3.55906477510122*^9, 3.55911580009848*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"Clear", "[", "\"\\"", "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"y", "[", "x_", "]"}], ":=", FractionBox[ SqrtBox[ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"a", "-", "x"}], ")"}]}], " ", RowBox[{"(", RowBox[{"a", "+", "x"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox["a", "4"], " ", SuperscriptBox["x", "2"]}], "-", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "A"}], "+", RowBox[{"a", " ", SuperscriptBox["x", "2"]}]}], ")"}], "2"]}], ")"}]}]], RowBox[{ SuperscriptBox["a", "3"], "+", RowBox[{"2", " ", "A"}], "-", RowBox[{"a", " ", SuperscriptBox["x", "2"]}]}]]}], ";"}], "\[IndentingNewLine]", RowBox[{"Integrate", "[", RowBox[{ RowBox[{ RowBox[{"y", "[", "x", "]"}], "^", "2"}], ",", RowBox[{"{", RowBox[{"x", ",", RowBox[{"-", "a"}], ",", "a"}], "}"}], ",", RowBox[{"Assumptions", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"a", ">", "0"}], ",", RowBox[{"A", ">", "0"}], ",", RowBox[{"Element", "[", RowBox[{"x", ",", "Reals"}], "]"}]}], "}"}]}]}], "]"}]}], "Input", CellChangeTimes->{{3.559115686592683*^9, 3.559115746637187*^9}, 3.559115786058456*^9}], Cell[BoxData[ FractionBox[ RowBox[{"2", " ", SuperscriptBox["a", "2"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", SuperscriptBox["a", "4"]}], "+", RowBox[{"a", " ", "A"}], "+", FractionBox[ RowBox[{"6", " ", SuperscriptBox["A", "2"], " ", RowBox[{"ArcTanh", "[", SqrtBox[ FractionBox[ SuperscriptBox["a", "3"], RowBox[{ SuperscriptBox["a", "3"], "+", RowBox[{"2", " ", "A"}]}]]], "]"}]}], SqrtBox[ RowBox[{ SuperscriptBox["a", "4"], "+", RowBox[{"2", " ", "a", " ", "A"}]}]]]}], ")"}]}], RowBox[{"3", " ", RowBox[{"(", RowBox[{ SuperscriptBox["a", "3"], "+", RowBox[{"2", " ", "A"}]}], ")"}]}]]], "Output", CellChangeTimes->{3.5591157559816036`*^9, 3.559115790598064*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"Clear", "[", "\"\\"", "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"G", "=", RowBox[{"6.67", " ", RowBox[{"10", "^", RowBox[{"-", "11"}]}]}]}], ";", RowBox[{"w", "=", RowBox[{"2", RowBox[{"Pi", "/", RowBox[{"(", RowBox[{"3600", " ", "17.2"}], ")"}]}]}]}], ";", RowBox[{"M", "=", RowBox[{"5.976", " ", RowBox[{"10", "^", "24"}]}]}], ";", RowBox[{"A", "=", RowBox[{"G", " ", RowBox[{"M", "/", RowBox[{"w", "^", "2"}]}]}]}], ";", RowBox[{"V", "=", RowBox[{"4", " ", "Pi", " ", RowBox[{ RowBox[{"r", "^", "3"}], "/", "3"}]}]}], ";", RowBox[{"r", "=", RowBox[{"6.371", " ", RowBox[{"10", "^", "6"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{"sol", "=", RowBox[{"NSolve", "[", RowBox[{ RowBox[{ RowBox[{"G", " ", "M", RowBox[{"(", RowBox[{ RowBox[{"1", "/", RowBox[{"r", "^", "2"}]}], "-", RowBox[{"1", "/", RowBox[{ RowBox[{"(", RowBox[{"r", "+", "h"}], ")"}], "^", "2"}]}]}], ")"}]}], "\[Equal]", RowBox[{ RowBox[{"(", RowBox[{"r", "+", "h"}], ")"}], RowBox[{"w", "^", "2"}]}]}], ",", "h"}], "]"}]}], "\[IndentingNewLine]", RowBox[{"h", "=", RowBox[{"h", "/.", RowBox[{"sol", "[", RowBox[{"[", "3", "]"}], "]"}]}]}]}], "Input", CellChangeTimes->{{3.5591146106276064`*^9, 3.559114758859065*^9}, { 3.559114992875073*^9, 3.5591150019386888`*^9}, {3.5591166409711466`*^9, 3.5591166684895945`*^9}, 3.5591167652097635`*^9, {3.5591168565635223`*^9, 3.55911686063513*^9}, {3.5591169401952686`*^9, 3.5591169918625584`*^9}, { 3.559117037367838*^9, 3.55911704410705*^9}, {3.5591172704322443`*^9, 3.5591172811338625`*^9}, {3.5591292185825796`*^9, 3.55912921867618*^9}, { 3.559129275725479*^9, 3.55912932211996*^9}, {3.559129358499223*^9, 3.559129359591225*^9}, 3.5591294113365154`*^9, {3.559129443316571*^9, 3.559129460913402*^9}, {3.559129503922677*^9, 3.559129532018326*^9}, { 3.559130108361331*^9, 3.559130130248169*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"h", "\[Rule]", "9.473118856891838`*^8"}], "}"}], ",", RowBox[{"{", RowBox[{"h", "\[Rule]", RowBox[{"-", "1.2720896289846623`*^7"}]}], "}"}], ",", RowBox[{"{", RowBox[{"h", "\[Rule]", "21459.130308216343`"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{ 3.5591160387164965`*^9, 3.559116559117804*^9, {3.5591166446839533`*^9, 3.5591166702211976`*^9}, 3.559116766410965*^9, 3.559116862538333*^9, 3.559116941833271*^9, {3.55911698704215*^9, 3.5591169929233603`*^9}, { 3.55911703884984*^9, 3.5591170457606525`*^9}, {3.559117271867447*^9, 3.5591172852210703`*^9}, 3.5591292208445835`*^9, 3.5591293852220697`*^9, { 3.559129435844158*^9, 3.5591294633782063`*^9}, {3.5591295073234825`*^9, 3.5591295189611034`*^9}, {3.5591300933541045`*^9, 3.5591301322293725`*^9}}], Cell[BoxData["21459.130308216343`"], "Output", CellChangeTimes->{ 3.5591160387164965`*^9, 3.559116559117804*^9, {3.5591166446839533`*^9, 3.5591166702211976`*^9}, 3.559116766410965*^9, 3.559116862538333*^9, 3.559116941833271*^9, {3.55911698704215*^9, 3.5591169929233603`*^9}, { 3.55911703884984*^9, 3.5591170457606525`*^9}, {3.559117271867447*^9, 3.5591172852210703`*^9}, 3.5591292208445835`*^9, 3.5591293852220697`*^9, { 3.559129435844158*^9, 3.5591294633782063`*^9}, {3.5591295073234825`*^9, 3.5591295189611034`*^9}, {3.5591300933541045`*^9, 3.5591301322293725`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"y", "[", "x_", "]"}], ":=", FractionBox[ SqrtBox[ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"a", "-", "x"}], ")"}]}], " ", RowBox[{"(", RowBox[{"a", "+", "x"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox["a", "4"], " ", SuperscriptBox["x", "2"]}], "-", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "A"}], "+", RowBox[{"a", " ", SuperscriptBox["x", "2"]}]}], ")"}], "2"]}], ")"}]}]], RowBox[{ SuperscriptBox["a", "3"], "+", RowBox[{"2", " ", "A"}], "-", RowBox[{"a", " ", SuperscriptBox["x", "2"]}]}]]}], ";", RowBox[{"b", "=", RowBox[{"y", "[", "0", "]"}]}], ";", RowBox[{"Raq", "=", RowBox[{"6.3782", " ", RowBox[{"10", "^", "6"}]}]}], ";", RowBox[{"Rpol", "=", RowBox[{"6.3568", " ", RowBox[{"10", "^", "6"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"sol", "=", RowBox[{"FindRoot", "[", RowBox[{ RowBox[{"V", "\[Equal]", RowBox[{"Pi", FractionBox[ RowBox[{"2", " ", SuperscriptBox["a", "2"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", SuperscriptBox["a", "4"]}], "+", RowBox[{"a", " ", "A"}], "+", FractionBox[ RowBox[{"6", " ", SuperscriptBox["A", "2"], " ", RowBox[{"ArcTanh", "[", SqrtBox[ FractionBox[ SuperscriptBox["a", "3"], RowBox[{ SuperscriptBox["a", "3"], "+", RowBox[{"2", " ", "A"}]}]]], "]"}]}], SqrtBox[ RowBox[{ SuperscriptBox["a", "4"], "+", RowBox[{"2", " ", "a", " ", "A"}]}]]]}], ")"}]}], RowBox[{"3", " ", RowBox[{"(", RowBox[{ SuperscriptBox["a", "3"], "+", RowBox[{"2", " ", "A"}]}], ")"}]}]]}]}], ",", RowBox[{"{", RowBox[{"a", ",", RowBox[{"r", "+", "h"}]}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"a", "=", RowBox[{"a", "/.", "sol"}]}], "\[IndentingNewLine]", RowBox[{"a", "-", "b"}], "\[IndentingNewLine]", RowBox[{"Raq", "-", "Rpol"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Pi", FractionBox[ RowBox[{"2", " ", SuperscriptBox["a", "2"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", SuperscriptBox["a", "4"]}], "+", RowBox[{"a", " ", "A"}], "+", FractionBox[ RowBox[{"6", " ", SuperscriptBox["A", "2"], " ", RowBox[{"ArcTanh", "[", SqrtBox[ FractionBox[ SuperscriptBox["a", "3"], RowBox[{ SuperscriptBox["a", "3"], "+", RowBox[{"2", " ", "A"}]}]]], "]"}]}], SqrtBox[ RowBox[{ SuperscriptBox["a", "4"], "+", RowBox[{"2", " ", "a", " ", "A"}]}]]]}], ")"}]}], RowBox[{"3", " ", RowBox[{"(", RowBox[{ SuperscriptBox["a", "3"], "+", RowBox[{"2", " ", "A"}]}], ")"}]}]]}], "//", "N"}], "\[IndentingNewLine]", "V"}], "Input", CellChangeTimes->{{3.5591158398005495`*^9, 3.5591158414073524`*^9}, { 3.559115872794607*^9, 3.5591158763514137`*^9}, {3.5591160672957463`*^9, 3.5591161034878097`*^9}, {3.559116139648672*^9, 3.559116226400424*^9}, 3.559116273153705*^9, {3.5591163762386847`*^9, 3.5591163966123204`*^9}, { 3.559116530678954*^9, 3.5591165455769806`*^9}, 3.55911668869163*^9, 3.559116775646181*^9, {3.559116875174355*^9, 3.5591169144084234`*^9}, { 3.5591169981337695`*^9, 3.5591169999901724`*^9}, {3.5591292924331083`*^9, 3.559129295069513*^9}, {3.55912936225883*^9, 3.559129370261644*^9}, { 3.559129426499742*^9, 3.5591294451261744`*^9}, {3.559129667114561*^9, 3.559129698002615*^9}, {3.559129747797902*^9, 3.559129750871107*^9}}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindRoot", "::", "lstol"}], "MessageName"], RowBox[{ ":", " "}], "\<\"The line search decreased the step size to within \ tolerance specified by AccuracyGoal and PrecisionGoal but was unable to find \ a sufficient decrease in the merit function. You may need more than \ \[NoBreak]\\!\\(MachinePrecision\\)\[NoBreak] digits of working precision to \ meet these tolerances. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/FindRoot\\\", ButtonNote -> \ \\\"FindRoot::lstol\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{ 3.5591162325312347`*^9, 3.5591162749321084`*^9, 3.5591163984063234`*^9, 3.559116561426608*^9, 3.5591166929660373`*^9, 3.559116777814585*^9, { 3.5591168777951593`*^9, 3.5591169156720257`*^9}, 3.5591170032193785`*^9, 3.5591170569458723`*^9, {3.559129371993247*^9, 3.5591293878584747`*^9}, { 3.559129429245347*^9, 3.5591294657962103`*^9}, {3.5591295118786907`*^9, 3.5591295208955064`*^9}, {3.559130098923314*^9, 3.5591301382665825`*^9}}], Cell[BoxData["6.385288452913401`*^6"], "Output", CellChangeTimes->{ 3.559116107060216*^9, 3.5591162325312347`*^9, 3.5591162749321084`*^9, 3.5591163984063234`*^9, 3.559116561426608*^9, 3.559116692981637*^9, 3.559116777830185*^9, {3.5591168778107595`*^9, 3.5591169156720257`*^9}, 3.5591170032193785`*^9, 3.5591170569458723`*^9, {3.5591293722428474`*^9, 3.5591293878740745`*^9}, {3.559129429245347*^9, 3.5591294658118105`*^9}, { 3.5591295118786907`*^9, 3.5591295208955064`*^9}, {3.5591300989389143`*^9, 3.5591301382665825`*^9}}], Cell[BoxData["21399.10482253693`"], "Output", CellChangeTimes->{ 3.559116107060216*^9, 3.5591162325312347`*^9, 3.5591162749321084`*^9, 3.5591163984063234`*^9, 3.559116561426608*^9, 3.559116692981637*^9, 3.559116777830185*^9, {3.5591168778107595`*^9, 3.5591169156720257`*^9}, 3.5591170032193785`*^9, 3.5591170569458723`*^9, {3.5591293722428474`*^9, 3.5591293878740745`*^9}, {3.559129429245347*^9, 3.5591294658118105`*^9}, { 3.5591295118786907`*^9, 3.5591295208955064`*^9}, {3.5591300989389143`*^9, 3.5591301382665825`*^9}}], Cell[BoxData["21400.`"], "Output", CellChangeTimes->{ 3.559116107060216*^9, 3.5591162325312347`*^9, 3.5591162749321084`*^9, 3.5591163984063234`*^9, 3.559116561426608*^9, 3.559116692981637*^9, 3.559116777830185*^9, {3.5591168778107595`*^9, 3.5591169156720257`*^9}, 3.5591170032193785`*^9, 3.5591170569458723`*^9, {3.5591293722428474`*^9, 3.5591293878740745`*^9}, {3.559129429245347*^9, 3.5591294658118105`*^9}, { 3.5591295118786907`*^9, 3.5591295208955064`*^9}, {3.5591300989389143`*^9, 3.5591301382665825`*^9}}], Cell[BoxData["1.0832069168457528`*^21"], "Output", CellChangeTimes->{ 3.559116107060216*^9, 3.5591162325312347`*^9, 3.5591162749321084`*^9, 3.5591163984063234`*^9, 3.559116561426608*^9, 3.559116692981637*^9, 3.559116777830185*^9, {3.5591168778107595`*^9, 3.5591169156720257`*^9}, 3.5591170032193785`*^9, 3.5591170569458723`*^9, {3.5591293722428474`*^9, 3.5591293878740745`*^9}, {3.559129429245347*^9, 3.5591294658118105`*^9}, { 3.5591295118786907`*^9, 3.5591295208955064`*^9}, {3.5591300989389143`*^9, 3.5591301382665825`*^9}}], Cell[BoxData["1.0832069168457536`*^21"], "Output", CellChangeTimes->{ 3.559116107060216*^9, 3.5591162325312347`*^9, 3.5591162749321084`*^9, 3.5591163984063234`*^9, 3.559116561426608*^9, 3.559116692981637*^9, 3.559116777830185*^9, {3.5591168778107595`*^9, 3.5591169156720257`*^9}, 3.5591170032193785`*^9, 3.5591170569458723`*^9, {3.5591293722428474`*^9, 3.5591293878740745`*^9}, {3.559129429245347*^9, 3.5591294658118105`*^9}, { 3.5591295118786907`*^9, 3.5591295208955064`*^9}, {3.5591300989389143`*^9, 3.5591301382665825`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"Curve", "=", RowBox[{"Plot", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"y", "[", "x", "]"}], ",", RowBox[{"-", RowBox[{"y", "[", "x", "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"x", ",", RowBox[{"-", "a"}], ",", "a"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", "Blue"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Ellipse", "=", RowBox[{"ParametricPlot", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"a", " ", RowBox[{"Cos", "[", "p", "]"}]}], ",", RowBox[{"b", " ", RowBox[{"Sin", "[", "p", "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"p", ",", "0", ",", RowBox[{"2", "Pi"}]}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", "Red"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{"Curve", ",", "Ellipse", ",", RowBox[{"AspectRatio", "\[Rule]", "Automatic"}]}], "]"}]}], "Input", CellChangeTimes->{{3.5590653376850014`*^9, 3.5590653397130046`*^9}, { 3.5590653724418616`*^9, 3.5590654461987906`*^9}, {3.5590654974292793`*^9, 3.559065714410058*^9}, {3.559065765562547*^9, 3.559065797074602*^9}, { 3.5590658300998597`*^9, 3.5590658374630723`*^9}, {3.559065974883712*^9, 3.5590660875783086`*^9}, 3.5591086372367916`*^9, {3.559108668561646*^9, 3.5591087628950105`*^9}, {3.5591087960918684`*^9, 3.5591088284775248`*^9}, 3.559115560092062*^9, 3.5591163185965843`*^9, {3.559116412789549*^9, 3.5591164482172103`*^9}, {3.55911659128506*^9, 3.559116625277519*^9}, { 3.5591167174892797`*^9, 3.5591167476441326`*^9}, {3.5591292697506685`*^9, 3.559129271107871*^9}}], Cell[BoxData[ GraphicsBox[{{{}, {}, {RGBColor[0, 0, 1], LineBox[CompressedData[" 1:eJw12Hc8Vf//AHCylbJXyCZ773Ffsi5Z4Z57Lw2lPkQICUWyypYUSosKlchu cd/4SGV8rMqKioiWImT0ffd4/H5/ncfz8TqPc97nPN7nNY7M/uBdBzcwMTHt Z2Zi+nvMEQo6n/KIjvaZNz0Kd7tEWsqRFSPq6Mh+a0+4pR4X7OF9fVW+io7U ftaVTvXwQWtGmvyPcjoKJ+/mz+sVApWNpDtNZXTUdmPrc00dMcg+O6+ZcYuO 3IS5/ulSkYBfbGW19BvYJ2n+MCQFzUz8aCGfjrhf0HYpdcuAcmybbUsuHQlV rkebfpCFzJXojuxsOprp5im3nZYD+q8Pr1XP0tGZJM07uZ0KMDdb+9XnOB39 fBTrcFJlO1D8D4drhtERm56Kndj37fDko9Tv1SA6qjguzfO5QgXOvjvDln+I jqxDQvW5pNVA+jVNosuTjqKXvSkCHRrg2rxGNtGjo9i01n3ardpQZ1nVzaFF R5EMWwELbh2QeHrIc0CVjk619SUZuujAdH33vhA5Olqdj3iw7ZUOxJXfOH6b n45ox4plPrzRhcp8m1v8czRU0++gBw/1QVj4t8r4Zxp6pyoq5vxDH06ev19Z Pk1DZ4eqOJtVDMA+U+SJ/TgNLYmojX4vMIDxhE+9sd00tHH1ScnwUUPYEpzB NFtOQ73zV0w+cBvDZ1VVhT1lNCRAPxxvRDKGF9Pt5J6bNLT5y0yJyDFjSNzP er7+Mg15l6t2x44aw2/PEwqJKTQ0eoybblxmAhNm/g6SB2mocZHFiqZhBo+4 bXJdJGjoCkkxJnnFAvKevW9oFqGhoPdt6hyKlhCeGDeqL0BDVe9PtNi7WoLm +mNFCW4a0g/ISeG/aQk353QeTi9SEe+lDzuOyJMg643029N9VCSRNTGw2kIC 39urSjWpVLRVZQdb7XeADTu4YsOTqSjSOf7IeWYruD4m1K8XT0WVR7SX9vJb waioZlxtFBVFax0watW1Akqaz5s6PypKvqie1njMCsghbWcf2lJRppsY+41f VqBhkv2pcQMVnW6J/Vftyw5Y7pS7+zyaQAvbDXpFa2yAdFiSuBZBoLC6uKce yAaS2UVYjoUSqG059FVylw0IWnJ7SR8m0B/Gk1ed0zagVfmd+ziNQJKvZGan JWzhUM7Tw/JGBNrV2pLhk2ALvRRC5dQCBS22Uy63ku3g3lhqmV4IBc3V3e4S aLIHw27Yz+HiifaYqdaxnnOEDRWVEYrggTIf1an6jTrBrcBgm3IJd9Tz+pYu N5crHBQMykpZc0N5Jf99e/XHFQ7dk9evmnNFT6+Sl+e63CDGmqKSdNUFvWMX Sj9RsAsqePtuWMQ5oQCGcMtma3eQfNC9pTjQES20pyckDbuDCU/gj1eqZLRF 6o/HlWAPeEOOTC/ms0Ma6XX/Fcx6gIKV/GXLF9aI79l8fdQBTyjf7Ho9OcMK hTDZtVu+9IT6x7Peh46TUF1PqaumNAXG7J0lVP3MUOb44seJYAooF5/3siHM kLzgk8u/QigQuvbm0l5bM7T+sqaII5QC7FUHxHLlzRCl+mKJUjgFNLZGCa+N m6Je4d1ZeyMpcPJLMW83zRTRrjdP34+jgFjOMmuYowka4Vl16simgPvwzS+P NIyQH3MQG3sFBaJ1s0oDJIxQRWqZHFMlBYrSog5IcBuhuGjqiWXs76bOgzGT hsi//eDrT1UUyCxcbCUVGqJci0xrVEeB9t0OhW1chmjilwebYyMFzMa/OfRO 6KPl61n3X3RQ4IDREFtCrz5ilWeRu9dJgdTsVoYuQx9xXj/Sl95FgUFSgcGF y/qIpyd02vE/CkTcsJKl7tJH3hzkSUYfBSr3X1gebdJDsg/NTZKHKKAwaVL6 6ZIuWjys0xU7RYGvpoOPvp3VRVJDG82tpynQkHO8cyFCFx0Z4p3j/EQBB1LN D2Z3XZT9Z1g5e4YCQZfVzcU36qJJG+XBC18oUO8m3etwQgfVa34hx/6kALmR bfWulzbyvn2or2WdAvxCNzdXkbWRaWUda+gfCowEWMk0GGojZhJ7zDYmAoLE Ttm2CmijMbUhq0hmAnLCl7JHXmohwz2d8jKsBAyrzCrwmGkhdbYgJ1cubI0b 10aMNFC2c2uqvAABEUMerrtENJDdTurFB9i8yZzM7QvqyOW5/QYLQQKsR4P3 V1epo5UF3xUPIQLup1kopKqpI6f4rLATIgTETQ/fMZRWQ4bDVQl3thKgWCRc l8Ohgmyvs/mckyOA4fTiEOfH7YhPOTeFW54Ar+UYkdjW7aiKJU8gATvbdTLS //R2VP5YbGeoAgHL69WmpBVldP/9dxeyEgEdXm7oy1clZGJQUTaoQsBRofQO +1cKaCHi6GCCNgFtT1LCZx8oIKWC92gOW8L3jERmhgLaWcE1t0cHx6viA/us FZBT7GU5Q10cd43atKdaHvn21Q291cPx1EM7w7LlkMrkRN5GIwK26vjOCwXK oZXGtG9HsEMGfQob7OTQ8xm/zO6/ceXdn9fWZVE7eeJGtjGO/+uefvaILGL8 c1+a25QAcSbouOIgg9QFhhKGzAkILrEIJynKoIZqp13aFgS0OptJfGCWQdvs /P3OYAdfMQxUfiiN2iqqArQtcdxEY1O1kjQim1l3nyQRMDJjkJ/NJYUGaeNb 5qwISBqYzZuMEUMB3kvXXOwIIEkPrK8ZiyHyrKFxGvbvw40HhX+JIrcWX502 7BCmc7p2QaJI5/cFaxN7AuiqBt2le0TQuwE0LU4mQOt0HHugpRBKPR6h2+hA wOxL/yOJK4LoJ+HQ+Q37trB7f2G9IPpnp80zaUf8/PcUijo1BVF7wNPhOGz2 1y/MNaUF0BzrljPGO/F61YSO/WDiQ6uJbz6lOxGQd3x9mPspLwruTT9Tg72r ecpKLooXCVnnJQ1jZ1aEy5S+34Is3jcUKTkT4MvWRA1o4UElg+E+DdhDl84p FpryoKH90huGsV20fOc7qjchpnx9ljVsEy+ubI2bG9HHP7cULV3w/n3g/mwu kQv52w3/eIydbKt4QXaNEy2u8OYNYq8ML+13P8aJhkzR5V/YH9mvrtcc5EAV LAVtGq4EPPGe1o+0ZUP+n6lFedg6Px6xlDWyorrYiw6V2CVnMnoGDVjRj35v 93bsnCqdI6ZKLGindGryL2w/zpib6xzM6I7SuTBHNwI+53zTHIpnQgpZY417 sN+ZtEYErPxhkAIWQo5i//DNvDn4cY0RYnyS6wL2SlTY4+yRVYZa8Ldvt7BZ s6i9dr0rjP7rGuZ12EIPZZnqni4zPgoIbRrAluriEDlSvcSYcmjM+ICt+OGz unzZImOTv2HEHLYRT71XTu4Co3yo8MTGXXi/yBaGklPnGWggNV8Em2x4OoUp 7idDg2OnhBw23cexPihgjsEy50A1wt4fodWl4POdwcuSymuF/UpY9cXD8q+M hvWrmo7YU6SFRaWCz4zQyNRad+ylw00KFxNnGEu/dxZ5YXNdSHFnDZlmAMfU wn5s8Sb306FeHxmZgrRKf2y1T5IV47YTjD3ZxT3B2OYC0yPOOu8ZEW/baMew nc2ruJ9KjjNGqC320diNXUshFP4RxqBLwdVY7HHZjA6dT68ZK1xkr3hsluMy ylsYfQyh1t7YJGyHiwEcr9W7GcYPDLjOYmdZ5DZmDbcxmn6H/07BLumMms3m f8qwYWTsSsPeZ3axmnWklFEmcpov/e96bTLHtIaySbw8LoZ//ZSU81x8uIY0 VrTY8vd88cjr5f1NzaTmnuj6VOyB7uO/tvK/JJWcfcP39/rt50uuaVf1kFru 8fYlY9+2FLD2HB4gpYrJrSRgJ87GTkeyDpHqr29JjcP2yZtJL1R/S7LnfnXq JPY1mbui+1TfkWIEI0aOY2cujs/MWX4gDYd9uxaKHdMp/DTBY5K0vQ9eBGIH Fu/MEvKfIoX1BXodwvaKivcpiflE4toU4r737/twadA1zpklRXM7VhHYSity b7wffyMl7Em+aYudNq3muTj1nXTm6qShOXZ8n50F6cccyT5KXEcX+2hZDG8v +zypwm+DlyS2X27+sjjfAqlPrCaF7+/7O1X9/sDWX6Tqq0bb2LBdPD/VLGgu kQ7N1cXO4P1oS2K9amGyTGrmqbAfwbZQ3XbmjPVv0sdLEVGd2GrMnjQx2irJ 3Thky33sJMp4rBQwgTb6vskPW3jURjirhwnUTqX4u2PfPnC3/I8PMyyOdvJZ YLeFRIyMx28AhU3uNnzY7GmbTIpbWUFIdJtRBf4+L/KF9gh4ssFp5cm5XGzF /Nd+iZNsoHTRfjEK2/ZWUd5Bdg7YemhsArCTm4wWlOy5wPgzr+NznC+Eba+m 57/hgnlqsk8Z9u0OFnkuf27YdfJMy1nstsGuXTMpG2GST+6wDTb7vG/lvZc8 oHnvONThfLUjt702X4kPjl5RFY7A+S+NOMHcE8wH7ITYhAN2n7iGE1cDH5iT D65JYR+4cX4iyp4ftl7m5v8X59eE+3sEvPwF4Jyragsndkv7fIjkXSHYtn4J gnC+3phe8sTzpxCEjzS3GGK7u9A5M02Fwa6L5xwT9oeBpmvrL4WBacfqpmyc /1k+pHaNfRYB+Y6WlNu4PuxYk1Yr0hCHhRX+8w9sCYhundzHESoOeht2a4Vj V6XduRBYKw6/Qz9qGmDLiuv9MTDbCgmkEs8GGwI2GNr1viRLQNJ8gGmtNa7/ IUciF3ylQPKwSeA5XL+WDHXK6aVS4HL8hKoDttafX++aZqUA3tTu3IB9LSPO MTV0GwhdYakNAQIsB8hZW0jS0Ch4UsIW10OpH977bzrJgMhAolUfrq+PPJtW KcEysPRr47Z4bMpDmTyuczIQVGe6Wws7M+7jy6B+GeA+fGo91YyAtc0h+sZe snC/Y63PCNfrEdUEzi4/OTA5OPFvBK7vlw+W3V9MUIBCy6zLJ3D/cLITucWU KoAnk5D8Zmxv/aH5DXhOlQmsn7uG+w0pto2mm4UV4TPS9EZaBBTdDHwmX6oI BieWA1c0CLjzQXPcrUMJ4gtFInaq4vX61PHdE1SBWJnUCR/cH11+3lWjbawC baMv10dk8f20p4j63SogNbDkTWBbMItebb6tAqW6MxNkGQKar0epDBqpwqej f2QUtxHwYsxsB/tuNVhbO2j3WBzXy90t4ftuaYDeXOA3go+A4+1hfqbPNWDl 3GnnRl4CBHXlvYW/aAA96PMvOWwXzqQdHXqakL7jgtyXzbj/qLLlN2rVhH9n 3s9EbCLgNcvk6zNHtSBGukHRj4MARPU7tjdbGxIPBP5mrFLgwobgys1dutBW YytCnqTAeEOa77O3urBwp2Ewc4ICKsGlonHfdEHvDel33wcKMIbfxc3x6kGC 1T5P7/cUmKnxcOt31wODY8GhB8coYPGPyY+CQT0QkLc1ow1SYPolm57CpD50 vXpX0Yb7eZMLhfWmq4bACk9+1j6gQImT4Pfdm40gVCQyrwnPC4Ls6cpx0kZw 9GFCZDueL75GnChotTaCOrrhj9flFCim0k84ZRjBgbw29KWMApskRC33ShmD etVi4/wNCrwtOt8WTzKBKjuhguBzFIivTB14nmgGJuPXWgSC8DyRtNNKYB8J Qk/brhfL4f7+rXuxcJgV9Msvv6jv9oT4E4+zHKqtITXgrtnkP55gbFTdOtlh CwGTvVmxvz0AjDcmPRAgQ3/L7fqIYx7wjPYzY0LXEXyYzX/0zLgD/22v64Fk J1D6EsG/z8Md7hvIWyiku8AKE5Nd2rVd4BtmkSPb7wrnxuKlsifc4NIX8Zcv ltygfzkh5AGLG1TF1J+VFHUH5Vy5MA0eV0iVLJyRk/eAbMuSigxRZ0gd08pp tvWElj28LUo3HYH8gOetIIUCSbPzYV9ZyTAt79pjRqdAf1nhl/fL9pCcn9Pq u4cCzn0H3v731R5aT4ncrTlEgU6mFJaCN/Zg4SwT6XEcz23BNgGDd+1BZ1aP PzefAr2roR/43exBXMHbVmCYAtbNz5wVLtrBbP6d+3z7CNDN+jbizmkL/J0s LLm+uC4XPLhm9NsGTJh3E8L+BHTnnTfg/2wDZ/23MG0NJWDe4mJacbcNyJuE u8knEiCnnRlkdtEG9gybzxuUEsC0ceGXjowN/CfZY+z9nYD1Ud7BU+rWsOim kvl2Hudham1lt6Q1SCUnvN+3TMBgdg0X32ZrCPyqn3ZwAxXm0kMb/vm6A7ia Lo0ECVLhmVe1gc/9HWC1z/fUaSMqTB2I1/RV3QG1RYutt+OoQHrPHnJa1Aok Ft9NDSdSwaUtinOSwwriHTu4+VKo4F5zR91wEcB5/rrriXNUGHTe7X17AOCj jcOISxEVgmMcn985ByD6qfDnUgsVbDwfRkuwA0RrWsk6ctCgiznyv+6nljCe oGYTt5EGHNS6RPc8S7B9I+xXu4UGT/4LPVkcYgn8p2fLt4nSYHhczPCWrCXc 7c01nlemgSNR80Is0QJGjk25XnGgwfZUtjesZuZg+ST91LcMGjyXOs3mFWkC N2Fk6kYODUoaWweYrE2A+5mqq0ceDdaknrQabDGB/p4X0g+v06A4ycrA/qYx HJ7ibI6vpsGlcFH/inYjyOVPYhUepMHD8oubzy4awHJef2D7KA0SDy1+KX5k AHsk5Qei39MgLO+ih0CMAWxXbrk1NkuDkcuj/K/+6EOjObNd2ToNVEpHp41X 9WDaLybFTJ4Opft4TkU91wGnrx1zX5XpIGm8V9wmVgeqwiToN9TpsGbxrPY/ bR2IiX2swm5Ih9GlkZMfc7WBP3f5ZTeZDhS39Muf7LUgUoysF+9MhxnrHZXa HzRh9Gp+oZ47HQpmG8LLtDShtMzwSIE3HTSG6+JfNKqDedOxzQeC6XB7mP3P 1ysqUGT9b4RQOB3q39rzOLVtB84XgmPPIunwzrJj2HZKGXr7qyvU4ulgteuZ afkWJTCis4iOJWNvqHVYElSEK2O74s6l0cG9akdX7GYF8JuZc/2VS4fY+Kb9 0a9koTMYHpYW0CEqjgsSr8mAzkK2jNdVOoS/e6rn5iYN+dHjKTzFdPDbK+SW 8VIS1v5o/mgqoQPf5ndknhJx2J90ih56jw5DhlOfC4REoJ27u1n+AR33+d1G Dwf5QT1bSvV1LR28vo7sDC/jhv/7Hwz//z/4fw0gCgE= "]]}, {RGBColor[0, 0, 1], LineBox[CompressedData[" 1:eJw12Hk0Vd37AHAyK2WeQmYyz/Nwt/mSKdxz76VBqZcIIaFIpjJLSkoTFSqR uYn74JXK8JoqU1RENClChr67tX6/v876rOesc/Y5a+/nefaW2Rey88AGJiam fcxMTH+vuULB51If0WGvedOjCPdLjKVcWTGijg4OW3siLPW4YDfv66vyVXRQ +1lXOtXDB62Z6fI/yukQQd7Fn98rBCobSXeayujQdmPrc00dMcg5M6+ZeYsO 7sJc/3SpSMAvtrJa+g3sE7QANCQFzUz8sHCRDtwvaDuVumVAOa7NriWPDkKV 6zGmH2QhayWmIyeHDjPdPOV203JA//XhteoZOpxO1ryT16kAc7O1X32P0eHn ozjHEyrbgRJwKEIznA5seir2Yt+3w5OPUr9Xg+lQcUya53OFCpx5d5rt4kE6 2ISG6XNJq4H0a5pElxcdYpZ9KAIdGuDWvEY20aNDXHrrXu1WbaizrOrm0KJD FMNOwIJbBySeHvQaUKXDyba+ZENXHZiu794bKkeH1fnIB9te6UB8+Y1jt/np QDtaLPPhjS5UXrS9xT9Hg5p+Rz30UB+EhX+rjH+mwTtVUTGXH/pw4tz9yvJp GpwZquJsVjEAhyyRJw7jNFgSURv9XmAA44mfeuO6abBx9UnJ8BFD2BKSyTRb ToPe+SsmH7iN4bOqqsLuMhoI0A8lGJGM4cV0O7nnJg02f5kpETlqDEn7WM/V X6aBT7lqd9yoMfz2Oq6QlEqD0aPcdOMyE5gwC3CUPECDxkUWK5qGGTzits1z laDBFZJibMqKBeQ/e9/QLEKD4Pdt6hyKlhCRFD+qL0CDqvfHWxzcLEFz/bGi BDcN9ANzU/lvWsLNOZ2H04tU4L30wfqwPAmy30i/PdVHBYnsiYHVFhL43V5V qkmjwlYVa7ba7wg2WHPFRaRQIcol4fA5Ziu4PibUr5dAhcrD2kt7+K1gVFQz vjaaCjFa+41ada2Aku77ps6fCikX1NMbj1oBObTtzEM7KmS5i7Hf+GUFGiY5 nxo3UOFUS9y/al+sYblT7u7zGAIWthv0itbYAumQJHEtkoDwuvinnmALKewi LEfDCGhbDnuV0mULgpbc3tKHCPjDePKqc9oWtCq/cx+jESD5SmZ2WsIODuY+ PSRvRMDO1pZM30Q76KUQKicXKLDYTrncSraHe2NpZXqhFJiru90l0OQAht1o H4erF+w2U61jPesEGyoqIxWRJ2Q9qlP1H3WGW0EhtuUSHtDz+pYuN5cbHBAM zk5dc4f8kv++vfrjBgfvyetXzbnB06vk5bkud4i1oagkX3WFd+xCGccLdkIF b98Ni3hnCGQIt2y28QDJB91bioOcYKE9IzF52ANMeIJ+vFIlwxapP55XQjzh DTkqo5jPHjQy6v4rmPUEBSv5y5YvbIDv2Xx99H4vKN/sdj0l0wpCmezbLV96 Qf3jWZ+Dx0hQ11PqpilNgTEHFwlVfzPIGl/8OBFCAeXic962hBnICz65/At/ d9jam0t77Mxg/WVNEUcYBdir9ovlyZsBpfpCiVIEBTS2RguvjZtCr/Cu7D1R FDjxpZi3m2YKtOvN0/fjKSCWu8wa7mQCIzyrzh05FPAYvvnlkYYR+DMHs7FX UCBGN7s0UMIIKtLK5JgqKVCUHr1fgtsI4mOox5exv5u6DMZOGkJA+4HXn6oo kFW42EoqNIQ8iywbqKNA+y7HwjYuQ5j45cnm1EgBs/Fvjr0T+rB8Pfv+iw4K 7DcaYkvs1QdWeRa5e50USMtpZegy9IHz+uG+jC4KDJIKDM5f1geenrBpp/8o EHnDSpa6Ux98OMiTjD4KVO47vzzapAeyD81NUoYooDBpUvrpki4sHtLpipui wFfTwUffzuiC1NBGc5tpCjTkHutciNSFw0O8c5yfKOBIqvnB7KELOX+GlXNm KBB8Wd1cfKMuTNoqD57/QoF6d+lex+M6UK/5hRz3kwLkRrbVu97a4HP7YF/L OgX4hW5uriJrg2llHWvYHwqMBFrJNBhqAzOJPXYbEwHBYiftWgW0YUxtyCqK mYDciKWckZdaYLi7U16GlYBhlVkFHjMtUGcLdnbjwta4cW3ESANyXFrT5AUI iBzydNspogH2O6gXHmDzpnAyty+og+tzhw0WggTYjIbsq65Sh5UFvxVPIQLu p1sopKmpg3NCdvhxEQLip4fvGOK8azhclXhnKwGKRcJ1uRwqYHedzfesHAEM 5xcHOT9uBz7lvFRueQK8l2NF4lq3QxVLvkAido7bZFTAqe1Q/lhsR5gCAcvr 1aakFWW4//67K1mJgA5vd/jyVQlMDCrKBlUIOCKU0eHwSgEWIo8MJmrjdf4k NWL2gQIoFbyHOWwJv9MSWZkKsKOCa263Do5XJQT12SiAc9xlOUNdHHeL3rS7 Wh78+uqG3urheNrBHeE5cqAyOZG/EeeBrTp+80JBcrDSmP7tMHbooG9hg70c PJ/xz+r+G1fe9XltXRbayRM3coxx/F+PjDOHZYHxz31pblMCxJlQxxVHGVAX GEocMicgpMQigqQoAw3Vzju1LQhodTGT+MAsA9vsA/xPY4dcMQxSfigNbRVV gdqWOG6isalaSRrIZjbdJ0gEjMwYXMzhkoJB2viWOSsCkgdm8ydjxSDQZ+ma qz0BJOmB9TVjMSDPGhqnY/8+1HhA+JcouLf46bRhhzKd1bUPFgWd3+dtTBwI oKsadJfuFoF3AzAtTiZA61Q8e5ClEKQdi9RtdCRg9mXA4aQVQfhJOHZ+w74t 7NFfWC8I/+ywfSbthL//nkJRp6YgtAc+HY7HZn/9wlxTWgDmWLecNt6Bx6sm dPQHEx+sJr35lOFMQP6x9WHup7wQ0ptxugZ7Z/OUlVw0LwjZ5CcPY2dVRMiU vt8CFu8bipRcCPBja6IGtvBAyWCEbwP20KWzioWmPDC0T3rDMLarlt98R/Um YLqoz7KGbeLNlaNxcyN8/HNL0dIVz98HHs/mkrggwH74x2PsFDvF87JrnLC4 wps/iL0yvLTP4ygnDJnC5V/YH9mvrtcc4IAKloI2DTcCnvhM60fZsUHAZ2pR PrbOj0csZY2sUBd3wbESu+R0Zs+gASv86PfxaMfOrdI5bKrEAjuk01J+Yftz xt5c52CGO0pnw53cCfic+01zKIEJFLLHGndjvzNpjQxc+cMgBS6EHsH+4Zd1 c/DjGiPU+ATXeeyV6PDHOSOrDLWQb99uYbNmU3vte1cY/dc1zOuwhR7KMtU9 XWZ8FBDaNIAt1cUhcrh6iTHl2Jj5AVvxw2d1+bJFxqYAw8g5bCOeeu/cvAVG +VDh8Y078XyRLQwjp80zYCDtogg22fBUKlP8T4YGxw4JOWy6r1N9cOAcg2XO kWqEvS9Sq0vB9zuDlyWN1wr7lbDqi4flXxkN61c1nbCnSAuLSgWfGWFRabUe 2EuHmhQuJM0wln7vKPLG5jqf6sEaOs1AHFML+7DFmzxOhXl/ZGQJ0ioDsNU+ SVaM200wducU94RgmwtMj7jovGdEvm2jHcV2Ma/ifio5zhihtjjEYDd2LYVS +EcYg64FV+Owx2UzO3Q+vWascJG9E7BZjskob2H0MYRae+OSsR0vBHK8Vu9m GD8w4DqDnW2R15g93MZo+h3xOxW7pDN6Nof/KcOWkbkzHXuv2YVq1pFSRpnI Kb6Mv+O1zRrTGsoh8fK4Gv71U1Luc/HhGtJY0WLL3/vFo66X9zc1k5p7YurT sAe6j/3ayv+SVHLmDd/f57efK7mmXdVDarnH25eCfdtSwMZreICUJia3koid NBs3HcU6RKq/viUtHts3fyajUP0tyYH71ckT2Ndk7oruVX1HihWMHDmGnbU4 PjNn+YE0HP7tWhh2bKfw00TPSdL2PvQiCDuoeEe2UMAUKbwvyPsgtnd0gm9J 7CcS16ZQjz1//4drg65x7iwphtupisBWWpF74/P4Gylxd8pNO+z0aTWvxanv pNNXJw3NsRP67C1IP+ZIDtHiOrrYR8pieXvZ50kV/hu8JbH98y4ui/MtkPrE alL5/v6/k9Xv92/9Raq+arSNDdvV61PNguYS6eBcXdwMno92JNarFibLpGae CocRbAvVbadP2/wmfbwUGd2JrcbsRROjrZI8jEO33MdOpozHSSEmpA3fN/lj C4/aCmf3MCG1k6kBHti3998t/+PLjBZHO/kssNtCI0fGEzYghU0etnzY7Omb TIpbWZGQ6DajCrw+L/CF9Qh4saFTypNzediKF1/7J02yIaULDovR2Ha3ivIP sHOgrQfHJhB2SpPRgpIDFzL+zOv0HOcLYburGRffcKF5aopvGfbtDhZ5rgBu tPPE6ZYz2G2DXTtnUjeiST65Q7bY7PN+lfde8iDNe8dQHc5X1nnttReV+NCR K6rCkTj/pRPHmXtC+BA7ITbhiN0nruHM1cCHzMkH1qSw9984NxHtwI+2Xubm /xfn18T7uwW8AwTQWTfVFk7slvb5UMm7Qmjb+iUUjPP1xoySJ14/hVDESHOL IbaHK50zy1QY2XfxnGXC/jDQdG39pTBisl7dlIPzP8uHtK6xzyJIvqMl9Tau D9Zr0mpFGuJoYYX/3AM7AmJaJ/dyhIkjvQ27tCKwq9LvnA+qFUe/wz5qGmDL iuv9MTDbihJJJV4NtgRsMLTvfUmWQMnzgaa1Nrj+hx6OWvCTQpKHTILO4vq1 ZKhTTi+VQq7Hjqs6Ymv9+fWuaVYKoTe1OzZgX8uMd0oL24aErrDUhiICLAfI 2VtI0qhR8ISEHa6HUj989t10lkEiA0lWfbi+PvJqWqWEyKClXxu3JWBTHsrk c52VQcF1pru0sLPiP74M7pdB3IdOrqeZEbC2OVTf2FsW3e9Y6zPC9XpENZGz y18OmRyY+DcS1/fLB8ruLyYqoELL7MvHcf9wohPcY0sVkBeTkPxmbB/9ofkN nQpIJqh+7hruN6TYNppuFlZEn0HTB7QIKLoZ9Ey+VBEZHF8OWtEg4M4HzXH3 DiWUUCgSuUMVj9e3ju+eoAqKk0mb8MX90eXnXTXaxiqobfTl+ogsfp/2FFG/ SwVJDSz5ENgWzKJXm2+roFLdmQmyDAHN16NVBo1U0acjf2QUtxHwYszMmn2X GlpbO2D/WBzXy10tEXtvaSC9uaBvBB8Bx9rD/U2fa6CVs6dcGnkJENSV9xH+ ooHowZ9/yWG7ciZbd+hpogzr83JfNuP+o8qO36hVE/07834mchMBr1kmX58+ ooVipRsU/TkIAKr/0T052ihpf9BvxioFzm8IqdzcpYvaauxEyJMUGG9I93v2 Vhct3GkYzJqggEpIqWj8N12k94b0u+8DBRjD7+LnePVQotVeL5/3FJip8XTv 99BDBkdDwg6MUcDiH5MfBYN6SEDezow2SIHpl2x6CpP6qOvVu4o23M+bnC+s N101RKzoyc/aBxQocRb8vmuzEQoTicpvwvsFQfYM5XhpI3TkYWJUO95ffI08 XtBqY4Tq6IY/XpdToJhKP+6caYT257fBlzIKbJIQtdwjZYzUqxYb529Q4G3R ubYEkgmqshcqCDlLgYTKtIHnSWbIZPxai0Aw3k8k77AS2EtCYafs1ovlcH// 1qNYONwK9csvv6jv9oKE44+zHattUFrgXbPJf7zA2Ki6dbLDDgVO9mbH/fYE ZLwx+YEAGfW33K6PPOoJz2g/Myd0nZAvs/mPnhkP4L/tfT2I7IyUvkTy7/X0 gPsG8hYKGa5ohYnJPv3aTvALt8iV7XdDZ8cSpHIm3OHSF/GXL5bcUf9yYugD Fneoiq0/IynqgZTz5MI1eNwgTbJwRk7eE+VYllRkirpA2phWbrOdF2rZzdui dNMJyA943gpSKCh5dj78KysZpuXdeszoFNRfVvjl/bIDpFzMbfXbTUEuffvf /vfVAVpPitytOUhBnUypLAVvHMDCRSbK8xgFFYXYBg7edQCdWT3+vIsU1Lsa 9oHf3QHEFXzsBIYpyKb5mYvCBXuYvXjnPt9eAulmfxvx4LQD/k4Wljw/AmUX PLhm9NsWTJh3EcIBBOrOP2fA/9kWzgRsYdoaRqB5iwvpxd22IG8S4S6fRCA5 7axgswu2sHvYfN6glEBMGxd+6cjYwn+SPcY+3wm0Pso7eFLdBhbdVbLezhOI nVpb2S1pA1Ipie/3LhNoMKeGi2+zDQR91U8/sIGK5jLCGv75ag1cTZdGggWp 6Jl3tYHvfWuw2ut38pQRFU3tT9D0U7WG2qLF1tvxVER6zx56StQKJBbfTQ0n UZFrWzTnJIcVJDh1cPOlUpFHzR11w0UELvPX3Y6fpaJBl10+twcQfLR1HHEt oqKQWKfnd84iEP1U+HOphYpsvR7GSLAjiNG0knXioKEu5qj/up9awniimm38 RhrioNYleeRbgt0bYf/aLTT05L+wE8WhlsB/arZ8mygNDY+LGd6StYS7vXnG 88o05ETUvBBLsoCRo1NuVxxpaHsa2xtWM3OwfJJx8lsmDT2XOsXmHWUCN9HI 1I1cGippbB1gsjEB7meqbp75NLQm9aTVYIsJ9Pe8kH54nYaKk60MHG4aw6Ep zuaEahq6FCEaUNFuBHn8yazCgzT0sPzC5jOLBrCc3x/UPkpDSQcXvxQ/MoDd kvIDMe9pKDz/gqdArAFsV265NTZLQyOXR/lf/dGHRnNm+7J1GlIpHZ02XtWD af/YVDN5Oirdy3My+rkOOH/tmPuqTEeSxnvEbeN0oCpcgn5DnY7WLJ7V/qet A7Fxj1XYDelodGnkxMc8beDPW37ZTaYjinvG5U8OWhAlRtZLcKGjGRvrSu0P mjB69WKhngcdFcw2RJRpaUJpmeHhAh860hiuS3jRqA7mTUc37w+ho9vD7H++ XlGBIpt/I4Ui6Kj+rQOPc9t24HwhOPYsio7eWXYM200pQ29/dYVaAh1Z7Xxm Wr5FCYzoLKJjKdgbah2XBBXhytjO+LPpdORRZd0Vt1kB/Gfm3H7l0VFcQtO+ mFey0BmCHpYW0FF0PBdKuiYDOgs5Mt5X6Sji3VM9d3dpuBgznspTTEf+e4Tc M19KwtofzR9NJXTEt/kdmadEHPYln6SH3aOjIcOpzwVCItDO3d0s/4COSLLd Rg8H+UE9R0r1dS0deX8d2RFRxg3/dx6M/v88+H+4j7sB "]]}}, {{}, {}, {RGBColor[1, 0, 0], LineBox[CompressedData[" 1:eJw122k4Ve/3MHCEkEwpQ2SWeZaZs8zjcTjzLopGUob6pkEkMw2ERERpkoQM JSqrCRVFaVJkSooQ0kD99++6nscb1+fax95n3/da617rBZXgcP/NfDw8PLsX 8fD877fR3cYLabcI4Pl/P3ovFTYnk36qHPfsCV8L5d4HSmkM6dxHlXxtfK8p 3TY6G8JIb/7D47jMsZ9ybPD9aSbpZSG7n93mG6Gw1weyLElnpXvfD74/Rnnb UXdcmnTc0RNreRynKH8+fXX70kDAkwc+ql6qM5SBKZG4BtJ7Ll/iy+Kbo0x0 rDI/THrp6kiPzoFflJrvZqHOpPmOxDgtuT9Pmd+ybhUPacfyIjv70n8UfHzB r+4mAVsFmfdqgRe0/dX4N5E2OtJk8qiaDw49+mwuSjpGOWrFCxV+6Ek8J+Nx g4BLv+I6P/AuhuVej1o+1BMwl3Vj9l2kEJz6kTm3nbT9knnNrn5hyA/IeLW7 jgD+wn2eFSgK5+3mkiZqCeAtLErLNBaDJd6H6zaSPrLvfHjYWXE4VBqXa1lD rt8blc3JFEnIvGntmn+dgINnv7hs+ykJX7MvRk9XE1CisarFpUoK4gtMA3Kq CBDjeTz3XUkaRIhJi9QKAsYbcvUl/62APXwz32uvEpAyvVT/Sb0MRCTKGPWU k+/nJ7c1bqcs8F3INpa+QgD1zIo33R/k4GjjcbblRQJuTlo/Nm9SgAcJ9drG FwiYNem61bRLEWoZnXFq50kPxR6y1V0FNzS/3pw5S8Ab9X1vNAqUYOwgZYtH EQHp1w6GtduowNKoo+vFCglQtV/VGpmnAolFi/ueFhDwK+Zej+S0CrA070Su OUWAc+Rwp/MVVRAT1ZxoPEFA1+ODabtk1eGxWOHJzGQCkmr4yz5Oa8IuVf3C x4kErJU5lDbEWg2r3aeXLBwmwFV0oam/YTWIiRjrecURkE1tKm+J04LiU9tj j0QT0G2yKNN+iQ50WNz7XLqFfP8I2Zt3lfUhN8z55aZNBJzIjE3xStCHoR8O pquCCais6RnuGNYHvXYH6v4AAqqOtPQ2XjEA0yL1JW/oBNgIJb5aLm8EKqsv 2uk4kPEns1biwjlj+DB1oibKloBtOiOvS14YQ4vIlss1VgTsj1uVnMdvAhfK hiWVzAiocMkI/G+LCXx3ff77tBYBPwXzoqe0TaExPKXPQIIAiQF/yb5KM3CX c7S/+5YLDRPKEQY3LGBBKXo9+xUXPqm7v7casgCb72vER7u4kDZtTHeQtISd NQ8Ufz3hgp7OFjPL7ZYQeC5K/OltLmwVoKd/XGUFI3+0l9qXcKFkX2ADb5I1 bLgUKsjdwIUTJ+6rWHrbwQFzVrP6Oi7o24k7h++1g8PVLqu/sLmg5TsjWnze DlaOhfAE+3Lh4MzeLZ/n7SB4rHpBwJ4Lr7Q4/CoV9qDNbJJVlueCVFecloQg BWJePXx0vp0Dmn+ObhAIBogZW9iIrRy4TfnmaRYN8JoTJfHmPgdMNxcKrD8C YCTbcnq6gQOTUqdeX6oHsF0cUTR4gQOdpvYbBJc4wvLDw2svx3DAp9Ig//t1 RxB/VZLarsGBppfbUzf9cwJOZnWEgTIHSqmeTYuknaExcEtmmjwHTLy1lM9o OYP4gGaMgTgHTqipbXrg5wzB5210TX+yITRfhv6q1BlSO/8oarexYV/I9KEl bi6gpSdF69zChiucxD+n0lwhw1/M4vt6Nmhlty51LXaF8Dz/ZFEuG+bXXr0z VusK9syDm/W92LCDNTKo+dEVzMLdA5QN2XDzwtksnzVu0PlPd/LpLAtSfQd0 1/e7wUtbgwlqLAvW/un3CjHxgM4nDS0le1jwW8g03t7NA8aL13qO7WRBhtS6 XNF1HiCvm8uKWM8C2R13HbOTPOC8qnarHoUF8GFhqecbD1hSN/VvLw8LzneU UjkxnmC4zr14ZwwTpjq9tZ7d8YLIiLz69N1MUPxrQQvq9AJJxfOpZ8OYUNx4 99zYoBcUxt/edHcdEzq7dtR9E/KGjhPlpk22TLDzVJEq8feG20uSPjX+YcCU bVBdzZA3qMwm1fT/x4CP8R+C4hZRQcPC7B9nBwMu2kvpGiynAuu0m+uTTQxg HPsvs1uTCsLlgs9K6AzQ/hyxSsyTCuY9QoISxgx4WO3QoZFJBesaixGer3Qw XNjR+FXeFzb8PXbGZi0dRAXLxYnlNGgLP+/u70eHL9Ejq0NX0aBV4Mf6jW50 +Pn1XP5/q2lQ9OJh/25TOszIacXFWtFAUarxMHMJHTTuVgpzAmhAjfJ5IH7L HxKTPXpWn6dBxpDk15hl/rCRJ65Cy8AP/jbe/fxJyB8cqA95i9f4gat66HWv v37wRaEmUMLBD/7kTdcJjPpBz8ZWrSFfP6BRh1/Y3PEDT3OeW76RfuDSl6Bc s9UPalzy681q/KD99+bpwXoaZIZ/zdEw9YdNMpF9/eU0qPO/OPLC2h9q7cu1 ekpoYPW2WOugoz/MNE713k2nwe6Kx2HNfv6w+PJ3I/9AGgher74nHOEPdT6e ykX8NGAzeJ9dKPcH9wqO2YPTVLA8MCI4u4oOfta3QksTqdCivFthQYMOz10z 1ffvoMKpwne9//TocOzgelsxByowLEODZ63psFrxv+rOfh+4LrnWoIRFh4Gy qA1PNHyA0l2iN32EDktPiOONS14gvPv8lfhZOvy33qTYNMsLNFYeyd35hw6B lLbOy/tJxzReYfEyYPciTc84Hy8I0bShyi5lwJ/HN5U/fPeE0r5F/kHqDFjL CjYNsvWEMwkPpYT9GXB/+7rq6sfucHI2e+RCGQMmee63B9e4Q1fB3h2z1xhQ qKw7ubTQHZI+mH13qGXAZt2TAf473MEvbsX0wzsM2MTvtC9H3B18+L5fy33B gKCWP+V8dDeQ6774w3GejKuEO4sru11g2xVzp6ceTOiJfKFyodEFXtns+vuB yoRklvDS3LMuZB16UDZKZ0JA2eHdwTtc4AbP8ONvAUxIzL7BruV3ASHuWvGq SCYM/NWbtTZxBgXuxW1LTzHhtf3dZUJpjnA6+9CM9kfyfhsSn7vudATO37dL ooeYsDA8Rj1Ed4RrnBufb38mr8uoVg4qknWsXOWb5RQT6jMedoVdBzjSKSU4 xMuCn/EznOPvKJDxWH/jUlUWRMue2n2Sag9Ho7x+nA9gAbEqpi9xpT082ib7 fWwDC67le77Y/tkOPt7bGmS4iawLp1ipaoftYDaUnlcayoKszMAMmzpbcNGU 9vSIZsHogSvScvI2YGYerph2nAVnihP4Tn6yhue9odPRJ1igYP/SWbTGGv5q lxgF5bLgqMafX0Ne1mA30v1B9TQLHu7MPeJ+0AqqTh6fCL7IAr/U/rTBXgvQ cH7L7G5kgcZxneY/VywgcKDHK/QOC26GK8mKRVtAZIHJmV/NLPj14dgGFXEL 0HxCZy5+xCLzaGPTIoc1sLJm7k/vMxaIx4R+vVdkBq3j30tjB8j3qz3pnhBi Bo/UFFc0DLHApFjmpL25Gci+W9kz/okF1McT5hefmsI60W1rXb+y4JKc1oz1 HxO47e6kWTNN3v+TdGsS2xjSawcjDyxiw+JIzvAidWOYMyqL8hNgQ9aS+bGD k0YQFLNyTmUxG8SKDo1tTjMCAVGftOsibDhqxisstscQipVXtmRJsqGYFZv/ RMwQbHj/e2uzjA2Si4Lux18ygK3tOmEfpdlgr/rU5c0bfRjOk5eUk2XDWGN7 vZGNHkiGbaW6r2KD+ve/9+te6ELAzDTnkRIbun6VG5qF6ULb+r6P9ips2Jun lqRYpAMjnju2Kqmz4dyVpJ7zf8k+ZmLE5Yw2G96seNh1ulkDHGqaeLJN2TA1 UE3c42jA0dNCt26aseFLWejTwUl12Hyq7usbczbsTnB8LqOiDh92qtwQsmTD Cv0PYm7xqpB4S2JEw5YNFtNRBW5yqvDaSm+dkh0bREXmup2qVaCCftltmT0b Utsn7hr2K8Nb7Z5DXxzYIGH1QLcqfhWslFDbw3JiQ/6847XIakWQf218Ss2Z Db/ytijq9SvAM+NZg6+kV6xXTT1CWQkf6983hLiyYbqQd8QkQh4MGi5UKLux wfbusl9dxXKg9S9seRfpmvLwBZ5/MvD2OJWi7sGG5HRR+nEDGYh1kV30kLTb +sOqsoErAKUdrDZ4sqGzv2JW7rY0JAZJ/ksiz91lergue2wZqPXuOSbhzYbz R7fxCyksg+qoisxc0ptzbN4P7peEVzslxo/5sMGXnfXl+CkJ0JbYPipKZcNX Y81zO8PEwcJVY10i6YqkuDO1l5bC6DFFl1nSTbMGHt8HlkCFndn59b5suKNm oaK1SgTiD3jtu0/aYLDwAosrBP9c3R4o09gwmHSmsvgZPyRVJDW2kJ5q8TFt FFkEvFH5WyX92HBIMn59pwsvyMmYpTFIpxod2fKQ9x/lWYujchbpsB722r7o P5SRxFq9FtIyN5NOT43PUUbU91/7QZqZaKmxsHGGosJIvqDkzwbjvpwffO8m KbU57ZKOpJ8LHF4pVPKFkmnvOh1Aesvv1k6+I4OU6LjPHrtI/youlp2ceEPh XXtt+WHSbRIXFC8df0gpupLNSif9WTxceANfeXOO4knRI6Qr/RvOhKZ2Ngta 1VmkkLbG1X6JVX3NY+zvrw+Qfu88etxz7FOzXDltPIT0gXsiad4835otrLqj /Ui/E365RC39e3NLb3K0CemsSw3F8tI/mpWoO8dESRs2SvCJkpVik/jJVx/J 96s/XGf2Y/VC87FDImuukb6udF9hjMODgd1vRf4jPX7HOO2NIh+ul5FmmpNe aXjq492BRfhB44n0BLne4jtGS89eEsCBOinPc6TXCMozDoYtRlN/nlkqaUdb 5Sh/Y2F0WVu+fJbcv0daM7HKP0TQZrH+1RzS4v3H/3y6JYo9QSU39f+3v/kT 2y7EiWG6t6pNMxkP84NC5/VVJdAxttPUm7Tww87K7k+SuPmOTSqNjC8RHzEf 1TYp3LDFbqiVjL+2KvMLIeXLMCT1V7016VPUwNSR8OW4Kd8RRcj4DfBS1c7/ KYuqhp93rHNng5Ldjau17+QwWzI+sYTMD+fa3InWJnn88cJ1WS+ZT3asw3Xv Dilg6rxZgYsLuf+iO7Sp/MqokibIeARs6M+qj7D5TxnnnDaZvKKwYanncnn1 T8pom/AqvY/M5883g6RHWlTwlZZ4eS+Z/5Sn14yMMtQwtq676pgVG4z2le8L mVFDqucznViyfrSdnh8oClDHdQGjP7dasCHmfnflnKEGempnntAh682LIzYb w15qYp7EzvvbjdhQeK2874DdamwoXP1Jn+yDBez5K5MursbDf9zWf9FnQ9Vc Z1/KXi2UvPCb6a9Lrtd/I41OCjpYYa8W2KvBhoty+bxqSTqY7SZgt4ush//e K7X/GdfBVb/HdvOpseGsJuV4/l1dfCRgMSWuzIbMzwbTh4P1saf4k8aI3P/2 z3bHAboRGtlgd7IoG+CUCqU6xQi18qyjkKzv1RUHhAcbjZCpzbPzhxAbzFZm e1urG+OqZuMkX/J8eHO1yLB+xhjV6a32dxZYwBMq5PYmxxQjz0q4HxtnQVpl 1HBY5xo0VVzvMt/KAgGu0N8AYQs0EPELppLn24YblvvdwQKb2Ao/8u+zwHTE JkC42gKZbWuTFMnzsav0p6hVpiW6jlRyXl1ngbqh5TJHH2u8Ubk+XZw8b4Me rBnCRGtc0a6+51ceC9gneeWsm6zRUFr80/scFuQ0+e1cpmuDHmtHmJnHWDD+ vUpkm5AtHklSOF4Wz4Ka4HStDLTDtVkBVhe3sKBfSORh2k87PCl29ZXWRhZI xNS3xBna4942t8UXyLlizC++2a/IHkVrIzXTOSyyH+x5lLHXAU3OwL7Fniw4 8t9e0wsnKfg+fvu6Z7osoBmchHtujgjzT7pPf2GCmm3IoE6wI14svCu59RMT soeLw4/FOGKS8UoHvQEmDLcs1nGtdsRVLxqizr9lAl/OXSJa3gk57+4Gercy 4bvxQdq7r06INXc5B84z4dywBm9Eugv2GzmcC+QwYc2B9yOrz7sgQasKDyX7 tduHfbre3nZB76NG7Aiyn7u+dTpXd8oFH6T17wt1ZkJDfqNkMtsVcfCHmawR EzI111yVUHPDvrT2O4QgE3xapWhede64I/HNumSy39xseCY49ak7Zp7M2GhN 9qNcp38H7w66Y7XTp8jRUgYUMMZvrpTywB2a1Vcs8xmQ77SnNGOnB4YO1led TmAA9cVkuLOmJ3oOW5j+xyb7291XiOlML9zR9PC/OLJf9vt2tmvLJS80e1bF l0f206Evz5S+uO2FzMl998om6KCTv8Im/4sXsiLXfmwepMPpFZyyF07eWBR3 /HbdEzr8q/V4wTPrjQJnlJ0aC+hwyEPdXMafikYSff0O5nTQ/uGbFBtCxWwq 1pka0OG4KMf6wyEqbt+iqaC6mg6vZzNuHqqkYunC/svfZOlweeuU93pRX6zN qJpznfeH4EoXuf8e+CLnwkbLH+gPbxd2TH45ScOgeyNdAR7+MNwYdFLsHA0v 0mpmLoE/vM/dmKJfQUP/tqmvY1b+MB4xbxdwn4bWS9I2btXxh4NJDRcSJmiY tvF27aol/pCrm3Rh1NUP5QMGn5964gcqPnF9/N/9UMCtsPWxOzlHHRFJePPH D1cRSgw+8AOdVQU3Lgn4o4NX/RlzKz8o+S3/2FLOH1ued+49ru0HvX4mfhoU f9yq/2aPhIgf3Mhuf6Zx1B/tRPyoHo9pcGuRjkitGh13neI4nnejweQpmtms Hh0Li1fEsh1owL/d6bzxGjo6e+28KWhBg3b6Lu9CdzrufbDKmKtJg5zCxFjL HXScatYK6yDnromKI5HpdXRckrbr+ly8L9iuLpdKcmRg4fp9aQNhVKAmn9s2 6MVA9+N1h99wqdC/t7DXhslAEbGR6w9dqcDb1KLUs5WBZ35DR5IyFegPL5d3 HWHgLG3pkaIXPvD595mc0m4GXqoaEjCy9IHkSgy8FczEq1KXEoT+eMHF9IZT 7O1M1E5dYz/5yQt0qHHF33Yxse2D8fCzLi+4S0kUFklioqNIk3j0FS/gT3IV +XGRic2Nz/Zv5npBfLPx6ZJRJpb/XtrRUucJuV/K7Q1DWbhIdJfQ020eoNy1 +N/7SBZ+9lHmUWJ4wJ92i5jEfSx83nXhU6iDB3z+kH6+OZWFp9pZeePLPSBo X3Hk+EUWmjPTtYvvuYOdE1dLZoCFPN3B2jflycqX+KGxkM7G79u/Nn556Apy DhXLYwk2PtdeGNpS5QoFs9u+coLYyLcoouF9gSv8d+Uf35+dbDw/lBhTF+4K 2lnHKybT2BhRrCytIOcKh0I2+ZTcYWNc696X90JdwG7/MC+/GgeZfjY+Z4Wd oaz/bnqZFgezMrqMPsw4wZf4J+YeBhzc3/nKRvqjEzBjeMf3WHHQ6Df7/a56 JxiqFLiV4ctB08g3L38FO8Hxf0vO5+/n4AbBaGmNRkdYNHu4yu4pB3c7Wh1T 2Aagljk/q97JQfmQPtsWf4CUgdtmgq84GB1yrn6HHYB+09r8+j4O/q21Vb4m BfD1jklvx3cOXvbRDVhookDAZDbjvQwXhe8cEH0nQQGer09NNwZw8Qm191nQ djs44Tgb9TGIi/4pDBVJWzuIdaXvYm3hYvTBOPmbonawZ3/Bdb1wLuZkC1iP V9iCYZGR0OFDXMyseXF5dsIGHv6nttBxlouXZCU1MqKsoTgtuPH1BfJ+oTe9 Gx2tgbgdsfptGRc5nC2zg1LWcPVpCP+jai42zkbzKNVYQXol9bobclH2MTPG btoShDMefh7t5SLLvEdIK8oCtuWx/8UMcFHH2OjqJ7AAJ42Q60KfuBgYyHv1 jKQFhMTwfBEe5+LdQuufv6vWAJtpkFDzm4tKtkEuWybMIVZ4NOmcNIG9ZpZP lm03g56dXnGjjgR6xsld0A00BikXnmtBLgTm2H0uijAyBiGKFf2FG4F6l+/6 VvMZQ/iRx3/PeRMYs1JKVPOSEQzwPG6bZhJoI5Xs3dRnCI9csiqJrQT2G6y6 qcGrD43SSryhIQQOO7kcML6gB4JNVs1R2wlM2HFdcY2HHpR8kkiKCCeQj+/b Wp0sXeDNjl+pGE3gOyUzB1TWgU3OpVLWSQROOjdVNdmvBkXdV4s+JhNY3Bqr GTSgCTutYW1cKoGzSoe//EvSBD1JdZsrGQQashwa9No14AZnTW1VFoHa0kdn Xdaqg/Fdy6cq2QTaMlyC7v5VA82ULW5Hcgg8luDxwvScGqTIfNjsk0dg/TW6 ptioKnSFNVsdKiQw8F/C0ow9KnDvaGFQUxGB4ZWXjPrkVEA3j+fL1BnyeRcj pQxvK0P0tP5Kj7MEDhVa77y3SBmCpOXmLl0g8G+3n9xuFUWoM1FfWXmRQI8f f5kftirAGcuDRZWXCGS9M4t3uLYSgvxXNhWWkevlMb92wloedh/fTayuIDC4 tcSIjyEDT8QuErOkVz5X+fuhaAV8eLrjV9M1AmuWzwXXjiwH5ZN8UlZVBO6S GSjxO0DOZbd2Z4yTvtem0Cn7cBmo01rOna4m0Dpr1rVHbBlIPt12Zeg6GQ/8 b6r9z0mCIF9IQVwNud7Xxf8uHpOAU+5XV0vXkt93M9g3mEuA/K5wOZ06Av9l Tmcpn1gKrN6rB8tJP48K6w7TEIUn2y7HatYT6GZaFXD9pggU39jMEr5B4Jq8 +HydvsUwZDOlF0n68umUPwFRguCwvqO4izRLxuNfhoAApL0YvqZ/k8D7+QLz tacWwfNOFeIw6YJPC+5vdfkgQH7n2eekoz7usvt1hweqZ24kyjYQODLCCNI5 /5eSZ/htEUF6WuDJNPH4D0VehVcpl3TpmRcfEyd/UrIuDnS3kZ5s//Xq8oof lLPqyVo/SS8oPE5ssZ2mOKR8XKF8i1x/Vtep/uBJynX4XAKkPe45JkxIfqXk U3MfrSPtq3ps7fLAYcqoYX9KJGllEZphGKOP8u3bo6E40qZ1m7Uvjr2g0BSc PqSQzuhxu1svjBSp5/TINNIh3duLMmQuNWeV/SxJJP1u+xG9A7lPm1OemYTt +581en2Sr79tXjzJ/2IraebCigj1owPN8ZRd3b6kd3yL9uCFz823DybsMia9 TGSS+WZmvPndV9sKUdJntxsqMF5ONfuyLh/uJ98v6XFXm8W1meZR8efTlaRf P6AnrUida773/S7P3v9d19MomQj63ey2IanMinSsfuP2ezYLzQZg+G2GXO+S Ww8VWkR4cJ9AR0cZ6bJryfpl0bwYen6XF4f0IZEPEilDfJg2rSp3jtxP8zD4 bnZHAMPWcELtSJ+Zbgd+3cVI3/eV3UXGh4/4uGBHnhAWTGcLjZLxIyb8Lcc3 cgkKzH18Ekp6Nx+3UrBXFGNP/NH9RMZfz+S2gXpPMcydk7vfRsZn0u32Ap4T Etj8MOi3Mem7plfoPu6SyGP+ryObjOfutPAvuX8lcdvTFzRnMt5pt5mL5Lcv w7I1Wxs5ZL4wYk1WjTqtwM8H/TfEkPlEbxw4kCgsg38FeRNOXSWQd9vFALln MriQdMXx1hUCda/MdRgQcggKN6+XkvmbJ6/A4ItQwKAMC58MMr/7se72RnNF fBJdv3X7eQI3JtZ/a/qtiITAL8Xl50hfam5lJinh37bvp8TI+rEnQu3CUh8V 3G+CtrdOE+h4VG9nQaUKfg/z3xhYQOaT7AdbVSlV9F7mZp1J1qNXT9xiVF+r 4oBzjXwBWc8CryocGNigjq/GaYJuiQTq37xi675rNUqdfCp04DCBzvaFG/d2 r8YXZpq5lw4ROHrcqvushRbeqU/w+RJDYOKZWe2Pf7Swdf6Z14r/CMx86p+z OFEH7YqezwttIlBDme45ekIfb9Xs7HsYRO7vvh96Z2f08cY+JeaB9WT98NY/ S2MZ4NUk3/lnBIFvrQOe58oZ4q6KA97KfgS6PHJU7/E1whiVHccDbMnv91Xp 2cpuY0x85/+4y4rA35THQeF8JqjdlpMHFmT+1wel3zE0wX7Tmq9LTMjn149F +KabIAysPUJfTeCm57jO2sEU63JPOFdKEGg2lr+gc9kMUwQNZ24vJc8v0XdN Jt1mONUwS3kkQqDa694N5nzmOOk3Vnufn8A0xgpd3XXmaNFkvzLqJ5fsH0uz X4uvwV2fTh8pIs/Xi5p10jp7LPBzW2aIaw8Xd9npz707Z4G5V8UmP7/m4vNv 8WLJzyzw+ZnpeflOLpqLzFU+1rLEu9ve65ve5+JP83cZoj2W2JEsUWF2kYvs 1HL5EHtrXFfpue/wVi66zsXHWvHY4aNImXDhjVy8tXdMp83QDncqVyxPC+Ti hsOiEbT1dsg9vGl/FJOLTTevH3W5Y4cnWjsfjTtycaZq8anWA/ZYZvHnVIoC F0v8U+Y1Zh1w59mHv/0ec3DJElZueASg+FOxuqkHHCzZ3NFemgK4MfeCT/pd DubcMfzddQZQ329O6WItB91C/yxTbQf86m/3X+oZDuJrnazD2o5YcCfP41cE 2b+NRQy49jvil/YN9nkSHKzi5NRd83ZGwt7MapcIBzUUrTc9CXZGlqRaqjs/ BwNSzoYO7HXGVfL62n1zbKyoMN7y77wzGrT7XW7qZeOnoyt2zMw7o+iGfjXr cjbGlBltfljugpr+Lje9HNhov/jHtTf8bii5Y/gK05KNbU8f5JjJu2Ftx4Qu 15iNNBdCMt3QDQ13C/Z7q7HR75beDRXCDZlZMz4jgmws/LVSYlGlG6LWrm+j T1k4tW5vVyrbHUt90xPMGCxcs5W/O7rUA5sn7d798mahmn+CdkG9ByYPXe1o cGHhgNxL9/o2D/TafidP24KF5w7t7e6Z8MB6Vf73rXIs/PJu4UKPjScm6P9s XNzLxKPLD6te6fJEF5Fr2dnkfPB36FtH8B8vHHcZWFVFMFHztHR53VJvlNam HW7xZ6KwcAjyKnvjareDbz85MtGWOZOe5OyNp0X6/46oMjHnhvcKm6PeqHNo 4fVsPwMnY82sCEUfdLWl5PqsY2Ayb2WQpyUV+Xy7j0vQGVg94jdo5kXFqUta D9s9yM/LyD5ZEUhFF9fML2ssGPict/hGcwIVTzBDQh9IMnCF9PLwqx1UVDBq ye54SEeTVVU1vcG+GNyWffS6Nh2Dgof/ODuSc2bVBpl7SnSUrzmYpOtFw1BB PdOny+n4V6tj51IGDWW3ius956VjyqbDrPubadjWdioq760//pH7ZjqXSsPM m+YvNVL98erXDz9jOmiY2FvkzOr3Q689Kw5PMP3we/w7m6jXfvjXcMVCZYAf 3jCUaktr98O260Z6YZv9cGgYKssb/DAsVrz49W7y85A+cSvLDzd3q5zbn+2H i98mhQqDHx67mrFn4pkfnolyfxNdSMOlxFdOios/7vqZvHHJCRpS/n6vv+VN zqFOjt0FKTQ8uGXq6CjdH/fxJu2piKIh9cONYqtgf+zPf/rzrDsNk1smlEtj /TFsdVRqzogvfuAflCyo98fX8687Gceo2F+06fISco5NanbuHt9LRamXsCuN XLd6wdn8uI1UpIy2VfIZ0TEtJOHAUXKfeitrokds6WhbstOAPuCDcgfkZPez 6Gj1sK0z3cwHrY4qu99LpeOobhzn2QsvvFK0RfDlKJ2sV/aJ6be9cM50u/Lt CfLvF6fLOFzywqsff8ecnaVjdul8Ws5+L0xVVP25lpeBBw8rmE4pe6Fba9Gz PFkGvo19CHvCPMl5PfyZtisDv35ZOv/yrzt6ncqprChioMFxk4Xnw+7onXR8 ILyUgRuD+AQfPXXHXvv/qHplDFy0s8rzdIE7Tmw8+Ti7loEL50U3T5u7Y7Jf 8WrFJwzs+WQvsmG7G3a/KlkuO8fAA9pa79d1uqD9znnF855M7L55/EjKDRc8 1P5f5BVfJr5LyNO6WuSCtOlz7VcYTNz87CBlMMQFtRIOFeQHMjG58fYtYT4X ZP4+uN4wiomXPoCFvbEz+gtd22x9ionbJ8tO7s1wxE1lL8we9jHRRGNCoT7K EdV/KQjiIBO9brurT3Ad8Y7cw/gbI0xc5BHUQddyxCQj1u6cCSZe6T9o9fkB 4AbafV2xf0yseJ3T4LhAwT0mVCkeRRbyaWVJdnXb47OW5dGvyTpRnH99yrnU Hnm4Rzo0OSwUoe4kqiLs8aswf23EWhamp/PPRSyxx7DI+xozQSycs4L6Aood Mj03H2vYycIC3Kp0uswGX9ddabFLYWFPl2vqxT2WmL6Bv0mzloXrPq0InLS3 xDcnNHP+1rOQXuvlaLrYEqmO41+eNbBwU8Hd/87nWeC2YkPJ9XdZeO+co/ua m2twJnyzjuFjFu47H1uk8dMMT2sVD+X2stA/hPt2cbMZNu929xToZ+G1z4Nq QylmmLD0qOrOQRa6R3A1MmTMkNi4WkvvMwtTgwMKc9aY4rb4B7u5UywkG6jd pbuMcfuqr4f38rFRatMhLQMbY+S+OXk/lp+N77eOjtXwGWOyYF1YLFlnT9/Z /7kqywjL43h+hIiw8Yb0lb6hQEOkTjp1T0mycZfclvfXxg1wk8onj/vL2Mij eHIhPMYAMeS06rHlbOwLFzd8k6+PgodN/4rLsfHXplU9Ti91ke+qMe2BEhv1 45VSHwfrYqFp9OsNKmxMMaUGeEzp4NJTo3d+qLKxh2bHZySmg8Efk66KabJx tl/+4Xk3LVyQDfjxTZeNtx+zyt0b1LHArOVzxBo2cgwf7hijqmORaKjucQvy XFr/QzFjSA1baj6cvUSeM4t3DlU3iKvhmjvK2i3WbLRkimT2blZBV/7W/GZ7 Nooqnr287o8yBuzO31VFnlPLpEsSXmUq467fpnEFFDb+6415dSBaCb/P9OgE OLIx74nsyyZVBZQ9/q290IWNAvpvF/b7rMTt/80aBrqy0cTelmWyVx7tDsnf k3NjY7/Ej61H22VxfKXAnjh3Nj6vf5Vh/FMGZYnIPG0P8n1u7ah5piqDETIe n9tJv40cHPgdvRzfBZ+SEfBiY9xfNjP1nDTO2NybzSO9L/5aiUT7Mlyv08Sn 4c1GpxKbPClVKUx7nXNZ34eNdzZIVRzxlkTHQ8ZQRtrVdUkkX7QE0qTkhZSo bJwX2FB9P28p6t8jlH+RlqYId+SNLcGJkJjIAF9yfy/mT20CEeQRPzjZSPrx 0TvJuieFMFPWpUiaxsbJcc2fX78I4vegO9FbSb93T6i66CCAp61fxtWS3ngs SX5dziI02bKr7g/pvn8fS0RHeTFZIV3Wzo+Na2Q8r92048FSYaGyvaRbe8YI VF9obvHp3XKNNA2spEQ6fjVHHPzn10s6aMtntnf0j+a1Z9eHCvmz0etbKG+a 8nRzypqZSj3Sl59ubo3Pnmi2MqhV9SJtKL+1fHT4c/Pf8bx7G0k/NL3H8rIc aL7qm398D+kRnXtfaMtfN286W5OWQNolu7FVZPf95uwLg1XppH/yp8Y6i1RQ fpxYLXCENO3L1Onc+eeULOm45BTSrwb1QnvZvZTbHl/NDpKuDmd2HX87TFn5 aKfEDtL3uE9+66SMUQqzlsqxSe+Ka4lfXzBF2a/VTrUhXbL/vbLG1AzF5kBt tRxpnX/GScNuPylXvz2B7+T7yjex3haf+UMZXKv47wHp8JF3bvTZvxQridrR LNKClpcP/tPgBdf2Ql4u6Td5GuJ+nXww0DrsKk/6SVbVeFEMPyhXnGnoJtf/ jY/KyOoXi+GlwsXVNqQ/njpjw79HGOapMv2fyP212VO2653sEljdqsR3lLSR uPd8VKAYNBhSD7SS8cEdH4MCHgm46lH6IoD0kYcRBnOdEnCN91v2OBlf+fqV N8/vloJF2zq958l4fKExNjXjsgwK/pW6HCC9/teXdIqMNNxZ/PHiDBm/f4bs l7U1LAc5HpnSV55slND++Tf/ryzcehmQzCbzw3Rw5sbdZ3LQ0uH7qIzMn4j8 Wfu+EnkI9XgR/oPMrwfr7L+JOCmAq4stTxyZf4aHhhniKUowI1ay0ofMTx61 vTz/dJWhOmhVRRSwcaf2891ROcrQJne17ASZz+pGxgyfzSrwTrVjCMn8V96w YPhcUA069gSHfCDrA9+d17c0wtUgPjKV+8qKjdf4s07vea0G2ruz6tvIerJw e/9v4cvqoNTpe7OYrD+BfpsjZ9w1wex4X7+QCbl+DSHfVKo1oV9iD6XbiMy3 WeVfHnKr4X2vrFihIRvnHpTR00ZXw88t1XOy+mz0CT0YdC9NG9aoNYo+X03G m+G2TJzSBuXLFmPhZP3jW6+fdIurA872DnYiGmwstt6zLl9bF+4S7uYGZL08 9/zcYqE2PYhT/Z1sqsBGVc9A9FhkBAHmv/jMxNlotX+9624bI8jMFDnGXErm n3q/ecEuI4hzyz4ZsYSNiZG/pd8MGoFzhJFq9mKy3pr0ycg/MIZsDSHT9H8s LFJ4xi5LMIXs/y42Ub6xkLtFRi6h0RSiFDdt/P2Vhcq0l27saVN4n8ufVTHK wmnDomvfgs2gezw1gWeYhd2bAj6/BnOosXZ8w+lhYUJGsUfy3zWgdCqj58Uj Fqp098X/M7eA3TpbHDUfsHBv/Uq9qB0WwDcRbBOFLKT1DOc5vbeA57ldn340 slAx7vGK3AZL2Fe8zqqmioUjf0OEhqOsIUi1OP7nKRYuITR6nvTZwZTh/bON m1goxg12vy5lD05zgbsOkuf5ngevfLNc7OHQuVP3rQJZePM6949luT2wXBQH z7JZmPj57Rb9PQ6wp9/W3MiThZm+lBe+SRQwuPf0o4sBC58ejY8qW+EIbt9T dq39zsRctVaHAkNH6D3iiju/MVGZ+So72d0RdqROlsV8YWJhw4MbfvsdQf71 qG/qABMXL092qv7gCON7vwuv72KidIr5ssxSJ/hw+saKJdVMfOm7cfk5XRew WiLUsymUiRNfnMd5XF3gKUf8W+hmJt5umy3mbnCB4eddXaEbmPikgRn9M8cF 3PXZMhwWE7ee3Cz3dcEFxN9b1H0GJgqERWcUdbjC1jaHrEEZJr5+5HScb7s7 fHpCHPO8y0B/3WGNY4nuMMRT1XK0gYFLOkIuSZ1xB933W/if1JD9Ydk5H8FO d+j0LT1mepmBDfYTnLw1HqDIWrf8diYDLe9e5yN4PEH4TGFIzQYGLrYx7zQ4 5gXZJ7f9ZzVPxzCeNV9YF70gopOnvpfsb1vFOLToO15gt69rayzZ/x61Vhst +eYFErl2mhUDdNwr55J1kuoNfNKelNZWOg7c4n4SWOoDwo48cT+y6ch3+faO RwlUuCAtpUHRpGOy42m/rwVU2PhIZqUFOdf0ekxVLr5OBcvIlXe0ZOn4rszr nuZHKszz/paaF6Yj86xq+oKNLxwS3G24cdwfrY6lSgp994VX+obTZjX++Pnl 9bDr5TT4Vzk8OWPtj5SgE0q/6miQM2PQ3Gbqj7LuEw9tmmmQJ7hLs0DPH2cO GA3WvqSB4eHBO/qr/LFIbJf8ngUaPBcOoSnz+uPvG74ME6of9DkXxJs98sPH +YvGD435gfMiapWGjx++qR30gx9+oGE90xXh4oclLpfWLfzzgy5h5qMbdn54 7WIosW2ZPxxJv69pZeCHCy3ydcLW/pDk+pNPRNwPmbwb30cn+8MaC4PO4Wc0 vGbL7J5aSQcieoNtNpWGkeGjeZ816CDWurj6igsNW1+s+vTWkA7FeSIJt21p OBwnL3/ViQ4LH3w9unVoKPQ1JEtsOx14uwIvtgvScFbGzKvtJh1Knok8ajng i0NP7Sr4aAxYbZa9NtKPiqNx95NcuAy4EhFzOd6eirI8604cDmZA+KM37DRd KrqZUey+7maA2M306VgBKna2H83ck88AN9+m+kUNPmiRs/dzVD8DbP2Ome5d 5YNDG9uPT4Yx4Yrf4G3TQS+84a/AHdrNBMHi/c3HnnvhYsbpjq4YJii5q/QM kHPUe6ZnffERJginaG/fc4qck1YUwq9yJrzUX7/Xy8cLGZ+/5paNMqFsx43A 6jpPVFrW8GVqAwsqk/V2HI33QL04swSLbSzoaM8+Kr3TA7u+/pcRHc6CCtui gFzCAz++kR36dJAF2xVnlRJNPTBO9FnC4QIWZP27ULeCnLPeiS6dKu9igYpR Uoi0mzt+S3gf1mjPhkDrfVtyBN1w0dz5czoubGiLYj5v/u6K4x1jnjlebHgc vrZ5pNcVFaU+jLE4bGC1v3bRuOGKbosavfIj2VAXeGj7+q2ueLLsvIZkKRtK 11r5vHzkguHPmVL3eTlQXB0nPBDvjAbbZ24dXMyB5N/esdt2OOOnpppLJks5 UKHiuOwLxxmP7VFzz5DlAPHcrLrX0BmJ4vnYOQMOSNRFuh774IQb99wvuEVw QK+vXm2JpRPm7CxNHKjggPtiXd8zo4AxtZFXfWs4cGm51X+tXYCK69WH6m6S z58tzx9vBDTuW31wx30OUA3C57SOAYY/KNdNe8MB26zDtzaZAL4XMzeO5OWC wuQy/h/7KOjvsOj3Lx8uvF+yTosg56yzSTqMZH8uxNvKqKa8t8OoJwezlrC5 sPnQ6ZvlFXbYonayimcDFwSfXZzs8bVD65rYgNxILvBW9S28zbbF+i73DVuz uXD5hLgiR94GDZ3thi50coFLSG1sULHEsLRRSbFuLuxu9q3QmLXAr9ctvoW/ 4UKa87GVR1otcNwgrFK5j3z+7erjbuEWmNKxw27NGBfuarOF9jetQUOFsldG AgQMd7a0v2KaY8aPE8dMhQioaZPa+0XLHNO3FX42XELAg23psz//mGH9R7VD yyUJoGb9zBU8Z4bpArNaOQoE+IwPlQyPm2KYuf2jQBMCdFcVJDQlmuDrvV+8 /pgRMO0anH6YbYLbztzblmVBwN2iR9LOOuT1TRGbr9gSYPvp0seGZ8b4K30n 97grAXMnh9Mj5Y1x10z8RXUuAfMT87fMyTmtLG4Xi1hLQKvshhLbRYbIjTHI SQsgYBkhOmlz2QBPOffVvwwi4Nxb99Xq3/WxCPPTFUMJMDGc2VSUpIfRYWLz e/cRYJT7sdi1XBvnpgbPxxwgYOdpqyZbF21c1XPEZu9BAmbHp9V1+7QwPXHS kB1Pfr/KefaUlBY+bbqmXptKgPl/E+Wa+zWxpvq2bnQ6ATf+dp0YWKaJZzeL K5gcIcAVTO7mVWigXmpQ/fHjBBw+L7f520d1jHA9dX80l4DYiO9Umpsaescn CIblEaBXu+TJ24+qeOv0HYNPp8j7r7VPDtivip4jXP17pwlwPtS/hlGhguvi NV7LniXg4YDsGO8yZcyQibMkzhHQPv4yDzlKKNj84kB2KQFbzz1q+y9jFa75 8P3i+AUCWDMfYu9NKqD7jEkS5QoBh+5ECM7cksP56NgxSjkBn0cXG4eNy+L7 bDFTq6sEPLt9/957JVl8sPH5OslrBAgecFpxOXEFuhtHOUyRPvmTy1x0Yzmu fGzJ01ZJwKSy9Q/2qDRue2G4ckM1ASF+yzjffJahpfHQbqXrBMRLhJ80OiSF euYPy1+RXnWPa7zjuiRSaEZVerWk39ZGv1ohgd90/A8+IX0yRBcWnxDDvybP NYPrCAjKfgqcH6LofFpXLLqeAEbgcZePd4TRIa/Wa5Y0q/Syk6yaENr8nNy0 /QYZv6zMT+4pgriiuYn7lvSMUWt51Fd+nNr+Q9vhJgGii15KnPRdhA4zsd2F pNcdt8quqeHFJ0+t1n0nfSv4UNUTGR5sHhK949BAQP2GjU9vz80337j48Xci aWBUbey5+6u5/mex1D3SN6v/qUyl/Gi2s7UUmCPNP7Gdj4c23exJze9Su0Xe r1CSR1h2spmTUrffnbSH/io9+ebR5vDG+IVNpK9sE4qTpw02K1jPEvtIi/91 69Vuftu8r0M4O4l0JP3FY/3mtub///9ycr47FIcb5O78H/ylBJg= "]]}}}, AspectRatio->Automatic, Axes->True, AxesOrigin->{0, 0}, ImageSize->{567., Automatic}, PlotRange->{{-7.006758589584157*^6, 7.006758589584157*^6}, {-6.113840864199892*^6, 6.113840864199892*^6}}, PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Scaled[0.02]}]], "Output", CellChangeTimes->{{3.559065395451902*^9, 3.559065415841137*^9}, 3.5590654478679934`*^9, {3.5590655658353987`*^9, 3.5590656079242725`*^9}, { 3.559065650902347*^9, 3.5590656733351865`*^9}, 3.5590657160012608`*^9, { 3.5590659936037445`*^9, 3.5590660084393706`*^9}, {3.5590660439294324`*^9, 3.55906608842071*^9}, {3.559108797683071*^9, 3.559108830396328*^9}, 3.5591155618704653`*^9, 3.559116631127529*^9, 3.559116696086043*^9, 3.559116783430595*^9, 3.559116924408041*^9, 3.559117016510601*^9, 3.5591170609862785`*^9}] }, Open ]] }, WindowSize->{1287, 827}, WindowMargins->{{175, Automatic}, {Automatic, 39}}, FrontEndVersion->"8.0 for Microsoft Windows (64-bit) (November 7, 2010)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[579, 22, 342, 9, 52, "Input"], Cell[924, 33, 75, 1, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[1036, 39, 1173, 38, 72, "Input"], Cell[2212, 79, 396, 10, 33, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[2645, 94, 847, 25, 67, "Input"], Cell[3495, 121, 1637, 49, 58, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[5169, 175, 1407, 43, 134, "Input"], Cell[6579, 220, 549, 15, 57, "Output"], Cell[7131, 237, 481, 13, 56, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[7649, 255, 621, 12, 52, "Input"], Cell[8273, 269, 795, 23, 51, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[9105, 297, 1477, 39, 96, "Input"], Cell[10585, 338, 1139, 32, 72, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[11761, 375, 1504, 41, 120, "Input"], Cell[13268, 418, 357, 5, 30, "Output"], Cell[13628, 425, 564, 14, 57, "Output"], Cell[14195, 441, 357, 5, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[14589, 451, 1432, 45, 102, "Input"], Cell[16024, 498, 847, 28, 113, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[16908, 531, 2122, 57, 92, "Input"], Cell[19033, 590, 875, 17, 33, "Output"], Cell[19911, 609, 593, 8, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[20541, 622, 4073, 117, 380, "Input"], Cell[24617, 741, 1107, 19, 39, "Message"], Cell[25727, 762, 551, 8, 30, "Output"], Cell[26281, 772, 548, 8, 30, "Output"], Cell[26832, 782, 537, 8, 30, "Output"], Cell[27372, 792, 553, 8, 30, "Output"], Cell[27928, 802, 553, 8, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[28518, 815, 1735, 42, 72, "Input"], Cell[30256, 859, 33743, 561, 510, "Output"] }, Open ]] } ] *) (* End of internal cache information *)