(* 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[     22755,        607]
NotebookOptionsPosition[     21444,        560]
NotebookOutlinePosition[     21791,        575]
CellTagsIndexPosition[     21748,        572]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["\<\
Objective function, Jacobian matrix (definitions for logistic curve fitting \
problem)\
\>", "Section",
 CellChangeTimes->{{3.7594730634801903`*^9, 3.759473065454303*^9}, {
  3.7594731330520473`*^9, 3.759473151192085*^9}, {3.7594731831699133`*^9, 
  3.7594732038230953`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"f", "[", "t_", "]"}], ":=", 
   FractionBox["x0", 
    RowBox[{
     RowBox[{
      RowBox[{"(", 
       RowBox[{"1", "-", 
        RowBox[{"x0", "/", "k"}]}], ")"}], 
      SuperscriptBox["E", 
       RowBox[{
        RowBox[{"-", "r"}], " ", "t"}]]}], "+", 
     RowBox[{"x0", "/", "k"}]}]]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"values", "=", 
   RowBox[{"{", 
    RowBox[{
    "0.04", ",", "0.35", ",", "1.05", ",", "1.9", ",", "1.98", ",", "1.97", 
     ",", "1.99", ",", "1.999"}], "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"resid", "[", "p_", "]"}], ":=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"f", "[", 
       RowBox[{"Range", "[", 
        RowBox[{"0", ",", "7", ",", "1"}], "]"}], "]"}], "-", "values"}], 
     ")"}], "/.", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"r", "\[Rule]", 
       RowBox[{"p", "[", 
        RowBox[{"[", "1", "]"}], "]"}]}], ",", 
      RowBox[{"x0", "\[Rule]", 
       RowBox[{"p", "[", 
        RowBox[{"[", "2", "]"}], "]"}]}], ",", 
      RowBox[{"k", "\[Rule]", 
       RowBox[{"p", "[", 
        RowBox[{"[", "3", "]"}], "]"}]}]}], "}"}]}]}], ";"}], "\n", 
 RowBox[{"plot1", "=", 
  RowBox[{"ListPlot", "[", 
   RowBox[{
    RowBox[{"Transpose", "[", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Range", "[", 
        RowBox[{"0", ",", "7", ",", "1"}], "]"}], ",", "values"}], "}"}], 
     "]"}], ",", 
    RowBox[{"PlotMarkers", "\[Rule]", "\[FilledCircle]"}]}], 
   "]"}]}]}], "Input",
 CellChangeTimes->{{3.728066367113092*^9, 3.728066407334959*^9}, {
  3.7280664448413467`*^9, 3.728066446226323*^9}, {3.728066649610647*^9, 
  3.728066649874896*^9}, {3.728066900033571*^9, 3.7280669001508837`*^9}, {
  3.7594732383160677`*^9, 3.7594732752461805`*^9}, {3.7594733153234725`*^9, 
  3.7594733225608864`*^9}, {3.7594736309259176`*^9, 3.7594736680960436`*^9}, {
  3.759473698218766*^9, 3.759473702176993*^9}, {3.760701698579077*^9, 
  3.76070173890305*^9}}],

