(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 11.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 158, 7] NotebookDataLength[ 45858, 1156] NotebookOptionsPosition[ 44172, 1097] NotebookOutlinePosition[ 44551, 1113] CellTagsIndexPosition[ 44508, 1110] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell["Gauss-Newton", "Section", CellChangeTimes->{{3.76118807524448*^9, 3.7611880882904496`*^9}, { 3.7624010936216574`*^9, 3.7624010939341407`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"GaussNewton", "[", RowBox[{"r_", ",", "p0_", ",", "tol_", ",", "stop_"}], "]"}], ":=", RowBox[{"(", "\[IndentingNewLine]", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", "s", "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"iter", "=", "0"}], ";", "\[IndentingNewLine]", RowBox[{"temp", "=", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Length", "[", "p0", "]"}], "\[Equal]", "1"}], ",", RowBox[{"{", "}"}], ",", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{ RowBox[{"Length", "[", "p0", "]"}], "-", "1"}]}], "]"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"\[Epsilon]", "=", SuperscriptBox["10", RowBox[{"-", "6"}]]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"Jacobi", "[", "p_", "]"}], ":=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{ FractionBox["1", "\[Epsilon]"], RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"r", "[", RowBox[{"p", "+", RowBox[{"Insert", "[", RowBox[{"temp", ",", "\[Epsilon]", ",", "k"}], "]"}]}], "]"}], "[", RowBox[{"[", "s", "]"}], "]"}], "-", RowBox[{ RowBox[{"r", "[", "p", "]"}], "[", RowBox[{"[", "s", "]"}], "]"}]}], ")"}]}], ",", RowBox[{"{", RowBox[{"s", ",", "1", ",", RowBox[{"Length", "[", RowBox[{"r", "[", "p", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"k", ",", "1", ",", RowBox[{"Length", "[", "p", "]"}]}], "}"}]}], "]"}], "//", "Simplify"}]}], ";", "\[IndentingNewLine]", RowBox[{"pOld", "=", "p0"}], ";", "\[IndentingNewLine]", RowBox[{"pNew", "=", RowBox[{"p0", "+", RowBox[{"Table", "[", RowBox[{"0.1", ",", RowBox[{"{", RowBox[{"k", ",", RowBox[{"Length", "[", "p0", "]"}]}], "}"}]}], "]"}]}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"While", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"Norm", "[", RowBox[{"pOld", "-", "pNew"}], "]"}], ">", "tol"}], " ", "&&", " ", RowBox[{"iter", "\[LessEqual]", "stop"}]}], " ", ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"pOld", "=", "pNew"}], ";", "\[IndentingNewLine]", RowBox[{"\[CapitalDelta]p", "=", RowBox[{"LeastSquares", "[", RowBox[{ RowBox[{"Jacobi", "[", "pOld", "]"}], ",", RowBox[{"-", RowBox[{"r", "[", "pOld", "]"}]}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"\[Alpha]", "=", "1"}], ";", "\[IndentingNewLine]", RowBox[{"pNew", "=", RowBox[{"pOld", "+", RowBox[{"\[Alpha]", " ", "\[CapitalDelta]p"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"While", "[", RowBox[{ RowBox[{ RowBox[{"Norm", "[", RowBox[{ RowBox[{"r", "[", "pNew", "]"}], ",", "2"}], "]"}], ">", RowBox[{"Norm", "[", RowBox[{ RowBox[{"r", "[", "pOld", "]"}], ",", "2"}], "]"}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Alpha]", "/=", "2"}], ";", "\[IndentingNewLine]", RowBox[{"pNew", "=", RowBox[{"pOld", "+", RowBox[{"\[Alpha]", " ", "\[CapitalDelta]p"}]}]}], ";"}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", RowBox[{"iter", "++"}], ";"}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", "pNew"}]}], "\[IndentingNewLine]", "]"}], "\[IndentingNewLine]", ")"}]}]], "Input", CellChangeTimes->{{3.7612461135414705`*^9, 3.761246124560611*^9}, 3.761247354462191*^9, {3.761250038662279*^9, 3.761250039772319*^9}, { 3.761830980013237*^9, 3.761830980169476*^9}, {3.761832130439988*^9, 3.7618321307299147`*^9}, {3.761833375799801*^9, 3.7618333892198925`*^9}, 3.761841169634926*^9, {3.76184125365444*^9, 3.7618412577335267`*^9}, 3.761842700401864*^9, {3.7618971859796963`*^9, 3.761897186111258*^9}, { 3.7618972307222824`*^9, 3.761897238991739*^9}, 3.7618976341651*^9, 3.761902167487347*^9, {3.7619022470659575`*^9, 3.761902308477112*^9}}] }, Closed]], Cell[CellGroupData[{ Cell["Diffusion parameter estimation ", "Section", CellChangeTimes->{{3.7618310283749423`*^9, 3.7618310371417418`*^9}, { 3.7623986889209166`*^9, 3.7623986890615306`*^9}, {3.7623989511670523`*^9, 3.762398957291571*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"HeatEquation", "[", RowBox[{"d_", ",", "T_", ",", "l_", ",", "h_"}], "]"}], ":=", RowBox[{"(", "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Tau]", "=", RowBox[{"Min", "[", RowBox[{ RowBox[{"Abs", "[", RowBox[{ RowBox[{"h", "^", "2"}], " ", "/", RowBox[{"(", RowBox[{"4", "d"}], ")"}]}], "]"}], ",", "0.01"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"n", "=", RowBox[{"Ceiling", "[", RowBox[{"l", "/", "h"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"m", "=", RowBox[{"Ceiling", "[", RowBox[{"T", "/", "\[Tau]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"t", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"j", "-", "1"}], ")"}], "\[Tau]"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", RowBox[{"m", "+", "1"}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"x", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"i", "-", "1"}], ")"}], "h"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", RowBox[{"n", "+", "1"}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"y", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", RowBox[{"n", "+", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", RowBox[{"m", "+", "1"}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"initial", " ", "condition"}], "*)"}], "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", " ", RowBox[{"n", "+", "1"}]}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"y", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "=", RowBox[{"Sin", "[", RowBox[{"2", " ", "Pi", " ", RowBox[{"(", RowBox[{"i", "-", "1"}], ")"}], "h"}], "]"}]}], ";"}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", "1"}], ",", RowBox[{"j", "\[LessEqual]", " ", "m"}], ",", RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "2"}], ",", RowBox[{"i", "\[LessEqual]", " ", "n"}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"y", "[", RowBox[{"[", RowBox[{"i", ",", RowBox[{"j", "+", "1"}]}], "]"}], "]"}], "=", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"1", "-", RowBox[{"2", "d", " ", RowBox[{"\[Tau]", "/", RowBox[{"h", "^", "2"}]}]}]}], ")"}], " ", RowBox[{"y", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}]}], "+", RowBox[{"d", " ", RowBox[{"\[Tau]", "/", RowBox[{"h", "^", "2"}]}], " ", RowBox[{"(", RowBox[{ RowBox[{"y", "[", RowBox[{"[", RowBox[{ RowBox[{"i", "-", "1"}], ",", "j"}], "]"}], "]"}], "+", RowBox[{"y", "[", RowBox[{"[", RowBox[{ RowBox[{"i", "+", "1"}], ",", "j"}], "]"}], "]"}]}], ")"}]}]}]}], ";"}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "[", RowBox[{"[", RowBox[{"1", ",", RowBox[{"j", "+", "1"}]}], "]"}], "]"}], "=", "0"}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "[", RowBox[{"[", RowBox[{ RowBox[{"n", "+", "1"}], ",", RowBox[{"j", "+", "1"}]}], "]"}], "]"}], "=", "0"}], ";"}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", "y"}], "\[IndentingNewLine]", ")"}]}]], "Input", CellChangeTimes->{{3.7623988078184967`*^9, 3.762398852049658*^9}, { 3.7623989444175434`*^9, 3.762398947182967*^9}, {3.7623990400199747`*^9, 3.7623990767673025`*^9}, {3.762401156273359*^9, 3.762401156617081*^9}, { 3.763096348792362*^9, 3.7630963512917747`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"res", "[", "D_", "]"}], ":=", RowBox[{"(", "\[IndentingNewLine]", RowBox[{ RowBox[{"d", "=", RowBox[{"D", "[", RowBox[{"[", "1", "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"y", "=", RowBox[{"HeatEquation", "[", RowBox[{"d", ",", "T", ",", "l", ",", "h"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"interpolation", "=", RowBox[{"Interpolation", "[", RowBox[{"Flatten", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"t", "[", RowBox[{"[", "j", "]"}], "]"}], ",", RowBox[{"y", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", RowBox[{"Length", "[", "x", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", RowBox[{"Length", "[", "t", "]"}]}], "}"}]}], "]"}], ",", "1"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"Join", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"interpolation", "[", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], ",", "1"}], "]"}], "-", RowBox[{"u1", "[", RowBox[{"[", "i", "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", RowBox[{"Length", "[", "x", "]"}]}], "}"}]}], "]"}], ",", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"interpolation", "[", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], ",", "2"}], "]"}], "-", RowBox[{"u2", "[", RowBox[{"[", "i", "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", RowBox[{"Length", "[", "x", "]"}]}], "}"}]}], "]"}]}], "]"}]}], "\[IndentingNewLine]", ")"}]}], ";"}]], "Input", CellChangeTimes->{{3.7623987102943296`*^9, 3.7623987111223927`*^9}, { 3.7624011724596834`*^9, 3.762401187661703*^9}, 3.7630963559312077`*^9, { 3.7630964007018657`*^9, 3.7630964012018695`*^9}, {3.7630964398644996`*^9, 3.763096498288415*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"T", "=", "2"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"l", "=", "1"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"h", "=", "0.2"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tol", "=", SuperscriptBox["10", RowBox[{"-", "5"}]]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"stop", "=", "100"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"u1", "=", RowBox[{"{", RowBox[{"0", ",", "0.6", ",", "0.45", ",", RowBox[{"-", "0.4"}], ",", RowBox[{"-", "0.6"}], ",", "0"}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"u2", "=", RowBox[{"{", RowBox[{"0", ",", "0.4", ",", "0.3", ",", RowBox[{"-", "0.27"}], ",", RowBox[{"-", "0.47"}], ",", "0"}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"diff", "=", RowBox[{"GaussNewton", "[", RowBox[{"res", ",", RowBox[{"{", "0.1", "}"}], ",", "tol", ",", "stop"}], "]"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"solution", "=", RowBox[{"HeatEquation", "[", RowBox[{"diff", ",", "T", ",", "l", ",", "h"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"x", "=", RowBox[{"Range", "[", RowBox[{"0", ",", "1", ",", "0.2"}], "]"}]}], ";"}]}], "Input", CellChangeTimes->{{3.7623988714076276`*^9, 3.762398878344618*^9}, { 3.762398965650378*^9, 3.7623990111314487`*^9}, {3.762399173446514*^9, 3.7623991822271233`*^9}, {3.7623998079523716`*^9, 3.7623998081398945`*^9}, 3.762399948098442*^9, {3.7624001068057003`*^9, 3.762400111617855*^9}, { 3.762401341621183*^9, 3.7624013698378286`*^9}, {3.763096495882801*^9, 3.7630965061303835`*^9}}], Cell[BoxData[ RowBox[{"{", "0.011445983169324409`", "}"}]], "Output", CellChangeTimes->{ 3.762400115305028*^9, {3.762401205334154*^9, 3.762401221004892*^9}, { 3.762401342433571*^9, 3.7624013477613244`*^9}, 3.7630963693342333`*^9, { 3.76309647351287*^9, 3.763096507879792*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Show", "[", RowBox[{ RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"Transpose", "[", RowBox[{"{", RowBox[{"x", ",", "u1"}], "}"}], "]"}], ",", RowBox[{"PlotStyle", "\[Rule]", " ", "Red"}], ",", RowBox[{"PlotRange", "\[Rule]", " ", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "1.1"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", "1"}], "}"}]}], "}"}]}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"Transpose", "[", RowBox[{"{", RowBox[{"x", ",", RowBox[{"Flatten", "[", RowBox[{"solution", "[", RowBox[{"[", RowBox[{"All", ",", "101"}], "]"}], "]"}], "]"}]}], "}"}], "]"}], ",", RowBox[{"PlotStyle", "\[Rule]", " ", "Black"}]}], "]"}]}], "\[IndentingNewLine]", "]"}]], "Input", CellChangeTimes->{{3.7618980666359816`*^9, 3.761898071325144*^9}, { 3.7618981979587746`*^9, 3.7618982149799385`*^9}, {3.762399155839685*^9, 3.762399165243984*^9}, {3.7623992174276934`*^9, 3.7623992214586515`*^9}, { 3.7623992514877167`*^9, 3.7623992536750574`*^9}, {3.7623996253875456`*^9, 3.7623996263562264`*^9}, 3.7624000144686613`*^9, {3.7624000654337068`*^9, 3.7624000694490395`*^9}, 3.762401318279126*^9}], Cell[BoxData[ GraphicsBox[{{{}, {{}, {RGBColor[1, 0, 0], PointSize[0.012833333333333334`], AbsoluteThickness[ 1.6], PointBox[{{0., 0.}, {0.2, 0.6}, {0.4, 0.45}, { 0.6000000000000001, -0.4}, {0.8, -0.6}, {1., 0.}}]}, {}}, {}, {}, {{}, {}}}, {{}, {{}, {GrayLevel[0], PointSize[0.012833333333333334`], AbsoluteThickness[1.6], PointBox[{{0., 0.}, {0.2, 0.6399183111777031}, {0.4, 0.3954912663312523}, {0.6000000000000001, -0.3954912663312523}, { 0.8, -0.6399183111777033}, {1., 0.}}]}, {}}, {}, {}, {{}, {}}}}, 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]], ImagePadding->All, Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({ (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[ Part[#, 1]], (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[ Part[#, 2]]}& ), "CopiedValueFunction" -> ({ (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[ Part[#, 1]], (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[ Part[#, 2]]}& )}}, PlotRange->{{0, 1.1}, {-1, 1}}, PlotRangeClipping->True, PlotRangePadding->{{0, 0}, {0, 0}}, Ticks->{Automatic, Automatic}]], "Output", CellChangeTimes->{{3.761897452103362*^9, 3.7618974805437927`*^9}, 3.7618976886141005`*^9, {3.7618979979381695`*^9, 3.761898080403367*^9}, { 3.761898215850362*^9, 3.7618982231425023`*^9}, {3.762399151402471*^9, 3.7623991843050594`*^9}, 3.7623992217555037`*^9, 3.7623992543780975`*^9, 3.762399627340476*^9, 3.7624000156403966`*^9, 3.7624001171799107`*^9, 3.762401222973516*^9, 3.7624012773757997`*^9, 3.762401318857173*^9, 3.763096372942726*^9, 3.7630965097855673`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Show", "[", RowBox[{ RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"Transpose", "[", RowBox[{"{", RowBox[{"x", ",", "u2"}], "}"}], "]"}], ",", RowBox[{"PlotStyle", "\[Rule]", " ", "Red"}], ",", RowBox[{"PlotRange", "\[Rule]", " ", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "1.1"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", "1"}], "}"}]}], "}"}]}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"Transpose", "[", RowBox[{"{", RowBox[{"x", ",", RowBox[{"Flatten", "[", RowBox[{"solution", "[", RowBox[{"[", RowBox[{"All", ",", "201"}], "]"}], "]"}], "]"}]}], "}"}], "]"}], ",", RowBox[{"PlotStyle", "\[Rule]", " ", "Black"}]}], "]"}]}], "\[IndentingNewLine]", "]"}]], "Input", CellChangeTimes->{{3.7618974595507283`*^9, 3.7618974635333962`*^9}, { 3.761897671487417*^9, 3.761897679117052*^9}, {3.7618982291953063`*^9, 3.7618982705674872`*^9}, {3.7624000288894835`*^9, 3.7624000598247366`*^9}, { 3.762401242597126*^9, 3.762401244643856*^9}, {3.7624012825785437`*^9, 3.7624013218257475`*^9}}], Cell[BoxData[ GraphicsBox[{{{}, {{}, {RGBColor[1, 0, 0], PointSize[0.012833333333333334`], AbsoluteThickness[ 1.6], PointBox[{{0., 0.}, {0.2, 0.4}, {0.4, 0.3}, { 0.6000000000000001, -0.27}, {0.8, -0.47}, {1., 0.}}]}, {}}, {}, {}, {{}, {}}}, {{}, {{}, {GrayLevel[0], PointSize[0.012833333333333334`], AbsoluteThickness[1.6], PointBox[{{0., 0.}, {0.2, 0.4305689913946605}, {0.4, 0.2661062711836612}, {0.6000000000000001, -0.2661062711836612}, { 0.8, -0.4305689913946606}, {1., 0.}}]}, {}}, {}, {}, {{}, {}}}}, 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]], ImagePadding->All, Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({ (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[ Part[#, 1]], (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[ Part[#, 2]]}& ), "CopiedValueFunction" -> ({ (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[ Part[#, 1]], (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[ Part[#, 2]]}& )}}, PlotRange->{{0, 1.1}, {-1, 1}}, PlotRangeClipping->True, PlotRangePadding->{{0, 0}, {0, 0}}, Ticks->{Automatic, Automatic}]], "Output", CellChangeTimes->{{3.7618982440563745`*^9, 3.7618982713701296`*^9}, { 3.762400032451686*^9, 3.762400060855863*^9}, 3.76240012056896*^9, { 3.762401224926481*^9, 3.76240124526877*^9}, {3.762401278828828*^9, 3.7624013221537943`*^9}, 3.763096374192433*^9, 3.7630965116288958`*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Monod model", "Section", CellChangeTimes->{{3.7618979671355076`*^9, 3.761897968437396*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"resid", "[", "p_", "]"}], ":=", RowBox[{"(", "\[IndentingNewLine]", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{"s", ",", "x"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"a", "=", RowBox[{"p", "[", RowBox[{"[", "1", "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"b", "=", RowBox[{"p", "[", RowBox[{"[", "2", "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"s0", "=", RowBox[{"p", "[", RowBox[{"[", "3", "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"x0", "=", RowBox[{"p", "[", RowBox[{"[", "4", "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"Res", "[", "t_", "]"}], "=", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"s", "[", "t", "]"}], ",", RowBox[{"x", "[", "t", "]"}]}], "}"}], "/.", RowBox[{ RowBox[{"NDSolve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{"s", "'"}], "[", "t", "]"}], "\[Equal]", RowBox[{"1", "-", RowBox[{"s", "[", "t", "]"}], "-", " ", RowBox[{"a", " ", RowBox[{ RowBox[{"s", "[", "t", "]"}], "/", RowBox[{"(", RowBox[{"b", "+", RowBox[{"s", "[", "t", "]"}]}], ")"}]}], " ", RowBox[{"x", "[", "t", "]"}]}]}]}], ",", RowBox[{ RowBox[{ RowBox[{"x", "'"}], "[", "t", "]"}], "\[Equal]", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"a", " ", RowBox[{ RowBox[{"s", "[", "t", "]"}], "/", RowBox[{"(", RowBox[{"b", "+", RowBox[{"s", "[", "t", "]"}]}], ")"}]}]}], "-", "1"}], ")"}], RowBox[{"x", "[", "t", "]"}]}]}], ",", RowBox[{ RowBox[{"s", "[", "0", "]"}], "\[Equal]", " ", "s0"}], ",", RowBox[{ RowBox[{"x", "[", "0", "]"}], "\[Equal]", " ", "x0"}]}], "}"}], ",", RowBox[{"{", RowBox[{"s", ",", "x"}], "}"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "2"}], "}"}]}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"Join", "[", RowBox[{ RowBox[{ RowBox[{"Res", "[", RowBox[{"Range", "[", RowBox[{"0", ",", "2", ",", "0.5"}], "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{ RowBox[{"Res", "[", RowBox[{"Range", "[", RowBox[{"0", ",", "2", ",", "0.5"}], "]"}], "]"}], "[", RowBox[{"[", "2", "]"}], "]"}]}], "]"}], "-", RowBox[{"{", RowBox[{ "1", ",", "1", ",", "0.87", ",", "0.52", ",", "0.04", ",", "0.01", ",", "0.04", ",", "0.14", ",", "0.48", ",", "0.965"}], "}"}]}]}]}], "\[IndentingNewLine]", "]"}], "\[IndentingNewLine]", ")"}]}]], "Input", CellChangeTimes->{{3.761901758066984*^9, 3.7619017763824043`*^9}, { 3.76190204301689*^9, 3.7619020459768076`*^9}, {3.762400146176508*^9, 3.7624001732370453`*^9}, {3.763096555540614*^9, 3.7630965774728193`*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"tol", "=", SuperscriptBox["10", RowBox[{"-", "6"}]]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"stop", "=", "50"}], ";"}], "\[IndentingNewLine]", RowBox[{"GaussNewton", "[", RowBox[{"resid", ",", RowBox[{"{", RowBox[{"1000", ",", "1000", ",", "1", ",", "0"}], "}"}], ",", "tol", ",", "stop"}], "]"}]}], "Input", CellChangeTimes->{{3.761901782980667*^9, 3.761901802429665*^9}, { 3.762400533039035*^9, 3.7624005431163783`*^9}, {3.76240065170337*^9, 3.7624006623276*^9}, {3.762400718026642*^9, 3.7624007181829147`*^9}, { 3.762400795130454*^9, 3.76240083583062*^9}, {3.762400966991584*^9, 3.762401050438691*^9}, {3.763096689805973*^9, 3.763096723751161*^9}}], Cell[BoxData[ TemplateBox[{ "NDSolve","ndsz", "\"At \\!\\(\\*RowBox[{\\\"t$\\\"}]\\) == \ \\!\\(\\*RowBox[{\\\"1.9023372950307762`\\\"}]\\), step size is effectively \ zero; singularity or stiff system suspected.\"",2,58,49,17248775219902964593, "Local"}, "MessageTemplate"]], "Message", "MSG", CellChangeTimes->{{3.7630966917740493`*^9, 3.763096724578864*^9}, { 3.7630967980928135`*^9, 3.7630968180724926`*^9}}], Cell[BoxData[ TemplateBox[{ "InterpolatingFunction","dmval", "\"Input value \\!\\(\\*RowBox[{\\\"{\\\", \\\"2.`\\\", \\\"}\\\"}]\\) \ lies outside the range of data in the interpolating function. Extrapolation \ will be used.\"",2,58,50,17248775219902964593,"Local"}, "MessageTemplate"]], "Message", "MSG", CellChangeTimes->{{3.7630966917740493`*^9, 3.763096724578864*^9}, { 3.7630967980928135`*^9, 3.7630968180881143`*^9}}], Cell[BoxData[ TemplateBox[{ "InterpolatingFunction","dmval", "\"Input value \\!\\(\\*RowBox[{\\\"{\\\", \\\"2.`\\\", \\\"}\\\"}]\\) \ lies outside the range of data in the interpolating function. Extrapolation \ will be used.\"",2,58,51,17248775219902964593,"Local"}, "MessageTemplate"]], "Message", "MSG", CellChangeTimes->{{3.7630966917740493`*^9, 3.763096724578864*^9}, { 3.7630967980928135`*^9, 3.763096818103737*^9}}], Cell[BoxData[ TemplateBox[{ "InterpolatingFunction","dmval", "\"Input value \\!\\(\\*RowBox[{\\\"{\\\", \\\"2.`\\\", \\\"}\\\"}]\\) \ lies outside the range of data in the interpolating function. Extrapolation \ will be used.\"",2,58,52,17248775219902964593,"Local"}, "MessageTemplate"]], "Message", "MSG", CellChangeTimes->{{3.7630966917740493`*^9, 3.763096724578864*^9}, { 3.7630967980928135`*^9, 3.7630968181193576`*^9}}], Cell[BoxData[ TemplateBox[{ "General","stop", "\"Further output of \\!\\(\\*StyleBox[RowBox[{\\\"InterpolatingFunction\\\ \", \\\"::\\\", \\\"dmval\\\"}], \\\"MessageName\\\"]\\) will be suppressed \ during this calculation.\"",2,58,53,17248775219902964593,"Local"}, "MessageTemplate"]], "Message", "MSG", CellChangeTimes->{{3.7630966917740493`*^9, 3.763096724578864*^9}, { 3.7630967980928135`*^9, 3.7630968181193576`*^9}}], Cell[BoxData[ TemplateBox[{ "NDSolve","ndsz", "\"At \\!\\(\\*RowBox[{\\\"t$\\\"}]\\) == \ \\!\\(\\*RowBox[{\\\"1.2046351143766452`\\\"}]\\), step size is effectively \ zero; singularity or stiff system suspected.\"",2,58,54,17248775219902964593, "Local"}, "MessageTemplate"]], "Message", "MSG", CellChangeTimes->{{3.7630966917740493`*^9, 3.763096724578864*^9}, { 3.7630967980928135`*^9, 3.7630968181662216`*^9}}], Cell[BoxData[ TemplateBox[{ "NDSolve","ndsz", "\"At \\!\\(\\*RowBox[{\\\"t$\\\"}]\\) == \ \\!\\(\\*RowBox[{\\\"1.2310392038236975`\\\"}]\\), step size is effectively \ zero; singularity or stiff system suspected.\"",2,58,55,17248775219902964593, "Local"}, "MessageTemplate"]], "Message", "MSG", CellChangeTimes->{{3.7630966917740493`*^9, 3.763096724578864*^9}, { 3.7630967980928135`*^9, 3.763096818244356*^9}}], Cell[BoxData[ TemplateBox[{ "General","stop", "\"Further output of \\!\\(\\*StyleBox[RowBox[{\\\"NDSolve\\\", \ \\\"::\\\", \\\"ndsz\\\"}], \\\"MessageName\\\"]\\) will be suppressed during \ this calculation.\"",2,58,56,17248775219902964593,"Local"}, "MessageTemplate"]], "Message", "MSG", CellChangeTimes->{{3.7630966917740493`*^9, 3.763096724578864*^9}, { 3.7630967980928135`*^9, 3.7630968182599497`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ "4.092624919861564`", ",", "0.11109174885391718`", ",", "1.012462819929313`", ",", "0.009373069256594067`"}], "}"}]], "Output", CellChangeTimes->{{3.762400547366068*^9, 3.7624005683187494`*^9}, 3.7624006575622683`*^9, 3.7624007206045637`*^9, 3.7624007985520396`*^9, { 3.762400828956066*^9, 3.7624008373773317`*^9}, {3.7624009624450293`*^9, 3.76240099923923*^9}, {3.7624010322212095`*^9, 3.7624010529697037`*^9}, 3.7630963842838116`*^9, {3.7630965715054407`*^9, 3.7630966067627773`*^9}, { 3.7630966953825755`*^9, 3.763096726203458*^9}, {3.7630967996705866`*^9, 3.763096819665868*^9}}] }, Open ]], Cell[BoxData[{ RowBox[{ RowBox[{"Clear", "[", RowBox[{"s", ",", "x"}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"a", "=", RowBox[{"pNew", "[", RowBox[{"[", "1", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"b", "=", RowBox[{"pNew", "[", RowBox[{"[", "2", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"s0", "=", RowBox[{"pNew", "[", RowBox[{"[", "3", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"x0", "=", RowBox[{"pNew", "[", RowBox[{"[", "4", "]"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"solution", "=", RowBox[{ RowBox[{"NDSolve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{"s", "'"}], "[", "t", "]"}], "\[Equal]", RowBox[{"1", "-", RowBox[{"s", "[", "t", "]"}], "-", " ", RowBox[{ FractionBox[ RowBox[{"a", " ", RowBox[{"s", "[", "t", "]"}]}], RowBox[{"b", "+", RowBox[{"s", "[", "t", "]"}]}]], " ", RowBox[{"x", "[", "t", "]"}]}]}]}], ",", RowBox[{ RowBox[{ RowBox[{"x", "'"}], "[", "t", "]"}], "\[Equal]", " ", RowBox[{ RowBox[{"(", RowBox[{ FractionBox[ RowBox[{"a", " ", RowBox[{"s", "[", "t", "]"}]}], RowBox[{"b", "+", RowBox[{"s", "[", "t", "]"}]}]], "-", "1"}], ")"}], RowBox[{"x", "[", "t", "]"}]}]}], ",", RowBox[{ RowBox[{"s", "[", "0", "]"}], "\[Equal]", " ", "s0"}], ",", RowBox[{ RowBox[{"x", "[", "0", "]"}], "\[Equal]", " ", "x0"}]}], "}"}], ",", RowBox[{"{", RowBox[{"x", ",", "s"}], "}"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "2"}], "}"}]}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"x", "[", "t_", "]"}], "=", RowBox[{ RowBox[{"x", "[", "t", "]"}], "/.", RowBox[{"solution", "[", RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"s", "[", "t_", "]"}], "=", RowBox[{ RowBox[{"s", "[", "t", "]"}], "/.", RowBox[{"solution", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}]}], "Input", CellChangeTimes->{{3.734356902660727*^9, 3.734356928536759*^9}, { 3.762400258652713*^9, 3.762400272417343*^9}, {3.7624005875205336`*^9, 3.7624005953480897`*^9}, {3.7630967793472195`*^9, 3.7630967840182076`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Show", "[", RowBox[{ RowBox[{"ListPlot", "[", RowBox[{"Transpose", "[", RowBox[{"{", RowBox[{ RowBox[{"Range", "[", RowBox[{"0", ",", "2", ",", "0.5"}], "]"}], ",", RowBox[{"{", RowBox[{"1", ",", "1", ",", "0.87", ",", "0.52", ",", "0.04"}], "}"}]}], "}"}], "]"}], "]"}], ",", RowBox[{"Plot", "[", RowBox[{ RowBox[{"s", "[", "i", "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "0", ",", "2"}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]}], "]"}]], "Input", CellChangeTimes->{{3.762400457310158*^9, 3.762400458028822*^9}, { 3.7624006257990036`*^9, 3.7624006399073534`*^9}}], Cell[BoxData[ GraphicsBox[{{{}, {{}, {RGBColor[0.368417, 0.506779, 0.709798], PointSize[ 0.012833333333333334`], AbsoluteThickness[1.6], PointBox[{{0., 1.}, {0.5, 1.}, {1., 0.87}, {1.5, 0.52}, {2., 0.04}}]}, {}}, {}, {}, {{}, {}}}, {{{}, {}, TagBox[ {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[1.], LineBox[CompressedData[" 1:eJwVkXk0lYkfxm/IkqJryVYjSySRVCakZ5Rkq8bSWEakrjJtopo2WRok1dCk Mqgs7732fcm4SpZhlCwJTbmUEC73XpO9e9/35/fHc57zOc853/Oc56t1JNDF X4xGo8Ut6v/u4D/S8WLUf2eQ3Q6z5VsEkNitfrZfzQYzb6K+TpgKsOVJYkqp mhtSHZh5bYvs+Lqn5JEaA6zOlpX3FtmdMfvsD7XzCHF6nqq2yJxNN3bfUYuE 4PQIZ91mAVKVD0THqiWgKhSlFpsEiJAu1b+pRoBU9tY4aCiAh2daoWCyDNE0 1dQEbQFkqQjvvd/VYOGknuU3NQHWOfNG6m814PNZ1W5lugAXW82RONWMC/VO YV2SAqiU0fuWLGnDrzs9vKzm+agY2N3StuMNKiQPuIV84qPuxCMNz/G3SB5k kk7/8DGtZyWff7oHxfeTlnIL+BCbz12pse495F9zdOXv8mFGHa5wLeiFvp9k VkUwHy2RW1sbN/Qjwyg81tuFD8147nw27yP2RCrqiQz54M7293WZDGDZhIxO jQQfxBj7pdzvnzHlsFt3STcPcmsy+YODg7i8SuKeRwYPt38pP5aqOwye/dps 90AeZn7IHQwI+wJG2NlCrc08nDdtKFtVPwJ5U4VztMkJ9Bc03OrQGsMrZnfX WmIC3068TTa4wsWTtRmNB3+cgKG1FHmMPY7KdE/25/lx7E3XOl6swsPVaScj 3T/GcXGNr/Ynbz7+2Gix/V+9cWwcPyjbs0+Az0+ib3wo4UJRXiDdsCBAK/+Y Sr4+FxI1yYU6KZMYqqr9hXw4hpNb2f/KOP6HIdd9/UXcUeTaH0l9PvMfquay T6qbjGKaq75y+MFXyL4x2vFd4AjKXfeKJ9hMIeznDQkBMV8QZhf1pGx4Cv4v t36yTh6GFsfhXUbUNN6ruVu0pgzB4vvIBDXTGYj3+LCCzgzC+7u2+SWdM1D2 87rhZvEZ0y/vDJyNmIVxbvDNttZPSJTwOe+iOwdPZvkqU9+PUOwrTxdrmMP9 l72T+X9xsPRahJt50DykfvdZnmb1HrRsn8IBhQXI78n3+vSmG8FJhXcU6hbQ MfqnqrJYJ54ZxkjTz3zDDGfKwTCqFXSH4I3X5YQI+ho0SHQ2IbVWXtaCLsTA 3IriV0+bYGKePzqpKMS5sL3GVHIT9q8fYfmpCRFeeK/+NqMJsZI+Wta6Qtgq aevYzjRCot5+Fc1CiHW3eOI/rm7EjKUWLdRfCI324VVl5xrQa9TeeblaCMaY dZq7ey2ecOnXxGuEMNi26aOjdS2OZLvq3akVQpUdnexiWItR3Z5LqY1ChL1j TzykXmBGo29NU7sQf21Oei2R/QL0ZePHlYaEaFNwmpsna2A7LCUsWCGC5pXr rdIXn0GGsGduXymCR9vyZu6BZ2jxu7W/TmExFyf2za9/BleOXFqXigiC664h CR+qcfitsq1QS4RTTiY6bbuqcbVOJ97eTIQqhzzTc2psFD2G7qCPCAq7bTfv GqnE1ytfpfz8RKizuD3n2VQJM/dMLueoCFG30ocfsSrBlpcv7QkQwf143c85 /pVoDO+zfhUsQu2ADat76Ck+HLnmWxItQlaNeGokrwKS+uyk8AIRjsbf7/XU Koe9eGCoqEgEPQXVTWKS5bjdr33kSqkIKSPGEj1jZaAnxhqcrxShQsxoMS7D ahmvyuN1IujvM0i/6FgGU+5c1/5uEc5sKVJJCS/FoUIz+hpSBFurjqsBsiWo 8JqTa6eRcCV6DhFfiyEnWbX8N3ESM1UaifTeYrzwtpIelSYh+2UqOCy/GNrL bKgyRRJ9lkqlGS7FGGY4TzgZkMjQvv5miFmEnSsVuZQhiUnp/mmtu0V4yH47 UmJMgjvnGn87pAh2Ch6DqltJXNNk24+7FiGnxufD4E4SkmH389okinBG/VRz iBuJ4/EM2zFOARr/Nmra5E7CLjSnS/hXATSD+A0DniRybuTu8H9QgPam4Bd2 vuRijy/J234sgOmFy0+VTpAguJuHG5vyMdsWzcwLIzG+7HqfVWMewqPSInpz SOxxUB+fEuZgxUUpB798Eumqz71iBnOQFHBaYbiQxDf3/F88WnJQ6miewS8j cZX2pSwvJQfD9I4GsRoSe98pqar/kAPHxzRpg04S9rHN9Vl3s/Eu7lh7fheJ 84aOOp0h2WBEtCSavlvc01Lj4JaAbIQyEg12cEgo/sO4eQ/ZKNpg4njgC4nH nu1ZTYIsrHrq+/uFbyQeNT5hpHtnISPr758WRCTymg+cnnPIgkmSoWYYjcJM kW3hNfMs2F+bLYxZSmF0NMBcSSULIbvjOpLlKZg2O4cee5uJgbbnSnU6FCz9 Do7Zu2cisFaXY6tHIbjPu55hlwlhSSzz1XoKOhtWsLLMM6H84CezLiMKQTEf qOzVmbDz5rmPfE8hbS7JLWaQhfyRNclyThQ21Bo3H73KQsyjqMN9+ylwG06G BZ1mgeEysa7AmcLjCPMtub4srK6uLtrnTqHDuSUn1YaF2LifG28fodAaXhoS vZKFALM/J5ddptAZHxjaXciEDZf29P1VCnP/eD04RTChmRoQkhNKIbLJMlb/ Tya6ZbZLOURSULk75L3sNyb2cLpX34yjcClzqfwlDya0I5XspFgU1tKDbBJl mCC3h6zoyaKwZ3WQjyONiX8nPr9h5VKgSW5UXjNLIN699JBt8eL9ICu++hAB aoPLuahqCjJx7OBf6wj0dsQ9Fu9c3Gey2dnnOoHK6NmjnV0UJn6qu9BxhUCC pa9BxjsK2mH2XYfOEXBkGpft6qNQcmDJiloGgapLr5sjRimw2WldWXYEHhht i3MeX+z3fMGDY00gaCDFTYtPgZei+HS9JYH1Tqf6X0xREMSb9IkZE5BY8paI n138V3qx6l19Ah/LLU8cXqDQ/TGEvk2LQPWJjE0mIgqBluqRPHUCDzVlpymK woM0+sZnSgT+B1Jzg2Y= "]]}, Annotation[#, "Charting`Private`Tag$416563#1"]& ]}, {}, {}}}, 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]], ImagePadding->All, Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({ (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[ Part[#, 1]], (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[ Part[#, 2]]}& ), "CopiedValueFunction" -> ({ (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[ Part[#, 1]], (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[ Part[#, 2]]}& )}}, PlotRange->{{0, 2.}, {0, 1.}}, PlotRangeClipping->True, PlotRangePadding->{{ Scaled[0.02], Scaled[0.02]}, { Scaled[0.02], Scaled[0.05]}}, Ticks->{Automatic, Automatic}]], "Output", CellChangeTimes->{{3.762400446795258*^9, 3.7624004583412848`*^9}, 3.762400600457054*^9, {3.7624006357513685`*^9, 3.7624006404854*^9}, { 3.7624007296976695`*^9, 3.7624007427748456`*^9}, 3.762400806254593*^9, 3.762400846875304*^9, {3.762400950242796*^9, 3.762400981256175*^9}, 3.762401011894562*^9, 3.762401060641034*^9, {3.7630967636790257`*^9, 3.7630968235712023`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Show", "[", RowBox[{ RowBox[{"ListPlot", "[", RowBox[{"Transpose", "[", RowBox[{"{", RowBox[{ RowBox[{"Range", "[", RowBox[{"0", ",", "2", ",", "0.5"}], "]"}], ",", RowBox[{"{", RowBox[{"0.01", ",", "0.04", ",", "0.14", ",", "0.48", ",", "0.965"}], "}"}]}], "}"}], "]"}], "]"}], ",", RowBox[{"Plot", "[", RowBox[{ RowBox[{"x", "[", "i", "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "0", ",", "2"}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}], ",", RowBox[{"PlotRange", "\[Rule]", " ", "All"}]}], "]"}]], "Input", CellChangeTimes->{ 3.7343601437902884`*^9, {3.734360293046259*^9, 3.7343603288445764`*^9}, { 3.762400479652272*^9, 3.762400487776699*^9}, {3.7624006174402337`*^9, 3.762400617596475*^9}}], Cell[BoxData[ GraphicsBox[{{{}, {{}, {RGBColor[0.368417, 0.506779, 0.709798], PointSize[ 0.012833333333333334`], AbsoluteThickness[1.6], PointBox[{{0., 0.01}, {0.5, 0.04}, {1., 0.14}, {1.5, 0.48}, {2., 0.965}}]}, {}}, {}, {}, {{}, {}}}, {{{}, {}, TagBox[ {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[1.], LineBox[CompressedData[" 1:eJwV0nk41PsXB/CxFN2yt6CEaKGUK0tU3pWSJVrGjUoNLi2UpdTcUqKhukrU lWTfKjeXSWgjsqTcq4tsLZhicjEzH6nsM9/vb35/nOc8r+c55zzPec7R9wna 5SfLYDDipPH/7OTX3/xiwM8Wig/NTM04kLfTDuZpbUbztB2cbTYcrM5ISi3W csOmk109AZs4cH7T8TBNyxdX+5kTJbs4cPcde35DKxTVelxO8AkOulZdsovV igK3xSM9tpSDzDnbL8ZoJeD+/o3JAeuiEKlYvPR3rVwsTa8wmvZLNDz2ZHG/ DpfgnuQUM23wEmbSkZ5bF1YinsVsZ6bGYPFO0l9zpRaUvSD07pVYsP+1RtKP eozJu3sYX43HvBK1bhmZRmg/O1mh53cDj3rsGhrXvYXTlGelXVYCqv3T5u8R toLt0cDbMvcWRpasVyk41gGDCb8IldzbkJ3IV52/+AOWjh/rvfksFZa01yNm YSdenrIy6RjOQEOU+b91xjzUhg0t5D/Jgm68YOJP8gm985x5m9tzIBjjdbeZ 9sA82W617oy7yB0s+1v5Wi+SDR91uBflQVnn3hCfz8eJi9fvW6Tm4+qR0oOZ hn14bD1Fr88pxOiGfP7h8//hMvHIS1tchFCz2pK5Nf1Qd+IxlWYVg1dYe6VZ fxDMnB+HvNxLMeXfmmJ0RoDBmsv1Dz4+xvKNCtTBMiEKstwqh9jPsDVb/1DR PAIeM8HsXtpzsHVYiz57DqHcNLhDN+QFVgh/mdnh8hU3bjh1t1tWQ0Plq2Lt 5Ff4szq9ez1qIV+ZwjVIHYbdrMRTStfqEGBe9n6G8zeIPBvYOyzqke/ok1kx +g2B+43WuGg3YESgrdqX+B019ftGPsxpRClzq1zC5h9QbvyiPcVsxnmH6IyS vh9QGvJXZxm3QL/L6V1O9AgmzQ8dTDRog41VVIKW2SiUWRnv2Gc74LmwcUKm ZRTjpl6fN3x+j5G/Y3uCI8cwyzuqt7SkE0nyB0J3GY6jvNdqblk8Dxrdpdmy teNQ0hT4Wzh+xrRzkW7WIRNYWrDkat60XjD+PMDtUZ+E5Rw95YynfBxP5saq V09i7XUrn9jMPjxffllRLXAK1v/VsAXh/VBzOr7igrIYbnn2LqrGg8isUplp oyaG976j7jI/D8LUumBgWEOMId/tlkNWg3Bd1n/XW0uMjU9m7ajaMoiY6Qf0 NxqK0VQ3xbD1HoR8jeNcho0YC+s1Xl6/NYjRtfqMcD8xeNkne/xkBOg0aWo5 XS5GBCGf378SIEOgdk6uUoy1eSt1PN4I4PMnc0lslRhxfPVnzW8FGDDs+C2z TowHf1jSj7oEGJ3frfOqSYygqHg5x+8CqP0kPDT7ixisRPOEaB0h7PsUxIVK ErwreHreOlCIGbmOd9aoSsDDTI7mCSEavK+4VqtLYHGp3eIbWwhml3JW2zwJ qoJI0s1IIbxa59iL9SXoX1/UnZIgRFi1QbyjpQQ2YUMl2U+FeJAOQ/4BCVJm H1ONoIX4fua7gre3BMIF8/OS5ESwdL8n6PpVgoZbGWX5CiKUqagUdxyWYJ5Z ys4qFRHqIro3/nNcgqQN/2jd1xXho8851sOLEmyXi1wYYivC9KVlyRGFEriM VUQ2sUVwlAsKlzyQ4KPHJx2rMBGu8hb5nCmWzt/qaZwULoJaUoxR6BMJMn98 SnaMFmHBjL1PDlVLoNl0QcX/hghmgvE213YJBF6sb9n5IuznWqrpUBKUP/er l/8gwqO948pNDAr241v/Gu8UQXn6s1kcOQoeroO+fTwRXniuVxxQpGD+uiyA +0WERT9tpks0KDxfeiRk/KsIfb47RduMKER06AUGTCewVdUQ0MspdKVv+YVW JLhV1tr/cCWFAR2XwNiZBA7qHnxNcwrDmjq7klQJ7lce+Mi3pRB0YTL3iDZB oPbR+rNuFFz12LcZJgR1L01erXKnsNWQ1W+zikA3ZKi2Zw+FjAVxH4J+Jmh6 dfyFA4vCOd7i1/UWBGYnTz+e7U8hb195jq0twVjjxTt/nadQFmhkvdqVICI6 K7LzPoUz621/jfAnUGIrOHkXULBetltR7ihB8uFj6n1cCrbLVfmRxwiKna1z hkooiDnZe0ODCfrUmmtlK6X9L1/LmpwicE5nKBq1UHAo7lr9PYLgXdzBpoI2 CtUhDqHrLhD4RjYkmb2jcOTUao8LHIJw3ySjdV0U+hSJofxFggfGps7b/6Mg v5Mb8ymGYO5j1rWTUxROe5u5yCYQ5OS93D0pkdYP7OjSvElgmrxc9zyDhtIf 4WdWJBI4nhvjXp5G4yDfa9QpieCsXVxzigqN1TFV2X6pBD2NFbOrDWisiehZ Pi+XIKjKsMt+CQ12+LLTMncIxA9j7vyzjMaWLyctB6Sek7jbss2ERuwOrlbR Pem9PIl7vxUNrmaOiUE+QUG/ToryNhqr0lQu3SwiuJwW7dXtSuOb7Hyb4IfS fXeJFhfupKFhNmelQzHBgvLyBy7uNPSVvMeGSghi4vbVXfWhsTvF8IvuE4LD lreHfzpNo2JSKVOhgmCzgPH4QxiNO6yA90+l1s08fPZ+OA0X5caXRyoJ2mes UXCKohG0e83b6hcEW7raF/weR0N2OOzE3hqCRVGzHRTu0thwqD528WsCas1Z pY48GtNztss9kvq9qPft3Xwa2wrfxG2uJ4h3L95vX0TDb9JR3vNvAtp414no chpOH4InQhsIOpvj0uVaaPADuJPMJoInF8d+bWmj8X3gedhrqRPWsoxy3tE4 3956ZG2z9B/urCzZ1E0jM6Msc8Fbgme/vamPHKCRfo0h+6aFINHEIm6nkEaW vt4Ji1aCkJ5UN/0hGvtl3dipUi/bdpT34geNygAX54NtBPIyrbnxYzT0Pk8a 1Ev9qXStv9ek9J6/cbKM2wnK/XNWmUpoeMWvfB0j9S3dmSM0TWPksl3qgNT/ A7/IkyU= "]]}, Annotation[#, "Charting`Private`Tag$416647#1"]& ]}, {}, {}}}, 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]], ImagePadding->All, Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({ (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[ Part[#, 1]], (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[ Part[#, 2]]}& ), "CopiedValueFunction" -> ({ (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[ Part[#, 1]], (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[ Part[#, 2]]}& )}}, PlotRange->All, PlotRangeClipping->True, PlotRangePadding->{{ Scaled[0.02], Scaled[0.02]}, { Scaled[0.02], Scaled[0.05]}}, Ticks->{Automatic, Automatic}]], "Output", CellChangeTimes->{ 3.7343601446809483`*^9, {3.7343603190941715`*^9, 3.734360329297714*^9}, 3.7619024427376475`*^9, 3.762400492245081*^9, {3.7624006082689867`*^9, 3.7624006469536753`*^9}, 3.7624007464776955`*^9, 3.762400851218746*^9, 3.762400983802864*^9, 3.762401014706875*^9, 3.7624010631252327`*^9, { 3.763096765553586*^9, 3.763096824602209*^9}}] }, Open ]] }, Open ]] }, WindowSize->{1280, 941}, WindowMargins->{{-8, Automatic}, {Automatic, -8}}, Magnification:>1.25 Inherited, FrontEndVersion->"11.0 for Microsoft Windows (64-bit) (September 21, 2016)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[580, 22, 151, 2, 86, "Section"], Cell[734, 26, 4573, 108, 724, "Input"] }, Closed]], Cell[CellGroupData[{ Cell[5344, 139, 224, 3, 68, "Section"], Cell[5571, 144, 4644, 128, 494, "Input"], Cell[10218, 274, 2441, 65, 158, "Input"], Cell[CellGroupData[{ Cell[12684, 343, 1692, 44, 254, "Input"], Cell[14379, 389, 286, 5, 37, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[14702, 399, 1381, 35, 86, "Input"], Cell[16086, 436, 2097, 41, 293, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[18220, 482, 1275, 34, 86, "Input"], Cell[19498, 518, 1894, 38, 293, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[21441, 562, 98, 1, 86, "Section"], Cell[21542, 565, 3524, 89, 326, "Input"], Cell[CellGroupData[{ Cell[25091, 658, 739, 16, 86, "Input"], Cell[25833, 676, 426, 9, 29, "Message"], Cell[26262, 687, 435, 8, 29, "Message"], Cell[26700, 697, 433, 8, 29, "Message"], Cell[27136, 707, 435, 8, 29, "Message"], Cell[27574, 717, 434, 8, 29, "Message"], Cell[28011, 727, 426, 9, 29, "Message"], Cell[28440, 738, 424, 9, 29, "Message"], Cell[28867, 749, 419, 8, 29, "Message"], Cell[29289, 759, 652, 11, 37, "Output"] }, Open ]], Cell[29956, 773, 2613, 76, 235, "Input"], Cell[CellGroupData[{ Cell[32594, 853, 712, 19, 37, "Input"], Cell[33309, 874, 4994, 95, 303, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[38340, 974, 852, 22, 62, "Input"], Cell[39195, 998, 4949, 95, 306, "Output"] }, Open ]] }, Open ]] } ] *)