(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 10.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 158, 7] NotebookDataLength[ 59481, 1274] NotebookOptionsPosition[ 58516, 1241] NotebookOutlinePosition[ 58860, 1256] CellTagsIndexPosition[ 58817, 1253] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"interc", "[", RowBox[{"c_", ",", "a_", ",", "b_", ",", "t_"}], "]"}], ":=", RowBox[{"Which", "[", RowBox[{ RowBox[{"a", "\[LessEqual]", "t", "\[LessEqual]", "b"}], ",", RowBox[{"c", "[", "t", "]"}], ",", RowBox[{"t", "<", "a"}], ",", RowBox[{"c", "[", "a", "]"}], ",", RowBox[{"b", "<", "t"}], ",", RowBox[{"c", "[", "b", "]"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"intertangentc", "[", RowBox[{"c_", ",", "a_", ",", "b_", ",", "t_"}], "]"}], ":=", RowBox[{"Which", "[", RowBox[{ RowBox[{"a", "\[LessEqual]", "t", "\[LessEqual]", "b"}], ",", RowBox[{ RowBox[{"c", "'"}], "[", "t", "]"}], ",", RowBox[{"t", "<", "a"}], ",", RowBox[{ RowBox[{"c", "'"}], "[", "a", "]"}], ",", RowBox[{"b", "<", "t"}], ",", RowBox[{ RowBox[{"c", "'"}], "[", "b", "]"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"trac", "[", RowBox[{"c_", ",", "a_", ",", "b_", ",", "oppt___"}], "]"}], ":=", RowBox[{"ParametricPlot", "[", RowBox[{ RowBox[{"c", "[", "t", "]"}], ",", RowBox[{"{", RowBox[{"t", ",", "a", ",", "b"}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"(", RowBox[{"PlotRange", "/.", RowBox[{"{", "oppt", "}"}]}], ")"}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"tracp", "[", RowBox[{"c_", ",", "a_", ",", "b_", ",", "t_", ",", "oppt___"}], "]"}], ":=", RowBox[{"Which", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{"a", "\[Equal]", "t"}], ",", "\[IndentingNewLine]", RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{"Point", "[", RowBox[{"c", "[", "a", "]"}], "]"}], "}"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"a", "<", "t", "\[LessEqual]", "b"}], ",", "\[IndentingNewLine]", RowBox[{"ParametricPlot", "[", RowBox[{ RowBox[{"c", "[", "s", "]"}], ",", RowBox[{"{", RowBox[{"s", ",", "a", ",", "t"}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"(", RowBox[{"PlotRange", "/.", RowBox[{"{", "oppt", "}"}]}], ")"}]}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"b", "<", "t"}], ",", "\[IndentingNewLine]", RowBox[{"ParametricPlot", "[", RowBox[{ RowBox[{"c", "[", "s", "]"}], ",", RowBox[{"{", RowBox[{"s", ",", "a", ",", "b"}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"(", RowBox[{"PlotRange", "/.", RowBox[{"{", "oppt", "}"}]}], ")"}]}]}], "]"}]}], "\[IndentingNewLine]", "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"tracpp", "[", RowBox[{"c_", ",", "a_", ",", "b_", ",", "t_", ",", "oppt___"}], "]"}], ":=", RowBox[{"Block", "[", RowBox[{ RowBox[{"{", RowBox[{"default", ",", "Vector", ",", "Point"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"default", "=", RowBox[{"{", RowBox[{ RowBox[{"Vector", "\[Rule]", "False"}], ",", RowBox[{"Point", "\[Rule]", "True"}]}], "}"}]}], ";", "\[IndentingNewLine]", RowBox[{"opptgraph", "=", RowBox[{"Select", "[", RowBox[{ RowBox[{"{", "oppt", "}"}], ",", RowBox[{ RowBox[{"!", RowBox[{"MemberQ", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{ "#", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}], "&"}], "/@", "default"}], ",", RowBox[{ "#", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}]}], "]"}]}], "&"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"default", "=", RowBox[{"Join", "[", RowBox[{ RowBox[{"{", "oppt", "}"}], ",", RowBox[{"Select", "[", RowBox[{"default", ",", RowBox[{ RowBox[{"!", RowBox[{"MemberQ", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{ "#", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}], "&"}], "/@", RowBox[{"{", "oppt", "}"}]}], ",", RowBox[{ "#", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}]}], "]"}]}], "&"}]}], "]"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"Show", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{"Join", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{"{", RowBox[{"tracp", "[", RowBox[{"c", ",", "a", ",", "b", ",", "t", ",", RowBox[{"Sequence", "@@", "opptgraph"}]}], "]"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"Point", "/.", "default"}], ",", RowBox[{"{", RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{ RowBox[{"PointSize", "[", ".02", "]"}], ",", RowBox[{"Point", "[", RowBox[{"interc", "[", RowBox[{"c", ",", "a", ",", "b", ",", "t"}], "]"}], "]"}]}], "}"}], "]"}], "}"}], ",", RowBox[{"{", "}"}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"Vector", "/.", "default"}], ",", RowBox[{"{", RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{"Arrow", "[", RowBox[{"{", RowBox[{ RowBox[{"interc", "[", RowBox[{"c", ",", "a", ",", "b", ",", "t"}], "]"}], ",", RowBox[{ RowBox[{"interc", "[", RowBox[{"c", ",", "a", ",", "b", ",", "t"}], "]"}], "+", RowBox[{"intertangentc", "[", RowBox[{"c", ",", "a", ",", "b", ",", "t"}], "]"}]}]}], "}"}], "]"}], "}"}], "]"}], "}"}], ",", RowBox[{"{", "}"}]}], "]"}]}], "\[IndentingNewLine]", "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sequence", "@@", "opptgraph"}]}], "]"}]}]}], "\[IndentingNewLine]", "]"}]}], ";"}]}], "Input", CellChangeTimes->{{3.6887451704019775`*^9, 3.688745193722209*^9}, { 3.6887452299165144`*^9, 3.6887452327055326`*^9}, {3.688745285792244*^9, 3.688745288523926*^9}, {3.688745333900796*^9, 3.6887453478273993`*^9}, { 3.6887454100645123`*^9, 3.6887454116329803`*^9}, {3.688785852569197*^9, 3.688785902043835*^9}, {3.688785987817837*^9, 3.6887860075161695`*^9}, { 3.6887861453495083`*^9, 3.6887861604720097`*^9}, {3.6887862265847187`*^9, 3.688786287004265*^9}, {3.688786404367365*^9, 3.6887865079588823`*^9}, { 3.688798269346449*^9, 3.6887982786460476`*^9}, 3.688798499044963*^9, 3.688798552824096*^9, {3.6887988889994135`*^9, 3.688798935265626*^9}, { 3.6887989717961426`*^9, 3.68879901642657*^9}, {3.6887990469947166`*^9, 3.6887990588821025`*^9}, {3.6887991197529697`*^9, 3.6887991204709225`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"interc3D", "[", RowBox[{"c_", ",", "a_", ",", "b_", ",", "t_"}], "]"}], ":=", RowBox[{"Which", "[", RowBox[{ RowBox[{"a", "\[LessEqual]", "t", "\[LessEqual]", "b"}], ",", RowBox[{"c", "[", "t", "]"}], ",", RowBox[{"t", "<", "a"}], ",", RowBox[{"c", "[", "a", "]"}], ",", RowBox[{"b", "<", "t"}], ",", RowBox[{"c", "[", "b", "]"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"trac3D", "[", RowBox[{"c_", ",", "a_", ",", "b_", ",", "oppt___"}], "]"}], ":=", RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"c", "[", "t", "]"}], ",", RowBox[{"{", RowBox[{"t", ",", "a", ",", "b"}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"(", RowBox[{"PlotRange", "/.", RowBox[{"{", "oppt", "}"}]}], ")"}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"tracp3D", "[", RowBox[{"c_", ",", "a_", ",", "b_", ",", "t_", ",", "oppt___"}], "]"}], ":=", RowBox[{"Which", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{"a", "\[Equal]", "t"}], ",", "\[IndentingNewLine]", RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{"Point", "[", RowBox[{"c", "[", "a", "]"}], "]"}], "}"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"a", "<", "t", "\[LessEqual]", "b"}], ",", "\[IndentingNewLine]", RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"c", "[", "s", "]"}], ",", RowBox[{"{", RowBox[{"s", ",", "a", ",", "t"}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"(", RowBox[{"PlotRange", "/.", RowBox[{"{", "oppt", "}"}]}], ")"}]}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"b", "<", "t"}], ",", "\[IndentingNewLine]", RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"c", "[", "s", "]"}], ",", RowBox[{"{", RowBox[{"s", ",", "a", ",", "b"}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"(", RowBox[{"PlotRange", "/.", RowBox[{"{", "oppt", "}"}]}], ")"}]}]}], "]"}]}], "\[IndentingNewLine]", "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"tracpp3D", "[", RowBox[{"c_", ",", "a_", ",", "b_", ",", "t_", ",", "oppt___"}], "]"}], ":=", RowBox[{"Show", "[", RowBox[{ RowBox[{"{", "\[IndentingNewLine]", RowBox[{ RowBox[{"tracp3D", "[", RowBox[{"c", ",", "a", ",", "b", ",", "t", ",", "oppt"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"PointSize", "[", ".02", "]"}], ",", RowBox[{"Point", "[", RowBox[{"interc3D", "[", RowBox[{"c", ",", "a", ",", "b", ",", "t"}], "]"}], "]"}]}], "}"}], "]"}]}], "\[IndentingNewLine]", "}"}], ",", "\[IndentingNewLine]", "oppt"}], "]"}]}], ";"}]}], "Input", CellChangeTimes->{{3.6887451704019775`*^9, 3.688745193722209*^9}, { 3.6887452299165144`*^9, 3.6887452327055326`*^9}, {3.688745285792244*^9, 3.688745288523926*^9}, {3.688745333900796*^9, 3.6887453478273993`*^9}, { 3.6887454100645123`*^9, 3.6887454116329803`*^9}, {3.688755988061946*^9, 3.688756042550669*^9}, {3.6887568857885895`*^9, 3.6887568954671183`*^9}}], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", "\[IndentingNewLine]"}]], "Input", CellChangeTimes->{{3.688745109042972*^9, 3.688745109182702*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"c", "[", "t_", "]"}], ":=", RowBox[{ RowBox[{"Cos", "[", RowBox[{"2", "t"}], "]"}], RowBox[{"{", RowBox[{ RowBox[{"Cos", "[", "t", "]"}], ",", RowBox[{"Sin", "[", "t", "]"}]}], "}"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"c1", "[", "t_", "]"}], ":=", RowBox[{"c", "[", RowBox[{"-", "t"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"a", "=", "0"}], ";", RowBox[{"b", "=", RowBox[{"2", "\[Pi]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"a1", "=", "0"}], ";", RowBox[{"b1", "=", RowBox[{"2", "\[Pi]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"np", "=", "80"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"xi", "=", RowBox[{"-", "2"}]}], ";", RowBox[{"xf", "=", "2"}], ";", RowBox[{"yi", "=", RowBox[{"-", "2"}]}], ";", RowBox[{"yf", "=", "2"}], ";"}], "\[IndentingNewLine]", RowBox[{"Animate", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{"GraphicsGrid", "[", RowBox[{"{", RowBox[{"{", "\[IndentingNewLine]", RowBox[{ RowBox[{"tracpp", "[", RowBox[{"c", ",", "a", ",", "b", ",", "t", ",", RowBox[{"Axes", "\[Rule]", "True"}], ",", RowBox[{"AxesOrigin", "\[Rule]", RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}]}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"xi", ",", "xf"}], "}"}], ",", RowBox[{"{", RowBox[{"yi", ",", "yf"}], "}"}]}], "}"}]}], ",", RowBox[{"ImageSize", "\[Rule]", "Medium"}], ",", RowBox[{"Vector", "\[Rule]", "True"}]}], "]"}], "\[IndentingNewLine]", ",", "\[IndentingNewLine]", RowBox[{"tracpp", "[", RowBox[{"c1", ",", "a1", ",", "b1", ",", "t", ",", RowBox[{"Axes", "\[Rule]", "True"}], ",", RowBox[{"AxesOrigin", "\[Rule]", RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}]}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"xi", ",", "xf"}], "}"}], ",", RowBox[{"{", RowBox[{"yi", ",", "yf"}], "}"}]}], "}"}]}], ",", RowBox[{"ImageSize", "\[Rule]", "Medium"}], ",", RowBox[{"Vector", "\[Rule]", "True"}]}], "]"}]}], "\[IndentingNewLine]", "}"}], "}"}], "]"}], "\[IndentingNewLine]", ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"t", ",", "a", ",", "b", ",", FractionBox[ RowBox[{"b", "-", "a"}], "np"]}], "}"}]}], "]"}]}], "Input", CellChangeTimes->CompressedData[" 1:eJwdzVkow3EAB/BtxQMlDCHlepDSSlaYa25KkhRSGJucaWwrR38NzRWJac7F 9iAzx7+wIbleFFLMVfLggVzlWB4w8//+Hr59Xr59v4Fl9bkSDovFCmYCxy1S oanlOdFzqysPrtkyCuBH0mMRPLfElEJD2htxVK+QQM5EfhXpiWT1UBAuaoDl cnsT1HSYKeiUxO6EUS82old38DD0UPpr4FB2tBbWBd7o4G7IjpH0XUJoKBYs mOD3DNtMfi4Vu/Av274HtS7Nx1AZKrTA+V7qDi46Tz3DaU3sO7wak35BxZyv DR6M9nHM2NuhneCAv4VYGzPNhe8NDwFQeHgbCnnmsAio7uZGwpTZNgF82lSn QruEnwUrHDNLIF/uLYa5Rr77OmNvkIiYvkVboYr39AXD9PpfuGINsEFKtuqw wXg/uE+MujZItxmrxyg57FlerDliPE3Oq4N+2+1qOJD1MwKDDMVTcImboYWb 6k/uGeNnpqsPNFqFjVB8opJBXf9FE4yvKmyBrQlsJcyh4oh0pdsjnHxVEf8B qwUzrg== "]], Cell[BoxData[ TagBox[ StyleBox[ DynamicModuleBox[{$CellContext`t$$ = Rational[3, 2] Pi, Typeset`show$$ = True, Typeset`bookmarkList$$ = {}, Typeset`bookmarkMode$$ = "Menu", Typeset`animator$$, Typeset`animvar$$ = 1, Typeset`name$$ = "\"untitled\"", Typeset`specs$$ = {{ Hold[$CellContext`t$$], 0, 2 Pi, Rational[1, 40] Pi}}, Typeset`size$$ = {626., {193., 197.}}, Typeset`update$$ = 0, Typeset`initDone$$, Typeset`skipInitDone$$ = True, $CellContext`t$439609$$ = 0}, DynamicBox[Manipulate`ManipulateBoxes[ 1, StandardForm, "Variables" :> {$CellContext`t$$ = 0}, "ControllerVariables" :> { Hold[$CellContext`t$$, $CellContext`t$439609$$, 0]}, "OtherVariables" :> { Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$, Typeset`animator$$, Typeset`animvar$$, Typeset`name$$, Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$, Typeset`skipInitDone$$}, "Body" :> GraphicsGrid[{{ $CellContext`tracpp[$CellContext`c, $CellContext`a, $CellContext`b, \ $CellContext`t$$, Axes -> True, AxesOrigin -> {0, 0}, PlotRange -> {{$CellContext`xi, $CellContext`xf}, \ {$CellContext`yi, $CellContext`yf}}, ImageSize -> Medium, $CellContext`Vector -> True], $CellContext`tracpp[$CellContext`c1, $CellContext`a1, \ $CellContext`b1, $CellContext`t$$, Axes -> True, AxesOrigin -> {0, 0}, PlotRange -> {{$CellContext`xi, $CellContext`xf}, \ {$CellContext`yi, $CellContext`yf}}, ImageSize -> Medium, $CellContext`Vector -> True]}}], "Specifications" :> {{$CellContext`t$$, 0, 2 Pi, Rational[1, 40] Pi, AppearanceElements -> { "ProgressSlider", "PlayPauseButton", "FasterSlowerButtons", "DirectionButton"}}}, "Options" :> { ControlType -> Animator, AppearanceElements -> None, DefaultBaseStyle -> "Animate", DefaultLabelStyle -> "AnimateLabel", SynchronousUpdating -> True, ShrinkingDelay -> 10.}, "DefaultOptions" :> {}], ImageSizeCache->{677., {230., 237.}}, SingleEvaluation->True], Deinitialization:>None, DynamicModuleValues:>{}, SynchronousInitialization->True, UndoTrackedVariables:>{Typeset`show$$, Typeset`bookmarkMode$$}, UnsavedVariables:>{Typeset`initDone$$}, UntrackedVariables:>{Typeset`size$$}], "Animate", Deployed->True, StripOnInput->False], Manipulate`InterpretManipulate[1]]], "Output", CellChangeTimes->{{3.688799125389207*^9, 3.688799140989805*^9}, { 3.6887991858947287`*^9, 3.6887991953039813`*^9}, 3.6887992429430695`*^9, 3.6888000768910666`*^9}] }, Open ]], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", "\[IndentingNewLine]"}]], "Input", CellChangeTimes->{{3.6887582164452095`*^9, 3.688758216592845*^9}}], Cell[BoxData[ RowBox[{"ContourPlot", "[", RowBox[{ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox["x", "2"], "+", SuperscriptBox["y", "2"]}], ")"}], "3"], "\[Equal]", SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox["x", "2"], "-", SuperscriptBox["y", "2"]}], ")"}], "2"]}], ",", RowBox[{"{", RowBox[{"x", ",", RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"y", ",", RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", RowBox[{"PlotPoints", "\[Rule]", "100"}]}], "]"}]], "Input", CellChangeTimes->{{3.6887513217261114`*^9, 3.6887513682747145`*^9}, { 3.6887527099414587`*^9, 3.688752718039686*^9}, {3.6887534984303904`*^9, 3.688753526223835*^9}, {3.6887535942521625`*^9, 3.6887536176230927`*^9}, { 3.688753681874908*^9, 3.688753689042772*^9}}], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]"}]], "Input", CellChangeTimes->{{3.6887557494310913`*^9, 3.6887557498194704`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"c", "[", "t_", "]"}], ":=", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"Cos", "[", RowBox[{"2", "t"}], "]"}], RowBox[{"Cos", "[", "t", "]"}], RowBox[{"Cos", "[", "t", "]"}]}], ",", RowBox[{ RowBox[{"Cos", "[", RowBox[{"2", "t"}], "]"}], RowBox[{"Cos", "[", "t", "]"}], RowBox[{"Sin", "[", "t", "]"}]}], ",", RowBox[{ RowBox[{"Cos", "[", RowBox[{"2", "t"}], "]"}], RowBox[{"Sin", "[", "t", "]"}]}]}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"c1", "[", "t_", "]"}], ":=", RowBox[{"c", "[", SuperscriptBox["t", "2"], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"a", "=", "0"}], ";", RowBox[{"b", "=", RowBox[{"2", "\[Pi]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"a1", "=", "0"}], ";", RowBox[{"b1", "=", SuperscriptBox[ RowBox[{"(", RowBox[{"2", "\[Pi]"}], ")"}], FractionBox["1", "2"]]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"np", "=", "100"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"xi", "=", RowBox[{"-", "2"}]}], ";", RowBox[{"xf", "=", "2"}], ";", RowBox[{"yi", "=", RowBox[{"-", "2"}]}], ";", RowBox[{"yf", "=", "2"}], ";", RowBox[{"zi", "=", RowBox[{"-", "2"}]}], ";", RowBox[{"zf", "=", "2"}], ";"}], "\[IndentingNewLine]", RowBox[{"Animate", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{"GraphicsGrid", "[", RowBox[{"{", RowBox[{"{", "\[IndentingNewLine]", RowBox[{"tracpp3D", "[", RowBox[{"c", ",", "a", ",", "b", ",", "t", ",", RowBox[{"Axes", "\[Rule]", "True"}], ",", RowBox[{"AxesOrigin", "\[Rule]", RowBox[{"{", RowBox[{"0", ",", "0", ",", "0"}], "}"}]}], ",", RowBox[{"Boxed", "\[Rule]", "False"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"xi", ",", "xf"}], "}"}], ",", RowBox[{"{", RowBox[{"yi", ",", "yf"}], "}"}], ",", RowBox[{"{", RowBox[{"zi", ",", "zf"}], "}"}]}], "}"}]}], ",", RowBox[{"ImageSize", "->", "Full"}]}], "]"}], "\[IndentingNewLine]", "}"}], "}"}], "]"}], "\[IndentingNewLine]", ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"t", ",", "a", ",", "b", ",", FractionBox[ RowBox[{"b", "-", "a"}], "np"]}], "}"}]}], "]"}]}], "Input", CellChangeTimes->CompressedData[" 1:eJwd0U0ow3EcBvBZXorawRAODIe1WitZsf8cZt6LHZaiHBim1LTMNpmahrwW Qs1L1lgSRlbYH00bFyW3v1CidrCYg2lNycye3+Hpc3n6fg9PQYdWqWazWKy8 eOAq0ydzDwUrMj0TTfAkWtcCv+RvrfCOkbbD3ZpP4rLDqIbsteYe0lPptZAq VulgpyE2CK2jtBmmyhPGYNlHlJg1WbQIMyz5VrigkNhgb8HTJvTxvU7S5/Bd sIvad8OfjQSa/Lk3+uCfInYJbRzTLbQIZAzcmzb74UHaehDareUh+LDSF4HG ndwovF6eYdO453Wlwtl8hqiR2rkwpAvwoOzmWQBFtLAELk1yS2HV9jAF38+X qmFMLW6A3cn1bVBsyO6CSqc4/TTudKGKWOtxheG46D0ChQ7HLzwK86LQrD9O Oov7On9FPDt80VxgF2NVP6yf8g9AzWOWCSaO5DzBucqUAKSUkm+49dxI/Adm fAbK "]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"trac3D", "[", RowBox[{"c", ",", "a", ",", "b", ",", RowBox[{"PlotRange", "->", "All"}]}], "]"}]], "Input", CellChangeTimes->{{3.688756249348799*^9, 3.6887562641777716`*^9}, 3.6887563088746758`*^9, {3.6887563825811653`*^9, 3.6887563827007933`*^9}}], Cell[BoxData[ Graphics3DBox[{{}, {}, {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[2], Line3DBox[CompressedData[" 1:eJwsmgc01//3x1VKWpIRhUpDw0gpUVxt4duwSiEjFSlZFRLKSmSVVGTvTfa4 9vp8ePuQvcfH5mPv8f/8zvl/zuec1zmP8z6v+7z3fc993Xte70M6Rop6GxkY GFo3MTDQ/wzX1tbpPxqI5G2r7j1nJ33//9cMyyvBg3S++vJOLrP3YxhY3Ph9 J32VSqW6BNO5EyXT0MnLCL5IbCP50Nftj9OyhOicT1bDwSnuLbQbbDYLoq+U X8xnXdZo8Lgx75CRpw38+hDiG0tfu+dzZtNXaeApeWLoZdYncDWaHYuir74v zm9OX6HBbUqvUWOME9yMCj0MsU4gfS6X88QSDbZl8QY6urlDmkOX1vuv7mDC K3PQb4EGPkGfyQ15nlBo7ZTlgJ4geszMoWuOBlIF4yJhad4w4KWz/3u6N/zs /fu2fYoGK12xm5XDfgD3n32mU+E/YPJRXnTkIA2y6+/FHv4QAId/unMt2gVA a/w5qQgqDWrD3uzwywyERY96i8HsQAhnT9N520MDzssd/NSUIEhYPvsyIC0I ZmhfXV1baZB/Ly8tNCIEJN9rL2nHhECAhJrdFxINeLbuOG70JQKOnqzs6PGM gCzz0ob+YBpcj37Z2COTAFcGipQO3E0AjtttPIz+NKD1Gr2t3JcIV7vdjx05 nAjsQh0V0z40KBiMrer7lghX+n+Ybf6TCAsz7QFyLjRQtekZVXFOgv1nuvSO eyfBqVei8k2vaDAuwXNu7kUKWLmodRVbpkDLltvlDiI0eBKeYfPtSDrwZ8Qw MkqnA8PwnkP17uPQr9yklH47D0gNWurcLnmglJUQafppHJz+cIhvLcsD6/3S g883I8z+2vKo6804WCzryJClEYJ3vZ+9aoswEDoxckRjHES0o38rcuaDisyj 33KS+ZC3+ujHY4FxeBlpnpLAUAB7fl10zDhVANn3pyeLksegyX3kt2dpIRTP pG27tLEIZvEif2L4KJTMj51eZikFwSBml4S3pVD+4dA502+joDts+Uzzcils 7VrJUYothbuZ6hu2fxwFjii14kcmpYA3D/gEdZWCdXBmXqb6KFxU32tyqK4U huJqx87LlkGFq2Xmvt2jcDzL2zXfvQz27jWVE2Urh097g0mHjUbgMOuEsfdM OVj9+OaW+KMCOJpvvz23exjyRWL28z8jwYt3qgKGKmTwN2Rw+b00BEo2fH/H nEjwZqDyI4cVGbbOZ1j09g1BatIjU/lIEnzS0AHNQDIwewiMcmcOgXjoI93s ARJc7Y35QRoig5iOEOtxrSEQzJBVfaxNBlGee1xb3lUBddeniQ1Rg7Cxrifu DlTBu8BWmo1dNSSMJk6HeA4Cw9W2JUa1KmghRCXbA6sh+0d93ynLQbj6cUj4 ukkVZDU2SVGwGlw9A/6Myw3C+h+27CchVRCU2vnMaqUauCY52m+MDIC6kkot B0M13HB+fanyNQG+S5/iPwgMQOlwkPNkcDXkxK88NeGvAY9xpbCAXQOQuZ7O dTK1GhzFcu9kCtcAv83O1ujZfkhdKOcfLK0GP5am7DLJGnDtKuKzK+oH89hf X6YHq4F61lhCUrEGVmT214pp9kOLXqCBjAABFj89wo0+1MCPT9lV8+5UEPlN c/7gSsCekba+UVINRKTIy3aYU6Hjql/+pA8BGeRqBa9/NfA3QepTyiMq6Gg/ G+sPJMBx6GPPclsNFFjmPz4lQIXkVN3PqikEbHfx+bdvrAasQTL/QW4fBHk4 jN+oJSAqejLg4zYK3LwzcCs3uA8sHi32lTQR0KYaK32NlQK2nV3de5z7wN6/ t6Kkg4Cmq7e+le+lwOeIMD0vxT6I2Km6QXaIgBxnuRSuIxToniruHO7vhayw R5uVFwhIuF/PXCtBgRt3nw8/I/dC3Ro0kpcIOO3xH5ObNAXW9LXqmpJ6Yfe/ CkbSCgHBLAxN269SYPW3ULfUj17446/seHuNAL0DOqq3b1JAgkvz4Y/3vdBx fY/q/XUC7iS91lCWp8C1OVcHqnYvJMH211Q6F01LLOe9Q4HGK9pRx2/S+Xe+ 2gU6FzcVMYtWpIA0i9sObcFeUJja8N6TzjkaPl5aUqEAYUue/craC4uBl40S 6HZ5yivWdqhR4EJXX0zyXA+IOMskKK4SwG4b69vziAK0yBBzcmsPWPFcuGy1 TMB6rfqgpSY9br7tP9vye8BI4YLgiUUCQj/5kOu1KDCU+/Jxb1gPuOkJWKjP ETDrEXBkWocC685iO7tdeoB/quAMxzQBNz/J1tQ+oUBDM+dwvVEPCHj80lOi EfAkYzTB7CkFOPbXKoRL0J9PtjbR6idA6z8F97XnFJB3/qj1ka8HmNebX4l0 E8B4+q1ivz4Fjn5bzlPd1AOhL/aesG8loG74fY27AQUuSZxMODzYDVy+02ma 9QSoKpqOLtH5Q5mlt0PkbngSx3Aaq+nxyV/5LPCCAtZ3VbQik7phZGNYXFgZ ASqqsm4cdM60fihBy6cb9ubuk9+TT0DKlT2txfR9WLuuxOyx6obJJ0ZHN2fQ 83M95ME5Oh/1+f4773E3xL779OhDAgEs3/6NPKXrJPonS59c64Ybs8y7PoQT 8LZuq5Ua3a/0ouNmm090A2ez02NGfwJ+/To7tuMZ3a87V4zEp7pAI/WbSZgz AcrV61xpuhT6+UCE6Nl3gfhEcrHUQwIEjyZtYnlIz4evrT8HnndB2dvZqjf/ ETDyfnb/lvsUkDXvu/H0vy4IKP1NEpchIIa5TSdXiQKpzA6dDzm7gFmb1Vjx CAGH9u0UMlSg56F+ndm9u53w0GlMd6K/GqLTNoSXS1Lg3Ks9Qcd/tMMF2es5 jKrV8Daqx//Vdgosa33jPiDQAj7tJJZf7WQYUvgtxPW+BhZingdmjTTD2Rqd +A35ZPA6IcK0+XUNGDEfz9VJbAbH+4vuA0FkOJbZxx+mWwORCje1WiWa4VJm mdTpJ2TQv8nFvX6rBmwMX5fb/tcE/oc044t7SPBs15EFbrYaCNwceln4XQNc efHZ7X15JVS0QFarEAEhher7HTtr4eseY+sEyXLw/+4l61NJBgsxwrlTuBxs /a+0LvUVQInj1oe5V0rhDSlQl1pdBmvvDnpDYgF4bR5ovchcCrfHZvovG5XB G/GtJdJWBRDV4RwiTZSApJP3woakUtjw6G6z0W76PpbbpGvVSqBN7K1ByoUS ep34U50ing+skTfNJfWLIUJS6YWNWiGQtDu4VmdyIPpB8C0140LAv/+xbiNl QsArHQ21+FQ4lyr2KlsxB+wTuhLOwS9QadQYMFvzhY69FSMZl4NgA0tCqG7/ D4hNYbnW1OIDc+/K4swCAmC/bdDZ3V7f6fkcfrA27RvoSSjGz6T7AZxOVpqa 8wDdp8I9iS89gGr+UuJThw/EMt/zrWizgfBbls/ytGyA4ffF7mVGR9B6FbeJ dP8jVkdVJFyL/4jvje4X/o1ywZH7LtecFkIQFG1ON+SHoYMOC7s0RyL6zAh3 +LWE4vnHubFhRDj6zsaGqg8kYsG518GKWWHoyejZtr0jApG3XP1EVhLefHRQ 198kAm3czzQ5L0ahfH/D41LtFPx9eP/BqNJo7DXfdyn3ZDz6pAozd2SkoVe9 ocXVD4lYtKkmeuRCCjacfPxgwT4XnSxbX5nXpKJLV3hu7YVsPHPs+89dxSV4 adu0vuS+NDzRU3C0KDkb+/UGu4XESvHKg9a9nTppGH3z2L0LgjkombbUbR9a in9iGLUqJtJQ5bH2x0jeXBRgPPG4waEM39+pNRRZT8e5LaV5ekt5eOLMot2Y XAVWbmgMyuvPxEFWTt0ihQJ06D2/iaeDjJTBwNlMag7maYTsvJxYjAUXxbxG d9fi/dZYztgduZgduc/ds7cYlXXI1K/ytXj+t9b4tzO5qBpHOuDAWYJ2Kto9 Uo61aLTUnXXbKhcHV5Op96xK8MkXky5crkXit62pDUMeKijFBOjLlOLBjXWn bgzUoadC2b2yljzMXcu++C+9DC2ze5TGsB4HPLKYQyfzsNMzjZvWV4YTpIfM B5br0f5D4t7PTIieUVURAazlyHdgh7XK+Qb8Y3CMz/g04iv9+qw8/XIUFkgL qI5tQLkTB2da3yAKCHeZdLJX4LYL6TW9PxvRJ3qrDp8jokWElZWXdAVG9T14 K1TfiDtzSHdNvyH2X9A4U/+sApmOjrS83d2E3z3NtjyMRxzl8a3eklmBaw8i G3Y4NeHr7e87ZRoRZ3YInGBSrcSRCUO3E2bNeJPNnSG7A3HiUhXth0UlLg5b +RnFNONXmSMLRlTE9zyyFwr9KrH1o1Jjck8zbnKS3e0+gZh0IV+ku7sSPWQy hUTutWDpd4XPyzOI8nis7t8mEvrxDr595tyCFNGkzxELiKrfL3Q/PkrCh2X8 8r+xBSf3Kl79toxY6LWvye4GCZ8OJEhXzrVg947Lw4WriPueOotLPiPhxZa6 oWmhVny7jn/F1hE5HBpsnZ1IGOxO6d2r14r8hy+2DtP5JOfMecMIEi4YjBSd 82vF4I1Xffvp/GeTYUdfKQntDz859F9dKwp8cTx+gs4DlwxWJ6kknA218dbY 1oZdLFWVKXS7bQsZ+70ZyRhD+5T09HIbnrgxgU50nTUvigrL+Ml4cjFzz/N3 bXhfskEunO5XUt0va3cZMr5xfXVbO6EN1aVeJO6eRTyoMds0pEHG1yenl+/1 t6EI6bd4OT1urwJO6zZYkpFt+E+yBG87ChncPlE9gmivrBSr/IOMRQ8IZi7l dhQ552lR2YloU11QxEiQsWLbw/m0gnY08/I2j0JEf15rLQGJKnx0UVXLxK8D OeI2jSykIpp3u4r4KFahbZxE7Ka6DvQ1faEaFEPXGS5r6f2iCgunOpWdmTvx 4Q+Ztxu/0+P2yWHDBb8q3PvDdr/Bm07c1GXALaiF+Hv3oQS+hSqEw2Y7CwW6 sEXF2UD/HuIu5WCB+Z3VaLAlzqz2Shd6n7CI23YF8Wrt8xuPDldjE5focItG F2ryPgj7cAjx+CE5oRKFaiRdKT9d592FryrrFpSa83DIdveGvb+q8ajAm0+P 1rvwYY922qxQHoZZGgaa8hH47uaNw3/+daPgYZu3vJx5WJ4s6NErSODs1Fep vbRu5El7BPmruWiw+yXRIklg5lmKsANzD16oJP17T8pFjbQ1PU0VAj+2famU lu7Bz2Ifxpa16HVnov/9X3sCrxn2PDgc0YPGdg51uuY5GGQcHfGGTGBs6blg hze96HT5x9lrSjnY/PNNYl0dgd4jyqaaXr2oHvajOfR0DrLbXmgKaCHQmjXe 8Fh8L75jUasfH8rGXamG2x36CfyUd7L1I7UXtVbVedpKsnGneGaX/wiBtyrP 1HJu6MMBjfwv94Kycfu/xOdHJwjs4INZX54+dH2kTLSpZqOuap6M7TyBZ+bS ip8q9aHav4u/505no6yNmpHSEoF1p47PJL3qQy291g+B27IxuCd65/cVAsPk o1WGP/fhnGXe27beLNQVfMl6fY3AcOsDMSxhfZhxb5t9bE4W1tdOmuiuE7j5 iGfaofw+jHtuGrznexbGSBwSmqDzyLenJQ+29uH1qTO1u19mIXeFpdQInf93 5P7Ctrk+/HttZU/UtSz8b5E9VJnOr/Y5pXTvpqLjrakn7fuzcPKev64g3W7W fiPh0FNUnBj1K0+YykRWzQd2JnSdWaa9pxRvULHYtUOavyITub4d2ixA9yvQ rP/ViBYVJbefe04zy8TmndIvWqcJVPryrrbjOxWDS0z5bstlov5DCYNWGoGS b/stLiRS0ZXNeUz5QCbu0PLaKkeP85XhILP3lVTknTDs0CnPwKVjByxMuwh0 MkgxLVyj4i/tb4yv/DLwUn7ke2H6exQSv89QzNWPL9lOKpx4nYFz3ZIfH9Df u+QF6kT8mX486aYnm8GZgf8Nf9qzVERgVKb+bpmn/cigOiD4WisdZcs47r39 QyC7sHODLLkfZdoP+YidSUdVR2n1eG8Cn9uE7HKn9uOPdhvhzE3peP7Uo2gV ZwL3ue6pz1nrx2OBv8e3haehwTOfrlkjAplfzB6tPj2AoU85grcOpqK9e8Py xwsEqjzr9irwHsADk1UJJ3X/ou5Kvdm2hGpUKtpyIkt5EFtGLjVUXU1Cyqou 6zGxKhQ+IT9kUjCE35TWT5awJ2GprFBIO3cVFuRef1dVP4QuHZyy2dREJOpH Cg6tk/GErtS/HcNDOBzTPhvllIh26wp375WTkT3z0V2JPcMoJLd+5O2PBNyU dTZLTYWMN0nBk9e0h/Hi7sH21ZhYDN1UU2P+kIS7ruy9w7g0jN+k/sZoLoXj yzNuwnr0c03tyFLbPo5RjNNgCk8OCkdaN9dT/9FylB4mKxcfHcWtV89yat0K RxFOxxyT4nLsEWqg/Xd+FJfL7/RTftL7oQ0T1+ZNynH/Es+2DfdHkSqld/kh hKK7ioCnBakMabeMxXt9RvE6V3yKeEAQfte5oz5nUIq3Y+s+W+0ewy3xrgrb HH+ig4moTOrzInxQRn1wcngMB3wYfa3mfbFY409dz/kibPqvjVtzbgyp0e/M Thr4Yi1139EyxiJ0/Hrkr9HGcdSbT5e9q/QDz7rbr6cFFOLE2x5OwX3jOP8l i+nYue+4IRKzZGsL8G/bMTYm2XHM7ZXyOyrpidvkL013nczHO7cMU874j6Ol 9ONDJwo80OWIiGXsEGIz1ddHImIc5ZnbzoTc8sA9b7q0ea0RFZ48XeJJGsck 3mfuLHlfsZ3wCvcNzsOop0zz74rH8UysWUud+Be8EPUjk38wB83ffcj8NDSO Yjsno24X2CPZY6PogmYmpm2dmPgnSEOtrQeeLTLYY/aj4uKdoxlIuW5g+FuM hgPqqy5jVz9hcHecUvG7DEy3+6ssdYmG6eYdKyV1dkiWd72y3zMdN+sYnGaT p2Gv0pzCe/EPOCXeVxWYm4r1cw32m5/RkCtnw7BekBmq9nru82JKxi5nMdbF HzRc77VQEvtqivbqf647uyfh+uzCkpE/Dbc3sZYf+GiCh6tq31/cm4SxzFT7 gmAaBuaxEfo2r1Hg+N2LXUcTkVNx8s5MLA0jjj7dyRSkjxMtlrnRtrHYGbEx +QfS8P41sbn6m49QLSzRbmtjKHb4drBUt9NwuI6bS+XnA/zEfIJRnhSCrK8D 0gW7aWg5HDG8tFMVw6JT48UwGF8fYCjW7KPhaxpzPLP2PUxyCJTbHR6I1b1B Y5eGaSi8a6pR2OY28r8t/2/vmQD0YX6yOjhKw797Dw6e/CSHHYWZGz7u9UeK XLnqcxoN3/kMzoDWVRQ88vNwc5cvXqntJbqmafjtxp/dnhIyGOk+a9lU7IM2 8YrU5lkalic5Luv9lUTX0KuBxpHfcGykViR8noZBZ6abPvaexnqj15d+vPJA G99epYwlGkqPG29UDDmBZtFHVV8Gu2JRru6xtWUavmy/dZLzDz/+gOS4U8+c ka/w34G9qzQ0sbgfuseIG3Xqfazfn7LHXyqJUpvWaNjNYavyDFnw6Yc/X6gd Nhig2mhdSOfHtt7ZlMqzAbnkW1muPn2LFxs9O5TXaWh7sfPo5paRvPsd7PuN Wl7ii0Pumll0rlolsiOPUpz3Iz38dCihhk+zqCtTdF7tVTwa+MAjb8PuPnU2 QRnoe8KWuUbnBV7Sh1VcWvOeK4oHVaw+gSf8q17tdM7/u3P1t/VcXh9zDTfX ggm8OktzcafzeEM5pjVnJjyoxgxzr62gjtM/iJ3OBzlHuXMF2LE8vdIvdsYO JpwONRnS9X8zSNTsLOfDVCvzj26nHWHMak3Elx6HDM/nxyvqjmIHUzSxzdAF pnh5Y3xWaCi0vWz7BqazSHw5qVT22hMuMJez76THOYl94fP5pWuYYRKyqNTz E34mM+zumqLhQaqK1S4XWdQm5tRjV3+DRtC15tAJej4/TGyf1VLAWZ4dfE5c f0D04LfeS+M0+lzgs5JtrIgXHhct3D4WBNbxnZT+QRpy8OoIZv5WRzuPaYtQ 1TAwvmtQwd9GQ50ezZYHuS9wLDHTUyQ0DkwHqZxXsmg48fdXFSXjJW7vo2zK jokHw3VrfelUGnrH+YcolxjhFgtdck5yApBnsh12JdKfPzF37/dWEyysLHnn 3JkIVTFX5U+G0fCmpYvFlcw3eOXpcWepwWSYFnhifM+NhqdPCYv41Noie2TI FuGhNKCq7Y01UKHh/uuPL+g52uGpabZ8/sfpYMFpL3b6Ng2VneZ6Q6Q+YkRY 2JGcf+nAONS7o+YGDeWWZo6qp33CCSWdwvq8DFhirB78I073V3SUVaDGER3a G56ZfM2C6NKRFcu9NMwzVJhVCHTDQ5ZCh8q480B8zzjTxZpx5Cu3ffyryAdX dK70mBgUgp6WzF2NY+Po571Zpuj2D0wMZHi27U8hiC1aMrvT67D1BX5bqdYf +Lp56yk1SiEYSHhY+O8aR66DT7Q5F3zxUTLnXW7xIvgmMN9+l17nOdN8hYZk fqOqq8uzXetFICkW+YS9eAy/+OzITdsciCMfLhnx25ZA6IxyyhG1MWy+FvlF +0wYlk33ig5eLodlU6UPJw1GcVvrbNcRjzC8SLu6Sf9VOXge2pNHfjCKDnvL xdTGw9A1mjXC4lc5KC90h926OYprPWJ1TLHhuD0wTvrmZDmQikZY/x0exVXT Op+DgpHInHti54efFXBsqmZEpHUERR6n/6dwKQaVbGPnU+oqgZEv/+KJqyNI 5ql2mS9IxNclr3jshsnA6LLg/3hxCBkoow+8ZhIx1HtQxZ+pCgy+ZtvF9g5h 68GzM4ICSSh21/ea2JEqmLOVq2yuGsKwuAvTaq5J2Kc2fL7mURXcOKyr0hA0 hGo/zbU1VJPx7m3TmOTiKmAp2D6jLjuELY0bI+4PpCAG66dkv6+GD7fsJKrd B3FZ64UHE/df9JhXZ5L3rIZXQtvqJd4NokSG3xFfub/o0zPCcDusGqb0Svg/ aw2i74FFNaW4v8jk1vQ2jlQNDjICSXWnB1HDoXv/0utUzDplw1TBSsCn/jdx fDUDONb+5flNWhqC2q0AkU90blTinMs8gGJVkiDOm477y4p9/3wloPbQhbWZ iX6M8AgcG5JLx21SbaN2Pwm4MzLsw97Uj+cjrx32CEvH4Kse/aFxBBgfuLZ7 Z3g/tjXX80Xfz8DKOy9M7xMEmBwYlRmS6Ufgq+mx/JSBamNX7Z80EPBN9/DY M4F+dOTl4FmPz8AoT57kvjYC2urHrWt29uNg+pLkvs2Z6BIrOCM1SIDMewNR rRYqyo+mXNWLy0SL1BmHewsEFNe2HLM1oSI7xy7y74ZMHAmmdc8sEbD/ubKU 8gMq1gpWPXmznokLF3kzt60SEJoqvJldmoq/X/PzTwtkYcrxqxu91wggsCgD D1PxUJIb1+67WShs8TPNbZ2ASVMfr/vMVLx1akCd9DYLHYsqKUt0bselWNU6 3ocOWvqMhwOycIv4zL0WOufMrQhR+NeHK+NZRzhKs/B9aJi4IJ0PRM3fjc7s Qz39hbqA0SwkrLodJuh6GIgVltk/fXhQOO94wZ5s/JXHc0VghYDVMwd5hOz7 8LfUGcm3F+jz0aC0Wc0iAY7a1YmK+n1Y9ixy6ufHbFzSZH1qOU3APctzOs/O 9qH08fXKjRHZ2OJJumxFI+D4lPWcKlcfPjtcULhQmY3cSzak8WECDp802nhu tRe7Aict7VlzkEGWfFG4i4AtYqPGuaW96JLrYsIploMjJQv3lpoJqKTw3H4e 04vXOufnhVRz8GDjTXOJOgKkA6S9N3r0osBm9iP9vjm4sOOmBmsxAZ/FmOS3 qPXS+9mX1je4cvGpdrXHqQACKoQe7Ls33IPC+g1qNedz8VcQQ7v3NwKmohr7 rat7sD03sapGORc/Khnrv/5MAMflYka/5B60SLDnEHbPRa4XgducjAk4ryR7 Ks6yB6XF+wOP0efZo/f3GZy+SAC11IZ6YlsPyvFsixXdm4d5pmIhMoIEyLEe Wx8Y68ZF8WOhoafzkP/1ueFSHgJMLwao/aR0423WgoPZ2nnIky6Zdn6lGngq nCsbfLtRRQqY3+Xk4QN/C7WYv9WgvtXItOhYN1oYmjOFKSI+Z9l8WHRHNQw0 CLb9kqHP9ebHj97VRnQvi/J5PFcFvtlys/5HuzCY0PDXNEJkGz+Xsa+rCmqE T2r93NaFZMM7vIGfEblFfDQuJFfB+ozk1HRkJ6qRk865pyEu3plf3XevCra3 VKcRRAcyc7DPC+cj7r5wKihDvApS+f/KP/3ZgVYV8VICFYj/OFlrqLz0euUY mDuh04E3tuvs5mlGJB85ozUxSAa2GqOGntl2/Ku1eGJqAvFK12cvigUZDIp9 RG7ub8dmVqFnorOISr8Ea+bUyVCow8Dg3NeGA6O7S/IWEPOz3/bGAxn2fr0+ jHFtKOEasO67jDhM45VaOkSGzI/HVGlv2nBktupf3ipiiJrZneZNZFhYfh7D KdOGMX88GUXWEfe2vv6qQCXBpvuN984zt6GxXv/dQTp3uB9ppFxKgtZUSLpd 24qXnP9ZD9F50biCwWQ4CZT5TUW1freiveebe6J0bnRGcP6wEwn+JV6TfPGk FVWl173/Z7dHOc5j9CkJ2EOcz74WasUTGsjiQdd5SpFIlr9BAs4vbA5Gcy3I FnbKM5TuV+segdbLR0nA9SrTRx9b8KtDT//SDGK4c7szsYkEBgMKNZrOLWg4 ntTtSY9bfXHcpdnuSuCS//Pzzr0W5N3tr/x8BHELc++7FKyEyxWWDy/ta0GT qvXF6g5EjpZXv5ctKsF1akZ7e2wz9mtlt6s3Ii7v3//VUbUS+O053o2bNWMI 9ajfUQLxblndodgzlWB1teQIIdWMb3SGjujmIO6qH1CpGKoAYT7XLheiCUOH f/xn5IZ4cfjts0dKFZAyX3F+fqoRdT9tON26FTGSbb5Qk70casZNOlclGrBj i+CVbzN5uFwsm944XAZuobNdwxsbsOVch8TbzjxsdUn3n8kvA8dtsuOtpHos ujwwm5ach9J2J8t4DMugaUy1jKRRjwzjH5Y33slDXrbP9hfySsEleT56zO4f 2uvPc2/QzMXgaI9KUcUSyLLbHOJNqsW9Px/8XLmQi+JLJpcuHiiBHdqGdooe tWgvqnxvZE8u8vXK7KCOFMPXN3KxvCq1qGb1TP1PUQ6uezPr89oXw8w/t+GW Dgo6pws+Bq4cfJLJOyqZWASLNpvZ16dqkGtfUN8ftyw8ST3KuDpfAK99yqyy L1Vj+inebzdC0pD/ROKNvT9zgDvhrEUFcxluvs7liJnh+Fl8eA+FIwoUNnkq votPxu2mOXZeY2EY5SCdrcoTCa8EnIfah5PQw+BygcaBMDxWzHHLmj8CuDxC 7BrofQCncRtHx/sQlD1T/j5HKAw+Vm+Q23EzAXdPhO9wPxCI9z02PTGSCQJD acMB5mtReNlXovLYVi9kMpia3fvXCy6+L1YOn/qFFnOSGZ2lLpCp0jkvs88F taMG+qYlvOBsGIsr81gq1Fz4qXa5IxltM2dYHPoyYMvVm2wP/qTBvC+bWpJ2 Cp6bljl4wDETekr3a0ffTQf+u5bfy/pS8Ga/UH3qsSzwqPQuVc7IgNyM23Lt Q3/R+4TIl39PsiH+8Hqrr3c2eOX6Rb6g9wX5EvPq6f9ywWG33rebOfnA8/W0 wbb+LPTKsB3cdbkAgstYhbUWSiHUOp8x40IB7pTZu+uMYwlsiMryz3tcBn9a bpF+mRSgCr/UttKSEvDdeETcvqwMlFnDiwRjC5C5+1ZEOWMpUD8Yzsv4lsNv 3dqDD/kKMXy+XYzLrhSOBtyKfCJdCZt3Uwpy1wpxO9eFHQ0mZWDSE0R+6VcF C7andMrii9Fe45K6tWQFdD84/m5avhbsrLVaREvK0Adbx5UfkCFHfqvPJvta yKK8/CkzVYY859Sz5ezJcFSE75VYTi0IXO817uMrx1cbq9ViEsjAGi27OitY B2/YypQG35SjaNLa/VjGKlBc0pIZZPkHbgf9HUmHKvAwZuVHRlTBvyC9a4yt 9dBm8VPtkU4lnrCu8HlBqYaWnaqnm9kbIEek7p2OcyVK0jZ5rUxXQ+tFuF5y uwFq9nofnY6rRP0oq8tzHAT8Jj961VXYANFh4y7x85V41E/J8PADAhTMrSVZ Yxvh63PGJCsHEqY7BrLp1hPQ134mt9quGd6afCgY+EzGuMKt42lqNVAUpORU mdMM+b2PcqQiyRil+Gji8tMaGHbtv1U33wwzix3lbKVk5G/ZsMXOpAbqkvft Y37VAmxkBUZphiocEWcw3+9SAx/ZDSymH7bCsckZButXVThpbE8pSK0BqQjz Z8Ln2+EnqcTg0plqlA8hvFrWaoBsI3Ar2rgdlLaEbXkoW40Jx1K/RG2hQGbX Hsvjce0w8dnTZEKjGpNklbiFdlFAYHv4nhNHOoA3wfSio1M1Xsm/LazIQ4Fn jS07j+/pBCbeK37y/6pRbjF+O+c5Cuy4db6/4E0XOGsqls8qEej5Xdvm+wMK CJkmcxl6dUG8zVT+Ny0CGQobxz8+okCCzn9f9sR3gfgnIUqmIYFWtOgGEU0K NEUwsypRuyC318P6+ycCD9/Qe1auQ4EY5oPfRhm6wVpoOF/+K4Gk7iypv08o cF3ssrgdTzeIuF42cPIlMHHq1sDjpxTgoXHk/lLqBnd3vPc1hsD/PIFj/n/3 /ocMWQ4YdcPnkMXran8JVNOMIvXoU+Dhf+za/i7dcKtSbUdUDoFvQxUzvxpQ ICytSp8jvBsKD363ty4m0N0/auF/9/5fjih2OxZ0g209eDSQCBwaiAv9371/ VrOqzWRbN3T8qdrzt5ZA3WsTTex0nvbFbZk+E8HgCa9JjmYCXyl5JRXS94l6 GLYzia0HOD8WsK10EFh294a8KJ3vrZB/uVmkB2yMX+s87SNwsup2kRZd52CK UO09uR5Ya/5WrjxE4LsXH8Tv0f1ylluY/q7XAyNbtpwtGyNw2kO8dOMzCvAy PXSh2PaAw/x7j9xJAonDhl/N9ehxDt/0H6NfD+wN/kYRnyWw6ZJ/QpQuBdT5 o6aF0nugWKRy8NwCgW5tLbcCtCkgHb7x/O3aHtgpal2ZvkTgpT9F1qqPKeD0 pYr0ZKwHlipphrkrBLbq2xsS6hSQvVtlaLy1F4bbbpTfWCMw9UiT+LaHFDjm Qa4xPdwLxw64NCit0+1qMW1kuk+BlRX/hBfSvbBVMtS7j85vbx8cK1GigMt9 3oIHar1ArddZXaBzXu52vpt3KUB7djpX0qwXHEPesHnSOevZhgwXBQpcnSt7 xOreCwYapo2xdLuaWne6PGQpwF3UpdwW1QsBBYSCwiqBiwl3s9SuUWDaxl7r d3EvyNtT3xktE0juemlIBQr8Lqs8f7uTbldyUodzkcDkraU8Yhcp8IGzxX1m sRd4ZfawX54jUFZy09T18xQoLhw8/JW9D9Q85H7Qpuj5fOIh235RCrz2P9Px 81YffDxifqtpmEDnba3NLMco0PCtsWx7Sh/wKd+91PKPwJcr/RMCOylQt2X3 vv2fqJDb4PUp1ZPAwP/ut4nX1IBzs3hynR8VVj8f/CDgRNfDFwFeZTVgW3WN 3SqNCuD0YonvPYGPzMjRyXk1cOi4G+/3ISpwHXM9GvCEwIy3/9RuxNZAPqks cPx2P/ApbX8xJEIgJ89J2zz7GniWoFC0xDEAXa3dg/8SqvHZLcMCvlM1YFy7 Z+7r6QG4tovlxCm/asxbqbs7yFsDiUyDdczyAyAwHOQ3S68PeSRbL+PdNbBe dutCss0AuDKEew7S6wmFe4+v8xAB5NpfN4sHBsBmuu7ZvU30OnNSRZRqRgDT lrUvIcmDQBERsHIVq8JX+TIu+mbVIGt89w+L6DD8+KT/nF2MhC97Lqky7ieD iELX5qabw/Dw4eftdSwk/LP3v+LkSRKUvpUcf685DLsj3Dk5RyrxweyBxZky EtQJ5LA/+jIMflXhmwUDK/HFsUhvcVMSZL5a9tXqHYaaA59PRTNW4sHc7NqK 4kq4p3Jx+q7rCNhy12W6pZej6lGqVLlyBZh8SFftxVHoKiw7f6eyBP2Odxh7 CpVCrIBClUPNKNQ/snyv/bsEZ3cw5ecsl4AVn8fYhu5R+DR2LofXsASPGHuo BVSUAOlVob7jhjHIf1m1/8eOEtz4SPqk25MScP/BcSBbZgy6PK5fDb5VjMEN ro58P4rhr6W5Zkz6GEjGuAaX/C3EUB5xjb+0Quip5522dhyHpq1zxmlReaip +rJAbSAPJN1OqjF4jcPou7E0HYk89EvSTwGzPLD+2i1x338cYOZchE95LkaT X4gWMuTBZ5WbVLuUcYgK73L3+993B+e964O5cmHwQ4lKWfs4MH55ukdxXzYe 5qZsX72aDXftz9fPn6CBmSjbrjmjdNxe4qVk4JIOH83y+cUCaaB+f3iSOTkB Q828OyYKEkDsJE/gZBgNGmcS34pjPNpKf05Kq4qHAwcaxD7G0KDEI0KVRorD ObmNzjzNcUCe5dXjSaPBv73v4nb0xeBNRt/z/8ZjoK7lzJWzFTTISr2zoMsS iaecW51uckZCbuC5JqcqGsw3uZnJcUeg7alXmtoHIsBd5iB7AoUGd0x63Zr4 w/HNvKzk2ePhEPTB9bZ5Mw2knvKrLImF4m3rRJbtEqHQyKkU2NlPAytoC89W CcRdO+8nd94PBFae32LlQzS4f0Wm3zXrD67/U3rMk/cHYmeLeT+N0qB0oXZI sMIPr6UKZ3CT/eD40THCaJIGH3jeTpyg+uKUqHz+zgFf+O9XpUvQNA22J2Y0 fJ/ywePmKu67Z3xgXEn2Q9gsDfS+CI5HrX/DW0LbV0gM36FY23Hh8CINlluv OA5yeeLnzEuz9tye0LFPxvzPEg0oXieHrsp9Rd6KCQVe+a+gvtmNc2SZBhsD 7Tqu6LjgrLKzda2OC6z6Jzcxr9IgZHBDZ5uFI74/eyXknKUjjK5ZZCzReUlx aMdBz48Yyt9IVfT8CKwjIim5azQou64WdvGuNRr2nm81vGsNDo2hpYrrNOjw adQpsTZHj56y0iZrc5CZFpzMpPOtdoLliyov0Ng5pmCT6guYFwwUnaZzCXTU YdC6j+ysEZMbtO7DYfGHTmt0PpvANX06CWBc95qoaBLgLRfJqVY6lwowuS7l pwv/HTs+I+OniwZ6oSZf/rfPx+dJjEHG8IGNTY8lyBjDyB+Yd9L5q++X5SNM LaDCs9a2zNQCdziKZWjR9XO5hcrMBNtC0+gNKa4QW0xk/2vrSPeXemTlTP0X e7gnzMXA7WqPz69p6MvS4zb45ttNP3VXiBArS72t4YquqnPuffT474tckjnD 7wFXRP/EpvJ7oBnzhXbFBRqYqwgHN3F4wbFdm2s9Ob1w3vTiWtgMDZhbeRVv r3yHzs42hd+r31G6PGbh2AgNhk/lee0o8YelHWfYb5X7o06dGk15kAbSy0sz omkBQI3g2cOTFYBb8/aACpUG6c/j1l/VBYLE6l/h0vpA/PXOY0dVB92vXVLC qkXBwBDio+RTGowrK9mvfxE0qC1ZMIwKCIdLzAJ62mHh6JPY5pEWSQPyonjv iZvxcOiEz9d15XhcvmOs4RVM1xn5iplyOwGuN3Zf4tdMwPoXCw8l/GjAFrz5 5+zJROjMKKc9P5uIlW4nYegrDdI21TGepb8oYviqzpOLSbgdO4TrjGnAr3f/ 9ZJECpxeZlJUkU3BPWlPNvQL0aBzQ4pXjmQ6cG+/uvpWOR1zstnc223G4Qrj z9VXQ3n0Oj2zKfE44syxCl5V03F47OClyqGFsEfj1OjLYMQo9eUm/6fjcBVF ngg3IEQf0njdujMfHfU32f5WGIfBbpnjD7TyQeDBXR1jx3y8mFKprcxJ5+02 k+o3C0A16p1chlkBxjhkbKAGjoGZos5pjy1F8Njmev7otSKcuSz0jOo9ChdT bdZarpQC050vOBRXisX9uQ6tH0ehUeWg9rnXpXDl3w1rlY5SZNS9WxBrPEo/ J2WKxf1Lge+iyNc3u8owN8VCtfv2KCi8kSeJzpXCp40pOxtelqHQkXobX6ZR 4OiT6RsIKYP8wsynC8fLkct6RNbQZARKHzJyV4yVg094cu9Fjwp85lOce1xk GAReWLKxiZOgv3mHqxEfGd84N0wNcQ/DcelgmbMPSGBscjztzjUySr+PEXFj HIbWZs/Tg+9IMO3EkUnRJ6Pt6IO/ls1D4FD1SG1XJgnYxwd3ZaaQUeTx09AF uyEQ3vteRFqMDPpmFR6br1bhljtuvWlVg7AlTceimKMKBEya3/HeqUbhCcF/ DIr/+357TM3lYzXMD7wQaNlRg8cn7ZSNxAfg+SH9qSWfavjSfTLy9b4a/O27 R7SYZwCuq6cHdERXwxx51+NSgRrslDuqwT/QD2H7m8irNdWgN+h0J/5yDXp4 VfBnWfaD4yXhHV6cBLR89v703KQGmaIbjug97occZvd/nscIsG7euThhXYO7 raevL17th8+tsjwbzxNwbvfJt+c/1+CU3QH5jh39cDKom+WsEgFyZ6o/zf+p wbJRrnNh/lS4veqgSzgS4HvgYrFcWQ0qHBzbEGJHhe3xthPyXgRsj7rz+G9N Df7n/CnPVY8Ka37uO8T9CXi0MWhjc3MN6nlUmRwVooLJ0P1TD5IIMNezNFIc qcHUjKgi8+w+sFr4dLygloAUcXv99W0UTFM6ytz+pw96Qopl4poJ2BUyVC65 h4Lbe3bsO/uxDwwNrgds6SLAcVu4pgg3BRkfeu+Kk+2DfvG4HIYRAhhYbihe PUZBg8Pp0ZRTffCCsvFWEI0AV69f13QEKXgv+r9P/bv6QPjbm7T4aQIubD/h I3GGgppXmNR7/vVCHmfyKMsSAT2Sh5x3S1GwgGzLUZneC7rsnL66KwRYCDAV sVyhYBLL3g8hv3rhkTRl6PAaAd89DnqV3qAg50j01CvrXmBQqG2+u04At7d3 sZg8BcWZTtcIavWCyEC8zjCdk97JnXx4h4KvHsSYtl/phQslIXYTdG4GM97n lSiYGiwqaXe0F9xSvp3XpfMNfJvrKlQpeH6VastFnwuUfHfYX6HbjWdKz2R7 SMF4lYXQkOEeCHwqZfSVrnOf0Lb97BoU7J4O2XS4qge+bWFguEn3y8iqpLLy MQWZV4WWfBJ6QET0n4zhPAEm3G/9xXUoeGG+l7ru2QNT0CW9Ro+PR4/Z84dP KNirsO+IplkP3JmYoB+HBEx2tm04+5Qe55qTe5JVe6Brn8Xf5/T46+N3mYJn FMwF050rF3pAS9NT6XI/AW9k76ww6lNQolzJ9NL+HvDKm19zpL9HAeLjzjU6 N3ARCjJdpc+J51koEi0EzIWa/hdrQEFtPu2ZoM5uUJXp7Fero+ehk4gPywsK Du/9r7KMPj+uK7BpDZAIyNvpXkifH3HjFwVbakg3tIotafYUEZDPej5+hr5P tUeq6bJDN+j7OmyVzyYA9n0794bOv5MbN2x73g0Pzt+yP5RMQFJ+37lYuh6B ySUnVrluEDpnuaIXSUAtx2k3n+cU1Dl9ZnLHrm54pr7fV8KbgJGExhwnPQrG pAR+3TDRBdmlHLwVTgQs1F9X+apLQY21At8JShe4KXtol1kRcFGgrlhWm4I8 P88+w+9doKLevZVVhwBx2q3kpkcUlMw25TjJ2wXTjZ12nQJ0ndbUZO3/KKiU OBi6Z70TxKUl3olzEfDR4tUm7lsU7FOamCIVdcLwgzTNzVsJ2BtavN/xGgVz 2NZKRP/rhDnRw3yM/dVQ/DxT3k6SgixqOg/aH3eA6kmHlu+/q0Ez0dxgkZ+C QQ/+65R3aAN22cQYlbEqkMq/LNpNrUEzuVLfs/VNcI9VY3G/MBn8+Rd9L56t QW6NkIqfv5qgcZ4Th3eSIZxjLsz9aA0GX+/etVurCVgCImYujZFAljdfxHdv DT5V+KsjMtIIBTcaZl9Gk0Cznt/ZjT43v469PSC3sRHcTJnS9hwgAXMCUSpm RaDKwnNZO7F6sA3tuTY5WQFezaflrxpU4yERA5aTcRToz1V9lHelDKKbZIXZ m0ho9ULB770FBXbkXw/U2lsGzlcUqg7GkfCXVF7d2HUKKNWKpkSNlML1rULH 4uxIONndtEu4swZyh4sHzniXgvR8Yer1EyT00X956gxbDWyXfX1frKMEMj0U Na1MKvHxzkYxlU1V0PGmbsTtSTEkcR+gt6rlOK69kdmNfm4VtZ4VShIvgKzD p1p2O5Tg85YPxIemdAhn/GFBvfkX/G7aV2FAFtoy14wdEEiH9NnmY/dLU8A7 In3g3YYs9IploeS9SYNvxNP1V9dSwPJdUs4ZnUw8+IIkscyVChMSRoa/IBnG HAb44g5mIFw3rvLST4b9vZqHDp1JBL5d+zjaPVPxrhLzJ/930XD55q74abZI kBCzOnnuRiJOGM1t5nnvBhS+A7ZeMa5w+VNNW8yaN272bomm7YrER6vKoZbZ 0chc7cbPupwCB40XIlfrI3HdczVD8W8M3qoaCHDJ/wtT8Ux6G35HIYfnX/7c 2Fi8upPP1csxFUb29sUO8cWgCSXsp4QfvU8KYOSRZU+HuRKhG5d2xSOnEXkl WTIJnyzkpPOfzwK3iszf70uSUG/rlX3+W1LRejJY8NStfMh/c0HByCYNP/N1 Cmio5+DPfsG2i5RS2JLphbaFaVhS/m+DfHUOTq8Kz727WAaR5hrRXozpOOQa 868UcvHnF7dMalgZmOsKeWU6pKPUyy5t7UN5+EWE+Yi8VTkczls4dvB9Bm6L em4w2YH4wwILRI5XAlvkg2GUyMLd24QX1wILcEBWUTPxcxVILH9MslrJwV0K br94q4uRYWHzgdUztXB4V8zVxQO5aNUpX6u9WowP0hkYLurXQrDlceNvV3Kx 4QSrqKpgCR7WqLD+FFALi/wdZccccvFUxtNfCp9LUNPVp/DMjjrQZJI6Xrue i4ayjR3qUIpRtlFu69Q6YOYa2LyJnIdH31Qom/8uQ8NcneilX/Uw6xjDiN15 uFnwxLRBSRmm5me7plDqIffN9c/hc3nY1tl3eAutDB+WhFe83toAzir9K8u8 iId4k4pPXi5HLccTOGfeADNxvDvadBBTwiKZebrLscxAwfnovUaI3z44o1yK uJ71Wns/SyWKCBxtSGVqhrI7j2cuE4g3x02W+M9WYr5kXq+9dDNsqH2/ZtCI +LY4c7BUtRJZ/0rN3DNvhqV7DyLMqYirH2Uiu/0qUSH5kWNfTzOYXWB1VBuh 6+mreXUfK/FABd+lKO4WaJB9RLOeQLxWI2f+rLsS33CyaRjcbYHUlzOLrguI TQqnpZSPkPDh2T+yXbkt8Em3edRiGVGtiJHh9HUSWr47e/PbTAtw87/eGbuK aKbmuzlUj4T+rgtcV0+1QtSvxs9H1hE/Pn5Tm+BAwp43A/dHtVtBS/yGVx+d 6/5JG1UMI+HV2EV/d99WcBdhVBuj82sC662uxSRU6z/uIEi0wvt/entk/ncf x+t4QqOXhEeWjCIKN7cBb2PJTDvd7nFLNq2iDWQkxZP+3LvUBiRHmwvldJ3h lRt34AEyPv92cXOzSRvMbx3ay7CIODVwLfuOFBkrWwrNHkS1wcjn/wjHWcTq SO1Vs4f0vvO4vhPR2QZW91ryHUYRh/Ru2lp5k9GyQmEtWKEdThsKWTMMIH6Z ODT4KIGMajfSfdY+tkOTVgkH0YVYbLZ6mqgko1GFn+TdzHbgfe84Y1iLqBwY VmLIUIXsLq47/x3tAEPmHyTpSsQVj/fXffZV4WUhZuPN6h3wTDohUL8AMY/D cNc9sSr0uz0ge8qrA3YcOmVJxCGK9aTPvntahc+YTqHiWgfcPbYhj8EG0WZ/ VpVTSRVq7v8UvbmhE17EL1//8RoxLvZQGndHFV5nK9QU3doFIrELhebaiA6f ohqOzlXh7c9XJe4f7YL8WDUn1suIOwt+scUcqUbv/xo6nB53wQTHMO/sfB6O /jW8pWlVjaM/xDrM67pAs7qffel6HvKzR/zjYCGQA6dMFDK7YV9FclCvYB7G JjJP8h4g0PtnYi+lvhuyDOaSRNny0CE+NuKXMIG3CaEQualukD/1KHdjRy4e dImpmZAnUM/LYzfHqR7YWBduHmuYi+LvDDmYbQnM/aafucWvB0JlaF8LjXKw UMHDL7OKQC9SwaVjVr1gsTz+4d9/OaiUV3xh4h+BS/buLEU+vbDp5o9d+qdy 8GZ+rYxvK4HaQpXX7yT3wp374tmHqNmorpBx4NIA/Vz8XUQVG+oFB3BNulWQ jeZ+hoXiowTqNJv9c2Psgw11IkvLftn4dOrv89AJAj2PW5yuP9AHIiJjUkxK 2UhWkzKtmSew4lv7RWHVPvBR/sn1RCgbXd000t2XCGQpkjUG4z4YDAplVGLK xuDftXm4QuARZb6D0q590Ljydrq+Kwsj99TYaK0RKCedv/1kRB9YPtjdOJ6Z hdUy26eN1wmcsEmb31jYB/peM0EBXlnooZzNNkPnA6kJEeVtdH771M1hgyxk NfnZ2U/nWkl2nZbzfWBHOZ1VcyULtXK7VOXo/NvOk7rce6gQ2TY6rrwvC/XU br/no9ttKq6eCxGkgoBTZPO7yUwMfHRIhd6H47WSPM19N6nwmovRSLo8E6P+ 7hviovt1zsHk8XttKjzk9n6ZZpaJaj93vGqYJvCudIvURh8qXE15XGAol4mj bbusOmgEGq0rhR1NpMKK1PA38oFM/KEvZ3V/hMCVc3yZZyupIOl01cujIgMH lrZF+XUR+Gsb+/KuNSoYHyIiFvwzsE4jYIt+C4GV8c//tu7tB6ODH9j3mGQg K80lPbyOwNPzWqe9RPvBXKilUos7A4dh/eubYnp/s/zFOftJP3DbXv8S+iQd z8nv2mAYQOAX2ZevFiv6wftNu2SPRDoe0Zog2X0ncO9rbn+t7n4gv2mdH9xF P9fc2F/zfCFwdzipN3mhH7hG3kzfz0hDJU6BwS4zAlejNm7kEhiA16u8ioNb 0zD9mNL7mCsERvgGOXV+HICQZ4MGYsF/8VCW51uzsmqs32ZieOn8IKwEfTUx yE9Cu3s/jQPfVqH6lwrOFZch8BS4se2fZxL65zFVFGlWIRGj+/t1wBAUFS2+ uKCbhE1NLr/Nrlch72wHUZkyRO/bTF3nNyehV+ItaVvWKmxYTo0QaRuCu9qP l97JJ+LBG0wpcqFkZEnOy0g7NQwXzHhNQxPi0PMXkbU9lYQJ389L+pYMw0R2 mEXxegReOIprHo4VeFn8/2o482gqv++PlwaERCIpJEWUQoMMbVPJUElESCIp mhSiInNIyBAiJEOSzPOwzTOPiz4yxuWauYYMIb63tX6/v/Zarz/O3ud99rPO PmedZxcmxAyMA92QfviTqARcjuta2HqjFqfNXq5MzI/DIp+W0itIQNRzymaQ rMVjqn4z7AwToHdZ1K7UKR53Pjjgn95Tg+OaETwzohPA0enpeY4hDjva2aSz jtbgRJlm54zVBMjIGofuEYtBFtPke4fzqnD/s30Z/QsTsD/QtfV2Wxge8mUZ qXUtxy3RNhZPBiahnNxp9Ph0GCZ0NX+5erEcVWWX31KmJ2F4LPr5YngoZsvo LrpxlmNe6mj6sbVJqGRXqPC5F4LDNQ2i5IQy7Hk673Zh9xQApwPPTa5gNGIY qo+pKsWqSfFlT7UpUGD5T66p1B+b34hJKIwhmqlGj4kmTEEwI6Gwu88FJU3/ UxmRzkWrK/YHleWpgBtkSpsPuOChtsVyzYoczMp+rdapQoVAp3hZBktnPD9C VtO/mIPZq9K/VS9TwaKvwZdhsxNWTOCwrlE26uygD8u6QYU3WcYBKaXPUVRv KUP8ZSaK3Uh/V2hHhZb2bYwXAq3w88v2+K2Jqfgoe5efQwIVGi8GpPz2fYxv tuoE3Tqairpe/N5x36iwz372nY7gIwyvry5eNElBJ4Vgmfg0KpjHMXZ+NLLE z4fc309YJSPBvjtMvIAKiRyzn/S8b6Ma32ymlXcikpU+ST5qoIJssv/zJv+r GLs5T00n8RPuEEoOdRymwmUlCxmmB1fwVE5sZV1SNL7uOtB7dIwKcb5lvF6k S+itLai2YhCFVDZWw6x/986nuQ/VOqigCsP+H9dMw/FnTc+ExgwV7G11W2sl lHGyRfCV4Z0w3GWiv2A4RwVbUqiCiJACtmWkJq/eDcHEcoeD5+apUH+CVF8Q cRqrfz28c+RBIDYr3Kd8XKJC5XEfpmQ6SSxKF+uqf/gO3wveeLR9mQo5Cg5D Q8RRjBWO1VLv88VrlvJsuitUoJz3tNINO4QnlcUCXMneePuYTc2zVSq4Kn8T 7b3Fjxltt801Bz0wUvlbwOO/VIgdfMx78SonHt5nVFVKccHXZ0ceKKxR4YOr 8RPtcibUDZdRJRodsMDaQW+Mxr0PejdKF/4tFvb8OvOZsMHwrexXLNepcMpi H6dcz3Dx4w1+Y4oxlqjX+/FaBY3X/d1a7TZbUvzuYyZ7lKouSmdY3J2ncRWT b0pmlIDicRsLG1VTBTjJxOHx7z5alq9uZaSjo3jXw2GKybXbUC3zOuUXjVeF LNu2Rv4ufhxqzUk6/AQOWmyihNL4ZkoGH5/JFqT7wS8sLPwckr25hURovGLE ifr6DDuejd5iyLTTGT7f/GYbSIs/J3DOw+7zAXzrX3Vx3w4vIAvfV5qi6bNp OkQ1wUwEf+zJZd69/S1Q3y2W/6LpKR1+98fdneLYPV2rT7j7w5N3hXqJNP0F V/STI7SksbOobIOnSxDcfSFZW7FIhfsjGsVJJqq4/sFJh/wsAlzQmm/fNBXS drN8PbeqgQqlOfnSFpEgT/Xx6pqkgs+nj+Iq9zTxlb/GYVaxaOgflEh9Ok4F T4alxFsftDEyjzHyBF8MXD9iel58iAorlnzfFe4YoXL/pS+8a3GwXml2lNxB y//+53fC6izxF5XkeDwpGeIVc4L486lQw7H/WtLYK7TqFZDS0s4GLaZwXw0j KihMFlEUPjvhXC41NbU2G956fTK+okvLQ9PVj0HazsjPmRWfczYHtrfpcIlr UmEuXsxLv9QFWa3bNucK5cJT5r2WzgpUCO6dF+wtc0e+Vs+8K/N5YL2ucd5c gApH5p8LHrrlg6w8RafMHYogMs7FS6dnCjT+vjTzmg7GmJK+06l7ysAt/PkP KakpUPmUBVvc3qPP9fXUj1pl4LRXjb5EdAou+sjms3CH4I7c48ZS3mUwefGX Mi/fFCycZbyQphKK6wKaNlZLZdBiZaZnvGUKmJi4c3KzP2BNFPlGe3M5mO4O F14jJsEstYJfvT0KnzGrXY98XAlc+wY6C25MgoXehuophzhk+iIuQTKrgV7Q fXXUcAKE+R1ENLri0Nv65qeP72pA6kRL2A+1CRiMfSLBeSYeW1JKWVYKa2DE jcpx/cwEJFkr3RuYj0dZBYNmM/Za0A4yb5rdNQFsqZa99NZfcKIqhmkopxYc apNqbRrHYe+aUkbgqyTMdmx9Nj1RB7Lpl6t6T47DbPnqQkNPKkZ7JLoELTfA 3kBFduvJUeiLDRJMYE7D41hPurCrEc5e+jZysWMUBras33SSSUOt5ZFE72ON IL7pVCFz5SjM2zbNHA6j7acZ9tezTBpBcyVs9UzEKIRrn+0S005H+43z7eHl jfBX0cNAV3UUrjG3sSVUZ+BpY2dFw0dNsD+5Z600cgRY3/ML589n4CAXEVLo 0gTVe75vW/ccgR430dCMA5l4J6GkPTa4CbKYxdIPPB2BCemvUldfZaLO7RIq Uz6NX3d0EVQZgfgte89wn8zCk54WpCN/m2Dl6+kFhqlheNRkx24Ylo1nveS2 utwnoPbeHVG308Nw45zkPc7KbJToaDmtaEfA/T90BZf5hyG4ezXy43Q28jTn vfVwJaCjxT5jC+MwdM6XSTGr5GCX6ebFD6EEtFQwh4p3DgGr5Ryz3EwOuieL GU0VEqDdxXHB98UQ9G90OS/Ok4vvqTFK9RUE9Lukn0w3HQLL36PPKcq5eKFf jF+0gYCKXdp0Zeq0OuoWy83XIbkY4XlL42kHATImbxfieIaA2dGnx/VMHu50 JBUJTBHA8UvOzS2fApP7C3LUbuVhwvSxcI5ZAiIY9GKYYyiwvLL+qcIzD6sl 7xV5zhOQSZnjdfGiwF0dlXVqSh4+dqdwvloi4HhepAjlMQWYjsVJ1P+XhwUS q+GLywScuku/cFKPAiphhlHaf/OwrHJOcXaVgGrWtTRroMAa0/AX7wP52CVK t/nRGi3OD0/9Ph2iQE7xSrqlaj5qtGl0PPh3P6xtj4UsFNi+MWT/wsN8VKke KZ6icW8X+adVvwdBPnDjdeGgfLyT1Zk+Q+NXfnN3F3cNwp619qINtLr94Tpn li2NJ4bZiMbT6vbcMbYsh558dMh2rnag+f3GSQ54mTgIQs8ZYz9uLMDQoKnJ rX8JuKTLrKzoPwgfIlymjQ4WYNQkvxj7CgH5Avh40XYQ1M3TN5RcoPHQer8o mg5D38lqETcG4duW+/QGvgV4YCddnzxNT96XyUSWyCBcPWb7515GAXoxRp5c mSQgSqiUg551EH79VFHY8LMAv+7vXdw/SkD4cze+83MDwG3yOz+frxDLV0J4 s3oIGGsUeu1bMADPvzyniCX++9+ffcc8ElDOKPH9svoAbLW5xypfV4hM60qt u3II0OJnOMZ5bAAKrnLY14wVopNQ64Z/7/mDs8+9aGAfgN6WbI7TIkXo4amm wRhGwN0UnReMnWS4fW3T7NWYIrwRJsX29y4BFy+InG8wJ4Nt9uD74eIilHQ0 dfY1oHGD7Iu71cnQUc2eQu4swvMy83R+FwkQUn8hdk2MDIyi8Q/o2YpxSCIT eo8RQBHLcfvyux8+WrSOttkUY++5J7pFk02Q5N9c892pH2YX1+Na+RB/1rme rbjUBEQao+DL932QOv0f55I44lEWwSoLqSZo5eLasMW+DzbkFAaZKSE6fq+N dt/fBC9WNg+6GvTB7Qg1E47biGo61CXu2UY4XddNus3fB2sxGy19IxAnhLzN an0awV/AXd/m6i+4xi5kWPMF0ZTLPZb3SSP43la+ucD1C2odnMdkMxF59Rnq x681gh0rZ6p5dy+UcW1pa6pBHN09U03lbYRZ0TmdfWa90N4SvddsBFHabqDz 25cGUH/W+37Apgei1tqElaYQv+6zjiS8GmCHUfuv+TM9tPNpzvrFOcR1W8UX Dy0aYNUyM3v1bzcs+9HNuS4iJoxco/NWa4Ba0enPS6XdoDx9VndgGfFJ1+5q AZEG4B/ZfnzMvRsu19S7WvxF1NipcfMMYwNoX/o4QFLthvCNoy1c64jbV1hM m4brYQvbMZVUlm4wCiz59O/d9XVKy31yZT188r4x4E7qApPJR5q9NL6qaMpt +7ke7hpXyV8N7gLda2yKv9doujX+pPd2qgcVwW0ju693wWu5yu7DNL88X8M7 eG/UQ8Fi3a32vV1Q9KrAy4EWZ/KLmxInz9TDN7GkU359nSAQJZs1uYC4YMWf 28hRDwYz6lT52E6YmkxvsptF9LnAv3OUWgcGv9Q5Js07oSEp48K+ScRAuYHJ 13V1IJJ+TixItBOKzLOlvv1CFHL/267lWAdFcramP9I7oL1a9XVBO+L8vpsu p/TqQPLaDmPvZx3Q6/o2c6gJUb/JnYgVr4PxlEs+srIdICvpx/QmH/G8rO0m wYFaUJwTtftQ+RNKnXPXaWUn5iXK7C2RrwXHn115aj/bIX+n2PWM2WJ0zBEb YB2uBrUjAwGTf3/A7r1bgnu7inHwUV/UbF41ZMjcoKaV/YCb97Z931FRjPIn ryda+lSDYK6r+4vXP+AGkZbi9q4YDbcXEQeOV0P3BJcRP9sPcHYKPPeGvxj3 T7xlKn9SBWIm/JnFgm3AqpsTyL29CI8yd3LYUSrAy+QKT9HFFmA6Yt8rMlyI urMCwQLpFVCW5ZixY2cLpFTqh4pjIdpmWJAfO1aAUEGSldVPEnCMJ7rz3C/E mgOFjqu7KoA5S1TH0JQEYYoyRip5BWhBkK+MypdDtJbnLYp9M8xbcKnK8edj Q/iV1TGPUhC41gwm2o3AJ/Jb8ylPNnrp9t6szyyAI+PDX1hOVYFQ8qO52P+y kPF3Q8aaYgEEnGh+ZdJYCfeKB7d1+Wdh39ZGvZrmfOgo1P04crsStO7Hw1W6 LFztjnA5OJ4H6V4HY94EVYBerwSnZUcGlljeFZTiyYVTRX4PFRbKYCLRtXFN Nw1P5a0PqlpkwrIE5QqvJ0LWua4ODdEElJNwY3zxKRGwd2Oh/t4MODHTpvZW /TXMl9xKvVDjgVUc9jEtPX6IL+40vax8A/KltQLbmd/gMS8Ryci+AMz6cXRE UtAfEkmUS3qX/ZD8TCfid1Mw3q86u98nMwhior/5cY4FIrk08M/DT+HIv86B odsiIdFf+4uhRQS6zOvpOfh+xrpQ5S6fii+QWisdw/MlHv3isu/XZX7HhgGF KtJENojkfD3nFp6BJ3rV2+kb85ClotxCODIHDk1pWpnvy8SkZIrYkGU+vvnp XvZSMxf2tO+nTkVm4gGTLQnljAXY5RVNz52bB7EU+8T+6Cx00jjr+kyxEO8L hZ9dCigE3gcvZ3dF56DsHpsrsQnFyBDFpEwnUwpjYonfi70K8AGXy62d0aV4 6ODx/L9PqiE5LzFNy6EU1T6uO7mTKnFOwqlCqLUa+rzeHQhLL8W0goqhBqYq XNBRN0yUrAE2WUOVsOFSfOZ7RibkXBWmNXkxp87XwE59V4NyzTJkFjS0y8ip wpOHjtJlOdSBW+fq7wS+crS/c3xt4/tq/Nk60j4Y0Qg3zUVSNqRUII9b0Fiu dC2GTVkdHmJogZOu29cifauRe/szw06BBoy/wnctRq4FDt0+0vEhrxr1TzeO rWk04LMmu/Mvn7RAwk97a6nBagwZD5qMtG3ANx1dbq+7W0C7h6pyTaoGde+L pWvWNKBgj3PHZForHDa0kNnWXYMKJ8WEn5k1YuFK+OY3t37AkMYWzuoddUge tdDJ8GjC3785D1k3/4QjTSOy7qI0f6uyywISzcipnDfVvrUDNgYLjC5daMC6 H/oJAWeb0SQ1x+OCXAdEyi2FbLjTgAd9F0xL1Zrxl8JDrktfO0CYD5XqIhtQ +qAPk6FJM7pzPlHrduuEGqkPe2SYGnFU2knSw7cZVd0DtP1luuHd+EtjaaIR KTVhpamdzShXcZhwedINh1lY2osojVisGBnAPtCMhXJ+RfaJ3TDz5JIwsdKI G80ZsyTHm/HTaHTbE64eOJoXwvThUBN+AA/2zOVmrHTnXHCZ64GRs1u31dg2 4fIZ1svfOUnowBsTHpP0C3g0v59S3khgP8PhmEBFEipEySyGk3/BaoRFpPMO AkWyIlao50nYcV81qIa5D3zJx+5f4CNwj8OCPbM6CcczVR9lK/VBl/7MdSsZ Aqu85o96aJHwvZar0IO0Psh75FAm8pDAW99jfB2NSEh/3tyT1NAHstu1ApLt CXwHnYd6jUlY4RcVKz5C27eN40/luhHonf3ddc2EhPGsygsDvP2wh2tZ53YY gYoh4RTnOyQkpYtHnZTuh8woZpctnwm83KrBTTEnocENuzFnnX7Q0K4aPPqN wHLhhoeM90gYQy57V/O4H6p9jhv0ZBJIL1G9NkXj1n8zLzH69IPNRYGKnUUE htyxIwdbkDBMdqJTKaEfEirlN5AqCFTvjZdco/EP1A3Tz8r64W31+gxHA4E+ O+OZhS1JKL7ztUpcTz+0SsU4/GohUCFWx3Mnjd8+yxpev0Qb56W01/EOAr+w 7S8roY3T6CcUO7aTDJn/fV2l/0Vgn+P6j2M07vXBcC/dMTIcqY4rNRsk8N7D p33GtDjzr0oXsamRYezm93LVUQINHb3Xrt4loeKGi/LcZmSoCNg9nTFJoMeB nouMNB3oBARu7nYiQ54d46nYGQJ1TH0nnMxImGWm2ro9nAxaAc889szT9PdR /5NvSsIpf51rq1lk2LtRtI1zicD5KfXAnFsk9H/U86mvmQxb84t2hi0TGPlO s9P2JgknSl8+Khgnw1JquMynVZo+I5+nlgxp8xL5YeezdQCM/9IpH1sj8Olr 7hk5fRIe3Of37Or+Aahc3CSktE6gXDvLpvO6JJSs0z+4Q3YAvDSd+7r/9Vep NDvPpk3CB0NB57ZeH4D819wv/tB465/HklqaJGxPyAmi2gzANULkvAeNtzVx kNU1SLjOyZtXHzAAVsP5Sq40vxuU816tqZDwm47k9Q8pAxDRqmM3TYvzXbKz tIUSCRM3W/MYNAzAvSqW8UbavFJPaij4nSUhU+KjeJbRAUjc+CKOh6YDl2Fa +YMzJHTPl007LzAI1PdPBRloeqrzXfqrK0Zbr5Ohjg1nB4F+cwNDwgSB2+p0 uG4Lk/DEm8oURYNBaNN8ap09TOCl/SFPRARo+n8s7FgLHIR9juuyQt0Elry/ nD/EQcI+pcch3zZToKibP+pqMYHh/io6Bxab8RX/pqYkfgrcXjHMtc8icO5h mTQ9tRlr+qo/h8pSwNG5T5ONlue/ojexJgw144VXexvPWFPgYqdBUnAIbd31 N93l+dGMfVd3O/uTKdAkeFppyZzAxTUW+Z7vzdgboKypVDAEW+snWHSnmnDH Bp2xZZ1m9PmsHzb3YwjUHxwy+t5N4+MWQ3mqzZh8oZ3TZ3oIZKZmV57XNaHB Jfs0KblmjAiyMbA+OAxGmv33bGKbMPHljgrbA83IEkwQu32HaecM9g9LV5rw RqxVxxLtOyBH1hnf1R8BFtZ4D5/ARhTmdRn/OdKEhhuHD4x2j0LA6zcqDJn1 2LDJbv94cgMuBOfetp0dBbmUsxIHguqR3Z10R8CrAT/za4sN0Y/BA0wVb3la j8sdopnLpg2otxh4977EGGw1VbeYEq/H/YxXI/13N+CV2BIys8cYDL88uXY0 sQ6z14xMbznU4yH3ZrAQGgcVNs12ebdaZLJM+cN5qg5l66V8t12ZAMYb26Fk ugrrUifo18yr0WbOkNfeZALm4wwfWVVU4eJ5A107qWqczEiWK386AaXmEp61 IVUo5usyHs9QjbW+FJPR4AlQc6o5ZiVXhVuOvH9+4EsVWqcYTD77OQG+j+hq 0twqcXXgmPu3gUo0WJB0s742CWtjCjMhmypQ6qjS74CLFWg2dbIv68QUjEcF lZ+uKcG6JJnrtjMl6LRp9mCq4hRYPIlwHXxTgnvYhfxUM0pwzjtFzlZzCo6g UZD85RJslWRq67AuQVGDnXvNLKeArbvUn/MrohY7hceSdm4xbi1Sr4mcgvLN R4y58oswyWgL3cRMEW7vuBT9cmUKnqfZJrljHjpwj4uvkfNw1nvXGd6tVFCK tasdlMvD/9jseaRv5iF1oqLHcwcVPt9JItUU5OIT58SOHV25mPx7NiNXkAri t3J8Dmfn4M/a4iMhpByMji2OdFSjQtveuajAuCy8fdmnXSE/C9Ni9Aw93lCB 0To2Vd4qDbm5Kv8bfpGGD95vjP0SQIUc5VSPvaOp6Fp4zrdxNhVvWaouvguj At234aE441S0lFykPrqXilpCEYZ58VR4s/WofsjH75jp9dN7KOk7GodyBOoV U8H6Itu43uWvSFf3vNlY7yu+fdcWk1pOBe626w9eSCeiacDfU48UE3HpymEm ooYKQ2SbMYGDX7AmsNqE+8gXFOB9H27UQoWnafR/dv2JQ9PizOtaa3EY9aJg oIpMBYYTEZ38EZ+Qr+SmSHPkJxQejOm8PESFyOkHweuu0djgYrn7lEc0bhwN S4ocpen2cMaXpz0S3afv02d1RqJ535Fvb6lUOGhh3R/+9AMWigW/6LP5gHGv uKqOztLibNI9tUshFHf9J+HZphiKrjYpvoG/qfCx7fj53dvfow/7iQd3Wd9j SmsSY+oSFcQOUxyG4t/h0bCerwEJ71Clgj75zr/+2HnPZH0pvnjujdiC9RDN 3tzQT16hgl/4kbbcl97o66p+UM3BG53UhKIF/1LhVYCJkslOD1T4Y2GUTbMJ IwMU0TVanFsYWaITnfG6JytvPc2OPHj6fYnGSzyl7rO+eo6qfeblkjT7ZXpk wXudCj9P2Xvf+GqFQ+MSl+/SbJzt0SIyjRfr72JM+noLXbZfz0+m2f/vN250 /O2/PuPFHv9n/wdSxRgf "]]}}, AutomaticImageSize->True, Axes->True, DisplayFunction->Identity, FaceGridsStyle->Automatic, ImageSize->{256.9463923152913, 440.5579603717025}, Method->{}, PlotRange->{All, All, All}, PlotRangePadding->{Automatic, Automatic, Automatic}, Ticks->{Automatic, Automatic, Automatic}, ViewPoint->{2.850000727794993, -1.8192858719532987`, 0.1330216812371001}, ViewVertical->{0.9380120195824551, -0.8505523701338612, 0.8224294158141447}]], "Output", CellChangeTimes->{3.688756264878182*^9, 3.6887563096055813`*^9, 3.6887563833107233`*^9, 3.6887567358373957`*^9}] }, Open ]], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", "\[IndentingNewLine]"}]], "Input", CellChangeTimes->{{3.6887578661990657`*^9, 3.6887578663494024`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"c", "[", "t_", "]"}], ":=", RowBox[{"{", RowBox[{ RowBox[{"Cos", "[", "t", "]"}], ",", RowBox[{"Sin", "[", "t", "]"}], ",", RowBox[{ FractionBox["1", "2"], "t"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"a", "=", "0"}], ";", RowBox[{"b", "=", RowBox[{"4", "\[Pi]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"a1", "=", "0"}], ";", RowBox[{"b1", "=", SuperscriptBox[ RowBox[{"(", RowBox[{"2", "\[Pi]"}], ")"}], FractionBox["1", "2"]]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"np", "=", "100"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"xi", "=", RowBox[{"-", "2"}]}], ";", RowBox[{"xf", "=", "2"}], ";", RowBox[{"yi", "=", RowBox[{"-", "2"}]}], ";", RowBox[{"yf", "=", "2"}], ";", RowBox[{"zi", "=", "0"}], ";", RowBox[{"zf", "=", "7"}], ";"}], "\[IndentingNewLine]", RowBox[{"Animate", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{"GraphicsGrid", "[", RowBox[{"{", RowBox[{"{", "\[IndentingNewLine]", RowBox[{"tracpp3D", "[", RowBox[{"c", ",", "a", ",", "b", ",", "t", ",", RowBox[{"Axes", "\[Rule]", "True"}], ",", RowBox[{"AxesOrigin", "\[Rule]", RowBox[{"{", RowBox[{"0", ",", "0", ",", "0"}], "}"}]}], ",", RowBox[{"Boxed", "\[Rule]", "False"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"xi", ",", "xf"}], "}"}], ",", RowBox[{"{", RowBox[{"yi", ",", "yf"}], "}"}], ",", RowBox[{"{", RowBox[{"zi", ",", "zf"}], "}"}]}], "}"}]}], ",", RowBox[{"ImageSize", "->", "Full"}], ",", RowBox[{"ViewPoint", "->", RowBox[{"{", RowBox[{"1.3", ",", " ", RowBox[{"-", "2.4"}], ",", " ", "2."}], "}"}]}]}], "]"}], "\[IndentingNewLine]", "}"}], "}"}], "]"}], "\[IndentingNewLine]", ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"t", ",", "a", ",", "b", ",", FractionBox[ RowBox[{"b", "-", "a"}], "np"]}], "}"}]}], "]"}]}], "Input", CellChangeTimes->CompressedData[" 1:eJwd0U0oA3AYBvBZPopyMITLNg5rtZSssDkM87GyHZaiphibUpM2M6GmIWYK oeYjayzkMytspjFclNwmU7PVDsQc0KJkZs//8PS7PL3v4WG2dUuVVAqFwogH LnnVAsdguCLLPd4Aj6N1TfCz8kUG77z8Vrhd805csOmUkLrc2El6cm035BXJ NbC9N9YPzSNOPUytTBiFpW9RYraxYA5mGuhmOCsps8Aupn8NXrA8u6SfzrJD BW/PAX9WE5zkz73uAv5JYpfQkj5wCw1sgRfumPQhuJ+2EoZWc/kH9C2qv6Bu Ky8KrxcmqU7c89hT4RTdS1TxrTT4oXlmQMFNgA0LnZxiOG+klUDh5hAPvp7O V8OYklsPO5JFLZDbm6OA0l1uxklcU76cWOu2R+BY4esX5Nhsv/AwwohCvfYo yRX3aeaK6DoIqs6wi07YA0UToT6oesgegInDuX44XZXyDHnSsm+4HhATzx8j 4vO49qCvAb5zhTLo0W40w3+QYhn2 "]], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]"}]], "Input", CellChangeTimes->{{3.6887586823687935`*^9, 3.6887586830235653`*^9}}], Cell[BoxData[ RowBox[{"g1", "=", RowBox[{"ContourPlot3D", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ SuperscriptBox["x", "2"], "+", SuperscriptBox["y", "2"]}], "\[Equal]", RowBox[{"4", "y"}]}], ",", RowBox[{ RowBox[{ SuperscriptBox["x", "2"], "+", SuperscriptBox["y", "2"], "+", SuperscriptBox["z", "2"]}], "\[Equal]", "16"}]}], "}"}], ",", RowBox[{"{", RowBox[{"x", ",", RowBox[{"-", "5"}], ",", "5"}], "}"}], ",", RowBox[{"{", RowBox[{"y", ",", RowBox[{"-", "5"}], ",", "5"}], "}"}], ",", RowBox[{"{", RowBox[{"z", ",", RowBox[{"-", "5"}], ",", "5"}], "}"}], ",", RowBox[{"Boxed", "\[Rule]", "False"}], ",", RowBox[{"AxesOrigin", "\[Rule]", RowBox[{"{", RowBox[{"0", ",", "0", ",", "0"}], "}"}]}]}], "]"}]}]], "Input", CellChangeTimes->{{3.6887586875324287`*^9, 3.6887588243061724`*^9}, { 3.688759200773491*^9, 3.688759201660813*^9}}], Cell[BoxData[ RowBox[{"g2", "=", RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"2", RowBox[{"Cos", "[", "t", "]"}]}], ",", RowBox[{"2", "+", RowBox[{"2", RowBox[{"Sin", "[", "t", "]"}]}]}], ",", RowBox[{"2", SqrtBox[ RowBox[{"2", "-", RowBox[{"2", RowBox[{"Sin", "[", "t", "]"}]}]}]]}]}], "}"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", RowBox[{"2", "\[Pi]"}]}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", "Red"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.6887592120408015`*^9, 3.688759249597642*^9}, { 3.6887592840082464`*^9, 3.6887593085218244`*^9}}], Cell[BoxData[ RowBox[{"Show", "[", RowBox[{"{", RowBox[{"g1", ",", "g2"}], "}"}], "]"}]], "Input", CellChangeTimes->{{3.688759262400465*^9, 3.688759268299056*^9}}] }, WindowSize->{958, 557}, WindowMargins->{{Automatic, 0}, {Automatic, 0}}, FrontEndVersion->"10.0 for Microsoft Windows (64-bit) (September 9, 2014)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[558, 20, 7335, 181, 492, "Input"], Cell[7896, 203, 3422, 85, 312, "Input"], Cell[11321, 290, 146, 2, 72, "Input"], Cell[CellGroupData[{ Cell[11492, 296, 3329, 88, 357, "Input"], Cell[14824, 386, 2673, 51, 484, "Output"] }, Open ]], Cell[17512, 440, 148, 2, 72, "Input"], Cell[17663, 444, 884, 24, 35, "Input"], Cell[18550, 470, 176, 3, 92, "Input"], Cell[18729, 475, 3069, 86, 309, "Input"], Cell[CellGroupData[{ Cell[21823, 565, 283, 5, 31, "Input"], Cell[22109, 572, 31316, 520, 456, "Output"] }, Open ]], Cell[53440, 1095, 150, 2, 72, "Input"], Cell[53593, 1099, 2779, 74, 301, "Input"], Cell[56375, 1175, 222, 3, 132, "Input"], Cell[56600, 1180, 1023, 30, 55, "Input"], Cell[57626, 1212, 711, 21, 39, "Input"], Cell[58340, 1235, 172, 4, 31, "Input"] } ] *) (* End of internal cache information *)