(* 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[     47441,       1144]
NotebookOptionsPosition[     46189,       1097]
NotebookOutlinePosition[     46537,       1112]
CellTagsIndexPosition[     46494,       1109]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["Problem 1: The 1D linear heat equation", "Section",
 CellChangeTimes->{{3.759309696794531*^9, 3.759309706357973*^9}, {
  3.759402411221094*^9, 3.7594024142211156`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"heatTransfer", "[", 
   RowBox[{"T_", ",", "l_", ",", "h_", ",", "d_"}], "]"}], ":=", 
  RowBox[{"(", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Tau]", "=", 
     RowBox[{
      RowBox[{"h", "^", "2"}], " ", "/", 
      RowBox[{"(", 
       RowBox[{"3", "d"}], ")"}]}]}], ";", "\[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"}], "]"}], "]"}], "=", "0"}], ";"}]}], 
     "\[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[{"\[Tau]", "+", 
            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"}]}], "]"}], "]"}], "=", 
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{"1", "-", 
            RowBox[{"2", "d", " ", 
             RowBox[{"\[Tau]", "/", 
              RowBox[{"h", "^", "2"}]}]}]}], ")"}], 
          RowBox[{"y", "[", 
           RowBox[{"[", 
            RowBox[{
             RowBox[{"n", "+", "1"}], ",", "j"}], "]"}], "]"}]}], "+", 
         RowBox[{"2", "d", " ", 
          RowBox[{"\[Tau]", "/", 
           RowBox[{"h", "^", "2"}]}], " ", 
          RowBox[{"y", "[", 
           RowBox[{"[", 
            RowBox[{"n", ",", "j"}], "]"}], "]"}]}], "+", "\[Tau]"}]}], 
       ";"}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", 
    RowBox[{"{", 
     RowBox[{"x", ",", "t", ",", "y"}], "}"}]}], "\[IndentingNewLine]", 
   ")"}]}]], "Input",
 CellChangeTimes->{{3.75930971766958*^9, 3.7593097526652794`*^9}, {
   3.7593107568896537`*^9, 3.7593107857953224`*^9}, {3.7593110158758144`*^9, 
   3.7593111867137566`*^9}, {3.759311383621805*^9, 3.759311461374673*^9}, {
   3.759311586839445*^9, 3.759311873329012*^9}, 3.7593122935935974`*^9, 
   3.7593126473744545`*^9, 3.759335337913406*^9, {3.7593423031290984`*^9, 
   3.759342309279076*^9}, 3.7594724962937155`*^9, {3.759472533746035*^9, 
   3.7594725669490366`*^9}, {3.7594726289858265`*^9, 
   3.7594726292669983`*^9}, {3.759472789262472*^9, 3.7594727897152605`*^9}, {
   3.759473407479626*^9, 3.7594734122139792`*^9}, {3.759473520070476*^9, 
   3.759473546842258*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"{", 
    RowBox[{"x", ",", "time", ",", "solution"}], "}"}], "=", 
   RowBox[{"heatTransfer", "[", 
    RowBox[{"1", ",", "1", ",", "0.1", ",", "1"}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.759334895315104*^9, 3.759334896349883*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Manipulate", "[", 
  RowBox[{
   RowBox[{"Quiet", "[", "\[IndentingNewLine]", 
    RowBox[{"ListPlot", "[", 
     RowBox[{
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{"x", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "1"}], "]"}], "]"}], ",", 
          RowBox[{"solution", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "j"}], "]"}], "]"}]}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", 
          RowBox[{"Length", "[", "x", "]"}]}], "}"}]}], "]"}], ",", 
      RowBox[{"PlotRange", "\[Rule]", " ", 
       RowBox[{"{", 
        RowBox[{"0", ",", "1"}], "}"}]}]}], "]"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"j", ",", "1", ",", 
     RowBox[{
      RowBox[{"Dimensions", "[", "solution", "]"}], "[", 
      RowBox[{"[", "2", "]"}], "]"}], ",", "1"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.7593128117950544`*^9, 3.7593128258842545`*^9}, {
   3.7593128900092087`*^9, 3.7593129133939*^9}, {3.759313301386424*^9, 
   3.7593133653026752`*^9}, {3.7593142191051555`*^9, 
   3.7593143078133802`*^9}, {3.759314492099165*^9, 3.7593145217656527`*^9}, {
   3.7593146479175887`*^9, 3.7593146578497086`*^9}, {3.7593147247041407`*^9, 
   3.7593147289637933`*^9}, {3.759314791654887*^9, 3.759314796749386*^9}, {
   3.759375303874225*^9, 3.759375312904493*^9}, 3.7594003316932087`*^9, {
   3.7594726399861975`*^9, 3.7594726484036083`*^9}, {3.7594726984193416`*^9, 
   3.7594727343399143`*^9}, {3.7594727955920725`*^9, 3.75947279574846*^9}, {
   3.7594728586431985`*^9, 3.759472859246338*^9}, {3.7594735522761717`*^9, 
   3.759473570246344*^9}, {3.762316055434261*^9, 3.762316055690096*^9}}],

Cell[BoxData[
 TagBox[
  StyleBox[
   DynamicModuleBox[{$CellContext`j$$ = 301, Typeset`show$$ = True, 
    Typeset`bookmarkList$$ = {}, Typeset`bookmarkMode$$ = "Menu", 
    Typeset`animator$$, Typeset`animvar$$ = 1, Typeset`name$$ = 
    "\"untitled\"", Typeset`specs$$ = {{
      Hold[$CellContext`j$$], 1, 301, 1}}, Typeset`size$$ = {
    360., {117., 121.}}, Typeset`update$$ = 0, Typeset`initDone$$, 
    Typeset`skipInitDone$$ = True, $CellContext`j$9575$$ = 0}, 
    DynamicBox[Manipulate`ManipulateBoxes[
     1, StandardForm, "Variables" :> {$CellContext`j$$ = 1}, 
      "ControllerVariables" :> {
        Hold[$CellContext`j$$, $CellContext`j$9575$$, 0]}, 
      "OtherVariables" :> {
       Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$, 
        Typeset`animator$$, Typeset`animvar$$, Typeset`name$$, 
        Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$,
         Typeset`skipInitDone$$}, "Body" :> Quiet[
        ListPlot[
         Table[{
           Part[$CellContext`x, $CellContext`i, 1], 
           
           Part[$CellContext`solution, $CellContext`i, $CellContext`j$$]}, \
{$CellContext`i, 1, 
           Length[$CellContext`x]}], PlotRange -> {0, 1}]], 
      "Specifications" :> {{$CellContext`j$$, 1, 301, 1}}, "Options" :> {}, 
      "DefaultOptions" :> {}],
     ImageSizeCache->{411., {163., 169.}},
     SingleEvaluation->True],
    Deinitialization:>None,
    DynamicModuleValues:>{},
    SynchronousInitialization->True,
    UndoTrackedVariables:>{Typeset`show$$, Typeset`bookmarkMode$$},
    UnsavedVariables:>{Typeset`initDone$$},
    UntrackedVariables:>{Typeset`size$$}], "Manipulate",
   Deployed->True,
   StripOnInput->False],
  Manipulate`InterpretManipulate[1]]], "Output",
 CellChangeTimes->{
  3.759400333082985*^9, {3.7594726340288296`*^9, 3.759472648903335*^9}, {
   3.759472703110239*^9, 3.759472734605303*^9}, {3.7594727929807196`*^9, 
   3.7594727960262165`*^9}, 3.759472859433609*^9, {3.7594735495897856`*^9, 
   3.759473570949327*^9}, 3.7594738559868855`*^9, 3.7610551595483036`*^9, {
   3.7623160581125154`*^9, 3.7623160625626183`*^9}, 3.762316279951446*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Problem 2: The 2D linear heat equation", "Section",
 CellChangeTimes->{{3.7593097103108053`*^9, 3.7593097146073627`*^9}, {
  3.7594024201306353`*^9, 3.7594024227109413`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"heatTransfer", "[", 
   RowBox[{"T_", ",", "l_", ",", "h1_", ",", "h2_", ",", "d_"}], "]"}], ":=", 
  RowBox[{"(", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Tau]", "=", 
     RowBox[{"1", " ", "/", 
      RowBox[{"(", 
       RowBox[{"2", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{"2", " ", "d"}], ")"}], "/", 
           RowBox[{"h1", "^", "2"}]}], "+", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"2", "d"}], ")"}], "/", 
           RowBox[{"(", 
            RowBox[{"h2", "^", "2"}], ")"}]}]}], ")"}]}], ")"}]}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"n", "=", 
     RowBox[{"Ceiling", "[", 
      RowBox[{"l", "/", "h1"}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"m", "=", 
     RowBox[{"Ceiling", "[", 
      RowBox[{"l", "/", "h2"}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"time", "=", 
     RowBox[{"Ceiling", "[", 
      RowBox[{"T", "/", "\[Tau]"}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"t", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"j", "-", "1"}], ")"}], "\[Tau]"}], ",", 
       RowBox[{"{", 
        RowBox[{"j", ",", "1", ",", 
         RowBox[{"time", "+", "1"}]}], "}"}]}], "]"}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"x1", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"i", "-", "1"}], ")"}], "h1"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", 
         RowBox[{"n", "+", "1"}]}], "}"}]}], "]"}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"x2", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"i", "-", "1"}], ")"}], "h2"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", 
         RowBox[{"m", "+", "1"}]}], "}"}]}], "]"}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"y", "=", 
     RowBox[{"Table", "[", 
      RowBox[{"0", ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", 
         RowBox[{"n", "+", "1"}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"j", ",", "1", ",", 
         RowBox[{"m", "+", "1"}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"k", ",", "1", ",", 
         RowBox[{"time", "+", "1"}]}], "}"}]}], "]"}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"(*", 
     RowBox[{"initial", " ", "condition"}], "*)"}], "\[IndentingNewLine]", 
    RowBox[{"For", "[", 
     RowBox[{
      RowBox[{"i", "=", "1"}], ",", 
      RowBox[{"i", "\[LessEqual]", " ", 
       RowBox[{"n", "+", "1"}]}], ",", 
      RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"For", "[", 
        RowBox[{
         RowBox[{"k", "=", "1"}], ",", 
         RowBox[{"k", "\[LessEqual]", " ", 
          RowBox[{"m", "+", "1"}]}], ",", 
         RowBox[{"k", "++"}], ",", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{
           RowBox[{"y", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "k", ",", "1"}], "]"}], "]"}], "=", "0"}], 
          ";"}]}], "\[IndentingNewLine]", "]"}], ";"}]}], 
     "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", 
    "\[IndentingNewLine]", 
    RowBox[{"For", "[", 
     RowBox[{
      RowBox[{"j", "=", "1"}], ",", 
      RowBox[{"j", "\[LessEqual]", " ", "time"}], ",", 
      RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"For", "[", 
        RowBox[{
         RowBox[{"k", "=", "2"}], ",", 
         RowBox[{"k", "\[LessEqual]", " ", "m"}], ",", 
         RowBox[{"k", "++"}], ",", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{"For", "[", 
           RowBox[{
            RowBox[{"i", "=", "2"}], ",", 
            RowBox[{"i", "\[LessEqual]", " ", "n"}], ",", 
            RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
            "\[IndentingNewLine]", 
            RowBox[{
             RowBox[{
              RowBox[{"y", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "k", ",", 
                 RowBox[{"j", "+", "1"}]}], "]"}], "]"}], "=", 
              RowBox[{"\[Tau]", "+", 
               RowBox[{"d", " ", 
                RowBox[{"\[Tau]", "/", 
                 RowBox[{"h1", "^", "2"}]}], 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"y", "[", 
                   RowBox[{"[", 
                    RowBox[{
                    RowBox[{"i", "+", "1"}], ",", "k", ",", "j"}], "]"}], 
                   "]"}], "+", 
                  RowBox[{"y", "[", 
                   RowBox[{"[", 
                    RowBox[{
                    RowBox[{"i", "-", "1"}], ",", "k", ",", "j"}], "]"}], 
                   "]"}]}], ")"}]}], "+", 
               RowBox[{"d", " ", 
                RowBox[{"\[Tau]", "/", 
                 RowBox[{"h2", "^", "2"}]}], 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"y", "[", 
                   RowBox[{"[", 
                    RowBox[{"i", ",", 
                    RowBox[{"k", "+", "1"}], ",", "j"}], "]"}], "]"}], "+", 
                  RowBox[{"y", "[", 
                   RowBox[{"[", 
                    RowBox[{"i", ",", 
                    RowBox[{"k", "-", "1"}], ",", "j"}], "]"}], "]"}]}], 
                 ")"}]}], "+", 
               RowBox[{
                RowBox[{"(", 
                 RowBox[{"1", "-", 
                  RowBox[{"2", " ", "d", " ", 
                   RowBox[{"\[Tau]", "/", 
                    RowBox[{"h1", "^", "2"}]}]}], "-", 
                  RowBox[{"2", "d", " ", 
                   RowBox[{"\[Tau]", "/", 
                    RowBox[{"h2", "^", "2"}]}]}]}], ")"}], " ", 
                RowBox[{"y", "[", 
                 RowBox[{"[", 
                  RowBox[{"i", ",", "k", ",", "j"}], "]"}], "]"}]}]}]}], 
             ";"}]}], "\[IndentingNewLine]", "]"}], ";"}]}], 
        "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", 
       "\[IndentingNewLine]", 
       RowBox[{"For", "[", 
        RowBox[{
         RowBox[{"i", "=", "2"}], ",", 
         RowBox[{"i", "\[LessEqual]", "n"}], ",", 
         RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{
           RowBox[{"y", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "1", ",", 
              RowBox[{"time", "+", "1"}]}], "]"}], "]"}], "=", "0"}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{
           RowBox[{"y", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", 
              RowBox[{"m", "+", "1"}], ",", 
              RowBox[{"time", "+", "1"}]}], "]"}], "]"}], "=", "0"}], ";"}]}],
         "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", 
       "\[IndentingNewLine]", 
       RowBox[{"For", "[", 
        RowBox[{
         RowBox[{"k", "=", "2"}], ",", 
         RowBox[{"k", "\[LessEqual]", "m"}], ",", 
         RowBox[{"k", "++"}], ",", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{
           RowBox[{"y", "[", 
            RowBox[{"[", 
             RowBox[{"1", ",", "k", ",", 
              RowBox[{"time", "+", "1"}]}], "]"}], "]"}], "=", "0"}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{
           RowBox[{"y", "[", 
            RowBox[{"[", 
             RowBox[{
              RowBox[{"n", "+", "1"}], ",", "k", ",", 
              RowBox[{"time", "+", "1"}]}], "]"}], "]"}], "=", "0"}], ";"}]}],
         "\[IndentingNewLine]", "\[IndentingNewLine]", "]"}], ";"}]}], 
     "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", 
    RowBox[{"{", 
     RowBox[{"x1", ",", "x2", ",", "t", ",", "y"}], "}"}]}], 
   "\[IndentingNewLine]", ")"}]}]], "Input",
 CellChangeTimes->{{3.759341407413829*^9, 3.7593414229877653`*^9}, {
   3.759341560612998*^9, 3.759341602362714*^9}, {3.7593417443448906`*^9, 
   3.759341793301532*^9}, {3.759341885131345*^9, 3.7593419659109697`*^9}, {
   3.7593421202397327`*^9, 3.759342145611512*^9}, 3.759342316998917*^9, {
   3.7593425290239897`*^9, 3.7593427428932123`*^9}, {3.7593427827529497`*^9, 
   3.7593428685326424`*^9}, {3.759342958696414*^9, 3.759342965076392*^9}, {
   3.759374741433628*^9, 3.759374747651861*^9}, {3.759375007815615*^9, 
   3.759375040785512*^9}, 3.75937510646729*^9, {3.759375629922736*^9, 
   3.7593756625055895`*^9}, {3.759375726361825*^9, 3.7593757301920004`*^9}, 
   3.759375899574464*^9, 3.7593759594451733`*^9, {3.759475876225706*^9, 
   3.7594759916839266`*^9}, {3.7594764710065527`*^9, 
   3.7594764786070013`*^9}, {3.7623160719504647`*^9, 3.7623160729623985`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"T", "=", "1"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"l", "=", "1"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"h1", "=", "0.1"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"h2", "=", "0.1"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"d", "=", "1"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"solution", "=", 
   RowBox[{"heatTransfer", "[", 
    RowBox[{"T", ",", "l", ",", "h1", ",", "h2", ",", "d"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"x1", "=", 
   RowBox[{"solution", "[", 
    RowBox[{"[", "1", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"x2", "=", 
   RowBox[{"solution", "[", 
    RowBox[{"[", "2", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"t", "=", 
   RowBox[{"solution", "[", 
    RowBox[{"[", "3", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"result", "=", 
   RowBox[{"solution", "[", 
    RowBox[{"[", "4", "]"}], "]"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.759375908651035*^9, 3.759375909061162*^9}, {
  3.7593759411609807`*^9, 3.759375985460458*^9}, {3.7594007418489323`*^9, 
  3.759400808311803*^9}, {3.759476640849617*^9, 3.759476641624323*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"First", " ", "approach", " ", "for", " ", "visualization"}], 
   "*)"}], "\[IndentingNewLine]", 
  RowBox[{"Manipulate", "[", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Quiet", "[", 
     RowBox[{"ListDensityPlot", "[", 
      RowBox[{
       RowBox[{"result", "[", 
        RowBox[{"[", 
         RowBox[{"All", ",", "All", ",", "j"}], "]"}], "]"}], ",", 
       RowBox[{"PlotLegends", "\[Rule]", " ", "True"}]}], "]"}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"j", ",", "1", ",", 
      RowBox[{
       RowBox[{"Dimensions", "[", "result", "]"}], "[", 
       RowBox[{"[", "3", "]"}], "]"}], ",", "1"}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.7594674067303867`*^9, 3.7594674094494205`*^9}, {
  3.759467440137745*^9, 3.759467446346943*^9}, {3.759474829683099*^9, 
  3.759474831580659*^9}, {3.7594749293269606`*^9, 3.7594749678077445`*^9}, {
  3.759476161248836*^9, 3.759476205381722*^9}, {3.7594762357657013`*^9, 
  3.759476248001685*^9}, {3.7594762966847696`*^9, 3.759476335578883*^9}, {
  3.75947638745156*^9, 3.759476388696495*^9}, {3.759476455994372*^9, 
  3.7594764585628624`*^9}, {3.7594766252705717`*^9, 3.759476629167714*^9}, {
  3.7594767012934566`*^9, 3.759476705972399*^9}, {3.75947682838948*^9, 
  3.7594768831758647`*^9}, {3.759477031382639*^9, 3.759477076018319*^9}, {
  3.761055228840185*^9, 3.7610552391363134`*^9}}],

Cell[BoxData[
 TagBox[
  StyleBox[
   DynamicModuleBox[{$CellContext`j$$ = 1, Typeset`show$$ = True, 
    Typeset`bookmarkList$$ = {}, Typeset`bookmarkMode$$ = "Menu", 
    Typeset`animator$$, Typeset`animvar$$ = 1, Typeset`name$$ = 
    "\"untitled\"", Typeset`specs$$ = {{
      Hold[$CellContext`j$$], 1, 801, 1}}, Typeset`size$$ = {
    420., {178., 183.}}, Typeset`update$$ = 0, Typeset`initDone$$, 
    Typeset`skipInitDone$$ = True, $CellContext`j$2815$$ = 0}, 
    DynamicBox[Manipulate`ManipulateBoxes[
     1, StandardForm, "Variables" :> {$CellContext`j$$ = 1}, 
      "ControllerVariables" :> {
        Hold[$CellContext`j$$, $CellContext`j$2815$$, 0]}, 
      "OtherVariables" :> {
       Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$, 
        Typeset`animator$$, Typeset`animvar$$, Typeset`name$$, 
        Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$,
         Typeset`skipInitDone$$}, "Body" :> Quiet[
        ListDensityPlot[
         Part[$CellContext`result, All, All, $CellContext`j$$], PlotLegends -> 
         True]], "Specifications" :> {{$CellContext`j$$, 1, 801, 1}}, 
      "Options" :> {}, "DefaultOptions" :> {}],
     ImageSizeCache->{471., {225., 231.}},
     SingleEvaluation->True],
    Deinitialization:>None,
    DynamicModuleValues:>{},
    SynchronousInitialization->True,
    UndoTrackedVariables:>{Typeset`show$$, Typeset`bookmarkMode$$},
    UnsavedVariables:>{Typeset`initDone$$},
    UntrackedVariables:>{Typeset`size$$}], "Manipulate",
   Deployed->True,
   StripOnInput->False],
  Manipulate`InterpretManipulate[1]]], "Output",
 CellChangeTimes->{{3.7594770764808393`*^9, 3.7594770970076275`*^9}, 
   3.7623160871661487`*^9, 3.76231629103191*^9, 3.7630962514868937`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"Second", " ", "approach", " ", "for", " ", "visualization"}], 
   "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"interpolate", "=", 
     RowBox[{"Interpolation", "[", 
      RowBox[{"Flatten", "[", 
       RowBox[{
        RowBox[{"Table", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"x1", "[", 
             RowBox[{"[", "i", "]"}], "]"}], ",", 
            RowBox[{"x2", "[", 
             RowBox[{"[", "j", "]"}], "]"}], ",", 
            RowBox[{
             RowBox[{"result", "[", 
              RowBox[{"[", 
               RowBox[{"All", ",", "All", ",", "2"}], "]"}], "]"}], "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "j"}], "]"}], "]"}]}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"i", ",", "1", ",", 
            RowBox[{"Length", "[", "x1", "]"}]}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"j", ",", "1", ",", 
            RowBox[{"Length", "[", "x2", "]"}]}], "}"}]}], "]"}], ",", "1"}], 
       "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"Manipulate", "[", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{"heat", "=", 
       RowBox[{"Interpolation", "[", 
        RowBox[{"Flatten", "[", 
         RowBox[{
          RowBox[{"Table", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"x1", "[", 
               RowBox[{"[", "i", "]"}], "]"}], ",", 
              RowBox[{"x2", "[", 
               RowBox[{"[", "k", "]"}], "]"}], ",", 
              RowBox[{
               RowBox[{"result", "[", 
                RowBox[{"[", 
                 RowBox[{"All", ",", "All", ",", "j"}], "]"}], "]"}], "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "k"}], "]"}], "]"}]}], "}"}], ",", 
            RowBox[{"{", 
             RowBox[{"i", ",", "1", ",", 
              RowBox[{"Length", "[", "x1", "]"}]}], "}"}], ",", 
            RowBox[{"{", 
             RowBox[{"k", ",", "1", ",", 
              RowBox[{"Length", "[", "x2", "]"}]}], "}"}]}], "]"}], ",", 
          "1"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"Quiet", "[", 
       RowBox[{"DensityPlot", "[", 
        RowBox[{
         RowBox[{"heat", "[", 
          RowBox[{"xx", ",", "yy"}], "]"}], ",", 
         RowBox[{"{", 
          RowBox[{"xx", ",", "0", ",", "1"}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{"yy", ",", "0", ",", "1"}], "}"}], ",", 
         RowBox[{"PlotLegends", "\[Rule]", " ", "True"}]}], "]"}], "]"}]}], 
     ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "1", ",", 
       RowBox[{
        RowBox[{"Dimensions", "[", "result", "]"}], "[", 
        RowBox[{"[", "3", "]"}], "]"}], ",", "1"}], "}"}]}], 
    "]"}]}]}]], "Input",
 CellChangeTimes->{{3.7593769099260097`*^9, 3.759376910985984*^9}, 
   3.7593770183049784`*^9, {3.7593772363728113`*^9, 3.759377256346032*^9}, 
   3.75937729614389*^9, {3.7593778698406305`*^9, 3.7593779615301046`*^9}, {
   3.7593780342196956`*^9, 3.759378041020095*^9}, {3.7593780791996527`*^9, 
   3.759378101193055*^9}, {3.7594003879727955`*^9, 3.759400388072689*^9}, {
   3.7594766471605024`*^9, 3.7594766477089715`*^9}, {3.759476757920019*^9, 
   3.759476762859883*^9}, {3.759477004713394*^9, 3.759477013861369*^9}, {
   3.761055206982399*^9, 3.761055217137909*^9}, {3.762316098372772*^9, 
   3.762316099292406*^9}}],

