(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 10.1' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 158, 7] NotebookDataLength[ 56743, 1152] NotebookOptionsPosition[ 55964, 1121] NotebookOutlinePosition[ 56299, 1136] CellTagsIndexPosition[ 56256, 1133] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[""], "Input", CellChangeTimes->{{3.682162298501871*^9, 3.682162325348839*^9}, { 3.682162368157855*^9, 3.682162407137453*^9}, {3.682162779811378*^9, 3.682162786995265*^9}, {3.682162913758072*^9, 3.682162914125402*^9}, { 3.682162981816029*^9, 3.682162982711321*^9}, {3.682163068288478*^9, 3.682163072168304*^9}, {3.682163117819425*^9, 3.682163123968315*^9}, { 3.6821631619993887`*^9, 3.6821631649896803`*^9}, {3.682163206185589*^9, 3.682163235744801*^9}, {3.6821632867642937`*^9, 3.682163318716198*^9}, { 3.682163432806164*^9, 3.682163434231059*^9}, {3.68216347636509*^9, 3.682163476646859*^9}, {3.682163853177684*^9, 3.682163855014164*^9}, { 3.682164136304613*^9, 3.6821641977257137`*^9}, {3.6821642975014677`*^9, 3.682164297838538*^9}, {3.682164536817321*^9, 3.6821645559348516`*^9}, { 3.682164616711033*^9, 3.68216461740899*^9}, {3.682164724612938*^9, 3.6821647339667683`*^9}, {3.682164796029511*^9, 3.682164798950366*^9}, { 3.682164946774913*^9, 3.682164950685629*^9}, {3.682164998122178*^9, 3.6821650163879547`*^9}, {3.6821650807292624`*^9, 3.68216509565979*^9}, { 3.682165263947599*^9, 3.682165264065371*^9}, {3.682165372163426*^9, 3.682165374372757*^9}, {3.682165882168922*^9, 3.682165887111968*^9}, { 3.682165952942607*^9, 3.682165960548774*^9}, {3.6821661923012*^9, 3.6821662556626463`*^9}, 3.6821663146379347`*^9, {3.682166662751055*^9, 3.682166680601881*^9}, {3.682166835793501*^9, 3.682166862621798*^9}, { 3.682167588818965*^9, 3.682167606024461*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"a", "=", "10"}], "\[IndentingNewLine]", RowBox[{"sol", " ", "=", " ", RowBox[{"NDSolve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{"x", "'"}], "[", "t", "]"}], " ", "\[Equal]", " ", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"1", " ", "-", " ", RowBox[{"1", "/", RowBox[{"x", "[", "t", "]"}]}]}], ")"}]}], RowBox[{"Sqrt", "[", RowBox[{"1", "/", RowBox[{"x", "[", "t", "]"}]}], "]"}]}]}], ",", " ", RowBox[{ RowBox[{ RowBox[{"\[Tau]", "'"}], "[", "t", "]"}], " ", "\[Equal]", " ", RowBox[{"1", " ", "-", " ", RowBox[{"1", "/", RowBox[{"x", "[", "t", "]"}]}]}]}], ",", " ", RowBox[{ RowBox[{"x", "[", "0", "]"}], " ", "\[Equal]", " ", "a"}], ",", " ", RowBox[{ RowBox[{"\[Tau]", "[", "0", "]"}], " ", "\[Equal]", " ", "0"}]}], "}"}], ",", " ", RowBox[{"{", RowBox[{"x", ",", " ", "\[Tau]"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"t", ",", " ", "0", ",", " ", "40"}], "}"}]}], "]"}]}]}], "Input",\ CellChangeTimes->{{3.682167124175424*^9, 3.682167191200797*^9}, { 3.682167235318005*^9, 3.682167425068624*^9}, {3.682167470482944*^9, 3.682167483552539*^9}, {3.682167540038033*^9, 3.682167564571343*^9}, { 3.6821675973718863`*^9, 3.682167601471653*^9}, {3.682167782005065*^9, 3.6821677855422497`*^9}, {3.682168295462328*^9, 3.682168298360331*^9}}], Cell[BoxData["10"], "Output", CellChangeTimes->{{3.6821672895654993`*^9, 3.682167420369385*^9}, { 3.6821674677069387`*^9, 3.682167483836611*^9}, {3.682167545705871*^9, 3.682167564944242*^9}, 3.6821676018038397`*^9, 3.68216796515742*^9, 3.682168298753285*^9, 3.682214940212914*^9, 3.682219716864674*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"x", "\[Rule]", TagBox[ TemplateBox[{RowBox[{ StyleBox[ TagBox["InterpolatingFunction", "SummaryHead"], "NonInterpretableSummary"], StyleBox["[", "NonInterpretableSummary"], DynamicModuleBox[{Typeset`open$$ = False}, PanelBox[ PaneSelectorBox[{False -> GridBox[{{ PaneBox[ ButtonBox[ DynamicBox[ FEPrivate`FrontEndResource[ "FEBitmaps", "SquarePlusIconMedium"]], ButtonFunction :> (Typeset`open$$ = True), Appearance -> None, Evaluator -> Automatic, Method -> "Preemptive"], Alignment -> {Center, Center}, ImageSize -> Dynamic[{ Automatic, 3.5 CurrentValue["FontCapHeight"]/ AbsoluteCurrentValue[Magnification]}]], GraphicsBox[{{{}, {}, { Directive[ Opacity[1.], RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1]], LineBox[CompressedData[" 1:eJwVx30wmwcAx/EQqzZeiighaBLvb60JQiWeX6h5qZdM1NucsZP2yhliOo2X ddXz1mnXuupQW0bXMquXUud2XqY7uplSNtHe9E2lupn2YVnjpWrpH9/73of9 UXaMRJtCofA1vX34VIpfx9DVgPkvyta3t23xQcnT4p8fLBDhF+UDTXJbELUb vZmD/xGnpO3beb62EKl2Pt3lSAHjN/LkUQUHvl7S7MTYd5A1nGhcmc+Bul+S t7VCQ9q3Dw7X6XHQIcwR3jM3wkmml8viNTaifl3VmXthgl38ydzlA2yEWS1U NfbtwUg631g9wQLDsK72SAwDUT8OOCQfYyHGjv4w0Y4Jbsh4arxqL1Lvi9pi R63BoVOtRjl7ofyhpLHMnYV9XgUGWWIbiD7jUddq2IgbOqokcqxhz552zZrn INnH1Nm23goeh6Ol6gg7LKgnmpo6mPjrBbX+cac9ZC8bHZSjltBPqD9957UD eqRhJz5fsoBETK+eTHfCef7Es6UNBpxsPk6Y6XGGwkWhnDRngF3YJiRZrsgp naTftzHHbQFLP6PMDYVmFW6ZjmaYnTxet/i7O7zmtMsVgj2wVnUclAn2w0nM 6SsVmGIxT5JYEeKBXOnxK8JDdFyUG2pRlz0wVfikeT3JBClmrbwK+btwWwnw 4IqNMZfXE7071BP2rVe4HklGyLfseiXf8sRQbP9UxwVDKPPCDOc7uNga/mOl Kl8ftJsxuvZxXlBJB4NmC2gIJyLPSWneWOmduZwetRNEV9/a8A1veN8z4ueG 7oCeYmj6p3gflEbT3CuidbAa1Oc/RuMhV3/CIthLG2fr2GUzXTyY1IYGjURQ EOxb9OHSEV9YKG50Ona/JroGQyypRn7oYhWJ96WvEyp1e/DYTT906r1J+qX2 FREj06VcPnYA3Q0HdWWuKkKycLch29gfFKZ75FdVq0TNpTOjkbf9YSdnCB0l K8SlrSdpf57g473O618HD5CE6FDrwwyWAKqGsPZNCUnYzaf9rTMjQIQzt1jE JYnrqWdm5cUB4CyrEsctScLggsj7Hw4BakZmr7UhSVTo63ICFQQc33/M26FL Esk887XxAGCgK/eTFipJaFc3O5TUAd0JZPodLZKg9xcblK8C8fuZPs8oJPEm xUj5ZbAQWq2ipAKNP73K1JF/IwTLpfqsv8ZjaY/onf8KIX7uM2+jcdswM3Mk KhCnNh457dZYvfGdzLMlEHev3arU1dim5blMvhmIZa2aZarG0+vNt+hxQVAU tb3U0jhtyHyz/PsgnCsKpby1p6loY3s7COdTKkspGv8P5ume2g== "]]}}}, { DisplayFunction -> Identity, AspectRatio -> 1, Axes -> {False, False}, AxesLabel -> {None, None}, AxesOrigin -> {0, 0}, DisplayFunction :> Identity, Frame -> {{True, True}, {True, True}}, FrameLabel -> {{None, None}, {None, None}}, FrameStyle -> Directive[ Opacity[0.5], Thickness[Tiny], RGBColor[0.368417, 0.506779, 0.709798]], FrameTicks -> {{None, None}, {None, None}}, GridLines -> {None, None}, GridLinesStyle -> Directive[ GrayLevel[0.5, 0.4]], ImageSize -> Dynamic[{ Automatic, 3.5 CurrentValue["FontCapHeight"]/ AbsoluteCurrentValue[Magnification]}], Method -> { "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> AbsolutePointSize[6], "ScalingFunctions" -> None}, PlotRange -> {{0., 40.}, {0., 9.999999807047462}}, PlotRangeClipping -> True, PlotRangePadding -> {{ Scaled[0.1], Scaled[0.1]}, { Scaled[0.1], Scaled[0.1]}}, Ticks -> {Automatic, Automatic}}], GridBox[{{ RowBox[{ TagBox["\"Domain: \"", "SummaryItemAnnotation"], "\[InvisibleSpace]", TagBox[ RowBox[{"{", RowBox[{"{", RowBox[{"0.`", ",", "40.`"}], "}"}], "}"}], "SummaryItem"]}]}, { RowBox[{ TagBox["\"Output: \"", "SummaryItemAnnotation"], "\[InvisibleSpace]", TagBox["\"scalar\"", "SummaryItem"]}]}}, GridBoxAlignment -> { "Columns" -> {{Left}}, "Rows" -> {{Automatic}}}, AutoDelete -> False, GridBoxItemSize -> { "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, GridBoxSpacings -> { "Columns" -> {{2}}, "Rows" -> {{Automatic}}}, BaseStyle -> { ShowStringCharacters -> False, NumberMarks -> False, PrintPrecision -> 3, ShowSyntaxStyles -> False}]}}, GridBoxAlignment -> {"Rows" -> {{Top}}}, AutoDelete -> False, GridBoxItemSize -> { "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, BaselinePosition -> {1, 1}], True -> GridBox[{{ PaneBox[ ButtonBox[ DynamicBox[ FEPrivate`FrontEndResource[ "FEBitmaps", "SquareMinusIconMedium"]], ButtonFunction :> (Typeset`open$$ = False), Appearance -> None, Evaluator -> Automatic, Method -> "Preemptive"], Alignment -> {Center, Center}, ImageSize -> Dynamic[{ Automatic, 3.5 CurrentValue["FontCapHeight"]/ AbsoluteCurrentValue[Magnification]}]], GraphicsBox[{{{}, {}, { Directive[ Opacity[1.], RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1]], LineBox[CompressedData[" 1:eJwVx30wmwcAx/EQqzZeiighaBLvb60JQiWeX6h5qZdM1NucsZP2yhliOo2X ddXz1mnXuupQW0bXMquXUud2XqY7uplSNtHe9E2lupn2YVnjpWrpH9/73of9 UXaMRJtCofA1vX34VIpfx9DVgPkvyta3t23xQcnT4p8fLBDhF+UDTXJbELUb vZmD/xGnpO3beb62EKl2Pt3lSAHjN/LkUQUHvl7S7MTYd5A1nGhcmc+Bul+S t7VCQ9q3Dw7X6XHQIcwR3jM3wkmml8viNTaifl3VmXthgl38ydzlA2yEWS1U NfbtwUg631g9wQLDsK72SAwDUT8OOCQfYyHGjv4w0Y4Jbsh4arxqL1Lvi9pi R63BoVOtRjl7ofyhpLHMnYV9XgUGWWIbiD7jUddq2IgbOqokcqxhz552zZrn INnH1Nm23goeh6Ol6gg7LKgnmpo6mPjrBbX+cac9ZC8bHZSjltBPqD9957UD eqRhJz5fsoBETK+eTHfCef7Es6UNBpxsPk6Y6XGGwkWhnDRngF3YJiRZrsgp naTftzHHbQFLP6PMDYVmFW6ZjmaYnTxet/i7O7zmtMsVgj2wVnUclAn2w0nM 6SsVmGIxT5JYEeKBXOnxK8JDdFyUG2pRlz0wVfikeT3JBClmrbwK+btwWwnw 4IqNMZfXE7071BP2rVe4HklGyLfseiXf8sRQbP9UxwVDKPPCDOc7uNga/mOl Kl8ftJsxuvZxXlBJB4NmC2gIJyLPSWneWOmduZwetRNEV9/a8A1veN8z4ueG 7oCeYmj6p3gflEbT3CuidbAa1Oc/RuMhV3/CIthLG2fr2GUzXTyY1IYGjURQ EOxb9OHSEV9YKG50Ona/JroGQyypRn7oYhWJ96WvEyp1e/DYTT906r1J+qX2 FREj06VcPnYA3Q0HdWWuKkKycLch29gfFKZ75FdVq0TNpTOjkbf9YSdnCB0l K8SlrSdpf57g473O618HD5CE6FDrwwyWAKqGsPZNCUnYzaf9rTMjQIQzt1jE JYnrqWdm5cUB4CyrEsctScLggsj7Hw4BakZmr7UhSVTo63ICFQQc33/M26FL Esk887XxAGCgK/eTFipJaFc3O5TUAd0JZPodLZKg9xcblK8C8fuZPs8oJPEm xUj5ZbAQWq2ipAKNP73K1JF/IwTLpfqsv8ZjaY/onf8KIX7uM2+jcdswM3Mk KhCnNh457dZYvfGdzLMlEHev3arU1dim5blMvhmIZa2aZarG0+vNt+hxQVAU tb3U0jhtyHyz/PsgnCsKpby1p6loY3s7COdTKkspGv8P5ume2g== "]]}}}, { DisplayFunction -> Identity, AspectRatio -> 1, Axes -> {False, False}, AxesLabel -> {None, None}, AxesOrigin -> {0, 0}, DisplayFunction :> Identity, Frame -> {{True, True}, {True, True}}, FrameLabel -> {{None, None}, {None, None}}, FrameStyle -> Directive[ Opacity[0.5], Thickness[Tiny], RGBColor[0.368417, 0.506779, 0.709798]], FrameTicks -> {{None, None}, {None, None}}, GridLines -> {None, None}, GridLinesStyle -> Directive[ GrayLevel[0.5, 0.4]], ImageSize -> Dynamic[{ Automatic, 3.5 CurrentValue["FontCapHeight"]/ AbsoluteCurrentValue[Magnification]}], Method -> { "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> AbsolutePointSize[6], "ScalingFunctions" -> None}, PlotRange -> {{0., 40.}, {0., 9.999999807047462}}, PlotRangeClipping -> True, PlotRangePadding -> {{ Scaled[0.1], Scaled[0.1]}, { Scaled[0.1], Scaled[0.1]}}, Ticks -> {Automatic, Automatic}}], GridBox[{{ RowBox[{ TagBox["\"Domain: \"", "SummaryItemAnnotation"], "\[InvisibleSpace]", TagBox[ RowBox[{"{", RowBox[{"{", RowBox[{"0.`", ",", "40.`"}], "}"}], "}"}], "SummaryItem"]}]}, { RowBox[{ TagBox["\"Output: \"", "SummaryItemAnnotation"], "\[InvisibleSpace]", TagBox["\"scalar\"", "SummaryItem"]}]}, { RowBox[{ TagBox["\"Order: \"", "SummaryItemAnnotation"], "\[InvisibleSpace]", TagBox["3", "SummaryItem"]}]}, { RowBox[{ TagBox["\"Method: \"", "SummaryItemAnnotation"], "\[InvisibleSpace]", TagBox["\"Hermite\"", "SummaryItem"]}]}, { RowBox[{ TagBox["\"Periodic: \"", "SummaryItemAnnotation"], "\[InvisibleSpace]", TagBox["False", "SummaryItem"]}]}}, GridBoxAlignment -> { "Columns" -> {{Left}}, "Rows" -> {{Automatic}}}, AutoDelete -> False, GridBoxItemSize -> { "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, GridBoxSpacings -> { "Columns" -> {{2}}, "Rows" -> {{Automatic}}}, BaseStyle -> { ShowStringCharacters -> False, NumberMarks -> False, PrintPrecision -> 3, ShowSyntaxStyles -> False}]}}, GridBoxAlignment -> {"Rows" -> {{Top}}}, AutoDelete -> False, GridBoxItemSize -> { "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, BaselinePosition -> {1, 1}]}, Dynamic[Typeset`open$$], ImageSize -> Automatic], BaselinePosition -> Baseline], DynamicModuleValues :> {}], StyleBox["]", "NonInterpretableSummary"]}]}, "CopyTag", DisplayFunction->(#& ), InterpretationFunction->("InterpolatingFunction[{{0., 40.}}, <>]"& )], False, Editable->False, SelectWithContents->True, Selectable->False]}], ",", RowBox[{"\[Tau]", "\[Rule]", TagBox[ TemplateBox[{RowBox[{ StyleBox[ TagBox["InterpolatingFunction", "SummaryHead"], "NonInterpretableSummary"], StyleBox["[", "NonInterpretableSummary"], DynamicModuleBox[{Typeset`open$$ = False}, PanelBox[ PaneSelectorBox[{False -> GridBox[{{ PaneBox[ ButtonBox[ DynamicBox[ FEPrivate`FrontEndResource[ "FEBitmaps", "SquarePlusIconMedium"]], ButtonFunction :> (Typeset`open$$ = True), Appearance -> None, Evaluator -> Automatic, Method -> "Preemptive"], Alignment -> {Center, Center}, ImageSize -> Dynamic[{ Automatic, 3.5 CurrentValue["FontCapHeight"]/ AbsoluteCurrentValue[Magnification]}]], GraphicsBox[{{{}, {}, { Directive[ Opacity[1.], RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1]], LineBox[CompressedData[" 1:eJwVx3swGwYAx/G0sSmNRxavSmhCKBUTSTQeJb/EVBmVS6majUaF4YzzWNsM 3Zx51Ki1K5eqW65WRzNL7TBblapTXcx109MyxqjoQ9kmKurRLv743vc+jKQs iWwngUA4aGj74X8k+Kl7m4IIVh9pxWe+C4ovelzY/9ecwEqtia++MysQ1K53 ZvS8FLAV7sWqUp1AvLLrsck+AogLS7+4rW8KfHk5WXHRb+Gmf8JrAY0Ifbcs d+s/U4wFOtDjabugFmYLx2wtUbFXN/1zPwlHfl02mlx6B802rZwSrSXCaHPn GrqscYnsqEmOp8DOXFGbIrHDU4eBL3des4aESZmKY1JRTLUYGXKww4lxsSr6 rgMu7K4piEq3h/b7ooZSTzqYcfObb7fQID7LJ65dZEAiSs5btt4LF8aIR+as E3QWFESr6GDHROXoI5igsO6RYmcZeLZEvPz3DRcsnLeZ15OdQTp+uWR40xXx /jrq5x5MyI5Squ4nu+ELdvf8wgcucHP85Phouzsy1Qt89zxXMD5TCf+le6C8 J7Oi5uo+DAbSSemlLCxaJbOm7rnh0f1PFU8eeKJjwiNpeswdDivq9+SBXnDN rKC3b+7Hk1xZXHkoG2uSglW9CQuXlOY7iItstPLCM+rMPZFg08IvV3pjYDxm hubyLiZz26MsDnMwXExYrGZ6Id++bVW5xYEnszHyxR42tLlh5rNqLvKUXtlj 9WyYdkiMXY7xMMRK3WgmeyNcEHk+x9QHg71RJo5V3hC0da31/eiD+lN9pbUb 3tj9sHfkduwBiOW39lQncbAc3BWgMeXjQH6iO32Ig2oFo3S0jY9DEXimYHIR 4luQuJDii4ZiRqxOzkVbT6g90dIPUquO0/29XKzofwjRdPjhz47b15tXuZDI jQlX0vyxoFqLkDnzIJv7vT6LHIBOXSLJDDxcrKu8GzkYAE3lA3HBUR7qtmak E2cOgks7WyxJ4EH8fstUOj0Q+yvMuqdP8sCclT43Gg2E88RPG+RUHlpPVD5S FgZBX3SOeuVjHswuiH1eOAmwfCiuLDWNh3KSsZPooQCv0m4yxw3+kG+79lsQ 8Hz4HxOLdB52VjW6FimAlq+zb7EMpnQXmpUtAyTqOjvA4NcJltqaECEcT0+a iww+1UQ1Un4rhCqmcSLEYI10mnJDJ4SfIKMp1GBVHzVj4IgIJ6U7rh42WL9+ Tc5pFiFDqqOGGezY/FSu3BDBJKUzZtsjrxrvUI4FgzXzTf62pb22G2XXg6El V5Zsm2MlXn/zJhgvibFfbft/9Z2Yrg== "]]}}}, { DisplayFunction -> Identity, AspectRatio -> 1, Axes -> {False, False}, AxesLabel -> {None, None}, AxesOrigin -> {0, 0}, DisplayFunction :> Identity, Frame -> {{True, True}, {True, True}}, FrameLabel -> {{None, None}, {None, None}}, FrameStyle -> Directive[ Opacity[0.5], Thickness[Tiny], RGBColor[0.368417, 0.506779, 0.709798]], FrameTicks -> {{None, None}, {None, None}}, GridLines -> {None, None}, GridLinesStyle -> Directive[ GrayLevel[0.5, 0.4]], ImageSize -> Dynamic[{ Automatic, 3.5 CurrentValue["FontCapHeight"]/ AbsoluteCurrentValue[Magnification]}], Method -> { "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> AbsolutePointSize[6], "ScalingFunctions" -> None}, PlotRange -> {{0., 40.}, {0., 20.415184756511696`}}, PlotRangeClipping -> True, PlotRangePadding -> {{ Scaled[0.1], Scaled[0.1]}, { Scaled[0.1], Scaled[0.1]}}, Ticks -> {Automatic, Automatic}}], GridBox[{{ RowBox[{ TagBox["\"Domain: \"", "SummaryItemAnnotation"], "\[InvisibleSpace]", TagBox[ RowBox[{"{", RowBox[{"{", RowBox[{"0.`", ",", "40.`"}], "}"}], "}"}], "SummaryItem"]}]}, { RowBox[{ TagBox["\"Output: \"", "SummaryItemAnnotation"], "\[InvisibleSpace]", TagBox["\"scalar\"", "SummaryItem"]}]}}, GridBoxAlignment -> { "Columns" -> {{Left}}, "Rows" -> {{Automatic}}}, AutoDelete -> False, GridBoxItemSize -> { "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, GridBoxSpacings -> { "Columns" -> {{2}}, "Rows" -> {{Automatic}}}, BaseStyle -> { ShowStringCharacters -> False, NumberMarks -> False, PrintPrecision -> 3, ShowSyntaxStyles -> False}]}}, GridBoxAlignment -> {"Rows" -> {{Top}}}, AutoDelete -> False, GridBoxItemSize -> { "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, BaselinePosition -> {1, 1}], True -> GridBox[{{ PaneBox[ ButtonBox[ DynamicBox[ FEPrivate`FrontEndResource[ "FEBitmaps", "SquareMinusIconMedium"]], ButtonFunction :> (Typeset`open$$ = False), Appearance -> None, Evaluator -> Automatic, Method -> "Preemptive"], Alignment -> {Center, Center}, ImageSize -> Dynamic[{ Automatic, 3.5 CurrentValue["FontCapHeight"]/ AbsoluteCurrentValue[Magnification]}]], GraphicsBox[{{{}, {}, { Directive[ Opacity[1.], RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1]], LineBox[CompressedData[" 1:eJwVx3swGwYAx/G0sSmNRxavSmhCKBUTSTQeJb/EVBmVS6majUaF4YzzWNsM 3Zx51Ki1K5eqW65WRzNL7TBblapTXcx109MyxqjoQ9kmKurRLv743vc+jKQs iWwngUA4aGj74X8k+Kl7m4IIVh9pxWe+C4ovelzY/9ecwEqtia++MysQ1K53 ZvS8FLAV7sWqUp1AvLLrsck+AogLS7+4rW8KfHk5WXHRb+Gmf8JrAY0Ifbcs d+s/U4wFOtDjabugFmYLx2wtUbFXN/1zPwlHfl02mlx6B802rZwSrSXCaHPn GrqscYnsqEmOp8DOXFGbIrHDU4eBL3des4aESZmKY1JRTLUYGXKww4lxsSr6 rgMu7K4piEq3h/b7ooZSTzqYcfObb7fQID7LJ65dZEAiSs5btt4LF8aIR+as E3QWFESr6GDHROXoI5igsO6RYmcZeLZEvPz3DRcsnLeZ15OdQTp+uWR40xXx /jrq5x5MyI5Squ4nu+ELdvf8wgcucHP85Phouzsy1Qt89zxXMD5TCf+le6C8 J7Oi5uo+DAbSSemlLCxaJbOm7rnh0f1PFU8eeKJjwiNpeswdDivq9+SBXnDN rKC3b+7Hk1xZXHkoG2uSglW9CQuXlOY7iItstPLCM+rMPZFg08IvV3pjYDxm hubyLiZz26MsDnMwXExYrGZ6Id++bVW5xYEnszHyxR42tLlh5rNqLvKUXtlj 9WyYdkiMXY7xMMRK3WgmeyNcEHk+x9QHg71RJo5V3hC0da31/eiD+lN9pbUb 3tj9sHfkduwBiOW39lQncbAc3BWgMeXjQH6iO32Ig2oFo3S0jY9DEXimYHIR 4luQuJDii4ZiRqxOzkVbT6g90dIPUquO0/29XKzofwjRdPjhz47b15tXuZDI jQlX0vyxoFqLkDnzIJv7vT6LHIBOXSLJDDxcrKu8GzkYAE3lA3HBUR7qtmak E2cOgks7WyxJ4EH8fstUOj0Q+yvMuqdP8sCclT43Gg2E88RPG+RUHlpPVD5S FgZBX3SOeuVjHswuiH1eOAmwfCiuLDWNh3KSsZPooQCv0m4yxw3+kG+79lsQ 8Hz4HxOLdB52VjW6FimAlq+zb7EMpnQXmpUtAyTqOjvA4NcJltqaECEcT0+a iww+1UQ1Un4rhCqmcSLEYI10mnJDJ4SfIKMp1GBVHzVj4IgIJ6U7rh42WL9+ Tc5pFiFDqqOGGezY/FSu3BDBJKUzZtsjrxrvUI4FgzXzTf62pb22G2XXg6El V5Zsm2MlXn/zJhgvibFfbft/9Z2Yrg== "]]}}}, { DisplayFunction -> Identity, AspectRatio -> 1, Axes -> {False, False}, AxesLabel -> {None, None}, AxesOrigin -> {0, 0}, DisplayFunction :> Identity, Frame -> {{True, True}, {True, True}}, FrameLabel -> {{None, None}, {None, None}}, FrameStyle -> Directive[ Opacity[0.5], Thickness[Tiny], RGBColor[0.368417, 0.506779, 0.709798]], FrameTicks -> {{None, None}, {None, None}}, GridLines -> {None, None}, GridLinesStyle -> Directive[ GrayLevel[0.5, 0.4]], ImageSize -> Dynamic[{ Automatic, 3.5 CurrentValue["FontCapHeight"]/ AbsoluteCurrentValue[Magnification]}], Method -> { "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> AbsolutePointSize[6], "ScalingFunctions" -> None}, PlotRange -> {{0., 40.}, {0., 20.415184756511696`}}, PlotRangeClipping -> True, PlotRangePadding -> {{ Scaled[0.1], Scaled[0.1]}, { Scaled[0.1], Scaled[0.1]}}, Ticks -> {Automatic, Automatic}}], GridBox[{{ RowBox[{ TagBox["\"Domain: \"", "SummaryItemAnnotation"], "\[InvisibleSpace]", TagBox[ RowBox[{"{", RowBox[{"{", RowBox[{"0.`", ",", "40.`"}], "}"}], "}"}], "SummaryItem"]}]}, { RowBox[{ TagBox["\"Output: \"", "SummaryItemAnnotation"], "\[InvisibleSpace]", TagBox["\"scalar\"", "SummaryItem"]}]}, { RowBox[{ TagBox["\"Order: \"", "SummaryItemAnnotation"], "\[InvisibleSpace]", TagBox["3", "SummaryItem"]}]}, { RowBox[{ TagBox["\"Method: \"", "SummaryItemAnnotation"], "\[InvisibleSpace]", TagBox["\"Hermite\"", "SummaryItem"]}]}, { RowBox[{ TagBox["\"Periodic: \"", "SummaryItemAnnotation"], "\[InvisibleSpace]", TagBox["False", "SummaryItem"]}]}}, GridBoxAlignment -> { "Columns" -> {{Left}}, "Rows" -> {{Automatic}}}, AutoDelete -> False, GridBoxItemSize -> { "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, GridBoxSpacings -> { "Columns" -> {{2}}, "Rows" -> {{Automatic}}}, BaseStyle -> { ShowStringCharacters -> False, NumberMarks -> False, PrintPrecision -> 3, ShowSyntaxStyles -> False}]}}, GridBoxAlignment -> {"Rows" -> {{Top}}}, AutoDelete -> False, GridBoxItemSize -> { "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, BaselinePosition -> {1, 1}]}, Dynamic[Typeset`open$$], ImageSize -> Automatic], BaselinePosition -> Baseline], DynamicModuleValues :> {}], StyleBox["]", "NonInterpretableSummary"]}]}, "CopyTag", DisplayFunction->(#& ), InterpretationFunction->("InterpolatingFunction[{{0., 40.}}, <>]"& )], False, Editable->False, SelectWithContents->True, Selectable->False]}]}], "}"}], "}"}]], "Output", CellChangeTimes->{{3.6821672895654993`*^9, 3.682167420369385*^9}, { 3.6821674677069387`*^9, 3.682167483836611*^9}, {3.682167545705871*^9, 3.682167564944242*^9}, 3.6821676018038397`*^9, 3.68216796515742*^9, 3.682168298753285*^9, 3.682214940212914*^9, 3.682219716897274*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Plot", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"x", "[", "t", "]"}], "/.", "sol"}], ",", " ", RowBox[{ RowBox[{"\[Tau]", "[", "t", "]"}], "/.", "sol"}]}], "}"}], ",", RowBox[{"{", RowBox[{"t", ",", " ", "0", ",", " ", "40.0"}], "}"}]}], "]"}]], "Input", CellChangeTimes->{{3.682167295878901*^9, 3.6821673077958517`*^9}, { 3.6821673447377443`*^9, 3.682167372464654*^9}, {3.68216741169899*^9, 3.682167422360126*^9}, {3.682167503608988*^9, 3.68216750750482*^9}, { 3.682167549827054*^9, 3.682167567207926*^9}, {3.682167686552191*^9, 3.682167692597928*^9}}], Cell[BoxData[ GraphicsBox[{{}, {}, {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[ 1.], LineBox[CompressedData[" 1:eJwd13k4Fd8bAHBLskdEsgslW/blbue994ZCJMnyJWv2bEmbUhGhkoQSLSqK okQhW2StpFJKlBKpqLHEtaTfmd9f83yeeeZdznnPPDMqPuFbd/JwcXElcnNx kde8wz5L9f2KGdk2Rpx//1ThktsoxVk5Ge0Qa7A1WVSF95tTopSVLyAByx3p oQuqsNfKdRmf8i0k9/vFj3ccVbB8Au4zSlXox53Aj3fHVeHOMaN4Qqkd1e5y afEYVIXP3iXSY0rvkaV9RmpxiyocdAjx+qn0HblUXQ+XPqUKgrMX20N155BC oveMi5gqRDbGnRyN4gX5ZVWuL86vBt3rGS/LSwQBLUkXOKK6GobsAuU/mYnD nnlXi9l8FciWTVS46iUJveYeTW3rVMA8qrKU95k0bDVPUBC+oQxNPRe3hK+V Ba7WdxUMFWXQP7nIkNmkAI9jbh/o26QEVsOUT3rOynCwhLVL1VwRvuk8sJga VQEJ5eUvY00V4OBCsmHNAVW4npS5/JepPPAvHwnapqkOLlq3J3SU5IDvOZ+V cO0aEKjh3fdaShbmO++Eno3UAAtxx3TnCRkYci5ZpaioCRLsHatWd60E5omw XW13tCAyISijslYaMhcrBdmeOuDhL1nckyYF/8k32fZN6sKxtEqemcgVINKg NbraSQ9+xy/GXDSVhJ3G+X9ctfUhoqo8wXJxOew9vG+r8h99aMpjxzx7Iw7H f2S87uw0AEvioPGMkhjYaordFMkxhB3dKy9IVYrAqTsmDbbeRtD1b616znYh 6HC8thBtbgzy9jB2vZYfGp6Hn9TmNYFa0yhodOCDsuCj8KzfBO7OXviWXM0D ijbILKzcFFr2Hl3pYcEF7H/u8ZdjzSDlpMZyBeYCupIi+533iBlQli6pCl6/ gBak3tkFxptBBO2772OFBfRA21FWP8UMBrhfCaXMzSMNt033GrPNoMVRNDOl fB6JVBh/Giozg9p1JW4umvOoO3gZVfu7GcTKP/W8oTCH/N7WT1Q6mUPBx4Ny gms5SPxzR+agizkUuQ+Uqa7ioEc/35gtczcHxgoDSWthDpLgHj3s62MOhR3L /ap/z6B6LRnhZeHmIJc2UaNZNYNkj0ao+p4wh2zrRwtOdjPopZbKNtFH5sCy 23ty/bFpFGuiPWNaZw4l7d2yfXumkQbTNMfnsTnkXV0ymxk0jQ5v3zzwsNUc uJv7uvQdppH20f2hPt3mcKx1rbywyjQ68fZlwsMxc3gXzay3afqDGEePVXgr U8Aloj6HKvEHzdByvWNVKbD/k/ZID/8fdJdTIZq9hgJ1D7YPHv07hVQjvu98 pk2BNXwuU4sjU0jAc4u0qTkFysY85vMaptBrmuJeka0U+Oicy4iKmkKBnCqz h/EUSH4b9WT/wCRSKX/99WUiBR705tx90TOJesPHzowmU0Chm+Jl+GISbf6m /F3lDAXc9+1YVK2bRAZvky6czKWAHDt76FbuJPp732nOu4ICpxQbPC3cJ9HZ 8PFHIt8ocK6FapM+NIF6UK5AyA8KNLewPy/2TyB5cSun9jEK/PeqvOXg2wlU eDf31/EpCiRWc9tda51ANeNWKv+4qXAh69x4YdEEGt59KXFcngq5T5ODYndP IK0Nm17bK1MhJFaZeBc6gSJWTCmVqFLh3EO1n2z/CTRfvqkqWJMKW85z/dvg OoHEp6d+DppS4e3Zl5/H0QSi7LNxeLuVCpe07xwYEZtAcRun84y2U8FItPqu mtAEeiJz9cdZVyrkdIrbhS2ZQPZV0/H2XlRw5KHsWTc7jvzmrj5o20UFtYzn u+mD4+h0LEfu0QkqcC9wOE+rxtGXIzeGL9dRIa7k/Ohs+DiK6cqLePKYCuuY 7RZLQ8eRkHLW3MgTKgSxezTlA8eRQX3iMoOnVOjU8N3u5zWO4v8GmDzpoYJm jJKfn8M4UjugmThCUGEunb9vhfE4qmxfLS46RQU7jTa9Mv1xZLtKLkd/hgqD 2msEXXTHUXSlcMmBv1Q4pHOct3btOGqeHn0jIkQDlWbX7T9kx5H/7lJ1fVUa SAcKGjzhGUdzjYWlTmtooKbYpPTlH4FOS1wxP7COBjOXuNL4/xLowb0zm5vW 02Bo96f7gTME4ici9zjRaHA5baXfsVEC3Qw1at7vRIN928Na+d4RiFajbX/J hQYGXZdOnHtDoC5h9feN/9HAfpH2dt1rAs0WSY0J+9CAw1btD+gkkPX36RWX wmigLUxfqtdMoJ87q/wak2gQ+ctX6+t9Aol0sNI0U2lwMTxl+OY9AmnrPqs6 e5oGr+6HWEWXEih0+uMy30wafD328aFiMYFGE5dU8V6jgbwm37GGfOxCe1GL OhqISv7WW5eB44m8N73zmAZ1PeFHt6TjeBE+PlLNNPDZWxV1KA3HM9vzcOgp DehxaQvfUgk01pbjnfieBu6mV8QWEwj06/tQResUDaa19Dep7CPQMrvwgfUc GnwbEWm4GEMg3TKO0Pl5bCvNcdk9BAo7IOwVyEMH9x+xdmuj8PNC+kKC4nSY K29JjQkl0G+t2B3WWnSgi4y2DnoTSPwMX3KZLh38mgb7rnkRSG/q9H1ZAzoE qd4bDPAkUETNVYGfZnRQNjbOmXcnEGHbWpZqSYcPXLVRW1ywd0nwP/emg+5G 35UF9gTydPgeNu9Hhzw5D9d7dgTqNGp4uy6QDkWnYmLqNxPo9vyugsQwOngX Cat+tSFQYHK7BRykwx37t6F2Gwk0cP1Iwv1zdFhJdV28zCSQ/Qnn0c/ZdJi0 rnZ6AwSqC9HdJn6RDh1Z28VFsXMN+lR3XaWDSZtwZjyDQC4Npk1rSujQ3Dul dpJKoBcffnHntNLhJOPrNT4TAjHqm4PaOuhwcVFfaYcxric/9+X0czr4C5Rm VxoRKCXI+qpjNx2I3ybBuw0JZDlzA0Q/00Gq89D8Pz2cX8Ij7sgcHdrq/y6G axNIZ9rwW8lfOhgqrqke0sL1vBey7+diwC9FU3kP7ANXKhUp/AwIFqts3KJJ IBPdFXWTkgwQHaqUsdEgUMmmpwv+OgwoX9v1NUmNQFof7+46rseAZa7FyjLY N6OyPl4zZICJ2o47t1QJlH/Ru2HAnAHnv/0q6lpNoMyxmXg3SwZs6jJ/paNC oNh0NWE7TwZEZrndUVHE868uFBvqwwCp0p9WzQoEiqn+PZqykwH+2gOzQdgR g9WdrSEMuL9gn/1AnkC+xg5nmfsYIPvlrb+7HD4v7w+tMk5nQN8Tlec9MgRq 3eWb4ngOxz/y9WAq9gaeTfOR2QwYMNORBGyG1or+kjwGnPhe8rt4JYH0DxVd 0ShiQPQDQ/HT0gSSUelZK9/EgP03aU/3r8D1V9Scp7QwoOMNtFOwJazzBV3b GZCuYnJiQRKfn927fma+YMD1zlZ0DPvfE95SsT4GvM50cMqSINBwoJ4J7x8G zEuvV/8kjud7eaHVTg4D/k4+aC3A5lQruLbOM0CoMFc9HFtQVDg2lQfBsLWd EA+29r2hx5LiCNYWD5roiRFo92zOZjUtBOKnBGVvi+L1zF/umaiLoP50qedh 7ASbExEj+gj6+zfQHbCz8vZk3DZDwOnjr+KIEKiaueW9kSWCG0X7jthhc6Xw 7bTwRnB10P6RpDDObxgbU+iHwLb/euEvIdxv30SSYCCCC6OSYh3YqroDRc/D EPj28m46hm35qppwikWwuaeeyREk0CnZiFj/bARMXQ3hWQGcv2n4VFsOgvYL 8n97sS+FelzWvITw91GYYy323Trrxl/XETTapq6Nx+72URfYW4agNpqwksSW Le7NSHqOYEOCZMgGfgJFvb70pbQLgScr87YOdvu8j/671wgSFMB4JfY+25+d Gr0I1oselvm5lEA9Y/OCHcMIPnFnimZhW59fBku5ACLfHp+Z4cPzdjsiIIcb YNunLMsh7PyGV6d1eQFsmOu7XmHPjWT1b18KcJFf50wJ9i1zpYOFIgCfPPtL A7EFPqx/sHEVQN/34lfDSwhk9Du9v08WwCDu6N4ebC/eqSWR8gDiGltt27Ar tSodLygB8HyrjC/CDoiF8e/qABKxQbKR2C0KDtqpBgBXpBp3LcWe0L/vqGQE IPRqzctZXgIpWkodvG8MUJDR898Y9p6w9+0fzABc7MyGurHV6r0DtBHAaQW5 mhvY8V5R+c9tAATiTsnZYU8tyufEbwawND02yML2z2tNN7cH4N2U/NIU27pX /uiNrQDC5QUKq7EltrV6HnLF/uccwOEh0BUreQWdAIBaZot7MbbEcMuKwUCA 248ibfOxExIiRS4EA0zq/na7gB3wuGV+SRiub/TeqyRsXWpkb180AL/uORV/ 7FrdluyTxwBoni+H1pL3n0eksRIA6pRdM5TJfCFySZzjAOoF7a6ryPg3I2L8 kgGOf5TSE8a2XS3nRDsDUGagP/ubm0C9UhESo3kAXXNnK2qxbctlhfIvA8RB 1NOH2HVbm7ldrgLQz4X8uYd9NU12vOk6QIdPYUIBdqBg84uLxQAxicrJZ7Bn FladsqkC4OyWHg/AJiZtkmerAbwLtoMv9siPQ8cLawCeSBpf3kHmf/f5EG8D wMEtIqlOZL7yW2GPWgDaNL5IbcBOCDXfotUNkC8q26mKfcg3xPb9GwCpknIu ZewYt7yNST0AbgNjSB47YCM3c7AXAL27OrgC21qtXf/iZ4CdsZxHS7HF+1wk hX/j59uvpI1y4fP7OkWsigCQTog/8B2bp6NGOGAC50+6HT6MPflQeUnTHwDb +9R9A9hvM0amDiwA/FT3WezGzrXZ9/aHIBO0gsWb67EzWUWvzgsz8feemUUt 9mnzvk5LUSakmtU9r8aOWwut+eJMGAjZM1OB7cMrUOm2kglyAXuLbmNrPMq+ 0KHGBNfM+JiL2G2REeX/rWECpWf72QvYgRqbXoyuZUI737LybOxbmXNLxLSY sEXqq8A5bM0IjwhHfSZER2+dPomtpa5q1UdnQs+w48s47Gcf5r13ISagghD2 YezQs92xi8CEel7n6ljs2/+OlylvYILPcb8H+7G1e0cUdlozwXJt+dBubN20 kqmx7Uyo8dWSDsLutEgSi3NhQiBtvDMAO2zBU1PcjQlXWitS/bFLg5Z7GXgw oSpXYLkftt6G6KcxvkxoeBW9yRPbYNY8/184zjcd7b8d+1WpRO2ZSCa8mXV2 dcKO8v/Zo7Ib9x8RZ78N+/7rXNENMUyobIux2IptVPJv34lYJhy/vcTZDtvY t9l+eTITwj2TBayw5a8PU56kMCF4T6+MJbm/Q/xr9p5kgoXfNm0Lsn5/64W+ NCYkV2W4s7H9g1/cvJnFxBv8dwBhn498zw3XmVBhEBNpjn24bG504gYTIroz Ss2w/Sbl3t0oZEJZCIswJfvbs6NEuJgJUzXr95tgt+8bdOu5x4RjVc43DbHn 4sbuh9Xj9dTbGq2LPdCw7LLKYyaYPYjr0cFu4dZL6W5kgqRkHY10RnyUF6WF CfRmR1Ftcj+SZkSWPsfrYazXsg7b4zS3f14vE8QzpNPVsTe8WO2wpQ+vv7Sj BGlN8Q003o9MGEaUTDXs6fREiaDPTLCe+XVZlZzPTOF6oxEm6CY7tqlgN+RK yTz7wwTZ8QUvRWy7vH4Vnxkm2Pwz4yb9Ie+GFofDhJJD1/MVyHiXjJHaAhO8 H8yNyGPrXHXyj+VhAf4njJfDzrmRWa4txoI8xUc/ZMj5LvCoaxRngfM9zVzS FQXqbS4SLLhy4Iwd6a7Cig8JUiyI1zJ6sBKbv+gNT78cC2wW089IY0eXSDmc 0mBBkrpz+ApyP0v73VQ1WXBHaUyXdFrpDb8qLRYY8dv/lsQuumu8b0iXBRNC oVGkP5c5XaYbs+BTX8NRCbK/h5ljY0wWvKs+1yiO3ffQYyaezQK3ospTpIMq 1bllLVjQYBrvSjqhqmKF5UYWcMqXTIphP3r0hnrJjgX06QQd0hoNUil2/7Eg 6NC7OlFs2qqd5/vcWcC9xPMcaYeo8oLgHSxI5A4KJn1AzbEp0ZsFmS5cMqSf nTizUBfAglxz4f0i2JEOwuHro1nwJZjpIIydWOR2qHYPC7hM2Lqkc3iLUm32 sqDxxD1h0k8qNt4MOMCC8o7OViFsadnEz5ePsKB7eGQD6eovXNvET7HAuDzB QZCcf+oWn0unWXD3UpUR6cFzlyO0z7Dg9RorGdIiVoxTGzNYcEYqYEAA27P4 YMuRCyx4tT8nmvSS6Gmz8RssCPwjfpsfW+a5hVVcIe7X0+UcaZ01mU6it3B9 0d9jSW9/Zxi17jYLVmkM2ZG+SYss9i5jQZqF5vRS7M18owqv6lhgt1LDgfTN kfiFow0sqHk5TyXN80zug34jC1ybpNaSrjxrcz69mQUl8zKLfNirlYslHJ6x 4KPIplLS05Rg/q73LFiecVmB9BZFnuG4DyxIVT8gQrqYO+fJ+n6cP/b0/BJs r7b2o2kDLEje6t1LusNp3YLdN3w/60IO6UsRI8TzKRZYSRWqkeY4HnlxaBrP u22jNGlHU5kSHQ4LdPhFBEkLLFqFnJpngX2/xi9e8n2WWjhky8OG93/3VZO2 KvDvfSrGBoHjc26k85MXKw8uZ4PC9lp70guhWdlakmzoFC3YQLrMsGVbqjQb /gh+0yUt/1j9hbUCG2ZUZfhIEx++NrVrsqEquO4hD/bEYdXH7tpsUCmjlJKe UvGp+63DBoO/XQWkOQGfKqX02UC5bp1Fmmvq/R1vUzb4eYbEkBYXeZE9y2aD 3FstKmmJUtHMkxZs+CKoZUR6xVbbs0pWuH5X0CEtc779pIU1GzzsbyuRVlZ7 cjR9CxvKmo/xkdajVYWs82DDvNpkNzf5/v80E1izgw3JDSYvSBsdM/G392KD TtrJdtJmbfe99viyYX3FrjrSsK3E6XEQG2o+z90kvSX0GnKNYQP/5MqjpLeK faGN7mXDcKtXLOltZcqUuP1suLC3ci9pF06e0Y1YNnBaMsJIeyWcXzd+jA1D j/+4kw7PPSWZnIbXX3gDjfQzbrXf1ulscM/5Y0paM6C6QySDDQcSHxiS/qr/ 7eiZLDZwZ+3U+n/8ViCy89gQ8F+UHGnWxOSzgmI2nLu2bJEL+7Jzys3AO2yY /ZAzR3qhRjlBs5QN71YbzpB+mLSZWlLGBgmTJIK0tkLhrYoqNkyMVA6Slt7o lviklQ2GRorPSO++Q3gntrPB6ItzO+kuiST6xqdsmFx6sYV0av/9qaedbLBY Rn9Mmmu3qO/rN2yIpI8+IP0jrwF9GWRDr9SXa6SteJ3lrg/h+HDjKunrgWPT O7+x4a9xxGXSnoayJd9/sEF9rcJF0t1tu+XHx9mQmtV1lrS+rhCnbJIN+hkv zpA+nXHldfQfNnRcfXOa9KYdz1M4HNzf+ekU0gVNvv7Vc2zwcVieTHrJujlm 7AIbvtYbJJH2Pn1GgbHIhuBnbsdJ102umf33jw2NnsnxpP8HBVLW2Q== "]]}, {RGBColor[0.880722, 0.611041, 0.142051], AbsoluteThickness[1.6], Opacity[ 1.], LineBox[CompressedData[" 1:eJwd13k4VV0UB2BDKiFTKlKmqEzhXuGOv3tvKNIkJaVMoUg0+kSRKaI0KqFo niilTAkpaaCigSglIiQlFOXb21/3eZ91zllr7b3O3c/R8ti0dJ2EmJhYrLiY GP1N2+Ux2tTrCk/MaNTXgfMXeekunawVmnH8Id2bO/+lxvDrHOI3a2qe4Nfr 3V5zPfUYf4ftyvFSmpf4lfPvTlqTep5vU47V/Rr5fPVDP/zen7zNv7aHGfld o5KfE/q1+HvKQ/5H96yJXRp1fJF25fUpJ17zdy7xc+vQaOf3j5fTWn+0lS/9 +2Slv/Ef/piNjMCB3f38oLLdCZ2bJfEz6Vv54wBxGJ89/OJWljQclBe7Vg+N RstCX/UPlgr45WSuYP9PFslqMVMz3JTx/nDQsJydIqw252VLPp2I9sDijgaH Cbj/5uTiTTPUcDXA+7QnJsE04R9v8vypQOOMSo2/qrBtZX0wWaGJEI7YT+ce dXwxum3d26kF8SnVPWezNbBzKI5RFKKD9LJ3iabGWhij2LZ+mb4uzqbJ6Xmd 0obUMylbmbt6EPvQUZzwSweDVdf8DwXNxNS0XIcz63TRsiJLddo0fXS+Wy/V +EQPgr0BGx9dM4DihBsabybOxNF/edKitUYIq2wrHRMyC6vU7y9o+GkMmy9S HQdK9SFbYtCp7WSCL1rXzGM0DLHOPPPXSkNTcGboxHoGGWHHruClmr9M8TR3 0qxr540R/fVwTVWVGRx2XzEJ+zYbC/TlL8qmMJD+O0xbI8EEidfmlCxwZ2Kb 2PiITaqmeOx4ZmirlTlEGVfuBF81RcmzTQmGknNgJjIPCTM2Q86GCDxtnIO1 j1gep2+YYZo93zLglgVerjii0qnLwLdSv7/771jghozZnT0zGLhnebwsu8AC l7JF4ZNnMeCm17Ogp8QCJQc/RwsMGcgUP+O+9ZkFdomk7ePMGJiRN3pfSKsF eqOE6t08Bkx0qxpjJltCPvVN9fgVDIiGV0eeCrXEvoHwU9uiGDgdr9YuGW6J lCXOF7uiGRhSebvQN5Jcn9a80yuWgduGjmqm8ZY4KX1y98J4Bma6zL9RlmwJ qQLGR+UkBmRzzT+05FiidBfsPE8yULthPNuw3RIXEiRzx98g+fuenE7qtESP 2Be/NTkMJEbEjf7VbYmvpp9fXr3JgM3xUS/v9lmCETcpy/o2A3kP/vguGmUF F/XaSJ9CBlI1vxzbrGEFMe5JGb8HDHi9vvcjz8kKSo/V5mx8y4DCx8dHm52t sEV9CWtXHQOFHa8sx6+2Qohbwdj99QwoiXfu8vSwguXdTxKXG8h6GkyWGb/J Chcrcj1rmhhQiwjU8dxrhT7Fu0kdbQy8MNBaJldoBcXgkubZAwyEzjHstyi2 wvnD8ddUf5P+BRYpHqVW+Oatw5b4w8Cu5Q5NdyqsEKjvn1E9yIBhxH/+HrVW YI49d3zVMAN7X7+IutNlBfUEfTvL0UzwIvbkumuykDzpU4+UMhP9nFT3UB0W 5u0flV5OfH0gVy5ZjwWDSE/NPROY0AlsX/fUkAWn7kWXfqswMXbt4okWViwM Kpu9rp/MRA1n2g7ZpSyUMRIkvKcx4TuQb3knkoWf71WLz8xkQutWzecXMSxs rp1pypvFRP2mrqTOOBZeFG9KfEPs8EWzXSuJBSmfmlFjDJgwex17IiGVhcl3 GFIrjZj4e9Ppj3suCxOrJ0g9NmXi0KaeQtkvLMQXKl9utWLiDT91rN9XFqbn RdlsYDGhrmDrVNnFwqPZd4o7iS9cT/0W3cvC3Hmq876xmSjqsdUaFmfjiN6H 1mYuE61b0mN61NnI2BWenCJgwmDu/JpFmmxcOHkqTVbIROCEXo0sHTbSOXrJ YcSDt+bnb9Bn40ph5obVIiYU+no7mi3Y8IhtKZazZoIVbL/k9VI2XOvaQzCP id3z+tKYy9kw3OK7JZO4fHLG10Mr2eBuywiWnM/Eovy+yEVubAgHCrPLiL3+ ZNx+tJGNps/Nf83tmdgfOjClcC8babWh0T0OTHwKP9d6qpgN8/1x5yqWMrH9 eVpgeSkbX1tLz0g7MjFO89iftnI23HqnD9oRm92LGW/2hI3rjXsLnhBH/vWZ U/6GjRsFyaceLGNieoh+TNt3NgKO3dBMXc5EXqW2glwvG/U3vO1qiReoTkkx 7Wfj7RuvjzIrmNiaJ5MV8peN/tMKq/8jftDX+Up2HAfFub3nHJyZ8N6SrWuq w4E3s/VEw0om/pRdyHbS46Bmqv9MWRfSj9Jpq5BZHPzNmGHBJr59I8nh/mwO elpy6pKJx3wP2ubE4WDTAuFS+1VMXPRnPvjPiYPsENtxqauZ4BQZLkp35qB/ MlurnPi5jG5d2SoO5GtqKzuIf19W6ZLx4KDJaLjGypUJu/a+CekBHDhHLWp7 TtyxLt+rLJYDVnRVT8caJmQfCw/o7+PA3mK3n8JaJgyNn+Yf2s/BlQW625nE /n3vx3se5WD2w3fCncSdMaPyJc9wYNaXVyfpRnxhkZx1MQdZ3fMGxcn/sKxs ncW1Ug78yqbKahIbBnp4qDzg4PTzS/e4xP6W2+60POHgR9d7pWDirkcp7jF1 HBxcEh7dRvytvSW3opeDRvPQ//I9mBi/cFPT7AEOrnnC8CWxcc7AuOODHCzI n3XwK3FAiIybrwQXyjfebFXzJPePMx0nrcBF4vnDhduJuw1C19gZcJHHeySl 7UXmNUkqLseYizHJ7ZhDbNK7/6aaGRcv6/zZdsSBRRljOyy5sMuMCgoi/r6g ImefDRefX6/vK6LeqDTmmTsXd7Injp6/jom1S9oDBr24CK8+u8OZuIpZ8nqW LxeH/NILfIivDm48HxPARcSP17ejiX3jKq2xkwtZ06Bj94ibzoZH3TzCherS Qg8Db/I+7F3R+TGZi61ew3qWxMV+xssUTnKxR39r+VziVLMGnY0ZXGjrnU5w JXYusbivl8VFiH6tayJx9btv4ikVXJyQNvNoJebde7D+0WMuQuPEO78TX81M fdH3jAv+9cNOg8Tx6+0yHGu5GPfQt1Lehwmb/nOQ+8iFZcNNrgVxsZLr7vA/ XDRvKLgQTmzUx/iS9ZeLSQNdK+KIU+vGLWoU4+F2xIaug8Qhp/Omscbw8M5o Q2Um8RzjCcU/lXnY6pnrcJ84a/6TIW8jHnQUx1QOExu8v74x2oSHSy5mz0b7 knnffOz9GQYPwW81i+SIM0+6lzRZ8XCr22TVFOKjXf2RLjY8yEwPP2tOHHpw uszCtTwkr+0+5k38W3dcqL8HD9PYVzz8ibcXdHfGr+Mh//Jz9c3Egc0FVRV+ PEyJ0XUNI/Y0X3JIEMzDSwMN5cPEdnVhquYHeehrTwwrJK7Y6BnveIQHDd5C mxLiuRLzB4OSeTAs2yT2gJhnMKExK42HzLeH5lcRm4ZdPj3zMg8rczYym4gn a72ZoX6fhzZV1QPi60n9uUXHWQ95UE/Kj5AiVrLLlF5ZyUPQfravNLHslo0d R6t5+OMarKhIPFwumS3fwEPXx09vNYhbfU3mSP7ioXvNpFVs4m7FC7brBniI XF3xiUc8UDB1ZcUgD/LC3DVCYmk5mdB9EnwcU/vAmE9seKOlVFmBD0nn/1Ys J97yO8VhugEfkROZZwOIQzMV18YY81ElUIgIIo6y3xvYZson5+G6ZVuJj6Vt O3zVko9Ji7vf/0dcIFhcx7ThI8TuakwUsVi81Dprdz4qdhuxjtP8jNDtF7z4 cB5kJKTQfht+xEr78nEoPfJ1KrGOcdPlZwF8bLzrsSyD2OZlwXenUD4GxqX0 XyJOVAsM9U7mo0CjqDuf5r/fmvgohY/L17RfFBKn+7ue0k/no+QD/+pd4uvF dmXfzvKhfW9gUSlxrYfu2B05fOgsLFv0iFjtSv3h2Gd85MhexivizTXpn7Kf 81Gk+Pnqa+LKQQ/TtzV8xERcmvCWOHhBR9XMej5cxii8qSd+0zUo/biVD+nt jTObiO2Oj8doMeDpC2PZr8Tbrwb6pIgDP7jbLDqIM0te7jeWBFzqFNd2Ev9p O9a4fDTQ5xBy5hvxJSuNnRdkgeh/kPhJPPbd7NvzVIGJ2R5ag8TM7oONDWqA YTlr0hCxm2TvqCB1IB9K4/4S5xnkOZ7QAGbhU+s/Yp9Q9LTrAi/G/AmT2MDE w6lLDPeZAV3m52dJE/8wvemowQQYmnE/qafZqOy8aQ6YZvcVjCPeFlBX+c4S MFBLEsoST7/n7mPIBxxNQszliSPdNmc+swdmBEy5PoG49596SqQDsH9flqcK sXdaxUGrRUDI6KYJE4nt6tUjzi0FPpxC4CRipWUVa8NWAq9MWmTUiE/bqk81 8gH2SQ02TKPx1ocTmn0BXUfZAA3iqKgg2RMbANvZJn+pfUofDo4KANaFPVbR IjZmB9U3bAWuv5Uy1SG+a/wwOWEPMPh9rMoMGn8WeEAYBbQ5WMZTn/abEjsQ DVi2v/1LHXUxcLtXHGD/ck/TTOIF2lOcOEmAiansEX3iepVApc40sr6Lt+YZ 0fgttXGZp4CikCmTjImLlz4Qd84Aaj6Gb6XOOKDWc/8sMDb/mMFsYl/pB9Un rwAPa9z2mRD3D6km2ucDxceHJzCIv/+0j/tdAPxsGO9G3fY1LPpCEdA9Y/pl 6vq3H8MkS4Dx5Hxi0ny3LgUUPgQS3zxeak7r97dabFALtJe5rLAgDvP0W1D3 CmjiM45Qb3dJmxf7hqzPA8Fzap954oLmemDLgXE2lnR/pleanvwIyKYlzLQi VmhwVpbpBl43rHnCIpauiZfP/w4wp+QMUUs8LpLx+QGc6XAzYhP/vKM56v4v Mk+n6hOoXx9u6w0ZApy0wuZyiFPtg19/lRYg7KracS7xUeHll8dlBNimrVpC vd+qocpGTgDriwe/UO+egYpMBQGYv/qZPGIPybF5LpMEeKPwpJJ6ZmHyicfT BUh7Kd/IJ34UFHhrlZ4AX00cB6l9Z86v7pwhQMsCS1UQXzr6Z5S8gQABwx2O 1PqBroGOpgK4qk56QG2gq2PbwBUgeNuGZAHx03eD7hv5Amivnn+d2v9Qbeg/ CBA3+t0j6qvD0TmacwVQblP7TW1Y3zZ1nZ0ADrYzlgvpvB3I6u1aLsCaGdqS IuIq61j53c6kP9ljk6gDhtbqK7gIsHbUYwPq7PWKbmaupH7HbY7UJnO3Ptnu KUBkcPcparPfVpnDmwQoH6tlPJf4ZbbS3aQgAS5J/eRRb/bueKO1RYC5BUmL qG/WpMrN3S5A7AnTQGpm1nDw3lABLjqrZlObez5YpBgnQL66xHRrYvWzrazy eAH0Df+aUUu0jNHbkSCA+bUHAuoqb7uhhgMCeP8qXUPtvaH64sVjAviW7zxG fTyoThxnBSiS2PaXelfOn84f5wTQ4R2WtiH2+jnl7bkLZL/W71WhNtm2Jkvm Crm+6LchdWVws8ubGwLMih3jQv1nd9fNgHsC/F4z7xp1U8n4U1qlpN4jKneo H4qbxNeWCXClqLSE+nDkZjfWQwF6dubVUBvG9suOfibAmRb9AWrX/eLeafUC 5L4MZdsSz63WXrK4QYCryutF1PoKczmS7wWQvcyyp+47GKO0/qMAtZsOuFDv Pypzj9kmQKaNdDB1SarK5Ke/BFjcuiuLemFao5ZHvwBiBjtvUb9LO2cwMCDA 3mCvgpHnpZvzpw8JEKQx/IDaKMPJO1RCiGrzVQ3UKeeO3jKUF2LAoVlqHp3v 867FZQpCnHB+Oo4697zuI2clIfKHz8lTP7+Q+y5KRYiqJwaq1GMuv5JonCLE u8FoA+qtWSpLEmcKMaXkhAO1RHaji46+EMMep5dQH8g+55VvIMRfmxQn6svX zYNbjIUI3+nrSv0xx+kU11wIhsUWf+qFd452dQmEuDrWey91wx3X/kiREH5P bu6jXp+nK65mLQSm9++njsrPnWAzT4iz8muPUhcWvmKnLxQiSvFOxkh/JSrx C1cJYd+ckUfNUV13vGG1EGKh8wupl2y+dX7DGiH4jW13qUOmO96PcReifb38 feqne5OGin2EqLz37yl10BKZTbO3CpHkufMDdcxll7C724TItL3ykTpF8vI+ +x1CMve1zdTlufMu+oQIEZ2k1EY9US3m46lwITb24Dt1wSexZQqJQhhb7Rim rmIv9kjfL8S4DpH4fOLmI6cCDZOEkFaWkaSWteUlzjssxCSVuNHUa6/sfBh+ QohzTdZy1KO29ln2nBNCMbRKlXryM2vb3ReEeCFlN4XaSO+ok9wlIZ72lqhT L3/L2DzrKumflaFBfZETdMU9R4jT7Ya61A5SnVNfFpP+xofNHom3RQ5FlAjB rak2oZZ4OuWdaZkQXsxpZtR5h+yPH3wgxCrxLCa1tuYVpSVPhTASP2lF3cfa MOZ5nRCLU3KF1IunSbTufifEsp4KEfUV8ZTy2Y1CyF19PZfa7VFlxIEmIf7t arehfuw0a2jhFyGWDrbYUacHtn1/1iuEz3TrpdQDjuHVYX1CyATpOFI7WkzO Mhog8/bq34jH/rP1Sxwk94+95kS9ed+FlgUSIlwc3eZMbXveu/6JvAgvv31Y Q50Z9y9vpyJx6eG11EP+x5INlEXwnWPtRp3DeLhs30QRHupluFOrl+pW200V wV2F50X9/d3n+5X6IrjdUF5P/WOXTulqQxF+dZ4cca+WR3G3kYh8h2lvGKnf 50OeiimJd+n7UYv11l1ztxDB79PMjdQKstXJv0Ui2L5rCKRWypY7mmAtIufD giDqCUsXHNKwFWHm4cIRTz5emWBtR+qvObKZWnN6ecTBxSLs9bfYSm3Cyfeb 5SrCfW/soDb70O9btEYEBaPUETP3zPFe5CaCwfO+EVs+uum2zVOE2rsXg6mx LMupdL0I0zqG/xvZT/8z/JXbRdBX2xVKvVT+E6dzhwhNCZUjXpajydr9nwht ihPCqJ0H0pjnQkX4vfr8iN2ijs/q2SPCCm7JLupNqYnKcQdEqHd+Ek6t71Pw WPYweb7psgjqz6ZfIpKOiTAHdSN2rsD35DQR2re83UMt/PHz6fkrImzPzooa 2b8izSj9bBGifdSjqe/EOrCzckRYphc7YsOpFy7l5otgKLMshnriPJeY8goR 9tx+Fkv9XCmWO++JCM7tunup9zXe7H1SJULLo50jFtsi51nzSoS+V5px1F/T Svifmkn+4CXxI+8rQy2r/asIXhWNI659tEW9p0eEoCSPfdTz1zyLHxgQwcWk acTFP/V+Dw+L8FFyRQL1/yfJRGk= "]]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->{True, True}, AxesLabel->{None, None}, AxesOrigin->{0, 0}, DisplayFunction->Identity, Frame->{{False, False}, {False, False}}, FrameLabel->{{None, None}, {None, None}}, FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}}, GridLines->{None, None}, GridLinesStyle->Directive[ GrayLevel[0.5, 0.4]], Method->{ "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> AbsolutePointSize[6], "ScalingFunctions" -> None}, PlotRange->{{0, 40.}, {0., 20.415184756511625`}}, PlotRangeClipping->True, PlotRangePadding->{{ Scaled[0.02], Scaled[0.02]}, { Scaled[0.05], Scaled[0.05]}}, Ticks->{Automatic, Automatic}]], "Output", CellChangeTimes->{{3.682167300038101*^9, 3.682167372779068*^9}, { 3.6821674120960503`*^9, 3.68216742262507*^9}, {3.682167485853795*^9, 3.682167507924296*^9}, {3.682167547646144*^9, 3.6821675676049423`*^9}, 3.682167693030156*^9, 3.682214941424644*^9, 3.682219720694622*^9}] }, Open ]], Cell[BoxData[{ RowBox[{ RowBox[{"enhance", "[", "x_", "]"}], " ", ":=", " ", RowBox[{"If", "[", RowBox[{ RowBox[{"x", "\[NotEqual]", " ", "0"}], ",", " ", RowBox[{"350", "+", RowBox[{ RowBox[{"(", RowBox[{"x", "-", "350"}], ")"}], "/", "5"}]}], ",", " ", "0"}], "]"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"basecolor", " ", "=", " ", "100"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Tau]start", " ", "=", "0"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"switch", "[", "t_", "]"}], " ", ":=", " ", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", "\[Tau]p", "}"}], ",", " ", RowBox[{ RowBox[{"\[Tau]p", " ", "=", " ", RowBox[{ RowBox[{"\[Tau]", "[", "t", "]"}], "/.", RowBox[{"sol", "[", RowBox[{"[", "1", "]"}], "]"}]}]}], ";", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"EvenQ", "[", RowBox[{"Floor", "[", RowBox[{"20", RowBox[{"(", RowBox[{"\[Tau]p", " ", "-", " ", "\[Tau]start"}], ")"}]}], "]"}], "]"}], " ", "||", " ", RowBox[{"(", RowBox[{ RowBox[{"enhance", "[", RowBox[{"color", "[", "t", "]"}], "]"}], " ", ">", " ", "1000"}], ")"}]}], ",", " ", "False", ",", " ", "True"}], "]"}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"switchalwayson", "[", "t_", "]"}], " ", ":=", " ", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", "\[Tau]p", "}"}], ",", " ", RowBox[{ RowBox[{"\[Tau]p", " ", "=", " ", RowBox[{ RowBox[{"\[Tau]", "[", "t", "]"}], "/.", RowBox[{"sol", "[", RowBox[{"[", "1", "]"}], "]"}]}]}], ";", RowBox[{"If", "[", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"enhance", "[", RowBox[{"color", "[", "t", "]"}], "]"}], " ", ">", " ", "1000"}], ")"}], ",", " ", "False", ",", " ", "True"}], "]"}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"color", "[", "t_", "]"}], " ", ":=", " ", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{"\[Tau]p", ",", " ", "rp", ",", "y"}], "}"}], ",", " ", RowBox[{ RowBox[{"\[Tau]p", " ", "=", " ", RowBox[{ RowBox[{"\[Tau]", "[", "t", "]"}], "/.", RowBox[{"sol", "[", RowBox[{"[", "1", "]"}], "]"}]}]}], ";", RowBox[{"rp", " ", "=", " ", RowBox[{ RowBox[{"x", "[", "t", "]"}], "/.", " ", RowBox[{"sol", "[", RowBox[{"[", "1", "]"}], "]"}]}]}], ";", " ", RowBox[{"basecolor", "*", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"(", RowBox[{"1", " ", "-", " ", RowBox[{"1", "/", "a"}]}], ")"}], "/", RowBox[{"(", RowBox[{"1", " ", "-", " ", RowBox[{"1", "/", "rp"}]}], ")"}]}], "]"}]}]}]}], "]"}]}]}], "Input",\ CellChangeTimes->CompressedData[" 1:eJxTTMoPSmViYGAQBmIQveyYXYZu/mtHgbJX+SD6ZoxXKYh2+vC7EkQ/l9Sv A9EzbsrWg+gfhwPbQLTdhttguq303AoQLfP/7kqwuleS60C0nq7nehCdXCe/ GUT/W2G/FUTztXDsB9EWC88fBNH3RNwugOjzNp5g+ty7p9dBtIuh3Q2w+Rk7 H4HoAxz9j8HmvpiQpQekC5+KFoBoq3OPHPWBtPbUBCcQ3bODxx9EX9qpEgyi Kw7OCAfz384G02Y+/dEg+vWVzWCaLepYhwvI/AufwTTTX60+EC2vcGICiJY4 bP3VDUg7/IgC0wDagJaq "]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Animate", "[", RowBox[{ RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{ RowBox[{"Polygon", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", RowBox[{"-", "1"}]}], "}"}], ",", " ", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", " ", "1"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1", ",", " ", "1"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1", ",", " ", RowBox[{"-", "1"}]}], "}"}]}], "}"}], "]"}], ",", RowBox[{ RowBox[{"ColorData", "[", "\"\\"", "]"}], "[", RowBox[{"enhance", "[", RowBox[{"color", "[", "x", "]"}], "]"}], "]"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{"switchalwayson", "[", "x", "]"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", " ", "0.1"}], "]"}], ",", RowBox[{"{", "}"}]}], "]"}]}], "}"}], "]"}], ",", " ", RowBox[{"{", RowBox[{"x", ",", " ", "23", ",", " ", "50"}], "}"}], ",", " ", RowBox[{"AnimationRepetitions", " ", "\[Rule]", " ", "1"}], ",", " ", RowBox[{"DefaultDuration", " ", "\[Rule]", " ", "30"}]}], "]"}]], "Input", CellChangeTimes->{{3.682167836371901*^9, 3.682167837944173*^9}, { 3.6821684056634912`*^9, 3.682168410426544*^9}, {3.682168453697583*^9, 3.682168454001737*^9}, {3.682168626637147*^9, 3.6821686291199017`*^9}, { 3.682168667661626*^9, 3.6821686684200497`*^9}, {3.6821687220398397`*^9, 3.682168723436467*^9}, {3.682169770518334*^9, 3.682169786926651*^9}, { 3.6821698578185997`*^9, 3.6821698603284483`*^9}, {3.682171446331304*^9, 3.682171471549302*^9}, {3.682215051805821*^9, 3.682215052499699*^9}, { 3.682215134097435*^9, 3.682215135330389*^9}, {3.682219872942745*^9, 3.682219873755604*^9}, {3.682219998107006*^9, 3.682219999885518*^9}}], Cell[BoxData[ TagBox[ StyleBox[ DynamicModuleBox[{$CellContext`x$$ = 50., Typeset`show$$ = True, Typeset`bookmarkList$$ = {}, Typeset`bookmarkMode$$ = "Menu", Typeset`animator$$, Typeset`animvar$$ = 1, Typeset`name$$ = "\"untitled\"", Typeset`specs$$ = {{ Hold[$CellContext`x$$], 23, 50}}, Typeset`size$$ = {360., {177., 182.}}, Typeset`update$$ = 0, Typeset`initDone$$, Typeset`skipInitDone$$ = True, $CellContext`x$310913$$ = 0}, DynamicBox[Manipulate`ManipulateBoxes[ 1, StandardForm, "Variables" :> {$CellContext`x$$ = 23}, "ControllerVariables" :> { Hold[$CellContext`x$$, $CellContext`x$310913$$, 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" :> Graphics[{ Polygon[{{-1, -1}, {-1, 1}, {1, 1}, {1, -1}}], ColorData["VisibleSpectrum"][ $CellContext`enhance[ $CellContext`color[$CellContext`x$$]]], If[ $CellContext`switchalwayson[$CellContext`x$$], Disk[{0, 0}, 0.1], {}]}], "Specifications" :> {{$CellContext`x$$, 23, 50, AnimationRepetitions -> 1, DefaultDuration -> 30, AppearanceElements -> { "ProgressSlider", "PlayPauseButton", "FasterSlowerButtons", "DirectionButton"}}}, "Options" :> { ControlType -> Animator, AppearanceElements -> None, DefaultBaseStyle -> "Animate", DefaultLabelStyle -> "AnimateLabel", SynchronousUpdating -> True, ShrinkingDelay -> 60.}, "DefaultOptions" :> {}], ImageSizeCache->{407., {212., 219.}}, 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.6821678277645683`*^9, 3.682167888318488*^9, 3.6821679738535357`*^9, { 3.682168337703006*^9, 3.6821683453064528`*^9}, 3.6821684107988462`*^9, 3.682168454311839*^9, 3.682168515921365*^9, 3.682168629566456*^9, { 3.682168668867856*^9, 3.682168680292962*^9}, 3.6821687242785807`*^9, 3.682168808239731*^9, 3.682169600726459*^9, 3.682169733284318*^9, 3.682169787802484*^9, 3.68216982331601*^9, 3.6821698607071943`*^9, 3.682171413483509*^9, {3.6821714489463797`*^9, 3.682171471986553*^9}, 3.682171546985385*^9, 3.6821715832368517`*^9, 3.682171615698099*^9, 3.682214944073044*^9, 3.682214986050104*^9, 3.682215052925708*^9, 3.6822169417374163`*^9, {3.6822198121474257`*^9, 3.6822198383696833`*^9}, 3.6822198741030483`*^9, {3.682219931636778*^9, 3.682219954016424*^9}, 3.6822200014695663`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Plot", "[", RowBox[{ RowBox[{"Sqrt", "[", RowBox[{"1", " ", "-", " ", RowBox[{"1", "/", "r"}]}], "]"}], ",", " ", RowBox[{"{", RowBox[{"r", ",", " ", "10", ",", " ", "1.1"}], "}"}]}], "]"}]], "Input", CellChangeTimes->{{3.682171848415207*^9, 3.682171856975987*^9}}], Cell[BoxData[ GraphicsBox[{{}, {}, {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[ 1.], LineBox[CompressedData[" 1:eJwVzns41Ikex/ExjLvcL2EMfohWqbVKtfp+k00r2omRS63cuyCKikbihHDc jlQuiWRDu22SLcmxkk2UEF2ObYsxZgYzv31ozTKM2Tl/fJ7P8/rvbRUe7xtF pVAox+X7/884lD+97LAANctV7vRyDlB5bY6r3i3Aua6IV8tGE5Ca4pd+iLMA fYHZN99aTIBYUzjcJFoAjdf2NrfsJmDamc72U1qEoLGOdLtNEzCUcb63bOMi aCR9oN1lTUAN3SOayF+EVL7KodmiCQD//mrXHRIQMyyJTtkE2HPH6fu8JeA6 KNjMUuaCTqK48miABMpTGdfGNLnAKWaUVcRJ4KfuxnzOai5kvzxRvFwugTtj 6RX5zlx4udMoo3NOAnaXWL5ZUVwI/vJQ+O66JfDTe/+vsG4u7HySOBZ6dwnW 3wiaMunjgiMzJySlbQla8isd+wa4sBJ3L/j24BJU9+9KNP2dCzcaFf00pUtw klSySpzlgsCyYeeg/zKksVTvaNEn4ZT2rE2gihTKjn0XaZ8wCbWxXRFh+lJg vve7V3FqEl71Xqo9xpBCZ5iyhQp7EuwzN1me2yyF4paItFdZk/BhgW1We1gK LU5+cSaVk+DBUdYT9kjl/bN7bnRPgsEv5rLzOStwQZo7363Pg6B1dcJhJQo6 pF3RUKnjwVxd2+MT6hQ8KjhAXarnQZ75UJ62DgUDj7VlzfzIg3aNFfs9ZhTk kRKPrvs8sJgOiOraSMFTqY12nk95wL2l8bEphIKTtI8RzeM8OG6ZOFDQSsEz EqM8R3M+XNB3b/aMVcC6ymz77Bw+bJ97MP79SQXs9wg3oebzYWHwC92kZAX0 rg83Zhfx4XihYUJNpgIuVZn7Rl7hQ6Da1LrFawrIU1a2Majjw3pZUePtfgX0 l+UlCDv4MDr9R42WExUz+g3Nyub4sLErpWh4loqO/1tqzfcVQN11R1b/31Q0 vcViB/kLwIj9yaRHSkV+W6q9TaAAJM7f1LapKWJz8UPv+98LoOsHnV9uWCni gEvSSusRAbBy60fj9ykix9NRNTJNAGeYb2w17ykibe6Y+u4GAfz344Z2jwQl 1Iz12dEqFsDItuARhTc0lK76d/2l5CmID6MrvQtQQY1lr6rM2Sm43OPLWOhU Rc26T8/P+E1DRaVK3oq1Oiqwmnr1GqbB2nW/Jv2mBiruHhN5L06Di6Hx2Ytq Wvg69LuAvq9n4GXbCZla4Src+3znDv+0GbDk20630nTw7EWHNVatcut695U8 1MHQmBIGbWoGwhnwIDlZF5PTSpm+q4Uw7TNzNdNRD4/EW97OcRfC/dDIrGt/ 6uFT7cfM9BghxKytMl1zVx+XGq1dE4qFcNSg6Id3YQZoN3L9G58mITRyfL4u JQzx7JzXMGtQCE5Dvz0//d4Q43pprRZCIdSdli6/yDXC9G8zE0ppIgimUN3W eRnj1AfZZBRDBPtv/3qyRWqMLPXwgnkXEZS6erV4dprg3LNNg5PeInhyrrTK 5/xqnFD5IlI5XARfTVxNTdloil/1HPitLEkEAscjh+/NmSKpViCuuCgCzpVZ T+16M6RbRV89VCaCFw+3Dn4KMEe6xuOVv+tFIG5nRrYb0/F1n8CK1iqCwxwG 2dFHx4bOkf1uz0SwWvj4Qn++BaogfT13RAS6RrraVHcG1ku4voJxESyIbd0b PzNQ51H5XSNSBDofL+D6YEs0l4n3NUtE0MMcGjVtscTQg7EykkZCtQuNXWBo hVsfdRws0SWBnDeJ33LKCoV/LJCp5iSE7lILSX5lhblrNrVstyPhs+qYr5uT Nd70tEiO2EBC44FL7pdzrfHFhkfMlC0kvM21y7bjWON2qVW+ujsJNmtrBtzc CfyyojpU4kXC5WCxvqCSwPOckAWZLwlMy8LE5CoCx+ndFoQfCZo822HVagKP ro2GXXJnnvAvcagl0Cw8Oj5f7qScFp2YBgIzxtLyTFgksB6c1CJbCOx23fOz kz8JBnp/Kv/VT6Czj96DoAASBt9lR2cOEPhspscuVe78KotnBkMEPn/D+891 uZUc9ma5jBC4RRoVMCH3X3BH8cwogZrlRdWxgSQMx8VQJDwCDxd02KYFkVDo rBiaJyCQu+1Htxq5vRYrfjWdJvB0jxGzS+6uzL60bSIC1V7ERyoHk9Bc6SBN /Szv79y7p1Du42FPDmrNE8i+ErC5Se61a4Laq8QE2kxN0V/LXducw+5YJLDf q2nU8AAJIcmWv+9dIjBlmde0WW7T7a1bPy3L++N/ygiS+60isyJ+hUDqvIIP W+6SXv6iTEagwgeBXpXc/wDxIEd6 "]]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->{True, True}, AxesLabel->{None, None}, AxesOrigin->{2., 0.5700000000000001}, DisplayFunction->Identity, Frame->{{False, False}, {False, False}}, FrameLabel->{{None, None}, {None, None}}, FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}}, GridLines->{None, None}, GridLinesStyle->Directive[ GrayLevel[0.5, 0.4]], Method->{ "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> AbsolutePointSize[6], "ScalingFunctions" -> None}, PlotRange->{{1.1, 10}, {0.5811907101741731, 0.9486832970932256}}, PlotRangeClipping->True, PlotRangePadding->{{ Scaled[0.02], Scaled[0.02]}, { Scaled[0.05], Scaled[0.05]}}, Ticks->{Automatic, Automatic}]], "Output", CellChangeTimes->{3.682171857297123*^9}] }, Open ]] }, WindowSize->{1920, 1033}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, FrontEndVersion->"10.1 for Linux x86 (64-bit) (March 23, 2015)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[558, 20, 1545, 21, 32, "Input"], Cell[CellGroupData[{ Cell[2128, 45, 1537, 39, 77, "Input"], Cell[3668, 86, 314, 4, 32, "Output"], Cell[3985, 92, 22895, 437, 125, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[26917, 534, 643, 15, 32, "Input"], Cell[27563, 551, 16234, 278, 238, "Output"] }, Open ]], Cell[43812, 832, 3380, 96, 165, "Input"], Cell[CellGroupData[{ Cell[47217, 932, 2004, 45, 55, "Input"], Cell[49224, 979, 3030, 58, 448, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[52291, 1042, 316, 8, 32, "Input"], Cell[52610, 1052, 3338, 66, 250, "Output"] }, Open ]] } ] *) (* End of internal cache information *)