(* 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[   3117696,      52046]
NotebookOptionsPosition[   3111850,      51838]
NotebookOutlinePosition[   3112250,      51855]
CellTagsIndexPosition[   3112207,      51852]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["Mathematica tools to study delay differential equations", "Title",
 CellChangeTimes->{{3.6803436809633603`*^9, 3.680343714873103*^9}, 
   3.6803437502124243`*^9}],

Cell[CellGroupData[{

Cell["Project 2: ", "Chapter",
 CellChangeTimes->{{3.6803438221000133`*^9, 3.680344133263152*^9}, {
  3.6803441726486826`*^9, 3.680344186598686*^9}, {3.680357826768152*^9, 
  3.6803578304716682`*^9}, {3.680358148394867*^9, 3.6803581490631094`*^9}}],

Cell[CellGroupData[{

Cell["Goals", "Section",
 CellChangeTimes->{{3.6803579202020936`*^9, 3.6803579436600676`*^9}, 
   3.680357995670724*^9}],

Cell["\<\
Introduction to elementary linear equations with constant delay the basic \
tools to study them\
\>", "Text",
 CellChangeTimes->{{3.680357948078904*^9, 3.680357992435274*^9}, {
  3.680358023497876*^9, 3.680358043527462*^9}, {3.680358079836583*^9, 
  3.6803580989843874`*^9}, {3.680358221239419*^9, 3.680358260155901*^9}, {
  3.680358298199155*^9, 3.680358308821535*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Task", "Section",
 CellChangeTimes->{{3.680358008909094*^9, 3.6803580173942356`*^9}}],

Cell[TextData[{
 "Assume that the Investigate the differential equation of population model \
with delay. ",
 Cell[BoxData[
  FormBox[
   RowBox[{" ", 
    RowBox[{
     RowBox[{
      RowBox[{"p", "'"}], 
      RowBox[{"(", "t", ")"}]}], "=", " ", 
     RowBox[{"-", 
      RowBox[{"p", "(", 
       RowBox[{"t", "-", "\[Tau]"}], ")"}]}]}]}], TraditionalForm]]],
 "."
}], "Text",
 CellChangeTimes->{{3.6803442710960474`*^9, 3.6803442742831683`*^9}, {
  3.6803443489543886`*^9, 3.6803444547939205`*^9}, {3.6803500450177865`*^9, 
  3.6803500559383497`*^9}, {3.680358451517172*^9, 3.6803584575353327`*^9}, {
  3.6803585231388903`*^9, 3.6803585278259573`*^9}, {3.6803586847121887`*^9, 
  3.680358693071994*^9}, {3.6810268064123735`*^9, 3.6810268077740493`*^9}, {
  3.681057983393797*^9, 3.6810580080495505`*^9}, {3.6810582030962543`*^9, 
  3.6810582030962543`*^9}}],

Cell["In Mathematica, it is written as", "Text",
 CellChangeTimes->{{3.6803442710960474`*^9, 3.6803442742831683`*^9}, {
  3.6803443489543886`*^9, 3.6803444547939205`*^9}, {3.6803500450177865`*^9, 
  3.6803500579078913`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"dde", "=", 
  RowBox[{
   RowBox[{
    RowBox[{"p", "'"}], "[", "t", "]"}], "\[Equal]", " ", 
   RowBox[{"-", 
    RowBox[{"p", "[", 
     RowBox[{"t", "-", "\[Tau]"}], "]"}]}]}]}]], "Input",
 CellChangeTimes->{{3.6803443850182686`*^9, 3.6803443950877733`*^9}, {
   3.6803444583121624`*^9, 3.680344462467309*^9}, {3.680358697275686*^9, 
   3.680358700193898*^9}, {3.680587900835311*^9, 3.6805879010823097`*^9}, {
   3.6807983244565153`*^9, 3.680798324841091*^9}, 3.6810580156505475`*^9, 
   3.6810582085290365`*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   SuperscriptBox["p", "\[Prime]",
    MultilineFunction->None], "[", "t", "]"}], "\[Equal]", 
  RowBox[{"-", 
   RowBox[{"p", "[", 
    RowBox[{"t", "-", "\[Tau]"}], "]"}]}]}]], "Output",
 CellChangeTimes->{3.680349898135892*^9, 3.6807983275879717`*^9, 
  3.681058016629506*^9, 3.681058209283497*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Equilibrium", "Section",
 CellChangeTimes->{{3.6803501264183836`*^9, 3.680350137457404*^9}}],

Cell["The concept of equilibria is the same as for ODE\[CloseCurlyQuote]s", \
"Text",
 CellChangeTimes->{{3.6803501384799824`*^9, 3.680350211256872*^9}, {
  3.680354211945733*^9, 3.6803543111930723`*^9}, {3.6803543700256233`*^9, 
  3.6803544032199736`*^9}, {3.6803591843719187`*^9, 3.680359221080619*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["ODE vs. DDE: observations before doing anything", "Section",
 CellChangeTimes->{{3.6803587244898353`*^9, 3.6803587282920046`*^9}, {
  3.680358765361582*^9, 3.680358784475958*^9}, {3.680358884817026*^9, 
  3.6803588947256002`*^9}, {3.680358977848401*^9, 3.6803590229408655`*^9}}],

Cell["Initial functions needed instead of initial values", "Text",
 CellChangeTimes->{{3.68035886817314*^9, 3.6803589508769894`*^9}, {
  3.6803590328500595`*^9, 3.6803590406717043`*^9}, {3.6803591007870708`*^9, 
  3.6803591214876275`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Solve the equation starting from an initial function :", "Section",
 CellChangeTimes->{{3.6803544120408106`*^9, 3.6803544620571012`*^9}, {
  3.68035923427151*^9, 3.680359235726017*^9}, {3.680359923198638*^9, 
  3.6803599251347837`*^9}}],

Cell[CellGroupData[{

Cell["The equation", "Subsection",
 CellChangeTimes->{{3.680359702224255*^9, 3.6803597131001215`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"Clear", "[", "k", "]"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"dde1", "=", 
  RowBox[{
   RowBox[{
    RowBox[{"p", "'"}], "[", "t", "]"}], "\[Equal]", 
   RowBox[{"-", 
    RowBox[{"p", "[", 
     RowBox[{"t", "-", "1"}], "]"}]}]}]}]}], "Input",
 CellChangeTimes->{{3.6803443850182686`*^9, 3.6803443950877733`*^9}, {
   3.6803444583121624`*^9, 3.680344462467309*^9}, {3.680354599297612*^9, 
   3.6803545996447678`*^9}, {3.6803596234358*^9, 3.6803596276640673`*^9}, {
   3.68035968027702*^9, 3.68035968417957*^9}, 3.68105802415473*^9, 
   3.6810582230183487`*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   SuperscriptBox["p", "\[Prime]",
    MultilineFunction->None], "[", "t", "]"}], "\[Equal]", 
  RowBox[{"-", 
   RowBox[{"p", "[", 
    RowBox[{
     RowBox[{"-", "1"}], "+", "t"}], "]"}]}]}]], "Output",
 CellChangeTimes->{3.6803544789295864`*^9, 3.68035463009929*^9, 
  3.680359690714684*^9, 3.6810580259167056`*^9, 3.681058223565312*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["The initial function", "Subsection",
 CellChangeTimes->{{3.6803596928488307`*^9, 3.680359721120282*^9}}],

Cell["Conditional definition is used!!!", "Text",
 CellChangeTimes->{{3.6803597266709075`*^9, 3.680359736852356*^9}, {
  3.6803597921130533`*^9, 3.6803597982300544`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"p", "[", 
   RowBox[{"t", "/;", 
    RowBox[{"t", "<=", "0"}]}], "]"}], "\[Equal]", "2"}]], "Input",
 CellChangeTimes->{{3.68035980692216*^9, 3.6803598073914785`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"p", "[", 
   RowBox[{"t", "/;", 
    RowBox[{"t", "\[LessEqual]", "0"}]}], "]"}], "\[Equal]", "2"}]], "Output",
 CellChangeTimes->{3.681058032216462*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Symbolic solutions", "Subsection",
 CellChangeTimes->{{3.6803631344558516`*^9, 3.680363144932556*^9}}],

Cell[CellGroupData[{

Cell["\[Tau] = 1, initial function= 1", "Subsubsection",
 CellChangeTimes->{{3.6803602955588417`*^9, 3.680360306681379*^9}, {
   3.680361966281682*^9, 3.6803619793520727`*^9}, 3.680798352504949*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"DSolveValue", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"p", "'"}], "[", "t", "]"}], "\[Equal]", 
      RowBox[{"-", 
       RowBox[{"p", "[", 
        RowBox[{"t", "-", "1"}], "]"}]}]}], ",", 
     RowBox[{
      RowBox[{"p", "[", 
       RowBox[{"t", "/;", 
        RowBox[{"t", "<=", "0"}]}], "]"}], "\[Equal]", "1"}]}], "}"}], ",", 
   RowBox[{"p", "[", "t", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"t", ",", "0", ",", "5"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.6803598973533936`*^9, 3.680359907073102*^9}, {
   3.680359961216595*^9, 3.6803599767446127`*^9}, {3.68036003661672*^9, 
   3.68036003729494*^9}, {3.6803600814055743`*^9, 3.680360099700844*^9}, 
   3.6803612852398825`*^9, 3.6810580366212864`*^9, 3.6810582284238896`*^9}],