Cell[BoxData[
 TagBox[
  StyleBox[
   DynamicModuleBox[{$CellContext`j$$ = 437, Typeset`show$$ = True, 
    Typeset`bookmarkList$$ = {}, Typeset`bookmarkMode$$ = "Menu", 
    Typeset`animator$$, Typeset`animvar$$ = 1, Typeset`name$$ = 
    "\"untitled\"", Typeset`specs$$ = {{
      Hold[$CellContext`j$$], 1, 801, 1}}, Typeset`size$$ = {
    426., {176., 182.}}, Typeset`update$$ = 0, Typeset`initDone$$, 
    Typeset`skipInitDone$$ = True, $CellContext`j$3193$$ = 0}, 
    DynamicBox[Manipulate`ManipulateBoxes[
     1, StandardForm, "Variables" :> {$CellContext`j$$ = 1}, 
      "ControllerVariables" :> {
        Hold[$CellContext`j$$, $CellContext`j$3193$$, 0]}, 
      "OtherVariables" :> {
       Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$, 
        Typeset`animator$$, Typeset`animvar$$, Typeset`name$$, 
        Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$,
         Typeset`skipInitDone$$}, "Body" :> ($CellContext`heat = Interpolation[
          Flatten[
           Table[{
             Part[$CellContext`x1, $CellContext`i], 
             Part[$CellContext`x2, $CellContext`k], 
             Part[
              
              Part[$CellContext`result, All, 
               All, $CellContext`j$$], $CellContext`i, $CellContext`k]}, \
{$CellContext`i, 1, 
             Length[$CellContext`x1]}, {$CellContext`k, 1, 
             Length[$CellContext`x2]}], 1]]; Quiet[
         DensityPlot[
          $CellContext`heat[$CellContext`xx, $CellContext`yy], \
{$CellContext`xx, 0, 1}, {$CellContext`yy, 0, 1}, PlotLegends -> True]]), 
      "Specifications" :> {{$CellContext`j$$, 1, 801, 1}}, "Options" :> {}, 
      "DefaultOptions" :> {}],
     ImageSizeCache->{477., {223., 229.}},
     SingleEvaluation->True],
    Deinitialization:>None,
    DynamicModuleValues:>{},
    SynchronousInitialization->True,
    UndoTrackedVariables:>{Typeset`show$$, Typeset`bookmarkMode$$},
    UnsavedVariables:>{Typeset`initDone$$},
    UntrackedVariables:>{Typeset`size$$}], "Manipulate",
   Deployed->True,
   StripOnInput->False],
  Manipulate`InterpretManipulate[1]]], "Output",
 CellChangeTimes->{{3.759378038203176*^9, 3.7593780416846027`*^9}, {
   3.75937808322948*^9, 3.7593781016893263`*^9}, 3.7594003890928674`*^9, 
   3.759400815323762*^9, 3.7594031029671125`*^9, 3.759475912365425*^9, 
   3.7594765083838816`*^9, 3.7594767637462435`*^9, {3.759477009384736*^9, 
   3.759477014692751*^9}, 3.7610551736566577`*^9, 3.7623161001169634`*^9, 
   3.762316294101735*^9, 3.7630962538300943`*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Problem 3: The logistic function - parameters estimation ", "Section",
 CellChangeTimes->{{3.7594024025511093`*^9, 3.75940244594067*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"y", "=", 
   RowBox[{"{", 
    RowBox[{
    "0.04", ",", "0.35", ",", "1.05", ",", "1.9", ",", "1.98", ",", "1.97", 
     ",", "1.99", ",", "1.999"}], "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"f", "[", "t_", "]"}], ":=", 
   FractionBox[
    SubscriptBox["x", "0"], 
    RowBox[{
     RowBox[{
      RowBox[{"(", 
       RowBox[{"1", "-", 
        FractionBox[
         SubscriptBox["x", "0"], "K"]}], ")"}], 
      RowBox[{"E", "^", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "r"}], " ", "t"}], ")"}]}]}], "+", 
     FractionBox[
      SubscriptBox["x", "0"], "K"]}]]}], ";"}]}], "Input",
 CellChangeTimes->{{3.7594029736370597`*^9, 3.759403081847352*^9}, {
  3.759403234856665*^9, 3.7594032466066093`*^9}, {3.762316389911727*^9, 
  3.7623163901857867`*^9}, {3.7623165468006864`*^9, 3.76231654905097*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"Clear", "[", 
   RowBox[{
    SubscriptBox["x", "0"], ",", "r", ",", "K"}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"error", "=", 
   RowBox[{
    RowBox[{"f", "[", 
     RowBox[{"Range", "[", 
      RowBox[{"0", ",", "7", ",", "1"}], "]"}], "]"}], "-", "y"}]}], ";", 
  RowBox[{"min", "=", 
   RowBox[{"NMinimize", "[", 
    RowBox[{
     RowBox[{"Sum", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"error", "[", 
         RowBox[{"[", "i", "]"}], "]"}], "^", "2"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", 
         RowBox[{"Length", "[", "error", "]"}]}], "}"}]}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["x", "0"], ",", "r", ",", "K"}], "}"}]}], "]"}]}], 
  ";"}], "\n", 
 RowBox[{
  RowBox[{"{", 
   RowBox[{"x0", ",", "r", ",", "K"}], "}"}], "=", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{
     SubscriptBox["x", "0"], "/.", " ", 
     RowBox[{"min", "[", 
      RowBox[{"[", "2", "]"}], "]"}]}], ",", 
    RowBox[{"r", "/.", " ", 
     RowBox[{"min", "[", 
      RowBox[{"[", "2", "]"}], "]"}]}], ",", 
    RowBox[{"K", "/.", " ", 
     RowBox[{"min", "[", 
      RowBox[{"[", "2", "]"}], "]"}]}]}], "}"}]}]}], "Input",
 CellChangeTimes->{{3.759403156356901*^9, 3.7594031823383665`*^9}, 
   3.759403268046238*^9, {3.7623163081618924`*^9, 3.7623163738915205`*^9}, {
   3.762316510791046*^9, 3.7623165173209352`*^9}}],

