(* Content-type: application/mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 6.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[     96957,       2229]
NotebookOptionsPosition[     93614,       2110]
NotebookOutlinePosition[     94021,       2128]
CellTagsIndexPosition[     93978,       2125]
WindowFrame->Normal
ContainsDynamic->False*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[TextData[StyleBox["APPENDIX B",
 FontWeight->"Bold"]], "Text",
 CellChangeTimes->{{3.3931674469262543`*^9, 3.393167452283793*^9}},
 TextAlignment->Left,
 FontSize->14],

Cell["\<\
Simultaneous determination of k and Tc with incremental n values.\
\>", "Text",
 CellChangeTimes->{{3.3932562596421833`*^9, 3.393256312572752*^9}},
 TextAlignment->Left,
 FontSize->16,
 FontWeight->"Bold"],

Cell["\<\
Programmed by: Mark D. Normand and Micha Peleg
Last modified: August 1, 2007\
\>", "Text",
 CellChangeTimes->{{3.393168181924315*^9, 3.393168198891725*^9}, {
   3.393241755676222*^9, 3.393241769788744*^9}, {3.393241858404467*^9, 
   3.393241870013986*^9}, 3.393327138150137*^9, 3.393773040473209*^9, {
   3.3938433404414883`*^9, 3.393843354875573*^9}, {3.393844380872176*^9, 
   3.393844382530843*^9}, {3.394188008788476*^9, 3.394188008958943*^9}, {
   3.394961051258519*^9, 3.394961052693719*^9}},
 TextAlignment->Left,
 FontSize->14]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
Assign the three parameters k, n and Tc their original values.  These are the \
values we are trying to recover.\
\>", "Text",
 CellChangeTimes->{{3.3931675076578712`*^9, 3.3931675129825773`*^9}, {
  3.393252746087364*^9, 3.393252752848744*^9}, {3.393328911293683*^9, 
  3.393328916410233*^9}},
 TextAlignment->Left,
 FontSize->14],

Cell[BoxData[
 RowBox[{
  RowBox[{"k", "=", "0.065"}], ";", 
  RowBox[{"n", "=", "0.8"}], ";", 
  RowBox[{"Tc", "=", "100."}], ";"}]], "Input",
 CellChangeTimes->{{3.393241591214342*^9, 3.393241619905489*^9}, {
   3.393242584467102*^9, 3.3932425879189653`*^9}, 3.393242672115157*^9, 
   3.393242705182542*^9},
 FontSize->14]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
Assign the final time values, tCycle1, tCycle2 and tCycle3, for the three \
survival curves.\
\>", "Text",
 CellChangeTimes->{{3.3933289214914227`*^9, 3.393328957210947*^9}},
 TextAlignment->Left,
 FontSize->14],

Cell[BoxData[
 RowBox[{
  RowBox[{"tCycle1", "=", "7."}], ";", 
  RowBox[{"tCycle2", "=", "7.5"}], ";", 
  RowBox[{"tCycle3", "=", "8."}], ";"}]], "Input",
 FontSize->14]
}, Open  ]],

Cell[CellGroupData[{

Cell["Define the three temperature profile functions.", "Text",
 CellChangeTimes->{3.393327174489436*^9},
 TextAlignment->Left,
 FontSize->14],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"T1", "[", "t_", "]"}], ":=", 
   RowBox[{"If", "[", 
    RowBox[{
     RowBox[{"t", "\[LessEqual]", " ", "4."}], ",", 
     RowBox[{"60.", "+", 
      RowBox[{"20.", "*", 
       RowBox[{"t", "^", ".7"}]}]}], ",", 
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{"t", "\[LessEqual]", "6."}], ",", 
       RowBox[{"60.", "+", 
        RowBox[{"20.", "*", 
         RowBox[{"4.", "^", ".7"}]}]}], ",", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"60.", "+", 
          RowBox[{"20.", "*", 
           RowBox[{"4.", "^", ".7"}]}]}], ")"}], "*", 
        RowBox[{"Exp", "[", 
         RowBox[{
          RowBox[{"-", ".9"}], "*", 
          RowBox[{"(", 
           RowBox[{"t", "-", "6."}], ")"}]}], "]"}]}]}], "]"}]}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.393773122681652*^9, 3.393773135377913*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"T2", "[", "t_", "]"}], ":=", 
   RowBox[{"If", "[", 
    RowBox[{
     RowBox[{"t", "\[LessEqual]", " ", "3."}], ",", 
     RowBox[{"60.", "+", 
      RowBox[{"30.", "*", 
       RowBox[{"t", "^", ".8"}]}]}], ",", 
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{"t", "\[LessEqual]", "6."}], ",", 
       RowBox[{"60.", "+", 
        RowBox[{"30.", "*", 
         RowBox[{"3.", "^", ".8"}]}]}], ",", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"60.", "+", 
          RowBox[{"30.", "*", 
           RowBox[{"3", "^", ".8"}]}]}], ")"}], "*", 
        RowBox[{"Exp", "[", 
         RowBox[{
          RowBox[{"-", ".7"}], "*", 
          RowBox[{"(", 
           RowBox[{"t", "-", "6."}], ")"}]}], "]"}]}]}], "]"}]}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{3.393773089732197*^9, 3.393773157210267*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"T3", "[", "t_", "]"}], ":=", 
   RowBox[{"If", "[", 
    RowBox[{
     RowBox[{"t", "\[LessEqual]", " ", "2.5"}], ",", 
     RowBox[{"60.", "+", 
      RowBox[{"40.", "*", 
       RowBox[{"t", "^", ".85"}]}]}], ",", 
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{"t", "\[LessEqual]", "6."}], ",", 
       RowBox[{"60.", "+", 
        RowBox[{"40.", "*", 
         RowBox[{"2.5", "^", ".85"}]}]}], ",", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"60.", "+", 
          RowBox[{"40.", "*", 
           RowBox[{"2.5", "^", ".85"}]}]}], ")"}], "*", 
        RowBox[{"Exp", "[", 
         RowBox[{
          RowBox[{"-", ".5"}], "*", 
          RowBox[{"(", 
           RowBox[{"t", "-", "6."}], ")"}]}], "]"}]}]}], "]"}]}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.393773092609137*^9, 3.39377311505787*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
Compute and save the three temperatures at the three tCycle times using the \
three temperature profile functions.\
\>", "Text",
 CellChangeTimes->{{3.393322237848249*^9, 3.3933222418990507`*^9}, {
   3.393325056163465*^9, 3.3933250652433043`*^9}, {3.393325526678794*^9, 
   3.393325542925816*^9}, {3.393326043888006*^9, 3.393326044481868*^9}, 
   3.393327194632271*^9, {3.393327561595056*^9, 3.393327569871024*^9}, {
   3.393773046546372*^9, 3.39377304796911*^9}},
 TextAlignment->Left,
 FontSize->14],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"tCycle1T1", "=", 
   RowBox[{"T1", "[", "tCycle1", "]"}]}], ";", 
  RowBox[{"Print", "[", 
   RowBox[{"\"\<T1[tCycle1]=\>\"", ",", "tCycle1T1"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.393322100432173*^9, 3.3933221278763723`*^9}, 
   3.39332554916901*^9, {3.393325607619711*^9, 3.393325612349854*^9}, 
   3.393326084597148*^9, {3.393327584308352*^9, 3.393327634860392*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"T1[tCycle1]=\"\>", "\[InvisibleSpace]", "45.85305487675453`"}],
  
  SequenceForm["T1[tCycle1]=", 45.85305487675453],
  Editable->False]], "Print",
 CellChangeTimes->{3.393320121581217*^9, 3.3933221852294617`*^9, 
  3.393323123591468*^9, 3.3933252890246973`*^9, 3.3933254396472607`*^9, 
  3.3933257653010883`*^9, 3.393326611064291*^9, 3.393327059591311*^9, 
  3.393327201570022*^9, 3.3933292484726152`*^9, 3.3933296520516663`*^9, 
  3.393773676103527*^9, 3.39377482917377*^9, 3.393775157462851*^9, 
  3.3937761979996357`*^9, 3.3938432909486427`*^9, 3.3938445217085533`*^9, 
  3.394188019496896*^9, 3.3941883450942993`*^9, 3.394188779957304*^9, 
  3.394189225299621*^9, 3.394189676157607*^9, 3.394189827876747*^9, 
  3.39418990707544*^9, 3.394190186204294*^9, 3.394194819754313*^9, 
  3.394961072366146*^9, 3.3949624869469023`*^9, 3.3949647243377132`*^9, 
  3.394966933020776*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"tCycle2T2", "=", 
   RowBox[{"T2", "[", "tCycle2", "]"}]}], ";", 
  RowBox[{"Print", "[", 
   RowBox[{"\"\<T2[tCycle2]=\>\"", ",", "tCycle2T2"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.39332215124137*^9, 3.3933221528608713`*^9}, {
   3.393325557500794*^9, 3.393325579892871*^9}, 3.393325619031274*^9, 
   3.39332609024548*^9, {3.393327586174065*^9, 3.39332763633394*^9}, {
   3.3933296279823647`*^9, 3.39332963268054*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"T2[tCycle2]=\"\>", "\[InvisibleSpace]", 
   "46.278126718300705`"}],
  SequenceForm["T2[tCycle2]=", 46.278126718300705`],
  Editable->False]], "Print",
 CellChangeTimes->{3.393320121694007*^9, 3.393322185357984*^9, 
  3.393323123750484*^9, 3.393325289185557*^9, 3.393325439912942*^9, 
  3.3933257654491453`*^9, 3.393326611145443*^9, 3.3933270598144608`*^9, 
  3.393327201651826*^9, 3.393329248584855*^9, 3.39332965217715*^9, 
  3.393773676200733*^9, 3.39377482927818*^9, 3.393775157600686*^9, 
  3.393776198085019*^9, 3.3938432910494003`*^9, 3.393844521792973*^9, 
  3.3941880195918922`*^9, 3.394188345214879*^9, 3.394188780092585*^9, 
  3.394189225433771*^9, 3.394189676279735*^9, 3.3941898279913263`*^9, 
  3.3941899071790667`*^9, 3.394190186317499*^9, 3.394194819830782*^9, 
  3.3949610724456663`*^9, 3.394962487068411*^9, 3.3949647243966417`*^9, 
  3.394966933100456*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"tCycle3T3", "=", 
   RowBox[{"T3", "[", "tCycle3", "]"}]}], ";", 
  RowBox[{"Print", "[", 
   RowBox[{"\"\<T3[tCycle3]=\>\"", ",", "tCycle3T3"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.393322170793523*^9, 3.393322172397147*^9}, {
   3.393325566670063*^9, 3.393325582959031*^9}, 3.393325626854475*^9, 
   3.3933261158418827`*^9, {3.393327587116534*^9, 3.393327637340363*^9}, {
   3.393329635081936*^9, 3.393329637480438*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"T3[tCycle3]=\"\>", "\[InvisibleSpace]", "54.1365314659016`"}],
  SequenceForm["T3[tCycle3]=", 54.1365314659016],
  Editable->False]], "Print",
 CellChangeTimes->{3.3933201217909107`*^9, 3.393322185505349*^9, 
  3.393323123897531*^9, 3.393325289315895*^9, 3.393325440102467*^9, 
  3.3933257655674067`*^9, 3.393326611225205*^9, 3.393327059898323*^9, 
  3.393327201743946*^9, 3.3933292486867743`*^9, 3.3933296522775106`*^9, 
  3.3937736762990932`*^9, 3.393774829395447*^9, 3.3937751577093563`*^9, 
  3.393776198165964*^9, 3.393843291145875*^9, 3.393844521873609*^9, 
  3.394188019684207*^9, 3.3941883453082027`*^9, 3.3941887802040777`*^9, 
  3.394189225539753*^9, 3.394189676423441*^9, 3.3941898280727262`*^9, 
  3.3941899072727547`*^9, 3.3941901864166594`*^9, 3.3941948198980503`*^9, 
  3.3949610725379066`*^9, 3.394962487187297*^9, 3.394964724463957*^9, 
  3.394966933191948*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Define the three log survival functions as function modules.", "Text",
 TextAlignment->Left,
 FontSize->14],

Cell[BoxData[
 RowBox[{
  RowBox[{"LS1", "[", 
   RowBox[{
    RowBox[{"k_", "?", "NumericQ"}], ",", 
    RowBox[{"n_", "?", "NumericQ"}], ",", 
    RowBox[{"Tc_", "?", "NumericQ"}], ",", 
    RowBox[{"tFinal_", "?", "NumericQ"}]}], "]"}], ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"b1", ",", "sfn1", ",", "t", ",", "T1", ",", "y1"}], "}"}], ",", 
    "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{"T1", "[", "t_", "]"}], ":=", 
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{"t", "\[LessEqual]", " ", "4."}], ",", 
        RowBox[{"60.", "+", 
         RowBox[{"20.", "*", 
          RowBox[{"t", "^", ".7"}]}]}], ",", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{"t", "\[LessEqual]", "6."}], ",", 
          RowBox[{"60.", "+", 
           RowBox[{"20.", "*", 
            RowBox[{"4.", "^", ".7"}]}]}], ",", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"60.", "+", 
             RowBox[{"20.", "*", 
              RowBox[{"4.", "^", ".7"}]}]}], ")"}], "*", 
           RowBox[{"Exp", "[", 
            RowBox[{
             RowBox[{"-", ".9"}], "*", 
             RowBox[{"(", 
              RowBox[{"t", "-", "6."}], ")"}]}], "]"}]}]}], "]"}]}], "]"}]}], 
     ";", 
     RowBox[{
      RowBox[{"b1", "[", "t_", "]"}], ":=", 
      RowBox[{"Log", "[", 
       RowBox[{"1.", "+", 
        RowBox[{"Exp", "[", 
         RowBox[{"k", "*", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"T1", "[", "t", "]"}], "-", "Tc"}], ")"}]}], "]"}]}], 
       "]"}]}], ";", 
     RowBox[{
      RowBox[{"sfn1", "[", "t_", "]"}], "=", 
      RowBox[{
       RowBox[{"NDSolve", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{
            RowBox[{
             RowBox[{"y1", "'"}], "[", "t", "]"}], "\[Equal]", 
            RowBox[{
             RowBox[{"-", 
              RowBox[{"b1", "[", "t", "]"}]}], "*", "n", "*", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{"-", 
                 RowBox[{"y1", "[", "t", "]"}]}], "/", 
                RowBox[{"b1", "[", "t", "]"}]}], ")"}], "^", 
              RowBox[{"(", 
               RowBox[{
                RowBox[{"(", 
                 RowBox[{"n", "-", "1."}], ")"}], "/", "n"}], ")"}]}]}]}], 
           ",", 
           RowBox[{
            RowBox[{"y1", "[", "0.", "]"}], "\[Equal]", 
            RowBox[{"-", ".00001"}]}]}], "}"}], ",", 
         RowBox[{"y1", "[", "t", "]"}], ",", 
         RowBox[{"{", 
          RowBox[{"t", ",", "0.001", ",", 
           RowBox[{"tCycle1", "+", "1."}]}], "}"}]}], "]"}], "[", 
       RowBox[{"[", 
        RowBox[{"1", ",", "1", ",", "2"}], "]"}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"sfn1", "[", "tFinal", "]"}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{3.393773219420247*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"LS2", "[", 
   RowBox[{
    RowBox[{"k_", "?", "NumericQ"}], ",", 
    RowBox[{"n_", "?", "NumericQ"}], ",", 
    RowBox[{"Tc_", "?", "NumericQ"}], ",", 
    RowBox[{"tFinal_", "?", "NumericQ"}]}], "]"}], ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"b2", ",", "sfn2", ",", "t", ",", "T2", ",", "y2"}], "}"}], ",", 
    "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{"T2", "[", "t_", "]"}], ":=", 
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{"t", "\[LessEqual]", " ", "3."}], ",", 
        RowBox[{"60.", "+", 
         RowBox[{"30.", "*", 
          RowBox[{"t", "^", ".8"}]}]}], ",", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{"t", "\[LessEqual]", "6."}], ",", 
          RowBox[{"60.", "+", 
           RowBox[{"30.", "*", 
            RowBox[{"3.", "^", ".8"}]}]}], ",", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"60.", "+", 
             RowBox[{"30.", "*", 
              RowBox[{"3", "^", ".8"}]}]}], ")"}], "*", 
           RowBox[{"Exp", "[", 
            RowBox[{
             RowBox[{"-", ".7"}], "*", 
             RowBox[{"(", 
              RowBox[{"t", "-", "6."}], ")"}]}], "]"}]}]}], "]"}]}], "]"}]}], 
     ";", 
     RowBox[{
      RowBox[{"b2", "[", "t_", "]"}], ":=", 
      RowBox[{"Log", "[", 
       RowBox[{"1.", "+", 
        RowBox[{"Exp", "[", 
         RowBox[{"k", "*", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"T2", "[", "t", "]"}], "-", "Tc"}], ")"}]}], "]"}]}], 
       "]"}]}], ";", 
     RowBox[{
      RowBox[{"sfn2", "[", "t_", "]"}], "=", 
      RowBox[{
       RowBox[{"NDSolve", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{
            RowBox[{
             RowBox[{"y2", "'"}], "[", "t", "]"}], "\[Equal]", 
            RowBox[{
             RowBox[{"-", 
              RowBox[{"b2", "[", "t", "]"}]}], "*", "n", "*", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{"-", 
                 RowBox[{"y2", "[", "t", "]"}]}], "/", 
                RowBox[{"b2", "[", "t", "]"}]}], ")"}], "^", 
              RowBox[{"(", 
               RowBox[{
                RowBox[{"(", 
                 RowBox[{"n", "-", "1."}], ")"}], "/", "n"}], ")"}]}]}]}], 
           ",", 
           RowBox[{
            RowBox[{"y2", "[", "0.", "]"}], "\[Equal]", 
            RowBox[{"-", ".00001"}]}]}], "}"}], ",", 
         RowBox[{"y2", "[", "t", "]"}], ",", 
         RowBox[{"{", 
          RowBox[{"t", ",", "0.001", ",", 
           RowBox[{"tCycle2", "+", "1."}]}], "}"}]}], "]"}], "[", 
       RowBox[{"[", 
        RowBox[{"1", ",", "1", ",", "2"}], "]"}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"sfn2", "[", "tFinal", "]"}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{3.393773245564724*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"LS3", "[", 
   RowBox[{
    RowBox[{"k_", "?", "NumericQ"}], ",", 
    RowBox[{"n_", "?", "NumericQ"}], ",", 
    RowBox[{"Tc_", "?", "NumericQ"}], ",", 
    RowBox[{"tFinal_", "?", "NumericQ"}]}], "]"}], ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"b3", ",", "sfn3", ",", "t", ",", "T3", ",", "y3"}], "}"}], ",", 
    "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{"T3", "[", "t_", "]"}], ":=", 
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{"t", "\[LessEqual]", " ", "2.5"}], ",", 
        RowBox[{"60.", "+", 
         RowBox[{"40.", "*", 
          RowBox[{"t", "^", ".85"}]}]}], ",", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{"t", "\[LessEqual]", "6."}], ",", 
          RowBox[{"60.", "+", 
           RowBox[{"40.", "*", 
            RowBox[{"2.5", "^", ".85"}]}]}], ",", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"60.", "+", 
             RowBox[{"40.", "*", 
              RowBox[{"2.5", "^", ".85"}]}]}], ")"}], "*", 
           RowBox[{"Exp", "[", 
            RowBox[{
             RowBox[{"-", ".5"}], "*", 
             RowBox[{"(", 
              RowBox[{"t", "-", "6."}], ")"}]}], "]"}]}]}], "]"}]}], "]"}]}], 
     ";", 
     RowBox[{
      RowBox[{"b3", "[", "t_", "]"}], ":=", 
      RowBox[{"Log", "[", 
       RowBox[{"1.", "+", 
        RowBox[{"Exp", "[", 
         RowBox[{"k", "*", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"T3", "[", "t", "]"}], "-", "Tc"}], ")"}]}], "]"}]}], 
       "]"}]}], ";", 
     RowBox[{
      RowBox[{"sfn3", "[", "t_", "]"}], "=", 
      RowBox[{
       RowBox[{"NDSolve", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{
            RowBox[{
             RowBox[{"y3", "'"}], "[", "t", "]"}], "\[Equal]", 
            RowBox[{
             RowBox[{"-", 
              RowBox[{"b3", "[", "t", "]"}]}], "*", "n", "*", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{"-", 
                 RowBox[{"y3", "[", "t", "]"}]}], "/", 
                RowBox[{"b3", "[", "t", "]"}]}], ")"}], "^", 
              RowBox[{"(", 
               RowBox[{
                RowBox[{"(", 
                 RowBox[{"n", "-", "1."}], ")"}], "/", "n"}], ")"}]}]}]}], 
           ",", 
           RowBox[{
            RowBox[{"y3", "[", "0.", "]"}], "\[Equal]", 
            RowBox[{"-", ".00001"}]}]}], "}"}], ",", 
         RowBox[{"y3", "[", "t", "]"}], ",", 
         RowBox[{"{", 
          RowBox[{"t", ",", "0.001", ",", 
           RowBox[{"tCycle3", "+", "1."}]}], "}"}]}], "]"}], "[", 
       RowBox[{"[", 
        RowBox[{"1", ",", "1", ",", "2"}], "]"}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"sfn3", "[", "tFinal", "]"}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{3.393773259076874*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
Compute and save the three final survival function values at each final time \
(tCycle) using the original values for the parameters k, n and Tc.  targetLS1 \
is the target survival value using the original values of all three parmeters \
k, n and Tc.\
\>", "Text",
 CellChangeTimes->{
  3.393240631786151*^9, {3.3932406926494513`*^9, 3.393240700414157*^9}, {
   3.393252830733748*^9, 3.393252835413404*^9}, 3.393328032193994*^9, {
   3.3933280834581413`*^9, 3.3933280838747587`*^9}, {3.393328301901437*^9, 
   3.393328325090005*^9}, {3.39332899578542*^9, 3.393329004051153*^9}},
 TextAlignment->Left,
 FontSize->14],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"lstCycle1", "=", 
   RowBox[{"LS1", "[", 
    RowBox[{"k", ",", "n", ",", "Tc", ",", "tCycle1"}], "]"}]}], ";", 
  RowBox[{"targetLS1", "=", "lstCycle1"}], ";", 
  RowBox[{"Print", "[", 
   RowBox[{"\"\<LS1[k,n,Tc,tCycle1]=\>\"", ",", "lstCycle1"}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.393328111310095*^9, 3.393328131458684*^9}, {
  3.39332827647143*^9, 3.393328282523467*^9}, {3.393328342399544*^9, 
  3.393328364344249*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"LS1[k,n,Tc,tCycle1]=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "3.428031736413475`"}]}],
  SequenceForm["LS1[k,n,Tc,tCycle1]=", -3.428031736413475],
  Editable->False]], "Print",
 CellChangeTimes->{3.393329248915043*^9, 3.3933296629005003`*^9, 
  3.393773684866783*^9, 3.393774829648288*^9, 3.39377515790489*^9, 
  3.393776198512752*^9, 3.393843291352399*^9, 3.393844522045969*^9, 
  3.394188019894277*^9, 3.39418834550069*^9, 3.394188780459243*^9, 
  3.394189225785534*^9, 3.3941896766435623`*^9, 3.3941898282395773`*^9, 
  3.394189907471603*^9, 3.39419018660972*^9, 3.3941948202041883`*^9, 
  3.3949610727347717`*^9, 3.3949624874282093`*^9, 3.394964724632956*^9, 
  3.3949669333968773`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"lstCycle2", "=", 
   RowBox[{"LS2", "[", 
    RowBox[{"k", ",", "n", ",", "Tc", ",", "tCycle2"}], "]"}]}], ";", 
  RowBox[{"Print", "[", 
   RowBox[{"\"\<LS2[k,n,Tc,tCycle2]=\>\"", ",", "lstCycle2"}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.393328167966692*^9, 3.393328183737896*^9}, {
  3.3933283727577753`*^9, 3.393328393870144*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"LS2[k,n,Tc,tCycle2]=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "7.117937915500087`"}]}],
  SequenceForm["LS2[k,n,Tc,tCycle2]=", -7.117937915500087],
  Editable->False]], "Print",
 CellChangeTimes->{3.3933292490380363`*^9, 3.393329663022585*^9, 
  3.393773684995525*^9, 3.393774829794098*^9, 3.393775158030463*^9, 
  3.3937761986192427`*^9, 3.393843291489208*^9, 3.393844522147708*^9, 
  3.3941880200108137`*^9, 3.394188345599552*^9, 3.394188780584734*^9, 
  3.394189225929885*^9, 3.3941896767783213`*^9, 3.394189828342917*^9, 
  3.394189907567644*^9, 3.394190186710937*^9, 3.3941948202920313`*^9, 
  3.394961072831262*^9, 3.394962487545512*^9, 3.394964724716858*^9, 
  3.394966933509622*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"lstCycle3", "=", 
   RowBox[{"LS3", "[", 
    RowBox[{"k", ",", "n", ",", "Tc", ",", "tCycle3"}], "]"}]}], ";", 
  RowBox[{"Print", "[", 
   RowBox[{"\"\<LS3[k,n,Tc,tCycle3]=\>\"", ",", "lstCycle3"}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.393328191203388*^9, 3.3933281978368673`*^9}, {
  3.393328399368285*^9, 3.393328404935911*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"LS3[k,n,Tc,tCycle3]=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "10.852269845844909`"}]}],
  SequenceForm["LS3[k,n,Tc,tCycle3]=", -10.852269845844909`],
  Editable->False]], "Print",
 CellChangeTimes->{3.393329249399015*^9, 3.39332966313656*^9, 
  3.3937736851272917`*^9, 3.393774829919445*^9, 3.393775158160285*^9, 
  3.393776198843849*^9, 3.393843291613603*^9, 3.393844522233078*^9, 
  3.39418802011979*^9, 3.394188345693172*^9, 3.394188780710678*^9, 
  3.39418922605256*^9, 3.394189676914523*^9, 3.394189828425954*^9, 
  3.394189907660692*^9, 3.3941901868486347`*^9, 3.394194820404652*^9, 
  3.394961072918508*^9, 3.394962487650503*^9, 3.394964724787868*^9, 
  3.394966933591548*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
Iterate values of the parameter n starting with n=nSt up to n=nEn with an \
increment of nInc.  Stop iterating when the absolute value of the difference \
between the target and trial values of LS1 is less than tol.  NOTICE THAT THE \
PROGRAM FINDS THE ORIGINAL GENERATION PARAMETER VALUES k=0.065, n=0.8 and \
Tc=100.\
\>", "Text",
 CellChangeTimes->{{3.3937759660918694`*^9, 3.393775972297305*^9}, 
   3.394189426455537*^9},
 TextAlignment->Left,
 FontSize->14],

Cell[BoxData[
 RowBox[{
  RowBox[{"nSt", "=", ".7"}], ";", 
  RowBox[{"nEn", "=", ".9"}], ";", 
  RowBox[{"nInc", "=", "0.02"}], ";"}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"tol", "=", "0.01"}], ";"}]], "Input",
 CellChangeTimes->{3.3932407866131372`*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Do", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"Clear", "[", 
     RowBox[{"k", ",", "Tc"}], "]"}], ";", 
    RowBox[{"frResult", "=", 
     RowBox[{"FindRoot", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{
          RowBox[{"LS3", "[", 
           RowBox[{"k", ",", "n", ",", "Tc", ",", "tCycle3"}], "]"}], "==", 
          "lstCycle3"}], ",", 
         RowBox[{
          RowBox[{"LS2", "[", 
           RowBox[{"k", ",", "n", ",", "Tc", ",", "tCycle2"}], "]"}], 
          "\[Equal]", "lstCycle2"}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{"k", ",", ".05", ",", ".06"}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{"Tc", ",", "80.", ",", "90."}], "}"}]}], "}"}], ",", 
       RowBox[{"MaxIterations", "\[Rule]", "2000"}]}], "]"}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"k", "=", 
     RowBox[{"frResult", "[", 
      RowBox[{"[", 
       RowBox[{"1", ",", "2"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"Tc", "=", 
     RowBox[{"frResult", "[", 
      RowBox[{"[", 
       RowBox[{"2", ",", "2"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"trialLS1", "=", 
     RowBox[{"LS1", "[", 
      RowBox[{"k", ",", "n", ",", "Tc", ",", "tCycle1"}], "]"}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"Print", "[", 
     RowBox[{
     "\"\<trialLS1= \>\"", ",", "trialLS1", ",", "\"\<   k= \>\"", ",", "k", 
      ",", "\"\<   n= \>\"", ",", "n", ",", "\"\<   Tc= \>\"", ",", "Tc"}], 
     "]"}], ";", "\[IndentingNewLine]", 
    RowBox[{"If", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"Abs", "[", 
        RowBox[{"targetLS1", "-", "trialLS1"}], "]"}], "<", "tol"}], ",", 
      RowBox[{"Break", "[", "]"}]}], "]"}]}], "\[IndentingNewLine]", ",", 
   RowBox[{"{", 
    RowBox[{"n", ",", "nSt", ",", "nEn", ",", "nInc"}], "}"}]}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.3932414854276457`*^9, 3.393241485667058*^9}, {
  3.393242748224288*^9, 3.393242758488509*^9}, {3.3932428181532717`*^9, 
  3.3932428321605043`*^9}, {3.394189774106435*^9, 3.3941897982279873`*^9}, {
  3.394189875842482*^9, 3.394189894037389*^9}, {3.394189997851327*^9, 
  3.394189998559903*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"trialLS1= \"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "3.280171973339469`"}], "\[InvisibleSpace]", "\<\"   k= \"\>",
    "\[InvisibleSpace]", "0.07587235182983353`", 
   "\[InvisibleSpace]", "\<\"   n= \"\>", "\[InvisibleSpace]", "0.7`", 
   "\[InvisibleSpace]", "\<\"   Tc= \"\>", "\[InvisibleSpace]", 
   "98.43870471409251`"}],
  SequenceForm[
  "trialLS1= ", -3.280171973339469, "   k= ", 0.07587235182983353, "   n= ", 
   0.7, "   Tc= ", 98.43870471409251],
  Editable->False]], "Print",
 CellChangeTimes->{3.393167885397414*^9, 3.393241317714065*^9, 
  3.393241438860729*^9, 3.393241682804503*^9, 3.393242913055911*^9, 
  3.3932501122223186`*^9, 3.3932507006513367`*^9, 3.393253094982565*^9, 
  3.393329250381813*^9, 3.3933296640171432`*^9, 3.393773686100819*^9, 
  3.393774839901455*^9, 3.3937751620406446`*^9, 3.3937761997900953`*^9, 
  3.393843292789316*^9, 3.39384452308354*^9, 3.394188021294117*^9, 
  3.394188346487739*^9, 3.394188781827223*^9, 3.394189227090362*^9, 
  3.3941896781017437`*^9, 3.394189829270961*^9, 3.394189908500795*^9, 
  3.3941901879394693`*^9, 3.3941948216981153`*^9, 3.394961073771452*^9, 
  3.394962488803829*^9, 3.39496472559095*^9, 3.39496693455931*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"trialLS1= \"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "3.3100547634345117`"}], 
   "\[InvisibleSpace]", "\<\"   k= \"\>", "\[InvisibleSpace]", 
   "0.07349653539604928`", "\[InvisibleSpace]", "\<\"   n= \"\>", 
   "\[InvisibleSpace]", "0.72`", "\[InvisibleSpace]", "\<\"   Tc= \"\>", 
   "\[InvisibleSpace]", "98.66670818824683`"}],
  SequenceForm[
  "trialLS1= ", -3.3100547634345117`, "   k= ", 0.07349653539604928, "   n= ",
    0.72, "   Tc= ", 98.66670818824683],
  Editable->False]], "Print",
 CellChangeTimes->{3.393167885397414*^9, 3.393241317714065*^9, 
  3.393241438860729*^9, 3.393241682804503*^9, 3.393242913055911*^9, 
  3.3932501122223186`*^9, 3.3932507006513367`*^9, 3.393253094982565*^9, 
  3.393329250381813*^9, 3.3933296640171432`*^9, 3.393773686100819*^9, 
  3.393774839901455*^9, 3.3937751620406446`*^9, 3.3937761997900953`*^9, 
  3.393843292789316*^9, 3.39384452308354*^9, 3.394188021294117*^9, 
  3.394188346487739*^9, 3.394188781827223*^9, 3.394189227090362*^9, 
  3.3941896781017437`*^9, 3.394189829270961*^9, 3.394189908500795*^9, 
  3.3941901879394693`*^9, 3.3941948216981153`*^9, 3.394961073771452*^9, 
  3.394962488803829*^9, 3.39496472559095*^9, 3.39496693536919*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"trialLS1= \"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "3.339841727705481`"}], "\[InvisibleSpace]", "\<\"   k= \"\>",
    "\[InvisibleSpace]", "0.07122527095895613`", 
   "\[InvisibleSpace]", "\<\"   n= \"\>", "\[InvisibleSpace]", "0.74`", 
   "\[InvisibleSpace]", "\<\"   Tc= \"\>", "\[InvisibleSpace]", 
   "98.93562795165332`"}],
  SequenceForm[
  "trialLS1= ", -3.339841727705481, "   k= ", 0.07122527095895613, "   n= ", 
   0.74, "   Tc= ", 98.93562795165332],
  Editable->False]], "Print",
 CellChangeTimes->{3.393167885397414*^9, 3.393241317714065*^9, 
  3.393241438860729*^9, 3.393241682804503*^9, 3.393242913055911*^9, 
  3.3932501122223186`*^9, 3.3932507006513367`*^9, 3.393253094982565*^9, 
  3.393329250381813*^9, 3.3933296640171432`*^9, 3.393773686100819*^9, 
  3.393774839901455*^9, 3.3937751620406446`*^9, 3.3937761997900953`*^9, 
  3.393843292789316*^9, 3.39384452308354*^9, 3.394188021294117*^9, 
  3.394188346487739*^9, 3.394188781827223*^9, 3.394189227090362*^9, 
  3.3941896781017437`*^9, 3.394189829270961*^9, 3.394189908500795*^9, 
  3.3941901879394693`*^9, 3.3941948216981153`*^9, 3.394961073771452*^9, 
  3.394962488803829*^9, 3.39496472559095*^9, 3.394966936009486*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"trialLS1= \"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "3.3694716472441204`"}], 
   "\[InvisibleSpace]", "\<\"   k= \"\>", "\[InvisibleSpace]", 
   "0.0690546665796205`", "\[InvisibleSpace]", "\<\"   n= \"\>", 
   "\[InvisibleSpace]", "0.76`", "\[InvisibleSpace]", "\<\"   Tc= \"\>", 
   "\[InvisibleSpace]", "99.24674606204337`"}],
  SequenceForm[
  "trialLS1= ", -3.3694716472441204`, "   k= ", 0.0690546665796205, "   n= ", 
   0.76, "   Tc= ", 99.24674606204337],
  Editable->False]], "Print",
 CellChangeTimes->{3.393167885397414*^9, 3.393241317714065*^9, 
  3.393241438860729*^9, 3.393241682804503*^9, 3.393242913055911*^9, 
  3.3932501122223186`*^9, 3.3932507006513367`*^9, 3.393253094982565*^9, 
  3.393329250381813*^9, 3.3933296640171432`*^9, 3.393773686100819*^9, 
  3.393774839901455*^9, 3.3937751620406446`*^9, 3.3937761997900953`*^9, 
  3.393843292789316*^9, 3.39384452308354*^9, 3.394188021294117*^9, 
  3.394188346487739*^9, 3.394188781827223*^9, 3.394189227090362*^9, 
  3.3941896781017437`*^9, 3.394189829270961*^9, 3.394189908500795*^9, 
  3.3941901879394693`*^9, 3.3941948216981153`*^9, 3.394961073771452*^9, 
  3.394962488803829*^9, 3.39496472559095*^9, 3.39496693697405*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"trialLS1= \"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "3.3988876180545873`"}], 
   "\[InvisibleSpace]", "\<\"   k= \"\>", "\[InvisibleSpace]", 
   "0.06698083410849545`", "\[InvisibleSpace]", "\<\"   n= \"\>", 
   "\[InvisibleSpace]", "0.7799999999999999`", 
   "\[InvisibleSpace]", "\<\"   Tc= \"\>", "\[InvisibleSpace]", 
   "99.60119271937853`"}],
  SequenceForm[
  "trialLS1= ", -3.3988876180545873`, "   k= ", 0.06698083410849545, "   n= ",
    0.7799999999999999, "   Tc= ", 99.60119271937853],
  Editable->False]], "Print",
 CellChangeTimes->{3.393167885397414*^9, 3.393241317714065*^9, 
  3.393241438860729*^9, 3.393241682804503*^9, 3.393242913055911*^9, 
  3.3932501122223186`*^9, 3.3932507006513367`*^9, 3.393253094982565*^9, 
  3.393329250381813*^9, 3.3933296640171432`*^9, 3.393773686100819*^9, 
  3.393774839901455*^9, 3.3937751620406446`*^9, 3.3937761997900953`*^9, 
  3.393843292789316*^9, 3.39384452308354*^9, 3.394188021294117*^9, 
  3.394188346487739*^9, 3.394188781827223*^9, 3.394189227090362*^9, 
  3.3941896781017437`*^9, 3.394189829270961*^9, 3.394189908500795*^9, 
  3.3941901879394693`*^9, 3.3941948216981153`*^9, 3.394961073771452*^9, 
  3.394962488803829*^9, 3.39496472559095*^9, 3.394966938397306*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"trialLS1= \"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "3.4280321466760553`"}], 
   "\[InvisibleSpace]", "\<\"   k= \"\>", "\[InvisibleSpace]", 
   "0.06499999402085757`", "\[InvisibleSpace]", "\<\"   n= \"\>", 
   "\[InvisibleSpace]", "0.7999999999999999`", 
   "\[InvisibleSpace]", "\<\"   Tc= \"\>", "\[InvisibleSpace]", 
   "99.99999648385807`"}],
  SequenceForm[
  "trialLS1= ", -3.4280321466760553`, "   k= ", 0.06499999402085757, "   n= ",
    0.7999999999999999, "   Tc= ", 99.99999648385807],
  Editable->False]], "Print",
 CellChangeTimes->{3.393167885397414*^9, 3.393241317714065*^9, 
  3.393241438860729*^9, 3.393241682804503*^9, 3.393242913055911*^9, 
  3.3932501122223186`*^9, 3.3932507006513367`*^9, 3.393253094982565*^9, 
  3.393329250381813*^9, 3.3933296640171432`*^9, 3.393773686100819*^9, 
  3.393774839901455*^9, 3.3937751620406446`*^9, 3.3937761997900953`*^9, 
  3.393843292789316*^9, 3.39384452308354*^9, 3.394188021294117*^9, 
  3.394188346487739*^9, 3.394188781827223*^9, 3.394189227090362*^9, 
  3.3941896781017437`*^9, 3.394189829270961*^9, 3.394189908500795*^9, 
  3.3941901879394693`*^9, 3.3941948216981153`*^9, 3.394961073771452*^9, 
  3.394962488803829*^9, 3.39496472559095*^9, 3.394966938886125*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Plot the three temperature profile functions separately.", "Text",
 CellChangeTimes->{{3.3931678109056807`*^9, 3.393167827311077*^9}, 
   3.3932409362764807`*^9, 3.393253017035926*^9, 3.39332715569305*^9, {
   3.3937732803530817`*^9, 3.3937732832844467`*^9}},
 TextAlignment->Left,
 FontSize->14],

Cell[BoxData[
 RowBox[{
  RowBox[{"tempCurveT1Plot", "=", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"T1", "[", "t", "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"t", ",", "0.", ",", "tCycle1"}], "}"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0.", ",", "10."}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"20.", ",", "160."}], "}"}]}], "}"}]}], ",", 
     RowBox[{"AspectRatio", "\[Rule]", "1"}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"Thickness", "[", "0.01", "]"}], ",", "LightGray", ",", 
         RowBox[{"Dashing", "[", "None", "]"}]}], "}"}], "}"}]}], ",", 
     RowBox[{"PlotLabel", "\[Rule]", "\"\<Temperature Profiles\>\""}], ",", 
     RowBox[{"Frame", "\[Rule]", "True"}], ",", 
     RowBox[{"RotateLabel", "\[Rule]", "True"}], ",", 
     RowBox[{"BaseStyle", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"FontFamily", "\[Rule]", "\"\<Helvetica\>\""}], ",", 
        RowBox[{"FontWeight", "\[Rule]", "\"\<Bold\>\""}], ",", 
        RowBox[{"FontSize", "\[Rule]", "14"}]}], "}"}]}], ",", 
     RowBox[{"FrameLabel", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"\"\<\>\"", ",", "\"\<Temperature (\[Degree]C)\>\""}], "}"}]}],
      ",", 
     RowBox[{"FrameTicks", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0", ",", "2", ",", "4", ",", "6", ",", "8", ",", "10"}], 
         "}"}], ",", 
        RowBox[{"{", 
         RowBox[{
         "20", ",", "40", ",", "60", ",", "80", ",", "100", ",", "120", ",", 
          "140", ",", "160"}], "}"}], ",", 
        RowBox[{"{", "}"}], ",", 
        RowBox[{"{", "}"}]}], "}"}]}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{
  3.393250047884212*^9, {3.393250439443324*^9, 3.393250463233591*^9}, {
   3.393250537805275*^9, 3.39325054746657*^9}, {3.393250635783272*^9, 
   3.393250648964038*^9}, {3.393329130805213*^9, 3.393329230622287*^9}, {
   3.39332926923064*^9, 3.393329297599441*^9}, 3.393329361557909*^9, 
   3.393329392420142*^9, {3.393329422559197*^9, 3.3933294510955667`*^9}, 
   3.3933297209700003`*^9, {3.393773322890156*^9, 3.393773352567934*^9}, {
   3.393773391502338*^9, 3.393773417731648*^9}, 3.393775082855269*^9, {
   3.39418784643503*^9, 3.394187883042354*^9}, {3.3941881522898397`*^9, 
   3.394188158929481*^9}, {3.394188628965201*^9, 3.394188633260126*^9}, {
   3.394188702198111*^9, 3.394188709717126*^9}, 3.3941887501024113`*^9, 
   3.394188837297305*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"tempCurveT2Plot", "=", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"T2", "[", "t", "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"t", ",", "0.", ",", "tCycle2"}], "}"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0.", ",", "10."}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"20.", ",", "160."}], "}"}]}], "}"}]}], ",", 
     RowBox[{"AspectRatio", "\[Rule]", "1"}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"Thickness", "[", "0.01", "]"}], ",", "Gray", ",", 
         RowBox[{"Dashing", "[", "None", "]"}]}], "}"}], "}"}]}], ",", 
     RowBox[{"PlotLabel", "\[Rule]", "\"\<Temperature Profiles\>\""}], ",", 
     RowBox[{"Frame", "\[Rule]", "True"}], ",", 
     RowBox[{"RotateLabel", "\[Rule]", "True"}], ",", 
     RowBox[{"BaseStyle", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"FontFamily", "\[Rule]", "\"\<Helvetica\>\""}], ",", 
        RowBox[{"FontWeight", "\[Rule]", "\"\<Bold\>\""}], ",", 
        RowBox[{"FontSize", "\[Rule]", "14"}]}], "}"}]}], ",", 
     RowBox[{"FrameLabel", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
       "\"\<Time (arbitrary units)\>\"", ",", 
        "\"\<Temperature (\[Degree]C)\>\""}], "}"}]}], ",", 
     RowBox[{"FrameTicks", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0", ",", "2", ",", "4", ",", "6", ",", "8", ",", "10"}], 
         "}"}], ",", 
        RowBox[{"{", 
         RowBox[{
         "20", ",", "40", ",", "60", ",", "80", ",", "100", ",", "120", ",", 
          "140", ",", "160"}], "}"}], ",", 
        RowBox[{"{", "}"}], ",", 
        RowBox[{"{", "}"}]}], "}"}]}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{
  3.393250047884212*^9, {3.393250439443324*^9, 3.393250463233591*^9}, {
   3.393250537805275*^9, 3.39325054746657*^9}, {3.393250635783272*^9, 
   3.393250648964038*^9}, {3.393329130805213*^9, 3.393329230622287*^9}, {
   3.39332926923064*^9, 3.393329297599441*^9}, 3.393329361557909*^9, 
   3.393329392420142*^9, {3.393329422559197*^9, 3.3933294510955667`*^9}, 
   3.3933297209700003`*^9, {3.393773328128667*^9, 3.393773331061606*^9}, {
   3.39377345772255*^9, 3.393773488524987*^9}, 3.393775086623687*^9, {
   3.394187890163568*^9, 3.394187902316348*^9}, 3.39418863886768*^9, {
   3.394188692519002*^9, 3.394188694324132*^9}, 3.394188753358597*^9, 
   3.394188835370954*^9, 3.3941892011974497`*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"tempCurveT3Plot", "=", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"T3", "[", "t", "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"t", ",", "0.", ",", "tCycle3"}], "}"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0.", ",", "10."}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"20.", ",", "160."}], "}"}]}], "}"}]}], ",", 
     RowBox[{"AspectRatio", "\[Rule]", "1"}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"Thickness", "[", "0.01", "]"}], ",", "Black", ",", 
         RowBox[{"Dashing", "[", "None", "]"}]}], "}"}], "}"}]}], ",", 
     RowBox[{"PlotLabel", "\[Rule]", "\"\<Temperature Profiles\>\""}], ",", 
     RowBox[{"Frame", "\[Rule]", "True"}], ",", 
     RowBox[{"RotateLabel", "\[Rule]", "True"}], ",", 
     RowBox[{"BaseStyle", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"FontFamily", "\[Rule]", "\"\<Helvetica\>\""}], ",", 
        RowBox[{"FontWeight", "\[Rule]", "\"\<Bold\>\""}], ",", 
        RowBox[{"FontSize", "\[Rule]", "14"}]}], "}"}]}], ",", 
     RowBox[{"FrameLabel", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
       "\"\<Time (arbitrary units)\>\"", ",", 
        "\"\<Temperature (\[Degree]C)\>\""}], "}"}]}], ",", 
     RowBox[{"FrameTicks", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0", ",", "2", ",", "4", ",", "6", ",", "8", ",", "10"}], 
         "}"}], ",", 
        RowBox[{"{", 
         RowBox[{
         "20", ",", "40", ",", "60", ",", "80", ",", "100", ",", "120", ",", 
          "140", ",", "160"}], "}"}], ",", 
        RowBox[{"{", "}"}], ",", 
        RowBox[{"{", "}"}]}], "}"}]}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{
  3.393250047884212*^9, {3.393250439443324*^9, 3.393250463233591*^9}, {
   3.393250537805275*^9, 3.39325054746657*^9}, {3.393250635783272*^9, 
   3.393250648964038*^9}, {3.393329130805213*^9, 3.393329230622287*^9}, {
   3.39332926923064*^9, 3.393329297599441*^9}, 3.393329361557909*^9, 
   3.393329392420142*^9, {3.393329422559197*^9, 3.3933294510955667`*^9}, 
   3.3933297209700003`*^9, {3.393773340625329*^9, 3.393773341397279*^9}, {
   3.393773496098831*^9, 3.393773516869606*^9}, 3.393775089063445*^9, 
   3.3941876269157963`*^9, {3.394187906299838*^9, 3.394187919891152*^9}, 
   3.394188163755245*^9, {3.394188644159461*^9, 3.3941886740831137`*^9}, {
   3.3941887554365*^9, 3.3941887683411837`*^9}, 3.394189194997942*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
ListPlot the temperatures at at times tCycle1, tCycle2 and tCycle3 for the \
three temperature profile functions.\
\>", "Text",
 CellChangeTimes->{{3.393168351042939*^9, 3.393168364232419*^9}, 
   3.3932537841910152`*^9, {3.393322372504961*^9, 3.3933224384617233`*^9}, 
   3.393322770438449*^9, 3.3933233004745197`*^9, {3.3933239438840647`*^9, 
   3.393323944731514*^9}, 3.39332503330643*^9, {3.393329514928504*^9, 
   3.3933295216675053`*^9}, 3.393844512283464*^9},
 TextAlignment->Left,
 FontSize->14],

Cell[BoxData[
 RowBox[{
  RowBox[{"tempPtsPlot", "=", 
   RowBox[{"ListPlot", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"tCycle1", ",", "tCycle1T1"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"tCycle2", ",", "tCycle2T2"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"tCycle3", ",", "tCycle3T3"}], "}"}]}], "}"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0.", ",", "10."}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"20.", ",", "160."}], "}"}]}], "}"}]}], ",", 
     RowBox[{"AspectRatio", "\[Rule]", "1"}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"Directive", "[", 
       RowBox[{
        RowBox[{"PointSize", "[", "Large", "]"}], ",", "Black"}], "]"}]}], 
     ",", 
     RowBox[{
     "PlotLabel", "\[Rule]", 
      "\"\<Temperatures at tCycle1, tCycle2 & tCycle3\>\""}], ",", 
     RowBox[{"Frame", "\[Rule]", "True"}], ",", 
     RowBox[{"RotateLabel", "\[Rule]", "True"}], ",", 
     RowBox[{"BaseStyle", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"FontFamily", "\[Rule]", "\"\<Helvetica\>\""}], ",", 
        RowBox[{"FontWeight", "\[Rule]", "\"\<Bold\>\""}], ",", 
        RowBox[{"FontSize", "\[Rule]", "14"}]}], "}"}]}], ",", 
     RowBox[{"FrameLabel", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
       "\"\<Time (arbitrary units)\>\"", ",", 
        "\"\<Temperature (\[Degree]C)\>\""}], "}"}]}], ",", 
     RowBox[{"FrameTicks", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0", ",", "2", ",", "4", ",", "6", ",", "8", ",", "10"}], 
         "}"}], ",", 
        RowBox[{"{", 
         RowBox[{
         "20", ",", "40", ",", "60", ",", "80", ",", "100", ",", "120", ",", 
          "140", ",", "160"}], "}"}], ",", 
        RowBox[{"{", "}"}], ",", 
        RowBox[{"{", "}"}]}], "}"}]}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{
  3.393251050701592*^9, {3.393320294224659*^9, 3.393320333211245*^9}, {
   3.393320571080579*^9, 3.393320590759631*^9}, {3.3933208225830812`*^9, 
   3.393320886103772*^9}, {3.393320919192649*^9, 3.3933209459323063`*^9}, {
   3.393320993422297*^9, 3.393321031966597*^9}, {3.393321480456262*^9, 
   3.3933215239096737`*^9}, {3.393321595255995*^9, 3.393321647106287*^9}, {
   3.393321707336195*^9, 3.39332171063087*^9}, {3.3933217576497927`*^9, 
   3.393321775498185*^9}, 3.3933218321220922`*^9, {3.393321874676762*^9, 
   3.393322013321027*^9}, {3.393322462851087*^9, 3.393322463364461*^9}, 
   3.3933227992668552`*^9, {3.393322838191132*^9, 3.393322875539483*^9}, {
   3.393322989348707*^9, 3.3933230303468647`*^9}, {3.393323166028068*^9, 
   3.3933231919176607`*^9}, 3.393323254576283*^9, {3.393323328414671*^9, 
   3.393323354375997*^9}, {3.393323489175621*^9, 3.39332349837894*^9}, {
   3.393323544312173*^9, 3.393323557908409*^9}, {3.3933237239305573`*^9, 
   3.393323751116066*^9}, 3.3933240451270103`*^9, {3.3933240967479343`*^9, 
   3.393324126899761*^9}, {3.393324719907443*^9, 3.39332473324393*^9}, 
   3.3933254311157*^9, {3.3933257483089733`*^9, 3.393325755298956*^9}, {
   3.393325812335999*^9, 3.393325843018764*^9}, {3.393326930078993*^9, 
   3.3933269324271517`*^9}, {3.3933270478841248`*^9, 3.393327054210685*^9}, {
   3.3933295264436274`*^9, 3.393329591902258*^9}, {3.3937734219003696`*^9, 
   3.393773422499653*^9}, 3.393775079519662*^9, {3.394187834728739*^9, 
   3.3941878618209333`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
Show the three temperature profile curves and the three final temperature \
points together on the same plot.\
\>", "Text",
 CellChangeTimes->{{3.393168351042939*^9, 3.393168364232419*^9}, 
   3.3932537841910152`*^9, 3.393322372504961*^9, {3.393323367039167*^9, 
   3.393323441740209*^9}, {3.3933296968342333`*^9, 3.39332970570168*^9}},
 TextAlignment->Left,
 FontSize->14],

Cell[BoxData[
 RowBox[{
  RowBox[{"tempProfilePlot", "=", 
   RowBox[{"Show", "[", 
    RowBox[{
    "tempCurveT1Plot", ",", "tempCurveT2Plot", ",", "tempCurveT3Plot", ",", 
     "tempPtsPlot", ",", 
     RowBox[{"ImageSize", "\[Rule]", 
      RowBox[{"Scaled", "[", "0.585", "]"}]}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.393323062018441*^9, 3.393323083173884*^9}, {
   3.393323272934407*^9, 3.393323282685367*^9}, {3.393323520755739*^9, 
   3.393323526992766*^9}, 3.393323796966511*^9, 3.393324168607872*^9, {
   3.3937735234755507`*^9, 3.393773537003829*^9}, 3.394962239746101*^9, 
   3.394964710025722*^9, {3.394964775659923*^9, 3.394964775966889*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Plot the three survival curves separately.", "Text",
 CellChangeTimes->{{3.3931678109056807`*^9, 3.393167827311077*^9}, {
  3.393167859756092*^9, 3.3931678691693497`*^9}, {3.3932530264807863`*^9, 
  3.3932530272597923`*^9}, {3.393773825536284*^9, 3.393773846728099*^9}},
 TextAlignment->Left,
 FontSize->14],

Cell[BoxData[
 RowBox[{
  RowBox[{"lstCurveLS1Plot", "=", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"LS1", "[", 
      RowBox[{".065", ",", ".8", ",", "100.", ",", "t"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"t", ",", "0.001", ",", "tCycle1"}], "}"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0.", ",", "10."}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"-", "12."}], ",", "0.5"}], "}"}]}], "}"}]}], ",", 
     RowBox[{"AspectRatio", "\[Rule]", "1"}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"{", " ", 
        RowBox[{
         RowBox[{"Thickness", "[", "0.01", "]"}], ",", "LightGray", ",", 
         RowBox[{"Dashing", "[", "None", "]"}]}], "}"}], "}"}]}], ",", 
     RowBox[{"PlotLabel", "\[Rule]", "\"\<Survival Curves\\nn<1\>\""}], ",", 
     RowBox[{"Frame", "\[Rule]", "True"}], ",", 
     RowBox[{"RotateLabel", "\[Rule]", "True"}], ",", 
     RowBox[{"BaseStyle", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"FontFamily", "\[Rule]", "\"\<Helvetica\>\""}], ",", 
        RowBox[{"FontWeight", "\[Rule]", "\"\<Bold\>\""}], ",", 
        RowBox[{"FontSize", "\[Rule]", "14"}]}], "}"}]}], ",", 
     RowBox[{"FrameLabel", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
       "\"\<Time (arbitrary units)\>\"", ",", 
        "\"\<\!\(\*SubscriptBox[\(Log\), \(10\)]\)S(t)\>\""}], "}"}]}], ",", 
     RowBox[{"FrameTicks", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0", ",", "2", ",", "4", ",", "6", ",", "8", ",", "10"}], 
         "}"}], ",", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"-", "12"}], ",", 
          RowBox[{"-", "10"}], ",", 
          RowBox[{"-", "8"}], ",", 
          RowBox[{"-", "6"}], ",", 
          RowBox[{"-", "4"}], ",", 
          RowBox[{"-", "2"}], ",", "0"}], "}"}], ",", 
        RowBox[{"{", "}"}], ",", 
        RowBox[{"{", "}"}]}], "}"}]}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.393090260110411*^9, 3.393090299577162*^9}, 
   3.39309034875788*^9, 3.3932363616796255`*^9, 3.393239825137187*^9, 
   3.39323991665901*^9, {3.39324338212543*^9, 3.3932434514121437`*^9}, {
   3.3932516147470293`*^9, 3.393251664445907*^9}, {3.3933440408495913`*^9, 
   3.3933440881418743`*^9}, {3.3933441635810003`*^9, 3.39334422746439*^9}, {
   3.3933443397075567`*^9, 3.393344340088715*^9}, {3.393771028920472*^9, 
   3.393771096174671*^9}, 3.3937713840630617`*^9, 3.393771612575638*^9, {
   3.393774331359338*^9, 3.393774346202752*^9}, 3.3937743923158092`*^9, {
   3.3937744274086037`*^9, 3.393774447461445*^9}, 3.393774520484353*^9, {
   3.393774587203412*^9, 3.393774604347769*^9}, {3.3937747999846907`*^9, 
   3.393774804605421*^9}, 3.3937749243186626`*^9, {3.393774961967884*^9, 
   3.393774991870867*^9}, 3.393775110389553*^9, 3.39418754066882*^9, 
   3.394187642164032*^9, {3.3941881873311863`*^9, 3.394188191059609*^9}, 
   3.394188307347417*^9, {3.3941891088362503`*^9, 3.39418914222752*^9}, 
   3.3941896528190517`*^9, 3.394189940555067*^9, {3.394190066276651*^9, 
   3.394190129877123*^9}, 3.394194874655121*^9, 3.394194955860907*^9, 
   3.394961033077262*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"lstCurveLS2Plot", "=", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"LS2", "[", 
      RowBox[{".065", ",", ".8", ",", "100.", ",", "t"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"t", ",", "0.001", ",", "tCycle2"}], "}"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0.", ",", "10."}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"-", "12."}], ",", "0.5"}], "}"}]}], "}"}]}], ",", 
     RowBox[{"AspectRatio", "\[Rule]", "1"}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"{", " ", 
        RowBox[{
         RowBox[{"Thickness", "[", "0.01", "]"}], ",", "Gray", ",", 
         RowBox[{"Dashing", "[", "None", "]"}]}], "}"}], "}"}]}], ",", 
     RowBox[{"PlotLabel", "\[Rule]", "\"\<Survival Curves\\nn<1\>\""}], ",", 
     RowBox[{"Frame", "\[Rule]", "True"}], ",", 
     RowBox[{"RotateLabel", "\[Rule]", "True"}], ",", 
     RowBox[{"BaseStyle", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"FontFamily", "\[Rule]", "\"\<Helvetica\>\""}], ",", 
        RowBox[{"FontWeight", "\[Rule]", "\"\<Bold\>\""}], ",", 
        RowBox[{"FontSize", "\[Rule]", "14"}]}], "}"}]}], ",", 
     RowBox[{"FrameLabel", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
       "\"\<Time (arbitrary units)\>\"", ",", 
        "\"\<\!\(\*SubscriptBox[\(Log\), \(10\)]\)S(t)\>\""}], "}"}]}], ",", 
     RowBox[{"FrameTicks", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0", ",", "2", ",", "4", ",", "6", ",", "8", ",", "10"}], 
         "}"}], ",", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"-", "12"}], ",", 
          RowBox[{"-", "10"}], ",", 
          RowBox[{"-", "8"}], ",", 
          RowBox[{"-", "6"}], ",", 
          RowBox[{"-", "4"}], ",", 
          RowBox[{"-", "2"}], ",", "0"}], "}"}], ",", 
        RowBox[{"{", "}"}], ",", 
        RowBox[{"{", "}"}]}], "}"}]}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.393090260110411*^9, 3.393090299577162*^9}, 
   3.39309034875788*^9, 3.3932363616796255`*^9, 3.393239825137187*^9, 
   3.39323991665901*^9, {3.39324338212543*^9, 3.3932434514121437`*^9}, {
   3.3932516147470293`*^9, 3.393251664445907*^9}, {3.3933440408495913`*^9, 
   3.3933440881418743`*^9}, {3.3933441635810003`*^9, 3.39334422746439*^9}, {
   3.3933443397075567`*^9, 3.393344340088715*^9}, {3.393771028920472*^9, 
   3.393771038960395*^9}, {3.3937711101583014`*^9, 3.393771131423704*^9}, 
   3.3937716159151897`*^9, {3.3937744305864763`*^9, 3.39377445693922*^9}, {
   3.393774523124353*^9, 3.393774544805518*^9}, {3.393774590680825*^9, 
   3.393774606248178*^9}, 3.3937749388117247`*^9, 3.393774993353723*^9, 
   3.393775121073962*^9, {3.3937751849309464`*^9, 3.393775185502006*^9}, {
   3.394187483146428*^9, 3.394187495687855*^9}, 3.3941876458626413`*^9, 
   3.3941883104606447`*^9, {3.394189111351191*^9, 3.394189152982156*^9}, 
   3.394189655301531*^9, 3.394189946975574*^9, {3.394190140517836*^9, 
   3.394190140864719*^9}, 3.3941948702186823`*^9, 3.3941949526774273`*^9, 
   3.394961036192305*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"lstCurveLS3Plot", "=", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"LS3", "[", 
      RowBox[{".065", ",", ".8", ",", "100.", ",", "t"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"t", ",", "0.001", ",", "tCycle3"}], "}"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0.", ",", "10."}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"-", "12."}], ",", "0.5"}], "}"}]}], "}"}]}], ",", 
     RowBox[{"AspectRatio", "\[Rule]", "1"}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"{", " ", 
        RowBox[{
         RowBox[{"Thickness", "[", "0.01", "]"}], ",", "Black", ",", 
         RowBox[{"Dashing", "[", "None", "]"}]}], "}"}], "}"}]}], ",", 
     RowBox[{"PlotLabel", "\[Rule]", "\"\<Survival Curves\\nn<1\>\""}], ",", 
     RowBox[{"Frame", "\[Rule]", "True"}], ",", 
     RowBox[{"RotateLabel", "\[Rule]", "True"}], ",", 
     RowBox[{"BaseStyle", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"FontFamily", "\[Rule]", "\"\<Helvetica\>\""}], ",", 
        RowBox[{"FontWeight", "\[Rule]", "\"\<Bold\>\""}], ",", 
        RowBox[{"FontSize", "\[Rule]", "14"}]}], "}"}]}], ",", 
     RowBox[{"FrameLabel", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
       "\"\<Time (arbitrary units)\>\"", ",", 
        "\"\<\!\(\*SubscriptBox[\(Log\), \(10\)]\)S(t)\>\""}], "}"}]}], ",", 
     RowBox[{"FrameTicks", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0", ",", "2", ",", "4", ",", "6", ",", "8", ",", "10"}], 
         "}"}], ",", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"-", "12"}], ",", 
          RowBox[{"-", "10"}], ",", 
          RowBox[{"-", "8"}], ",", 
          RowBox[{"-", "6"}], ",", 
          RowBox[{"-", "4"}], ",", 
          RowBox[{"-", "2"}], ",", "0"}], "}"}], ",", 
        RowBox[{"{", "}"}], ",", 
        RowBox[{"{", "}"}]}], "}"}]}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.393090260110411*^9, 3.393090299577162*^9}, 
   3.39309034875788*^9, 3.3932363616796255`*^9, 3.393239825137187*^9, 
   3.39323991665901*^9, {3.39324338212543*^9, 3.3932434514121437`*^9}, {
   3.3932516147470293`*^9, 3.393251664445907*^9}, {3.3933440408495913`*^9, 
   3.3933440881418743`*^9}, {3.3933441635810003`*^9, 3.39334422746439*^9}, {
   3.3933443397075567`*^9, 3.393344340088715*^9}, {3.393771028920472*^9, 
   3.393771029628439*^9}, {3.393771150698284*^9, 3.39377118332802*^9}, {
   3.393771372516172*^9, 3.393771377822114*^9}, 3.393771619026457*^9, {
   3.39377443412571*^9, 3.393774460342142*^9}, {3.393774525851582*^9, 
   3.393774552811377*^9}, {3.3937745933092403`*^9, 3.393774607421873*^9}, 
   3.393774808161901*^9, 3.3937749418603086`*^9, 3.393774994879775*^9, 
   3.393775124754537*^9, {3.393775188127206*^9, 3.3937751886836653`*^9}, {
   3.394187442481225*^9, 3.394187475625634*^9}, 3.394187655358603*^9, 
   3.394188195830495*^9, 3.394188319026825*^9, {3.394189113709844*^9, 
   3.394189156445352*^9}, 3.394189658300911*^9, 3.394189951335783*^9, {
   3.3941901447699623`*^9, 3.394190145151984*^9}, 3.394194806095833*^9, 
   3.394194951156131*^9, 3.3949610422574997`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
ListPlot the survival values at times tCycle1, tCycle2 and tCycle3 for the \
three survival curves.\
\>", "Text",
 CellChangeTimes->{{3.393168351042939*^9, 3.393168364232419*^9}, 
   3.3932537841910152`*^9, {3.393322372504961*^9, 3.3933224384617233`*^9}, 
   3.393322770438449*^9, 3.3933233004745197`*^9, {3.393323859054771*^9, 
   3.393323888274333*^9}, {3.393326404691359*^9, 3.393326415582448*^9}, 
   3.393326465735125*^9, {3.393338826352191*^9, 3.3933388739886723`*^9}, {
   3.393339390263722*^9, 3.393339393939906*^9}, 3.3933394758287697`*^9, 
   3.3933419167867002`*^9, {3.3933437706360283`*^9, 3.393343802397408*^9}, {
   3.393344334096384*^9, 3.393344335366448*^9}, {3.393773944570592*^9, 
   3.3937739523266*^9}},
 TextAlignment->Left,
 FontSize->14],

Cell[BoxData[
 RowBox[{
  RowBox[{"lstPtsPlot", "=", 
   RowBox[{"ListPlot", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"tCycle1", ",", "lstCycle1"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"tCycle2", ",", "lstCycle2"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"tCycle3", ",", "lstCycle3"}], "}"}]}], "}"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0.", ",", "10."}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"-", "12."}], ",", "0.5"}], "}"}]}], "}"}]}], ",", 
     RowBox[{"AspectRatio", "\[Rule]", "1"}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"Directive", "[", 
       RowBox[{
        RowBox[{"PointSize", "[", "Large", "]"}], ",", "Black"}], "]"}]}], 
     ",", 
     RowBox[{
     "PlotLabel", "\[Rule]", 
      "\"\<Survival at tCycle1, tCycle2 & tCycle3\>\""}], ",", 
     RowBox[{"Frame", "\[Rule]", "True"}], ",", 
     RowBox[{"RotateLabel", "\[Rule]", "True"}], ",", 
     RowBox[{"BaseStyle", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"FontFamily", "\[Rule]", "\"\<Helvetica\>\""}], ",", 
        RowBox[{"FontWeight", "\[Rule]", "\"\<Bold\>\""}], ",", 
        RowBox[{"FontSize", "\[Rule]", "14"}]}], "}"}]}], ",", 
     RowBox[{"FrameLabel", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
       "\"\<Time (arbitrary units)\>\"", ",", 
        "\"\<\!\(\*SubscriptBox[\(Log\), \(10\)]\)S(t)\>\""}], "}"}]}], ",", 
     RowBox[{"FrameTicks", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0", ",", "2", ",", "4", ",", "6", ",", "8", ",", "10"}], 
         "}"}], ",", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"-", "12"}], ",", 
          RowBox[{"-", "10"}], ",", 
          RowBox[{"-", "8"}], ",", 
          RowBox[{"-", "6"}], ",", 
          RowBox[{"-", "4"}], ",", 
          RowBox[{"-", "2"}], ",", "0"}], "}"}], ",", 
        RowBox[{"{", "}"}], ",", 
        RowBox[{"{", "}"}]}], "}"}]}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{
  3.393251050701592*^9, {3.393320294224659*^9, 3.393320333211245*^9}, {
   3.393320571080579*^9, 3.393320590759631*^9}, {3.3933208225830812`*^9, 
   3.393320886103772*^9}, {3.393320919192649*^9, 3.3933209459323063`*^9}, {
   3.393320993422297*^9, 3.393321031966597*^9}, {3.393321480456262*^9, 
   3.3933215239096737`*^9}, {3.393321595255995*^9, 3.393321647106287*^9}, {
   3.393321707336195*^9, 3.39332171063087*^9}, {3.3933217576497927`*^9, 
   3.393321775498185*^9}, 3.3933218321220922`*^9, {3.393321874676762*^9, 
   3.393322013321027*^9}, {3.393322462851087*^9, 3.393322463364461*^9}, 
   3.3933227992668552`*^9, {3.393322838191132*^9, 3.393322875539483*^9}, {
   3.393322989348707*^9, 3.3933230303468647`*^9}, {3.393323166028068*^9, 
   3.3933231919176607`*^9}, 3.393323254576283*^9, {3.393323328414671*^9, 
   3.393323354375997*^9}, {3.393323489175621*^9, 3.39332349837894*^9}, {
   3.393323544312173*^9, 3.393323557908409*^9}, {3.3933237239305573`*^9, 
   3.393323751116066*^9}, {3.393323894850007*^9, 3.39332392933667*^9}, {
   3.393324198458137*^9, 3.3933242203481817`*^9}, {3.393324253002262*^9, 
   3.393324299245779*^9}, {3.3933244200930223`*^9, 3.3933244387682858`*^9}, 
   3.393324554558779*^9, {3.393325728692493*^9, 3.393325735147439*^9}, {
   3.393325871387785*^9, 3.393325876316062*^9}, 3.393327824896139*^9, 
   3.393338900893297*^9, {3.393338934394993*^9, 3.393338937474901*^9}, {
   3.393338977946499*^9, 3.393338983905766*^9}, {3.393339369878024*^9, 
   3.393339400826179*^9}, 3.393339490325576*^9, 3.39333966498777*^9, {
   3.3933438154187803`*^9, 3.393343906927309*^9}, 3.39334395010356*^9, {
   3.393344286656887*^9, 3.393344328046926*^9}, {3.3933443997857103`*^9, 
   3.393344400917284*^9}, 3.393344737862132*^9, 3.393771771283959*^9, {
   3.39377399242447*^9, 3.393774095999634*^9}, {3.3937748864857903`*^9, 
   3.393774914319977*^9}, {3.3937749730203257`*^9, 3.3937749735143757`*^9}, 
   3.393775094110031*^9, 3.394187532342841*^9, 3.39418763602553*^9, 
   3.394194807556344*^9, 3.394194946696479*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
Show the three survival curves and the final survival values at times \
tCycle1, tCycle2 and tCycle3 together on the same plot.\
\>", "Text",
 CellChangeTimes->{{3.393168351042939*^9, 3.393168364232419*^9}, 
   3.3932537841910152`*^9, 3.393322372504961*^9, {3.393323367039167*^9, 
   3.393323441740209*^9}, {3.39332435574153*^9, 3.3933243770256643`*^9}, {
   3.393326446391835*^9, 3.3933264740635233`*^9}, {3.3933395940344677`*^9, 
   3.39333960323083*^9}, {3.393344235442288*^9, 3.393344262681156*^9}, {
   3.3937744804632387`*^9, 3.3937744870941687`*^9}},
 TextAlignment->Left,
 FontSize->14],

Cell[BoxData[
 RowBox[{
  RowBox[{"survivalPlot", "=", 
   RowBox[{"Show", "[", 
    RowBox[{
    "lstCurveLS1Plot", ",", "lstCurveLS2Plot", ",", "lstCurveLS3Plot", ",", 
     "lstPtsPlot", ",", 
     RowBox[{"ImageSize", "\[Rule]", 
      RowBox[{"Scaled", "[", "0.6", "]"}]}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.393323062018441*^9, 3.393323083173884*^9}, {
   3.393323272934407*^9, 3.393323282685367*^9}, {3.393323520755739*^9, 
   3.393323526992766*^9}, 3.393323796966511*^9, 3.393324168607872*^9, {
   3.393324385406657*^9, 3.393324411709361*^9}, 3.393326584402074*^9, {
   3.393344349354158*^9, 3.393344356586041*^9}, {3.393771193393984*^9, 
   3.3937712066024857`*^9}, {3.393774467386056*^9, 3.393774517386146*^9}, 
   3.3949622552267923`*^9, {3.3949647154073257`*^9, 3.394964718386483*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
Show the temperature and the survival plots together in the same graphics \
image.\
\>", "Text",
 CellChangeTimes->{{3.393168351042939*^9, 3.393168364232419*^9}, 
   3.3932537841910152`*^9, 3.393322372504961*^9, {3.393323367039167*^9, 
   3.393323441740209*^9}, {3.39332435574153*^9, 3.3933243770256643`*^9}, {
   3.393326446391835*^9, 3.3933264740635233`*^9}, {3.394961631708231*^9, 
   3.394961672094633*^9}, {3.3949626643722353`*^9, 3.3949626667600803`*^9}},
 TextAlignment->Left,
 FontSize->14],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"figure4", "=", 
  RowBox[{"GraphicsArray", "[", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", "tempProfilePlot", "}"}], ",", 
     RowBox[{"{", "survivalPlot", "}"}]}], "}"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.394961405706451*^9, 3.39496143058958*^9}, {
  3.394961466930292*^9, 3.394961475094851*^9}, {3.3949615279010687`*^9, 
  3.39496153087962*^9}, {3.3949669262847147`*^9, 3.394966928181616*^9}}],

Cell[BoxData[
 GraphicsBox[{{}, {InsetBox[
     GraphicsBox[{{{}, {}, 
        {GrayLevel[0.85], Thickness[0.01], Dashing[None], 
         LineBox[CompressedData["
1:eJxF0nk4lHsbB/DBFCKP3YQkSyVyWt4ild8tPFI8M3aziJEWy1SkRZ2Q5SQq
soxTqQ6RpU0KvZSlQnOYkhYGkTepV8TIsWR7n/cf88d93dfn+l2/+3td93W7
HvJfvpdCoUiTlfT/bh1rllB+2zsbCMqiTINnGxCXfzNG81cuEAbllpF56SiU
f9pkQZIzEJYBc6qsanSSW5Ayd3kHEDu/ah5sBnSGHzj9+9scINhZJWzbcRQ/
NVcc0f4XEMFLktXUBSiRmxFwrOcaEBF3GwrrBOgS32+f5fEmINJbGruhHqUJ
TGHT1C4gsq+pVNZkoYypMe2NkX8DcXeqOl2qHV3hJr1aH98ARF1RnUVQC8rm
s69qOvYAIX4eZTtGRTn3a+U0ksi8aeu4P2hsdFOw4qi60A/oCxPTnHfQUd6U
mK7G6AK6ds7HJ1u+oCLuWRllr06g2+juSGm3QQ/4zI3x6weAnlD1eXX0bVRi
qDBAfdsE9EtOmnKTWajkflVOXPgdoF+uOwGexeihwAiLLQsGelHUs9J6Hiqd
+vHfmM39QG80Ks5+xUYV3LhrUfAVGNSW5Ap2Narle0s/yQ8FBpI22pnfiGr/
k2URpWcHjO0vlaTwK+iZ+acQmwwaMPBzK9UcVNGzhgPv62NqgOE8Q9dMd0fP
JyIKmncrAYNTvKRkxVFUz7pK9KrfAUaE67Lq5QdQo17XVYUzfcC4P6Ml/yeO
WtxKLt81NwcXjYC4M7f0UYdV5GFxtSm4VD9O1Y+noZ5fu4Kbz0qBq+c3Isor
HPWVLbBwTS8E1x8fF79IWI6+6/fk0E/zwS3c5kKOwxga6g/OO9ZUDW4zNxO+
hWxCI4RwwuRpMLhHXkn/qeuKRqPGeR4xNHAfOSWz+0UYGvt0/B4/4jN4hEiZ
qbR+QpO0cyi9mwEe7arbfas80JQ9b5Q3ugE8bU1PKIwao5nswqOcrRvBs2Tv
zu4Ttmju/TFcViwCLy2h+3PtRJDSdRIH8AjwivUt8/5iANI7driUbA4Cr/5X
Q9sz8kHGT5HZ9cQLvB0jFbOoLKAmn3qp0fYOvO8FNNhQN8CCooddpncMgKkg
Lb8+7Qcs/JjJcqLNAZNXH2IfcgDkFvUeYH7MBKbQemtpVz7IG/7Sw5kUYBkr
7vWPfgiL3CuttvEmgBXzz/otgUagEGo3q7a5FliiSH88VQEUEyzdPvjYA3sN
I1+eOQCLy1qxXG0msM+apshZy4DSG+sjjIvWwG7H7K9fSQClDwurjgeFA7sj
/KbJjAEotb+Wv44XAbtTNKxhzgWlXt8b3+c0gN2d603jVoHSeHTT2UMDwP5i
yf0tqQywpc+Nq1yuAnvYh3VObhiwIIc2M/UJ4Cy8ShkwYwF2SNnQVWwOHNlZ
RhtGvh9pOxghDACOPDfl/L0WwE4HUhviXwNHcWU3Z8gHsJTEdf5jt4Cj8rBi
ZGQAsLKmxMttHsDRfRnydAn5vzL9fc2jRODora4o6QoCrMZH/2tKDXCWnZ/I
1f8OmGCwfKOjGXAMXOzii6wA+6j0ubmCCpyVIt0s5lLAelrNx/lWwDGx6q/7
6yBgfTdO6oUdBs7qrEKZ3caAfd9fZ0+Q+aazXt53jgM2vFY5ZHUHcNb4/vMy
NwGw0Ql2mqwycMxr4zx3k/mTNbf+/RkHzlpDGapYD7DZBPGn6t+Bsy4utDUs
A5SlCq7HCkWkhx+5vTMCZaeG0pCvxhJn9jV5UkLnrWIkT2hdDJP4iVRMV6DV
vFUduxTSbwRL/G7Ix+ZMxLzVgj9wfxbTJZ7IrsiwujBv9Uvb67eqpM1bY9lL
4xY7aYlvxxC+BsS8NW1VTcTO/hI3/8EMe9s1b63A5i2T+QYSJx/qEy3sl7gU
u/R4X4nEDd37V2nyJW6/vzWz/qTEvVor6ddDJB6MUpU9tlvikQ2ye/y9JB7r
m64mGBJPwjjPzlniqSsjOlY7JZ4ZHRCsxedNoxBfj6+wBc5v05sORxQATTr6
SLjoHnk/mhYma7WBRn0wF3aBBhx1C1Ur4TagyWto8X6S9z9xwbHPVQg0lQ6H
gNp0YNcmcVcL9YFmuD/fldMHbPdGWwvBFNDw6P1r0hKAFb1dkJJXCrTzD/p6
KULw7vd5YummC0uM+i97rloGHo05nSrWE7DkoeyN8iIRuG09lfq0pRK0rZmH
MFEouFDFOmh9CGi/jxbz3JvBWdgarjZ+DXT8zzdUmP4JO91xJweVPaAzddBF
kb4SHKivXmOuQ6SjM7NftJE+vMpvhEo6tcvCKpG0SqzF872ky4P9jX+Ag4z7
tvDepaAzTYkvny4DBynRC7ODXNKpFXuKcMBnv2gnmQeDzky5YeWCfYBPTPd3
3LxLWhDEOEUjfT28IcyZdEfxl+G/SQPF4jY5f5ZirdxpDvh4rE6e3WLSjt77
SsYAH5MPeOSXTLrjgqpfPOA/NVHWGx7pwbf5HzYBPlLe08T+F+jMUbS37fpG
mhkf9b2AtFH+AQsnwMVZb47IPiXNq61SUgd82PBUQqQm6Wg5t9h6wIfqTFov
jpJOJb5NniC9T1QUXki6rFOttxPwH0XbOvIekxYYFLIukt41srmX3OdcR6B1
MwA+OFig032Y9OD9d/YjpJN9zySPk54bC6zMJb2OFjT063/0h+v+
          
          "]]}}, {{}, {}, 
        {GrayLevel[0.5], Thickness[0.01], Dashing[None], 
         LineBox[CompressedData["
1:eJxF1Xk0lOsfAHBZZihkNzRDC7JvIbr0frPNnWhmZBfmLUQRkpQ1JFRIEtGP
pLkXIUuSutpsJSIi2qhEssu+zn3vH7+ZP57znM85zznf5Xyf5zkYcGSbFxcX
Fze2cv7b98ZdVmho7psHKhfP2Rs4V+RI775QQm4XUEkN63wFGUhQL+FIuyQA
dVdXapHWKBJhLiIg+kkKqOTyY6JrTkhsL0/PXXoIUF1MbmrqbEQSkll3W+4f
Bqrv13hvFUEkyWwlclzKGqgRHXKnW02RtJ6lNOkJUaDmkBZ4LDch13UXwgtv
xQO1hKkX3rQHyUye8zKirwL1n87hnYw65H9m04Yu94eB+sHGzWXUAMnvmQ2g
1vYDTTTFMtpcHmHynfKUlswGmrySrzPrCfKX7m+nfn87oKlf774fGIcUJE+Z
Bsq/BhqZd5mV1ISUmI1JXo2tBloYQ+bQ1Riksmf6y4iCJdA+UzKKzVDkvn+F
E5XOBNqw8rHW15eRKr7AzsrIDUCbme3PiwLkge7ky7DuWqDzldwKjSUiD5PH
KzYl6ACdWNYQ7tWE1Jr9ilMfIQKd4rmtcWg7UtczdX13zxLQc2+/ylutR+pV
ClB+3RqgM01eS0b5IPXh7qq9ySFAL+rmfecoizTItTwPNf0N9Kq5C3Yq+kij
19+jj0tHgN6yQd7y7Wbk5Yybqcn5z0Cft+19EWyIvBFunjTTeg42++3aGgUo
yLv8DmaBWT3YjMj2Gt2KR76oHxl61E8EWy6iAQvJRwbM8HcZfhJgy9r1s+i1
EDJsT/3HUicN7HG7Umr+vIqMi0w/177xAeyn32+UXutCpmW4288xL4PDcElf
bWwAMpNz7eJAujQ4vvumby5sjMzNd+6l2CWAU1FfcGMwDlmMbWbcNjoFzkmt
nU+NDJHlAf+w5UlzcHH1bSysCUPWDGUP2xjHwCGF458DT5gjrGdRVVn39sCh
pfeD9LcjsEHoVfT0oh+4PhkvXeVOAG57jT4RtSVwu3LS9vuYHvBkZOSVFBuD
u4u1vYhTK/DWOaD7xaaBIakkcy33J+DwQQ5o0g9gfDnJ215aADhB7e3qI/7A
+OrkmnF3HXAi4+NL5AVgDCDpXrcbACd77EIaz0Zg/BJiGgcMA04DraoP1QbG
bKHFfY+tgLOliSt5RQLKN9oSbAeAu6XWPrZXClDSZns3riDA3fl1s+ZmHqBy
MQKZuzYBrrDAO25JBVD53/fsM7B8KrezSFXGgG7vbvaPTAdco6wmTdUDUOWb
lpE6woAb5U+qlCoHVG8Hw/eBKuANBilnJimAWinvJC18BLwxf7ycDhFQ6+qy
uP51wO9Tq2sImgCUaqFkvTsJ8FZBxiLz6YDaHO4X/p4OeAZLs2itD1DH7Mih
Yz6Avygj/lkwCFBP/uDhn1sBn2JMj6WaA+p1fTA29E/AX2MkK6dKA+q9w7xJ
YwDwOX/jTovXAnrcpN2dtAr4St0l4S04QAODJKmR2PlP1n2gkg1oRI8PiRIG
+K+BskO+foBGeTQG78MBfvCaY1IpAui5KdHar1mAH6lO150UAzQmwnolJRPw
kx/eftD+Aeh5/jN6ZTTAz64JnQuqBjQuPdWPicVb3rpfseoioPFbbxS0NwOe
ZRbfMu8KaELJlZHkt8DPe7Q+yFAT86B1cF4Q8BNmqjchZznWiC5kWtZz7JCd
/N7eieNzVaf4PcbYFthGenylQYbjugRuaew9+L83TpanKL0RY1vQwIBelkhk
W8hbXMRI7Qbbwldn9PfJvmZ783ut1DriHrZFYg8H5/MIsy2q5GwsHyPH8cuD
1a6KAmyLuZ0pkzHh1CvOtb9erC+I4xsU+6SZSbYlNCjZTVpkjtsWVv6Yecq2
pFfHR111f7alBJy7PCdMOM4XOkBQsGBbenfd2UvRyhwnzm759ruI44+Kz4y8
VDkeyC1OcuHEkx6XzuynH+R4ukuYv4XTP+n51DhdiyGOFyeekXmsOF6xXj70
vJzjtWL9wCgJtglcAoEX/gjF3LzWCk1A4PYuzlrE5j96uvwD8hIIvI1D96pN
AQ0jGkXJXgeCQKxrr44goH6hTX2bjYAgtkJRU8Lur41DhWGnABAUx7d1iL4C
dEvGs3X5o0BQtnIdbNcAVGIUXsj9AQTVoszlpDRAhU3jUjKVgaB1VFgBj80n
1xK92sUUCIZ9yyGrk8B4s7k6vzAZCFZv3xGHpYDhSTvLY50IhJNVcT5PvcA9
o7vEIiYHCMEWZ0ZOTID7RSJuu904EELe+/qRzoJ7uE9iUtdjIIQvHgyIuATu
qPjKesIbIMQbbw3ZUw7uKlV0s1Y9IGTX156vXgW3pwRKrwoChPrO2dyyNHAd
dclR7B8CGclpz668Z3DIUT+Zr9obZB5qLbAKZ8Ep4qE2re0oyB7IFihW9AAH
gmbVXxUDIDtVaXBNdRps64PTO+TGYEvi3LcVZg3YtNNFempVgIj1YIe5C9Al
cW0H8saB+FzzUcKJLWDdIndPaHkKSI4HKvw2SgJloWDI96kykFZkG5K1NYG8
O3jkiJw+Zk1nm0EFzIqukc4qmPdNSGXLAtmgu71IE4f5mHQeLx9mg0e6GUuY
Hx6v7PkMZL3FrE0sHyCt2on0nrsEZN2Iy91LwUBaS3bb0TYEZI0LJqWC2pjz
poex/5Gsod/Rdeor5qr4UsNOIKsPeQtytWH+VG6Q/wQzOb04fAhI66o8lJB0
IKsJrBanfsfcVBQgZwpk5RRWKc4f88e9u94ZYkaY3UIGmMc7FxOxenZOWQmu
rwCJJbEa81sGs012Sfp/9qClN04BWUnC8q7CK8whP5zDsPwUG2ZtSHcwXwqV
08LyUzx9B//DHXNFfkHWSyAr9PAwW+cxN+72o9ZiTnhw6Wci5t5WHZ5KzIY+
0S9OYB5D52oKgbzjFzHekQ8za+7xiVzMWR03mcX/AvWaYIM=
          "]]}}, {{}, {}, 
        {GrayLevel[0], Thickness[0.01], Dashing[None], 
         LineBox[CompressedData["
1:eJxF03k41GsbB3BmyISM3cyEOdZC1mMSpd8tJIYZE1knmZ8lFEXKUo44dMqS
OnlR2eskJPvypiQzqVdIKp3KkqVspSQ69vN7/5n547nu63Pd1/29r+e6nmf/
MVQtQEhICIed7P/X3cmTrS8mAhyAIYTbSWovRfzJQ8cp1iLAUP2eepRigUSS
a3NzCn4FhtmT1TOOJshv2lEJmspjwGA097e/V0GSyHk1b7rowPBPdf0w54lc
oK7/pZKmCYyYlGXrzDHkojYn289hBRjpU/mgbYtkkv988NPfHRj1Fm/srrQg
WbncRoo1DxhPnVWS/L8gV6k/qiw1jIDxXj7KfD4Tydd2K04cIQBjvatsoz4J
uUlOt/7DcACYe8siGvSYyF/eWq2kgkxgun9rUbNtRm7ltuwqIzoAM1ha1yTO
GymlfqV1zjQAM2UMvwC/IHe1WVtlKy4Cs0uh2vCjK1JHTjnHVCgEZ+0gy/EC
FKm758bL7joNzga0y+8/5yH13hq4D8lu4Exb7EujGyINuc3x4fOS4GwbFJg4
vID8l/o5OvNVNDj7fSR1RV9GWrQZIe+usMC5ID3QbMIK4ZEvbKklJABLttw0
v/gkwjtec/6ynS6wKI7tW/JmEN6T95PHknuBpTYZlf+6BXl8yqB8G14TWMby
2WeqWpD23lcGN1fbgcU01SsMVUD+l0alZc1JAiulNE5nyhDpXqu3jh3IAdba
38MheHHk1ey3yiJPD9g/mH/60ZIsMqRqx05YlwDXkKnkN0QhZCyVOzHD0oID
/YVjjVbByOSPHMS8OwrcfSo8Gts7ka/CwoMydUrgcX/zuVN64cj3Xc87rUca
wEt/r+LJtQ7kx2fPzXkhbeCdUREWvFiG/AytIFpqzwP7lTpe7ScBWdLaLF2+
iAMfA7UnxEQfZOWuS4Pv+nY4dNTBmkLEIeu67d2mSZ3g6yx1gLDFGYRmShen
L1HB98fyb7gNJYAzej049HMaOOnkR6hZNOC2O7RlHG4ATkYpF4r6Abfz4a3d
b84C50/zkps7fAG3tzQ0r1ERONmeJI2oE4Bjn1nxPrUHOEVZV5WrLwHuvCa5
78c14DRIajTdqAXcUIRr1+w+4AwOF1oF5wJubHxHHEcOOMNhot/lGwA3xVbe
1jMInNElpq/PZcDN246mVkUCZ1I2/brdY8BLkCPsjxUAZ86qKeNMHeC3t2bw
viwAKnat5w+iEeAvSnU0T90E1GjHPvNBzFdCNJjZYYCa1Hs+/80L8DmPT4/a
mAFKM2EhPtmAv3FGXzyvA1ALvdnue6WAb5rM8HSeBdRGZenSy3DAt9iMf169
B+jeazsHzscCnlsIZ8uSALVX9LNt5AG+2232tggJUIaU//mnBMCP8lwWGncD
6r5qfC8O2zdJvZPiTwDUK/omd64E8DOnRVWlewFlz61t/dYK+Lm+g9UPcgE9
dMxyluMP+EXjBtuQAEA5U/5a+m2AX0uXeqtkAKh/wPFO8xgQwU8cPsr9B9AA
mV9jlKZAZKN1q9DxNoGlCkiZKmkCa0rGjvTUCGweM2D8+1uBMwwDjzpHCjw2
9rVERYpv0aS2cgNrJ743aBYG1gcL8sXig8bKYYFvAlexjmciyBNXe5HKuzvE
t0SszvVXIct8SzYHrolZ6fEtpVR9WK3Yjm8ie1j9kEoU39JObyQOfN4m8BKP
8I0ox7dMrkbl9YcifMvSN75tNj0i8EJLIWpfzrdc1idpyvQ83/I2mXg7l3qB
P9UVJO+L51sh+cbfDgNqfCuakDOTzCwEfim8+FyVzbeSrnmYQeRtgROOj6R3
CPYplenYCb/1Ebh3uPzkeLvA45YDlTNLAi9fI04uGAq8ppVlbKbIN0moWjk6
agVQv1ixBNMBIOF3FT9sGgX0YBG66cI9IIk+1dmwiL13946c0FErIBFcKhkW
VYCy5svEF58DSXyI9p/T2YA6qmc51fcCibiwR2PVH1BIWAmq0gGSbMLTI7sd
sP9Sd1Ku0gVICpsYNfHGgJpOtG98nQAkiqa3ldAaoDpuM/Ha2Lzm/kgfkRxA
Ze3t9S3CgaT9gqbkh/UlS2yOkClA0mHO93CxfDExYY2UE0AycIqyTjQCzvKz
+y/HvYG0Y1/s1vWnwPkYoj/ldwtIjrvj55Z+AqdpKW4n1xRIkbqpF2bdgMMJ
kKgZjgMSD1dcPygHvo02TkEfYoEcbu6aECcKPk3h3OwvVkCRHan0KRUCtt85
0c53KFDaZNMrr0WAF7H2tbCUKmwOkyGc13sA7gv20lM2z0BZKnFbhIMMuG0l
LpfUrIPy/ZiXPfr94OqltZtw4yKo+IxpIvlVsL9zauDqNB1U1lRktlceBJb/
5AGZI7+Dall/QlHyWWAuhtFOoMVAdZS6rjJRAU7pA6slRpFAXXqyQ0qGBXQN
3YpDOjjM/ZtxoilAVx99ObdAw/xtdX6Rhzl32ap5BqjLFO7AyHbMREeCTjTm
MGZFLQXov8wtgPs/QF1ROOzkOgJ01QehbDYC1FW/rLScSKBTnD/Rnu/DHB1z
Nq0Cs7jGpw0emNPZkWc/AZ3MRRmL+Zgb1dnBHphpw/RNVkBdE7+7zcIS6CTK
+HSdKubqx8/6NwBd4aNSl1cA5vbbD3sQzPkcQkQb5veptbxozO5lKzomQF0X
ZV2/MwV0+Q6rTUFYf91j4EhcN9DlKk/2rWL3WQ9tPRQhhjmI19eriznxhksg
YFaXz6CFYL4TvJNRA3TZzAadomOYHzka7pnGzCBI3ArD3GeoYaaJWczrlsk0
5mlZJT0foMu0VizarGFen5egZmOOxYm86/wXng4D2A==
          "]]}}, {{}, 
        {GrayLevel[0], PointSize[Large], 
         PointBox[{{7., 45.85305487675453}, {7.5, 46.278126718300705`}, {8., 
          54.1365314659016}}]}, {}}},
      AspectRatio->1,
      Axes->True,
      AxesOrigin->{0, 20.},
      BaseStyle->{
       FontFamily -> "Helvetica", FontWeight -> "Bold", FontSize -> 14},
      Frame->True,
      FrameLabel->{
        FormBox["\"\"", TraditionalForm], 
        FormBox["\"Temperature (\[Degree]C)\"", TraditionalForm]},
      FrameTicks->{{{0, 
          FormBox["0", TraditionalForm]}, {2, 
          FormBox["2", TraditionalForm]}, {4, 
          FormBox["4", TraditionalForm]}, {6, 
          FormBox["6", TraditionalForm]}, {8, 
          FormBox["8", TraditionalForm]}, {10, 
          FormBox["10", TraditionalForm]}}, {{20, 
          FormBox["20", TraditionalForm]}, {40, 
          FormBox["40", TraditionalForm]}, {60, 
          FormBox["60", TraditionalForm]}, {80, 
          FormBox["80", TraditionalForm]}, {100, 
          FormBox["100", TraditionalForm]}, {120, 
          FormBox["120", TraditionalForm]}, {140, 
          FormBox["140", TraditionalForm]}, {160, 
          FormBox["160", TraditionalForm]}}, {}, {}},
      ImageSize->Scaled[0.585],
      PlotLabel->FormBox["\"Temperature Profiles\"", TraditionalForm],
      PlotRange->{{0., 10.}, {20., 160.}},
      PlotRangeClipping->True,
      PlotRangePadding->{Automatic, Automatic}], {168., -182.4}, 
     ImageScaled[{0.5, 0.5}], {312., 316.}], InsetBox[
     GraphicsBox[{{{}, {}, 
        {GrayLevel[0.85], Thickness[0.01], Dashing[None], 
         LineBox[CompressedData["
1:eJwV1Hk0VesbB/DjqDh7c4vI3mcTRdKgMktpP+kmc6KUMoSS4USTNAh1IxmS
IZSkFKmE+CHltN+igXsroYk6GS5uaUAZivi9/thrr89ae+33eZ/vsx6nYO9Z
23k8Hh8/jRNv1tG7Wy/T4RRnIFPpdtPUlvV/PAcle2hza54of1wecJcNKQz8
sSD2GOckXL3eLK+F/StnqvywRjHn7uGxTLy3mI2ROZz+0cGP89ogfW/SgQVs
vOi/Wc1h6pxPd6xJuHE7m5zZek4S+Jjza/MVmq2dyaY+S/vwmp/DBdqfHjUq
j2XTpey1np8N53YWP++0qJRiM30rC9FjY273kkxa6q0lm0PVeq1LcuZCuG9j
RbYse+u1eU2cFs39lejZYTLoyd4vDx2c+yWMSzZUCX59RJtt8PF3Cktw53JS
83OR6Sm2xbqtUfpPU67A2j46Nc6ebZdqsrhTVc+Vbp9jkt3pwnaJvcTlAwxX
WVaTs7r1Ktujo/TUwjyKQ2pfbD+9z2C/9e0eO335BFczOaLQ1W0N2++0TCtg
mx/3ZKViSYBtJPvjL8fpjeMM90/qPxqFiwzZwY47Suo657j65o8Z/2nbsj+Z
csmtZh2uMVRZPjdrJztidcKIC3TlXi0P0VQ21GR/XynIKoxz4d4eeShxiUbs
+JttzzsNvLl3D9NP79C7D1JqlYLwMjWu1d76/oGIBOBbX9ttBVu4dudOMT2a
DtJePiuMEi9x/7o+H8xzi4dJpx+npOTZcV0ldJv2rn9h8o0BB+XIce6jTHbk
4sPPYIrkaZKsJsX13BbvcJMUgKzc1HOmRtLc15iDMd5dUSCYw8jGEwNc79V8
pe87WSA21I875ttz/Su6Gqi6QiD3OG5ZKrjJfa8njlpslwa5k0uPeDRLuIGd
es+E06RAvuLp3IG9IdzQqsE3VUp68EfDfJOieiXu57SS9Tkzg2Aan+r+ckiO
G2GyFCbFvINpGevSZ1UPcb81X7odnn4LFLSChru/dXLjvDF5X5gNClVeF7oi
ohEvYvnHaPYkKNqYDu3vW4Ck/M/v+2r+Nyi+XG2gLaxEfDv/vgB7MUwXaZfK
X/2FpNl9QtOCBpg+/LnzhAqLJqmf5Eo1p4JSUoyb6ScTNHmqWK9m7n+grCEb
FZatjyZ/uXxeW70VlG+Ikz7IzUBTHiqJrI/egxmr1rQffHcTyRTE+vPyBDCj
vrwzL9EKyR5LM6hZ2wwq/iO6ZG4aEvhF9F0q6AEV0QaFsv89RwIRnysvuA4q
QYWDnjWySBAcFfRglASVvd5cafshJNgfG9wiugYqYXVOburuSHA8eQzVzQOV
xLOhBWe1kODihQ0yqo9ApczkgX18CRK8ufBHm5s6qFQk5Q1l9iDBOzVZ5xj8
fWVP3KUb+PsPmR0nVr4FlXvZ6wdqcT1dZ1nblk+g8mRK94UpB5HgR3L8jxJN
UGl+Jfc1EhAxLSI6a5cxULx9GxP3PkWElcW/P3kzgOLnPpeuO4UI27BlRgtC
gJr02vLALAdEOBT8iiSuAiUwM9r6rB4R63+M78s/BJQib7r+vCZEeM3SfTvw
CijNuGdNH94hIvSXRri7DVBzxJbWxlmIOFRyY3x5AVBzv4rFCR6IOOJTbfAd
n79wXUGeWSsijhU1n9/gDJSRSmzomQ5ExPelbJsiD5TlZUuh7SdEXJzX0uVz
ASirJvGpnBuIyNkrHFkXDpTtFCP+TxEirpT11nVkALXWX7Mn7ysirulk37XY
BtSmxVJiqT5E3KqucR2RBsr/jtizYhgR9yXhLn2NQAV+Nmz6oxIR1U98BOme
QAWpFVhtP4SIhwX1Z45rALUnMlNPcRQRtduPeERHAnXY8iA/iIeI+vTsOyEB
QMU3GOZqySBCQhLq/QRQiQZS891cEPHhkXbn3b1AJaU8LUrNQ0Tb4cazq28B
lbbB985kS0T8+1rVTzcYqOw3Gc+6jiOiZ63OiLMfUCWS0eFrfEQMOSW4FvgA
VcbWHmlzQsSwuP2RwVagKi6mSdM5iPildVum2guoKu/F8rErETHayeh6dgL1
sNNrdmAkIvn65rqxSkC96XlopzuGSDkBo370BFAtdkkvfB0QKa9/XuuFBVDv
b3q4XLiAyD82+Xf4TgeqPWjIa6o5IqedpVLC9wDV0z8vpDcMkUrffvdZ7gZq
bDghu+QnIlV/XyxmjIGWct2s1WOFSLU+taWes4GWvqN9XTMDkTPbksa/ygEt
e5j7X6opIjUq415sHgFaYayvNmQ/IrUsX8o5ioCeM3nDd5MfiJwvKdYJtwR6
rnmZRcU/iFyQ0b2W1wn0vP1KySZXELnQccgqTwVo3a5GPRNnROriHJV6gDZ+
4hhsXILIJSFZzVuWA20VZ9NjtAuRRgmdlTHaQNtUXzcrw/UZ65yvNr0OtN2o
INZIA/t+wrtifaDXiZ7MM3yOSJOvRkf9xEBvdli9w3ARIs0MDCqqU4EWKbAd
+p8RuWKXUY7cSaCDrLP1S2qwe0fHsp4Avevo2DH984hkd63lEeVA7+sXz9K3
RSQE5Dmjt0AfaTLz1LuOSAuH890jWUAnphu2LPFD5OrXBrJy+Pyk56nziwGR
lhZXDLMGgU6V+X5oCYVd0L6vQQJ0xoESevFjRK458krBPBjoS5sXb1ysjUhr
hRHF/TFAl6rNa9TtRKTdlLmH5pcAXWbjVqSD67HzsHY5OBfoitDEOM312GUH
Ug1xP+++GFhFPUCkvRepatQIdE00KpW+iEiHm9aP5v8G+mXvxuS3bohcp7rn
+hPc7zdqsUFNSti+1/e99Ae62UZs8+xv7GK67BmuT5Krya8xQ6QT5BNed4Du
dv26p0gFkc6bkNQrXaCHq6OcjjcgcoNvlXsDD+hfvbcXRcRi5+6q+zQV6N9q
PcShldgdEVa+vSCUOuD4ILgYkS7upQtXHgchoauqvzkBkRttcp8NtYOQSStV
XILzdJVdkFeaD8IV/q0vWjhEuqk5Pb/wDYTwyO/e7SnYDj4ra9NAuGp2741U
O+zwz8H5t0G45h0vyu4NtiTYfKMjCNc5aJiKexHpnu69aekYCLfpe2ZfnI1I
j343G+M2EPqe6o4Pw3l5MqVNng9A6NcTfHBTIfafWcqsBIQ7r0Q4KyzDPjPv
V/dmEIbOyJ5yDPd7q5799iEDEJ78Jdm5LRqRXi5VuobqIIzf6LsZ8Hx7Haxa
pDgOwlOlX9eoKWKf3xruuheEKaIxjZfZ2K3n8vhaIDz/Qa3J8jYivbcNwOTX
ICyqdl+m8xGRPh5Lykcvg/Bl3DvZL/h+26m34XIfQPiaG0QKQuwlaY7rL4Lw
7Q+FUMP/sK1U87Px/SVulp1hUdj7T/yMyARh98JbD8h72PXdkZr4/8NPTxye
vwSRvmHmK3pOgXBEKkffHu8P39P1qpKXIPxtVPVxN67fN9ewskkdGKms/o2V
/thPQ9hJgcAQQe6GVjmI3MG0NS89AYzqVP0vO/D87ChuyFl0D5iZq+yvxLZj
P7hTLcMHRiPUb0shzntH0+GqtRnAaH3Iqh20xx46+WBrLjC6xbJ5Ubh/fsta
O+ruAMM6SjwvD2OLX5TzDwLjnXzyRet7RPpfM+N7XgNm28DPs114Hv1vz+CX
SAOzw9XfqwfvI/9H9XfXPgVGpGHVN+iE3dFRP1MLmJBCWUXyESIDVGPb34wC
E1Mb7WSM8w2ImW/++RYwcQsHhctJ7DNMykVXYBJOb29fiech4FJvk78yMMkb
V++xw/MScNfKRKAKTGbnpGSvNOzPJZ6LjgNTxDvWEB+ByECb0iVeOsDc8uk/
l/QO22VST9QGYEofe3un43oDvVb82rMLmNunVvZfwvswMFTUSwIwD1SlFCvw
eYE55dDXDcxLkwindkfsgUVFfp7AjAQdbjSdhUhRnN6yujrsyPvilLQJ7zut
yGCnyFz9JjfhitYtstgVKWG5g9jxK/b1fwBmlHdDe3oddoJtkMVl7OTmA9/2
YCduu3T5NDC/K5bOzMP7TpSSKr83Hrs2QsBfOuFX7eH+2C01391xv0SpdGnq
YmDGeI5PlDImnG3dhvMcs/bbfUyEfeb6Gjk57Jb0Gg9l7HTur9Zo7C/vCyvx
PhBl8OmqMmDGeZpnlXkT/vNG/m5sraKdf3+acG11Ec5/fOdjlaV4XkVnG2ui
cD7jkfJSZwywz80w68LzN57s1NObP+FNN10nY5dL0NXkCb+PDsHzMF6rdUNa
BjtT/cvMEeyWgDMeYRP2XtvyG/tLcURl34SvFBY4Y48PBij7TrhbcHrN/wGe
SAb4
          "]]}}, {{}, {}, 
        {GrayLevel[0.5], Thickness[0.01], Dashing[None], 
         LineBox[CompressedData["
1:eJwV1Xk4FlsYAHBJJBEh3/fNfEqEZEklyjKvimSJZElcSSpL9rW6oYRSV1FJ
opItiixJJHOyZatsRSFboWwpZIt7/DXP75nznHPeZd4xdbeXOM7BwcHJwbGE
b/FJmdj/kj5bcIHeshCqt+X+Q8pleVnBTU4LWi/5+0DClnWU308ehQkRX9r0
q8iZ4J4M6mLvo+8bEkJp27dF/l9zNanLd85Gq8vI0fbbvEbWUzpUpPF+dZPs
etphxN7J19OJutFxuvJVDz/trEqkbd00S8WoHSP0UrvoU7WiLdkr1KjYG4bu
TS65tLtFRBuK9qAS9q1lDEwcor3kvbdPWW2jkg7bq0RwHaD9hW5J6Yw4UXmM
j6d+inrTocsYP2PMxKjSKcFLSeJc9E12o9MugUaqbF/55VO6yfRN71ZwSlGn
yu4GRKi476Zv1og1fxdeoMq1uq68QRfoW/57A9ruclMVodmRQ8c46JjmoN3J
g8VUlbDJLZW0OfpO3GQ97+mD1FvF64lVyuP0fRs/1Yur91JNuZf+8Xv4i07p
F4cIvXGqqTXTXHLiD53K2Foo2lBLNc03GjbozdGp+zym9yc1Us0GpMamUW46
9ckxVxLpUM292axudSad5iPveFbqGvVRpKXVoAnodD6PDy86eqhWf5mDEkv/
o59YvbNX8JSmOtQr9WodZOmcktnuF02yVMfRTXm5Ial0Tvu/Fa7tKlRHWJT4
naQNdM6MqMbPukSqo8F6/GTPejpX5b9NWrpq1JcTYw+47Eg6NzOup/qrCNUZ
Tf7RtBGk8x5qXD57NpHqGt/R2JI+R+ffW7X+jnEY1at//fy+xF66sDdtVaNz
EtVr2zpksCGZLlpxf82GM+ZUr7eEpWG6A12kXNBWYLSd6r2bq2Cc/Y0uCjaY
07cIpHqHmj6Z0/30y7Uc08JPsqmv18SUj3YM0sX2OYJHRSyob00J3WeYv2l6
adlpXcqZGrB5rJ0RvZQuPVSuMm38gRpwvyNr4itOl7qLlr/TMKEGLoQLTB7a
SZeG+i8QxzOpgUf27bvEPejSXLe4+5bj1MA4y+9TRjtdxr/hVDPTgPp+5fKT
5WXP6LKK12eklrpSP4oc15yYcKArdv+WNbBUooYZUoOShyrpKj8tyQCKRQ3L
Pwq3VOqhq25oK3MkRFDD2vKSV7gX6KrswxcdZLWoYeet1r/y1eiqH/Vyf3yK
qOFi7bpSoXS62q6qRJeZRo3Y2WbZ1Vyla0ymdxZluFGjj2I972mY0nW70t2e
19ZQY2To3Mpro3T99u+pppUl1Jj07b7xkXi6Xk9QJiP/JjW2Ob2hw0iPrrfa
VXtr+RA1tudtWtaKh3T92brrU0q/qDFX0YPGoWZ0/WuvjLbSamrsVXJG9JkX
dIMxf+dYAYv6ZVtmxTh+gW70ThV6bDhO/U7keC6lLkI3915+eFf8FTUxI/f1
TcMV+vMSVXs7oz5qKoyvYuPAJrprxwCfsqM5NdMXzxL7waZ7s59vrStTo/6q
b334MO8R3X+5bfOcUii1UHq+ffCvDT2kVRtu9XEBlqyiKrruKNM/b3w6WuuV
DZwWdyblvKvp3+HnpLQ4pmHp7T9oV64CPVm06YsmzxRwlWXuedC3k56eVZXK
5R0A7uUlEl8O29B/30tcCP6XBB416cKj5TKIo/ldVd25QFge/CO93I0bLeXv
3lXQuBl4y3K3TpzNRcu2z2xQd3oBK342BJwK4kE8iiHxqnfGYaWqoWhpqSfi
3ZPW2MPKAn5HcVfb2RDEpxBte+1gEghEj9tWx9ojfq7MjIcnRmBVy4Z0YV8L
JDDra+yvrg+CF1osGpKmkaBD5V+7GRqEpKtDLyuYIiHrtCN+wxUg9ObVeZmB
KLRav6g0sMsZVv/TYfX9XjISll54cum/9yDMkdHJGi9FIsJE77C8NQjHJmu7
6zxDIsNiYp+H1UFE4b5HRtVVJIpSPI4YNYDIez1D80F3tCZpfC6sIg5ET7CX
SfFKIzG/lfpPR0RgDe9956OV/YhhWWJVFZQEax4aPSkNOoeYsulb/xYZgpjq
3+IVTtOI+fvTG75zIiD22dr/9XgwYr27YNA5WwFi7YX3NtjGIFbDU11OIh/E
OsXKL1U9RqzmdlkRrcMg1tssZBz/EbHaVJLZd16D2LBhRoeOPGL19//r/PYu
MDg0O6ZutSBiie7jZc+EgMEZz2W9MIQILjd67TtVYHDNyNHOnIjgiRK2ZNwC
Bu9z/zAteUTwl9+lIn4AQ1hJSLjvAiKY/T/v76gGhvS6PUqqiohQ3jAwGdoA
DNlA5+jEPYjYOmGXnckHDLn265N8VojYTkc8U8gHhlJsR3HnRURo7BGyLcgF
xg4hf/3wT4jYK5DhZtQNDEPOjOOt+L11tsjnoFRgeH4TjD/dhojzB6tPNV0G
ho/vcM5DD0SE5NMzE5bA8OeuqazlRkTY6nMcOWXA+Ff24hi5BRERxadS2AeA
EeYyrVNyCRE3mgT0aA1gxP3qHl26DREpq29fjBEBRkIIzbWpBhFpUl3HDKKB
8UAknml2BBHpW3xKdEqAkbLdYndKBCIyqQd/Jm4A4+npmti93YjIXz84r6sD
jHKOPO2rkYgoF4Hv7/F+lVHXLfKlEFHRuVzXhRcY1etdXToKEfEm5a6cYzMw
3u2RiVH8ioha6V0t8RPAaA2/+6NhJyIaJt3/YSwFxrDAxZui/Yj40joQIlkH
zDWked99QMTE54g0IhKYjMSdSg4JiPgjozTq2w1MQlo8YOMMIqY8MnLye4G5
Tqmf91k+ImZmdl5uOQhMud0B8jVyiFgYCGjccg+YWs5xXhPCiORVXWHKkAYm
jAa+LPJA5AozHp4oXWDu9jnGFfQOkXxumo0fOYGpF7gpZsUlRPLHeXbnZQLT
NKq4cO0cIoUaE44dFwXmiRdfFvT7EMl8s0HWiheYTppleoK7EckqCp/30wOm
S2na9Q8PEEk8Pm+kVQhMjzqP9XbWiGSHGV1yeQ3MM11LdHzrESmxtv+LhzUw
r/FIXnlQiEiZ/g8NfbnALDB3FJu8isitwzf6PM8AszBr9sPeZkRui+a967wK
mMXLr0XfYSNSZdtPBZY4MNHLAn71LERu9xPfqysMzJr13JzBeH+1PvvZg5bA
bB9JGuYTQaTm+aST8znA7NRTzbCxQaSWQF19AQCzO7H2ZGYKdlxCUEkFMPvM
fvWYbEcklaW93/8EMEdfarfcOoRI7Yof8y8UgcUZ3lm6Nh6ROoUhxnfHgcXV
7R3k8RWRulti4qUGgcWjzqNRqoCdcf5kuwOw+EYUCxxKELk3ZsDj2DtgiZqd
e5Lehch99nojTw2AJStB3N4mhUjDioLb9d+BZVRkeUrvCSJN6z9732wElolu
YtshH0QeVLb+rTMNLNPGIX1HTezrX04UcgLLcvC8bDiun5nhT/2Xp4FlT2Z8
rRxDpHlh5LiQJLD8Auf/0cHxHDoxumWKBtZpPr06swXsF0+//doKrLO3b6g7
VCHSijdO8F4JsIKzZZkXrbAf/c69cxxYEd0HPpT9i8jDHZO9tpXAur87Zb/2
a0T+IynM35oCrMT3o68ORGA7OkXIpQEr2Wan/FFT7Ew+DQVLYKX7vF8RjPNj
u00tdedGYOWlTFUibkQe2enZsukTsKqWG1Jahog8KtNjae0HrJ91vzfv/IhI
h+j4ydcrgfVbYVj4I97f4bWHjI4PsCYi+/54aWGPNtFPDwBrxuRzSXoqIo8b
mBRPvgJiaXOpEcMP+6/2zGPsNe1RTpNrEHnSaKougB8IpuYVo6gabD+L1pud
QBD3QpUVArHvWax0HwJinV3AlAPu75PD6qm7xoHY+NUurOk5Ih3DXpVc9QNC
Y0gpKccCkU7prkNzwUBQRhvDjVZg1xld9TQCQjtL0nkA19tppP3Bg1ogdD3E
tqyTRqTz1pWKPwSAMP49jyInsQsqmZpSQNjPvu04FYNIl7z9HxJXAnGZ75SY
zAdEuvo/31eJz7ti56poZ4Id5eintxaI//LddGLrsDN6bDJSgIg+4uHDW47d
ofRImgXE3Tyf+sFcRLpR8DPnIxBZ1ucinkZiT2cPMtKAyM4OTBzA93UXFC9d
xQFE3rKgwnWh2NI2mpz4vIKn5/ujcL3dTfd2HC4G4vXSsD3eLthpw3+l84Fo
ehw5r6qHSA/9lmPtD4H4yHFd1B2f72HrXfdzEohW8yj5RxS2Z1/Tj1NAtC/c
sGbi/vOINT3DjYDoO3i7YBb3u0evoY36aiCmZh940ksQ6eklsUxrJZDiBjnf
dIsQ6eUbcCSYBnJdk1OWJs6/V2hX3G1RICWtJfy34fi8YqSyzPuAlHWO4l2f
jv1ce5j/BpBbLnvIz8dij7clb7wIpO4bRZ/nuB+8XTqv/pIHct/+Ps1ML+wz
yoUv1YA0+HiPO9kN+9LxzQJdQJr08cdGH8dO9lxV6AekNddQsZsZdpuUalM/
kO67HnNJb0Gkz24VwU4ZID1r7N+yFbGNi32n/wXSx5QVIyKHbbMxUF4dyNNH
I2Q5JbB98iO1XIEMCXYy+CKInZyw7205kLdLZG7cGsFe6LMuDwAS7UiWXJaB
SN+UD5loEMjSGD67S1HYWV3n7skCWf7L6x5fAHbBZ7vK5UBWP9ZmCuliV8dz
1GoD2UR+WUX2YA9Zc18IBLKPQ2x2G4FIP8UfL0Lngc1XeaXJHs9Xv3RlZ7GT
wBZY/1vwmzd2tn2LrymwBYMO73c8jF0Q7hgbBWxRNdlqN1nsigzHbcXAFk+v
oE/jevt1H54ZvAps5asLj6//RaT/miBbkAW25UHvEBrn0z8gZHkzE9hWPVHe
HLif/QPfSHKlANvG86n9Lg3siwsZKirAPnptECr4sKP2O6rNAdulzv5vzSPs
9ONBmjeBHaxj6v8Bx+Pfcn3uHV6frqrsMmSOyAD5Cz7K94H9uNLYWqEYWzlk
x2ZRYGeZu+q74/4K2H52R0w0sPO802XHfmHDAbU6fWC/yl737Q/+HgIO3tZV
0gB240ZBO048/wIC4vf6vwD2LGvUjInvF/CKMyE+HFsxxDVvfNESiWFXsHeJ
he/PXDQkfXXCdqKKLoovOjBpDsc7W3B93SjHomeiEkaBPWe2ebi8Arvk13qE
/fc/r3CP/dioi630ATuRO5GPe9Hzo+U52M/iilJo7NdsVB+A3VY63LZ50da2
RnLAnpcTNtNdveiPYR112JXPJAg8LwNK6yLx/5I9/1lv5/PF+EoH9YIvYQ+3
mx3Yi122YslyXmAviHKFh71Y9F7X4CrsY2bDY3GLLt1GDGD7fee5arroru4x
wI44JyHDt+iFSKoeOyfFzAbPh4ByjR69aeyKHa5/VBZtHbG0DLv1bXjUyKLP
bFaTwR46miifuujY5j5N7IWJojdHFl0QICD+P9g85Vs=
          "]]}}, {{}, {}, 
        {GrayLevel[0], Thickness[0.01], Dashing[None], 
         LineBox[CompressedData["
1:eJwV1nk8VF8bAHBp05R9GzNmJpn7ixbK1oLOKS0UkpSyhEiylBRZSlJkX6JI
kaXIkpBCyT0YomTfXoqSPVsiO+/x1/18/zjnPvd5nvOcq3/lvNQFLi4ubi6u
FXuWn0Dv/ALdZKMLqbAw7Mtn/gvYH7XIfJjsR2qFnL2Ur3UJuO5MvFgUZ0ue
fK9XdaylCvgSBzgnacak+bawBq2CfBC41SMoaGmEtPzPMyQ6ugeEUdu2WDoV
k9Ybmso2b/IBD1lC47FnDEi71/cUW68pgudP5taV3NAjXcPNtxR37gO5RldE
T62pJn2LrEwpHQyQO6ciZ6LhRd6/eDHObvQqePtk4YillyLpJ1DLrrF9Bd59
D3S7OhNN+lsNaHOmvUGBecr3kCELMkioVqmx/S8osv7xoqJuggzjvNIQzt8P
OHtS6UzdYPLRzXCaQX0M4ITVZE8aU8lHv/SJ6bJhwOn7p1l1MZGMOhpqqJ91
CpQ90nBx93pHRtPuvpHhlwDl49/rmnM6yMfIbc2R2IugMkPIL1RcjowlMl8y
DRNBDdNjcqmzlkyUy5XrlbgBapwTg5qHjMhEu1Xm/67LgpqqSulX091kYsoL
lYgmUVCrHiiTyTNPJkm1debs3Axqf7eeNtkkQD6nsfOi6PKgXtMpO99gN5lM
e3Enb5YHNK1IsXHM9yfTzg2EaERWg6Yths6VGvlkWkJwh0WqEWg6ufauVHUv
mdbtm8zbSANNzy/G1vccINPtbuyPdZ0HzUc2NygLzpMZt8xhRwUXaAl6qT5j
d5nMzH4VUC3bBtrEU4W9mfpk9usQ0TwXUdAGz25sS/Uis5tDzlj0FYK2S+vk
FJUyyewF0lNBIxi0vb+k2X1sPZmj/fZopl4faDeV9Tx8g0PmDE18MTGxB98S
UwcodSpkrspfeYdYCujcnlYUcY9Ovht3KBXl3Qs6T4W4Rc0eIPNY/NuVex1A
5y0npaeOtmTesd4d+uHyoPPrnrQXJvlkXrJyQ/ZGN/BDxeqYcUw6mW9qoef2
Ohn8aPq7y+GaJlnQLP5Nw/oP6BIV5A9ne5OF4x7GJ86cBl2y3MPCZo/Ij//t
ctot1Qm61CY+P4xJJT8aybfLZaSDLqtWnxiBWvJjSfMtW8droCv32VziAoss
iuKaZp65CX4ZyPXlNhWSiLh1NXH+KeiO0CZbff+RxX+8H6z3MADdL0fNRmMf
kSWC3QocQUPQXfSAa02uClmi4GDFoE+D7oEWqPTThSy5NrUmi98M9EALFK46
SZZM/95Z6jUOeoZdirXGJkgOxcXpc0cc6DsSX/LB6A9ZriNsHjrBAH2mGufr
r4aR5U4Z1tJajaDPqZd7wG8HWR591aV2ry/oi9uqIf7WkSzv9tLyuSQH+ibe
ll7nHyM/eeYmhNUMgP6Ez5xtnBGy4kPUDblDlWBgbqIsTm6I/HIGXrrxfhwM
ChDl6RVc5Bevgny55kIwyD79qcBClPySanqYsj4ADGrnVzZG7CO/zEk/GJZX
B4OxN6spU2FkVXwW3ZWnD/wGq5pvICXy69+uFQyZFjDkLdJ34oQHWWuaf3WL
rw0YenSo/9xgGFl7PbfSRKMJDKW5DNjffU7WBlbwfAgwAUP1rb/vv60ia98r
Xl94R4Bh6aejRTQmWUdlBTrE3wbDZezpLV0lZF3rmUFlgXNglEd53arrFLLB
NisoofIyGGWKVSX3p5ENvnWzOh+fgVHF6TAt02NkQ5KAwcaDSmDUrFA85HAQ
2dDezWdcrwlG8zT+o0rwk426A/YbXQzAmKJMwsNOCtm01/2gx1IOGOtSKUq9
J0C2yMjQApeUwDjUn66tppPtFvJO3x+qgfFjAuW8K7eQ7fcslbgrcsD4qeoI
rd27yfaUd0/6EobAuN3R7cWJp8n2kcKHz7m7wPij/edzXMLIb56HX9alPQfj
Q/JVkaxV5PckM14oth/8jd4Qb3hlhPzBveC8cZ02mBgtP9IhiMiuKKSY8msl
mJj5cuih3zqyK2Pdqf70MDC5slZDe1Gf7EIW/sXEbjAp3rbvw2A32TUAdpkE
T4NJMKL8mLOW/KWaJR5qdgVMholJn3LRJX/9stN2VxkA/3ZcWKxq+0b27Dsf
mBE4Dqaurn77IWmG7JdmtH/6FQVmZUxdlFIKyNHD2zoO3vwJ5rPvnLo1y0v+
vSc2+vAWD1ja/oPfZ+4POdVyzrxISRNyjW0NEcnLJBdYJw3/eodBbgUd6ero
V2iF9ru9jXdK4cpQ965mBUe0Kuzss9ZtrnBV72F+dek5tHaqveLSmZ1wDTF1
bINaLVrP94JCKTODa71ys2T5bBDf8TAUs8If8nDOsDPXHkUCiXGnORZckCK9
tsQrfwwJKTa8UO23hvgMfs0pOYeE3940tawTgxsK7/dr/xVDomIFZvH9wpCP
+mmb0rdwJBZ6LPCgVCbkNx26fW/4NBIfsCWfcwaggK5hPG3LOJLYzevjc5gB
BeaAKPXDE0SLA1+PfhGCgk8fh38xG0D0fr2+vlg2FDrmLs3+EY0YWlU0Rpc+
FJqStZPj3YqY8bVjPLMtUDjq2PvfOpaI9f3i3K36fVDk4ATr+KUFtDE6+eDe
KQ0o0scfeDFGDUkpl5nb61yDoj4zmfIuZkiq+hrj9LdGKKZwqSdOrhttcsrS
t37sDMUadNhojg9JU170bBPQhuJ33rkGJD1G0pm843q76qD4vf0D0zEBSDrL
oXjzf7JQ3PersXW4O5LO+Xor/U8vFA/s3nfgtgmSzgsiq3aGQ/FHQitnjVlI
uphiU1IcAsXTHEJthV4g6UZuPcukYCjewE7WupONpGe6g51Y8lC8KUs8zzUR
Sc8foF6ox/u1qvoTVx4g6cVn4XaRz6H4d337FabXEXulkc67s5egeJ/XzoI9
uxB7Q/UYl6EvFJ9r/ygz9hGxWbnfD9n4Qir7QSOPWSViHzQq1C5hQOpmQa8T
oraIfdh74qJ4BaTKhm17/IWC2JovI+2aSiBVLthnq8pRxNYeElFsnoTU3X4q
uhuqEPuUwdHjR+5Cqvat6IiCasS+kK694XkkpF63MWUJNyD2vTsOKvP/g1SX
vnUXK3F8vkbOnxtiINXtwtus2yKI7beTkJGWh1RPS779gwaIHdjitdKcD1L9
TEkr1IzYETwm90bcIPXpSak0h1bEjp83t4wvglSOeq/yp++I/UGVu9VjEVLL
x9PqqCOIXTgS4Mm3G1IrU67Y23IhdtGzKHrQLKTWCEwnbWAjdvGsZqv4O0j9
X/d6AT17xK54cCKwJQBSRwPkB1tnEbvpuOCdjCtQgtbs9myQitjDG4RPProA
JRiB+1RVtyD2yPUO5a8boMRGyN0cpIbYo+3uJylLUIJIDeKVM0fs8UT5J3MK
UGLHzTh3pxTEnhatmdwfCiUObyo1mFNBBHfUhMKGIChx9QovZcNJRIhRqQoi
slDi2mQsZ3suIsS39tu4F0EJl5tynsdFEEHdJySn+wtKePgfHw9vRATN1K1C
tRhK+CQ9aBc/jQimX1UTZQxKPG6lZmwyQsR/ESz99f1QomQ/W3e3OSKU726L
dTwLJTif3vKcRYhQuZBHTzgNJT7pHirxkELErsON7dHfoESV8UWVj12I2LOq
boexKpRodk5jaVghQu3K4uN0fygxmCY/pnsRERqr8jKNWJAmLKz6wOoyIo4r
Rvv61UCaaE33MymIPXdAau1KSKMGBad3CCFCr/j2IxM+SGOu7iwzfIcIfa0i
P24uSJOZvjt9dBERp7Tuc7wyIW3f96pzO4MRYcwjWZ9/G9L2xzjbjpzDLlic
yD4CaQcNmTfSdyDCxPZhAViCNK1ax3A2zofpp/65ALzeoFSEQ5VEhPk1U2aZ
MKTZvjSVXcxAxIXQna91AaQ5XFij8sETEdbbNX3U10Oao9RrDVc97M+KGlfu
QZrzkxUmfyYQYcMVdJTyG9K8gl+E9qghwtbMu+iyPqQ9chqe/FKFiMv9Yusl
4yGtRN2zOPo3Im7sUDlZuBrSOMPnU0XssdP4HL7g/T/FHg4LHUGE66Z5d4Eh
SKta4DPzHUeEm9BGQy0/SGspjF24Po0I956siJjPkDas+nHPiVWI8LTcF7wJ
53P0d8LGKl/sms12wdyQNv7Eh0dzLSJu71UfNLCBtKk57RZAQYQX70zuUQNI
537/zXm7ACLupDf1LxlDOnXPXPY6nI97Jbbnje0gnTbQGX0vDhE+Ut/+BFdC
OuNxqRcXC9srojlzDNI3zQTqTm1ChK9qQ0mnN6RvL6AN9cgg4n4qWtk7COka
u/fKFCshIsDqkvx0NqRfVnGNd9NGRKhp9wraEUh3zJdJLsPvC00y4pEVgvRr
e1pfCQxh9+tWD+J4XNV3F7wMQkSYY750uBake2tM17d8RUS4R2LC52JIf3TC
bbWKLiIi3J/mfgyG9CJ7t0t/jyMiSlrFM46EdDQi6wg2YhvtXM2QgPTSK/9z
CRjDDuPkPxOE9Ipre+5uCsOef+8P1kJ6vdvM0xO1iIiumeDuOwjpvffdqrNw
vWOsQF/MFJTkS3RXdNBHRKyWOvIzgZKCRgeCDCKxXa25OgagpIgQpUe1CTtF
fnWOP5SUuPMkinIGEXErk/d1/4GSbPPC2RQz7Hz1uBWzUHIvc6n4J65/vED2
XbdMKKnWXEGryMTeq+Y/9RxKgpCwa69HsS2vfT/gCCUPLkoRnk7YuYPa+reh
5PHvB/zpbohIMKC49/NDyQtPfPRO4fon3n2780MGlLQ5qZOqVoH98kykSx6U
tFsvxs1eh12VUxAoCyUd3ZPfjOP8JokE3zV8CyU9zlaIh0Zgxw8lJIVDyXAq
paMyHhHPM77abnaHkoWRYQ7qBYhIjpyt8TgBJUkmu3vyKPbr5z9+yEDJ4pQC
k4w27M+MuR2mULL8fZc2fRZ70TPttjWUrP2ptGV2LyJSLqyMOI6/r2d7S/+7
QkS8lLXxd9kBGbyfWFby+Hc99cGanvIcyBDQe9Peh/OfmqJ2j1MKGUJtmifj
fmAXalMawyFDfPiqBh8Xds9Mo2g8ZGwSLpUawvMkbVdU7LlLkKFibt35sgS7
UYwhGAAZ52YyjaXKEJH+d9RCmQkZFi63XkXj85OxKtjq633IsPyrw8V3EFuU
r+S6HWTYDP1+PsvBVi70pfRChtMPmZGGUmznzxKzSpDhW57k7VOMPdbsWpgA
GZkR0akDHxHxqmH+gyuOJ1vg4rzZPPbPCxZnf0DGm+Bdus14PrwaybMaeAIZ
+b5Nf0vw92euo8JfS5BR4i6s/uQDtlr95bEpyGg2D67TwfnOjNcZeNkHGYvb
vaez3yLi9Rlj85bNkMn1wza21BrbSqb4Mh9krozQP9BExb7yU2BVI2TyzEgF
Tbtj3xfeuu46ZApzkBTA+Xmd661xUgEyZY0XtKo+IyJr/Rv+yOOQaRBwI6av
AzvjfEX1N8g0VDcDM6HYb9vf7CqDzLNjh3vW4/VZ5EFFWhRkmp0Wk5dLwq7r
drEzgUxbqdxSZxvsfxcS1upC5u38P0PceH5lq91XivoAmWl99upMPI+yizYl
2klB5iv7jSu41mCXaaqFu0Pm67EGTteyq85nb+mDzNxZVe3kZbe5OJ9KgkyS
l2Iit+x/dPefopDZoPDSYx+uX85W+9V370Lm3K3ugnPc2CHpl7UNIHORO/oW
XHZk2HhwJGRx3dfev2nZMQ5JEZ2QtTrsTUXvspM3FE8DyOJP8m69vAK78Jet
yRPIkq7cOOWJ+yunf1H1RgJkaYuYKsfi++eNquZl0ReQpVsjq013xIaPgrKV
IOuE/+T56J/Yh74d3U1A1uml4NAI3D9v9HR7W2iQZTH4ccDfH9uq3m9FL2Td
KGbGOuP7+k2QBfIYh6wkhx+rtbdgN81585ZDVrJMBr3yKXab0YqvuZD18per
whFe7I7slmefICvzrIDZfjzP3vTp5FGGIavgICzYhfvjzfSuZ3+0IauGluAg
jeuVS+34E4f3m/10vnHWEPtUlc3HV9jfDjEaDZbdvytSBfvPZutMfN5yT69a
vbMOsuZov6fPay5b9WmJN/ZlZ/rn3ctO+6iKPS96/3wM7rdcQ7/SxzieBcv0
sT2t2Gc1ep0TsV1D9gg3LtssRIeFHex4Z7h22R5KwynYecrCCRXLzvEMoELW
IqVo17p8bCOplQpq2Nm1t/4XveyFeh0t7PKc8jd43uYaS1gX/sRuf8gfHLZs
pZl/XpC1tPpswn6/ZdtLVohgn+nipLouu+1krw62Qxnv3evLnuxIacP2Tjlt
iuuXayJoo4jzt5Rh1y9os2zNG+m/sIt1d/y2WrblQlQPdvMOtzLzZd/21s3A
/i1U8sxk2U9W1x7BXppc73Z22e/8JD/9H02XaX8=
          "]]}}, {{}, 
        {GrayLevel[0], PointSize[Large], 
         PointBox[{{7., -3.428031736413475}, {7.5, -7.117937915500087}, {
          8., -10.852269845844909`}}]}, {}}},
      AspectRatio->1,
      Axes->True,
      AxesOrigin->{0, 0},
      BaseStyle->{
       FontFamily -> "Helvetica", FontWeight -> "Bold", FontSize -> 14},
      Frame->True,
      FrameLabel->{
        FormBox["\"Time (arbitrary units)\"", TraditionalForm], 
        FormBox[
        "\"\\!\\(\\*SubscriptBox[\\(Log\\), \\(10\\)]\\)S(t)\"", 
         TraditionalForm]},
      FrameTicks->{{{0, 
          FormBox["0", TraditionalForm]}, {2, 
          FormBox["2", TraditionalForm]}, {4, 
          FormBox["4", TraditionalForm]}, {6, 
          FormBox["6", TraditionalForm]}, {8, 
          FormBox["8", TraditionalForm]}, {10, 
          FormBox["10", TraditionalForm]}}, {{-12, 
          FormBox[
           RowBox[{"-", "12"}], TraditionalForm]}, {-10, 
          FormBox[
           RowBox[{"-", "10"}], TraditionalForm]}, {-8, 
          FormBox[
           RowBox[{"-", "8"}], TraditionalForm]}, {-6, 
          FormBox[
           RowBox[{"-", "6"}], TraditionalForm]}, {-4, 
          FormBox[
           RowBox[{"-", "4"}], TraditionalForm]}, {-2, 
          FormBox[
           RowBox[{"-", "2"}], TraditionalForm]}, {0, 
          FormBox["0", TraditionalForm]}}, {}, {}},
      ImageSize->Scaled[0.6],
      PlotLabel->FormBox["\"Survival Curves\\nn<1\"", TraditionalForm],
      PlotRange->{{0., 10.}, {-12., 0.5}},
      PlotRangeClipping->True,
      PlotRangePadding->{Automatic, Automatic}], {168., -547.2}, 
     ImageScaled[{0.5, 0.5}], {320., 342.}]}, {}},
  ContentSelectable->True,
  ImageSize->{348., 741.6},
  PlotRangePadding->{6, 5}]], "Output",
 CellChangeTimes->{3.3949624439747677`*^9, 3.3949625032368507`*^9, 
  3.39496474036182*^9, 3.394964796881201*^9, 3.394966950556263*^9}]
}, Open  ]]
}, Open  ]]
},
CellGrouping->Manual,
WindowSize->{638, 727},
WindowMargins->{{1, Automatic}, {Automatic, 0}},
PrintingCopies->1,
PrintingPageRange->{5, 5},
FrontEndVersion->"6.0 for Mac OS X PowerPC (32-bit) (April 20, 2007)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[590, 23, 172, 4, 28, "Text"],
Cell[765, 29, 215, 6, 30, "Text"],
Cell[983, 37, 545, 11, 45, "Text"]
}, Open  ]],
Cell[CellGroupData[{
Cell[1565, 53, 341, 8, 45, "Text"],
Cell[1909, 63, 324, 8, 29, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[2270, 76, 221, 6, 28, "Text"],
Cell[2494, 84, 170, 5, 29, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[2701, 94, 142, 3, 28, "Text"],
Cell[2846, 99, 877, 27, 67, "Input"],
Cell[3726, 128, 874, 27, 67, "Input"],
Cell[4603, 157, 882, 27, 67, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[5522, 189, 512, 10, 45, "Text"],
Cell[CellGroupData[{
Cell[6059, 203, 410, 8, 48, "Input"],
Cell[6472, 213, 928, 16, 22, "Print"]
}, Open  ]],
Cell[CellGroupData[{
Cell[7437, 234, 459, 9, 48, "Input"],
Cell[7899, 245, 924, 16, 22, "Print"]
}, Open  ]],
Cell[CellGroupData[{
Cell[8860, 266, 461, 9, 48, "Input"],
Cell[9324, 277, 927, 15, 22, "Print"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[10300, 298, 113, 2, 28, "Text"],
Cell[10416, 302, 2897, 84, 215, "Input"],
Cell[13316, 388, 2896, 84, 215, "Input"],
Cell[16215, 474, 2903, 84, 215, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[19155, 563, 626, 12, 62, "Text"],
Cell[CellGroupData[{
Cell[19806, 579, 470, 11, 48, "Input"],
Cell[20279, 592, 743, 13, 22, "Print"]
}, Open  ]],
Cell[CellGroupData[{
Cell[21059, 610, 375, 9, 48, "Input"],
Cell[21437, 621, 742, 13, 22, "Print"]
}, Open  ]],
Cell[CellGroupData[{
Cell[22216, 639, 375, 9, 48, "Input"],
Cell[22594, 650, 736, 13, 22, "Print"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[23379, 669, 473, 10, 79, "Text"],
Cell[23855, 681, 145, 4, 29, "Input"],
Cell[24003, 687, 116, 3, 29, "Input"],
Cell[CellGroupData[{
Cell[24144, 694, 2233, 57, 228, "Input"],
Cell[26380, 753, 1241, 21, 22, "Print"],
Cell[27624, 776, 1246, 21, 22, "Print"],
Cell[28873, 799, 1244, 21, 22, "Print"],
Cell[30120, 822, 1244, 21, 22, "Print"],
Cell[31367, 845, 1279, 22, 22, "Print"],
Cell[32649, 869, 1279, 22, 22, "Print"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[33977, 897, 302, 5, 28, "Text"],
Cell[34282, 904, 2579, 58, 174, "Input"],
Cell[36864, 964, 2552, 58, 192, "Input"],
Cell[39419, 1024, 2581, 58, 192, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[42037, 1087, 513, 10, 45, "Text"],
Cell[42553, 1099, 3497, 75, 228, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[46087, 1179, 383, 8, 45, "Text"],
Cell[46473, 1189, 668, 13, 48, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[47178, 1207, 313, 5, 28, "Text"],
Cell[47494, 1214, 3253, 70, 212, "Input"],
Cell[50750, 1286, 3181, 69, 212, "Input"],
Cell[53934, 1357, 3267, 70, 212, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[57238, 1432, 770, 14, 28, "Text"],
Cell[58011, 1448, 4171, 87, 230, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[62219, 1540, 604, 11, 45, "Text"],
Cell[62826, 1553, 812, 15, 48, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[63675, 1573, 508, 10, 28, "Text"],
Cell[CellGroupData[{
Cell[64208, 1587, 433, 9, 29, "Input"],
Cell[64644, 1598, 28942, 508, 757, "Output"]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)