(* 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[ 50843, 1065] NotebookOptionsPosition[ 49438, 1012] NotebookOutlinePosition[ 49793, 1028] CellTagsIndexPosition[ 49750, 1025] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell["Solving the differential equation for problem 1", "Subsubsection", CellChangeTimes->{{3.720050704220381*^9, 3.720050717000328*^9}}], Cell[BoxData[ RowBox[{"Clear", "[", RowBox[{ "xmin", ",", " ", "xmax", ",", " ", "xmid", ",", " ", "psol1", ",", " ", "psol2", ",", " ", "y1", ",", " ", "y2", ",", " ", "dy1", ",", " ", "dy2", ",", " ", "W"}], "]"}]], "Input", CellChangeTimes->{{3.7201059559018793`*^9, 3.7201059857524223`*^9}, { 3.7201066406677523`*^9, 3.720106651090966*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"xmin", "=", " ", "0.0001"}], "\[IndentingNewLine]", RowBox[{"xmax", " ", "=", " ", "20"}], "\[IndentingNewLine]", RowBox[{"xmid", " ", "=", " ", "3"}]}], "Input", CellChangeTimes->{{3.720052641014125*^9, 3.720052675764423*^9}, { 3.720054916905231*^9, 3.7200549189350357`*^9}, {3.72005497371109*^9, 3.720054974238084*^9}, {3.7200550607880697`*^9, 3.72005506167378*^9}, { 3.720095654622414*^9, 3.720095677767705*^9}, 3.720098595536298*^9}], Cell[BoxData["0.0001`"], "Output", CellChangeTimes->{3.7201062118640757`*^9, 3.720106748431492*^9}], Cell[BoxData["20"], "Output", CellChangeTimes->{3.7201062118640757`*^9, 3.720106748446653*^9}], Cell[BoxData["3"], "Output", CellChangeTimes->{3.7201062118640757`*^9, 3.720106748457711*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Define the weight", "Subsubsection", CellChangeTimes->{{3.72011693271446*^9, 3.7201169362857437`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{ RowBox[{"w", "[", "x_", "]"}], " ", "=", " ", RowBox[{ RowBox[{"x", "^", "3"}], "/", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"x", "^", "2"}], " ", "+", " ", "1"}], ")"}], "^", "2"}]}]}]}]], "Input", CellChangeTimes->{{3.720116940770988*^9, 3.720116973523247*^9}}], Cell[BoxData[ FractionBox[ SuperscriptBox["x", "3"], SuperscriptBox[ RowBox[{"(", RowBox[{"1", "+", SuperscriptBox["x", "2"]}], ")"}], "2"]]], "Output", CellChangeTimes->{3.720116975345306*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Construct the solution for small x", "Subsubsection", CellChangeTimes->{{3.720107847920246*^9, 3.720107856024756*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"y1small", "[", RowBox[{"e_", ",", " ", RowBox[{"x_", "/;", " ", RowBox[{ "0", "\[LessEqual]", " ", "x", " ", "\[LessEqual]", " ", "1"}]}]}], "]"}], " ", ":=", " ", RowBox[{"1", "-", " ", RowBox[{ RowBox[{"e", "/", "8"}], " ", RowBox[{"x", "^", "2"}]}]}]}]], "Input", CellChangeTimes->{{3.720107253660522*^9, 3.7201073662203712`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"dy1small", "[", RowBox[{"e_", ",", " ", RowBox[{"x_", "/;", " ", RowBox[{ "0", " ", "\[LessEqual]", " ", "x", " ", "\[LessEqual]", " ", "1"}]}]}], "]"}], " ", ":=", " ", RowBox[{ RowBox[{ RowBox[{"-", "e"}], "/", "4"}], " ", "x"}]}]], "Input", CellChangeTimes->{{3.720107484133718*^9, 3.720107533928644*^9}, { 3.720107568631857*^9, 3.720107575039998*^9}, {3.720107625377116*^9, 3.720107669947968*^9}}] }, Open ]], Cell[CellGroupData[{ Cell["Construct the solution for x up to xmid", "Subsubsection", CellChangeTimes->{{3.720107870319697*^9, 3.720107904078603*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"psol1", " ", "=", RowBox[{"ParametricNDSolve", "[", RowBox[{ RowBox[{"{", " ", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"D", "[", " ", RowBox[{ RowBox[{ RowBox[{"x", "^", "3"}], " ", RowBox[{"D", "[", RowBox[{ RowBox[{"y", "[", "x", "]"}], ",", "x"}], "]"}]}], ",", "x"}], "]"}]}], " ", "\[Equal]", " ", RowBox[{"ee", " ", RowBox[{"w", "[", "x", "]"}], " ", RowBox[{"y", "[", "x", "]"}]}]}], ",", RowBox[{ RowBox[{"y", "[", "xmin", "]"}], "\[Equal]", RowBox[{"y1small", "[", RowBox[{"ee", ",", "xmin"}], "]"}]}], ",", RowBox[{ RowBox[{ RowBox[{"y", "'"}], "[", "xmin", "]"}], "\[Equal]", RowBox[{"dy1small", "[", RowBox[{"ee", ",", "xmin"}], "]"}]}]}], "}"}], ",", "y", ",", RowBox[{"{", RowBox[{"x", ",", "xmin", ",", " ", "xmid"}], "}"}], ",", " ", RowBox[{"{", "ee", "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.7200952488082113`*^9, 3.7200952837933073`*^9}, { 3.7200956412868834`*^9, 3.720095644538021*^9}, {3.720095686586073*^9, 3.720095702664434*^9}, 3.720095783292328*^9, {3.7200989136019382`*^9, 3.720098923663498*^9}, 3.72009903852103*^9, {3.720107472100485*^9, 3.7201074809706583`*^9}, {3.720107541353344*^9, 3.720107558241611*^9}, { 3.7201169795977297`*^9, 3.7201170008951607`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"y", "\[Rule]", TagBox[ TemplateBox[{RowBox[{ StyleBox[ TagBox["ParametricFunction", "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[{{}, {}, { GrayLevel[0.55], AbsoluteThickness[1.5], Opacity[1.], LineBox[{{9.09090909090909*^-8, 9.090909090884856*^-8}, { 0.0874367476365131, 0.08197166587636243}, { 0.18222810297558026`, 0.1358255024317191}, { 0.27073779576926765`, 0.13315455086535802`}, { 0.3575112797365835, 0.08143316377978886}, { 0.4516394623155443, -0.007347854733426912}, { 0.5394859823491253, -0.08816932746972758}, { 0.6346872009943513, -0.13787329990112937`}, { 0.7281522108132057, -0.12816574147938434`}, { 0.8153355580866803, -0.07001063326902106}, { 0.9098736039718, 0.02121659675083076}, { 0.9999999090909091, 0.10024804094746914`}}]}, { GrayLevel[0.55], AbsoluteThickness[1.5], Opacity[1.], LineBox[{{9.09090909090909*^-8, 9.090909090880112*^-8}, { 0.0874367476365131, 0.08090369567458766}, { 0.18222810297558026`, 0.12746559200130417`}, { 0.27073779576926765`, 0.1116266937044405}, { 0.3575112797365835, 0.046963725126600256`}, { 0.4516394623155443, -0.04509050859182233}, { 0.5394859823491253, -0.11115559892361665`}, { 0.6346872009943513, -0.12642564594664163`}, { 0.7281522108132057, -0.07740196037964171}, { 0.8153355580866803, 0.004179083230780074}, { 0.9098736039718, 0.08937493944143977}, { 0.9999999090909091, 0.12834702174618903`}}]}, { GrayLevel[0.55], AbsoluteThickness[1.5], Opacity[1.], LineBox[{{9.09090909090909*^-8, 9.090909090875369*^-8}, { 0.0874367476365131, 0.07984396759533968}, { 0.18222810297558026`, 0.11939549826670079`}, { 0.27073779576926765`, 0.09184904274707059}, { 0.3575112797365835, 0.017913117173780694`}, { 0.4516394623155443, -0.07109962300031061}, { 0.5394859823491253, -0.11717052830989395`}, { 0.6346872009943513, -0.09879571259704382}, { 0.7281522108132057, -0.022693471356141894`}, { 0.8153355580866803, 0.06128703302641453}, { 0.9098736039718, 0.11605625481609519`}, { 0.9999999090909091, 0.10415981267620744`}}]}, { GrayLevel[0.55], AbsoluteThickness[1.5], Opacity[1.], LineBox[{{9.09090909090909*^-8, 9.090909090870287*^-8}, { 0.0874367476365131, 0.07879243663407876}, { 0.18222810297558026`, 0.1116081822210312}, { 0.27073779576926765`, 0.07372482040530459}, { 0.3575112797365835, -0.006267252216927747}, { 0.4516394623155443, -0.08743897872291118}, { 0.5394859823491253, -0.11109134853087661`}, { 0.6346872009943513, -0.06371549528296311}, { 0.7281522108132057, 0.02544576624054603}, { 0.8153355580866803, 0.09464282937855747}, { 0.9098736039718, 0.10731931496853601`}, { 0.9999999090909091, 0.0516787429232188}}]}, { GrayLevel[0.55], AbsoluteThickness[1.5], Opacity[1.], LineBox[{{9.09090909090909*^-8, 9.090909090865543*^-8}, { 0.0874367476365131, 0.07774905743802465}, { 0.18222810297558026`, 0.10409673959866513`}, { 0.27073779576926765`, 0.05716133039337165}, { 0.3575112797365835, -0.026084496726167965`}, { 0.4516394623155443, -0.09591154799603378}, { 0.5394859823491253, -0.09685594593277226}, { 0.6346872009943513, -0.02730376120671303}, { 0.7281522108132057, 0.061646715198529656`}, { 0.8153355580866803, 0.1046634148537508}, {0.9098736039718, 0.07509702498429204}, { 0.9999999090909091, -0.006536873471333553}}]}, { GrayLevel[0.55], AbsoluteThickness[1.5], Opacity[1.], LineBox[{{9.09090909090909*^-8, 9.0909090908608*^-8}, { 0.0874367476365131, 0.07671378556951527}, { 0.18222810297558026`, 0.09685439581479163}, { 0.27073779576926765`, 0.042069810325321415`}, { 0.3575112797365835, -0.042006785469922474`}, { 0.4516394623155443, -0.09808544986431629}, { 0.5394859823491253, -0.07760167026270944}, { 0.6346872009943513, 0.0063643674428021085`}, { 0.7281522108132057, 0.08405712128907425}, { 0.8153355580866803, 0.09555213349985489}, { 0.9098736039718, 0.03202960165577907}, { 0.9999999090909091, -0.054402034659985464`}}]}}, AspectRatio -> 1, Axes -> {False, False}, AxesLabel -> {None, None}, AxesOrigin -> {0, 0}, Background -> GrayLevel[0.93], BaseStyle -> {FontFamily -> "Arial"}, DisplayFunction -> Identity, Frame -> {{True, True}, {True, True}}, FrameLabel -> {{None, None}, {None, None}}, FrameStyle -> Directive[ Thickness[Tiny], GrayLevel[0.7]], FrameTicks -> {{None, None}, {None, None}}, GridLines -> {None, None}, GridLinesStyle -> Directive[ GrayLevel[0.5, 0.4]], ImageSize -> {Automatic, Dynamic[ 3.5 (CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[ Magnification])]}, LabelStyle -> {FontFamily -> "Arial"}, Method -> {"ScalingFunctions" -> None}, PlotRange -> {All, All}, PlotRangeClipping -> True, PlotRangePadding -> {Automatic, Automatic}, Ticks -> {Automatic, Automatic}], GridBox[{{ RowBox[{ TagBox["\"Expression: \"", "SummaryItemAnnotation"], "\[InvisibleSpace]", TagBox[ TagBox["y", HoldForm], "SummaryItem"]}]}, { RowBox[{ TagBox["\"Parameters: \"", "SummaryItemAnnotation"], "\[InvisibleSpace]", TagBox[ RowBox[{"{", "ee", "}"}], "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[{{}, {}, { GrayLevel[0.55], AbsoluteThickness[1.5], Opacity[1.], LineBox[{{9.09090909090909*^-8, 9.090909090884856*^-8}, { 0.0874367476365131, 0.08197166587636243}, { 0.18222810297558026`, 0.1358255024317191}, { 0.27073779576926765`, 0.13315455086535802`}, { 0.3575112797365835, 0.08143316377978886}, { 0.4516394623155443, -0.007347854733426912}, { 0.5394859823491253, -0.08816932746972758}, { 0.6346872009943513, -0.13787329990112937`}, { 0.7281522108132057, -0.12816574147938434`}, { 0.8153355580866803, -0.07001063326902106}, { 0.9098736039718, 0.02121659675083076}, { 0.9999999090909091, 0.10024804094746914`}}]}, { GrayLevel[0.55], AbsoluteThickness[1.5], Opacity[1.], LineBox[{{9.09090909090909*^-8, 9.090909090880112*^-8}, { 0.0874367476365131, 0.08090369567458766}, { 0.18222810297558026`, 0.12746559200130417`}, { 0.27073779576926765`, 0.1116266937044405}, { 0.3575112797365835, 0.046963725126600256`}, { 0.4516394623155443, -0.04509050859182233}, { 0.5394859823491253, -0.11115559892361665`}, { 0.6346872009943513, -0.12642564594664163`}, { 0.7281522108132057, -0.07740196037964171}, { 0.8153355580866803, 0.004179083230780074}, { 0.9098736039718, 0.08937493944143977}, { 0.9999999090909091, 0.12834702174618903`}}]}, { GrayLevel[0.55], AbsoluteThickness[1.5], Opacity[1.], LineBox[{{9.09090909090909*^-8, 9.090909090875369*^-8}, { 0.0874367476365131, 0.07984396759533968}, { 0.18222810297558026`, 0.11939549826670079`}, { 0.27073779576926765`, 0.09184904274707059}, { 0.3575112797365835, 0.017913117173780694`}, { 0.4516394623155443, -0.07109962300031061}, { 0.5394859823491253, -0.11717052830989395`}, { 0.6346872009943513, -0.09879571259704382}, { 0.7281522108132057, -0.022693471356141894`}, { 0.8153355580866803, 0.06128703302641453}, { 0.9098736039718, 0.11605625481609519`}, { 0.9999999090909091, 0.10415981267620744`}}]}, { GrayLevel[0.55], AbsoluteThickness[1.5], Opacity[1.], LineBox[{{9.09090909090909*^-8, 9.090909090870287*^-8}, { 0.0874367476365131, 0.07879243663407876}, { 0.18222810297558026`, 0.1116081822210312}, { 0.27073779576926765`, 0.07372482040530459}, { 0.3575112797365835, -0.006267252216927747}, { 0.4516394623155443, -0.08743897872291118}, { 0.5394859823491253, -0.11109134853087661`}, { 0.6346872009943513, -0.06371549528296311}, { 0.7281522108132057, 0.02544576624054603}, { 0.8153355580866803, 0.09464282937855747}, { 0.9098736039718, 0.10731931496853601`}, { 0.9999999090909091, 0.0516787429232188}}]}, { GrayLevel[0.55], AbsoluteThickness[1.5], Opacity[1.], LineBox[{{9.09090909090909*^-8, 9.090909090865543*^-8}, { 0.0874367476365131, 0.07774905743802465}, { 0.18222810297558026`, 0.10409673959866513`}, { 0.27073779576926765`, 0.05716133039337165}, { 0.3575112797365835, -0.026084496726167965`}, { 0.4516394623155443, -0.09591154799603378}, { 0.5394859823491253, -0.09685594593277226}, { 0.6346872009943513, -0.02730376120671303}, { 0.7281522108132057, 0.061646715198529656`}, { 0.8153355580866803, 0.1046634148537508}, {0.9098736039718, 0.07509702498429204}, { 0.9999999090909091, -0.006536873471333553}}]}, { GrayLevel[0.55], AbsoluteThickness[1.5], Opacity[1.], LineBox[{{9.09090909090909*^-8, 9.0909090908608*^-8}, { 0.0874367476365131, 0.07671378556951527}, { 0.18222810297558026`, 0.09685439581479163}, { 0.27073779576926765`, 0.042069810325321415`}, { 0.3575112797365835, -0.042006785469922474`}, { 0.4516394623155443, -0.09808544986431629}, { 0.5394859823491253, -0.07760167026270944}, { 0.6346872009943513, 0.0063643674428021085`}, { 0.7281522108132057, 0.08405712128907425}, { 0.8153355580866803, 0.09555213349985489}, { 0.9098736039718, 0.03202960165577907}, { 0.9999999090909091, -0.054402034659985464`}}]}}, AspectRatio -> 1, Axes -> {False, False}, AxesLabel -> {None, None}, AxesOrigin -> {0, 0}, Background -> GrayLevel[0.93], BaseStyle -> {FontFamily -> "Arial"}, DisplayFunction -> Identity, Frame -> {{True, True}, {True, True}}, FrameLabel -> {{None, None}, {None, None}}, FrameStyle -> Directive[ Thickness[Tiny], GrayLevel[0.7]], FrameTicks -> {{None, None}, {None, None}}, GridLines -> {None, None}, GridLinesStyle -> Directive[ GrayLevel[0.5, 0.4]], ImageSize -> {Automatic, Dynamic[ 3.5 (CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[ Magnification])]}, LabelStyle -> {FontFamily -> "Arial"}, Method -> {"ScalingFunctions" -> None}, PlotRange -> {All, All}, PlotRangeClipping -> True, PlotRangePadding -> {Automatic, Automatic}, Ticks -> {Automatic, Automatic}], GridBox[{{ RowBox[{ TagBox["\"Expression: \"", "SummaryItemAnnotation"], "\[InvisibleSpace]", TagBox[ TagBox["y", HoldForm], "SummaryItem"]}]}, { RowBox[{ TagBox["\"Parameters: \"", "SummaryItemAnnotation"], "\[InvisibleSpace]", TagBox[ RowBox[{"{", "ee", "}"}], "SummaryItem"]}]}, { RowBox[{ TagBox["\"Generator: \"", "SummaryItemAnnotation"], "\[InvisibleSpace]", TagBox["ParametricNDSolve", "SummaryItem"]}]}, { RowBox[{ TagBox[ "\"Dependent variables: \"", "SummaryItemAnnotation"], "\[InvisibleSpace]", TagBox[ RowBox[{"{", "y", "}"}], "SummaryItem"]}]}, { RowBox[{ TagBox[ "\"Independent variables: \"", "SummaryItemAnnotation"], "\[InvisibleSpace]", TagBox[ RowBox[{"{", "x", "}"}], "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->("ParametricFunction[<>]"& )], False, Editable->False, SelectWithContents->True, Selectable->False]}], "}"}]], "Output", CellChangeTimes->{ 3.720099039881722*^9, 3.7200996109196663`*^9, 3.720101145616954*^9, 3.720106215749918*^9, 3.720106319874859*^9, {3.72010674398343*^9, 3.720106751190228*^9}, 3.720107585818474*^9, 3.720107630054943*^9, 3.720107673464197*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"y1", "[", RowBox[{ RowBox[{"e_", "?", "NumericQ"}], ",", " ", RowBox[{"x_", "/;", " ", RowBox[{ "xmin", " ", "\[LessEqual]", " ", "x", " ", "\[LessEqual]", " ", "xmid"}]}]}], "]"}], " ", ":=", " ", RowBox[{"Evaluate", "[", " ", RowBox[{ RowBox[{ RowBox[{"y", "[", "e", "]"}], "[", "x", "]"}], " ", "/.", " ", "psol1"}], "]"}]}]], "Input", CellChangeTimes->{{3.720097165595976*^9, 3.720097168471047*^9}, { 3.7200989261482964`*^9, 3.720098931456727*^9}, {3.7200991616252604`*^9, 3.720099168206476*^9}, {3.720099616009536*^9, 3.720099617023347*^9}, { 3.720101087722439*^9, 3.7201010911480494`*^9}, {3.720102195350438*^9, 3.7201022298547707`*^9}, {3.720106013385662*^9, 3.7201060543521976`*^9}, { 3.720106140083528*^9, 3.720106142441667*^9}, {3.720106665207788*^9, 3.720106665666418*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"dy1", "[", RowBox[{ RowBox[{"e_", "?", "NumericQ"}], ",", " ", RowBox[{"x_", "/;", " ", RowBox[{ "xmin", " ", "\[LessEqual]", " ", "x", " ", "\[LessEqual]", " ", "xmid"}]}]}], "]"}], " ", ":=", " ", RowBox[{"Evaluate", "[", " ", RowBox[{ RowBox[{ RowBox[{ RowBox[{"y", "[", "e", "]"}], "'"}], "[", "x", "]"}], "/.", " ", "psol1"}], "]"}]}]], "Input", CellChangeTimes->{{3.72009740913022*^9, 3.720097415971609*^9}, { 3.720097715356227*^9, 3.7200977165887623`*^9}, {3.72009893436171*^9, 3.7200989397125597`*^9}, {3.720099277252261*^9, 3.720099278477332*^9}, { 3.720099622787568*^9, 3.720099651991693*^9}, {3.720101093292836*^9, 3.7201011086500807`*^9}, {3.720102237588152*^9, 3.720102249141653*^9}, { 3.720106048175606*^9, 3.720106068609885*^9}, {3.7201061459302673`*^9, 3.720106149442648*^9}, {3.720106667705306*^9, 3.720106668226614*^9}, { 3.720106761298559*^9, 3.7201067618317757`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell["Construct the solution for x from xmid up to xmax", "Subsubsection", CellChangeTimes->{{3.720107887511181*^9, 3.720107911855659*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"psol2", " ", "=", " ", RowBox[{"ParametricNDSolve", "[", RowBox[{ RowBox[{"{", " ", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"D", "[", " ", RowBox[{ RowBox[{ RowBox[{"x", "^", "3"}], " ", RowBox[{"D", "[", RowBox[{ RowBox[{"y", "[", "x", "]"}], ",", "x"}], "]"}]}], ",", "x"}], "]"}]}], " ", "\[Equal]", " ", RowBox[{"ee", " ", RowBox[{"w", "[", "x", "]"}], " ", RowBox[{"y", "[", "x", "]"}]}]}], ",", RowBox[{ RowBox[{"y", "[", "xmax", "]"}], "\[Equal]", RowBox[{"1", "/", RowBox[{"xmax", "^", "2"}]}]}], ",", RowBox[{ RowBox[{ RowBox[{"y", "'"}], "[", "xmax", "]"}], "\[Equal]", RowBox[{ RowBox[{"-", "2"}], "/", RowBox[{"xmax", "^", "3"}]}]}]}], "}"}], ",", "y", ",", RowBox[{"{", RowBox[{"x", ",", "xmax", ",", " ", "xmid"}], "}"}], ",", " ", RowBox[{"{", "ee", "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.720052578461878*^9, 3.720052705806698*^9}, { 3.7200543324088078`*^9, 3.7200543789490757`*^9}, {3.72005468030716*^9, 3.7200547178462067`*^9}, {3.7200937107074947`*^9, 3.720093744108863*^9}, { 3.720093844023532*^9, 3.7200938478000603`*^9}, {3.720094217436778*^9, 3.720094245875267*^9}, {3.720094637756028*^9, 3.720094639253981*^9}, { 3.720095272663909*^9, 3.720095281857132*^9}, {3.720095670567429*^9, 3.72009567570378*^9}, {3.720096813044058*^9, 3.720096813185864*^9}, { 3.720098604314726*^9, 3.7200986062761393`*^9}, {3.720098882770576*^9, 3.720098882929544*^9}, {3.720117008846217*^9, 3.720117017658704*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"y", "\[Rule]", TagBox[ TemplateBox[{RowBox[{ StyleBox[ TagBox["ParametricFunction", "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[{{}, {}, { GrayLevel[0.55], AbsoluteThickness[1.5], Opacity[1.], LineBox[{{9.09090909090909*^-8, 9.090909090884856*^-8}, { 0.0874367476365131, 0.08197166587636243}, { 0.18222810297558026`, 0.1358255024317191}, { 0.27073779576926765`, 0.13315455086535802`}, { 0.3575112797365835, 0.08143316377978886}, { 0.4516394623155443, -0.007347854733426912}, { 0.5394859823491253, -0.08816932746972758}, { 0.6346872009943513, -0.13787329990112937`}, { 0.7281522108132057, -0.12816574147938434`}, { 0.8153355580866803, -0.07001063326902106}, { 0.9098736039718, 0.02121659675083076}, { 0.9999999090909091, 0.10024804094746914`}}]}, { GrayLevel[0.55], AbsoluteThickness[1.5], Opacity[1.], LineBox[{{9.09090909090909*^-8, 9.090909090880112*^-8}, { 0.0874367476365131, 0.08090369567458766}, { 0.18222810297558026`, 0.12746559200130417`}, { 0.27073779576926765`, 0.1116266937044405}, { 0.3575112797365835, 0.046963725126600256`}, { 0.4516394623155443, -0.04509050859182233}, { 0.5394859823491253, -0.11115559892361665`}, { 0.6346872009943513, -0.12642564594664163`}, { 0.7281522108132057, -0.07740196037964171}, { 0.8153355580866803, 0.004179083230780074}, { 0.9098736039718, 0.08937493944143977}, { 0.9999999090909091, 0.12834702174618903`}}]}, { GrayLevel[0.55], AbsoluteThickness[1.5], Opacity[1.], LineBox[{{9.09090909090909*^-8, 9.090909090875369*^-8}, { 0.0874367476365131, 0.07984396759533968}, { 0.18222810297558026`, 0.11939549826670079`}, { 0.27073779576926765`, 0.09184904274707059}, { 0.3575112797365835, 0.017913117173780694`}, { 0.4516394623155443, -0.07109962300031061}, { 0.5394859823491253, -0.11717052830989395`}, { 0.6346872009943513, -0.09879571259704382}, { 0.7281522108132057, -0.022693471356141894`}, { 0.8153355580866803, 0.06128703302641453}, { 0.9098736039718, 0.11605625481609519`}, { 0.9999999090909091, 0.10415981267620744`}}]}, { GrayLevel[0.55], AbsoluteThickness[1.5], Opacity[1.], LineBox[{{9.09090909090909*^-8, 9.090909090870287*^-8}, { 0.0874367476365131, 0.07879243663407876}, { 0.18222810297558026`, 0.1116081822210312}, { 0.27073779576926765`, 0.07372482040530459}, { 0.3575112797365835, -0.006267252216927747}, { 0.4516394623155443, -0.08743897872291118}, { 0.5394859823491253, -0.11109134853087661`}, { 0.6346872009943513, -0.06371549528296311}, { 0.7281522108132057, 0.02544576624054603}, { 0.8153355580866803, 0.09464282937855747}, { 0.9098736039718, 0.10731931496853601`}, { 0.9999999090909091, 0.0516787429232188}}]}, { GrayLevel[0.55], AbsoluteThickness[1.5], Opacity[1.], LineBox[{{9.09090909090909*^-8, 9.090909090865543*^-8}, { 0.0874367476365131, 0.07774905743802465}, { 0.18222810297558026`, 0.10409673959866513`}, { 0.27073779576926765`, 0.05716133039337165}, { 0.3575112797365835, -0.026084496726167965`}, { 0.4516394623155443, -0.09591154799603378}, { 0.5394859823491253, -0.09685594593277226}, { 0.6346872009943513, -0.02730376120671303}, { 0.7281522108132057, 0.061646715198529656`}, { 0.8153355580866803, 0.1046634148537508}, {0.9098736039718, 0.07509702498429204}, { 0.9999999090909091, -0.006536873471333553}}]}, { GrayLevel[0.55], AbsoluteThickness[1.5], Opacity[1.], LineBox[{{9.09090909090909*^-8, 9.0909090908608*^-8}, { 0.0874367476365131, 0.07671378556951527}, { 0.18222810297558026`, 0.09685439581479163}, { 0.27073779576926765`, 0.042069810325321415`}, { 0.3575112797365835, -0.042006785469922474`}, { 0.4516394623155443, -0.09808544986431629}, { 0.5394859823491253, -0.07760167026270944}, { 0.6346872009943513, 0.0063643674428021085`}, { 0.7281522108132057, 0.08405712128907425}, { 0.8153355580866803, 0.09555213349985489}, { 0.9098736039718, 0.03202960165577907}, { 0.9999999090909091, -0.054402034659985464`}}]}}, AspectRatio -> 1, Axes -> {False, False}, AxesLabel -> {None, None}, AxesOrigin -> {0, 0}, Background -> GrayLevel[0.93], BaseStyle -> {FontFamily -> "Arial"}, DisplayFunction -> Identity, Frame -> {{True, True}, {True, True}}, FrameLabel -> {{None, None}, {None, None}}, FrameStyle -> Directive[ Thickness[Tiny], GrayLevel[0.7]], FrameTicks -> {{None, None}, {None, None}}, GridLines -> {None, None}, GridLinesStyle -> Directive[ GrayLevel[0.5, 0.4]], ImageSize -> {Automatic, Dynamic[ 3.5 (CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[ Magnification])]}, LabelStyle -> {FontFamily -> "Arial"}, Method -> {"ScalingFunctions" -> None}, PlotRange -> {All, All}, PlotRangeClipping -> True, PlotRangePadding -> {Automatic, Automatic}, Ticks -> {Automatic, Automatic}], GridBox[{{ RowBox[{ TagBox["\"Expression: \"", "SummaryItemAnnotation"], "\[InvisibleSpace]", TagBox[ TagBox["y", HoldForm], "SummaryItem"]}]}, { RowBox[{ TagBox["\"Parameters: \"", "SummaryItemAnnotation"], "\[InvisibleSpace]", TagBox[ RowBox[{"{", "ee", "}"}], "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[{{}, {}, { GrayLevel[0.55], AbsoluteThickness[1.5], Opacity[1.], LineBox[{{9.09090909090909*^-8, 9.090909090884856*^-8}, { 0.0874367476365131, 0.08197166587636243}, { 0.18222810297558026`, 0.1358255024317191}, { 0.27073779576926765`, 0.13315455086535802`}, { 0.3575112797365835, 0.08143316377978886}, { 0.4516394623155443, -0.007347854733426912}, { 0.5394859823491253, -0.08816932746972758}, { 0.6346872009943513, -0.13787329990112937`}, { 0.7281522108132057, -0.12816574147938434`}, { 0.8153355580866803, -0.07001063326902106}, { 0.9098736039718, 0.02121659675083076}, { 0.9999999090909091, 0.10024804094746914`}}]}, { GrayLevel[0.55], AbsoluteThickness[1.5], Opacity[1.], LineBox[{{9.09090909090909*^-8, 9.090909090880112*^-8}, { 0.0874367476365131, 0.08090369567458766}, { 0.18222810297558026`, 0.12746559200130417`}, { 0.27073779576926765`, 0.1116266937044405}, { 0.3575112797365835, 0.046963725126600256`}, { 0.4516394623155443, -0.04509050859182233}, { 0.5394859823491253, -0.11115559892361665`}, { 0.6346872009943513, -0.12642564594664163`}, { 0.7281522108132057, -0.07740196037964171}, { 0.8153355580866803, 0.004179083230780074}, { 0.9098736039718, 0.08937493944143977}, { 0.9999999090909091, 0.12834702174618903`}}]}, { GrayLevel[0.55], AbsoluteThickness[1.5], Opacity[1.], LineBox[{{9.09090909090909*^-8, 9.090909090875369*^-8}, { 0.0874367476365131, 0.07984396759533968}, { 0.18222810297558026`, 0.11939549826670079`}, { 0.27073779576926765`, 0.09184904274707059}, { 0.3575112797365835, 0.017913117173780694`}, { 0.4516394623155443, -0.07109962300031061}, { 0.5394859823491253, -0.11717052830989395`}, { 0.6346872009943513, -0.09879571259704382}, { 0.7281522108132057, -0.022693471356141894`}, { 0.8153355580866803, 0.06128703302641453}, { 0.9098736039718, 0.11605625481609519`}, { 0.9999999090909091, 0.10415981267620744`}}]}, { GrayLevel[0.55], AbsoluteThickness[1.5], Opacity[1.], LineBox[{{9.09090909090909*^-8, 9.090909090870287*^-8}, { 0.0874367476365131, 0.07879243663407876}, { 0.18222810297558026`, 0.1116081822210312}, { 0.27073779576926765`, 0.07372482040530459}, { 0.3575112797365835, -0.006267252216927747}, { 0.4516394623155443, -0.08743897872291118}, { 0.5394859823491253, -0.11109134853087661`}, { 0.6346872009943513, -0.06371549528296311}, { 0.7281522108132057, 0.02544576624054603}, { 0.8153355580866803, 0.09464282937855747}, { 0.9098736039718, 0.10731931496853601`}, { 0.9999999090909091, 0.0516787429232188}}]}, { GrayLevel[0.55], AbsoluteThickness[1.5], Opacity[1.], LineBox[{{9.09090909090909*^-8, 9.090909090865543*^-8}, { 0.0874367476365131, 0.07774905743802465}, { 0.18222810297558026`, 0.10409673959866513`}, { 0.27073779576926765`, 0.05716133039337165}, { 0.3575112797365835, -0.026084496726167965`}, { 0.4516394623155443, -0.09591154799603378}, { 0.5394859823491253, -0.09685594593277226}, { 0.6346872009943513, -0.02730376120671303}, { 0.7281522108132057, 0.061646715198529656`}, { 0.8153355580866803, 0.1046634148537508}, {0.9098736039718, 0.07509702498429204}, { 0.9999999090909091, -0.006536873471333553}}]}, { GrayLevel[0.55], AbsoluteThickness[1.5], Opacity[1.], LineBox[{{9.09090909090909*^-8, 9.0909090908608*^-8}, { 0.0874367476365131, 0.07671378556951527}, { 0.18222810297558026`, 0.09685439581479163}, { 0.27073779576926765`, 0.042069810325321415`}, { 0.3575112797365835, -0.042006785469922474`}, { 0.4516394623155443, -0.09808544986431629}, { 0.5394859823491253, -0.07760167026270944}, { 0.6346872009943513, 0.0063643674428021085`}, { 0.7281522108132057, 0.08405712128907425}, { 0.8153355580866803, 0.09555213349985489}, { 0.9098736039718, 0.03202960165577907}, { 0.9999999090909091, -0.054402034659985464`}}]}}, AspectRatio -> 1, Axes -> {False, False}, AxesLabel -> {None, None}, AxesOrigin -> {0, 0}, Background -> GrayLevel[0.93], BaseStyle -> {FontFamily -> "Arial"}, DisplayFunction -> Identity, Frame -> {{True, True}, {True, True}}, FrameLabel -> {{None, None}, {None, None}}, FrameStyle -> Directive[ Thickness[Tiny], GrayLevel[0.7]], FrameTicks -> {{None, None}, {None, None}}, GridLines -> {None, None}, GridLinesStyle -> Directive[ GrayLevel[0.5, 0.4]], ImageSize -> {Automatic, Dynamic[ 3.5 (CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[ Magnification])]}, LabelStyle -> {FontFamily -> "Arial"}, Method -> {"ScalingFunctions" -> None}, PlotRange -> {All, All}, PlotRangeClipping -> True, PlotRangePadding -> {Automatic, Automatic}, Ticks -> {Automatic, Automatic}], GridBox[{{ RowBox[{ TagBox["\"Expression: \"", "SummaryItemAnnotation"], "\[InvisibleSpace]", TagBox[ TagBox["y", HoldForm], "SummaryItem"]}]}, { RowBox[{ TagBox["\"Parameters: \"", "SummaryItemAnnotation"], "\[InvisibleSpace]", TagBox[ RowBox[{"{", "ee", "}"}], "SummaryItem"]}]}, { RowBox[{ TagBox["\"Generator: \"", "SummaryItemAnnotation"], "\[InvisibleSpace]", TagBox["ParametricNDSolve", "SummaryItem"]}]}, { RowBox[{ TagBox[ "\"Dependent variables: \"", "SummaryItemAnnotation"], "\[InvisibleSpace]", TagBox[ RowBox[{"{", "y", "}"}], "SummaryItem"]}]}, { RowBox[{ TagBox[ "\"Independent variables: \"", "SummaryItemAnnotation"], "\[InvisibleSpace]", TagBox[ RowBox[{"{", "x", "}"}], "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->("ParametricFunction[<>]"& )], False, Editable->False, SelectWithContents->True, Selectable->False]}], "}"}]], "Output", CellChangeTimes->{ 3.7200527068456707`*^9, 3.720054382423677*^9, 3.7200545366356497`*^9, 3.7200546509519*^9, 3.720054682857225*^9, 3.720054718595395*^9, 3.720054922859538*^9, 3.720054977767272*^9, 3.7200550665191107`*^9, 3.7200937450842943`*^9, 3.7200937776728354`*^9, {3.720093820070238*^9, 3.72009384870443*^9}, {3.72009422128966*^9, 3.720094248986588*^9}, 3.7200946407894697`*^9, 3.720095274127344*^9, 3.72009550331115*^9, 3.720095609034442*^9, 3.720095720757526*^9, 3.720096406540399*^9, { 3.720096799390588*^9, 3.7200968140334272`*^9}, 3.720098375824091*^9, 3.720098607769125*^9, {3.720098875147534*^9, 3.720098894089326*^9}, 3.7200996037414017`*^9, 3.720106227248007*^9, 3.720106769909294*^9, { 3.72011701873006*^9, 3.7201170287110167`*^9}}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"y2", "[", RowBox[{ RowBox[{"e_", "?", "NumericQ"}], ",", " ", RowBox[{"x_", "/;", " ", RowBox[{ "xmid", "\[LessEqual]", " ", "x", " ", "\[LessEqual]", " ", "xmax"}]}]}], "]"}], " ", ":=", " ", RowBox[{"Evaluate", "[", " ", RowBox[{ RowBox[{ RowBox[{"y", "[", "e", "]"}], "[", "x", "]"}], " ", "/.", " ", "psol2"}], "]"}]}]], "Input", CellChangeTimes->{{3.720097058177556*^9, 3.720097118598814*^9}, { 3.720098898095475*^9, 3.7200989019125433`*^9}, {3.72009914139384*^9, 3.7200991458556623`*^9}, 3.720099592235539*^9, {3.720101067774225*^9, 3.720101074443303*^9}, {3.720101116292568*^9, 3.720101133252392*^9}, { 3.7201016737588663`*^9, 3.720101703428467*^9}, {3.7201059310010567`*^9, 3.720105935569765*^9}, {3.720106077573945*^9, 3.7201061157020617`*^9}, { 3.720106670824196*^9, 3.7201066712178297`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"dy2", "[", RowBox[{ RowBox[{"e_", "?", "NumericQ"}], ",", " ", RowBox[{"x_", "/;", " ", RowBox[{ "xmid", " ", "\[LessEqual]", " ", "x", " ", "\[LessEqual]", " ", "xmax"}]}]}], "]"}], " ", ":=", " ", RowBox[{"Evaluate", "[", " ", RowBox[{ RowBox[{ RowBox[{ RowBox[{"y", "[", "e", "]"}], "'"}], "[", "x", "]"}], " ", "/.", " ", "psol2"}], "]"}]}]], "Input", CellChangeTimes->{{3.720097381189123*^9, 3.720097396271584*^9}, { 3.720097722693634*^9, 3.720097723071767*^9}, {3.720098905651449*^9, 3.720098909801659*^9}, {3.7200992910765944`*^9, 3.7200992928929367`*^9}, { 3.720099598019903*^9, 3.720099598912471*^9}, {3.720099655327845*^9, 3.720099657799645*^9}, {3.7201010772534857`*^9, 3.720101080715218*^9}, { 3.7201011181867323`*^9, 3.720101119543161*^9}, {3.72010593863481*^9, 3.720105939311638*^9}, {3.7201061191448812`*^9, 3.7201061294384613`*^9}, { 3.7201062329623537`*^9, 3.7201062336194*^9}, {3.720106673104115*^9, 3.720106675521823*^9}}] }, Open ]] }, WindowSize->{808, 755}, WindowMargins->{{183, Automatic}, {39, Automatic}}, FrontEndVersion->"10.0 for Mac OS X x86 (32-bit, 64-bit Kernel) (June 27, \ 2014)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[580, 22, 138, 1, 35, "Subsubsection"], Cell[721, 25, 361, 7, 28, "Input"], Cell[CellGroupData[{ Cell[1107, 36, 477, 7, 63, "Input"], Cell[1587, 45, 100, 1, 28, "Output"], Cell[1690, 48, 95, 1, 28, "Output"], Cell[1788, 51, 94, 1, 28, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[1931, 58, 109, 1, 35, "Subsubsection"], Cell[CellGroupData[{ Cell[2065, 63, 362, 11, 46, "Input"], Cell[2430, 76, 214, 7, 57, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[2693, 89, 125, 1, 35, "Subsubsection"], Cell[2821, 92, 409, 12, 28, "Input"], Cell[3233, 106, 476, 13, 28, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[3746, 124, 130, 1, 35, "Subsubsection"], Cell[CellGroupData[{ Cell[3901, 129, 1469, 36, 46, "Input"], Cell[5373, 167, 18863, 344, 65, "Output"] }, Open ]], Cell[24251, 514, 886, 20, 28, "Input"], Cell[25140, 536, 998, 22, 28, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[26175, 563, 140, 1, 35, "Subsubsection"], Cell[CellGroupData[{ Cell[26340, 568, 1705, 40, 63, "Input"], Cell[28048, 610, 19394, 351, 65, "Output"] }, Open ]], Cell[47457, 964, 907, 20, 28, "Input"], Cell[48367, 986, 1055, 23, 28, "Input"] }, Open ]] } ] *) (* End of internal cache information *)