Cell[BoxData[
 TemplateBox[{
  "Clear","ssym",
   "\"\\!\\(\\*SubscriptBox[\\\"x\\\", \\\"0\\\"]\\) is not a symbol or a \
string.\"",2,30,5,17248772297139696545,"Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{{3.76231635150203*^9, 3.7623163743122034`*^9}, {
  3.762316491451322*^9, 3.7623165504810085`*^9}, {3.763096257719806*^9, 
  3.763096265452365*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "0.04378363518774445`", ",", "2.0105916623258993`", ",", 
   "2.006564273109243`"}], "}"}]], "Output",
 CellChangeTimes->{
  3.759403506998983*^9, 3.759403548314893*^9, 3.761055185327669*^9, {
   3.762316123481413*^9, 3.76231616207134*^9}, 3.76231621129591*^9, {
   3.762316318931595*^9, 3.762316375248827*^9}, {3.762316492244912*^9, 
   3.7623165505247374`*^9}, {3.7630962587039537`*^9, 3.763096265483627*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{
   RowBox[{"ListPlot", "[", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"i", "-", "1"}], ",", 
        RowBox[{"y", "[", 
         RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", 
        RowBox[{"Length", "[", "y", "]"}]}], "}"}]}], "]"}], "]"}], ",", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"f", "[", "t", "]"}], "/.", " ", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"r", "\[Rule]", " ", "r"}], ",", 
        RowBox[{
         SubscriptBox["x", "0"], "\[Rule]", "x0"}], ",", 
        RowBox[{"K", "\[Rule]", " ", "K"}]}], "}"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"t", ",", "0", ",", "8"}], "}"}]}], "]"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", " ", "All"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.7594037715238657`*^9, 3.759403849633113*^9}, {
  3.759403902092827*^9, 3.7594039052932177`*^9}}],

Cell[BoxData[
 GraphicsBox[{{{}, {{}, 
     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
      0.012833333333333334`], AbsoluteThickness[1.6], 
      PointBox[{{0., 0.04}, {1., 0.35}, {2., 1.05}, {3., 1.9}, {4., 1.98}, {
       5., 1.97}, {6., 1.99}, {7., 
       1.999}}]}, {}}, {}, {}, {{}, {}}}, {{{}, {}, 
     TagBox[
      {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], 
       Opacity[1.], LineBox[CompressedData["
1:eJwV1Hk0FFwUAPCRtbEb68wwWSpki6KEd9tElohC2SlSKmtS1kKkqChJKEKW
lJBsIVtlF7LWl7ILFdnne/PXO79z3nnv3vvevZKOF46e2kAgEFaYCATGevjU
WHvVeJjOdhHNK61uGYhlP/niVzFnRPJguXUnPAOppSQkvRbzRgkzBRo+TzOQ
QXNPwWOx62j16uDvJ30ZyML5X8VdsTgUS+c/FG+YiQaVI/bfEktH9zlrN7Fp
ZKFUoSPhUWKFyP2AzmelrdnoSvOD2uOlhcimWl/W81A2Oh72bYP0ySLkERKi
NeCSjbjmPYLKkorRQtrZY5ufZyO/z3F+U7S3yLffVSpWMQeZ3utzM95SiRr+
/hc6o52LQjheb40Uq0XVvF8NxAJfIPlDd24UXKxFZh/+pBdlvEAdYRfG+xtq
Ud5pFdfw1hdIhlkhR8m3Dk1zS8rXS+aj+rV0pc8d9Sg1+U9cXWM+4vpzfyct
+gNSmH9vGhj7EiUM+u8vWm9BS35m6zzrr9BeqlX60LFWZOHyxH9YqABNnNBg
5chrRa83ONt2KxYgianLTnmKbWgxXn9SxLYAXTijvpkw0Iakrbk5eaoKkIBz
ftazXR1osD3u2Ejka2Rp9SR/du4zkkVb843Ui9BMvznxKbkLmbzUMrQ2KULh
thynzPZ3IX6n7RGRbkXotdMFcnFcF3pSz2wFKUWI57x22BWNbjR/j7RBmliM
6kJ7T7AF9CCPwS/uvT+LkWouPxuVow/dEVCTWSwqQR+U6uybVfoQ0SAonPK5
BNm/8isLtOpDu5dzPjr8LkG3ir9e/C+7D613j0dYKL9FY1V5/RlG/Sjni239
ZPZblNyl/2r7vQH0VsNHrierFHHSQ6wPSXxFY7YfrSqqylHCGQtbBd2vaNO5
a6pd38uRzGcFe/7zX1FT0wcuEmsF0nne49Rf8RVdbIy1ZDpcgd4+v1/0Rukb
Gu0sM8vtqkDTZlPlBi3fUEVFnZn130p0LOvBJy/u7ygkakZXOKsKJdmEGp6W
/I7IsbWxGWVV6LuAe7Plzu/IIkRvt2FrFboQsK9V2+Y70hVO+9uzUIUij053
sOd9R377AzX/6FajitV9fY8Mh1G4P5Px4clqtNn019j7mz/QhaXsBLaD7xE5
6JCmy5MfyNdflc/E5j3iy0u9SXzzA/GQlZ8W+bxHyxxmSqbff6D5Y0cDuzLe
o7aqN15Du36irKF7C27EWuSvErK++PMn2qh/c5m5pxa18AoKKsEoCm7gRVnh
9ei9tvup9mOjqF7a0CDraT0qOVtf7H12FGURUmsbKutRWoOfVdn9URRxYoPb
qX/1yC9oMFlvehRlHtK4KOHWgKRmMuScEsdQnr3cXPKxRnSpZTdK+DuOBDJr
WuMPfkQ8SkU/P3BMoPLWttA7Th9R+i2V6BXqBPIJ6w7MCfmI2g239toenECl
XduyNCs/om2fBL02359AbbN+lMZdn9DXhpmMAo1J5EcKNChQa0KHqp/xNPtP
IVnjyDOWe1vQ0CbJovXbU2hUq+WmrlML8g5OOqmSNoX2mkg/OhbWgp6ge1n3
Pk2hyK2lHd0fWtBKRfC+E9RplM1UNORi3oryS0/6jlZMo3iuNcKkRhsSKeQf
YmKaQeX+ZT922rcj6+gV8x/CM4jL3Ut83LsdPXX++aleYQal3JYlv4psRwpC
b9/etJxBTvLZt4Jet6N9vvbxQi9nkH+eXpIIRwdST7LIiFOeRf4nGkR1ijqQ
g22GTdbQLLKVMXvxgvYZFX/f39Sq9Rt9+NBculmjB3XrhczdOvQbGf+zOzNo
0oMWXlQKGx79jchverVz3HqQxpXdDo0uv9Hx0eBzkck9qISkvFB19zdyDAw+
ZcL+Bb09QN5UMPYbjeeyhnANfkHlGbNecXF/0Mv0NyLmD/tQjdtjitXUX2TY
2ip64tgQmt+izZvn/g8Zrv+SLqcMow1LOXyUzctoRi1To3x5FKnT7YvNXqwi
UYJUXuvUFGq6vqOlXp6OPHtfHClfmEW02Mml578I0DdQujKe8gdN/vs61KWy
AUro+dT+mXmUPlH2kec2M5ReIFy1ll1E5x8ckXa7xwzo5mrjaZVFpHFg+Epd
AjMQF/Y9OrtrEX18TFS8msYMncStH5z1FtGciVXMeAkzXKFr0YRcF5HOmwWz
2mFmgGG3a+cyFlHvte2Dl3ezwFl5ObYZiSXEI5458+MHC6xNrjQp0pcQHL82
ODzOAglha8oOrMvIM8bu0/dfLNBvrvouhnMZdTOJZn5bZIH7Hy1PDYkso+SR
G7YDXKzAnzI3uldlGSm9dGvu3MEKdsaW8VL2y8hov1JOzXVWKImY0j5avowC
rxIfVkeygtichMR09TJ6WTQSXnWbFbwHe42vNSwjQdkUx8oEVpjaMZj5pGMZ
DXDxkUtzWeFgV6djxtgyOtc9d+PVZ1aYfZjt0E5aQdFnik6nyrBBheime29P
r6AJyLvBKc8GaXZvrymdW0GHRJ9l+yqzgZz6WbVkjxW0oSHul6EmG7D8PtHn
fnUFXZbx9l06wgZ8W/fFtNxZQae/qoYfvcIGq4+eTzwoW0F7j718xtLBBqIO
4ZM2nKsoWSGr4UIPG9QWBXCo8K2iVebU8b4BNrhrF6yzJriKSgpiFF+NssE2
wzamSIlVpMR3ochmjQ16O8z2nN6+iqifFOuKZdnhjo+liubxVbQAOT9cg9jB
6/G+I7ZJq6g8ffLRpevs0Luwet80dRWFciiYhUdil32eg/RVxNOWW512jx0o
7PGTpNxVtMXuRfJQJjv4736oebNsFR0PeGVp3sYOJduC06p7V1HR2zdNSJID
JktCsu1Ja8hffPGa8RYO0DyeFN8vvIYgZNcem20c4Df6Y+IoeQ016b197r+T
A7ZOrQypS66hHz2l4UX6HNA1Kqf5WWkNCS5UwDZPDnhX/kHaRX8NeavWFgq/
54BvksR/hKtr6PSOTal3GjmgtkYm7UjQGrJUv3qTq4UD2kPfByeGrqE9mjsc
N/RygH+X7ZBM5BrasC+d99cvDohqeXiN6/4aijG97lZH3ggnFVIvyeevoeyL
BzZ5e26Er5Lp/SVDayjJM5Vz5tJGUGG7ZRH33xq65b26cCZgI1i73PU492MN
XfQrbLa7sRHOF2y/JTCxhjSCZa4YJG+EhWzZd+rza6g2hqVb6uNGeGsuPvUf
cR19fVF7s12SCA7xd1n71NZR2zTLlOpWIpRuEZE6oL6OqhQOGsYpECH5B79d
zq51lJpdx2WhQQS945t1PbTXkWNG/a0BQyKktgWT+3XX0cjjxpjRS0S4bsju
LGe1jqajm+6tNRMhRUjFqvbKOhr6xPXXtpMITgFNfF8C1lEL0ci86gsRMt+q
7x8PWkcvbjQLXhsmwq7zSJrl+jq6eL0lnmOJCFoHdSjU6HX090rbA5IMJyjV
Gmv3PFpHK26fH8n5c4Jk9E6v6LfryFzd5o5SECeUjEXzqpeto1ymkXC165yQ
dphYM1i+jmwT/nlo3+aEOb0zHlur1lF1HVnP9CknkBII29Pr11HEJoe/fh85
wWyxeuxA5zoidU8bNpK5IEYx667K5Do698R3XzONC0bffshPnVpHtefoGh0y
XOBHy6fx/lpHvsz80gNKXODUMbc+MruOelV2LM3s4wL2RI3p0IV1lHLT/5nI
WS6I2tMo4M5ER9v2sq+fLuMCr9jLocLCdPQyO2frlyou4KngYr8gQkc7BE1M
9eu5oOWCbmedKB1pjSakb2vngi18/YfPU+jI+Ja8wewIF7jJWcxkb2LMWaOE
ywLc4Ml7tvkznrtvPeJUo89wg/7koNeiNh3p9O+yXjvPDebMwRXiiI5qDgyG
nffmhjy6a/ReoKOPIlt6TYK4QVBmb/X1fXQ0UPEmSCieG1DEcsmSLh0RiP1N
yVXcMAyTZ1KO0NGhp5Iur0R4wMo2uuynHY43gWdLPZUHAtSdI1rs6WjT7ZUf
fZI8cP2dyK5iBzpavtzlyKLAA07ZykahTnSUZxppawE8cCXq+RluFzoiMc0d
W3PlASPDgmuz7nT01a7qgH4pD+guqJb89qejpmN5zLbveGAjd3BI6RU6KjFI
rPGs5YF0kc3HQ67S0V0Nb0hq4YEdln82cwXS0UFeWa2Z7zxQUfU+UDiEjnIq
Y9TiOXkh9u9ZIj2cji6J20n9Z80LzynFC1vu0pGjicrxAQde+OPCOtKCbXyN
KarnNC/cCb9B9L1HR1vG0+aaL/KC/VzfQnUcHX0uHKsqvc4LryaNfE0e0JGK
oaddXC4vzCXT+vSS6GjMPyzp0CovLK6lmOQ9o6OJ49tcnxD4YJNc3mPtDDqy
nyVtoDHzQWdt8p5mbAPpnzvJHHzgYePiPpFJR1KREY/5BPhgbvxHJTWbjtrM
m86tbeaDN/OuHxxe0JHC1DHOHiM+iK1xbfEsxvUwN3wXZcIHTfZuakvY+yv2
eemY8UG07KeRwDd0dOK2cv8zSz6wytfSulFCR5EqG3O8nfhAf1N66d1SOhrx
Lj8scJkPtrvJNEZU0lHqmlSUUTofcDmKsCzW4ftOi+kwZeL4H//rP1NPR29a
eH8XPueDX5lN1H7slpRVK2o+H6T/dR4oa6Cjtb3dcpNv+eBBT4Sr7wd8f3jk
hxutfBCSanqyrQm/N+8sR+0yH6g7Cb6idtJRR0VUhP4aH7jY3Hvmi3333Ga2
Vjof+NSYyLdhC3w8wdzHwg/erzNHQj5jh9WtzfDwA+u9yN3fuuiIbyXxN0Wa
HwKz1cpvfKEjrtGDA14G/JDlZlccPYj/T/y3E0tG/MCbKJPfix194EpvoAk/
UNOIeluG6Ijz6cvuyGP8UCrLa12JTbShtKfY8cPpoUKxia90xNE5W/fJix/K
/SOOqHynI5Z3j/Klk/jho1hEtv8IHfVpHdD7mswPHqEzTDXY+aVT3xKf8AO/
85dFjlGcf7G2gEAmP/geVWqIx36V99WbUMAPs896/2SP0ZFtkrTmYAM/7GiJ
XH4/gd/zck7t/T/8ELM3iL3wFx3d/mdmc3SBH5bLbLRmsZ18Vue5l/B+v6c0
hRk64vEw2hq2zg/5taQ/T7FPuc5Eem0UAL++M73Rs7heFmpHTGgCoHCghHD0
Nx2d3VHWu9FAADY3R11Pn6ejqrau3BIjAeC1nnTsxhZynw1yMRGA846d6+wL
dPTu2eYtdccEIOFewkU3bJJIrGewvQCQ/vgcVPhHR2VLzpz/fATgGnkDPW0R
1/Mdj/aPVAEo3XZy5tQK/u8n5fjupWFf994Ui130b//w3gwB+F6gJVaKbady
OTIlRwBGYpiduFfp6PWT4a4TxQKwKHvV7hX2yesl59s/CYChqMuemTXcz/qO
qZULAjC21VRNikCAHQqSOmhJAOxMJr6oYlfyfOt/tyIAxWk1Z/Zjt3XaiFQT
SDAZ8MDIGfuvjdXt90QS2HAnbnqGreNpcrVRggQWw4MbpZnw/kQdyw5dEggM
Cd/l20AAq4C1+aP6JNC6LNwqgf3drvxepwEJpF+YLitg/5XRbP1sQoKXKx0U
fWyx/J26PSdI4PvV6VYwtuN7hR0D50nQHRWgPYU9P0nmG7lPglvXtB3fMBMg
7vFoQ9JDEsxpr8bXYqsdKQwySyKBzPO+N+3YHgVGM1VPSLCn9Xj9JPa0X3DL
o1wSyLKP8tFYCPCTdSTatIYEg7pdF0Oxu2ivNlZMkyCIR2xJi5UA3u0B1V6z
JAiJkBLXwyZdO3xZ/g8JJBxIqmbYpiPDYw8WSVCUuKp0Brs5V7jRk1kQjP9t
uhSHXbf7arismCAIaS+ajmEXmR1ijjsgCN8iPxaHsRHghkFJGPWQILiq3Q24
jW29X47jmb4gVCZr736AvUGNk7vYWBA8yi7FZGEbC7QKf7EUBEWdtZKP2GNt
x+Wp7oLQ7KvOxc1OALLxadP0eEHYuqzFG4U9fbCnUyFBEOQnrKPvYFdp6x0v
ShSEr6UEwkPs04ryJ+tTBCGml/lTJvZrrl/OY88FgVq7f64W2/CTj59CpSB8
f9VzYw07SC8spXBEEIY3N6m6chCgw9bVWmVcEOhBDr7nsWV8DMRyJwUhbEyp
wBv74xP+e2mzgmDoliAQgi248jj87rIg3IZXlx5iZ+UVu1/kEYKnfgr9jdht
AmN7FNSFYIwvUldqIwGk5D4tZu4SAt6Ns0+3YvugF0XSe4TA/EneogI2+Zy3
MgWEYKLAKFYD27mOIE08LATVW9g8DLEXL4kRx2yE4MyIDd0bW3Lw8Je0MCEY
6OkJq8aOqtmTf+wGPu8zS3A99u9MhXD2m0KwL+GK1yfsGk+eHWdjheDvLTfd
Lsb5HB0xqo+EYH673uMx7Cw1q0PVL4WgITBHmIdIAKUol+KhfiF46UfwMse+
f8HyVuyQEIR7sGZaYtPN9Z33/ScEu5jsuqyx22jbBDJGhABOK0mewvYo+uXu
PicE4vf4fX2wC7/5bF5lE4arM46UeGxNjetxYqrCcOQSPawNu9tRTP3xDmF4
fXGqqxPb6/aLnk0awtD+JU+yBzv35xeyrJYw6I3V5A5iS8QpPlHXFYZOzYY7
E9gb5rrzzE4IwzaU94iZkwCfsuUaboUKgwTRoF0V27W70pU/TBjeF3s924nN
ssGcGB8hDK7+XL67sHWsAo2SooXhkmIqlw72S/bOzuz7wvDIsolXHzvO+eq3
hmxhWLZh87TFtpZoXWLqFIaHRjbvw7F5xDyLyrqEIfR+hkMkdjVJyMPnizBE
/rRevYm9ZePJsbFBYcjgkdl6B3vmz0h365gw3HBNO/oIO/TD+uvH68Jg8imC
NR87y1vp/G55EWATyV7uxD5xoV3uj4IILPNGKXdjc7t5/8xVFoHM7tMOX7A9
7UqtN+0UgQaNreUD2HsO6xqwgwj0reww/ondQrOV7TomAvvkCWPz2PMfb32/
GCICKfMUbiEuAgxwdrmuXhMBI+lXH4SxawypMxHhIuBiz3tdFDu2JXsl+aYI
7D2i+puCLd/RINgcLwI7ZOwLpbBte5kOyWWLADpY3aGMXT/qk/NfhwiYipCi
9LBzZSu2u3eJQPy3dZbD2PfOsJQs9ojAV1f2QANs+8m7tbyDImBh5udmjL38
K39Qe1QEepWoiubYygvjvIkrIjBEu3TSDjuBxdb7qIwoCG7WPOeNve4ksrN9
iyh07BvM8sF2ft82f0ROFNLO/h32xVYJ3X/JSEkU3j4IP3YZu5Eu66+3SxSa
hi9sCcT+t/wnSMdQFIx6l4MjsI/PRd6U9xaF/XLJmx9il5vsN3zuKwqvHc13
JGJLvVzlkr0sCrq8JXsfYf86f+H25kBRSLG+bfkYO2zK/M6mG6KQGpvu+wS7
cJT2QDhJFIrXh3yeY5OGip5uqBWFrjGHmyXYt7kCAt7Ui0LIkn34W2zingOW
5z6IwuIjzuBS7A0JHdzdLaIg89HtQjn2b9MZv+e9ooBehOpUYbfVyR4xmRGF
v59TUhqwo/MerTwmi0GEssnlLmyOAcfuo+JiIKwmYNSNfZ0o/4p9kxgUrhbS
erCvuJSc9tgsBr/nn1Z/wXaT7Go/oCIGYQmHF/ux9eJ5siYPikFNlYDqd2zW
gBDz3R5iQN/j8d80I7+5qaV3XmIQHrbzxi9G/U9bpuj6ikHisy1KM9i1R5Qn
zK6IgXH8vUuzjPulBoPOh4lBCbWI8AebpXFXTtpDMVjsW5r6x8iXNMvEUyMG
p7LtdzBxE6Ag4mRmXK0YSC/n1jDsuFpvSGkQg73kgyYbsGt+Pn4g2yQG8yes
zzBjh5YYKO7vFoOXUegeKzaTbaal34QYWLkk1m/EJjy3zf8uQIamhUuv+LHr
p3rOpwiRoa+VZacAdrSKqZK1KBmKFAffMCxSsj+vS5wMlebeZSRspQbZnEZZ
Mp7vFqVC2NY/f2e80CHDLbaqB2LYbyUjkv3dyOCrdq5nE3bgKYKthjsZ9O6r
6EpiH3h+WfzvBTJoEWwKGW5TOZfk7kMGE/m521LYY8g00T6YDE624toy2MK2
lPu698kwN7EzcCu2Z2L+LYEaMiQdKEpQxH6zmhhRU0uGbSmH/jK8ahse4tlA
hg/JtkeUsCOkbXw6mnC8MuIsythJuUSbuz1koOz75aKCXVd5SkFgmgx7H7uL
qWGLDlM+8YtS4KUBVXYXts1B9rpqMgU+DXv7MPw083elhzgFVA+fqmFY4dyH
gnYpCtwu4D25G3vvvG/iHUUK/EdNjNDEdmPvOMO/nwJ3zFebtbArtt3g4D9P
AfeQSNJebPXPIs13L1KAeiXAhOGXVzPvCHpRwP+LzC2GnzbXk0X9KDCiE8O6
D/vGBVZFWigF6mUGfzN89HWoqeJ9CuRI5r87gD2qGZCoX0kBrcKnMnrY9sNc
dp+qKBDwpdCS4b6bSdJG7ykgFRUWzXDTQFmuSSMFNIac/jBcELhcadFBgRPK
4pX62FdrLg2fGqFAvkXjIQNs/sOeCqE8VPiXS1Eyxk6t5uXU5KeCaoW5LcMq
u/PG50hUUBM1vM2wsexYhoMYFbYIXppmOIrNVnKvDBVK2xqyjmCzvNcXJmhS
YZ9WKskUe2GPJCHwFBXoSt7lZthhryuHdrpSoUie8pNhwW3WFdNuVOitjuI2
x1YjP/C3uUiF0ZlaG4YvLnLNa12hwi6vvBWGxwsXx1diqXCZO0bxOPaAYlvn
5XIqhBwNPmuJnTLJH8D8jgrXTgTfZNjxudmWW9VUsMixyGF4XKbHL7WeCq4L
0RMML1CGxBvaqBDYJOJixagHccpF8CcVNC/1WZ7A1h1hX33BLQ7Mo6c3WWNv
TNd/totPHIgP8vcw3ORw07hGQByYYruOM2w2yPOkS0Qc5LhLoxm2/yykuyop
Do4r4wsMX6mRjtVXF4cgec46G8b/SEYyP2zFYXTJS98O+4//H3YHB3FI991n
z7C6RebkoJM4HA/+68twGS/v6x5XcbBuZE5nuD54aO8nT3F4z9S5ynC/Y4Bd
Qbg40PcefW6Pzba1LDH4hTi4Pyj/5YCtz3whcO2lOOzYpbPBkTEfvko5+r8W
h+ZrmUIM8ydEyXmXiMPruO1aDFM3nihxqRGH+Z9bbjCsOrnYZdwtDgf5X1Kc
GP2Ur84vvi4OT0OYVZ2xi08s8rQRJCDcnnqAYR62Uq5rzBLwsIZ6nOEqa22O
cQ4JyI37fJlhKeIBeiFJAoBVporhEWfTaUM5CXBQstU7hX2efO7DVXMJsEgW
0TvNyLdOsUHZQgL2tNtYMEzzmKn9biUBm61iTjPc1uBZpWcnATU9hdcZVvW5
/EbQTQLG4rmqGP7XGv4sNwjH85FNzQU7OOxJyEC2BBQ6DhBdsbkvsR92yJPA
fftDmOFEV3eBkXwJmHYelGL4tcHutJlCCeiceqLJ8Ah/e+2GdxLwdPrSGYYN
kgkccp0SsJHgWsuw8Bu72z4rEmAo1H3hDHZaVt3x5TUJuGlH8WdYJXEbLYhA
gxOvzK4zrB/wL/8GKw0E8x4nMHx1f0z7I14axHxMqGL4e2ulYI00DTIjdvO4
YeeNiT/iMaTBf0djUxm+8TjMfsiYBsqiVtkMOx+d3vzClAbuJMFChqnl5S+N
LGhwcsm+geGomJP10Y40cE6ynGbYVf3hHPEyDfoKWtXPMt7nuqAeewYNLlw8
WcXw+q6r3D1ZNMhJut3IcO/0cEdGDg28NYvbGI61eG2j+4oGVZzD3ximyx/1
CiunwdOGQsI5Rr+2xyQzd9JA+7yuDsMl4f+cOrvw+TKNBxmO22Mnl/aFBnq6
yIhhg2dKhfuGaBDCwWPNcKlf84eQcRrYSCtfZvi+4s4Y0ykaFHXaBDPs8T3J
XHKGBsXOoREMyxqe+1r1lwYFDnnxDLMwfU6P/UeDf+Ovkxj+VrTHzX6ZBkE3
X6YxXO6WpqyyRoM95mnZDD+gcc7T6TTgPhr9iuH/AWFksEA=
        "]]},
      Annotation[#, "Charting`Private`Tag$5049#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.7594038331834087`*^9, 3.7594038500230904`*^9}, {
   3.759403903747207*^9, 3.7594039056631474`*^9}, 3.761055188327463*^9, 
   3.762316219611923*^9, 3.7623163276648397`*^9, {3.7623165282611027`*^9, 
   3.762316551730707*^9}, 3.7630962670145226`*^9}]
}, Open  ]]
}, Open  ]]
},
WindowSize->{1280, 941},
WindowMargins->{{-8, Automatic}, {Automatic, -8}},
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, 174, 2, 70, "Section"],
Cell[757, 26, 5463, 145, 392, "Input"],
Cell[6223, 173, 289, 7, 30, "Input"],
Cell[CellGroupData[{
Cell[6537, 184, 1732, 37, 50, "Input"],
Cell[8272, 223, 2153, 44, 349, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[10474, 273, 180, 2, 70, "Section"],
Cell[10657, 277, 8652, 220, 715, "Input"],
Cell[19312, 499, 1259, 34, 202, "Input"],
Cell[CellGroupData[{
Cell[20596, 537, 1407, 28, 69, "Input"],
Cell[22006, 567, 1756, 35, 473, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[23799, 607, 3455, 85, 107, "Input"],
Cell[27257, 694, 2537, 51, 469, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[29843, 751, 143, 1, 70, "Section"],
Cell[29989, 754, 891, 26, 78, "Input"],
Cell[CellGroupData[{
Cell[30905, 784, 1436, 44, 69, "Input"],
Cell[32344, 830, 379, 8, 23, "Message"],
Cell[32726, 840, 453, 9, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[33216, 854, 996, 28, 30, "Input"],
Cell[34215, 884, 11946, 209, 243, "Output"]
}, Open  ]]
}, Open  ]]
}
]
*)