Cell[BoxData[
 GraphicsBox[{{}, {
    {RGBColor[0.368417, 0.506779, 0.709798], AbsolutePointSize[6], 
     AbsoluteThickness[1.6], 
     GeometricTransformationBox[InsetBox["\[FilledCircle]", {0., 0.}], {{{0., 
      0.04}}, {{1., 0.35}}, {{2., 1.05}}, {{3., 1.9}}, {{4., 1.98}}, {{5., 
      1.97}}, {{6., 1.99}}, {{7., 1.999}}}]}, {}}, {}, {}, {{}, {}}},
  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, 7.}, {0, 1.999}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.02], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{{3.759473658228479*^9, 3.759473704297114*^9}, 
   3.759475052852432*^9, 3.760701598175826*^9, {3.7607017357119455`*^9, 
   3.76070173955529*^9}}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"J", "[", "p_", "]"}], ":=", 
    RowBox[{
     RowBox[{
      RowBox[{"Transpose", "[", 
       RowBox[{"ArrayFlatten", "[", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"D", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"f", "[", 
              RowBox[{"Range", "[", 
               RowBox[{"0", ",", "7", ",", "1"}], "]"}], "]"}], "-", 
             "values"}], ",", "r"}], "]"}], ",", 
          RowBox[{"D", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"f", "[", 
              RowBox[{"Range", "[", 
               RowBox[{"0", ",", "7", ",", "1"}], "]"}], "]"}], "-", 
             "values"}], ",", "x0"}], "]"}], ",", "\[IndentingNewLine]", 
          RowBox[{"D", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"f", "[", 
              RowBox[{"Range", "[", 
               RowBox[{"0", ",", "7", ",", "1"}], "]"}], "]"}], "-", 
             "values"}], ",", "k"}], "]"}]}], "}"}], "]"}], "]"}], "/.", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"r", "\[Rule]", 
         RowBox[{"p", "[", 
          RowBox[{"[", "1", "]"}], "]"}]}], ",", 
        RowBox[{"x0", "\[Rule]", 
         RowBox[{"p", "[", 
          RowBox[{"[", "2", "]"}], "]"}]}], ",", 
        RowBox[{"k", "\[Rule]", 
         RowBox[{"p", "[", 
          RowBox[{"[", "3", "]"}], "]"}]}]}], "}"}]}], "//", "N"}]}], ";"}], 
  "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{{3.7607019703698344`*^9, 3.7607021266620474`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Levenberg--Marquardt", "Section",
 CellChangeTimes->{{3.7594765365202017`*^9, 3.7594765418795085`*^9}, {
  3.759478853243374*^9, 3.7594788908933663`*^9}, {3.760701358540709*^9, 
  3.760701366438812*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"pOld", "=", 
   RowBox[{"{", 
    RowBox[{"10", ",", "10", ",", "10"}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"pNew", "=", 
   RowBox[{"{", 
    RowBox[{"1", ",", "1", ",", "1"}], "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"tol", "=", 
   SuperscriptBox["10", 
    RowBox[{"-", "6"}]]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"iter", "=", "0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"maxIter", "=", "100"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Mu]", "=", "1"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"lastStepAccepted", "=", "True"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"residuals", "=", 
   RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"While", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{
       RowBox[{"Norm", "[", 
        RowBox[{"pOld", "-", "pNew"}], "]"}], ">", "tol"}], "||", 
      RowBox[{"lastStepAccepted", "\[Equal]", "False"}]}], ")"}], " ", "&&", 
    " ", 
    RowBox[{"iter", "<", "maxIter"}]}], ",", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"pOld", "=", "pNew"}], ";", "\[IndentingNewLine]", 
    StyleBox[
     RowBox[{"\[CapitalDelta]p", "=", 
      RowBox[{"LeastSquares", "[", 
       RowBox[{
        RowBox[{"Flatten", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"J", "[", "pOld", "]"}], ",", 
            RowBox[{
             SqrtBox["\[Mu]"], 
             RowBox[{"IdentityMatrix", "[", "3", "]"}]}]}], "}"}], ",", "1"}],
          "]"}], ",", 
        RowBox[{"Flatten", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"-", 
             RowBox[{"resid", "[", "pOld", "]"}]}], ",", " ", 
            RowBox[{"Table", "[", 
             RowBox[{"0", ",", 
              RowBox[{"{", "3", "}"}]}], "]"}]}], "}"}], ",", "1"}], "]"}]}], 
       "]"}]}],
     Background->RGBColor[1, 1, 0]], 
    StyleBox[";",
     Background->RGBColor[1, 1, 0]], "\[IndentingNewLine]", 
    RowBox[{"If", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"Norm", "[", 
        RowBox[{
         RowBox[{"resid", "[", 
          RowBox[{"pOld", "+", "\[CapitalDelta]p"}], "]"}], ",", "2"}], "]"}],
        "<", 
       RowBox[{"Norm", "[", 
        RowBox[{
         RowBox[{"resid", "[", "pOld", "]"}], ",", "2"}], "]"}]}], ",", 
      "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"pNew", "=", 
        RowBox[{"pOld", "+", " ", "\[CapitalDelta]p"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"\[Mu]", "/=", "3"}], ";", "\[IndentingNewLine]", 
       RowBox[{"lastStepAccepted", "=", "True"}]}], ",", 
      "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"\[Mu]", "*=", "2"}], ";", "\[IndentingNewLine]", 
       RowBox[{"lastStepAccepted", "=", "False"}]}]}], "\[IndentingNewLine]", 
     "]"}], ";", "\[IndentingNewLine]", 
    RowBox[{"Print", "[", "lastStepAccepted", "]"}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"AppendTo", "[", 
     RowBox[{"residuals", ",", " ", 
      RowBox[{"Norm", "[", 
       RowBox[{
        RowBox[{"resid", "[", 
         RowBox[{
          RowBox[{"pNew", "[", 
           RowBox[{"[", "1", "]"}], "]"}], ",", 
          RowBox[{"pNew", "[", 
           RowBox[{"[", "2", "]"}], "]"}], ",", 
          RowBox[{"pNew", "[", 
           RowBox[{"[", "3", "]"}], "]"}]}], "]"}], ",", "2"}], "]"}]}], 
     "]"}], ";", "\[IndentingNewLine]", 
    RowBox[{"iter", "++"}]}]}], "\[IndentingNewLine]", 
  "]"}], "\[IndentingNewLine]", "pNew", "\[IndentingNewLine]", "iter"}], \