Cell[BoxData[
 TagBox[GridBox[{
    {"\[Piecewise]", GridBox[{
       {"1", 
        RowBox[{"t", "\[LessEqual]", "0"}]},
       {
        RowBox[{"1", "-", "t"}], 
        RowBox[{"0", "<", "t", "\[LessEqual]", "1"}]},
       {
        RowBox[{
         FractionBox["1", "2"], " ", 
         RowBox[{"(", 
          RowBox[{"3", "-", 
           RowBox[{"4", " ", "t"}], "+", 
           SuperscriptBox["t", "2"]}], ")"}]}], 
        RowBox[{"1", "<", "t", "\[LessEqual]", "2"}]},
       {
        RowBox[{
         FractionBox["1", "6"], " ", 
         RowBox[{"(", 
          RowBox[{"17", "-", 
           RowBox[{"24", " ", "t"}], "+", 
           RowBox[{"9", " ", 
            SuperscriptBox["t", "2"]}], "-", 
           SuperscriptBox["t", "3"]}], ")"}]}], 
        RowBox[{"2", "<", "t", "\[LessEqual]", "3"}]},
       {
        RowBox[{
         FractionBox["1", "24"], " ", 
         RowBox[{"(", 
          RowBox[{"149", "-", 
           RowBox[{"204", " ", "t"}], "+", 
           RowBox[{"90", " ", 
            SuperscriptBox["t", "2"]}], "-", 
           RowBox[{"16", " ", 
            SuperscriptBox["t", "3"]}], "+", 
           SuperscriptBox["t", "4"]}], ")"}]}], 
        RowBox[{"3", "<", "t", "\[LessEqual]", "4"}]},
       {
        RowBox[{
         FractionBox["1", "120"], " ", 
         RowBox[{"(", 
          RowBox[{"1769", "-", 
           RowBox[{"2300", " ", "t"}], "+", 
           RowBox[{"1090", " ", 
            SuperscriptBox["t", "2"]}], "-", 
           RowBox[{"240", " ", 
            SuperscriptBox["t", "3"]}], "+", 
           RowBox[{"25", " ", 
            SuperscriptBox["t", "4"]}], "-", 
           SuperscriptBox["t", "5"]}], ")"}]}], 
        RowBox[{"4", "<", "t", "\[LessEqual]", "5"}]},
       {"Indeterminate", 
        TagBox["True",
         "PiecewiseDefault",
         AutoDelete->True]}
      },
      AllowedDimensions->{2, Automatic},
      Editable->True,
      GridBoxAlignment->{
       "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
        "RowsIndexed" -> {}},
      GridBoxItemSize->{
       "Columns" -> {{Automatic}}, "ColumnsIndexed" -> {}, "Rows" -> {{1.}}, 
        "RowsIndexed" -> {}},
      GridBoxSpacings->{"Columns" -> {
          Offset[0.27999999999999997`], {
           Offset[0.84]}, 
          Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
          Offset[0.2], {
           Offset[0.4]}, 
          Offset[0.2]}, "RowsIndexed" -> {}},
      Selectable->True]}
   },
   GridBoxAlignment->{
    "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
     "RowsIndexed" -> {}},
   GridBoxItemSize->{
    "Columns" -> {{Automatic}}, "ColumnsIndexed" -> {}, "Rows" -> {{1.}}, 
     "RowsIndexed" -> {}},
   GridBoxSpacings->{"Columns" -> {
       Offset[0.27999999999999997`], {
        Offset[0.35]}, 
       Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
       Offset[0.2], {
        Offset[0.4]}, 
       Offset[0.2]}, "RowsIndexed" -> {}}],
  "Piecewise",
  DeleteWithContents->True,
  Editable->False,
  SelectWithContents->True,
  Selectable->False]], "Output",
 CellChangeTimes->{{3.680359965672501*^9, 3.6803599778146496`*^9}, 
   3.6803600383494806`*^9, {3.680360087626328*^9, 3.6803601007713165`*^9}, 
   3.680361285940071*^9, 3.681058039062127*^9, 3.6810582291082215`*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["\[Tau] = \[Pi]/2, initial function= sin(t)", "Subsubsection",
 CellChangeTimes->{{3.6803602955588417`*^9, 3.680360306681379*^9}, {
   3.680361966281682*^9, 3.6803619793520727`*^9}, {3.6803620172960167`*^9, 
   3.6803620393115377`*^9}, 3.6807983559436245`*^9, {3.6810583388695803`*^9, 
   3.6810583392383595`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"DSolveValue", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"p", "'"}], "[", "t", "]"}], "\[Equal]", 
      RowBox[{"-", 
       RowBox[{"p", "[", 
        RowBox[{"t", "-", 
         RowBox[{"Pi", "/", "2"}]}], "]"}]}]}], ",", 
     RowBox[{
      RowBox[{"p", "[", 
       RowBox[{"t", "/;", 
        RowBox[{"t", "<=", "0"}]}], "]"}], "\[Equal]", 
      RowBox[{"Sin", "[", "t", "]"}]}]}], "}"}], ",", 
   RowBox[{"p", "[", "t", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"t", ",", "0", ",", 
     RowBox[{"2", "Pi"}]}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.6803598973533936`*^9, 3.680359907073102*^9}, {
   3.680359961216595*^9, 3.6803599767446127`*^9}, {3.68036003661672*^9, 
   3.68036003729494*^9}, {3.6803600814055743`*^9, 3.680360099700844*^9}, {
   3.6803608652238503`*^9, 3.6803608716162405`*^9}, {3.6803609077200108`*^9, 
   3.6803609098532667`*^9}, {3.680361037095152*^9, 3.6803610376916866`*^9}, {
   3.680361140132572*^9, 3.68036116569905*^9}, {3.681026946936861*^9, 
   3.6810269474413414`*^9}, {3.681058043442069*^9, 3.6810580552942886`*^9}, 
   3.6810582336265373`*^9}],

Cell[BoxData[
 TagBox[GridBox[{
    {"\[Piecewise]", GridBox[{
       {
        RowBox[{"Sin", "[", "t", "]"}], 
        RowBox[{
         RowBox[{"t", "\[LessEqual]", "0"}], "||", 
         RowBox[{"0", "<", "t", "\[LessEqual]", 
          FractionBox["\[Pi]", "2"]}], "||", 
         RowBox[{
          FractionBox["\[Pi]", "2"], "<", "t", "\[LessEqual]", "\[Pi]"}], "||", 
         RowBox[{"\[Pi]", "<", "t", "\[LessEqual]", 
          FractionBox[
           RowBox[{"3", " ", "\[Pi]"}], "2"]}], "||", 
         RowBox[{
          FractionBox[
           RowBox[{"3", " ", "\[Pi]"}], "2"], "<", "t", "\[LessEqual]", 
          RowBox[{"2", " ", "\[Pi]"}]}]}]},
       {"Indeterminate", 
        TagBox["True",
         "PiecewiseDefault",
         AutoDelete->True]}
      },
      AllowedDimensions->{2, Automatic},
      Editable->True,
      GridBoxAlignment->{
       "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
        "RowsIndexed" -> {}},
      GridBoxItemSize->{
       "Columns" -> {{Automatic}}, "ColumnsIndexed" -> {}, "Rows" -> {{1.}}, 
        "RowsIndexed" -> {}},
      GridBoxSpacings->{"Columns" -> {
          Offset[0.27999999999999997`], {
           Offset[0.84]}, 
          Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
          Offset[0.2], {
           Offset[0.4]}, 
          Offset[0.2]}, "RowsIndexed" -> {}},
      Selectable->True]}
   },
   GridBoxAlignment->{
    "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
     "RowsIndexed" -> {}},
   GridBoxItemSize->{
    "Columns" -> {{Automatic}}, "ColumnsIndexed" -> {}, "Rows" -> {{1.}}, 
     "RowsIndexed" -> {}},
   GridBoxSpacings->{"Columns" -> {
       Offset[0.27999999999999997`], {
        Offset[0.35]}, 
       Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
       Offset[0.2], {
        Offset[0.4]}, 
       Offset[0.2]}, "RowsIndexed" -> {}}],
  "Piecewise",
  DeleteWithContents->True,
  Editable->False,
  SelectWithContents->True,
  Selectable->False]], "Output",
 CellChangeTimes->{
  3.6803611672164717`*^9, 3.6810269505802717`*^9, {3.681058044312237*^9, 
   3.681058056532352*^9}, 3.6810582342590666`*^9, 3.6810582678234544`*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["\[Tau] = \[Pi], initial function= sin(t)", "Subsubsection",
 CellChangeTimes->{{3.6803602955588417`*^9, 3.680360306681379*^9}, {
   3.680361966281682*^9, 3.6803619793520727`*^9}, {3.6803620172960167`*^9, 
   3.6803620393115377`*^9}, 3.6807983559436245`*^9, {3.681058334418941*^9, 
   3.681058342872754*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"DSolveValue", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"p", "'"}], "[", "t", "]"}], "\[Equal]", 
      RowBox[{"-", 
       RowBox[{"p", "[", 
        RowBox[{"t", "-", "Pi"}], "]"}]}]}], ",", 
     RowBox[{
      RowBox[{"p", "[", 
       RowBox[{"t", "/;", 
        RowBox[{"t", "<=", "0"}]}], "]"}], "\[Equal]", 
      RowBox[{"Sin", "[", "t", "]"}]}]}], "}"}], ",", 
   RowBox[{"p", "[", "t", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"t", ",", "0", ",", 
     RowBox[{"2", "Pi"}]}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.6803598973533936`*^9, 3.680359907073102*^9}, {
   3.680359961216595*^9, 3.6803599767446127`*^9}, {3.68036003661672*^9, 
   3.68036003729494*^9}, {3.6803600814055743`*^9, 3.680360099700844*^9}, {
   3.6803608652238503`*^9, 3.6803608716162405`*^9}, {3.6803609077200108`*^9, 
   3.6803609098532667`*^9}, {3.680361037095152*^9, 3.6803610376916866`*^9}, {
   3.680361140132572*^9, 3.68036116569905*^9}, {3.681026946936861*^9, 
   3.6810269474413414`*^9}, {3.681058043442069*^9, 3.6810580552942886`*^9}, 
   3.6810582336265373`*^9, {3.6810583466772656`*^9, 3.6810583469304895`*^9}}],

Cell[BoxData[
 TagBox[GridBox[{
    {"\[Piecewise]", GridBox[{
       {
        RowBox[{"Sin", "[", "t", "]"}], 
        RowBox[{"t", "\[LessEqual]", "0"}]},
       {
        RowBox[{"1", "-", 
         RowBox[{"Cos", "[", "t", "]"}]}], 
        RowBox[{"0", "<", "t", "\[LessEqual]", "\[Pi]"}]},
       {
        RowBox[{"2", "+", "\[Pi]", "-", "t", "-", 
         RowBox[{"Sin", "[", "t", "]"}]}], 
        RowBox[{"\[Pi]", "<", "t", "\[LessEqual]", 
         RowBox[{"2", " ", "\[Pi]"}]}]},
       {"Indeterminate", 
        TagBox["True",
         "PiecewiseDefault",
         AutoDelete->True]}
      },
      AllowedDimensions->{2, Automatic},
      Editable->True,
      GridBoxAlignment->{
       "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
        "RowsIndexed" -> {}},
      GridBoxItemSize->{
       "Columns" -> {{Automatic}}, "ColumnsIndexed" -> {}, "Rows" -> {{1.}}, 
        "RowsIndexed" -> {}},
      GridBoxSpacings->{"Columns" -> {
          Offset[0.27999999999999997`], {
           Offset[0.84]}, 
          Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
          Offset[0.2], {
           Offset[0.4]}, 
          Offset[0.2]}, "RowsIndexed" -> {}},
      Selectable->True]}
   },
   GridBoxAlignment->{
    "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
     "RowsIndexed" -> {}},
   GridBoxItemSize->{
    "Columns" -> {{Automatic}}, "ColumnsIndexed" -> {}, "Rows" -> {{1.}}, 
     "RowsIndexed" -> {}},
   GridBoxSpacings->{"Columns" -> {
       Offset[0.27999999999999997`], {
        Offset[0.35]}, 
       Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
       Offset[0.2], {
        Offset[0.4]}, 
       Offset[0.2]}, "RowsIndexed" -> {}}],
  "Piecewise",
  DeleteWithContents->True,
  Editable->False,
  SelectWithContents->True,
  Selectable->False]], "Output",
 CellChangeTimes->{
  3.6803611672164717`*^9, 3.6810269505802717`*^9, {3.681058044312237*^9, 
   3.681058056532352*^9}, 3.6810582342590666`*^9, 3.6810582678234544`*^9, 
   3.6810583478998575`*^9}]
}, Open  ]]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["NumericalNDSolve, NDSolveValue)", "Subsection",
 CellChangeTimes->{{3.680354483285797*^9, 3.680354493894696*^9}, {
  3.6803545633893867`*^9, 3.68035456972301*^9}, {3.6803623507009544`*^9, 
  3.6803623692740774`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"NDSolveValue", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"p", "'"}], "[", "t", "]"}], "\[Equal]", 
      RowBox[{"-", 
       RowBox[{"p", "[", 
        RowBox[{"t", "-", "1"}], "]"}]}]}], ",", 
     RowBox[{
      RowBox[{"p", "[", 
       RowBox[{"t", "/;", 
        RowBox[{"t", "<=", "0"}]}], "]"}], "\[Equal]", "1"}]}], "}"}], ",", 
   RowBox[{"p", "[", "t", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"t", ",", "0", ",", "5"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.6803598973533936`*^9, 3.680359907073102*^9}, {
   3.680359961216595*^9, 3.6803599767446127`*^9}, {3.68036003661672*^9, 
   3.68036003729494*^9}, {3.6803600814055743`*^9, 3.680360099700844*^9}, 
   3.6803612852398825`*^9, 3.6803624110495253`*^9, 3.6810580689419837`*^9, 
   3.6810582774303913`*^9}],

Cell[BoxData[
 RowBox[{
  TagBox[
   TemplateBox[{RowBox[{
       StyleBox[
        TagBox["InterpolatingFunction", "SummaryHead"], 
        "NonInterpretableSummary"], 
       StyleBox["[", "NonInterpretableSummary"], 
       DynamicModuleBox[{Typeset`open$$ = False}, 
        PanelBox[
         PaneSelectorBox[{False -> GridBox[{{
               PaneBox[
                ButtonBox[
                 DynamicBox[
                  FEPrivate`FrontEndResource[
                  "FEBitmaps", "SquarePlusIconMedium"]], 
                 ButtonFunction :> (Typeset`open$$ = True), Appearance -> 
                 None, Evaluator -> Automatic, Method -> "Preemptive"], 
                Alignment -> {Center, Center}, ImageSize -> 
                Dynamic[{
                  Automatic, 
                   3.5 (CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
                    Magnification])}]], 
               GraphicsBox[{{{{}, {}, {
                    Directive[
                    Opacity[1.], 
                    RGBColor[0.368417, 0.506779, 0.709798], 
                    AbsoluteThickness[1]], 
                    LineBox[CompressedData["
1:eJwB0QMu/CFib1JlAgAAADwAAAACAAAAStFeOa+/dj6NoYDS///vP1t+433D
3bQ/NZBDkEdk7T8/lvq2a77FP3BaQRJlkOo/UPQI4wkn0D/XhXsOe+znPzg0
cm9ZVNU/5OVGSNNV5T/3vGVz/fHaP4QhTUYBh+I/r0FwQdgX4D+ifB99T9Df
P07J8gTc7uI/ZG0a9kci2j9JHOSKnbjlP27HN+rEjtQ/GA6ZlmZS6D+gx5ul
ZbbOP1IkE95ZJOs/uG6zh5huwz9g2VCrVMbtP/40eaVazbE/5ax+nYYt8D8+
SbdW7KKGv1B/NwP4k/E/1HlEzsQAuL8lIdIrbeLyP0aRw/bq/sS/MFVPcvdM
9D9wxhi3A8jNv+nuA5rgsPU/fmdw9CK30r8MWJqEzfz2P/UrQCPN2NW/ZVMT
jc9k+D9EFVfAqMTYvygebljVtPk/vtD7w30M278he6tB8CD7P8SX2bOkCN2/
xz0gDGqG/D9HRAyml33ev9fPdpnn0/0/fxYJwABp378d9K9Eej3/P3bQRyGG
7d+/53NlWYhHAEDq/udxH/bfv5WgDgED7QBAJhYaVYSW379eFqk3iKABQBK2
/csbxN6/3HO0Tw9IAkBIVOI0E6Xdv3Uasfag/QJAa4cNAsYV3L/lc0kO4q8D
QKnbxL9AP9q/CrVSByVWBEAU62Up20vYv0o/TY9yCgVAL0sTnKz21b8/sbj4
wbIFQIFNqGuZntO/C9a/0sBXBkAUGlFxYTPRv/JDuDvKCgdA5WUW9qP1zL+O
mSGG1bEHQHSawaNbxse/RTh8X+tmCEDS6zyWqR7Cv7G+RxoDEAlANEovwzbD
ub/0965FyrUJQFNS8dgHZq+/UnoHAJxpCkCoYYaijpeUv2Xk0JtvEQtADITG
EwX9kD+Tl4vGTccLQJ5vcLVupqs/l/3hYdt5DEBUa/ftkYW2P1BLqd5qIA1A
RQmQRRLuvT8k4mHqBNUNQFc9/O2CkMI/rWCL16B9DkBqSunxgILFPw2SUDXs
Ig9AluG+gtr9xz+IDAciQtYPQAZ/Ut/zN8o/XDcX+Mw+EEDZdz4lz9rLPwKN
oyZ+mRBAIZ/8J+YlzT8TvH0Nh/IQQAndgaS18M0//l4Q5ZBFEUAK08EnAkTO
P3emGwSgnxFAV3jSmdAvzj/KYd8TsPMRQCeeL0LXus0/q8Eba8VOEkDH5SmL
mtjMP/f6pXoyqBJAYCQAyUmeyz8eqOh6oPsSQNmnKwzwLso/0vmjwhNWE0B+
SYnCclnIP2G/F/uHqhNAwqk3Z3Jpxj8yFFD6//8TQAhGqeRDRMQ/lF/epw==

                    "]]}}}, {}}, {
                DisplayFunction -> Identity, PlotRangePadding -> {{
                    Scaled[0.1], 
                    Scaled[0.1]}, {
                    Scaled[0.1], 
                    Scaled[0.1]}}, PlotRangeClipping -> True, ImagePadding -> 
                 All, DisplayFunction -> Identity, AspectRatio -> 1, 
                 Axes -> {False, False}, AxesLabel -> {None, None}, 
                 AxesOrigin -> {0, 0}, DisplayFunction :> Identity, 
                 Frame -> {{True, True}, {True, True}}, 
                 FrameLabel -> {{None, None}, {None, None}}, FrameStyle -> 
                 Directive[
                   Opacity[0.5], 
                   Thickness[Tiny], 
                   RGBColor[0.368417, 0.506779, 0.709798]], 
                 FrameTicks -> {{None, None}, {None, None}}, 
                 GridLines -> {None, None}, GridLinesStyle -> Directive[
                   GrayLevel[0.5, 0.4]], ImageSize -> 
                 Dynamic[{
                   Automatic, 
                    3.5 (CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
                    Magnification])}], 
                 Method -> {
                  "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> 
                   AbsolutePointSize[6], "ScalingFunctions" -> None}, 
                 PlotRange -> {{0., 5.}, {-0.49939714549689873`, 
                  0.9999999152542373}}, PlotRangeClipping -> True, 
                 PlotRangePadding -> {{
                    Scaled[0.1], 
                    Scaled[0.1]}, {
                    Scaled[0.1], 
                    Scaled[0.1]}}, Ticks -> {Automatic, Automatic}}], 
               GridBox[{{
                  RowBox[{
                    TagBox["\"Domain: \"", "SummaryItemAnnotation"], 
                    "\[InvisibleSpace]", 
                    TagBox[
                    RowBox[{"{", 
                    RowBox[{"{", 
                    RowBox[{"0.`", ",", "5.`"}], "}"}], "}"}], 
                    "SummaryItem"]}]}, {
                  RowBox[{
                    TagBox["\"Output: \"", "SummaryItemAnnotation"], 
                    "\[InvisibleSpace]", 
                    TagBox["\"scalar\"", "SummaryItem"]}]}}, 
                GridBoxAlignment -> {
                 "Columns" -> {{Left}}, "Rows" -> {{Automatic}}}, AutoDelete -> 
                False, GridBoxItemSize -> {
                 "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, 
                GridBoxSpacings -> {
                 "Columns" -> {{2}}, "Rows" -> {{Automatic}}}, 
                BaseStyle -> {
                 ShowStringCharacters -> False, NumberMarks -> False, 
                  PrintPrecision -> 3, ShowSyntaxStyles -> False}]}}, 
             GridBoxAlignment -> {"Rows" -> {{Top}}}, AutoDelete -> False, 
             GridBoxItemSize -> {
              "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, 
             BaselinePosition -> {1, 1}], True -> GridBox[{{
               PaneBox[
                ButtonBox[
                 DynamicBox[
                  FEPrivate`FrontEndResource[
                  "FEBitmaps", "SquareMinusIconMedium"]], 
                 ButtonFunction :> (Typeset`open$$ = False), Appearance -> 
                 None, Evaluator -> Automatic, Method -> "Preemptive"], 
                Alignment -> {Center, Center}, ImageSize -> 
                Dynamic[{
                  Automatic, 
                   3.5 (CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
                    Magnification])}]], 
               GraphicsBox[{{{{}, {}, {
                    Directive[
                    Opacity[1.], 
                    RGBColor[0.368417, 0.506779, 0.709798], 
                    AbsoluteThickness[1]], 
                    LineBox[CompressedData["
1:eJwB0QMu/CFib1JlAgAAADwAAAACAAAAStFeOa+/dj6NoYDS///vP1t+433D
3bQ/NZBDkEdk7T8/lvq2a77FP3BaQRJlkOo/UPQI4wkn0D/XhXsOe+znPzg0
cm9ZVNU/5OVGSNNV5T/3vGVz/fHaP4QhTUYBh+I/r0FwQdgX4D+ifB99T9Df
P07J8gTc7uI/ZG0a9kci2j9JHOSKnbjlP27HN+rEjtQ/GA6ZlmZS6D+gx5ul
ZbbOP1IkE95ZJOs/uG6zh5huwz9g2VCrVMbtP/40eaVazbE/5ax+nYYt8D8+
SbdW7KKGv1B/NwP4k/E/1HlEzsQAuL8lIdIrbeLyP0aRw/bq/sS/MFVPcvdM
9D9wxhi3A8jNv+nuA5rgsPU/fmdw9CK30r8MWJqEzfz2P/UrQCPN2NW/ZVMT
jc9k+D9EFVfAqMTYvygebljVtPk/vtD7w30M278he6tB8CD7P8SX2bOkCN2/
xz0gDGqG/D9HRAyml33ev9fPdpnn0/0/fxYJwABp378d9K9Eej3/P3bQRyGG
7d+/53NlWYhHAEDq/udxH/bfv5WgDgED7QBAJhYaVYSW379eFqk3iKABQBK2
/csbxN6/3HO0Tw9IAkBIVOI0E6Xdv3Uasfag/QJAa4cNAsYV3L/lc0kO4q8D
QKnbxL9AP9q/CrVSByVWBEAU62Up20vYv0o/TY9yCgVAL0sTnKz21b8/sbj4
wbIFQIFNqGuZntO/C9a/0sBXBkAUGlFxYTPRv/JDuDvKCgdA5WUW9qP1zL+O
mSGG1bEHQHSawaNbxse/RTh8X+tmCEDS6zyWqR7Cv7G+RxoDEAlANEovwzbD
ub/0965FyrUJQFNS8dgHZq+/UnoHAJxpCkCoYYaijpeUv2Xk0JtvEQtADITG
EwX9kD+Tl4vGTccLQJ5vcLVupqs/l/3hYdt5DEBUa/ftkYW2P1BLqd5qIA1A
RQmQRRLuvT8k4mHqBNUNQFc9/O2CkMI/rWCL16B9DkBqSunxgILFPw2SUDXs
Ig9AluG+gtr9xz+IDAciQtYPQAZ/Ut/zN8o/XDcX+Mw+EEDZdz4lz9rLPwKN
oyZ+mRBAIZ/8J+YlzT8TvH0Nh/IQQAndgaS18M0//l4Q5ZBFEUAK08EnAkTO
P3emGwSgnxFAV3jSmdAvzj/KYd8TsPMRQCeeL0LXus0/q8Eba8VOEkDH5SmL
mtjMP/f6pXoyqBJAYCQAyUmeyz8eqOh6oPsSQNmnKwzwLso/0vmjwhNWE0B+
SYnCclnIP2G/F/uHqhNAwqk3Z3Jpxj8yFFD6//8TQAhGqeRDRMQ/lF/epw==

                    "]]}}}, {}}, {
                DisplayFunction -> Identity, PlotRangePadding -> {{
                    Scaled[0.1], 
                    Scaled[0.1]}, {
                    Scaled[0.1], 
                    Scaled[0.1]}}, PlotRangeClipping -> True, ImagePadding -> 
                 All, DisplayFunction -> Identity, AspectRatio -> 1, 
                 Axes -> {False, False}, AxesLabel -> {None, None}, 
                 AxesOrigin -> {0, 0}, DisplayFunction :> Identity, 
                 Frame -> {{True, True}, {True, True}}, 
                 FrameLabel -> {{None, None}, {None, None}}, FrameStyle -> 
                 Directive[
                   Opacity[0.5], 
                   Thickness[Tiny], 
                   RGBColor[0.368417, 0.506779, 0.709798]], 
                 FrameTicks -> {{None, None}, {None, None}}, 
                 GridLines -> {None, None}, GridLinesStyle -> Directive[
                   GrayLevel[0.5, 0.4]], ImageSize -> 
                 Dynamic[{
                   Automatic, 
                    3.5 (CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
                    Magnification])}], 
                 Method -> {
                  "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> 
                   AbsolutePointSize[6], "ScalingFunctions" -> None}, 
                 PlotRange -> {{0., 5.}, {-0.49939714549689873`, 
                  0.9999999152542373}}, PlotRangeClipping -> True, 
                 PlotRangePadding -> {{
                    Scaled[0.1], 
                    Scaled[0.1]}, {
                    Scaled[0.1], 
                    Scaled[0.1]}}, Ticks -> {Automatic, Automatic}}], 
               GridBox[{{
                  RowBox[{
                    TagBox["\"Domain: \"", "SummaryItemAnnotation"], 
                    "\[InvisibleSpace]", 
                    TagBox[
                    RowBox[{"{", 
                    RowBox[{"{", 
                    RowBox[{"0.`", ",", "5.`"}], "}"}], "}"}], 
                    "SummaryItem"]}]}, {
                  RowBox[{
                    TagBox["\"Output: \"", "SummaryItemAnnotation"], 
                    "\[InvisibleSpace]", 
                    TagBox["\"scalar\"", "SummaryItem"]}]}, {
                  RowBox[{
                    TagBox["\"Order: \"", "SummaryItemAnnotation"], 
                    "\[InvisibleSpace]", 
                    TagBox["3", "SummaryItem"]}]}, {
                  RowBox[{
                    TagBox["\"Method: \"", "SummaryItemAnnotation"], 
                    "\[InvisibleSpace]", 
                    TagBox["\"Cubic Hermite\"", "SummaryItem"]}]}, {
                  RowBox[{
                    TagBox["\"Periodic: \"", "SummaryItemAnnotation"], 
                    "\[InvisibleSpace]", 
                    TagBox["False", "SummaryItem"]}]}}, 
                GridBoxAlignment -> {
                 "Columns" -> {{Left}}, "Rows" -> {{Automatic}}}, AutoDelete -> 
                False, GridBoxItemSize -> {
                 "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, 
                GridBoxSpacings -> {
                 "Columns" -> {{2}}, "Rows" -> {{Automatic}}}, 
                BaseStyle -> {
                 ShowStringCharacters -> False, NumberMarks -> False, 
                  PrintPrecision -> 3, ShowSyntaxStyles -> False}]}}, 
             GridBoxAlignment -> {"Rows" -> {{Top}}}, AutoDelete -> False, 
             GridBoxItemSize -> {
              "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, 
             BaselinePosition -> {1, 1}]}, 
          Dynamic[Typeset`open$$], ImageSize -> Automatic], BaselinePosition -> 
         Baseline], DynamicModuleValues :> {}], 
       StyleBox["]", "NonInterpretableSummary"]}]},
    "CopyTag",
    DisplayFunction->(#& ),
    InterpretationFunction->("InterpolatingFunction[{{0., 5.}}, <>]"& )],
   False,
   Editable->False,
   SelectWithContents->True,
   Selectable->False], "[", "t", "]"}]], "Output",
 CellChangeTimes->{3.680362412687416*^9, 3.681058076381729*^9, 
  3.6810582782772713`*^9, 3.6810583571111174`*^9}]
}, Open  ]]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Visualization 1: Dependence on the delay", "Section",
 CellChangeTimes->{{3.68035690544081*^9, 3.680356927980993*^9}, {
  3.6803601842416887`*^9, 3.6803601959493046`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Manipulate", "[", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"DynamicModule", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"p", ",", "ddesol"}], "}"}], ",", 
     RowBox[{
      RowBox[{
       RowBox[{"ddesol", "[", "t_", "]"}], ":=", 
       RowBox[{"NDSolveValue", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{
            RowBox[{
             RowBox[{"p", "'"}], "[", "t", "]"}], "\[Equal]", 
            RowBox[{"-", 
             RowBox[{"p", "[", 
              RowBox[{"t", "-", "\[Tau]"}], "]"}]}]}], ",", 
           RowBox[{
            RowBox[{"p", "[", 
             RowBox[{"t", "/;", 
              RowBox[{"t", "<=", "0"}]}], "]"}], "\[Equal]", "1"}]}], "}"}], 
         ",", 
         RowBox[{"p", "[", "t", "]"}], ",", 
         RowBox[{"{", 
          RowBox[{"t", ",", 
           RowBox[{"-", "\[Tau]"}], ",", "T"}], "}"}]}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"Plot", "[", 
       RowBox[{
        RowBox[{"Evaluate", "[", 
         RowBox[{"ddesol", "[", "t", "]"}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"t", ",", 
          RowBox[{"-", "\[Tau]"}], ",", "T"}], "}"}], ",", 
        RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]}]}], 
    "\[IndentingNewLine]", "]"}], ",", "\[IndentingNewLine]", 
   RowBox[{"{", 
    RowBox[{"\[Tau]", ",", "0", ",", "6"}], "}"}], ",", "\[IndentingNewLine]", 
   RowBox[{"{", 
    RowBox[{"T", ",", "1", ",", "20"}], "}"}]}], "\[IndentingNewLine]", 
  "]"}]], "Input",
 CellChangeTimes->{{3.680357254315432*^9, 3.680357422411628*^9}, {
   3.680357498154586*^9, 3.680357546547616*^9}, {3.680360233612874*^9, 
   3.6803602507150707`*^9}, {3.6803604173059483`*^9, 
   3.6803605489615545`*^9}, {3.6803605802914715`*^9, 
   3.6803605813606644`*^9}, {3.680361062668417*^9, 3.68036108114921*^9}, 
   3.6803611298789377`*^9, {3.680362910159089*^9, 3.6803629145617795`*^9}, {
   3.6807983627476025`*^9, 3.680798372932252*^9}, 3.681058081199033*^9, 
   3.681058292221365*^9}],

Cell[BoxData[
 TagBox[
  StyleBox[
   DynamicModuleBox[{$CellContext`T$$ = 12.78, $CellContext`\[Tau]$$ = 0., 
    Typeset`show$$ = True, Typeset`bookmarkList$$ = {}, 
    Typeset`bookmarkMode$$ = "Menu", Typeset`animator$$, Typeset`animvar$$ = 
    1, Typeset`name$$ = "\"untitled\"", Typeset`specs$$ = {{
      Hold[$CellContext`\[Tau]$$], 0, 6}, {
      Hold[$CellContext`T$$], 1, 20}}, Typeset`size$$ = {648., {194., 203.}}, 
    Typeset`update$$ = 0, Typeset`initDone$$, Typeset`skipInitDone$$ = 
    True, $CellContext`\[Tau]$60088$$ = 0, $CellContext`T$60089$$ = 0}, 
    DynamicBox[Manipulate`ManipulateBoxes[
     1, StandardForm, 
      "Variables" :> {$CellContext`T$$ = 1, $CellContext`\[Tau]$$ = 0}, 
      "ControllerVariables" :> {
        Hold[$CellContext`\[Tau]$$, $CellContext`\[Tau]$60088$$, 0], 
        Hold[$CellContext`T$$, $CellContext`T$60089$$, 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" :> 
      DynamicModule[{$CellContext`p, $CellContext`ddesol}, $CellContext`ddesol[
           Pattern[$CellContext`t$, 
            Blank[]]] := 
         NDSolveValue[{
           Derivative[
             1][$CellContext`p][$CellContext`t$] == \
-$CellContext`p[$CellContext`t$ - $CellContext`\[Tau]$$], $CellContext`p[
              Condition[$CellContext`t$, $CellContext`t$ <= 0]] == 1}, 
           $CellContext`p[$CellContext`t$], {$CellContext`t$, -$CellContext`\
\[Tau]$$, $CellContext`T$$}]; Plot[
          Evaluate[
           $CellContext`ddesol[$CellContext`t]], {$CellContext`t, \
-$CellContext`\[Tau]$$, $CellContext`T$$}, PlotRange -> All]], 
      "Specifications" :> {{$CellContext`\[Tau]$$, 0, 6}, {$CellContext`T$$, 
         1, 20}}, "Options" :> {}, "DefaultOptions" :> {}],
     ImageSizeCache->{730., {309., 320.}},
     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.6803611311801453`*^9, 3.68036291501567*^9, 
  3.6807983737663126`*^9, 3.681058083032409*^9, 3.681058133782401*^9, 
  3.6810582928994813`*^9, 3.681058384461447*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Visualization 1I: Dependence on the initial function", "Section",
 CellChangeTimes->{{3.68035690544081*^9, 3.680356927980993*^9}, {
   3.6803601842416887`*^9, 3.6803601959493046`*^9}, {3.680362485722199*^9, 
   3.680362489008337*^9}, {3.680362733777731*^9, 3.6803627356943054`*^9}, 
   3.6803633328722396`*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Manipulate", "[", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"DynamicModule", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"p", ",", "ddesol"}], "}"}], ",", 
     RowBox[{
      RowBox[{
       RowBox[{"ddesol", "[", "t_", "]"}], ":=", 
       RowBox[{"NDSolveValue", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{
            RowBox[{
             RowBox[{"p", "'"}], "[", "t", "]"}], "\[Equal]", 
            RowBox[{"-", 
             RowBox[{"p", "[", 
              RowBox[{"t", "-", "\[Tau]"}], "]"}]}]}], ",", 
           RowBox[{
            RowBox[{"p", "[", 
             RowBox[{"t", "/;", 
              RowBox[{"t", "<=", "0"}]}], "]"}], "\[Equal]", "ifunc"}]}], 
          "}"}], ",", 
         RowBox[{"p", "[", "t", "]"}], ",", 
         RowBox[{"{", 
          RowBox[{"t", ",", 
           RowBox[{"-", "\[Tau]"}], ",", "T"}], "}"}]}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"Plot", "[", 
       RowBox[{
        RowBox[{"Evaluate", "[", 
         RowBox[{"ddesol", "[", "t", "]"}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"t", ",", 
          RowBox[{"-", "\[Tau]"}], ",", "T"}], "}"}], ",", 
        RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]}]}], 
    "\[IndentingNewLine]", "]"}], ",", "\[IndentingNewLine]", 
   RowBox[{"{", 
    RowBox[{"\[Tau]", ",", "0", ",", "6"}], "}"}], ",", "\[IndentingNewLine]", 
   RowBox[{"{", 
    RowBox[{"ifunc", ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"UnitStep", "[", "t", "]"}], ",", "1", ",", 
       RowBox[{"Cos", "[", "t", "]"}]}], "}"}]}], "}"}], ",", 
   "\[IndentingNewLine]", 
   RowBox[{"{", 
    RowBox[{"T", ",", "1", ",", "20"}], "}"}]}], "\[IndentingNewLine]", 
  "]"}]], "Input",
 CellChangeTimes->{{3.680357254315432*^9, 3.680357422411628*^9}, {
   3.680357498154586*^9, 3.680357546547616*^9}, {3.680360233612874*^9, 
   3.6803602507150707`*^9}, {3.6803604173059483`*^9, 
   3.6803605489615545`*^9}, {3.6803605802914715`*^9, 
   3.6803605813606644`*^9}, {3.680361062668417*^9, 3.68036108114921*^9}, 
   3.6803611298789377`*^9, {3.68036274458715*^9, 3.680362799654569*^9}, {
   3.6803628367741537`*^9, 3.6803628689530897`*^9}, 3.6803629621526423`*^9, {
   3.6807984007555904`*^9, 3.680798417729432*^9}, 3.6810581407402067`*^9, 
   3.6810584231748643`*^9}],

Cell[BoxData[
 TagBox[
  StyleBox[
   DynamicModuleBox[{$CellContext`ifunc$$ = 
    UnitStep[$CellContext`t], $CellContext`T$$ = 1, $CellContext`\[Tau]$$ = 0,
     Typeset`show$$ = True, Typeset`bookmarkList$$ = {}, 
    Typeset`bookmarkMode$$ = "Menu", Typeset`animator$$, Typeset`animvar$$ = 
    1, Typeset`name$$ = "\"untitled\"", Typeset`specs$$ = {{
      Hold[$CellContext`\[Tau]$$], 0, 6}, {
      Hold[$CellContext`ifunc$$], {
       UnitStep[$CellContext`t], 1, 
       Cos[$CellContext`t]}}, {
      Hold[$CellContext`T$$], 1, 20}}, Typeset`size$$ = {648., {194., 203.}}, 
    Typeset`update$$ = 0, Typeset`initDone$$, Typeset`skipInitDone$$ = 
    True, $CellContext`\[Tau]$68491$$ = 0, $CellContext`ifunc$68492$$ = 
    0, $CellContext`T$68493$$ = 0}, 
    DynamicBox[Manipulate`ManipulateBoxes[
     1, StandardForm, 
      "Variables" :> {$CellContext`ifunc$$ = 
        UnitStep[$CellContext`t], $CellContext`T$$ = 1, $CellContext`\[Tau]$$ = 
        0}, "ControllerVariables" :> {
        Hold[$CellContext`\[Tau]$$, $CellContext`\[Tau]$68491$$, 0], 
        Hold[$CellContext`ifunc$$, $CellContext`ifunc$68492$$, 0], 
        Hold[$CellContext`T$$, $CellContext`T$68493$$, 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" :> 
      DynamicModule[{$CellContext`p, $CellContext`ddesol}, $CellContext`ddesol[
           Pattern[$CellContext`t$, 
            Blank[]]] := 
         NDSolveValue[{
           Derivative[
             1][$CellContext`p][$CellContext`t$] == \
-$CellContext`p[$CellContext`t$ - $CellContext`\[Tau]$$], $CellContext`p[
              
              Condition[$CellContext`t$, $CellContext`t$ <= 
               0]] == $CellContext`ifunc$$}, 
           $CellContext`p[$CellContext`t$], {$CellContext`t$, -$CellContext`\
\[Tau]$$, $CellContext`T$$}]; Plot[
          Evaluate[
           $CellContext`ddesol[$CellContext`t]], {$CellContext`t, \
-$CellContext`\[Tau]$$, $CellContext`T$$}, PlotRange -> All]], 
      "Specifications" :> {{$CellContext`\[Tau]$$, 0, 
         6}, {$CellContext`ifunc$$, {
          UnitStep[$CellContext`t], 1, 
          Cos[$CellContext`t]}}, {$CellContext`T$$, 1, 20}}, "Options" :> {}, 
      "DefaultOptions" :> {}],
     ImageSizeCache->{730., {309., 320.}},
     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.6810270700685654`*^9, 3.6810271061670284`*^9, {3.681058399341807*^9, 
   3.681058423737665*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["ODE vs. DDE: more observations", "Section",
 CellChangeTimes->{{3.6803587244898353`*^9, 3.6803587282920046`*^9}, {
  3.680358765361582*^9, 3.680358784475958*^9}, {3.680358884817026*^9, 
  3.6803588947256002`*^9}, {3.680358977848401*^9, 3.6803590229408655`*^9}, {
  3.6803591259608994`*^9, 3.6803591312480264`*^9}}],

Cell["\<\
Solutions form very complex semiflow of infinite dimension
Solutions may not be differentiable, may not be continued backwards
 \
\>", "Text",
 CellChangeTimes->{{3.68035886817314*^9, 3.6803589508769894`*^9}, {
  3.6803590328500595`*^9, 3.6803590406717043`*^9}, {3.6803631887550755`*^9, 
  3.680363189624318*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
 "An important observation, time transformation: ",
 Cell[BoxData[
  FormBox[
   RowBox[{"t", "=", "s\[Tau]"}], TraditionalForm]]]
}], "Section",
 CellChangeTimes->{{3.681194325747265*^9, 3.68119440651764*^9}}],

Cell[BoxData[
 FormBox[
  RowBox[{
   RowBox[{
    FractionBox["d", "ds"], 
    RowBox[{"p", "(", 
     RowBox[{"s", " ", "\[Tau]"}], ")"}]}], "\[LongEqual]", 
   RowBox[{"\[Tau]", " ", 
    RowBox[{
     SuperscriptBox["p", "\[Prime]",
      MultilineFunction->None], "(", 
     RowBox[{"s", " ", "\[Tau]"}], ")"}]}], "\[LongEqual]", 
   RowBox[{"k", " ", "\[Tau]", " ", 
    RowBox[{"p", "(", 
     RowBox[{
      RowBox[{"s", " ", "\[Tau]"}], "-", "\[Tau]"}], ")"}]}], "\[LongEqual]", 
   RowBox[{"k", " ", "\[Tau]", " ", 
    RowBox[{"p", "(", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"s", "-", "1"}], ")"}], " ", "\[Tau]"}], ")"}]}]}], 
  TraditionalForm]], "Text",
 CellChangeTimes->{{3.6811934849679556`*^9, 3.6811934867225122`*^9}, {
   3.681193527075882*^9, 3.6811936389340053`*^9}, {3.6811936797409697`*^9, 
   3.6811936988558645`*^9}, {3.681193771805807*^9, 3.6811937801270113`*^9}, 
   3.6811938946426897`*^9, {3.681194445794169*^9, 3.6811944570339875`*^9}}],