"Input",
 CellChangeTimes->{{3.759473300406619*^9, 3.759473306635976*^9}, {
   3.759473800247602*^9, 3.759473943452793*^9}, {3.7594740291855626`*^9, 
   3.759474072685051*^9}, {3.7594743015848684`*^9, 3.75947432678431*^9}, 
   3.7594743592391663`*^9, {3.7594755821109686`*^9, 3.7594755821489706`*^9}, 
   3.7594758255377045`*^9, 3.759475899847536*^9, {3.7594760551194167`*^9, 
   3.759476123442028*^9}, {3.759477702632553*^9, 3.759477702832565*^9}, {
   3.759477966100357*^9, 3.7594780342692566`*^9}, {3.759478213494303*^9, 
   3.759478225450987*^9}, {3.759478915173755*^9, 3.7594789185849504`*^9}, {
   3.759478956409114*^9, 3.759478958983261*^9}, {3.760701402995467*^9, 
   3.760701486860176*^9}, {3.760701520522996*^9, 3.7607016851414194`*^9}, {
   3.7607018006164417`*^9, 3.7607019077331324`*^9}, {3.760702146816066*^9, 
   3.760702147313735*^9}, {3.7607021819824743`*^9, 3.7607022563803425`*^9}, {
   3.760702309037472*^9, 3.7607025565082335`*^9}, {3.7607025944738903`*^9, 
   3.760702609245425*^9}, {3.76070285457364*^9, 3.7607029601134963`*^9}}],

Cell[CellGroupData[{

Cell[BoxData["True"], "Print",
 CellChangeTimes->{{3.7607029093666396`*^9, 3.760702960913358*^9}}],

Cell[BoxData["True"], "Print",
 CellChangeTimes->{{3.7607029093666396`*^9, 3.7607029609163485`*^9}}],

Cell[BoxData["True"], "Print",
 CellChangeTimes->{{3.7607029093666396`*^9, 3.760702960924327*^9}}],

Cell[BoxData["True"], "Print",
 CellChangeTimes->{{3.7607029093666396`*^9, 3.76070296092832*^9}}],

Cell[BoxData["True"], "Print",
 CellChangeTimes->{{3.7607029093666396`*^9, 3.7607029609343023`*^9}}],

Cell[BoxData["True"], "Print",
 CellChangeTimes->{{3.7607029093666396`*^9, 3.760702960941282*^9}}],

Cell[BoxData["True"], "Print",
 CellChangeTimes->{{3.7607029093666396`*^9, 3.7607029609462686`*^9}}],

Cell[BoxData["True"], "Print",
 CellChangeTimes->{{3.7607029093666396`*^9, 3.7607029609512596`*^9}}],

Cell[BoxData["True"], "Print",
 CellChangeTimes->{{3.7607029093666396`*^9, 3.760702960957243*^9}}],