Cell["Introduce q(s) := p(s \[Tau])", "Text",
 CellChangeTimes->{{3.681193783879429*^9, 3.681193807841154*^9}, 
   3.681193911418721*^9, 3.6811944457962093`*^9}],

Cell[BoxData[
 FormBox[
  RowBox[{
   FractionBox[
    RowBox[{"d", " ", 
     RowBox[{"q", "(", "s", ")"}]}], "ds"], "\[LongEqual]", 
   RowBox[{"k", " ", "\[Tau]", " ", 
    RowBox[{"q", "(", 
     RowBox[{"s", "-", "1"}], ")"}]}]}], TraditionalForm]], "Text",
 CellChangeTimes->{{3.681193818599145*^9, 3.6811938570523295`*^9}, 
   3.681193911418721*^9, 3.6811944457962093`*^9}],

Cell[CellGroupData[{

Cell["Dependence of q(s) on the delay", "Subsection",
 CellChangeTimes->{{3.68035690544081*^9, 3.680356927980993*^9}, {
  3.6803601842416887`*^9, 3.6803601959493046`*^9}, {3.681194631734655*^9, 
  3.6811946410451183`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Manipulate", "[", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"DynamicModule", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"q", ",", "ddesol"}], "}"}], ",", 
     RowBox[{
      RowBox[{
       RowBox[{"ddesol", "[", "s_", "]"}], ":=", 
       RowBox[{"NDSolveValue", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{
            RowBox[{
             RowBox[{"q", "'"}], "[", "s", "]"}], "\[Equal]", 
            RowBox[{
             RowBox[{"-", "\[Tau]"}], " ", 
             RowBox[{"q", "[", 
              RowBox[{"s", "-", "1"}], "]"}]}]}], ",", 
           RowBox[{
            RowBox[{"q", "[", 
             RowBox[{"s", "/;", 
              RowBox[{"s", "<=", "0"}]}], "]"}], "\[Equal]", "1"}]}], "}"}], 
         ",", 
         RowBox[{"q", "[", "s", "]"}], ",", 
         RowBox[{"{", 
          RowBox[{"s", ",", 
           RowBox[{"-", "1"}], ",", "S"}], "}"}]}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"Plot", "[", 
       RowBox[{
        RowBox[{"Evaluate", "[", 
         RowBox[{"ddesol", "[", "s", "]"}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"s", ",", 
          RowBox[{"-", "1"}], ",", "S"}], "}"}], ",", 
        RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]}]}], 
    "\[IndentingNewLine]", "]"}], ",", "\[IndentingNewLine]", 
   RowBox[{"{", 
    RowBox[{"\[Tau]", ",", "0", ",", "6"}], "}"}], ",", "\[IndentingNewLine]", 
   RowBox[{"{", 
    RowBox[{"S", ",", "1", ",", "20"}], "}"}]}], "\[IndentingNewLine]", 
  "]"}]], "Input",
 CellChangeTimes->{{3.680357254315432*^9, 3.680357422411628*^9}, {
   3.680357498154586*^9, 3.680357546547616*^9}, {3.680360233612874*^9, 
   3.6803602507150707`*^9}, {3.6803604173059483`*^9, 
   3.6803605489615545`*^9}, {3.6803605802914715`*^9, 
   3.6803605813606644`*^9}, {3.680361062668417*^9, 3.68036108114921*^9}, 
   3.6803611298789377`*^9, {3.680362910159089*^9, 3.6803629145617795`*^9}, {
   3.6807983627476025`*^9, 3.680798372932252*^9}, 3.681058081199033*^9, 
   3.681058292221365*^9, {3.681194644623201*^9, 3.681194784255621*^9}}],

Cell[BoxData[
 TagBox[
  StyleBox[
   DynamicModuleBox[{$CellContext`S$$ = 10.38, $CellContext`\[Tau]$$ = 0., 
    Typeset`show$$ = True, Typeset`bookmarkList$$ = {}, 
    Typeset`bookmarkMode$$ = "Menu", Typeset`animator$$, Typeset`animvar$$ = 
    1, Typeset`name$$ = "\"untitled\"", Typeset`specs$$ = {{
      Hold[$CellContext`\[Tau]$$], 0, 6}, {
      Hold[$CellContext`S$$], 1, 20}}, Typeset`size$$ = {648., {199., 208.}}, 
    Typeset`update$$ = 0, Typeset`initDone$$, Typeset`skipInitDone$$ = 
    True, $CellContext`\[Tau]$38054$$ = 0, $CellContext`S$38055$$ = 0}, 
    DynamicBox[Manipulate`ManipulateBoxes[
     1, StandardForm, 
      "Variables" :> {$CellContext`S$$ = 1, $CellContext`\[Tau]$$ = 0}, 
      "ControllerVariables" :> {
        Hold[$CellContext`\[Tau]$$, $CellContext`\[Tau]$38054$$, 0], 
        Hold[$CellContext`S$$, $CellContext`S$38055$$, 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" :> 
      DynamicModule[{$CellContext`q, $CellContext`ddesol}, $CellContext`ddesol[
           Pattern[$CellContext`s$, 
            Blank[]]] := 
         NDSolveValue[{
           Derivative[
             1][$CellContext`q][$CellContext`s$] == (-$CellContext`\[Tau]$$) \
$CellContext`q[$CellContext`s$ - 1], $CellContext`q[
              Condition[$CellContext`s$, $CellContext`s$ <= 0]] == 1}, 
           $CellContext`q[$CellContext`s$], {$CellContext`s$, -1, \
$CellContext`S$$}]; Plot[
          Evaluate[
           $CellContext`ddesol[$CellContext`s]], {$CellContext`s, -1, \
$CellContext`S$$}, PlotRange -> All]], 
      "Specifications" :> {{$CellContext`\[Tau]$$, 0, 6}, {$CellContext`S$$, 
         1, 20}}, "Options" :> {}, "DefaultOptions" :> {}],
     ImageSizeCache->{730., {291., 302.}},
     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.6803611311801453`*^9, 3.68036291501567*^9, 3.6807983737663126`*^9, 
   3.681058083032409*^9, 3.681058133782401*^9, 3.6810582928994813`*^9, 
   3.681058384461447*^9, {3.681194718860942*^9, 3.6811947254336195`*^9}, {
   3.681194758550372*^9, 3.6811947858714285`*^9}}]
}, {2}]]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
The structure of the solutions of linear equations with delay\
\>", "Section",
 CellChangeTimes->{{3.680363783018197*^9, 3.6803638257275133`*^9}}],

Cell[BoxData[
 FormBox[
  RowBox[{
  "Try", " ", "to", " ", "find", " ", "the", " ", "equations", " ", "in", " ",
    "the", " ", "form", " ", 
   RowBox[{
    SuperscriptBox["\[ExponentialE]", 
     RowBox[{"\[Lambda]", " ", "t"}]], ".", " ", "Substitute"}], " ", "into", 
   " ", "the", " ", 
   RowBox[{"equation", ":"}]}], TraditionalForm]], "Text",
 CellChangeTimes->{{3.6803639020543585`*^9, 3.6803639924210625`*^9}, {
  3.6805881890257673`*^9, 3.680588200620825*^9}}],

Cell[CellGroupData[{

Cell[BoxData["dde1"], "Input",
 CellChangeTimes->{{3.6803639981076736`*^9, 3.680364005517214*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   SuperscriptBox["p", "\[Prime]",
    MultilineFunction->None], "[", "t", "]"}], "\[Equal]", 
  RowBox[{"-", 
   RowBox[{"p", "[", 
    RowBox[{
     RowBox[{"-", "1"}], "+", "t"}], "]"}]}]}]], "Output",
 CellChangeTimes->{3.68036400589564*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"p1", "[", "t_", "]"}], "=", 
   SuperscriptBox["\[ExponentialE]", 
    RowBox[{"\[Lambda]", " ", "t"}]]}], ";"}]], "Input",
 CellChangeTimes->{{3.680364032174918*^9, 3.680364047134936*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"(", 
  RowBox[{
   RowBox[{
    RowBox[{"p1", "'"}], "[", "t", "]"}], "\[Equal]", 
   RowBox[{"-", 
    RowBox[{"p1", "[", 
     RowBox[{"t", "-", "1"}], "]"}]}]}], ")"}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   SuperscriptBox["\[ExponentialE]", 
    RowBox[{"t", " ", "\[Lambda]"}]], " ", "\[Lambda]"}], "\[Equal]", 
  RowBox[{"-", 
   SuperscriptBox["\[ExponentialE]", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "1"}], "+", "t"}], ")"}], " ", 
     "\[Lambda]"}]]}]}]], "Output",
 CellChangeTimes->{3.6805882989629946`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Simplify", "[", 
  RowBox[{"Thread", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"p1", "'"}], "[", "t", "]"}], "\[Equal]", 
       RowBox[{"-", 
        RowBox[{"p1", "[", 
         RowBox[{"t", "-", "1"}], "]"}]}]}], ")"}], "/", 
     RowBox[{"p1", "[", "t", "]"}]}], ",", "Equal"}], "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.6805883242949743`*^9, 3.680588326596652*^9}, {
  3.6805884322034807`*^9, 3.6805885246614494`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   SuperscriptBox["\[ExponentialE]", 
    RowBox[{"-", "\[Lambda]"}]], "+", "\[Lambda]"}], "\[Equal]", 
  "0"}]], "Output",
 CellChangeTimes->{{3.680588513890236*^9, 3.680588525784316*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"cheqn", "=", 
  RowBox[{"Simplify", "[", 
   RowBox[{"Thread", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{"p1", "'"}], "[", "t", "]"}], "\[Equal]", 
        RowBox[{"-", 
         RowBox[{"p1", "[", 
          RowBox[{"t", "-", "1"}], "]"}]}]}], ")"}], "/", 
      RowBox[{"p1", "[", "t", "]"}]}], ",", "Equal"}], "]"}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.680364051458361*^9, 3.6803641350394707`*^9}, {
  3.680364485551271*^9, 3.6803645186763525`*^9}, {3.6803645487479706`*^9, 
  3.6803646337351265`*^9}, {3.6803646839318485`*^9, 3.6803648784891567`*^9}, {
  3.6803649192123604`*^9, 3.6803649224980917`*^9}, {3.680588545855546*^9, 
  3.680588551259328*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   SuperscriptBox["\[ExponentialE]", 
    RowBox[{"-", "\[Lambda]"}]], "+", "\[Lambda]"}], "\[Equal]", 
  "0"}]], "Output",
 CellChangeTimes->{{3.6803640715488267`*^9, 3.680364135808674*^9}, 
   3.6803644452568417`*^9, 3.680364487069353*^9, 3.680364519823472*^9, 
   3.6803645514043255`*^9, {3.680364582532942*^9, 3.6803646344139333`*^9}, {
   3.680364688455904*^9, 3.68036475263435*^9}, {3.680364787150921*^9, 
   3.680364806780751*^9}, {3.6803648515092373`*^9, 3.680364879341175*^9}, 
   3.680364924382612*^9, 3.6805885957252913`*^9}]
}, Open  ]],