Cell[BoxData["True"], "Print",
 CellChangeTimes->{{3.7607029093666396`*^9, 3.760702960964222*^9}}],

Cell[BoxData["True"], "Print",
 CellChangeTimes->{{3.7607029093666396`*^9, 3.7607029609731956`*^9}}],

Cell[BoxData["True"], "Print",
 CellChangeTimes->{{3.7607029093666396`*^9, 3.7607029609781833`*^9}}],

Cell[BoxData["True"], "Print",
 CellChangeTimes->{{3.7607029093666396`*^9, 3.7607029609831705`*^9}}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "2.010591707366255`", ",", "0.043783632288488035`", ",", 
   "2.0065642690193806`"}], "}"}]], "Output",
 CellChangeTimes->{
  3.759473407737502*^9, {3.759473801072649*^9, 3.759473944574857*^9}, 
   3.759474073650106*^9, {3.759474302499921*^9, 3.759474327189333*^9}, 
   3.7594743597741966`*^9, 3.7594745315609674`*^9, 3.7594755829460163`*^9, 
   3.7594758259057255`*^9, 3.7594759003665657`*^9, {3.759476047110959*^9, 
   3.759476123933056*^9}, 3.759477703743617*^9, {3.7594779700235815`*^9, 
   3.75947800390752*^9}, 3.7594782166814857`*^9, 3.759478919052977*^9, {
   3.7607021386156564`*^9, 3.7607021490001383`*^9}, {3.760702185589113*^9, 
   3.760702212078062*^9}, {3.760702247848171*^9, 3.760702257196788*^9}, {
   3.7607023114984145`*^9, 3.760702332559862*^9}, {3.7607024468762627`*^9, 
   3.760702481311214*^9}, {3.760702518226778*^9, 3.7607025571520634`*^9}, {
   3.7607025958714952`*^9, 3.7607026215981264`*^9}, {3.760702855601885*^9, 
   3.7607028767224555`*^9}, {3.7607029094583406`*^9, 3.760702960989154*^9}}],