Cell["\<\
This is the characteristic equation. The roots are the characteristic roots\
\>", "Text",
 CellChangeTimes->{{3.680364886125332*^9, 3.680364913924396*^9}, {
  3.680365003109096*^9, 3.6803650277995515`*^9}}],

Cell["DO not use Solve, since :", "Text",
 CellChangeTimes->{{3.680365074909512*^9, 3.680365101228593*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Solve", "[", 
  RowBox[{"cheqn", ",", "\[Lambda]"}], "]"}]], "Input",
 CellChangeTimes->{{3.6803651119325385`*^9, 3.68036511502139*^9}}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Solve", "::", "ifun"}], "MessageName"], ":", 
  " ", "\<\"Inverse functions are being used by \[NoBreak]\\!\\(\\*RowBox[{\\\
\"Solve\\\"}]\\)\[NoBreak], so some solutions may not be found; use Reduce \
for complete solution information. \
\\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", \
ButtonFrame->None, ButtonData:>\\\"paclet:ref/message/Solve/ifun\\\", \
ButtonNote -> \\\"Solve::ifun\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{3.6803651162884097`*^9, 3.6805886011995907`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"{", 
   RowBox[{"\[Lambda]", "\[Rule]", 
    RowBox[{"ProductLog", "[", 
     RowBox[{"-", "1"}], "]"}]}], "}"}], "}"}]], "Output",
 CellChangeTimes->{3.6803651162894106`*^9, 3.6805886011995907`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Reduce", "[", 
  RowBox[{"cheqn", ",", "\[Lambda]"}], "]"}]], "Input",
 CellChangeTimes->{{3.6803650315542855`*^9, 3.6803650593644743`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"C", "[", "1", "]"}], "\[Element]", "Integers"}], "&&", 
  RowBox[{"\[Lambda]", "\[Equal]", 
   RowBox[{"ProductLog", "[", 
    RowBox[{
     RowBox[{"C", "[", "1", "]"}], ",", 
     RowBox[{"-", "1"}]}], "]"}]}]}]], "Output",
 CellChangeTimes->{{3.680365044856118*^9, 3.680365060918049*^9}, 
   3.6805886129984903`*^9}]
}, Open  ]],

Cell[TextData[{
 "Infinitely many solutions exist in this form: ",
 Cell[BoxData[
  FormBox[
   RowBox[{
    SubscriptBox["\[Lambda]", "i"], ",", " ", 
    RowBox[{"i", "\[Element]", 
     RowBox[{"Integers", "!"}]}]}], TraditionalForm]]],
 "\n",
 Cell[BoxData[
  FormBox[
   RowBox[{"\[Sum]", " ", 
    RowBox[{
     SubscriptBox["c", "\[ImaginaryI]"], " ", 
     SuperscriptBox["\[ExponentialE]", 
      RowBox[{
       SubscriptBox["\[Lambda]", "i"], "t"}]]}]}], TraditionalForm]]],
 " is a solution."
}], "Text",
 CellChangeTimes->{{3.680365153511257*^9, 3.6803652262458315`*^9}, {
  3.680365450167611*^9, 3.680365514299721*^9}, {3.680366636232228*^9, 
  3.6803666523204083`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"scharsol10", "=", 
  RowBox[{
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"ProductLog", "[", 
      RowBox[{"i", ",", 
       RowBox[{"-", "1"}]}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", 
       RowBox[{"-", "10"}], ",", "10"}], "}"}]}], "]"}], "//", 
   "N"}]}]], "Input",
 CellChangeTimes->{{3.6803652340023866`*^9, 3.68036529991215*^9}, {
  3.680365523765436*^9, 3.680365531237012*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{
    RowBox[{"-", "4.0637417027918294`"}], "-", 
    RowBox[{"58.0495734344775`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{
    RowBox[{"-", "3.9495227424225288`"}], "-", 
    RowBox[{"51.7601220040207`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{
    RowBox[{"-", "3.8205543078136768`"}], "-", 
    RowBox[{"45.46926540371086`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{
    RowBox[{"-", "3.672450068709818`"}], "-", 
    RowBox[{"39.176440021735246`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{
    RowBox[{"-", "3.4985152121541034`"}], "-", 
    RowBox[{"32.880721480068914`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{
    RowBox[{"-", "3.287768611544094`"}], "-", 
    RowBox[{"26.580471499359145`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{
    RowBox[{"-", "3.0202397081645014`"}], "-", 
    RowBox[{"20.272457641615222`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{
    RowBox[{"-", "2.6531919740386973`"}], "-", 
    RowBox[{"13.949208334533214`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{
    RowBox[{"-", "2.062277729598284`"}], "-", 
    RowBox[{"7.588631178472513`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{
    RowBox[{"-", "0.31813150520476413`"}], "-", 
    RowBox[{"1.3372357014306895`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{
    RowBox[{"-", "0.3181315052047641`"}], "+", 
    RowBox[{"1.3372357014306895`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{
    RowBox[{"-", "2.062277729598284`"}], "+", 
    RowBox[{"7.588631178472513`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{
    RowBox[{"-", "2.6531919740386973`"}], "+", 
    RowBox[{"13.949208334533214`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{
    RowBox[{"-", "3.0202397081645014`"}], "+", 
    RowBox[{"20.272457641615222`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{
    RowBox[{"-", "3.287768611544094`"}], "+", 
    RowBox[{"26.580471499359145`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{
    RowBox[{"-", "3.4985152121541034`"}], "+", 
    RowBox[{"32.880721480068914`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{
    RowBox[{"-", "3.672450068709818`"}], "+", 
    RowBox[{"39.176440021735246`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{
    RowBox[{"-", "3.8205543078136768`"}], "+", 
    RowBox[{"45.46926540371086`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{
    RowBox[{"-", "3.9495227424225288`"}], "+", 
    RowBox[{"51.7601220040207`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{
    RowBox[{"-", "4.0637417027918294`"}], "+", 
    RowBox[{"58.0495734344775`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{
    RowBox[{"-", "4.1662424475284165`"}], "+", 
    RowBox[{"64.33798412035904`", " ", "\[ImaginaryI]"}]}]}], "}"}]], "Output",
 CellChangeTimes->{{3.6803652721337214`*^9, 3.680365300908415*^9}, 
   3.680365564705014*^9, {3.6804327742925873`*^9, 3.680432802182106*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"scharsol10", "/.", 
  RowBox[{
   RowBox[{"z_", "?", "NumericQ"}], ":>", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"Re", "[", "z", "]"}], ",", 
     RowBox[{"Im", "[", "z", "]"}]}], "}"}]}]}]], "Input",
 CellChangeTimes->{
  3.6803657528961573`*^9, {3.6803659068488693`*^9, 3.68036590915143*^9}, {
   3.6803659445514946`*^9, 3.680365949509065*^9}, {3.6804327796592093`*^9, 
   3.680432782379562*^9}, {3.6804328247696495`*^9, 3.680432851136113*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "4.0637417027918294`"}], ",", 
     RowBox[{"-", "58.0495734344775`"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "3.9495227424225288`"}], ",", 
     RowBox[{"-", "51.7601220040207`"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "3.8205543078136768`"}], ",", 
     RowBox[{"-", "45.46926540371086`"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "3.672450068709818`"}], ",", 
     RowBox[{"-", "39.176440021735246`"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "3.4985152121541034`"}], ",", 
     RowBox[{"-", "32.880721480068914`"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "3.287768611544094`"}], ",", 
     RowBox[{"-", "26.580471499359145`"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "3.0202397081645014`"}], ",", 
     RowBox[{"-", "20.272457641615222`"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "2.6531919740386973`"}], ",", 
     RowBox[{"-", "13.949208334533214`"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "2.062277729598284`"}], ",", 
     RowBox[{"-", "7.588631178472513`"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "0.31813150520476413`"}], ",", 
     RowBox[{"-", "1.3372357014306895`"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "0.3181315052047641`"}], ",", "1.3372357014306895`"}], 
    "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "2.062277729598284`"}], ",", "7.588631178472513`"}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "2.6531919740386973`"}], ",", "13.949208334533214`"}], 
    "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "3.0202397081645014`"}], ",", "20.272457641615222`"}], 
    "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "3.287768611544094`"}], ",", "26.580471499359145`"}], "}"}],
    ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "3.4985152121541034`"}], ",", "32.880721480068914`"}], 
    "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "3.672450068709818`"}], ",", "39.176440021735246`"}], "}"}],
    ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "3.8205543078136768`"}], ",", "45.46926540371086`"}], "}"}],
    ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "3.9495227424225288`"}], ",", "51.7601220040207`"}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "4.0637417027918294`"}], ",", "58.0495734344775`"}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "4.1662424475284165`"}], ",", "64.33798412035904`"}], 
    "}"}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.680365704962778*^9, 3.680365753928321*^9, 3.680365910013301*^9, 
   3.68036595030929*^9, {3.6804327837495623`*^9, 3.680432852444702*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{"scharsol10", "/.", 
   RowBox[{
    RowBox[{"z_", "?", "NumericQ"}], ":>", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"Re", "[", "z", "]"}], ",", 
      RowBox[{"Im", "[", "z", "]"}]}], "}"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.6803652913503947`*^9, 3.6803653786896744`*^9}, {
   3.680365558184767*^9, 3.680365558523387*^9}, {3.6803657219421597`*^9, 
   3.680365741706524*^9}, 3.680365959980899*^9, 3.6804328740571175`*^9}],

Cell[BoxData[
 GraphicsBox[{{}, {{}, 
    {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.012833333333333334`],
      AbsoluteThickness[1.6], PointBox[CompressedData["
1:eJxTTMoPSmViYGAQBWIQnd/u3OjqKHAg/ZV0TgSb74E5UjqJ82fwH5gi17bW
+6HngZlmT4/VTeA70M3V+lhgj9uBHZlZDbqJvAeq2nWnTZ3mfMBBVK376x/u
A+8eO1TvLnA4UKtrvTDSl+vAZ8Zdx2dOsTrgo7s0v1iT40Csx5YT+11NDvDx
bNm2x5L1QNCl5kX/nmofyFGLL++uZzggdGKz4qEouQPWC1fbucRc2X+R2ds5
MPXrfitUvj2aegc08xzQ7HNAc48Dmnsd0PzjgOZfB7TwcEALL4d9sWdLrVcJ
HJj1QLVjoWiAAwBi8Zq3
      "]]}, {}}, {}, {}, {{}, {}}},
  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->{{-4.1662424475284165`, 0}, {-58.0495734344775, 
    64.33798412035904}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{3.680365379690632*^9, 3.6803655742978897`*^9, 
  3.6803657427228193`*^9, 3.6803659611491337`*^9, 3.6804328756699066`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["For general delay :", "Subsection",
 CellChangeTimes->{{3.680366010963182*^9, 3.6803660275704627`*^9}}],

Cell[CellGroupData[{

Cell[BoxData["dde"], "Input",
 CellChangeTimes->{{3.6807984680092287`*^9, 3.680798468355998*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   SuperscriptBox["p", "\[Prime]",
    MultilineFunction->None], "[", "t", "]"}], "\[Equal]", 
  RowBox[{"-", 
   RowBox[{"p", "[", 
    RowBox[{"t", "-", "\[Tau]"}], "]"}]}]}]], "Output",
 CellChangeTimes->{3.6807984688412333`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"p", "[", "t_", "]"}], "=", 
   SuperscriptBox["\[ExponentialE]", 
    RowBox[{"\[Lambda]", " ", "t"}]]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"cheqn", "=", 
  RowBox[{"Simplify", "[", 
   RowBox[{"Thread", "[", 
    RowBox[{
     RowBox[{"dde", "/", 
      RowBox[{"p", "[", "t", "]"}]}], ",", "Equal"}], "]"}], 
   "]"}]}]}], "Input",
 CellChangeTimes->{{3.680364051458361*^9, 3.6803641350394707`*^9}, {
  3.680364485551271*^9, 3.6803645186763525`*^9}, {3.6803645487479706`*^9, 
  3.6803646337351265`*^9}, {3.6803646839318485`*^9, 3.6803648784891567`*^9}, {
  3.6803649192123604`*^9, 3.6803649224980917`*^9}, {3.680366065274536*^9, 
  3.6803660984392357`*^9}, {3.680588646196726*^9, 3.68058865296996*^9}, {
  3.680798473529209*^9, 3.6807984873826885`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   SuperscriptBox["\[ExponentialE]", 
    RowBox[{
     RowBox[{"-", "\[Lambda]"}], " ", "\[Tau]"}]], "+", "\[Lambda]"}], 
  "\[Equal]", "0"}]], "Output",
 CellChangeTimes->{{3.6803640715488267`*^9, 3.680364135808674*^9}, 
   3.6803644452568417`*^9, 3.680364487069353*^9, 3.680364519823472*^9, 
   3.6803645514043255`*^9, {3.680364582532942*^9, 3.6803646344139333`*^9}, {
   3.680364688455904*^9, 3.68036475263435*^9}, {3.680364787150921*^9, 
   3.680364806780751*^9}, {3.6803648515092373`*^9, 3.680364879341175*^9}, 
   3.680364924382612*^9, 3.68036610110987*^9, 3.6805886550030093`*^9, {
   3.6807984753454995`*^9, 3.6807984882203317`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Reduce", "[", 
  RowBox[{"cheqn", ",", "\[Lambda]"}], "]"}]], "Input",
 CellChangeTimes->{{3.6803650315542855`*^9, 3.6803650593644743`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(", 
   RowBox[{
    RowBox[{"s", "\[Equal]", "0"}], "&&", 
    RowBox[{"\[Lambda]", "\[Equal]", 
     RowBox[{"-", "1"}]}]}], ")"}], "||", 
  RowBox[{"(", 
   RowBox[{
    RowBox[{
     RowBox[{"C", "[", "1", "]"}], "\[Element]", "Integers"}], "&&", 
    RowBox[{"s", "\[NotEqual]", "0"}], "&&", 
    RowBox[{"\[Lambda]", "\[Equal]", 
     FractionBox[
      RowBox[{"ProductLog", "[", 
       RowBox[{
        RowBox[{"C", "[", "1", "]"}], ",", 
        RowBox[{"-", "s"}]}], "]"}], "s"]}]}], ")"}]}]], "Output",
 CellChangeTimes->{3.6803661254722557`*^9, 3.680588659339037*^9}]
}, Open  ]],

Cell["Different structures for s = 0 and s > 0", "Text",
 CellChangeTimes->{{3.6803661610866146`*^9, 3.680366198534652*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Manipulate", "[", 
  RowBox[{
   RowBox[{"Column", "[", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"ListPlot", "[", 
       RowBox[{
        RowBox[{"tmp", "=", 
         RowBox[{"Map", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"Through", "[", 
             RowBox[{
              RowBox[{"{", 
               RowBox[{"Re", ",", "Im"}], "}"}], "[", "#", "]"}], "]"}], 
            "&"}], ",", 
           RowBox[{
            RowBox[{"Table", "[", 
             RowBox[{
              FractionBox[
               RowBox[{"ProductLog", "[", 
                RowBox[{"i", ",", 
                 RowBox[{"-", "\[Tau]"}]}], "]"}], "\[Tau]"], ",", 
              RowBox[{"{", 
               RowBox[{"i", ",", 
                RowBox[{"-", "n"}], ",", "n", ",", "1"}], "}"}]}], "]"}], "//",
             "N"}]}], "]"}]}], ",", 
        RowBox[{"ImageSize", "\[Rule]", 
         RowBox[{"{", 
          RowBox[{"300", ",", "200"}], "}"}]}], ",", 
        RowBox[{"PlotRange", "\[Rule]", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{
             RowBox[{"-", "50"}], ",", "10"}], "}"}], ",", 
           RowBox[{"{", 
            RowBox[{
             RowBox[{"-", "100"}], ",", "100"}], "}"}]}], "}"}]}], ",", 
        RowBox[{"AxesOrigin", "\[Rule]", 
         RowBox[{"{", 
          RowBox[{"0", ",", "0"}], "}"}]}]}], "]"}], ",", 
      "\[IndentingNewLine]", 
      RowBox[{"Row", "[", 
       RowBox[{"{", 
        RowBox[{
        "\"\<Max(Re(\!\(\*SubscriptBox[\(\[Lambda]\), \(i\)]\)))=\>\"", ",", 
         " ", 
         RowBox[{"Max", "[", 
          RowBox[{"tmp", "[", 
           RowBox[{"[", 
            RowBox[{"All", ",", "1"}], "]"}], "]"}], "]"}]}], "}"}], "]"}]}], 
     "}"}], "\[IndentingNewLine]", "]"}], ",", "\[IndentingNewLine]", 
   RowBox[{"{", 
    RowBox[{"\[Tau]", ",", "0.01", ",", "6"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"n", ",", "0", ",", "10", ",", "1"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.6803662182537413`*^9, 3.680366235985077*^9}, {
  3.6803662680882816`*^9, 3.68036632256427*^9}, {3.68036637772369*^9, 
  3.6803663907540164`*^9}, {3.680366426900672*^9, 3.680366440109706*^9}, {
  3.6803664864053564`*^9, 3.680366518421238*^9}, {3.6803669214175043`*^9, 
  3.680366987382147*^9}, {3.6803670220955076`*^9, 3.6803670513623962`*^9}, {
  3.680367330191185*^9, 3.6803673481999316`*^9}, {3.6803674013106565`*^9, 
  3.68036745287344*^9}, {3.6803680112880955`*^9, 3.6803680116503563`*^9}, {
  3.680368051239793*^9, 3.680368055898922*^9}, {3.6807985372572627`*^9, 
  3.680798547376212*^9}, {3.680798588074071*^9, 3.68079862602701*^9}, {
  3.6807986884909678`*^9, 3.6807986897394443`*^9}}],

Cell[BoxData[
 TagBox[
  StyleBox[
   DynamicModuleBox[{$CellContext`n$$ = 10, $CellContext`\[Tau]$$ = 0.01, 
    Typeset`show$$ = True, Typeset`bookmarkList$$ = {}, 
    Typeset`bookmarkMode$$ = "Menu", Typeset`animator$$, Typeset`animvar$$ = 
    1, Typeset`name$$ = "\"untitled\"", Typeset`specs$$ = {{
      Hold[$CellContext`\[Tau]$$], 0.01, 6}, {
      Hold[$CellContext`n$$], 0, 10, 1}}, Typeset`size$$ = {
    660., {232., 245.}}, Typeset`update$$ = 0, Typeset`initDone$$, 
    Typeset`skipInitDone$$ = True, $CellContext`\[Tau]$92672$$ = 
    0, $CellContext`n$92673$$ = 0}, 
    DynamicBox[Manipulate`ManipulateBoxes[
     1, StandardForm, 
      "Variables" :> {$CellContext`n$$ = 0, $CellContext`\[Tau]$$ = 0.01}, 
      "ControllerVariables" :> {
        Hold[$CellContext`\[Tau]$$, $CellContext`\[Tau]$92672$$, 0], 
        Hold[$CellContext`n$$, $CellContext`n$92673$$, 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" :> Column[{
         ListPlot[$CellContext`tmp = Map[Through[
             {Re, Im}[#]]& , 
            N[
             Table[
             ProductLog[$CellContext`i, -$CellContext`\[Tau]$$]/$CellContext`\
\[Tau]$$, {$CellContext`i, -$CellContext`n$$, $CellContext`n$$, 1}]]], 
          ImageSize -> {300, 200}, PlotRange -> {{-50, 10}, {-100, 100}}, 
          AxesOrigin -> {0, 0}], 
         Row[{"Max(Re(\!\(\*SubscriptBox[\(\[Lambda]\), \(i\)]\)))=", 
           Max[
            Part[$CellContext`tmp, All, 1]]}]}], 
      "Specifications" :> {{$CellContext`\[Tau]$$, 0.01, 
         6}, {$CellContext`n$$, 0, 10, 1}}, "Options" :> {}, 
      "DefaultOptions" :> {}],
     ImageSizeCache->{758., {369., 382.}},
     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.680366309155188*^9, 3.6803663299534636`*^9}, {
   3.680366429232878*^9, 3.680366441996362*^9}, {3.6803664897848406`*^9, 
   3.6803665373532114`*^9}, 3.6803670420295486`*^9, 3.6803673495008984`*^9, {
   3.680367414467495*^9, 3.6803674535224943`*^9}, 3.680368016006957*^9, 
   3.6803680584808607`*^9, 3.6807985037707367`*^9, 3.6807985495773106`*^9, {
   3.680798598886139*^9, 3.680798626759019*^9}, 3.6807986906101885`*^9}]
}, Open  ]]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["NOTE", "Section",
 CellChangeTimes->{{3.6823996793515368`*^9, 3.6823996898574977`*^9}}],

Cell[TextData[{
 "Oscillation appears in the solutions if ",
 Cell[BoxData[
  FormBox[
   RowBox[{"\[Tau]", ">", 
    FractionBox["1", "\[ExponentialE]"]}], TraditionalForm]],
  FormatType->"TraditionalForm"],
 ", i.e., all characteristic roots are complex!"
}], "Text",
 CellChangeTimes->{{3.6823996913924026`*^9, 3.682399814870842*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["TASK : ", "Section",
 CellChangeTimes->{{3.6803676674222765`*^9, 3.6803677053610244`*^9}}],

Cell["\<\
Find real solutions using the characteristic roots, and investigate the \
asymptotic behavior of the solutions.\
\>", "Text",
 CellChangeTimes->{{3.6803676674222765`*^9, 3.6803677509152946`*^9}, {
  3.6803690691964426`*^9, 3.680369120885888*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Localize the characteristic roots by visualization", "Section",
 CellChangeTimes->{{3.6804465700762625`*^9, 3.680446603798191*^9}}],

Cell[TextData[{
 StyleBox["Idea:",
  FontWeight->"Bold"],
 " Visualize the Absolute Value of the characteristic quasipolinomial: "
}], "Text",
 CellChangeTimes->{{3.680446688964653*^9, 3.680446823990629*^9}, {
  3.680446864483438*^9, 3.680446869658433*^9}}],

Cell[TextData[{
 Cell[BoxData[
  RowBox[{
   SuperscriptBox["\[ExponentialE]", 
    RowBox[{
     RowBox[{"-", "\[Tau]"}], " ", "\[Lambda]"}]], "+", "\[Lambda]"}]],
  CellChangeTimes->{{3.6803640715488267`*^9, 3.680364135808674*^9}, 
    3.6803644452568417`*^9, 3.680364487069353*^9, 3.680364519823472*^9, 
    3.6803645514043255`*^9, {3.680364582532942*^9, 3.6803646344139333`*^9}, {
    3.680364688455904*^9, 3.68036475263435*^9}, {3.680364787150921*^9, 
    3.680364806780751*^9}, {3.6803648515092373`*^9, 3.680364879341175*^9}, 
    3.680364924382612*^9, 3.6804342885816507`*^9}],
 ", and from here"
}], "Text",
 CellChangeTimes->{{3.680446688964653*^9, 3.680446823990629*^9}, {
  3.680446864483438*^9, 3.680446961411192*^9}, {3.680447889241246*^9, 
  3.6804478958909836`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Chf", "[", 
   RowBox[{"\[Lambda]_", ",", "\[Tau]_"}], "]"}], ":=", 
  RowBox[{"\[Lambda]", "+", " ", 
   SuperscriptBox["\[ExponentialE]", 
    RowBox[{
     RowBox[{"-", "\[Lambda]"}], " ", "\[Tau]"}]]}]}]], "Input",
 CellChangeTimes->{{3.680446688964653*^9, 3.680446823990629*^9}, {
  3.680446864483438*^9, 3.680446993105179*^9}, {3.680447036417987*^9, 
  3.6804470375173883`*^9}, {3.6804479182163677`*^9, 3.6804479305507092`*^9}}],

Cell["\<\
This method will be quite important in the case of several delays when the \
characteristic equation cannot be solved symbolically.\
\>", "Text",
 CellChangeTimes->{{3.6804476864748335`*^9, 3.680447754890559*^9}}],

Cell[CellGroupData[{

Cell["A special case", "Subsection",
 CellChangeTimes->{{3.6804466582796087`*^9, 3.6804466721373463`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot3D", "[", 
  RowBox[{
   RowBox[{"Abs", "[", 
    RowBox[{"Chf", "[", 
     RowBox[{
      RowBox[{"\[Mu]", "+", " ", 
       RowBox[{"I", " ", "\[Nu]"}]}], ",", "1"}], "]"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"\[Mu]", ",", 
     RowBox[{"-", "5"}], ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"\[Nu]", ",", 
     RowBox[{"-", "40"}], ",", "40"}], "}"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"0", ",", "20"}], "}"}]}], ",", 
   RowBox[{"PlotPoints", "\[Rule]", "100"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.680447072809614*^9, 3.680447217617277*^9}, {
   3.680447274473708*^9, 3.680447323053584*^9}, {3.6804474130227613`*^9, 
   3.680447494664062*^9}, {3.6804475557969604`*^9, 3.680447639468896*^9}, {
   3.680447959740569*^9, 3.6804480233889637`*^9}, 3.680799183792692*^9}],

Cell[BoxData[
 Graphics3DBox[GraphicsComplex3DBox[CompressedData["
1:eJxk3XWUlWXf9nGlBUS6Y4IZJmFmNindoLRgEALSJY2kgIQSClIiIUiXIogi
[encoded data, 744840 characters]
uXQcxsn8/gu5XHxA7nt/vQ76xb0r8LzcvJBbR9Ny/h96SBiQ
   "], {{
     {RGBColor[0.880722, 0.611041, 0.142051], EdgeForm[None], Specularity[
       GrayLevel[1], 3], 
      StyleBox[GraphicsGroup3DBox[{Polygon3DBox[CompressedData["
1:eJxMnXd4FsXXht/QO4QQeu8CERGRTkIKKEiRGorS8oZeVJqgtAAh2JAiCihC
[encoded data, 383700 characters]
hv9QrA==
          "]], Polygon3DBox[CompressedData["
1:eJws3WXcVkUTB+DzKN0g3d3doICiIHaL+aoIioqIYHeDga2IYgJ2B9hioSIW
[encoded data, 97920 characters]
+YXPIgw7cXHtkT3ax0NcDfpbecbkXnuHryL/B7B6YS4=
          "]], Polygon3DBox[CompressedData["
1:eJwl13mcz3UeB/D3t3ENw7iPDHKEKDMMGVckV+XKmY5N0jS0Rrtqa+uxbZtx
5XhgS22Oh9lhkSyFVNTy2FIr0bWyi0dEB1Y5Wh5q2/b5XX+8H7zfv8/7/Xq9
r8/nO03HTR5WfFlEbCDlyZP5SWTXS2IHmbmI3jGJsoIkcr+KaHUiYvDUiF9f
lcTcBkkUdI548LOIhf+KmEO/qmYSn98TUfvriE1HnGucxFS2joV82U6z/Zbt
E/G2nImo9n7Elssj3qa3qpPExEERj38ZceKomM6tIk/D/JJMwiW3RhJd8KuC
W2bFJBo2iiigz+G7LCOiBJf/HI+YS1/D901+fyGf3hEx0vmX2Gv+N6IG6XV3
xMv0zngPWhlx4FTEP8kGfkvII/h+xbfNs/jBu6l2Eg0mRNzBPuBzWCk3Uu9Y
RBnO2UMj2ou3nK0hDsNORlQtFaN+EvtbJDE8M4lyFZK4pVbE7qZJXM12Q5Ya
fRRRXS0+VINZfBf+FLGAZIyKeBFm38uT+FtORH0Y9+PXj74+N4nNtZI4rR6P
j8OZfob/CBjVYbwO40MYGdck0WhMxPbySSyuJoc2SfTuoi5w7m2SxPTDEfOa
JdGT7/hDen06Yirp3xk3tiK2mfRZ5Ea2v7KV4lrC91UxrhfrTnmsqpREe/E7
9YxoBuNcB7/Pjfj9jeKbm4/5rW6f/H92WqrpWbam7ZLogdf5KhF3fhGx3Gws
FWcc/g/UjRiB2xdsf5RTWv+0D1/zPaq2ZbAr6eei282JWI31I7eyOa2axPtd
I3bq79vkFefmmd+nzU6+WtURZ6N52auGT5n18Xo35ZuID+jbyyXx/ciIfWr/
d1KAyzR1zS+O2Krn6fyne1Ag35V+ayvf6vL90S5kyvforRG/k09RXdzoM+HO
g9sBbi24i+E+i8+1cu0g1qjR5sP5kzg374Gr2b0oj6Iq9kXflunbHLF2wHlB
Tdtk4Z4bkWnXZpiDkTjniFlR7PYwlqpFI3q++h2prV/qd0j9fgNzg97l7Iu4
Te9mss/X82+rJ9Eav2Nq2pmscG6/Hi0xT/8W7xReOfqbDbuZPL+wJ8dIR7Z8
NVgrdo4Zv2ax+ZNrvSuTqKZvVckgGPvUagGcWf7/U3oW/kt4THBuUv2IQvYi
dVjp3FZ13mKez8Atgb+NXtgQTkP9dS7n24j76OnMpLMzXbwO+lLiXCmM7AER
3dSjK+mM3/ViFprdaeZ2OhkAt4o9fKtCxBhzcSvJor8s34t29hu4H3Q3+3q5
gu9dMGrBXVMtop38dzuX6MU556rp1VLnfmCf3CdiGfunevQJfS4u0/i9oc4p
fi+xNqvrIfwPkm6t3DH6mO5cuntvOpfuWG/n5rCVsG1jS3fxuoJ075O4H/+N
bE1Sm1gf4bJIjb7D5YQeNZTvbPYTcjpO5pv/mnK7yH+v3XnOPPS/LaJHnjly
nxSqwavmewtZS2/D/5x/j5u7lu0iJpq9C3yf53sb38p8b2K7mJveK/YVbhn8
7+lL1Xug/XkNx4N7Iio1SuIes7HafnRQk27mprW6LNDzD8WrnW2m+rG1xZnt
A7ZabHXZmrEN06cf9tq/9M5Ri3R+XpJbwFjvTViTF/GZ80PSfRsfMdTbUgan
iRzewSf05ix+r+H5Ln022aGvp9g2sc2hT3Z+iR2sJ/5YMbqLUV5truN7wblt
zlWg34P3m73Nst8f0uN8mM/gk2HWBuFQQR+rw50kjz3yKHGuq3fueTO6yz60
WBZR2R1V1X4W2ctm/CeK85O9O+GuaE5/z0yPNhd9zHNxy4hd9O2wa3vXpsGo
aDZ30EvTe+2qiMf091HyopoM6XLpDr5RrV5I90Hdd+PUwIw9lt4j7om1fi8V
Z9NNEUO8H93Yn2TLFbOf2X0Y5+/km0c/r+d/kMcYPS+n5330fLt828l3hBhV
xa4Js8T/N7rXN6vfLD1pzN6J/Wt3wwwztLGt2e8kJvsCvw+Auw7+cn4LYdyt
5jnX2yE92Es/rlaj1aqKWtVXq2NqtQnH1urSiu9qWLvwfB2X1rgMFydL7Bow
/6GGg3FtoG9r1XY/vakZ3K0mxe7Tq32fNKO/L2Y9sZ+yFzfDOgpzIq4V8HvD
ufZ6/40ab4X7Lt93yJR8e6JXLeBtw7VYvrv1dzEuU/kO5XuBb/oGpG/BEjvU
H95yb8oGdZ4pVjon6bwc4J/B9gTbHjmf0a+P+0Z0Eu8Z85Inl3n4zyVPert/
TN/u/EvvW/rOZeB+GT3sxyZ3wivk5HsRl9GLccnBpQ+MUnmcl8cFdVjsTl+N
c3mzdJ7e0H5l2bPqOJ3kkw1zuT5dtPOV1LOXXe7CVgt2oR71kOcD8mrE9gne
gw9futt/UffSm5xpPpp7o3Z4w5Y4N8q+bFbXO9W4uVm4boX7VPzd/n/WLJ4m
O9tEjIM9Vq3Kecfv8672bm5nvN0V9Wm2b6h1ej5KrabDbAy7VOwhuBzEqbc8
7nauKNt9iE++e3EsfTLsZc49on55ONXF7Vq+XeRUUW5fynEuLgfSbz1v5xF9
m1TZXZxy04+R8qlh7rer/afyz+K7Es8pYvY1q23Z/my/7qrhOwBuOzKRrYpz
Zc790rmezl3JluM+eE68QeKtrxPR0dyeglGHbZdajROj2LmWdnuJGmSahYVX
8KU/Sk6zFeHXZqw9ptcX77AzA2GuFS9PvHlpf+U6RV3q2uM98h+vpndNlJea
DlW/aelbgfMwuBm43Cy3nXALyDBcH8J5Ne7V0vsfv1HO1XRup3OHnLnFvGWZ
uxn62djOvKW3z8Pddy9fuAf4X8F3S9on+zIBx9xMs+s7qpRPnm+D0bi0lluh
/t56re8S3/q99eQGMoa+zr3zsLfkW/v7c1wqi1dAuvu9B6ljd1bJ95DendSr
BmSSb+aD9HVwD5upIyS3zPe33+rblaPmfb4eLyDrnalpR7bYl9fIBTtTm77B
Xo+033l29Q5vy8vObVWD0eZpFEm/jxfhMcLsvGu3GjlzbKB7w+4Mh73cHfe0
eM+Qbs6tkO9z/sZY47shX95l9O/l9DP72Mtevs5/o5jN5XvWffsr9Tukji3o
ffToT97LRWLluEsy9TL9myH922G4WB/b4xf4dtCTgfr2oLqvUquqcG/Xj8fN
wS3infJWLeT7P5XQxxg=
          "]]}],
       Lighting->{{"Ambient", 
          RGBColor[0.30100577, 0.22414668499999998`, 0.090484535]}, {
         "Directional", 
          RGBColor[0.2642166, 0.18331229999999998`, 0.04261530000000001], 
          ImageScaled[{0, 2, 2}]}, {"Directional", 
          RGBColor[0.2642166, 0.18331229999999998`, 0.04261530000000001], 
          ImageScaled[{2, 2, 2}]}, {"Directional", 
          RGBColor[0.2642166, 0.18331229999999998`, 0.04261530000000001], 
          ImageScaled[{2, 0, 2}]}}]}, {}, 
     {GrayLevel[1], EdgeForm[None], 
      StyleBox[GraphicsGroup3DBox[{Polygon3DBox[CompressedData["
1:eJw0nQd4FUXbhs+hS+8d6SCdAKF3kE4ghSIkhKaAHEFIKJKcgICKCpKQEwki
[encoded data, 586920 characters]
dSaEp5rOQ4QW4AracxFeAVciviC7pj7rh/8HsaHxzw==
          "]], Polygon3DBox[CompressedData["
1:eJws3XccjtUXAPDnFULZe69foaSUMiuSVUb2Hg3KatESQksZEYqUlBktFKXS
[encoded data, 133320 characters]
pGWlkZ8LjdjfItMVeeHISwb/H0I3hFU=
          "]], Polygon3DBox[CompressedData["
1:eJwl13m4j9UWB/D1mmfnGJPQEcpVpm7mEyrzcCQUEk6TipNuN5W6SZmlSaZE
KBzXUEoXjXqKjI0UDYRGGlTmwv28T398n/Pb33dNe621194nK/f2nnkFImIv
FIb+FybRv3fEzEJJ9CydxP0XJHGgWRKZxZPYUDiJTpkRfbKSyK+VRHncm7g/
MiJuwJ0kV6hEEjuKJPHn2RG31kxiZsMkcltGrGKrc8kkijdP4pXGSeRNirjw
aEQ96HlREpOqJ/HJvohxByM2LYh4sUoS37L3ZbEkxvPxW5mIw+cmcRBXjt/3
cDliuZbf73G7yU3EnSaX4F5l7/AfEX9AhwYRI8smsVrM3eieJndh+YjF5L6j
+xXdCbhjdE/ysZbuN+LaD4/Z20u+bZeX0X0i1lZIoqN9/FQniX10d9Ed4/se
uvvoTrbPI80j5v4QMSdFGge/T/LbXl7GVIwoIS/P0F3F/suwXCzvimW83NxJ
fgduO0wl8wOUpfsuH9n225ncYLEMVqNZanQVfxPUaLzfb+VH1DkVURvKiusb
ul/4O47uQfH9Ir575CGnE5/ysgxG8/eG+LrzcYZcY7Es52M/3c/pjsV9S/cH
ujPJbjkWsRluLx0x0LcfyWXQXe/35eLrQXek+EaK72jlJPqKb634JjZJouhr
EY3sqyFcqeaH6U6j21teTshLVXmZIpYyuLVpzuW+HXtHyE3H9SF3hlwWueO4
rbincFefFdEb16l2EptxY3E/VJYrXEX25ttHU/Y2lIx4zz52i2eAuOZX5E+c
U8Q7H9cLNxuXh8vDHTtPz9Jdiq+jh2vVtQc1OqLHK+GuaRJxPntFcH/gKuJ6
4rJwH4llPG7tWUnMeVYu/plEXXJr9M9qWFk0ogl7K8g1JleoaBKHzonYJuaL
yb1B70040D1iTz39frYeKpjEE1CiTsQG64H8fFooYoaYZ0J9etlwTbUkblSr
yrsjispflhw8Zx/d8R/LwSf0cvkJeZuGmw6X0msEg9RsINSsGfFm/SSGycHY
GhHl2C8PF5J5z99qB+juiXh8cMRG693ns18qiQp+l+0aMUzMR+n3Y2ey/fWz
5/Pp1knj87svFM+JOE3ulFpeTKa0HPStFvGx2P7CFcF9ppZFq0aMwG33fb29
l9fr5WCb7/nk6uuniyCdKavMlq/JNtfbzWD1lIhDuC/sY6Scf0Iv7xx50I+7
9GDugIhc+y2fkcRcsexhr6Nc7SC35CNzranekMsDztSP0LVgRHV1nkPuGT5n
wwjfr7WvaXLQUQ42OGO7u0R0ZG8+rjPuPdw3uO64MDvyxb5V3gvVMl/EFvrg
Kb3b0X4f1+OV7OFdfHWxTBPLGnVuq27ZfJX/ytzjtz+st5dMujPodqH7HN26
dMeTW/+nuQJVzPUWbFzHX0VntwKkc3KA9cFK6n9YLWFsR2danNnsbWTvYfa+
dI4uZW8vPx/I5/vQ1jz92vpycpvIjSG3j9wV5DqyWfrjiFtORgyB18Xxe1pj
fXybPFeR5zesZ+jfGgXsDw7okenWU8meTMwC2I+boN5N5O4IfzN+JwP9yI1q
4A5pFfFTZhIvl0/iQev+4h4o/tP2cQpurJpEAzVZrB8XQRmX2wz5u469go0i
jofzDXV9G+QcbtYbXeV6IXu101qKpb2z1g4uWudOEkst/HO+Pw93mAe96dXE
zbOeD3nuh7fNjYb4pex1YW8B/jwyF8D/9Pwr8KRvlZyNWWRnmTULxTBYTx01
g2rpjWV6Y4g89b0yohqZ6riluJtw3XCVccHef9P+F9cmPfytXj6tJyfg7sXl
4DbgWpH7TK062U9/fTOTTIJbXiKN0ZwkV0gsfczHLtdEfKc3voVVYr9TXleK
/6A8D8qN+Lf10+wcV4dj0OXTiHz1niy+39k7YV7tmR/R+BJnF7cE15T+vFsj
zvCRgR+lD6rh1pHNtM7VQ5N36DG9nA3TxNdfje9gfzjs6qzXrdvJW+FbIn4+
o+7wsDiqlNEzl9rncecAGujlG9h7gr1mbDWFGeydNP+m6tMO+nSS85GhT6fj
W/jeHFL5Fc759XQf9bslLv2WyqTxpHG1gvTbMnJncBNLpPPE3SV/W+X5PP33
n/bO25GId6Cf/hvgnZD/D3PeGbwENnnT7JDnt92TW9jbDCPt425yHeq5h05E
PAI7yX2sF5Y140/u56Z3UqMkHiO7Sm+/DDX7O+e+zZTLRd4ZFb9X27r2a/0o
/TLmRWJeBTxm/TO/D5eLeI3PV+E+tuqIcZlaL4Vb23iHyeehGkkUGyQWtkvq
ybP0SlX5W6AfmrGx1f2xzTnqbzblmPMf2HdXZ6o8uan6bJg6DIV31OVJ67Jk
27gTy7F1SM4Kkzssz4v8flr+jp8XsYit13EDcc1xHd03DXFFxdNNXobIyc2w
QszF2WvBXiZ7v6b3JnuP4O+T35HQt7h9i3NJOVyHiLvVYgQsIrMQJvk9ERb3
c7/ph2z3cjH9ea68VXW+stJ7EDfErKqOm6JfV5tXW8Q3jr888S3K8vYSXz2y
mWqTAaPEfDNuLt0TdJemdVO/9epYgVwPOdouV8Pk7KQ38XC/b0rnFGy8mZx4
B5CpafZlwUvdyJF5xVnPNVc+F9sV+mSbc3U7Hx/qu+58jOErz/qk+BaLb7b4
TsvpMrE8iX9NLZeQu1q+qohjFz/NnKly5ngm5MvJ+3p6pXPcgm53+/jA+k89
Xpi9T+Uocf8Od2bK6oN5+qCJ/K7TB2+nbx/cnbiCuDy5X43rIcaD6X3OZxm1
6iUH+c768IXOOB/POvuHxTaF76rqejZMvNr5cMaHOut3ketjNveGJ8gEH//i
owAfQ/hYyccO+2uP2+xeyzBPF+IG0a1jZtSGYg2df/ksae78Klc90npY1xPL
T/ZWwTzYxN61zuIt3kjt2En/B5r6l76F3eLc790Q8lvE77+8zfZZf0n3lPVw
8pv847TfefkFl8NeSbma7f38Env99Mwh931ZvdGN37XeAb/Sf9y5LUo/YfcX
67bOSNWNEe35bAe5zurZ6jRQ7laIeZf6PeDs/8pHDz5KFUnvmog1fBzCXYkr
jVuMexP3Fa6GvMwRX5Y7r7281IeB6lCZzWPqUoz9q/TTGn0ywf07Hhak94Me
2sbehPT/N/f5ADW/RP895M4ZRn6BfVSmWwBeYGe5fex9K+IVb6lVdCZ7aCx0
tp6HWembiW622dIKdq7xvwa5r8XXqVj6blVffZ1P5i6yK82kF+HuCuohhjZ8
3IEfDv305GXq/bm+vMrvlvyu0q9fWBck97NYMuztYnssYK+t+Ci1zf8Y4ugB
La2Hqnsd34fy+xv5B/ToTvwnfF6nBwZk/v3/1jlsl67u/+HfvHEgm+xo995q
74263jMXwEPWz4rrhH0dh16tPaXL/f0mfADfVS6WuGcayN/TaT3Mw+owwrto
kpj34p6Xgyv032DcaFxNs+5DMzOf3GJ4lJ0ryV3n781iHqrXelgPZvu2St7U
4p4Eb/j+oB5aL653YZr7ZlT6JiTb8Ev3hO+PWLeQp4ugin48C+bK4zg5/Ulv
5Y+PaGMfraG0HursrB5zj+fxezz9H0au1uuNsX4Pk5955t9Uuundmt6xTaD+
HHeLM52hnq31y/X6ZbZ+qcjnvc7C+eZeZX574o46Cw3I3U9uKLl5uErkdrPX
QV628zvKmd6pX0aL/X65P2hvB+Ae+zllTtTzP2xbs7MNPI7rpu7NZ3lnOFen
oXE6c/n4ppd+Yq+U/sjRAzn2sUV/jZPHX/R4DzWaSz+H7e6wjq8fnbv/4C/3
Xmrn2xXQmr2D7F0mR20hzdnPcjfCTH4hvZvgr+yIgurY0D5C/73PRwe6DaxP
yM2r5l5Vd8oB75iT1oX1zGnvqqPqfQRK071Mrlek94y8TJSfVnz+HxBNVMM=

          "]]}],
       Lighting->{{"Ambient", 
          GrayLevel[0.8]}}]}, {}, {}}, {
     {GrayLevel[0], Line3DBox[CompressedData["
1:eJwl22W4FWXbgOHZdHd3d4h0l3R3l4KKgAgWKIgoqSIgUtIoIC3d3d3d3d39
[encoded data, 7440 characters]
/wFyVzQi
       "]], Line3DBox[CompressedData["
1:eJwd2nf8T+X/x/H3ByUSksie2SIhZc8kGtp777QHSkv2jKZU0k57p5Iy2n0z
[encoded data, 7680 characters]
q//eWCdth8/g20obp/hYvqyj3/O69vN746L//z/r/wGZ7lIf
       "]], Line3DBox[CompressedData["
1:eJwl1MdSVFEUheHDi2g5M/sa5gECkkOTocnQ5Cg5mrNlleUrOPAZFHVkRCkQ
aJIEnTrwO+Vg8++11m3uvWfv7mOJpszGjBDCUX9e/Q1hbSOE50TXZgjp7RDO
r4TwNh1CzX4IWWshbMi7ZT1qkV/Hv87f5Pfy+tQ7/sZWCPn8bX4/b0C95x/n
f4ifWw2hYDeEeizEzOUQijFJF2EDXuM1YpPaPpIRsuhmfYvK1rdim8rRt2OJ
z3XgmwPPxCujU3QpdmICc/kV2EWXYzdWYh6/GnvoKuzFGtxx33xZnb6PV4v9
WI8F/BPeZ4Au9q6DeJIu5J/CofhM/GFM06dVkewMbqqR+GzyG9i4F8JZXrH8
HG6pUX6ZfCyegbxcv2UmCWymW9R4fP54zvxWegcr6Db9C3NsxwnXVPN2ZR30
L6yiU/pONSmvo/f5XfQB1tLd+h41JW+gf/Nfm/cFO5Gk+2QrzrpJv4oDdIXn
/+O6bfpZfC/ZU1yQ9cd3kB/Kt+RP4jnKH+O8rFcl5HvyTfmjOAP5Q5yTLZnF
MN3h7Idwlpd23QzW7YTwwHXT+nb5oHxKvyGfxFr5ffmEvk0+IB/Xr8vHsEZ+
Tz6qb40zjnPRr8lHsFp+Vz4czz3ugnxI/1M+iFXyO/LmuDuyZef0Q1Xyb8f5
xd3ivzST7/wlVSG7Fc8o7qHsG++rKuffjGce95T/hfdZJfgL/GTcZ/4n3kdV
xp/np3hXzSYln6M76Cv0rL49/n/nd5meid8b+pJ+Wt+qv6gvinsVz/vQO2Kh
2vH9LcBkvI/ZrNvLvLiH/Fys55fy1/g5dJqfjbX8Ev5q3HXz/ImL7r8ur5Y1
ukdSNRz+/935ByMvmIE=
       "]], Line3DBox[CompressedData["
1:eJwl1GWUFkQUBuAPyZUu6e7uRkBJAeku6e4QAWkMUBpRQIUFlmbZoFO6u7tD
AQFpWJRnjz+eeX/MnJl7z8yZLO37NOgdIxAIZDLMigoEfuFnZvITM5jONEIs
miqnMJlJTORHfmAC4/me7/iWbxjHWMYwmlGMZATD+ZphDGUIXzGYLxnEQCo4
e4DsTz/60ofe9KInPehON7rShc50oiMdaE872vIFbWhNK1rSguY0oylNaEwj
GtKA+tSjLnX4nNrUoiafUYPqVKMqVVigh8ryUz6hEhWje+NjylOOspShNKUo
SQmKU4yiFKEwhShIAfKTj7zkITe5yEkOspONrGQhM5nISAbSk460pCE1qfiI
lKQgOclIGl2bfpLIxCQiIQmIz4cEEY+4xCE2sYjJB8QgwH9vA4F/eUcUb3nD
a17xkhc85xlP+YcnPOYRf/OQB9znL/7kHne5w21ucZNgtd+Q17nGVa5wmUtc
5ALnOcdZznCaU5zkBMc5xlGOcJhDHOQA+9nHXvYwz7m75S52soPt/ME2trKF
zWxiIxtYzzrWsobVRBJBOGGsIpSVrGA5y1jKEhaziBAWsoD5BKkpWHZ+ZH/3
sZpIIghnkfkwuYpQVrKC5SxjKUtYHL2WEBZGv3fmE8w85vI7v/Erc5gd9f+/
8x6NgadP
       "]], Line3DBox[CompressedData["
1:eJwl1Edz1VUcxvGTlVs2stThHaCi0mxIFNnQlu6cAYSQhDRSgVTSKAHSG6mk
AoKyAXkHoCNtxyA1kHpTCC0wfs64+Oa5v+9z/3dOzjn3rvglfXtaUgjhZ38+
eB9C3+MQlsZDeIctD0I4fD+EazMhvDf3685PhzCXCOGcrNW94f/Uv5W9+lF+
Vj8ia/Sv+csPQ+jR/SGH+Wr+Ff+7uZu/JDc/D6GKf8lfNHfxv8mBj5PCT7pK
3SZ5gTsT1yF/NB/hf5AXuQXrKJywNl3/pOe5eW4OBfwo38df4Ge5BPL5Eb6X
P8/PcNPI44fjuvlz/BQ3iVx+iB+0rm7dqG6CH8cB3WBcOz/CP+eeIYcfiOvm
h/kx7imy+bN8Jz/EP+EeI4vv5zv4Qf4R9xCZfB/fzg/wD7h/kcH38me5ex8l
hR3Orc17eri/ncXdsRC2c61cN/cXd4fbxvU5j07PTXO3ua1ch/mKz72Kdq9f
PvV5MsdcYb/f2e9W8yLfIrP5cn6JX25fmrkXuiaZpSvTvdU1xjPiG2QmXxrv
D18fz4nP4PbPQpbEO6JbdB9eyzl9Op+uT5PF8a7wL/QLKHnmrqIY//g/Ur3n
piz12ZOeLZPj3r/TcxOyPJ4pXxF78y5+Kt6nuBd8ZUzzbn4m3s14N/hZrz/0
Pybkr7pT9u80TqIB9djDN8lmNKIVLdjLt8l21MV9RwdS+C7ZjTPoRQ/28ZPO
ZQoTSGAGqXync8xznjfMydZzXabxBVwh8nEQRUjnD8nDyEUJirGfL5VlSHGH
KmQ5MvgjshL7+GpZhUy+RtYilT8mjyKLPy5PII0/KeuQHfcnfn+s73Rch65B
1iNH1yib4jr4FtmMA/xKe/sJqu35Z/JT5PKr5Oeo4b+UXyCPXy3XoJZfJ9ci
n18vv8JR/hv5NQr4b+V3WGZdx3Tfe70BhbqNsj+eg/09ruuK3ytzMl+kn5MD
5hO6hPtQJwfN8/xB/YIcivch/sboT8XvAndItyhvjf3/O/sfX6sK1w==
       "]], Line3DBox[CompressedData["
1:eJwl1Hecz3UcB/Dvz0olKlsyQrYkFBXKSEJ2OBSibBUa9miYDdFElD2a9t57
73l3nH137s664+70/D388fy+Ht8/vo/H6/3+fL/fop36NOsdCoIgwqVKWhBU
5jkq8SwVeYYKlKccZSlDaUpRkqcpwdzCoaC4LMZTFKUIhSnEkxTkCQqQn3zk
JQ+5yUVOHucxHiUH2XmEbDzMQzxIVh4gC5nJREYyECIId9XrXmoQpJNGKne5
QwrJ3OYWN7nBdZJIJIFrxBNHLFe5wmUucZELnCeGc5wlmigiOcNpTnEyfF8o
FJyQxznGUY5wmEMc5AD72cde9rCbXexkB9vZxla2sJlNbGQD61nHWtawmlWs
ZAXLWcZSJsbYabp9kplMZCQDIQLu2Wk6aaRylzukkMxtbnGTG1wniUQSuEY8
ccRylUXO6Iq8zCUucoHzxHCOs0QTRSRnOM0pTnKC4xzjKEc4zCEOcoD97GMv
e9jNLnayg+1sYytbwu+vXpvlJjaygfWsYy1rWM0qVrKC5SxjKUtYzH/8yz/8
zV/8GZ6bhSxgPvPC3xFzmM0sZvIHvzOD6fzGtPBz+k2VU/iVX/iZn/iRH5jM
JL5nIt/xLd/wNRMYzzjGMobRfMWXfMHnjGIkIxjOMIYyhMEMCs+gz0D5GZ/y
CR8zgP704yM+5AP60ofe9KInPehON97nPbrShXfpTCc68g5v04H2tCOCtrSh
NW/Ripa0oDnNaEoT3qQxjWjIGzTgderzGvWoS53wGZmxtnyVV6hFTWrwMi/x
ItWpxgs8T9W0+//Z/wF+HhUx
       "]], Line3DBox[CompressedData["
1:eJwl1VdbFVcUBuB9EEF6ERVUrHf5E/kB6cbeRU25SZ40TdEUe0XsvSaxpGEX
K/beuyAiYO/Gjgp558nFYs16vz3n7HNm5tC24PMOn8VCCHf9uVEfwstrISTd
CqGoMoREvZGqZQl6QzWJN9CLW8dCvP5CFqfHVKEs6GuuhPCcr9brb4Ywkdfp
F6+HsPxGCOX6oPshfPEwhMF6C/MKnng7hIKqECrM3/Iv5V+p7xy3ZK2i/ViT
673zHU/3ut/LvrbmBz3fmjY8yZrW+gz5EP6NfKjeSn7BntrKZsp+ZGXm1rxc
32aexQdZf9HchlfopXw2H8wvRefzSj3Z+8zh7bzeZXM7XqWvsr8U2VxZe1k1
ay+r0av0y6r4XghXzEOte+b1n6tqvuxOCAeq7Y2/YLWqhh9/EMJBvlw+wHd0
x7mHzOu8zkDzXfNh80/Oe+mcV+oeW20vK5xzRPaz7DWvU/dlf/Cj/Bdez4L3
eMD/5Mf4MB5jceoh/4sf58N5Axav/uV/8xN8BG/IEtQj/g8/yUfyRNZIPebF
/A37OiUbJUviyeqJbKXsNB/NU1iqehp9p/wMH8PTWLp6Fn0+fpaP5RksUz3n
a/g5Po5nsWz1gq/l5/l43pjlqFq+jl/gE/ib9taEN1UvZetlZbKJsmYsV73i
G3g5L+R5rLl6zUv4RT6Jt2AtVR3fyCt4Ec9n9Wyg52Ke++Qj/ZLsKqvUU62Z
zz/ml83Xonsnut72liZbIPtEVs2uy2r0dL6Qf8qvmG/wq3oGX8RXuleumW/y
LLbdNc/UF8tWyR7Z02r9lrwx3yHP1pfI1/DH8rX6bXkO3ynfpfq7/y65R1u4
159En1c/H30O1tzxNOfn6edYc3ZWzzVP5U31hOh71fNkTaLnReWoKfLs6Hcn
ulZ6rjxLr1KZarI8XY+Prr3eTJ4WPYuqUg3wPZzxXv3sr4E1y/zGNLWmf/Rd
ywv00/K+8jh5batYWGpNE2v6ySqs6aufsqaPNTFrfpfnyPtEv2fy3vpJeW95
kP8mbyzvxcvlPfUT8l7yevOv8mx5D15m7q4fly/hi1WWbL9r1I0f44vYQpXJ
9/Gu/ChfwOarDe6HDNleWRfZEdk8Plel8z28Mz/M57DZKo3v5p34IT6LzVSp
fBfvyA/yGWy6SuE7+Yf8AJ/GpqpkvoN34Pv5FDZZJfHtvMTePpDtkxXxSaqR
rFT2Pt/LC9nE6Lefb+Pv8T18AhuvEvhW/i7fzcexsaoh38Lf4Tv5GDZa7XK8
mb/NtzsexUaqHY438Y329Jas1PzUvftEDXHdNsqWej638hHWD1fbomeZb9Z7
un5bomfevEnvYY63h/XRa5q7R/eYucRxTO8W3VP6BnOq61zmONl7d+XDvHYw
p/BkdcFxl+j+8D9zvfXnzZ3NdeZ15nPmTubX5rXRd+O+71j1///r/wAcyLwr

       "]], Line3DBox[CompressedData["
1:eJwl1VOQXVkABdDX49hOp9Oxbdu2bXTHtm3btm3btoYZezLJrFfzse6uuh/n
7HtQN7xVZK2IkEAg8Nbj3odA4C53uM2BsJDALXmTG1znGle5wmUucZELnOcc
ZznDaU5xkhMc5xhHOcJhDnEwOBf72cde4ph7j9zNLnayg+1sYytb2MwmNrKB
9axjLWtYzSpWsoLlLGMpS1jMIhaygPnMYy5zmM0sZhJXpxlyOtOYyhQmM4mJ
TGA84xjLGEYzipGMYDjDGMoQBjOIgQygP/3oSx9604ue9KA73Ygkgq5E062L
7EwnOtKB9rSjLW1oTSta0oLmNKMpTWhMIxrSgPrUoy51qE0talKD6lSjKlWo
TCUqUoHylKMsZShNKUqyR+cSsjjFKEoRClOIghQgP/nISx5yk4uc5CA72chK
FjKTiYxkID3pSEsaUpOKcFISRgpCSU4ykpKExCQiIQmIT7zgWSAOsYlFTGIQ
PbgXRCUKX/A5n/Epn/AxHxFCgA/vA4H3/Ms7/uFv/uJP/uB3fuNXfuFnfuJH
fuB73vId3xJqbb+RX/MVX/KG17ziJS94zjOe8oTHPOIhD7jPPe5yh9vc4iY3
uM41rnKFy1ziIhc4zznOcobT7NLtlDzJCY5zjKMc4TCHOMgB9rOPvexhd3AM
drKD7WxjK1vYzCY2soH1rGMta1jNKlayguUsY2lwTN2WyMUsYiELmM885jKH
2cxiJjOYzjSmMoXJTGIiExjPOMYyhtGMYiQjGM4whjKEwQxiIAOC+5oiJNBf
9qMvfehNL3rSg+50I5IIutKFznSiIx1oTzva0obWtKIlLWhOM5rShMY0oiEN
SG6N6st61KUOtalFTWpQnWpUpQqVqURFKlCecpSlDKUpRUlKUJxiFKUIhSlE
QQqQn3zkJQ+5yUXO4FnQL4fMTjaykoXMZCIjGUhPOtKShtSkIpyUhJGC0OA3
k4ykJCExiUhIAuITL/jO3HFlHGITi5jEIDrRiEoU3r1xL93/N7zmFS95wXOe
8ZQnPOYRD3nA/Q///6//A8Tm0Rs=
       "]], Line3DBox[CompressedData["
1:eJwl1WV8lWUcxvH7UNLdKdYLBem2xfaVSShlF9gtrdi+UUG6tzEWjLEx2Nhg
zShbxCCUZkE3fp+PL/773dfvus9z6rl3Oo8Z98DYWAih3J+SWCzM32e9P4Tr
DoXQsVMsDN4dwoPyrIoQJlWFMBsX2FPBdbHnDv1D1nP4yfq5uFBfyVWZh808
bujeEBbxJdalpthMsX8Yv5gvkzeZjWYqP5xfwnf1HNebR1xnxcEQHuOX8t24
bK+vOz6qS9M9rovTHZWHmL07Pd71RvBdD4cwEvdxcXtcU47H/XICdrFvGR6Q
E3G5OWidhBePhJCMOypDOMTleN5QHkIKd0mXijH5sK46rpCrYRrWwCN8LVwp
18R0vAzL+Tq4Sq6NGVgXK/j6mCnXw9XYACv5hpgld/P612AjuYpvjHs6xsJa
rrsuG5twR3VNMUfuwa/DZvIxvjnmyj35PGwhH+db4nq5F78BW8kn+NaYL/fm
C7CNfJJvi4VyH74I28mn+PZYLPflc31mJdYduNO6jlgq99NtxE7yGf5yLJP7
85uws3yWvwI3ywP4LXilfI6/CrfKA/ltONp3fMx3Pwq/PhDC91y8e6KH+6Sn
GcN/w//AJ/C9uDyvrTc+oftW96Numa4P19c8yU/nf+IT+X5cf/MUP99983N0
v/ADuM3ync7E07oF1oui82CS9AP1W6zv0j8T3ffWS80Sk6wfpN9qvc3cYH3R
+8tyX8fJa/A5jxnjva21Ho0XzHnzFf+sWc2Pku93/cEef876Sz6T32QyzEju
Pv3t+rPWX+hX8WUm3Yzg7tXfpj9j/bl+Jb/RlJotPocOu1wPb43eL35mTxmW
RPuwvb4Eb9F38tmWWn9qTzEW21OE7ewpwJvtKcRP9PlYpN+AbfV5eJN+PX6s
z8VC/Tpso8/GG/U5OE2/Fgv0a7C1PgvP+W4+0p3F1XIrPhPPyB/ypzFDbsmv
wlPydudoqu6kdTrXQrcST8hT+OOYJjfnV+AxeTJ/FFPlZnwKVsmT+EpMlpvy
SVghT+TLo/tGbsIn4hF5An84uv/kxnwCHpLH8wcx3+cZzzXSxeEB7gPdflwq
N+SX4D75fX4vLpYb8IvwX/k9/h9cKNfnF+Ae+V1+N86X6/HzcJf8Dr8T58p1
+Tn4t/w2/xcWeE2zuTq6a5zHP7lBzuMfeLVcm5+l3yG/5TG/R///uO34pvwb
ztTX4r7DX+U3+F9whlyTz/e9vs5lOSev4Qa5Br8eX43uff4VzJOr84Ve0zrr
XJOpy8FsU033sn0Z3DiMyWOjMyC/hEF+EdPlF/CSM/h8dAbkNJMa/c44F/c4
I6ej3yQuhUuOfpvku/kiz33KOlWXxC83w+WT5oRJ4Ye4ZiI/VD5uhpl4+bzn
m+68j/cbOAPn2buMv9Z9noAX9DP5CfqJVf//Zv8HzZXKkA==
       "]], Line3DBox[CompressedData["
1:eJwV1WMQJUcUBtB5a9u2bdu2bdu2bdu2zWxsO9lN1jbPq6nTX9X8mdt9u+6k
79C3QZ9QEAT3LEdDoeAIhznEQQ6wn33sJUgbCvbI3exiJzvYzja2soXNbGIj
G1jPOtayhtWsYiUrWM4ylrKExSxiIQuYTzTfnifnMofZzGImM5jONKYyhclM
YiITGM84xjKG0YxiJCMYzjCGMoTBDGIgA+hPP/rSh95cTxMKesme9KA73ehK
FzrTiY50oD3taEsbWtOKlrSgOc1oShMa04iGNKA+9ahLHWpTi5rUoDo31FRN
VqUKlalERSpQnnKUpQylKUVJSlCcYhSlCIUpREEKkJ985CUPuclFTnKQnWxk
JQuZyURGLuhhBpmedKQlDalJRUpSkJxkJCUJiUlEQhIQn3jEJQ6xiUVMYhCd
aEQlCpGJREQiEMISvA9CwTve8obXvOIlL3jOeTU/k095wmMe8ZAH3Oced7nD
bW5xk//5jxtc51/+4W/+4k/+4Hd+41d+4Wd+4kd+4Hu+41u+4Wu+4ku+4HM+
41M+4WOu8RFXucJlLnGRC+F9cY6znOE0pzjJCY5zjKMc4TCHOMgB9rOPvexh
N7vYyQ62s42tbGEzm9jIBtaHz8NdXSfXsobVrGIlK1jOMpayhMUsYiELmM88
5jKH2cxiJjOYzjSmMoXJTGIiExjPOMYyhtGMYmT47NU2Qg5nGEMZwmAGMZAB
9KcffelDb3rRkx50pxtd6UJnOtGRDrSnHW1pQ2ta0ZIWNKcZTWlC43C/3MlG
siENqE896lKH2tSiJjWoTjWqUoXKVKIiFShPOcpShtKUoiQlKE4xilKEwhSi
IAXIT77we3XllXnITS5ykoPsZCMrWchMJjKSgfSkIy1pSE0qUpKC5CQjKUlI
TCIS8lSvEsj4xCMucYhNLGISg+hEIypRiEwkIhKB8GMwBO/fB8E73vKG17zi
JS94zjOe8oTHPOIhD7jPPe5yh5PO5ba8xU2W3w+Ca+bPVa5wmUtcDM9IznOO
s5zhNKc4yQmOc4zwP/sD2LsLjw==
       "]], Line3DBox[CompressedData["
1:eJwl1Md3VWUYxeH3EnovoYqJqHQIPdT03hOwoIKEUAUFy9+gYwEFsXcRTKFD
KEEZOHCB0mzAwLIEy6KMmNCesxzs7Hf/9rn3fDnfd+6opo0LN6QiotGfzLRU
tF+JuHs1ouzfiNrfI1aad12PeO1mxFf8iP4eFn9HrOLN2Ja/Io7iKawTrcZb
8DfxY3ga9m1mKjrzNbpW3TbdcV0XrCutxdvwt/AO/Kz5DGX8FvE29g5lmqdg
75qz+Hv8Iez2HxF3aJT5Lr9HtzJS8bAcf0aMvWZt/NKNiEew8XIneRxP4xP4
o/gk3lmeyLvwyXx0ck/eVc7i3fhUPgafxrvLC/6L6MGny2PxGbynnIP34jPl
cfgsfsW6emO5uj48Gxuvm837ynl4Pz5HnoDP5f3lfHwAnydPxOfzgXIBPihZ
gzwJz+HpciE+mH/gGXXzfN/n3/Mh2Efm7uYe9KH5B55uf4bqPpHL7X0FfWw+
oxuGf2auxKroU/NMz384/oW5Gquhz82z8CZnJZt/Kdfi58x1ztIOebZ5Ls2h
nXKd/rx5srXXmxsow/fu0j3IX3EuLlvDy3wEH07tupG6Mp+7JB+Sh/FSeSg/
LD+gL5EvygfkIbxYHswPyiP0RfKv8j55hfWlm/ebv/MsfjH/THvkJt0g817z
T/xHajMvxwead5sv8PPUYm7EB5hbzef4Wepw9pbhJ3ih+9banwZ7VU/HsGd1
x3mBrkZXh9fSEWyp7ijP11XravBsa6zmh/El+naep6/SV+GVdBB7RneI5+oq
dRV4Oe3HntYd4Dm6Cl0ZXkp7sad0+/gCXbmuBC+mPdji5P+wx/N1ZboivJB2
657Udejm6Up1BXg+NXsWT+hOWXeb60p0eXguvWRvH9e14sX4DJ+dTi/ij+Et
eBE+DZtKG/FF+Db3adat4iuT3yf3WIivMC+n7fIG1zZgW5PetY18Gb2A1+On
rae//TnpN2drct7/iXhe1w/rS9/gbyTvCL4e74P1pq/xzck7ha/De2E96QS+
CX8O22mtWc71Kde8jp3mJc53afK+JOt0zQ7XvHrz/9/e+7NAxsw=
       "]], Line3DBox[CompressedData["
1:eJwV1GV411UYBuDffxu5kYPRA0Q6LEoQFUwEEZEOSYsYKSEoBiEooWAAKggS
UwmlS1BRYXTH6NG1YDQ47n24z3M+nfNc1znXW7pL72ZxoSAIOlmiw0NBfvKR
lzzkJhdRRJKTHGQnG1nJQgThhBEiICPMvmQo+F/e4y53uM0tbnKD61wjnauk
kUoKyVzhMpe4yAXOc46znOE0p0jiJCc4zjGOcoTDVNAnUR7iIAfYzz72sofd
7GInO9jONrayhc0ksImN/Me//MMG/uYv/mQ96/iDtaxhNatYyQqWs4ylLGEx
v/Mbi1iYea/eC+R8fuUXfiaeecxlDrP5iVnM5EdmMJ0f+J7vmMZUpvAt3/A1
XzGZSXzJF0xkAuMZx+d8xljG8CmjMzvGhoJRciQj+ISP+YgPGc4HvM8whvIe
QxjMIAbyLgPoTz/60ofexNGLnvSgO+/wNm/xJm/Qja50oTOd6MjrdKA97WhL
G1rTipa0oDmv0YxXacorNOFlGtOIl2jIi7zA8zxHgvd6Vj5DA+rzNE/xJPV4
grrU4XFqU4ua1KA6j/Eoj/AwD1GNqlShMpWoSAXKU46yPEgZHqA0pdikV0kZ
SwmKU4yiFKEwhYihIAWIJj/5yEsecpOLKCLJSQ6yk42sZCGCcMIIEZARMhO4
x13usFGv2/IWN7nBda6RzlXSSCWFZK5wmUtc5ALnOcdZznCaUyRxkhMczzzb
Xz0mj3KEwyRyiIMcYD/72MsedrOLnexgO9vYyhY2k0B8chCUMv9iKUFxilGU
IhSmEDEUpACZs/c+U5u/ZA==
       "]], Line3DBox[CompressedData["
1:eJwl1NdSFUEUheGmyrcwAaKYURAEVBRRwEDwHHJQQYKI4cbwCOrrmKuMZDCL
gJmkZS7Dnd75dXmxWXv9q2fOTO8e0tpPJU6mhBBS/clbkBIGP4XwZGlKWPgl
hMS3EDrmQzj9OYTLP0O4+TGEYfki2WJ1Br+K38FH8SXYUnUWv4734ffxVOwc
1vU7hKy5EKZ/WfchhHc0lc7QDXjyRwhpfIKm043YMlruXhn0z5KUkI0t169Q
OfpMulJt0q+iNa5dTWfdMxer49fwtXQtrad5eCNdxzfQ9bSJbsabaRZf9d0z
0Rb+r9/Nl7XqN2LVsmx6kC/AD9Ec/gC+ie7xvIexQlkb3cvnxveS59F9fDu+
RX6E7uc340l5Pq3gO/Ct8k5ayRfEd5MXxmfju/Bt8m5azW+J7yjfGp+DP4oX
yXvokDmMm2kC3yYf5evNtoZvoCN8rb5OFckf8E14A99M7/ON+ia1Xf6Yb8Vb
+IP0UfT6bvN9qj9KD/FlcW/kk3E/6BNZ0vWd8TpVKm/Be61P4A/lB2hH/F21
O85JfkxejY/Jq+gR2Rczqae74oyt6bGmUjZsTQVtlyVVibxWPo+38ZnO4nv9
fmtefQ1hjp/lV9DD8uV0jt8nfymf5af5DPqCX0Zn+L3x3Orf6tNpmnqn34NP
69/Es+8MjqkmcynHT5jFa3zCLEbxEdUoK5Mdl72K3xc2pBrwUrwXf4kPYgOq
Ht+NH8Nf4P1Yn6rDd+FTWI9vtkQ/GWeh36mf0H+2Z918Mf+c79Lv0I/rO/WD
9Jm6YR922rca+zbAP/AN99Mp71+MT9J7/Bh+l16zfgc+gd/mR/Bb9Eo8L/g4
nnSv5/S1ZxiSX5JdNLPz6sLv//97/gFRP4MT
       "]], Line3DBox[CompressedData["
1:eJwV1GW0lkUUBtD3o6UF6bz30t3dpQjSIgJSgoR0d6cSgjRIN0iIoKDSJSEo
Kd3d3bDvjz3Pj1mzzpyZNRPWrEOt9qEgCNIbUkQLBclJRlKS8BGJScTRdKHg
Q5mQBMQnHnGJQ2w+IBYxiUF0ohGVKIQIeBc1FLzlDa95xUte8JxnPOUJFdR9
LB/xkAfc5x53ucNtbnGTG1znGle5wmUucZELnOccZznDaU5xkv85wXGOcZQj
HOY//uUQB/mHA+xnH3v5mz3sjlxj77vkTnawnW1sZQub2cRf/MkfbGQDv/Mb
61nHr6zlF9awmlWs5GdWsJxlLGUJi1nEQhYwn3nMZQ6zmcVPzGQG05nGVKYw
mUlMpKR+fpQTGM8PjGMsYxjN93zHKEYyguEMYyhDGMwgBjKA/vSjL33oTS96
0oPudKMrXehMJzrSgfa0oy3f0obWtKIl39CC5pH3p4evZTOa0oTGNOIrGtKA
+nxJPb6gLp9Th9rUoiY1qE41PqMqVfiUynzCxxxSt5KsSAXKU46ylKE0pSLP
mBIUpxhFKUJhClGQAuQnH3nJQ25ykZMcZCcbWclCZjKRkQxEEE4C+wqTa+4F
QWFvsyAFyE8+jpnPK/OQm1zkJAfZyUZWspCZTGQkAxGEE0Z60pGWNKQmFSmJ
/HveA2ONkQ8=
       "]]}, 
     {GrayLevel[0.2], Line3DBox[CompressedData["
1:eJwV0ldaE2EARuEZCxZ6RykhBAgEJBASIJSAvTeW4ALE3gu4BNfgEtyFF67C
KxXpHfX9L77nzDmTTDKTJB8szD+Moyj6YaVkHP0sxNEv+2Tz6Tia0e7jEr+H
0/wuLvI7OMVv40d+Cyf5TfzAb2CRX8f3/BpO8Kqi69s77ao2rl3Bt/wyjvFL
+IZfxAK/gK/5eczzOXzFZ3GUl/Aln8Ecr3T9JXuhTWsj2hQ+55M4zIv4jE9g
lo/jUz6GQ7yAT3gez/FRfMxzOMhH8BEvm4yjCp81oA1rC1oWM3wI5zLej/18
EGf5APbxDJZ4P6Z5X3jmPI29vDc8c/7V9cutR+sJz13rxm6eCs+dd2GKJ8Nz
553YxRPhufMOTIb7dPx30H2G14T7xEOeynm9rfichN7pOGkHziWwI0xvw3bb
11vxbHC9Bc/Ynt6MTdamN2Cj7er1WGeteg3W2o5ejVXhWnoFVtq2Xo6nw3X1
b77XScenbMu5E1hmLc4dw+O2qR/FI9asRxjbhv5vxH1bk553z+vaIT+wxvD7
amvaPt+zBi2nrWq7fMfqtXbfY9FW9G1ty+rC/8Fr/2ibfMNqw/9BW9bW+ZrV
aFntt7bKV6w6/E+0LxPez5fts+Pv9h+4fXfU
       "]], Line3DBox[CompressedData["
1:eJwV0+laTlEchvG9zfM8DxGZ55DMQ0VFSspQaaLMlTILyVAJSZKkOAKzk3Dh
THzgCPzWh/3ez3Ovevd617r+iZW1+TVxFEWzfLxPjKOB5XG0elUc9WN1chy9
46rwrb6K78NT+gB/Et/oK/lerNT7+Qp8ra/ge7A8/D9fhq/05Xw3lup9/Al8
qS/ju7AkfC//e3McFcsvuKXWOrFI77V2HJ/rS/gOPBbeyR/FZ/pi/ike0Xv4
QnyiL+LbsSDshT+Mj/Ukvg3z9W7+ELbqC/kWzAt75HNxnn21eRZYe2TtINdl
LQcf6on8AzwQ9s7vx/v6fL4Zs/VOPgvv6fP4JswMv4nfh3f1BP4O7tU7+Ay8
rc/lGzE9/FY+DW/pc/h/9pXg2cM9tbYbb1qbbe0G7grnwO/E6/os/hru0Nv5
7XhVn8lfwW3hfPjL8gyuAbdybVy9PJ27hFvCeXF18jSuFr/bx2a+ha/Rp/IX
MZV7xF2Qp3DncVM4O+6cPJk7iyncA+6MPIk7jRvDWXLV8kSuCjdwzdwpeQJ3
Etdzf73/Hl+pj+crMJlv4srlcVwZrgvnzZXKY7kTuJa7w5XIY7hiXBPOnyuS
R3PHw6xwjdwxeRR3NMxJuA/ul/cf0UfyhWFOwl3wBfII7nCYEe4Gly8P5w6F
+Qh3w+XJw7jcMBvcNe6gPJTLCTMR7oo7IA/h9od54K5wP707Wx/MZ4V54C/z
mfIgbl+YBa6B2yvHXEaYA66eS5cjLg2bU9wx921DHH311MlfktwxtnrPn1T3
K3/mLuInvIAf8Tx+wHPY5+9+eP4DlSSb5g==
       "]], Line3DBox[CompressedData["
1:eJwV01Vzk0EYQOEE10Jxd3dti7VQtLi7++Du7vViFUrR4u7ufs/fwN15cnFy
9j2bSTY7X2pOnj94XjAQCER4eV4rGEirHwxkG1L5mTmFT7QNBpL5qTmJj5sT
+Yk5gY+Z4/mxeQfnmrfzI/M2Pmreyg/NpdoFA1k4om3RHmibeRPuW2/kDbhn
vZ7X4a71Wl6DO9areRVuW4f7rEysNN8yr+DluGm9jJfihvUSnu53zcQMXNem
8jRMwhRMxjV9Ak/EWIzHOFzVR/MYjEBY+2BgFI/EFXvDeDgGYyiG4LI+kAeh
HwagPy7pfbgveqE34nBR78E90RXd0Q0X9C4ci8PurTPH4LzeiaNxSG/nHgo5
V0fzOXvtuQMOhvY4Cmf1CI7EAb0tt8EZvRW3Ro7eklvgtN6Mm2O/3pSb4JTe
iBsjW2/IDXBSr8f1sU+vywWdKcrZTtirba6DLHu1uCaO69W5BjL1alwVx/TK
XAUZeiWuiFy9PFdAul6Oy+KoXprLYK9eisNxRC/BJbFHD+PiobvUI52ruPMV
M++2V5SLhO7TXiEujF16QS4Quks9H+fHTj0v5wndpR76EwWRpv/DX+Tov/kP
UvELP7Ff/84/kIJv+Ips/QsXc6YIZ0u2/oxP2GfvA39EEt7jHbL0xZ7vROtF
nGleyAnmBZxhns/x5nmcbp7LO8xzeK+5pO/KwHZttrZHm8XbzG/xBru1+238
N7Er9KxGeVYxJNJzj53aJb6INOt79u4i1foOp/v8WKSYY81dEee9tzhZ68xd
0Eu7yUlaNMeEvke7wYnaS36FHtp1TtCe8wt0165xvPaMX+M/PS3Kuw==
       "]], Line3DBox[CompressedData["
1:eJwV1GWUlVUUBuB7UVCQ7u5SRBGQsrsLuxW76O7u7u7uzoEJumPo7s6hhpTn
/Njr+d79njXruzP3TKFqNapWj0YikTwpIpGSRaKRK+WjkSQTUzoauVMuGrlt
Stgn85ZZYn+TN0xx+2u8bhbbX2WSKWZ/mef9nCtcpLvEdRWikYssqj/Pc/oL
XKg/x7X6syyiP82z+jNcoD/FNfqTLKzfXjkaOe75jDMnON+ZY1ztzFEWcuYw
T+uPcJ7+EFfpD7Kgfj9P6Q9wrn4fV+r3soB+N0/q93COfhdX6Hcyf3gHntDv
4Gx9IpfrtzGffguP67dyln4zE/SbmFe/gcf0Gznc55npzHrP8c6sYx5n1vCo
M2s5Q7+acfpVzK1fwSP6lZyuX85YfQJz6eN4WB/PafpYLtMvY059DA/pl3Kq
fgmX6hczh34hD+oXcYp+AWP085ldP5cH9PM4WX/L55jjubvdbGZzZib3y7M4
yZkZ7CZPZ1b9VO6Tp3Gifgq7ypOZRT+Re+VJnKCfwC7yeGbWj+UeeRzH68ew
szyamfQjuVsexXH6EewkD2dG/VDukodxrD7ZZxjiuaPdYGZwZiB3yoM4xpkB
7CD3Z3p9X+6Q+3G0vg/by72ZTt+T2+VeHKXvwXZyd6bVd2Wi3I0j9V3YVu7M
x/QduU3uxBH6Dmwjt2cafVtulduF75K+DVvLrXnT50ntTEvPW+xacZgzLdhK
bs5H9U25WW7GofombCk35iP6htwkN+IQfQO2kOszlb4uN8r1OFhfh83l2kyp
r8kNci0O0tdgM7k6H9b/y/Xyfxyo/4dN5b/5kD7R5/jT8zq7vzjAmT/YRP6d
KZz5lWvl39hf/wsby9UY1f/ENfLP7Kf/kY3kHxjRf8fV8vfsq/+WDeVv+H/h
aOQrrpK/Zh/9l2wgf8H7+s+4Uv6cvfWfsr5clff0H3OF/En4G/g8vZz5yHM9
uw9515n3uVz+gD3177Gu/C7v6N9mgvwOe+jfYh35Td7Wv854+Y1wD/Wvsbb8
Km/pX2ac/Eq4h/qXWEt+kcn65xkrvxDuof451pSrhO+SvhIrh++I9++ir+i5
grmhK89nw/2zL8ey5rp9GT4T7p3903zKXLN/kqXDXbMvxSfMVfuSfDzcMfsS
LG6S7IuyWLhb9kVY2FyxTyjjd2aqlPI/23vFe44zl3XLGGsq65YyxlyyX8wl
ppL9Ii40F+3nc4GpaD+Pc80F+9mcYyrYX/b7uGTO21/gxXD3vFMh71PQnLPP
zwLhHtvnY97wO/JuZ3W5PecJd1iXiznNGfvszBHurn02ZjWn7TMzS7iv9pmY
0ZyyT88M4Z7ap2Nac9I+Tfi/Ee6nfepwv80J+1ThLod7aZ8y3Dtz3H6b90oR
7lm4k7pouA/mmO5+Wd9909T+Hu+ao/a3ecc0Cf/3mWyO2N/gTTO+kr89r5lB
npN41SR6fgD7jps1
       "]], Line3DBox[CompressedData["
1:eJwV1FV8VnUYwPFzFEEx7/UGC1EkVMpuAUUFuxXr8/HCvLKD7u7u2EY3Y/TY
xoABY2wDxoAxYEO6bL/PxW/xfXbO+3/P+9nToPMXnT5PkyRZ68vFRmnyTcs0
+VqX/LxcNc3T5NMmaXJDszS5XhfYMh3jn/Dr2LU6z5bqKP+YX8Ou1jm2REf4
R7w+u0pn2WJV8w/5layezrBFOsw784o2aVKXX6HTfKGqzD4wq8Mu1ym2QIf4
+/wyluokm6+D/D2esP+apskJNk8H+Lv8X/aPfmdzVcnf4X+zv3SczdF+/jb/
k/2hfc5Wy7NUYfaW2SV+UTUsU/v4m/wCO69jLEN7+Rv8HDuro2y29vDX+Rl2
WkfYLJXz1/gpdlLVbKbK+Kv8RJxfh1l757pJpWavmB3ntaoym6Hd/GVeE+fR
ITZdJfwlfjReVwfZNO3inXh13F8H2FQV8468Ku6hSjZFO/mL/GD8rfazds5z
o3aYvWBWGa4Ks8nazp+Pzzuea8QmqYh34HvZHu1lE7WNP8fLWVnM2ARt5c/y
UrZb5Wy8tvD2vITtimvYOBXyAudqZ1bMd6qUj9Vms7Z8B9se92NjVMCf4UVs
m0rYaOXzp/lWtiVeh41SHn+KF7LNKmYjtYk/yQtYfrw2G6Fc/kS8Z+fK45u0
gw/XRrPHzXLZxjgXG6YN/DG+ga1XERuq9fxRvo6tjfOyIVrHH+Fr2GptZYO1
lj/Mc9iqeB9skNbwh3g2WxnP29kK+UCtNnvQbAVfHu+RDVAOf4AvY0tVwPpr
Fb+fL2GL472zfsrmbfgitlB5rK9W8tZ8AZsfz4P10Qreis9jc+OzdqZc3jt2
l1lLszk8K54V6xW7i7fgmSxDG1jP2F38Pj6bzYpnyHrE7uL38plshtax7rG7
+D18OpsWz5bVa+j17JjmfCqbojW8rXPVV5ZZM7PJfFI8d7O6rql1r0yzHL9/
Zfd+1ipNvvS9Ywt7z3yW701dN9E1E+Jz8Xd1eIZrmvDxbJyyebfYn+53Nx/L
xsTnxbrG7uSN+Wg2KnZO4zTpEruT5zvfXWYj+Yj4/zf7Lfan2Z18OBsWe4f/
GvuTN+JD2ZDYC/yX2J/8Dj6YDYqdw3+O/ckb8oFsQOwL/lPsTn4778/6xc7h
P8bu5LfxMufqy/vELjH7Ifan2a1mvVmv2Dv8+9if/Bbek/WI/cK/i/3Jb+bd
Wbf4XFqnybexP3kD3pV1if8Z/j8Y/m9q
       "]], Line3DBox[CompressedData["
1:eJwV0kVzFVEQBtAZ3C2B4BAgBBICcWSJOxR7WLDC3d3dCe7uXgUFfwB3CE5w
K9wdziy+19On33v3Vk3H9xzYbUAYBEGcj64pYZCdEAaxyWEQI0fSwyCnYRg0
Sw2DplKOvUoKgz18KW/CGktZ/pLv5kt4NsuSMvwF38UX80yWIaX5c76TL+Lp
LE1K8Wd8B1/IHzQNg1TeSEqaPTXbbrbArCFLkRL8Cd/G5/MGLFmK88d8K5/H
k1h9KcYf8S18Lq/HEqUof8g38zm8LkuQIvwB38Rn8zqstuS5W2GzPLONZrPM
avF4KcTv8w18Jq/JakhBfo+v5zN4dVZNCvC7fB2fzquyKpKf3+Fr+TRemVWS
fPw2X8On8oosTkLew71qymqzKWYVeHkJzG75zSo+mceyGPnHbspKPomXi96n
/GU3ZAWfyLPsxR99rixnE6L3yX7rr8syNp5lsF/6a/LFPXL4uOgd85/sqixl
Y1ka+6G/IkvYGJbKvusvy2I2mjVi3/SXZBEbFb1/9lV/URaykSyFfdFfkAXs
ovNHRDvBP7PzMp8PZ8msi30/x+axYdGOsM7sLJvLhrL6rBM7w+awIdHOsI7s
NJvNBrNE1oGdYrPYIHbU+d+ahEF7fpLP5AOj3fLdduwTO8wGsATWln1kh1j/
aNdYG/aBHWT9WG3Wmr1nB1jfaO9YK/aO7Wd92DFnf3V2S/6W7+O9ebzvtmBv
2F7WK9pP1py9ju6YHQYnMu0rO65m17KnnrPUamqmWlXt7v8zPFfxnK5WVtPU
SmqqWlHNcX6cmqv+B4Czs/o=
       "]], Line3DBox[CompressedData["
1:eJwV0VVTVlEAheF97A5KEAv1yhlnEMXA/A1e6L03Nna3iHfqb7C7A+yYAVss
QpESCRFUEBXbh4s171rvnn2+M99JmbNkZmYUQmiT3HFRyJGd46NwZ1QUZoyO
wm3cmRaFDsOjkI237On8TdxhR3wW3rCn8ddxux34bXjNnspfxa32v5QobMEr
9hT+Mm62//KbMNeezOfgRvsPnzopCq1yictwdhE3OPvtbD1esCfx53Gd/Ytf
i+fsifxZXGP/5FfjGXsCfxpX2T/4lXjKHs+fxBV2G78cT9jp/HFcZn/nl+IX
73SMG+fsKC7hvjnLxCP2WP4wLra/8ovwkJ3GH8SFdiu/AA/YY/j9OL/92fw8
3Gen8ntxrt3CZ/s++XYW5rV/B2zmC/CS75eV7j/Cv97vMz94YhSGT4jCCEmR
T9wgbqg+TIbIRy6ZG6QPlmRp4gZySfpASZRGLolL0AdIvHzgErlYPU5iJCEj
Cg38AL6f3V/6ynsugeut95FeUs/Fcz30ntJd6rg4rqveTbpILRfLddI7S0ep
4WK4SO8gof2d/PY7vmlkFKqxEd/iB6zCBqzE91iB9ViOe/xfXzyrTK/j3mAt
lmINvsZ3+AqrsQR3u9PiTrH+livCKizESnyJFfgCy/E57nKn2Z1nehn3FN9g
AZbiE3yNj/EVPsISnO3eQ71Yf4BFeB8L8R6+xLv4AvPxOc5yJ09/5Pf+A1If
ubY=
       "]], Line3DBox[CompressedData["
1:eJwVz0VTVXEAxuFzPoG6sQsFC3EcZ+xu3SqIQae02J1gkdLYhT2jbuzuxO5e
+AFcuvNx8c49v+d/z7nnRqQVzSoMgyD4YyVDw6AxKgwih4dBqeucvmEQEx0G
DWz3sDBYqAfqet2ss3W0rtNNOksP0LW6UWfq/rpGN+gM3U/v0u1GeQ5LZ31Z
NavTabqPrtK1OlVH6Updo1N0pK7Qu3Sy7q3LdbVO0r10mW7rN6pYIotgO1ml
TtA99Q5doRfoHnq7LtfzdXe9TZfpebqb3qqLY8Jgru6qS/UinTkyDNpYF1bC
ili873TWW3ShnqM76c26QMfpjnqTztexuoPeqPP0bN1eb9C5eoDn/7bzg8Pg
zKAwWM9nus5xds7nWTvN1/EZrhfyX0PC4Ked4mv5dJ7Nf7DvdpKv4dN4Fv/G
vtoJvppP5Zn8C/tsx/kqPoVn8E/so730XsecrXQ22Vm6sw/8vbXwFXwST+Pv
2Fs7ypfziTyVv2Gv7QhfxifwFP7q//PtMF/Kx/Nk/oK12iG+hI/jSfw5e2YH
+S/v9XdEGIx1lujsKX9iB5wtds8YnsAfs0e2nxfz0XwBf8ge2D6+iM9n9/U9
28uK2Dx2V9+xPayQzWW39S3bzQrYbe8Sz2+yG9bM8/kcdl1fsyaWx+LYVX3F
Glkui2WX9SVrYDmsxX+7qC9YPWvV/wAZvaWQ
       "]], Line3DBox[CompressedData["
1:eJwV1NVCVVEUQNFz7O7uTmyxu1vAbsXuF7t+wPoPA7sTsLu7u7sLx36Y3DXn
As5mc6Fq6tyUOXEURc19eJwYR9NrxNGKFnFUsXYcNagXR9P4cl6BJ/CpfBkv
z+vzKXwpL8fr8cl8CS/L6/JJfDEvw+vwibxwmzhapJXWamup2kJeitfiE/gC
XpLX5OP5fF6C1+Dj+CtnLc6r87G8mLmaeYw5uXUcFeVV+WhexFzFPCo821zZ
PNL8yPdIbhlHSaikjdCyPGeO/hAPUFEfrv/TZ/P7uIcK+jD9rz6LJ3nmXa93
UN5uqN0fu5n8Nm6hnD5E/63P4DdxA2X1wfovfTq/jmsoow/Sf+rT+FVcQWk9
Rf+hT+WXcQml9GT9uz6FJzpTPnd90VzSLsnum91kfgHnUUIfqH/VJ/FzOBvu
Vh+gf9En8jM4jWJ6f/2znspP4SSK6v30T/oEfgLHUUTvq3/Ux/NjyOtMzZ2t
sF0fuw924/SjOIJCem/9vT6WZyIDBfVe+jt9DE/HYRTQe+pv9dH8EA4iv95D
f6OP4gewH/n07vprfSTfh73hXHo3PafzvbIboe3BbuSx62r3Uh/Od2Encutd
9Bf6ML4D25FL76w/14fybdiKnHon/Zk+hG/BZuTQO+pP9cF8E9KQXe+g53Cm
tu5skLYRG5DNrr3dE1+TwtdjHWK9XThruEvv94i3DWcMd8iz6sZRm3C2cHf8
H28dzhTujP/lrXh+z3wa7kv7o7UMzwr3xH/zFjyP+Zc5MbxPvb5M8Ls1zwt/
L/6uZmGzfy7NtLRWfr9oah7oZ+mMTN6Er/V528wZaMzX8K3mdDTiq/kW82E0
5KvC9zUfQgO+km8yH0QCv+j1P4Z6ykY=
       "]]}, 
     {GrayLevel[0.2], 
      Line3DBox[{73407, 68462, 69110, 68461, 70846, 69107, 70845, 80200, 
       68460, 70127, 68459, 70844, 69089, 70843, 68451, 69088, 68450, 72140, 
       73622}], Line3DBox[CompressedData["
1:eJwV0ltoz2Ecx/HnyRzaNNtCZmytTRZiU9IWI+fDhUMpCrMZUg4JF+aCUCZh
bixujIwo2swFan+iHKIYN0YO5VAsUyjnw+t38e39/Xx+z/P5fX/P7yms2bhw
QwwhfFfX82M4pT5lxVA6MIbBmTFk5cSQh2V0Me7IiKEEd+L73Bg+qHp9KW8f
zi+PoYuXqxroSv5h/JYXw1e5M+RM4hXJnYwz6dnY0T+Go9bN0x9L1ts/VP1R
TXQV/wT+kPNTzmL7VvBK5VTjEvovvxbX0eX89Vhifb6MXhUxNNu/mXcaD/WM
4Txuoy9g2uAYCq0rs76NXjM0hkvJd9H7rRmLl+nxeAWfmreB/xLb6SP6FL6i
G/Vv8AbdpL+JPbNjqDPbZ+dzgjfPfCdxO6+PZ7vwLL2Ifw43ef/FITG0qNH+
SStu5m1RLcn3WteK6fbutfeifimvDevpTP5vs7bTVfwUTqSLfOOyQc7KHOPk
LtfXqMme3bKm1trb+IfOktEg6w69mn83+Zd0Dv+Bfi3vIdaZ6Zr52lWFzBTu
5u1RTzzfYF0nHrd3gL1N+JTeyn+Gy5Mz4b3W33NW7/A+ZvknxeYdhI/oj/zH
yb1ztt36bnxCzzZrJ87Cf/w5+JxegC+SO4klcorlvKF794vhLR7o5f4m76G7
cJU58szXZpZs3kHz5eAluoCfwny6kV+A1+jhuFT+F/tH6b/iCO8Zldwl/ElP
4P9K7rR5h8npsG8ir1lOJT6ip+NKOT36xjBXn4ZT7J+qMt3ddHohPwOn8crk
Z+ureTk4Uu5HOTX0VbkrsZsewz/jn1QVxPAfuKiRzg==
       "]], 
      Line3DBox[{73411, 67849, 66995, 72872, 73258, 72871, 72873, 73060, 
       69251, 70214, 69250, 72699, 69889, 72698, 69249, 70209, 69248, 73059, 
       73058, 73057, 73056, 73594}], 
      Line3DBox[{73414, 67933, 67086, 80214, 67922, 67075, 73589}], 
      Line3DBox[{73415, 73341, 67123, 67958, 67110, 67946, 67098, 67934, 
       67087, 80215, 67923, 67076, 67912, 67065, 67902, 67054, 67892, 67041, 
       73590}], 
      Line3DBox[{73417, 71420, 69493, 70910, 68635, 69501, 68636, 72264, 
       71989, 72273, 68644, 69519, 68645, 70918, 69939, 71445, 80377, 69520, 
       70356, 69521, 71446, 69529, 72563, 68655, 69536, 68656, 73418}], 
      Line3DBox[{73421, 73379, 71621, 73372, 73631}], 
      Line3DBox[{73424, 68483, 69163, 67749, 80395, 66886, 72172, 68492, 
       69175, 70601, 70602, 72734, 71738, 68499, 73047, 72483, 72484, 72855, 
       72179, 68501, 69194, 67772, 66912, 72487, 69206, 70185, 73337, 73598}],
       Line3DBox[CompressedData["
1:eJwVkE0ow2Ecx7+/g79pXhq1vMzbTe5TkmIuIksKua0lb/23cdtljYi8xEV5
u1FKucrV5GYph1kiDi47unBZ1j4O3z7P83l+z+956YwmJuImKUGe20wLtaZg
venFY1qtIL2mZVw/7jxguiB7+LFGU5gUqDtm3t1i2qRuiLq3BtM7GWc95jdt
4Zda6U0czliEV+yZYn2aqMp0w7yPHqfUDtPDi7vDnTHPcodRXCV7PeQRPw9d
+szBDmrz/3eAX7AHfsMBWIQj0HFM+/SZZOxjHIEB6MIumIRBuA5DcJveB4zD
jD+4wwrn37L/BDeD+8QlcRn+457M4go4v890TU2ceVOzKec1paj74z9KxOW9
l/xJkdo1/A/n/JI09SHekqXXIHygxy7ukLXqOtMGtU+4DHc4wtfgdnCvuFi7
qQyDkEET
       "]], Line3DBox[CompressedData["
1:eJwNjztIAmAUhc8dlHKwSBKhhzSp1KiDEKgt4RBBRLREUDj5Gl1ERQ3CSKiG
WiOCqLXdwa22BsFGF8mCBlvUoL7h8B3Of/7LvStH+Z2cSdpDb8umvNu0Pmf6
WDIN0JnDFCS/XTQFYH/KdE2W4q1Od4Nu12N6R9s+U8ZrOiHfJP+k84Ue6K/x
95EZq/CPGc9krmlTC55GTfP4F7wfdmAI9mAEfsMYHMMkdDhN5/zZxc/iD+EC
TMMALMAwrMIEbMIteAMP+LePv8Mfwwv2e8Jn8Zf4EbsX2X3CPb8ozU333DQm
r5AP6fygEv04t7xyUwy2mdEgu+LNNWOq0S2T5/ymfyawNKs=
       "]], 
      Line3DBox[{73428, 68002, 67170, 73100, 73101, 73102, 73103, 69365, 
       70267, 69366, 72702, 69912, 72703, 69367, 70271, 69368, 73104, 72906, 
       72904, 73278, 72905, 73583}], 
      Line3DBox[{73429, 69402, 70291, 70660, 70661, 67244, 68590, 69416, 
       68075, 67254, 72235, 68599, 69430, 70666, 70667, 67264, 70904, 68606, 
       72543, 71800, 71801, 71802, 70906, 68608, 69447, 72924, 72925, 67281, 
       80396, 72249, 68615, 69454, 68105, 67291, 73571}], 
      Line3DBox[{73430, 69412, 70299, 72911, 72912, 67253, 68598, 69429, 
       72917, 72918, 67263, 72245, 68605, 72542, 71797, 71798, 71799, 70905, 
       68607, 69446, 72922, 72923, 67280, 80279, 68614, 69453, 73344, 73572}],
       Line3DBox[{73431, 73391, 67352, 72935, 72934, 69531, 68651, 70922, 
       71820, 70702, 70701, 70351, 69511, 80284, 71435, 67337, 70693, 70692, 
       69510, 68640, 72266, 67327, 73520}], 
      Line3DBox[{73433, 70952, 71853, 70733, 70732, 70399, 69604, 71842, 
       72748, 70727, 70726, 69603, 68692, 73663}], 
      Line3DBox[{73434, 73386, 71911, 73376, 73679}], 
      Line3DBox[{73435, 68825, 70791, 70790, 66619, 67546, 70506, 68811, 
       73450}], Line3DBox[{73436, 72402, 72401, 72996, 72399, 72400, 73634}], 
      Line3DBox[{73437, 68924, 70809, 70808, 66690, 67595, 71054, 70534, 
       68919, 72994, 72993, 73541}], Line3DBox[CompressedData["
1:eJwV0s9Pz3EcwPH3OytSK/1i6YfFCiE30RBlM7Y2S7n4VR0cRCFkFnWy2czF
/B4XmwP5cfCj8gewoVxy4SDq4lSqzfwYj8/hued3r9fr/Xq9vq/vt6yts7Ej
hhBG8aU0hoqsGHoyYhgrjOErvqFKbILHkV4TwyKexBRq5GaKYpjGLLYsjGEr
ynJjqOc6bFcznB/DCBp8nvVuBkX4hd/YK37H3J8LYmj0Zjeq9GjmJvwRP8AH
US3eyi34J36YJzfE0M61cvtK9MNK+zzQs1vvlMUxRCwxay7Pw0Xxj3YaxSH1
q9SvxmXxDPn+4hjm81JvsjgbV+Q+q/+ElJwYjpl5HDvN7eITmCN+kqfsdIob
5Lr5NNLkerkv+Y7i9/XrNPso2t2/gydS1XI3zuClmgFkejvIQ3hht2co8Hs8
57y0GNrSYyjlDq7ks1zNPXqcw3n0og9HzL5uh2vI0/cG30SX+G0esfstHjPr
tRuO8xvOcYNc5OG7WAHno9yNirkEf8Xr3bEO2zDlVj+Se2W7v3wpKtSXJbdF
qviY3k1qm7EHF+zRb/4HezzkQjs+4se4JPeUn6BEfIAHcVX8FQ+hXM9JPSt5
hRnLUZn8l3kt1ovP2mkGy/R468073NNjmN/jrju12KUVO9RXe7cOa5L/kLtv
9HkTdsnV8uakd/JfyYxhv9i03/M/hip2wQ==
       "]], 
      Line3DBox[{73441, 69500, 68634, 72263, 71988, 72272, 68642, 69516, 
       68643, 70917, 69938, 71443, 80376, 69517, 70355, 69518, 71444, 69528, 
       72562, 68653, 69535, 68654, 72280, 71992, 73511}], 
      Line3DBox[CompressedData["
1:eJwV0l9M11UYx/FzYgmGSBD4B6VS9CKtFvxwiG46m23lqo3CEdOLpiawFf5b
2ebCjQpK5h8Wm0PWdESTC0bd5HQzmVzY3BTn/GmlG4KxsKlFxebUab2+F5/v
9/u8n+d8nud7zpm3YctbjTGE8JzHjadj6M6NoWxaDPmZMVQ8EUN5QQxX8mOY
PTOGdGUMs7yz5Url2mbF0KG+uiiGtykDX4yPzY1hNx4LmdK9KTEswBvnxPAR
/oH3+zSBF+NFepzUowGrp5t4Ib5ZbYFZhsXTxNXiHHGnOdJYBnZdrz5rr86I
4Vc6iz+Yah3Po/hlLE2D+CSei6/iU27eFJ3A7+BpPi/iO/zTdurHx/Cp6jv5
DPI4TT34NTwDP4DXmvcdOoTnLYvhYtKbz53sGPZjPyUe4nFxs/hH8ZBezdb2
8jtKH+M/4Fv1DWrf5PcGbcP78ftP2avpMdzW8xbV4d/id/FxfIj3OqwLG7E3
k/hO/j28v6EquXa5v/Gt+BfFMbTSK3gL3uK70vcu3xXuQCUtpeexLdiEfL11
Hby+omfwjXiDec/p/573JirAa/Db+mxQP2mWfer3Upbc63J/yL0r14btoYeP
x/Ayvt76AV4r/fcK+hdfimfnxXAMHxe/IC41V4rK6Co2H/vLfDU8P+f3GQ3h
M/ExvarwpoSZ5ZPkHsjlyPUknvbydzqOPYaNqu/Ef8NuUB9+LyuGZeZpx9ea
sZq68X/wYfWr+ZfLp6gLv4XvUvu9M2kXj4o/FPcm52+GT7FfsCPuwHJr68y0
ObkD+AX8Z55l+HpsHTXiZ/BXkz3mE9zZ/2b7xk/hX/OpwtPmvUS1+DE8zWch
nzU8XkvuAP4dfh6fg88z77O0Bs91b3vlyvi0mXMldli8SNwsXiI+KN6jV37y
v/xStBjfjw/yzMRLsPlUgrfiA3gejxHzXqcivAn/ks+jHGfpn16iJ/GdyX7a
n5PW3HXmWXwyk3sjt01uQv2fzvx/lP6qhQ==
       "]], 
      Line3DBox[{73443, 69600, 68691, 70940, 69950, 71490, 69601, 70398, 
       69602, 71491, 69621, 70951, 73362, 73667}], 
      Line3DBox[{73447, 69699, 70960, 70959, 73158, 72614, 72615, 69714, 
       72350, 72968, 72349, 73526}], 
      Line3DBox[{73448, 69815, 70979, 70978, 67535, 68297, 70785, 69807, 
       73464}], Line3DBox[{73449, 72651, 72650, 73174, 72648, 72649, 73682}], 
      Line3DBox[{73451, 73387, 68812, 70507, 67547, 66620, 70792, 70793, 
       68826, 73397, 73636}], 
      Line3DBox[{73452, 73388, 69025, 70834, 70833, 71113, 70831, 70832, 
       69021, 70830, 70829, 68999, 70826, 70827, 73489}], 
      Line3DBox[{73458, 70300, 69424, 69921, 71365, 71366, 70304, 72705, 
       72704, 73182, 72706, 72707, 70314, 73117, 73116, 72921, 73538}], 
      Line3DBox[{73463, 69698, 70753, 68251, 67471, 70964, 71894, 70965, 
       69713, 72967, 73350, 73677}], 
      Line3DBox[{73465, 73394, 69808, 70786, 68298, 67536, 70980, 70981, 
       69816, 73403, 73684}], Line3DBox[CompressedData["
1:eJwV0klIllEUxvF7G7VBLcTK1GalJHURklpQ0EQLo4gKKhdGm6JJWuSiwBZJ
QbppWqQuCgpKxWjWoHFRUBmUEg0LpTKqRUqi0PT7Fof/eZ57z7nnvfedVblv
w94YQugV/3JiWDAphrQJMVSkxpCKRdNiWJ4ZQ7J8PS8Jl9FJpTEUWFssH8Fb
ay1i5fQYVsg7M2Iow4uTYyjFl/RcvQfTYxgeH0MJbwiL1OfrMx/76Xz+T9yq
T678nrpsPKNPFnbQOfpcw1T6JD8F2+hM/qeSGAac8VWPWmt9mKv3PNErH2Nv
D86h85w7E9/T65z3Dlfj35QYVmI3vQq78Iee3/l9+Ir+Ju/EKeozxDP5J95T
7MFG83TjL/McNEeXvIGXbsbz+Ibew3+NR7NjqBHtWTHcFR3iBb/YezzHSWrq
1TyR71TzGFN5j3AH/RCr1d9Wt0hNhre5Jb+T6GWtwp52nKimFbfSLXhMz2Re
s3wL7ypW6dOibqE+zdgq9vMuW9tozyWsUTdWXTU20eX8G76zUX6IN8raA/fQ
QN/H0/jRvXXIT8nb8QNdL3+Lt+hi93wTC/EEfyFeowuwDZPd8UzvNRKv0ofs
uYJ1o/1jWEVfwOvmOJC4F3PsxsVm24WV9J+0GDZjBV3I346b6GH+kHOa1G9L
vB3+dlaWGBTn6HL+Wew075rE/6FuKc7WZwmuoAf0GdSn3r4yXh32q58qxnmT
PqxNfGPi7vEL/VkckefxDif+Szw+I4b/7bmGfA==
       "]], 
      Line3DBox[{73474, 73395, 70138, 71179, 71178, 69862, 71160, 80329, 
       71161, 70126, 71159, 71158, 69859, 69087, 70122, 73033, 73542}], 
      Line3DBox[{73475, 70182, 69188, 69877, 71229, 71230, 70174, 72688, 
       72687, 73181, 72685, 72686, 70169, 73046, 73045, 73044, 73603}], 
      Line3DBox[{73476, 71364, 70303, 70673, 68089, 67271, 71384, 73575}], 
      Line3DBox[{73479, 70362, 69526, 69940, 71433, 71434, 70350, 71432, 
       80374, 71431, 69937, 71429, 71430, 70346, 73396, 73559}], 
      Line3DBox[{73483, 70550, 71688, 71095, 69004, 70079, 71693, 71694, 
       71695, 70836, 68413, 69027, 73398, 73652}], 
      Line3DBox[{73484, 67695, 66822, 71156, 71157, 70125, 70580, 80345, 
       67701, 66829, 73367, 73619}], 
      Line3DBox[{73485, 67757, 66895, 71227, 71228, 70173, 70606, 73486}], 
      Line3DBox[{73488, 68357, 72435, 68358, 72434, 68921, 70812, 73352, 
       73639}], 
      Line3DBox[{73490, 73353, 70828, 69845, 71093, 69001, 70078, 69002, 
       71094, 69022, 70835, 68412, 69026, 73651}], 
      Line3DBox[{73492, 73364, 70961, 69701, 72616, 68748, 72629, 68749, 
       73675}], 
      Line3DBox[{73493, 71034, 71035, 70034, 70533, 67594, 66689, 71052, 
       71053, 73494}], 
      Line3DBox[{73495, 71089, 71090, 70076, 70555, 67645, 66758, 73648}], 
      Line3DBox[{73497, 71200, 69867, 70858, 68474, 69132, 68473, 72157, 
       71958, 72151, 68468, 69116, 68467, 70854, 69109, 71173, 80202, 69101, 
       70133, 69100, 71172, 69861, 72684, 68457, 69099, 68456, 73628}], 
      Line3DBox[{73500, 71692, 71691, 66759, 67646, 70556, 70077, 71092, 
       71091, 66752, 67640, 73658}], 
      Line3DBox[{73501, 73355, 71698, 72728, 71697, 71696, 70080, 69006, 
       71096, 71689, 73654}], 
      Line3DBox[{73502, 67406, 71485, 71486, 70396, 70730, 73348, 73670}], 
      Line3DBox[{73503, 71854, 70735, 70734, 70400, 69606, 71845, 72749, 
       71844, 71843, 73664}], 
      Line3DBox[{73505, 66814, 72142, 68453, 69092, 70575, 70576, 71721, 
       71163, 69093, 70129, 70583, 80347, 70584, 71726, 70848, 68464, 69112, 
       72836, 72837, 66839, 73615}], 
      Line3DBox[{73506, 66815, 72143, 72144, 69094, 71722, 67696, 66823, 
       71164, 71165, 70130, 70585, 80348, 67702, 66831, 70849, 70850, 69113, 
       72838, 67709, 66840, 72153, 73614}], 
      Line3DBox[{73508, 66876, 68484, 69164, 72852, 72853, 80388, 66887, 
       72173, 68493, 69176, 70603, 70604, 72735, 71739, 68500, 72485, 70607, 
       70608, 66903, 70859, 68502, 69195, 72856, 72857, 66913, 68508, 69208, 
       67780, 66924, 71247, 69209, 70188, 67789, 73597}], 
      Line3DBox[{73509, 71962, 72482, 69190, 70183, 69189, 71243, 69878, 
       71231, 69174, 70175, 69173, 72690, 69874, 72689, 69172, 70170, 69171, 
       66885, 80207, 73400, 73604}], 
      Line3DBox[{73510, 67328, 72267, 68641, 69512, 70694, 70695, 71815, 
       71436, 80285, 69513, 70352, 70703, 70704, 71821, 70923, 68652, 69532, 
       72936, 72937, 67353, 73553}], 
      Line3DBox[{73518, 72262, 72261, 71987, 72270, 72271, 69515, 70916, 
       70915, 71617, 71439, 80375, 71440, 70354, 71442, 71441, 73129, 72560, 
       72561, 69534, 72279, 72278, 71991, 73392, 73551}], 
      Line3DBox[{73521, 72277, 68650, 69530, 68649, 70921, 69527, 70920, 
       68648, 71454, 68647, 80358, 70919, 69509, 70914, 68639, 69508, 68638, 
       73560}], Line3DBox[CompressedData["
1:eJwV0k1IlFEUxvF7oZnSxgILKc3KpdQsU9QUoo2RoE0Q4S6h0kIHWhhSQUFS
UFJEREktEgs0bFW2iBZDCX1ooi0cog+jNkI0SZsKoX4uHs49/3Pe55z3cqs6
spmeGELIUcOWGI5WxHCEPiRjKC6OoXVtDA3lMdRTSX0Ms3jEP22KoVntzIYY
TtNL/G9RDKech0pjeF0WwyvK4b/w1PoYavWn+WynJ/h3fIZPNV5YHcMY9hVL
6r3OowIfwt5jl+Sj/EboJjaNfamzud5+tf12ztBltQm1Fb6d53lW/lQ+ZU7E
fuov0En8Ef5vXQxLa8yyUzll8TH8hf4FfJrHYWwYe471mTWPt2OD2LgdTmDv
sFbsKnZcPivfLT+/fCeVMdQ49zlXi91ipXhILBUPiFf4rHTe67yUiGGXuCjW
it/EbWLn5hiOURflsSqsyDlj3j3zJrEyrEU+xe82lsNS2B4s6+56aByL2DX1
Ev+colHs96oYHvjHHe5xQC3vnubortqi2kd3dRFP6k/QLXwB78VG3NOA/LN8
mEeTeV1mvbFHp3hObU4tz6NmueabXuwtlpbv09NG3dgE1q7+Z6OdqAN7hg3y
bcGn7DRJB/HH+AzPrTyavL1GasMf4jud6/QXePygZjzhDd9Xu8ErrdZv70b8
DnbB+/8Plw9yhg==
       "]], 
      Line3DBox[{73537, 72903, 69364, 70270, 69363, 69911, 69362, 70266, 
       69361, 73099, 73098, 73543}], 
      Line3DBox[{73539, 68138, 67329, 72268, 72269, 69514, 70696, 68144, 
       67338, 71437, 80339, 71438, 70353, 70705, 68150, 67345, 71822, 71823, 
       69533, 72938, 68156, 67354, 73552}], 
      Line3DBox[{73546, 73257, 72869, 72870, 69247, 70213, 69246, 69888, 
       69245, 70208, 69244, 73055, 73593}], Line3DBox[CompressedData["
1:eJwV1F9IkGcUx/HncWSuNDNj6bLSzUnDonS7yBYro1VjXphBf4xoWtkwpgVu
g20ENVwZsUUjsFVGBeuii11ssVFbKTSDIiidg2LOSsqwsk0pFqzV5734+Xi+
53d+z/F9XyyobaxqiCGE2X5UTIthx9QYjmbGUJIew96cGKbnxfDGxBhKqWdC
DN/otZXFMGlSDLv1W2gX7cOLzFS9HEPK6Bj6zd2ibfgU/Jn5L/kKnPnUKWur
3tOXYmjGc7EcOoNvxB/iA6kxTBwTwwZ1toxe9Vj1cnW6utUe3VgKVoGNwrbJ
KnJvtqwJ9J28BXrd8pr0xmEZNB/7b6y/y74/y+g200Uz8Uf4Vt7X1Gm8o6lV
ToFeu5zjyZ1YpK/xfHzl5BgO4OPnxnDlRfdj9+V8hXWq09R31NvVv6hfUPep
G91T6J4G54dUIec6/p5zC++T7Bj+pYfjYvgdH3RnHf4Ye0TlfJfxerOvyhnB
+jyXYefHdus1d8zOlWb+TnKoEf8Tb54SwyK82T5vee/t5uc552DDev9QPe9V
3v0yivEh8w9oE34J34d/4O5NVEcXsU77rPf7K/Lu8Q5SDf+IvTr095ip0X+f
1lE79puZNcn3YeYu/wCtNXNKr4V/eJTnh92mavxH/CP73cFn2H+WvX8wO9t5
DRvSe0AreE/yfiHjMt5v/hZV4seTbDt9rlft7tW0iuZ5puftc9uz/snMjeSZ
Uiv/OfwmfhIv45tDm91zGl9uNt8OVc5ltIf/FP4t72BaDJ+qv1fvTTLVTeoT
yTdhhx3YEbNtNNduDXob7fUJ/oe7e6gEX49X40vkL6Z3qDZ5x3L+z9XnPyTj
IFXiB/Eu+y7k65bRRYVyFuotlVOOL6D5VJ68T/53ZSyz0wk7lmC7sbexNvXr
6u3qN9UXZHVSlrw8vFRemZwOrJ1S8fF4AV6Kn8V+pUzsMxl/2TdXzi677qSn
Gd4pPos3T32Gt9cep52Ps3y7+qmyZurPoGIawurNHPa/6zlwhbOX
       "]], 
      Line3DBox[{73557, 71453, 67344, 68149, 70700, 70349, 71428, 80338, 
       71427, 67336, 68143, 73558}], 
      Line3DBox[{73573, 73406, 80278, 67279, 69443, 70315, 69442, 72709, 
       69926, 72708, 69428, 70305, 69427, 71368, 69922, 71367, 69426, 70301, 
       69425, 72541, 71977, 73578}], 
      Line3DBox[{73612, 69129, 68472, 72156, 71957, 72150, 68466, 69115, 
       68465, 70853, 69108, 71171, 80201, 69098, 70132, 69097, 71170, 69860, 
       72683, 68455, 69096, 68454, 72147, 71952, 73627}], 
      Line3DBox[{73613, 72155, 72154, 71956, 72148, 72149, 69114, 70852, 
       70851, 71180, 71168, 80330, 71169, 70131, 71167, 71166, 72682, 72680, 
       72681, 69095, 72146, 72145, 71951, 73389, 73626}], 
      Line3DBox[{73616, 66838, 72835, 72834, 69111, 68463, 70847, 66830, 
       70582, 80346, 70581, 70128, 69091, 71162, 71720, 70574, 70573, 69090, 
       68452, 72141, 66813, 73399, 73623}], 
      Line3DBox[{73640, 68925, 70811, 70810, 72433, 72431, 72432, 68920, 
       72077, 72076, 73645}], 
      Line3DBox[{73660, 71841, 68198, 67407, 71487, 71488, 70397, 70731, 
       68202, 67413, 73385, 73669}], 
      Line3DBox[{73661, 73393, 69599, 70939, 70938, 71489, 70945, 70946, 
       69619, 70948, 72945, 70947, 69620, 70949, 70950, 73668}], 
      Line3DBox[{73671, 70752, 68250, 67470, 71541, 71542, 70446, 70756, 
       73349, 73678}]}, {}, {}}},
   VertexNormals->CompressedData["
1:eJzs3Hk0Vl/8OPoiIhSlQYVSIhJCZcg7EUkZmzQglQaR0ECURKSSuRIyz/M8
[encoded data, 976980 characters]
4v8XqfY3JA==
    "]],
  Axes->True,
  AxesLabel->{None, None, None},
  BoxRatios->{1, 1, 0.4},
  DisplayFunction->Identity,
  FaceGridsStyle->Automatic,
  ImageSize->{369.4170859008568, 136.64280487285498`},
  ImageSizeRaw->Automatic,
  Method->{"DefaultBoundaryStyle" -> Directive[
      GrayLevel[0.3]], "RotationControl" -> "Globe"},
  PlotRange->{{-5, 0}, {-40, 40}, {0, 20}},
  PlotRangePadding->{
    Scaled[0.02], 
    Scaled[0.02], Automatic},
  Ticks->{Automatic, Automatic, Automatic},
  ViewPoint->{-2.1341362834353457`, -2.6105670838001074`, 
   0.28355215517461685`},
  ViewVertical->{0., 0., 1.}]], "Output",
 CellChangeTimes->{{3.6804479644985075`*^9, 3.6804480267771797`*^9}, 
   3.6807991878059683`*^9}]
}, Open  ]]
}, Open  ]]
}, Open  ]]
}, Open  ]]
}, Open  ]]
},
WindowToolbars->"EditBar",
WindowSize->{1218, 677},
WindowMargins->{{-8, Automatic}, {Automatic, -8}},
Magnification:>2.2 Inherited,
FrontEndVersion->"11.0 for Microsoft Windows (64-bit) (July 28, 2016)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 168, 2, 366, "Title"],
Cell[CellGroupData[{
Cell[773, 28, 248, 3, 158, "Chapter"],
Cell[CellGroupData[{
Cell[1046, 35, 120, 2, 152, "Section"],
Cell[1169, 39, 380, 7, 65, "Text"]
}, Open  ]],
Cell[CellGroupData[{
Cell[1586, 51, 91, 1, 152, "Section"],
Cell[1680, 54, 862, 21, 106, "Text"],
Cell[2545, 77, 224, 3, 65, "Text"],
Cell[CellGroupData[{
Cell[2794, 84, 539, 12, 66, "Input"],
Cell[3336, 98, 336, 9, 66, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[3721, 113, 98, 1, 152, "Section"],
Cell[3822, 116, 305, 4, 65, "Text"]
}, Open  ]],
Cell[CellGroupData[{
Cell[4164, 125, 284, 3, 152, "Section"],
Cell[4451, 130, 239, 3, 65, "Text"]
}, Open  ]],
Cell[CellGroupData[{
Cell[4727, 138, 242, 3, 152, "Section"],
Cell[CellGroupData[{
Cell[4994, 145, 102, 1, 106, "Subsection"],
Cell[CellGroupData[{
Cell[5121, 150, 602, 14, 108, "Input"],
Cell[5726, 166, 375, 10, 66, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[6150, 182, 110, 1, 106, "Subsection"],
Cell[6263, 185, 170, 2, 65, "Text"],
Cell[CellGroupData[{
Cell[6458, 191, 202, 5, 66, "Input"],
Cell[6663, 198, 188, 5, 66, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[6900, 209, 108, 1, 106, "Subsection"],
Cell[CellGroupData[{
Cell[7033, 214, 198, 2, 83, "Subsubsection"],
Cell[CellGroupData[{
Cell[7256, 220, 822, 21, 66, "Input"],
Cell[8081, 243, 3356, 94, 339, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[11486, 343, 319, 4, 83, "Subsubsection"],
Cell[CellGroupData[{
Cell[11830, 351, 1166, 28, 108, "Input"],
Cell[12999, 381, 2227, 60, 112, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[15275, 447, 313, 4, 83, "Subsubsection"],
Cell[CellGroupData[{
Cell[15613, 455, 1186, 27, 108, "Input"],
Cell[16802, 484, 2084, 58, 171, "Output"]
}, Open  ]]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[18947, 549, 222, 3, 106, "Subsection"],
Cell[CellGroupData[{
Cell[19194, 556, 851, 22, 66, "Input"],
Cell[20048, 580, 11762, 232, 147, "Output"]
}, Open  ]]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[31871, 819, 177, 2, 152, "Section"],
Cell[CellGroupData[{
Cell[32073, 825, 2039, 51, 399, "Input"],
Cell[34115, 878, 2573, 50, 804, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[36737, 934, 316, 4, 152, "Section"],
Cell[CellGroupData[{
Cell[37078, 942, 2351, 59, 441, "Input"],
Cell[39432, 1003, 2983, 62, 804, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[42464, 1071, 320, 4, 152, "Section"],
Cell[42787, 1077, 323, 7, 146, "Text"]
}, Open  ]],
Cell[CellGroupData[{
Cell[43147, 1089, 226, 6, 152, "Section"],
Cell[43376, 1097, 979, 25, 106, "Text"],
Cell[44358, 1124, 161, 2, 65, "Text"],
Cell[44522, 1128, 380, 10, 106, "Text"],
Cell[CellGroupData[{
Cell[44927, 1142, 222, 3, 106, "Subsection"],
Cell[CellGroupData[{
Cell[45174, 1149, 2109, 52, 1293, "Input"],
Cell[47286, 1203, 2642, 52, 760, "Output"]
}, {2}]]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[49986, 1262, 156, 3, 241, "Section"],
Cell[50145, 1267, 474, 11, 69, "Text"],
Cell[CellGroupData[{
Cell[50644, 1282, 98, 1, 66, "Input"],
Cell[50745, 1285, 280, 9, 66, "Output"]
}, Open  ]],
Cell[51040, 1297, 236, 6, 66, "Input"],
Cell[CellGroupData[{
Cell[51301, 1307, 205, 7, 66, "Input"],
Cell[51509, 1316, 376, 12, 66, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[51922, 1333, 506, 14, 66, "Input"],
Cell[52431, 1349, 224, 6, 66, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[52692, 1360, 751, 19, 66, "Input"],
Cell[53446, 1381, 571, 11, 66, "Output"]
}, Open  ]],
Cell[54032, 1395, 216, 4, 65, "Text"],
Cell[54251, 1401, 107, 1, 65, "Text"],
Cell[CellGroupData[{
Cell[54383, 1406, 160, 3, 66, "Input"],
Cell[54546, 1411, 565, 10, 95, "Message"],
Cell[55114, 1423, 239, 6, 66, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[55390, 1434, 164, 3, 66, "Input"],
Cell[55557, 1439, 366, 10, 66, "Output"]
}, Open  ]],
Cell[55938, 1452, 684, 21, 121, "Text"],
Cell[CellGroupData[{
Cell[56647, 1477, 438, 13, 66, "Input"],
Cell[57088, 1492, 2844, 67, 316, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[59969, 1564, 474, 11, 66, "Input"],
Cell[60446, 1577, 2892, 89, 316, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[63375, 1671, 487, 11, 66, "Input"],
Cell[63865, 1684, 1841, 42, 524, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[65743, 1731, 109, 1, 106, "Subsection"],
Cell[CellGroupData[{
Cell[65877, 1736, 97, 1, 66, "Input"],
Cell[65977, 1739, 267, 8, 66, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[66281, 1752, 809, 18, 108, "Input"],
Cell[67093, 1772, 677, 13, 66, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[67807, 1790, 164, 3, 66, "Input"],
Cell[67974, 1795, 615, 18, 92, "Output"]
}, Open  ]],
Cell[68604, 1816, 124, 1, 65, "Text"],
Cell[CellGroupData[{
Cell[68753, 1821, 2732, 67, 691, "Input"],
Cell[71488, 1890, 2722, 53, 786, "Output"]
}, Open  ]]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[74271, 1950, 93, 1, 153, "Section"],
Cell[74367, 1953, 338, 9, 125, "Text"]
}, Open  ]],
Cell[CellGroupData[{
Cell[74742, 1967, 96, 1, 153, "Section"],
Cell[74841, 1970, 256, 5, 107, "Text"]
}, Open  ]],
Cell[CellGroupData[{
Cell[75134, 1980, 137, 1, 243, "Section"],
Cell[75274, 1983, 257, 6, 66, "Text"],
Cell[75534, 1991, 781, 16, 67, "Text"],
Cell[76318, 2009, 469, 10, 67, "Input"],
Cell[76790, 2021, 223, 4, 107, "Text"],
Cell[CellGroupData[{
Cell[77038, 2029, 106, 1, 105, "Subsection"],
Cell[CellGroupData[{
Cell[77169, 2034, 860, 21, 111, "Input"],
Cell[78032, 2057, 3033754, 49774, 334, "Output"]
}, Open  ]]
}, Open  ]]
}, Open  ]]
}, Open  ]]
}, Open  ]]
}
]
*)