Cell[BoxData["13"], "Output",
 CellChangeTimes->{
  3.759473407737502*^9, {3.759473801072649*^9, 3.759473944574857*^9}, 
   3.759474073650106*^9, {3.759474302499921*^9, 3.759474327189333*^9}, 
   3.7594743597741966`*^9, 3.7594745315609674`*^9, 3.7594755829460163`*^9, 
   3.7594758259057255`*^9, 3.7594759003665657`*^9, {3.759476047110959*^9, 
   3.759476123933056*^9}, 3.759477703743617*^9, {3.7594779700235815`*^9, 
   3.75947800390752*^9}, 3.7594782166814857`*^9, 3.759478919052977*^9, {
   3.7607021386156564`*^9, 3.7607021490001383`*^9}, {3.760702185589113*^9, 
   3.760702212078062*^9}, {3.760702247848171*^9, 3.760702257196788*^9}, {
   3.7607023114984145`*^9, 3.760702332559862*^9}, {3.7607024468762627`*^9, 
   3.760702481311214*^9}, {3.760702518226778*^9, 3.7607025571520634`*^9}, {
   3.7607025958714952`*^9, 3.7607026215981264`*^9}, {3.760702855601885*^9, 
   3.7607028767224555`*^9}, {3.7607029094583406`*^9, 3.7607029609941387`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ListLogPlot", "[", 
  RowBox[{"residuals", ",", 
   RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.7594740757102237`*^9, 3.7594741577289147`*^9}}],

Cell[BoxData[
 GraphicsBox[{{}, {{}, 
    {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.012833333333333334`],
      AbsoluteThickness[1.6], 
     PointBox[{{1., -0.007317785203721156}, {2., -0.5198691106979378}, {
      3., -1.4646252850260084`}, {4., -1.717139690637691}, {
      5., -2.0307126779475557`}, {6., -2.033057736288337}, {
      7., -2.0331235070406595`}, {8., -2.033125545616949}, {
      9., -2.033125624172501}, {10., -2.0331256272897194`}, {
      11., -2.0331256274161262`}, {12., -2.033125627421278}, {
      13., -2.0331256274214877`}}]}, {}}, {}, {}, {{}, {}}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0.8125, -2.180160067582455},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{
     Charting`ScaledTicks[{Log, Exp}], 
     Charting`ScaledFrameTicks[{Log, Exp}]}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Part[{{Identity, Identity}, {Log, Exp}}, 1, 2][#]& )[
         Part[#, 1]], 
        (Part[{{Identity, Identity}, {Log, Exp}}, 2, 2][#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Part[{{Identity, Identity}, {Log, Exp}}, 1, 2][#]& )[
         Part[#, 1]], 
        (Part[{{Identity, Identity}, {Log, Exp}}, 2, 2][#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{{1., 13.}, {-2.1456705075446973`, -0.007317785203721156}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.02], 
     Scaled[0.05]}},
  Ticks->FrontEndValueCache[{Automatic, 
     Charting`ScaledTicks[{Log, Exp}]}, {Automatic, {{-1.6094379124341003`, 
       FormBox[
        TagBox[
         InterpretationBox["\"0.2\"", 0.2, AutoDelete -> True], NumberForm[#, {
           DirectedInfinity[1], 1.}]& ], TraditionalForm], {0.01, 0.}, {
        AbsoluteThickness[0.1]}}, {-0.6931471805599453, 
       FormBox[
        TagBox[
         InterpretationBox["\"0.5\"", 0.5, AutoDelete -> True], NumberForm[#, {
           DirectedInfinity[1], 1.}]& ], TraditionalForm], {0.01, 0.}, {
        AbsoluteThickness[0.1]}}, {0., 
       FormBox["1", TraditionalForm], {0.01, 0.}, {
        AbsoluteThickness[0.1]}}, {-2.3025850929940455`, 
       FormBox[
        InterpretationBox[
         StyleBox[
          
          GraphicsBox[{}, ImageSize -> {0., 0.}, BaselinePosition -> 
           Baseline], "CacheGraphics" -> False], 
         Spacer[{0., 0.}]], TraditionalForm], {0.005, 0.}, {
        AbsoluteThickness[0.1]}}, {-1.8971199848858813`, 
       FormBox[
        InterpretationBox[
         StyleBox[
          
          GraphicsBox[{}, ImageSize -> {0., 0.}, BaselinePosition -> 
           Baseline], "CacheGraphics" -> False], 
         Spacer[{0., 0.}]], TraditionalForm], {0.005, 0.}, {
        AbsoluteThickness[0.1]}}, {-1.2039728043259361`, 
       FormBox[
        InterpretationBox[
         StyleBox[
          
          GraphicsBox[{}, ImageSize -> {0., 0.}, BaselinePosition -> 
           Baseline], "CacheGraphics" -> False], 
         Spacer[{0., 0.}]], TraditionalForm], {0.005, 0.}, {
        AbsoluteThickness[0.1]}}, {-0.916290731874155, 
       FormBox[
        InterpretationBox[
         StyleBox[
          
          GraphicsBox[{}, ImageSize -> {0., 0.}, BaselinePosition -> 
           Baseline], "CacheGraphics" -> False], 
         Spacer[{0., 0.}]], TraditionalForm], {0.005, 0.}, {
        AbsoluteThickness[0.1]}}, {-0.5108256237659907, 
       FormBox[
        InterpretationBox[
         StyleBox[
          
          GraphicsBox[{}, ImageSize -> {0., 0.}, BaselinePosition -> 
           Baseline], "CacheGraphics" -> False], 
         Spacer[{0., 0.}]], TraditionalForm], {0.005, 0.}, {
        AbsoluteThickness[0.1]}}, {-0.35667494393873245`, 
       FormBox[
        InterpretationBox[
         StyleBox[
          
          GraphicsBox[{}, ImageSize -> {0., 0.}, BaselinePosition -> 
           Baseline], "CacheGraphics" -> False], 
         Spacer[{0., 0.}]], TraditionalForm], {0.005, 0.}, {
        AbsoluteThickness[0.1]}}, {-0.2231435513142097, 
       FormBox[
        InterpretationBox[
         StyleBox[
          
          GraphicsBox[{}, ImageSize -> {0., 0.}, BaselinePosition -> 
           Baseline], "CacheGraphics" -> False], 
         Spacer[{0., 0.}]], TraditionalForm], {0.005, 0.}, {
        AbsoluteThickness[0.1]}}, {-0.10536051565782628`, 
       FormBox[
        InterpretationBox[
         StyleBox[
          
          GraphicsBox[{}, ImageSize -> {0., 0.}, BaselinePosition -> 
           Baseline], "CacheGraphics" -> False], 
         Spacer[{0., 0.}]], TraditionalForm], {0.005, 0.}, {
        AbsoluteThickness[0.1]}}, {0.09531017980432493, 
       FormBox[
        InterpretationBox[
         StyleBox[
          
          GraphicsBox[{}, ImageSize -> {0., 0.}, BaselinePosition -> 
           Baseline], "CacheGraphics" -> False], 
         Spacer[{0., 0.}]], TraditionalForm], {0.005, 0.}, {
        AbsoluteThickness[0.1]}}, {0.1823215567939546, 
       FormBox[
        InterpretationBox[
         StyleBox[
          
          GraphicsBox[{}, ImageSize -> {0., 0.}, BaselinePosition -> 
           Baseline], "CacheGraphics" -> False], 
         Spacer[{0., 0.}]], TraditionalForm], {0.005, 0.}, {
        AbsoluteThickness[0.1]}}, {0.26236426446749106`, 
       FormBox[
        InterpretationBox[
         StyleBox[
          
          GraphicsBox[{}, ImageSize -> {0., 0.}, BaselinePosition -> 
           Baseline], "CacheGraphics" -> False], 
         Spacer[{0., 0.}]], TraditionalForm], {0.005, 0.}, {
        AbsoluteThickness[0.1]}}, {0.3364722366212129, 
       FormBox[
        InterpretationBox[
         StyleBox[
          
          GraphicsBox[{}, ImageSize -> {0., 0.}, BaselinePosition -> 
           Baseline], "CacheGraphics" -> False], 
         Spacer[{0., 0.}]], TraditionalForm], {0.005, 0.}, {
        AbsoluteThickness[0.1]}}, {0.4054651081081644, 
       FormBox[
        InterpretationBox[
         StyleBox[
          
          GraphicsBox[{}, ImageSize -> {0., 0.}, BaselinePosition -> 
           Baseline], "CacheGraphics" -> False], 
         Spacer[{0., 0.}]], TraditionalForm], {0.005, 0.}, {
        AbsoluteThickness[0.1]}}}}]]], "Output",
 CellChangeTimes->{{3.7594740822916*^9, 3.7594741585249605`*^9}, {
   3.759474304568039*^9, 3.7594743292414503`*^9}, 3.759474361653304*^9, 
   3.7594758279318414`*^9, 3.759475901993658*^9, 3.75947892050706*^9, 
   3.7607026398427973`*^9}]
}, Open  ]]
}, Open  ]]
},
WindowSize->{1370, 729},
WindowMargins->{{-8, Automatic}, {Automatic, 0}},
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, 285, 6, 70, "Section"],
Cell[CellGroupData[{
Cell[890, 32, 2023, 58, 112, "Input"],
Cell[2916, 92, 1634, 38, 243, "Output"]
}, Open  ]],
Cell[4565, 133, 1555, 44, 69, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[6157, 182, 209, 3, 70, "Section"],
Cell[CellGroupData[{
Cell[6391, 189, 4695, 118, 500, "Input"],
Cell[CellGroupData[{
Cell[11111, 311, 98, 1, 22, "Print"],
Cell[11212, 314, 100, 1, 22, "Print"],
Cell[11315, 317, 98, 1, 22, "Print"],
Cell[11416, 320, 97, 1, 22, "Print"],
Cell[11516, 323, 100, 1, 22, "Print"],
Cell[11619, 326, 98, 1, 22, "Print"],
Cell[11720, 329, 100, 1, 22, "Print"],
Cell[11823, 332, 100, 1, 22, "Print"],
Cell[11926, 335, 98, 1, 22, "Print"],
Cell[12027, 338, 98, 1, 22, "Print"],
Cell[12128, 341, 100, 1, 22, "Print"],
Cell[12231, 344, 100, 1, 22, "Print"],
Cell[12334, 347, 100, 1, 22, "Print"]
}, Open  ]],
Cell[12449, 351, 1061, 17, 30, "Output"],
Cell[13513, 370, 947, 13, 87, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[14497, 388, 205, 4, 30, "Input"],
Cell[14705, 394, 6711, 162, 243, "Output"]
}, Open  ]]
}, Open  ]]
}
]
*)