(* 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[     74619,       1741]
NotebookOptionsPosition[     71584,       1628]
NotebookOutlinePosition[     72000,       1646]
CellTagsIndexPosition[     71957,       1643]
WindowFrame->Normal
ContainsDynamic->False*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["APPENDIX A", "Text",
 CellChangeTimes->{{3.3931681445083933`*^9, 3.393168149137631*^9}},
 TextAlignment->Left,
 FontSize->16,
 FontWeight->"Bold"],

Cell["Simultaneous determination of k, n and Tc.", "Text",
 CellChangeTimes->{{3.393253902157148*^9, 3.3932539041373873`*^9}, {
   3.393253955422089*^9, 3.393253988587735*^9}, 3.393256321009637*^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.393251160115095*^9, 3.39332011478198*^9, {3.393339984565557*^9, 
   3.393339985027905*^9}, {3.393341267796788*^9, 3.3933412735904903`*^9}, 
   3.393768131682469*^9, 3.393844075005616*^9, {3.394190265187067*^9, 
   3.394190265395576*^9}, {3.3949606069721947`*^9, 3.394960608655511*^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.393253624737135*^9, 3.3932536441172047`*^9}, {
  3.393325187213808*^9, 3.393325205517309*^9}, {3.393328605242693*^9, 
  3.393328606694022*^9}},
 TextAlignment->Left,
 FontSize->14],

Cell[BoxData[
 RowBox[{
  RowBox[{"k", "=", "0.05"}], ";", 
  RowBox[{"n", "=", "0.9"}], ";", 
  RowBox[{"Tc", "=", "80."}], ";"}]], "Input",
 CellChangeTimes->{3.393241936908412*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
Assign the three final time survival values at the original values of the \
three parameters.\
\>", "Text",
 CellChangeTimes->{{3.393327709047627*^9, 3.393327709592849*^9}, {
   3.393328629989909*^9, 3.3933286575575323`*^9}, 3.393341347327615*^9},
 TextAlignment->Left,
 FontSize->14],

Cell[BoxData[{
 RowBox[{
  RowBox[{"lstCycle1", "=", 
   RowBox[{"-", "3.56086"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"lstCycle2", "=", 
   RowBox[{"-", "6.34801"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"lstCycle3", "=", 
   RowBox[{"-", "9.73642"}]}], ";"}]}], "Input"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Assign the final time value, tCycle.", "Text",
 CellChangeTimes->{{3.393341325722166*^9, 3.3933413350988073`*^9}},
 TextAlignment->Left,
 FontSize->14],

Cell[BoxData[
 RowBox[{
  RowBox[{"tCycle", "=", "8."}], ";"}]], "Input",
 CellChangeTimes->{3.393238132173972*^9},
 FontSize->14]
}, Open  ]],

Cell[CellGroupData[{

Cell["Define the three temperature profile functions.", "Text",
 CellChangeTimes->{{3.393168351042939*^9, 3.393168364232419*^9}, 
   3.3932537841910152`*^9, 3.393322299976026*^9, {3.3933223521656*^9, 
   3.393322353083847*^9}},
 TextAlignment->Left,
 FontSize->14],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"T1", "[", "t_", "]"}], ":=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"25.", "+", 
      RowBox[{"20.", "*", 
       RowBox[{"t", "^", ".70"}]}]}], ")"}], "/", 
    RowBox[{"(", 
     RowBox[{"1", "+", 
      RowBox[{"Exp", "[", 
       RowBox[{"1", "*", 
        RowBox[{"(", 
         RowBox[{"t", "-", "8.5"}], ")"}]}], "]"}]}], ")"}]}]}], 
  ";"}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"T2", "[", "t_", "]"}], ":=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"35.", "+", 
      RowBox[{"25.", "*", 
       RowBox[{"t", "^", ".65"}]}]}], ")"}], "/", 
    RowBox[{"(", 
     RowBox[{"1", "+", 
      RowBox[{"Exp", "[", 
       RowBox[{"1", "*", 
        RowBox[{"(", 
         RowBox[{"t", "-", "8.5"}], ")"}]}], "]"}]}], ")"}]}]}], 
  ";"}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"T3", "[", "t_", "]"}], ":=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"45.", "+", 
      RowBox[{"30.", "*", 
       RowBox[{"t", "^", ".60"}]}]}], ")"}], "/", 
    RowBox[{"(", 
     RowBox[{"1", "+", 
      RowBox[{"Exp", "[", 
       RowBox[{"1", "*", 
        RowBox[{"(", 
         RowBox[{"t", "-", "8.5"}], ")"}]}], "]"}]}], ")"}]}]}], 
  ";"}]], "Input"]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
Compute the three temperatures at time tCycle 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},
 TextAlignment->Left,
 FontSize->14],

Cell[CellGroupData[{

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

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"T1[tCycle]=\"\>", "\[InvisibleSpace]", "68.93231469526532`"}],
  SequenceForm["T1[tCycle]=", 68.93231469526532],
  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.393327837129407*^9, 3.3933286890721703`*^9, 
  3.393339018778117*^9, 3.393341301635853*^9, 3.393341787060822*^9, 
  3.393768137702388*^9, 3.393768238721755*^9, 3.393844144482864*^9, 
  3.3941906364993677`*^9, 3.39419464112774*^9, 3.394195024128025*^9, 
  3.394203324320195*^9, 3.3949596669301558`*^9, 3.3949600952184153`*^9, 
  3.39496020975946*^9, 3.394960241518434*^9, 3.394960613826387*^9, 
  3.3949614851621857`*^9, 3.3949616906047497`*^9, 3.394963491262267*^9, 
  3.394963878605187*^9, 3.39496487173556*^9, 3.394967173245634*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"tCycleT2", "=", 
   RowBox[{"T2", "[", "tCycle", "]"}]}], ";", 
  RowBox[{"Print", "[", 
   RowBox[{"\"\<T2[tCycle]=\>\"", ",", "tCycleT2"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.39332215124137*^9, 3.3933221528608713`*^9}, {
   3.393325557500794*^9, 3.393325579892871*^9}, 3.393325619031274*^9, 
   3.39332609024548*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"T2[tCycle]=\"\>", "\[InvisibleSpace]", "81.91168472067831`"}],
  SequenceForm["T2[tCycle]=", 81.91168472067831],
  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.393327837210575*^9, 3.393328689155407*^9, 
  3.393339018860034*^9, 3.3933413017423553`*^9, 3.3933417871582623`*^9, 
  3.393768137807526*^9, 3.3937682388145523`*^9, 3.393844144584133*^9, 
  3.394190636580682*^9, 3.394194641235832*^9, 3.394195024225338*^9, 
  3.39420332441868*^9, 3.394959667009986*^9, 3.39496009528869*^9, 
  3.394960209849423*^9, 3.3949602415982933`*^9, 3.394960613907619*^9, 
  3.394961485227672*^9, 3.394961690674549*^9, 3.394963491342585*^9, 
  3.394963878662035*^9, 3.394964871825856*^9, 3.394967173331924*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"tCycleT3", "=", 
   RowBox[{"T3", "[", "tCycle", "]"}]}], ";", 
  RowBox[{"Print", "[", 
   RowBox[{"\"\<T3[tCycle]=\>\"", ",", "tCycleT3"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.393322170793523*^9, 3.393322172397147*^9}, {
   3.393325566670063*^9, 3.393325582959031*^9}, 3.393325626854475*^9, 
   3.3933261158418827`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"T3[tCycle]=\"\>", "\[InvisibleSpace]", "93.03654847288784`"}],
  SequenceForm["T3[tCycle]=", 93.03654847288784],
  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.3933278373035583`*^9, 3.3933286892366247`*^9, 
  3.3933390189409037`*^9, 3.39334130184315*^9, 3.3933417872473917`*^9, 
  3.393768137893032*^9, 3.393768238893093*^9, 3.393844144690846*^9, 
  3.394190636662025*^9, 3.394194641302989*^9, 3.394195024319271*^9, 
  3.394203324498796*^9, 3.3949596670867367`*^9, 3.394960095355447*^9, 
  3.39496020992784*^9, 3.394960241677791*^9, 3.394960613985651*^9, 
  3.3949614852860394`*^9, 3.394961690743464*^9, 3.394963491425802*^9, 
  3.394963878714937*^9, 3.394964871894867*^9, 3.394967173411784*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Define the three survival functions as function modules.", "Text",
 CellChangeTimes->{{3.3933277119529543`*^9, 3.393327712248871*^9}},
 TextAlignment->Left,
 FontSize->14],

Cell[BoxData[
 RowBox[{
  RowBox[{"LS1", "[", 
   RowBox[{
    RowBox[{"k_", "?", "NumericQ"}], ",", 
    RowBox[{"n_", "?", "NumericQ"}], ",", 
    RowBox[{"Tc_", "?", "NumericQ"}], ",", 
    RowBox[{"tc_", "?", "NumericQ"}]}], "]"}], ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"b1", ",", "sfn1", ",", "t", ",", "T1", ",", "y1"}], "}"}], ",", 
    "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{"T1", "[", "t_", "]"}], ":=", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"25.", "+", 
         RowBox[{"20.", "*", 
          RowBox[{"t", "^", ".70"}]}]}], ")"}], "/", 
       RowBox[{"(", 
        RowBox[{"1", "+", 
         RowBox[{"Exp", "[", 
          RowBox[{"1", "*", 
           RowBox[{"(", 
            RowBox[{"t", "-", "8.5"}], ")"}]}], "]"}]}], ")"}]}]}], ";", 
     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[{"tCycle", "+", "1."}]}], "}"}]}], "]"}], "[", 
       RowBox[{"[", 
        RowBox[{"1", ",", "1", ",", "2"}], "]"}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"sfn1", "[", "tc", "]"}]}]}], "]"}]}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"LS2", "[", 
   RowBox[{
    RowBox[{"k_", "?", "NumericQ"}], ",", 
    RowBox[{"n_", "?", "NumericQ"}], ",", 
    RowBox[{"Tc_", "?", "NumericQ"}], ",", 
    RowBox[{"tc_", "?", "NumericQ"}]}], "]"}], ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"b2", ",", "sfn2", ",", "t", ",", "T2", ",", "y2"}], "}"}], ",", 
    "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{"T2", "[", "t_", "]"}], ":=", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"35.", "+", 
         RowBox[{"25.", "*", 
          RowBox[{"t", "^", ".65"}]}]}], ")"}], "/", 
       RowBox[{"(", 
        RowBox[{"1", "+", 
         RowBox[{"Exp", "[", 
          RowBox[{"1", "*", 
           RowBox[{"(", 
            RowBox[{"t", "-", "8.5"}], ")"}]}], "]"}]}], ")"}]}]}], ";", 
     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[{"tCycle", "+", "1."}]}], "}"}]}], "]"}], "[", 
       RowBox[{"[", 
        RowBox[{"1", ",", "1", ",", "2"}], "]"}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"sfn2", "[", "tc", "]"}]}]}], "]"}]}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"LS3", "[", 
   RowBox[{
    RowBox[{"k_", "?", "NumericQ"}], ",", 
    RowBox[{"n_", "?", "NumericQ"}], ",", 
    RowBox[{"Tc_", "?", "NumericQ"}], ",", 
    RowBox[{"tc_", "?", "NumericQ"}]}], "]"}], ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"b3", ",", "sfn3", ",", "t", ",", "T3", ",", "y3"}], "}"}], ",", 
    "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{"T3", "[", "t_", "]"}], ":=", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"45.", "+", 
         RowBox[{"30.", "*", 
          RowBox[{"t", "^", ".60"}]}]}], ")"}], "/", 
       RowBox[{"(", 
        RowBox[{"1", "+", 
         RowBox[{"Exp", "[", 
          RowBox[{"1", "*", 
           RowBox[{"(", 
            RowBox[{"t", "-", "8.5"}], ")"}]}], "]"}]}], ")"}]}]}], ";", 
     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[{"tCycle", "+", "1."}]}], "}"}]}], "]"}], "[", 
       RowBox[{"[", 
        RowBox[{"1", ",", "1", ",", "2"}], "]"}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"sfn3", "[", "tc", "]"}]}]}], "]"}]}]], "Input"]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
Compute the three survival values at time tCycle using the original values \
for the parameters k, n and Tc.\
\>", "Text",
 CellChangeTimes->{{3.393322237848249*^9, 3.3933222418990507`*^9}, {
  3.393325056163465*^9, 3.3933250652433043`*^9}, {3.393327714849176*^9, 
  3.393327715408972*^9}, {3.3933286670715427`*^9, 3.39332866829493*^9}},
 TextAlignment->Left,
 FontSize->14],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Print", "[", 
  RowBox[{"\"\<LS1[k,n,Tc,tCycle]=\>\"", ",", 
   RowBox[{"LS1", "[", 
    RowBox[{"k", ",", "n", ",", "Tc", ",", "tCycle"}], "]"}]}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.393319966469637*^9, 3.393319998943075*^9}, {
   3.393320064182743*^9, 3.393320064479691*^9}, {3.393322124789736*^9, 
   3.393322131563163*^9}, 3.3933256619833193`*^9, {3.393328440806755*^9, 
   3.3933284624275837`*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"LS1[k,n,Tc,tCycle]=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "3.560857877297127`"}]}],
  SequenceForm["LS1[k,n,Tc,tCycle]=", -3.560857877297127],
  Editable->False]], "Print",
 CellChangeTimes->{3.393320121581217*^9, 3.3933221852294617`*^9, 
  3.393323123591468*^9, 3.3933252890246973`*^9, 3.3933254396472607`*^9, 
  3.393325779579335*^9, 3.3933266114096413`*^9, 3.39332706028478*^9, 
  3.3933272019418087`*^9, 3.393327837503456*^9, 3.393328689457782*^9, 
  3.393339019298667*^9, 3.393341302034075*^9, 3.393341787434205*^9, 
  3.3937681382295427`*^9, 3.393768239077735*^9, 3.39384414491912*^9, 
  3.3941906369021053`*^9, 3.3941946415387163`*^9, 3.394195024627583*^9, 
  3.394203324677058*^9, 3.39495966741*^9, 3.394960095508278*^9, 
  3.3949602100975933`*^9, 3.3949602418481903`*^9, 3.3949606141582813`*^9, 
  3.3949614854364567`*^9, 3.3949616908886423`*^9, 3.394963491602689*^9, 
  3.394963878847652*^9, 3.394964872037624*^9, 3.3949671735885477`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Print", "[", 
  RowBox[{"\"\<LS2[k,n,Tc,tCycle]=\>\"", ",", 
   RowBox[{"LS2", "[", 
    RowBox[{"k", ",", "n", ",", "Tc", ",", "tCycle"}], "]"}]}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.39332007300701*^9, 3.3933200828250637`*^9}, 
   3.393322161484642*^9, 3.393325693290319*^9, {3.393328469035684*^9, 
   3.393328472514614*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"LS2[k,n,Tc,tCycle]=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "6.348013272924269`"}]}],
  SequenceForm["LS2[k,n,Tc,tCycle]=", -6.348013272924269],
  Editable->False]], "Print",
 CellChangeTimes->{3.393320121694007*^9, 3.393322185357984*^9, 
  3.393323123750484*^9, 3.393325289185557*^9, 3.393325439912942*^9, 
  3.393325779705103*^9, 3.3933266115143642`*^9, 3.393327060391836*^9, 
  3.393327202047213*^9, 3.393327837595304*^9, 3.3933286895429773`*^9, 
  3.393339019404811*^9, 3.393341302127349*^9, 3.393341787541155*^9, 
  3.393768138319903*^9, 3.393768239170231*^9, 3.3938441450409327`*^9, 
  3.3941906370296392`*^9, 3.3941946416214952`*^9, 3.394195024712863*^9, 
  3.394203324782199*^9, 3.394959667502877*^9, 3.3949600955835457`*^9, 
  3.394960210178719*^9, 3.394960241929119*^9, 3.394960614237028*^9, 
  3.394961485499854*^9, 3.394961690967836*^9, 3.3949634916920757`*^9, 
  3.39496387891539*^9, 3.394964872116544*^9, 3.394967173667007*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Print", "[", 
  RowBox[{"\"\<LS3[k,n,Tc,tCycle]=\>\"", ",", 
   RowBox[{"LS3", "[", 
    RowBox[{"k", ",", "n", ",", "Tc", ",", "tCycle"}], "]"}]}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.393320087098538*^9, 3.3933200985041313`*^9}, 
   3.393322181469367*^9, 3.393325697880682*^9, {3.393328477223815*^9, 
   3.393328481340178*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"LS3[k,n,Tc,tCycle]=\"\>", "\[InvisibleSpace]", 
   RowBox[{"-", "9.73642130875272`"}]}],
  SequenceForm["LS3[k,n,Tc,tCycle]=", -9.73642130875272],
  Editable->False]], "Print",
 CellChangeTimes->{3.3933201217909107`*^9, 3.393322185505349*^9, 
  3.393323123897531*^9, 3.393325289315895*^9, 3.393325440102467*^9, 
  3.3933257798152037`*^9, 3.393326611610318*^9, 3.393327060510661*^9, 
  3.393327202142065*^9, 3.393327837688336*^9, 3.3933286896359453`*^9, 
  3.3933390196553507`*^9, 3.393341302222069*^9, 3.393341787845975*^9, 
  3.3937681385185337`*^9, 3.39376823926608*^9, 3.3938441453338413`*^9, 
  3.394190637141371*^9, 3.394194641745226*^9, 3.39419502480593*^9, 
  3.394203324884741*^9, 3.3949596675799513`*^9, 3.3949600956660433`*^9, 
  3.394960210271525*^9, 3.39496024200951*^9, 3.394960614317299*^9, 
  3.3949614855660067`*^9, 3.3949616910374537`*^9, 3.3949634917748537`*^9, 
  3.394963878981738*^9, 3.394964872186202*^9, 3.394967173746498*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
Find values of the three parameters k, n and Tc that simultaneously solve the \
three survival equations at time tCycle.  Give two initial guesses for each \
parameter in order to avoid the use of derivatives.  NOTICE THAT THE \
SOLUTIONS ARE ALMOST IDENTICAL TO THE ORIGINAL PARAMETER VALUES k=0.05, n=0.9 \
and Tc=80. (The slight discrepancy is due to the rounding of lstCycle1, \
lstCycle2 and lstCycle3.)\
\>", "Text",
 CellChangeTimes->{{3.3932382688789*^9, 3.393238269494989*^9}, {
   3.393243093902789*^9, 3.393243125209425*^9}, 3.393253598496396*^9, 
   3.393253749723592*^9, {3.3933252176594467`*^9, 3.393325233383534*^9}, {
   3.3933262342727222`*^9, 3.3933262354459248`*^9}, {3.3933277198356867`*^9, 
   3.3933277210012293`*^9}, {3.393328674593845*^9, 3.393328676071*^9}, 
   3.3933416822146597`*^9, {3.393844063195136*^9, 3.393844065700686*^9}, {
   3.393844337472574*^9, 3.393844340865069*^9}},
 TextAlignment->Left,
 FontSize->14],

Cell[BoxData[
 RowBox[{"Clear", "[", 
  RowBox[{"k", ",", "n", ",", "Tc"}], "]"}]], "Input",
 CellChangeTimes->{{3.393243616755053*^9, 3.3932436170077763`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"FindRoot", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{"LS1", "[", 
       RowBox[{"k", ",", "n", ",", "Tc", ",", "tCycle"}], "]"}], "==", 
      "lstCycle1"}], ",", 
     RowBox[{
      RowBox[{"LS2", "[", 
       RowBox[{"k", ",", "n", ",", "Tc", ",", "tCycle"}], "]"}], "\[Equal]", 
      "lstCycle2"}], ",", 
     RowBox[{
      RowBox[{"LS3", "[", 
       RowBox[{"k", ",", "n", ",", "Tc", ",", "tCycle"}], "]"}], "\[Equal]", 
      "lstCycle3"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"k", ",", ".04", ",", ".045"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"n", ",", ".6", ",", ".7"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"Tc", ",", "50.", ",", "55."}], "}"}]}], "}"}], ",", 
   RowBox[{"MaxIterations", "\[Rule]", "500"}], ",", 
   RowBox[{"AccuracyGoal", "\[Rule]", "6"}], ",", 
   RowBox[{"PrecisionGoal", "\[Rule]", "6"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.393241157232133*^9, 3.393241203779924*^9}, 
   3.393241235337063*^9, {3.39496020179797*^9, 3.394960202273522*^9}, 
   3.3949602349821663`*^9, {3.394960310774068*^9, 3.394960330891343*^9}, {
   3.3949604170068083`*^9, 3.394960526402029*^9}, {3.394960583267733*^9, 
   3.39496058513557*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"k", "\[Rule]", "0.04999963804016134`"}], ",", 
   RowBox[{"n", "\[Rule]", "0.9000069226611157`"}], ",", 
   RowBox[{"Tc", "\[Rule]", "80.00036498347139`"}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.3931687123151217`*^9, 3.393238834655732*^9, 3.393241119796772*^9, 
   3.393241175743135*^9, 3.39324121643988*^9, 3.393241250957178*^9, 
   3.393243950369433*^9, 3.393244009966853*^9, 3.393251086417383*^9, 
   3.393253839954664*^9, 3.393254028134224*^9, 3.3933201326291523`*^9, 
   3.393323139280815*^9, 3.3933253001919394`*^9, 3.393325450955555*^9, 
   3.39332580582519*^9, 3.393326622768568*^9, 3.393327071513707*^9, 
   3.393327213914505*^9, 3.393327849208995*^9, 3.39332870059485*^9, 
   3.39333903129743*^9, 3.3933413142447166`*^9, 3.39334179920826*^9, 
   3.3937681496552353`*^9, 3.393768249599525*^9, 3.39384415654799*^9, 
   3.394190648138094*^9, 3.394194666285651*^9, 3.394195046188807*^9, 
   3.39420333674893*^9, 3.394959678713756*^9, 3.394960106751211*^9, {
   3.3949602205629396`*^9, 3.39496024803575*^9}, 3.394960345479535*^9, 
   3.394960465748316*^9, 3.394960501958528*^9, 3.394960533554636*^9, {
   3.394960590600037*^9, 3.3949606158457203`*^9}, 3.394961487120377*^9, 
   3.394961692780471*^9, 3.3949634933330717`*^9, 3.394963880579673*^9, 
   3.394964873920156*^9, 3.3949671756677628`*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Plot the three temperature profile curves.", "Text",
 CellChangeTimes->{{3.393168351042939*^9, 3.393168364232419*^9}, 
   3.3932537841910152`*^9, 3.393322372504961*^9, {3.393323367039167*^9, 
   3.393323380899116*^9}},
 TextAlignment->Left,
 FontSize->14],

Cell[BoxData[
 RowBox[{
  RowBox[{"tempCurvePlot", "=", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"T1", "[", "t", "]"}], ",", 
       RowBox[{"T2", "[", "t", "]"}], ",", 
       RowBox[{"T3", "[", "t", "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"t", ",", "0.", ",", "8."}], "}"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0.", ",", "10."}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"0.", ",", "140."}], "}"}]}], "}"}]}], ",", 
     RowBox[{"AspectRatio", "\[Rule]", "1"}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{"Thickness", "[", "0.01", "]"}], ",", "LightGray", ",", 
          RowBox[{"Dashing", "[", "None", "]"}]}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"Thickness", "[", "0.01", "]"}], ",", "Gray", ",", 
          RowBox[{"Dashing", "[", "None", "]"}]}], "}"}], ",", 
        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[{"\"\<\>\"", ",", "\"\<Temperature (\[Degree]C)\>\""}], "}"}]}],
      ",", 
     RowBox[{"FrameTicks", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0", ",", "2", ",", "4", ",", "6", ",", "8", ",", "10"}], 
         "}"}], ",", 
        RowBox[{"{", 
         RowBox[{
         "0", ",", "20", ",", "40", ",", "60", ",", "80", ",", "100", ",", 
          "120", ",", "140"}], "}"}], ",", 
        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.393323043264168*^9, 3.393323051596662*^9}, {
   3.393323223959093*^9, 3.3933232250892553`*^9}, 3.393323513728169*^9, {
   3.393323758811049*^9, 3.393323777192898*^9}, {3.393324056535119*^9, 
   3.393324057701312*^9}, {3.393324139126712*^9, 3.393324141726007*^9}, 
   3.3933254343125553`*^9, {3.393325824794448*^9, 3.393325852804388*^9}, {
   3.393768110411003*^9, 3.393768111385027*^9}, {3.394190300473057*^9, 
   3.394190335350337*^9}, 3.394190368730852*^9, {3.394190450076414*^9, 
   3.394190473476864*^9}, 3.394961754246566*^9, 3.394963574587223*^9, {
   3.394963642870945*^9, 3.394963643558456*^9}, {3.3949636827064857`*^9, 
   3.394963726043089*^9}, 3.3949637714028263`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
ListPlot the temperatures at time tCycle 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},
 TextAlignment->Left,
 FontSize->14],

Cell[BoxData[
 RowBox[{
  RowBox[{"tempPtsPlot", "=", 
   RowBox[{"ListPlot", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"tCycle", ",", "tCycleT1"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"tCycle", ",", "tCycleT2"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"tCycle", ",", "tCycleT3"}], "}"}]}], "}"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0.", ",", "10."}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"0.", ",", "140."}], "}"}]}], "}"}]}], ",", 
     RowBox[{"AspectRatio", "\[Rule]", "1"}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"Directive", "[", 
       RowBox[{
        RowBox[{"PointSize", "[", "Large", "]"}], ",", "Black"}], "]"}]}], 
     ",", 
     RowBox[{"PlotLabel", "\[Rule]", "\"\<Temperatures at tCycle\>\""}], ",", 
     
     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[{
         "0", ",", "20", ",", "40", ",", "60", ",", "80", ",", "100", ",", 
          "120", ",", "140"}], "}"}], ",", 
        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.393768222004148*^9, 3.394190366265501*^9, 3.39496158748779*^9, {
   3.394961738913563*^9, 3.394961759929484*^9}, 3.394963576300029*^9, {
   3.394963645815774*^9, 3.394963647553526*^9}, {3.39496368655053*^9, 
   3.3949637343070593`*^9}, 3.394963777353129*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
Show the three temperature profile curves and the temperature points at time \
tCycle 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}},
 TextAlignment->Left,
 FontSize->14],

Cell[BoxData[
 RowBox[{
  RowBox[{"tempProfilePlot", "=", 
   RowBox[{"Show", "[", 
    RowBox[{"tempCurvePlot", ",", "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.394961591747574*^9, {3.394963488390458*^9, 3.394963513027947*^9}, 
   3.394963583640977*^9, 3.3949637799479723`*^9, 3.394963819192431*^9, {
   3.3949638694659758`*^9, 3.394963896715249*^9}, 3.394963933351191*^9, {
   3.394963975960422*^9, 3.394963998985455*^9}, 3.3949640868480263`*^9, 
   3.394964134425975*^9, {3.39496417563127*^9, 3.394964175765367*^9}, {
   3.394964211335552*^9, 3.394964249307542*^9}, {3.39496483079762*^9, 
   3.3949648309829473`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Plot the three survival curves.", "Text",
 CellChangeTimes->{{3.393168351042939*^9, 3.393168364232419*^9}, {
   3.393168404854501*^9, 3.393168414243094*^9}, 3.393326428654821*^9},
 TextAlignment->Left,
 FontSize->14],

Cell[BoxData[
 RowBox[{
  RowBox[{"lsCurvePlot", "=", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"LS1", "[", 
        RowBox[{".05", ",", ".9", ",", "80.", ",", "t"}], "]"}], ",", 
       RowBox[{"LS2", "[", 
        RowBox[{".05", ",", ".9", ",", "80.", ",", "t"}], "]"}], ",", 
       RowBox[{"LS3", "[", 
        RowBox[{".05", ",", ".9", ",", "80.", ",", "t"}], "]"}]}], "}"}], ",",
      
     RowBox[{"{", 
      RowBox[{"t", ",", "0.001", ",", "tCycle"}], "}"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0.", ",", "10."}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"-", "10."}], ",", "0.5"}], "}"}]}], "}"}]}], ",", 
     RowBox[{"AspectRatio", "\[Rule]", "1"}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{"Thickness", "[", "0.01", "]"}], ",", "LightGray", ",", 
          RowBox[{"Dashing", "[", "None", "]"}]}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"Thickness", "[", "0.01", "]"}], ",", "Gray", ",", 
          RowBox[{"Dashing", "[", "None", "]"}]}], "}"}], ",", 
        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[{"-", "10"}], ",", 
          RowBox[{"-", "8"}], ",", 
          RowBox[{"-", "6"}], ",", 
          RowBox[{"-", "4"}], ",", 
          RowBox[{"-", "2"}], ",", "0"}], "}"}], ",", 
        RowBox[{"{", "}"}], ",", 
        RowBox[{"{", "}"}]}], "}"}]}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.393238418107068*^9, 3.3932384316360903`*^9}, {
   3.3932510563551903`*^9, 3.393251064837891*^9}, {3.393320379397503*^9, 
   3.393320459627843*^9}, {3.393320519548664*^9, 3.393320529808571*^9}, {
   3.393321062039523*^9, 3.393321160769993*^9}, {3.393321206645903*^9, 
   3.393321213195038*^9}, {3.393321248567235*^9, 3.393321273553196*^9}, {
   3.393324310926723*^9, 3.393324342146593*^9}, 3.3933244574107733`*^9, {
   3.393324507246973*^9, 3.393324526172906*^9}, 3.3933245598079777`*^9, {
   3.393325881760735*^9, 3.393325899679166*^9}, {3.3933412888037357`*^9, 
   3.39334129069744*^9}, {3.393768120150341*^9, 3.393768122735518*^9}, {
   3.394190482251485*^9, 3.394190550353314*^9}, 3.39419463420961*^9, 
   3.394195014291848*^9, 3.394959644464952*^9, 3.3949599014110518`*^9, {
   3.394960038378326*^9, 3.394960072180572*^9}, 3.394963566424345*^9, {
   3.3949636343492403`*^9, 3.394963637450425*^9}, 3.394963706440468*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
ListPlot the survival values at time tCycle 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.393338863451004*^9, 3.3933388660281363`*^9}, 
   3.393341400159108*^9, 3.39334175038418*^9},
 TextAlignment->Left,
 FontSize->14],

Cell[BoxData[
 RowBox[{
  RowBox[{"lsPtsPlot", "=", 
   RowBox[{"ListPlot", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"tCycle", ",", "lstCycle1"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"tCycle", ",", "lstCycle2"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"tCycle", ",", "lstCycle3"}], "}"}]}], "}"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0.", ",", "10."}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"-", "10."}], ",", "0.5"}], "}"}]}], "}"}]}], ",", 
     RowBox[{"AspectRatio", "\[Rule]", "1"}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"Directive", "[", 
       RowBox[{
        RowBox[{"PointSize", "[", "Large", "]"}], ",", "Black"}], "]"}]}], 
     ",", 
     RowBox[{"PlotLabel", "\[Rule]", "\"\<Survival at tCycle\>\""}], ",", 
     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[{
       "\"\<\>\"", ",", "\"\<\!\(\*SubscriptBox[\(Log\), \(10\)]\)S(t)\>\""}],
        "}"}]}], ",", 
     RowBox[{"FrameTicks", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0", ",", "2", ",", "4", ",", "6", ",", "8", ",", "10"}], 
         "}"}], ",", 
        RowBox[{"{", 
         RowBox[{
          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.393768228139571*^9, 3.3941904114955177`*^9, 3.3941906287898073`*^9, 
   3.394194631612575*^9, 3.394195018966229*^9, 3.3949635698982353`*^9, {
   3.394963629567026*^9, 3.394963630359771*^9}, 3.394963701913322*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
Show the three survival curves and the survival values at time tCycle \
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}},
 TextAlignment->Left,
 FontSize->14],

Cell[BoxData[
 RowBox[{
  RowBox[{"survivalPlot", "=", 
   RowBox[{"Show", "[", 
    RowBox[{"lsCurvePlot", ",", "lsPtsPlot", ",", 
     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.3949598014125843`*^9, 3.3949598046823587`*^9}, {3.3949598370204353`*^9, 
   3.394959840801425*^9}, 3.3949600883031683`*^9, 3.394961601593752*^9, {
   3.394963453636385*^9, 3.394963461787159*^9}, {3.3949635043276463`*^9, 
   3.394963505347754*^9}, 3.3949635588109303`*^9, 3.3949640978855953`*^9, {
   3.3949641401671247`*^9, 3.394964198978652*^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.394963269309606*^9, 3.394963271125407*^9}},
 TextAlignment->Left,
 FontSize->14],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"figure3part", "=", 
  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.394964015767812*^9, 3.3949640233060083`*^9}, {
  3.3949640595347357`*^9, 3.394964094181012*^9}, {3.394967162408063*^9, 
  3.39496716636302*^9}}],

Cell[BoxData[
 GraphicsBox[{{}, {InsetBox[
     GraphicsBox[{{{}, {}, 
        {GrayLevel[0.85], Thickness[0.01], Dashing[None], 
         LineBox[CompressedData["
1:eJwV13k8VPsbB3BbESLbzJkZmrJdCRGRJecRshbDnDZLXUVRiBQqa1xlqUgI
2RKVLNmpq7K3kCXkZksoy89u7Pp9+2te79eZ8z0z5/s8n+ccK3eHnY5sbGwc
bGzskn8+dUNH37b+cjSF/Rvv53ucB/CzlP5LVH0u0DQdHbzlz8S9KEXJCamq
oHkt6ElUgSHuL+sdJC0+BJqPpQjzJhYeQnlU2NVkBpoNbZxVQdL4bfrvJxKR
0qA52sbXO8bC78j+HX/GdA20trQfPiY2iMdSYv5dPHsMtAxddtWrFuBxyTVl
VP1a0Prbj6XANow/pM8XHJBSBq0bT8L2VhThKbJHM4IHeUCr4LE4vsKNZ1Ki
9MP29IK2kHsX7xwv/sRG5i2WGgvach6+u1Us8azkKp3ngqagjfOtD3CX48/o
U/s+TZaC9gU5uwcpZDxPliEnnHsHtN8kd/zwH8WLKeH/WIilgY7tN+65N854
ceXR2vim66DjnmaTUamJl9hIcQyEHgWdYC0F33lBvDT5VYDHAj/oZM18UHUK
wivoEz6xX3xAZ/LhFqspc7xK9ojLf/cZcMC7sXdwwhWvpdz+q4gnCHS9Vd//
lBPCay8V3oo2kgfdsJrysYJuvLbh26h7aBvoJugGWV5m4XVXlXIUOKVBt5zl
xH+sBa9v+6KUuV4PusvXRaJLo/D3kfR9cXP8gHvvOvqxWB9v3ijRv9abAKAa
NTfT+BH/MjOdn37iOOgJixgJrX7A+7cb2Qb95oODYXWLAnRtfCii5tckQwb0
2//ePzCJ4aPzCbhmszcY7rgjLh3SjU+xs/cJFZPhkMGI/cTdk/iszudP+oOl
YOTRTpiWmeDzEydoj1yqwThyO7/Luhu+6JoreEB2AUxuRxW3ewjiKzK0bTnL
HGB6X3itm5WBr+VZl57+rQ5mN/VLcVFO/Ld8fbNayCcwP2985bd2KrBNPlse
v0eHwyc1TPgVhYFDuaOvf3EcjqgfDmq7mQqcdwQ+vBrLBAu6lIxpzgpwDQ1N
ZUsIgMUGlIhKasFm6TSnEudIsDzbxdfmFAvcAeeHcoAFli3FdtuCuICnhlRc
u9cLGKp0tcmqh8C7szWiNq8fGEmWCr3OosB3bVfSF5dVYPym6s7T0oH/ldMG
t95usDoxspRgMg8C5JfndmYYgVVFlb2pgwgI2n6XPCXhDdbbntyXZ2yFbYe7
+IgJBbA+ozf1U90ctq3U8kwLioB1RVVGjecdEEqWyk96wwVMbh2vLYrZIGy2
pfuV2gVgHiPiOXw1QZhVleZgkgPM9JgwySECROJGtlHHF4A56Ph8a7sbiBrE
chpZlwAhpZurl1MEoiPFqaHGAUCcWD1p9TMZxEIffzXt3QnEbd+VrD5tELtF
5GtXhwMRrhkwMhQNYhHc/yhkzyIvBS/bNoLYvQtqAp7VQEReYQ/8ayuIJe29
28ZzFoi7575eHS8AsYJ3RidUnwLxQK0gvU0XxL71lZ4LVwIizZnMSKkBsb57
Ucfd44BIJze1+iB/P3jWhPkbuW5Z9cQQiI1kCe+mtwDxeEfulfxXIDbj6jpZ
4gHEk883vyqKAGnzmsyVgWIgcoSzKv6OAdJeSmyIhhYQJQqFqda9QFJ7pE5x
6ESuWlK+KQskjR1deXcuA1Fq+drRfh5IOnLUr8MvgCjzfPPqtCOQDqmnK8Zu
B6KiIIdmPwykk9a5nTMcQFSJ9JyzmwZScFTDrucfgai/8v38X9lAChU6X9V5
HnnRn6NfHUhhD3iZHFxANPhEkzoigBSZfNj/hC4Qjb6UR+tKQIp73t7OXQDE
h8svPuYzgfS8YeDG2Vggmo+k2d/oAlI7+/Ln7bZAdDzsiddTBlLHqxJPigAQ
nSIyzzYOA+nrFQ8xkXfIEdUbHDeA1Pvrl80WWSC6ri1zOHoA6dfnzpG5KSC6
GcZp1A9A2nhUtPE+GIieXtnW3ttAZj/mnl6rBkSv2e2I/7UCmUtot0HVT+QK
ra+xukDeEpIRXmwORF9MWPmkL5BFL0STUsWAGNB56zPcDGR5TTclL7SfP9xt
41qvAllhbler+0nkthvL8Z1A3vNi2MtlKxBDaiO0x/1A3rfDtvKUJ/JitFSj
PpCB2+yQqTYQI174LwEzIBMdcvb0ZiBGTV+MYUtADro0GPVhAYhJZ8mR4z1A
vrl2voOvBzk/++X8BpD/CZsSN0f1NrkQrNYRCuSIlNXnzfeAmArQ2P7cBcgP
mkQb2uSBmA6PL2DNAfn5LiP2bnsgZgMnpnS2AvlFcZMJ1QD5dfahjHgg54P1
PZvdyIuFB4d1gFx87PT2niUg5i6Uf2flALkq1Fd7IAaIebOA1zlvgdz2PefK
z3ogFkZzsncXALnDVeVfOVQ/LKrJ0wx0P7uWyrmc7yObqkrVhQC5V7A2Zvw0
cvbP/YvVQP55oDd/cgWIxZMfyiy4gTz6/szingHk0FeneB4BeYI5pnupATl/
/nt/CpBnLrA+zaL6WOLQ+2voLpBXEwVHWXuQU/b5at8A8oZMnIoGCbnmTJXX
VsDYXor7eK8h//yq5TUH2KbGXZuXPwCxrBQTzPoBmCDroNQ66s/lIv3NCvsA
Ew5qvKBritz2sPCqBmBiW48UBqggz6iOP3AHjCpto8e2AcSKgkpuuRlg0lZe
9lwJyA9NBOJlAJNt3Uc+8+d44b6t12cA22Wx0FKD1l/54Hx0rQEwpcPe+sHK
yEvPQrz1AFP+pLH2HfX7qkBxXcw3wFRNl0oB/d5VaW6+O4mA7Te+Jve7Edk8
4FS8J2BaDVqDp1BerJ7K/KAYCNgBw9WkKnQ/Vj0FZ0vqATuo7yfgh+7vamzp
vv4qwAyqDzT0oP1ZzZx+4aUEmBFsBOkoIBe7JwyEAGauGzC3sojcqhlug44f
+RdyT6J+W+27OdY4BRhDh92pAtXX6riQd1MGYNaV1Tsocsis733jhYAd1Qzu
9rkDxBr7fNP9OMCOl+vHdM0j85vMSI0CZqPBZa6B+mFNbIA1ehswu5K6TXFv
kbdXRNIzATutFlrFQv26JttuN9UNmEPRIR8iEllJieJpAZjjXm6VYpSva2rv
L+R+BuxcQcO46HFkrSfs8T2Auey5lXm5Chl/e0v0I2AX80zs2qWRDWg5i+i4
uyIveS/K6zWjUiqnHWAeOR9aoqeRTe/6/yoDzEs+4vbMUWTz7DDjaMCuPjPX
t3yNfHi9v3AEMF85/rUCSeQjMWrZOwC7ntVUKnjrj8/ujUX14S9zx91t8s/3
L4d+rQUsMNNCrpmJbPbmvxo6YDelBAcVK5FNjHuj0X6FprckRe1APiSg5SkC
2K0d0cyJUOSDgo9uigEWnmolYDaGfMC8sG0WsKjtwg05lsgadbKnbQG7m9we
xFuKrOKTSXUDLIYWq+Uijizv1DWArh+byJz7k29rknddb8YCFk8Ry5X7hUyZ
3VbTAlgyKZ7+E82Htc1nlb48BSzlwbFuIyraz9XLbN+aAEsXxWKyApGnKu/8
LxewLKHETY4o/1Y7ZiPXqIA9vXeyqrYQuWHQuuIWYDmCNB9pDLliSwjnS8AK
+JPHBoeQE9aaJf0BK4ywyzyI+mk1rFH3kChgJbzb7dJR3q9eeV8bGwFYJXdq
y+kbyBaub8pQvdZyZJT0iaD+WFbeaXINsPoyPz9+VF8r/YnDTkGAvb943ECb
F7lOdPK+OWDNnVvbEtA8WonaJy6K6rP7mc8UA9XLijBlrkIFsB5769JAN9S/
07PTXwwB6xdR8stD93e5ueu9Gtq/oRtD/HzfkcNqiIwBwKYsLeTq0HxdmhTr
G18CbHbTrql59P+X6vTKGeyAzVdylUmiflxKuvJbLhKwFalKw4A3yAYsCXkp
oHAtyjho5qG8ijS+CyJAoT5aTcpF9bww3lGPcpsiYdXp0MODXMR8P6cClB3c
L+V5Q5B9OwiNRqDIeDiVn7uOzNl96EUMUJQNWjt2uqC85RvwejsJFMPxpwIP
jIGYGQ/vtggBinFacEdNLXLGpIV6AFDMCLvkWUA+YW5zgQcolm9EdltooPyv
ZvlafwOKTYy/8Z/5OHWblfTZFCgemswgv01ATIwY/M8zHyiX3yxPzKD6mLiP
8bp2A+XqodTjTii/JnB2p3eCQLluPbbHAl1vPFbM07oFKKFuAb2S6PePqe4f
WakAysPHOfsb0bz/aSjR4q4IlGR5y0wdtN7ID3mrhNdASSlgCb5E640EGp2j
qwHlcRWMJKL1hsteftsYAEru185Y15toPtJ0N3uqA6VagH1aBO3vQPGPqRx/
oNTGZtuEofob0P+8xfAzUBpo5vVrDkD0t9ZtwdD/+yQX/2gY5VXfWItl8xJQ
OvUVzcpRP/QKnSQbJQJlzOd4tj07eh7YM258WhqowoP59s/YgGgpiOgdjgWq
WKpb1wTKyxZRrm03eIGK2SoxlNH+fb5qIRv2CKgSHbn65WT0fKJhq1bID1S5
xhy5BpTHH3MMHSa7gaqblzU7jPKz/sT0XSo63+VaStjOCiAqP8fEWzUB1VXD
ns3RH5k/Is1aH6ju8xK+z9D8rTA28dAMAeoV1+SLyqj+yl4L5RYVATXwdCID
R/1bHNt4scIJqHGH4mj2aL7l8brp7wKgVgtH5Sei54k0zuS9G+5Ara0OS1wr
AyJV/Enuo3CgNngG/WOP5kmKmqd8Fjr/U9sVO6kwIJJPW8U3nwFqV4wd74tk
IBKeyWrpTwP1f8IKjlV1QNzbeMcdHwM0TOgj7Qfq78AyIZlWa6BR39VxG9oB
EZBdaRsyAjQJj6rZ7Awg/OOpJqo6QJNsLXx/Ec2XGx5rDS3VQFOMTry6oAeE
D6a3rkEBmr6Qcxs3qmd32etU4UyguQnx3Nr9LxC2RpLv9uwA2qWYJvIXdSDQ
e2GSxCegXRaNzr6B8uLEesWiDSfQfDBqPZpDxLH8PGpqBtCC6Yqb3FEeW2/K
tFowAlqcIiO4yAoIE9+IS/3zQKsySfTTQs9z+wLOiWjsAHGBYAVPUyEgtiwo
c7u4griQFp/guzQgeBw/8rzjAXGRmdFcTTSfN39hFW1VBXHK31m/5I4AwZm3
MdrJBHFpvZ32m8OBuc74aHrCG8S12EgmNej9Y9oKC8OTQdzRn22nNguYHdrv
hBv9QPy1b3uLYg8wE7NeXG0OAPE3nH5C46nATPjvZJwNL4i/i5K1fuYAzHj+
18fTvEC8Ps2nS3oUmPddEw2MzUG8pUGij4bWi6S/NhYtAvFh0XMTvELA9DP0
HxZLAImteWvcI4eAaR9ncJDtEkjYD0njKQXAlDD2yTwyCxIbEkLq+XZgHbVS
bxtPhe3Pe4LSQwPB6tvPWK/KYaCbCyRJ/MoFq522EqcDXwN9pWG/gBADGAfT
t4nefYjcQ+PYFI4cCqXs2sjT6wvLtcjOHsIe3UBfpdb0Dqojq7SsEbbIbha5
RVRg6FVHxZpGAn1N7Nxh5iAw4MdyJyke6Otn4iIT0PsiLpvyPAAdX/fxDYzM
Reatz3OoR46y9QocAYbu5KSruzdymaSt83HksgMSvPeBvsGbp6B1ANmo6+DF
UeSXdR97NgPjwPmNbH495Pqnb1pw5MOS9rJPkL9FFNX6IKsYtgQ0Af33JkbS
izFg6CxHjMxVIR/vveDXjBwu+o37GrLr21Oe3Mhu6hcPWCEHP7Z2AmSrY+Wv
k5BfOGsfKUSmPryZ/Ge9d+Z7Do4DQ3ujkoPNAblzj5SGNPLgf9qF/yKPC5N3
2yPXryiXqCD/XuCjxyPnUL5vIf8fnZQrFQ==
          "]]}, 
        {GrayLevel[0.5], Thickness[0.01], Dashing[None], 
         LineBox[CompressedData["
1:eJwV13k4VG8UB3A7Ibu5d8YyEaIQFf0suQfJLjKWuYoKJWWpiJRsTWQL2YpS
SMlatqKUPUmkEFGypY02yRK/t7/m+Tz3ufe9c973fM/MLr/9cp5sbGwcbGzs
8v8+DVgfH7+Y8rQAwzBt/yOzXYQH9Z0/zZgLDB8dXv4+EkcEUCuyM3M2g+E4
30KqQwxxRikoQkF6HIw4TkqnPOwgzlKv3O3vtAQjmdNZ25RciPP0lRsy8Qpg
pK1V2nfLi0hU2pfhbrEERtbE+dM1uUQqNeXhnIcTGAXl9svmjBHp2U01NONm
MIo7pDphX0xcov8q37ZWA4yu+kyf6s8lrio55kaO8oFRY6qzJPgR+dQE4+iN
w2DMjWnGtXcRN1wUH+M5qWAs0dTEmpklCrLr9W8LW4CxguwZcydxopA+o/Vs
uhqMDX9zv1YfI0qV7JTFShLBOKR0UWvsPlFJjT23U/IaGI8JWO5xJIjKWsfm
jM5TYDxzcF6gJYeoclnLMcJyBOPFPIMXP94R1dl1YUdnBWG7eG3ThVdpxH36
l+DUV8Gw3SjhnE6+MFGvZOM9eNEOtmclMPoDColm6vl1FXwRYEJYqvk45hDN
/ndjkk3Xg4l5r0XQ1SGiue3NRz9WD5gwdlY52WYRLSfUi1Q5FcDkEOVlvkcp
0drzSj3/byuYJM9RF7YJEO3xdK30n4Jg8t5AvTq1n3i+XGUcMpwJO4L82FcS
NxGvvn8ru850BtPYGLUHweNELy2Xn/FmBkxT1q35Snckeo3tPflcz4Hppccz
4tPTRG9qNc3PvQJMb03YiRnjRN/WU6xtvkJg2vojqY8MJ/pDeciByBYw41AJ
9ihQJwb5pHlEizeBWaD9xrXLBcQ7WdPdESsCYG57fHdFug4xHtc0NW2nCBZD
6VRl7Vni469MQud5EFh5xslkvN1HzLCzvxWtxMC69oWtf3Q08UO/65nxaDXs
xGRorbkXiV9fmFJXvBvBVq0863dLPzHnUyK8TWkWbOsZnAqTgcSCopRI0TwH
2DmwnT7uTyWWSu2r965og91QxvCSz1ViZX3r8y1nn8GuPeZJHG4KwDZdOP85
iQ67xvOjYtVigUOj9+27uc9g72b8dDbnFnAmCj2t+5QP9iMn5/zvFAHX+PjM
TRkhYLh0HQzHyoBH4dqBqkPxwOgYlF0WqwLeMK/xIvgNDkYGmudSmoCviVLZ
vCkAHEpGSA69KeCXexHXXPoOHOlel0SsZEEgRCXrlfciOMbdPzFRUweCdQeW
eQ03gOPXyeVvafdBCLtzUC7XFJwcaFXrL70G4d3v5d1kgsDprtHHXb9GQMS6
X8Dhiyo483u25DIOgshCM983YXFwPqCi21mgCaLZa8uyHnGBc81gv83LZBCz
XDVQt+UwMAVuuh99LwZiv+uv7TcvAqZzKt1ymQDx9EkR2udZYF6e7/P5ngAS
21M5Te2rgDnk+YRcfgUSk5U5LLMwIMVOX9q6XQMkPmrs9nLtAVLcSKrAnQES
n4uploEKQErwKueLbgKJb3kXRa8/BVLywvtDxABILKScu7IoASSe1fahmAck
hfx9qktuASl7RfuBgjlIaqvqTYl1A7nehEeQuxgkWXmvLYblgASpPVER/iAZ
41Cm1xiL3NXnre0FknG851Rv/gDSMKrC8YYRSCYd3iJ0rBFIo48ZCwEAklmb
LvTweQBpUsp5ErsLkuUNpszNaD1LqtTazgWQfPO2+mCsOpAOJzhVkgZA8m1S
grNfOpCOFJfSoAmQfG/kYc5YQa5eXXM9HSQnC8Q20NH7Of08vZieDJLffXym
q44CSR4409g2BhSeJcXAkUog3VQbYy00gLKJmnp2qy6QXluH9MY6gbLlijZ1
fx/y/WNWvmFA2bqmvzTxOJCHdH0eXg8Air4y7fVEMZDeuhczOBKBskP7ulqq
LJBHNqd8b18NFNK+pO87B5D+fLLEhYNAiUxoU7ndAWSQju2JHZxAYYl61fd5
IRes+Pm5ACU6jZ/BwQVksAT/ltueQInPtj7DNED+ojkVqAuU9NsvX/KWAxmS
2rDLvgMot9tGTnukAhnaerks2BIoL9nnu2R3AxlVxL3mTRRQeuuqjlGFgDzL
v7RPSxgorwOPSoo3IHsrPV3WAsrw1JTLKiUgWUrigTUpQJnq6pv8OQPkuVS4
7IruX75SsdweCeR58wSV8fOAsTv5XW/egpzNsrkYCRiX6Ibt9R+Qpwt/vR4G
bNXZ3NhKKyBjk4IoR6wAkzicTMmRBDLuqURfkAZg63V81QPQfiZKrPDG8AOm
+lPlhR+JTMY7heQCtrF4IsB7NXIO4QTdgGmt2V3rdgzIC0oSxZvqAQNeyx0W
ekAmqfjHr+8AzKiB5+P2aWQvur6TIGAmIY3xxDXkG98UVukBZjGt06OF6pks
+3uwbxVgDr3KrvTnQKbwxOXJfQPMOXGcnRaObEB0FkYD5mKWky+5GTlQrMRW
BbC9Dyif+DORR8TaujkBO5zHHTi7H8iLRf6cLdOARfiPJjydBTIN0/3wfAiw
qCWvXoEh5G3TsotXADsXPSNthc53mntfRd57wOKuLt5+noRcsiN0TztgaZ0S
bT3rgUzXTjceegdYpnOWsLgIsuMnKv4fYJfH5Rzt55BP+Dt3jACWs6g23tuM
XCGWNfcTsNsqpuwDrkBmyOe1F4sDVlzZaU7bjqwfHleH1i8D+ySXDcgO8fVU
R8AqnfbKDv1BZmn+usgNWPXo5AGZt8hZT/97SQHsvu+REtcW5DtJPoXlgNWz
TuqNpCAP3lDWzAKsQWQlSi4Y+euPryqxgDVls57udwMyk81HzsQLsCcVKc7j
qsjyEwyzQsA6CPyaohjyptHkrBeAdT69+sFzHtmIv/bTZcB63hcFfmhF3jNo
rfwWsF4fzYfKqJ8yvcN/Fi0C1v/nHtehi8gn7Nh36QA2eNbA8nYIcoSFRG4K
YMPCzSmf9yLHHX626ihg77IsBlV3IKdWtbx+Atio0gs5HzXkK+uumSkANn7X
yasU5Vlm/uP5gAbAPmwbLpteQC6K8NtXBtjHdve5jSPId7xy5t4A9oXxycC/
Dbn6hAqPOmDTI/6sO6XItYW3DAsA+37497MfqP8y6znqTPwB+zkXKrH5FPLj
8JWUecB+R3G6BKDzktmoIJYiCti80PncKrN//pyctADY4mXhj783/nM/tRnV
c1kxXXMrBblh0szMEHC2O9LBQUv/nofXyu4FnFM/99G9MeSHvoIaPIBzP1Hh
mUd5nXl/4o21JeB89mU2uigfMqvCmtWcAOd/p5V2KgO5XM8d5Re+2rtu6MEZ
5Nu4Q/puwIV/G639i/I2M4+yVXQacLGIJ4cNLJCztuRr3QZccrXN3TBN5JTj
OttuAo5lvpp/jCPHvEixzgOcpuBiyLaMfGan5VEtwKXL3scYjiMf+8KWNQA4
Xe9gVyTKu8wDN60LJQBX2BXgyoX6KdMyWHb/OsCVXmhh7v/u1x+YyKEDrrJz
trvp3/uo7Ts7Ugi4unWQcaQGssBIt/IawDWebV16j/I+Y76/3r4I8M0Wf6oB
1Sdj8idflwHg/5mFKK88Qa5P08ZSAddt0x11Q/Mio1CK3IuubzNZzKpH+5WR
0rzn2wzgRsahQqHo/GS4h55dagPcyiDs58K/fluIGHokBbjNQyghUZ6mv5Vg
vZEF3E6f/cB91M/pj9tEbE8D7qgTORCciBwVa/ZiJ+B7t7Dqf6M8TVs+8zg/
FvD9FTuCHeKR33B/TUL189zEq1mJ5l9adeGf6VDAvTfG5B+vRz5sqiybBHjA
+rjz3x2BTG3vsBX7CPiJQitj2wfIOeqve8UBP6ksuFQuj3y8+D+Pq4CfUUz0
80X5l0r98fvRNsBj1iQzvrBQPrnceGJRA3hszi4hy0/IKp1ub+8AniAr1lZk
i/JsTqBRFr1vilSqrrc0cmKHm1k+4NmUDPoHNA+TKyxNDHQAv5rmNGBKQw6R
PZVlAvh1CTylAOVlMnAfMbUFvED0MrcnyvukdsGfGzIBLxfM/jSK9v9C11wP
jva7mSO36q04kAnlkyFJpwBvrQkNFUT1Sth/TWm4B/D2I87b9fiRxfVj1h0B
/Hnf6p5MNC/jj0r5/ioAfKAweMYO1SuO7j69VxHwGdudyi1onseYlEeGoOs/
uFVmfqH1oyeYf5Q8AP9Vy1Ujj/Y/OuqthuBGwBfW1pqEPULz7GHAh/zXQOWa
U9yvg/qXJd8Z0qAOVNqVxawStD8RVXXrRbqBKrOrb/8QH7LRqyf5DKCu4b2z
nv8skOHPG9aNnQaq4tED9w6i/g4bE1FyGwOqxvYXvXLeaP7+PSm2eAuoJp9v
CaWh/j45z7Nqvg+oR3UYEaHcQB7lZU2o9AL1+KP5L99Rff2DTTYs9QD1xI4c
5wPovPt9+MZncxaop+w/bdyJ8t23mWdYqxyoLN+wYXm0/hE/tZDjWUC9lFf0
3xM0zw4mc6SP6QK1UYj9mziq1+6Vr8y2v0BtTr3pEn0X2dp+sMoGqG1SVq1L
KH9cLmX2dK4F6jPljCsT6PyS6g9YAwNA7TNWs7yH5ruTpdu5y8eA+inY+aYr
O5B21veiHwsDTWy0zLWQDcjtx4eBVxxokjm+/V9QXhvHap7UuQU0fLe6nQaq
n1GOvJvHMaDJ9JYY38PQ773m14rmukBTflKk3Ib6T3+FZTB4D2gGpQU/JlA/
aWke1+JeBzTvkKvRcveBVDgrQpN5BTSfra5sniif1h75u0+mEmh+v2ROFqJ5
Js8IknlUCLRAn+wjGmj/19BXz7t2AS1872U7YgpIqdw6K5Yd0NJ3pEu5ovwU
c78R1JsDtEaxhDL0vZjLT8wtp2SB1twYfXmpBph/o9tquV4Cre1YxDnXJWAu
mXT/qR0F2rOewD1ro4E5/zCSZvkFaP0pe/iLs4H5K2sXM7oYaF/FVD3rW4D5
SWRHt1EDSOGiHVJjODBfDh02yC4GKV9RvpgND4GZO/E0dEIVpPxTOrFX2sC8
vrfge60HSB2XSL55uhyYOYPrNBu0QCoYp7V25gIzu120d0AepCLpatx+McBM
Tz/3hyULUulqdpEVu4B5fqaoxyUMpOrNL4fqfgCmX2ab5qQ3SAtFqh6zEAWm
drV/qwsHSIvqCgg3XAPmlrLES5VvQFr8+8cSHQ1gbipYFhe9BdLUfQVTyjbA
VE/qbpSNBWkFQzlXnlhgKtnqV1nWg7QuG8W8iQuYkqFfq2QOgLTnGTY5vd/g
/COlJSZ+GqQfnHzZrTYEzjdTCz7TpUHGdVyBuFoOznzjZ7LvfgaZZRlR7bI9
6P9P8g7DkOsge3so4jorHByD2xejGnmBbiWUJTNVAg4X9XllnBuAvtD2n5Co
HTAi+zZ2lyYiD0lxcMciV5qoPg9D/vZ3dr4ZOYXUtSkH+iKtaXhUG9km6tAJ
I2TfnSUVNGBEtPbgiXNAX5I8aM0YBUZ4zd46a/S8v+7p8ZkBwDiTaSyr8xU5
+GR4fAnyCZvCUCfkhN0B4ZPIDOYkg0Sukd99yBlZxC+gKQPoy/ylqrrbgBEa
kz7+Ej1/+U5LxxAPME4H94srFCG33nrUTSA7Dt/bn4f8Jq6iORh5y7iw9Q+g
r3DbZRV/AsapbzM9m7mQnYcPhz5HPsRZ7EFD9nnsdowX2ZT3nc4ocmSe/QFA
VhRoap9ELj6kZ3MXGCHvxXLpIsgNVhuNPiM/ksz64Ivct3HtVgXkK7ipfwjy
ZzFsgyvyKVp+XRfyyqwAPQOZKV3RIfg/tAvwxg==
          "]]}, 
        {GrayLevel[0], Thickness[0.01], Dashing[None], 
         LineBox[CompressedData["
1:eJwd13c8Ve8fAHAkfCky7j33XiuhhFA0RM4nyt5k3nto8LVHIamsIkmyUpmR
yCY7pewSkowy+4qMCg0z8Xv6/XVf79c95znPc57PeI6p5ykxByYmJhYmJmbp
v79qYdMv3k456IJm0C5oydLEz1BHvWgarKA5iake3+SO+1DLUu6mK4KW0KBE
hfwIHrjzfIiE0Dhoads35Snm41epqY/7O/RAy8NaREKhAL8uuvFQOEoCtGJy
ee96iODRO0/eOa27Blolxua3ZFzwBGrcs6UzlqA1PXzEtGQOT0xprKJpNIE2
61mjjTds+D3RXyVHxBVAW3ik/YeaMp620yIzdIwDtA3eJ6XoK+NZ1Jsa1+SH
QfshX6xnvAD+0FbyBSU9AbRrfExNN67i2Sl1qnk8uqDdQb9f61aG54rO7W+f
rQTtnxNiC/mn8aKdJlJ8hdGgo5r80ojnOl5OjQw3It0HnRdvLKZdXfDyJxZN
dzougk6Xtens/AG8wlac5WOYBeiMltkssNjjlSm1Qd4LW0Dnz6iWbPIgXiP6
1T+hxx90Dx1vfNU7j9ftNHQZiDcB3ewL4QF4F95Evb6rjCME9Dx6ElUKK/Am
r8cRsVrSoHcxo/B9Gzve1Do47RnWDXoRVutHn5DwZj+5fNlNEqCXEaFuSArE
W7p75LL+tIBe9+J23HEZfxUluj/x5xbQ3+d0NdZiEe9cr9AIGL4L+hP9+1WT
P+E93+eLM6ytwFDaOXtfMhnvpWVymg/OgaGihWbNhiTeq2HmwEGEg6Gq6mZ2
rid4b0IlzfN0GRga/JRQapPH+w5eDDviwQ2G3ryDb+dV8f7LbDYfQpvBsFK0
+FXAMD7AIcTGW7APjFQ+BRVkKuGjIlr0kA0uMBYyGEgIccbHbzROzZpIgvG6
UNmIjzU+/esurtx5HkwGNU47aNPwOWbmEd5yDEzv0ZvtL3fgP1TftGuMVYKZ
dVV9o3gn/uurtWCqSwOYc51z8lyexpfcC3mO7FwA81wmyYDGXnxVUnBb/goL
nNBXyJ3npuJrRWaV9hsH4MSI1sWjv1vxDemWTqWr7WDhZCA/mfsRmGZzV77E
iILFcg/vu6uawKLQOzK69AUsLx7wEiQNwKZo7rbamSywYuWW5FKQBtbx8bkc
YW6wCt/1KaZGEdgk7jtWOEeB1e+pVM8if2APchrPh0WwPv9xXP9hHXA0ksub
9vmA9ZQMuycrH3CKvb3RVDQKNoS0Bn/UNuAK2J3c4/IbbDp1aiNVp2FLreM6
+1EZsFUMETnhwwPcWOm/YplaYJvqqaiV1QU89P922AmfB9u1O+l3vKNhm0E/
14mvskC3kNq9riIL21abOOZ5+IH+uF8gJisWeFPEi5OfswKDtaHlBf0j8On9
86FWyRUYVjzRh7fuAL7FuvundPKB8SB/bk2kFPgTP2+jfVkAxjhHQETgYxA4
lrBJy6wCiB3HHHOr7EBASyisWi8CCHGhzPuufCCg85BjtwYd+depyAYJEDCs
3MqlyAKEZKa7IY8/CFi9J7/hMwRCavkR72wACLgLSVt0TQEhfzNQmXUQBG5n
mTga0oBQtdwUJf4TBD6Xp4dpBwFhblsdYC0EAtMKdCeiG4gTPF5WJ4xB4EsB
Vc9XArkho1fHEwTmH8TzZrQBYbnLqzDXHARW48JTfwsAYT1TOSS3E0jcXu6V
hY+AsLMsmKhtB9IBWZUpvi4gnF/5Pb4SDKRDOcuvd+8AwsV1vHMoAUgq4hVF
uC8QrlubrJZGgASC8j5uaH5uRn05zelA0uUUZ2o+A4RHB9TXRwKJMcVJPb8M
xLmCOS4zKSCFPXivOywGxKU/nUahN4AUcaJYpSESiMtJ2/MkvYB0gz1cNucH
EIH7X4wKzwEpxlWJ+2wDEEEu005LK0BK3nermwONH9LeNRXZAKSSei1rRTT/
cI9H9murQBocqfw3Ug6IW7s+M9moAmkk5qaVZyJyoZqLaBuQ/lM/o2O+AUTM
3vfDLJ+B9DmbT0YUrTf2ULdlOzOQvru7z1Z4AxF/pH1mCAMy25qk78dyIO6I
Ryqf+gHkfdSEqwcPA5Fm6+z/8gCQlVIPUE/1Ib9YN07zBvLB7f1F0eeASJcc
G09mA7KqFO39RAHyt3vfhVOArHkgY0+CCBAZPpUOe1SAbGNW2PcdxUPWSSJq
shPI9H4jN6E05NrW4Od/gGxn851ZWwWIh6SA/0SKgHzmlJJMug9yy/PoqhYg
e3g9CTb4DETO9iqpB+ZADr3ZujvvNRC5j/nVsjOAHMbrVNfnBEQe06r+Mx4g
X7vNac7Cimy00/yuE5CjUgwCrdWQp/XmJeuBnJj37h17CRAFAo1ZZej5ea0f
L51JAKJIG2Ox4ARyoV4ob4wC8g3x5UfXgFz8Rjy7tgO5w8tMbDuQy/sc3/Kz
AVFs1Nx20hrIdePfJOv9gSjRs36XxgHkd8wrb0RQPjwW7nOzUgByb23FWSo3
srU6DYsF8ntfbxJ/PXJ8U8GVp0Aenpqy/WcnEGVsW3I49IE89abv88855NGk
9S0RQP5yPe76bCYQ5fwlw82OQJ7VMJSdNkfW/PC18hGQf1Y2nx2pRs47/YTP
DsjrqWXrr0KBqHDUXlBWBozZ0jOjSQk57su+sGXAWHlljtVNItelDC+JAPbP
1czIcn0gKvllk42RBVxjyekk5JJ5F3IPYJikfs29l8h9TD2btAGjjnLQEwKQ
V4VGg58DJmoalHl9BIgq3Pd7aDVg0soecj4oXqvq0hLYDwAm+3P3W08b5MH1
NJE7gMkXTPi4bEVedvTu/wrY/u30J3ZngaiW0wsf4wTs4ACFYYPyu1qrtUtw
GrDD8T1M5ij+qu11frRxAQbsepq6KF6qoxndrjhg6vVs08dmkR/M3v4UAtjx
gIYo/D5y9RV84BpgurPK3ftRPFSPVJo6SwJmkLPgK1+JPG/2QuINYMYnSynS
zkDUMP0UyYkC7ESvFCHaiSy6v9JJHTCr6HFmWjCybP9PN17AbLXTs0iKyMoB
Ug0fASNYbLV4UHzWHBfsH0PztX9KnuG8i2z0NBSjAnbat/vmZl1ka9tHQUqA
OcpHy2/8QT653M7TD5jTtE73KorHGud497gHgLk+2Oy7cArZS2ZrtgVgHvR6
yjyqZzW+z2XrBgDzJl+qnWlFDjDSk3QH7FzXIWIC7UfN5YGp5QrA/K7/Yv4o
ixx8Mlz2A2AXNIqzBtD+1ISMZd5A41/846rVG4McynhWHgdYYNXOmS71v/93
WwozARbiNXazbeHv/fjstsOAXVlz6uUa+v/4x/bYABZ+bU5Iv+Hv86PdVRIA
u87ve/pmLrIfnWlrEWA30n7ndf4d31s93qIFsGjpkB88fsiuR5/luQEWU8mh
bMxAPmNNulQHWPzR6ODYY8j0iL07nAG73SHQ2i2NbNaZ3fUesLtWyTz825B1
Zbek7QMsaVzMwmwJGU9fDrkEWKpnTmrC3/UpSTOZVgGW/nvPeG8T8q5XNW1j
gGWGl8uS85BpgU+uigP2kO/wOctY5C3ad1Z6AcvbrcX8gUDx8I2q5uEDWEF5
hw4Nzad6SHRv2y/AisEsxlYGue2QRbosYOWW9iJDy3/jLTNCMAawyrHPjsIj
f+PxR51cEmA1Hm6FRDPyBXP/WwaA1YVdUPkYh6yn+7oZ7X/9to0rYqheVO8b
uYnOlVhjSljbKTtkaqh4nilgL8virMbR/lV9WtpujPK1+79838kWZPeowBX0
Pnvd9z6TQvW2yiCeLzMMsP7lalbneGTZbGxpFLBhnqa4L/YoPyfnGit+AzZ5
ZLh4dhXZbHPh/XuATb86vST/EVnm9PC9rYB9NZ9R80LxVMnyWsAPrf+762L7
D1QPK4oejzr8BOx3Es/0ojyqL6uDo2Uon9YlE/ceJCO/9jui6A8UplIh//Nr
yMlk+ff6QNn8cjfbCurH5QcdqPlBQOFZVBf/g/pZ2Zl9w7a2QJEw9SFYUX6U
Rss+i6cDZefb/djpdWRjWUWdCaDsNlroakTXl/IqadZ2A0XO4LxGKKrHJTH0
dPtqoBzSDpDaQPWo+ApL4D0zoOirBf1cRfFQqHx8LKIJKIbPoNAG1d+CL6mO
WCFQTFSZHWtQvBakLI4GygPFQjn0g380EPlrRWxjK0CxVwqrW0T1N69EzGkp
HCg+0jeuf7dA/eZzndpKLVD8cvU1jJ8iX2cpLWMBygWpLWsl6LyQI63OoGsD
JVAy2tMD1Z9s59LCpBmgRGyPNf8ahvrb0DHjEHegpJDviE6i/ppxa8Am9StQ
mlgyK0b4gUja87pP7C1QWqouX96C5nPvQfdhJW6gvHKzOqbCiUy5nU1ZBkpn
39buu6h/3lk3OXlUDCgfcv3nTND54faTqF/m00CZMzaSakb9PWbhbmfCY6DS
Un8nF0ah88LZHTtkE4EqbNp3aogDiLC6Kq6sXUDdzl4qzXkV+Z8lNf5IoEp6
O1b/exGIKynn1Qs/A1Xh2NteMRcggksdpXMGgHr8yyPu29pABCSF2pwqB6q3
snnI5c1AuDu07X7nBdRzz1e+fkfrc7udTjabBKqfZrqVI9o/1yaVAicHoF40
m5E3QvnoIty2yXQeqGEeQcM70PMdG+TScmuBeu9B/qGXqD7bvSnN2ZQP1AZu
5nl+tF7jVa+SPD+gNiXk2F57jMxibTkdDNRWQf2WNVQfjTh+b78lAtR2qTup
E2h/DficNIefArVPY49eNep/OqQh/EoPUGf8rXIIZiBgpLkl2gxofGPFRC4T
On/6XqCgcWikdI/+ryi/5FQsw+YngEahy5kooPcnu/FQ5u0NoAn3FmpUY0BI
X7kxkogs9TJfqhWdvyT9z2/ZOwk0taLsHxMo3gTJmnpmdKC5BKRdE6sBYrPw
sNXxdaC5HySYHAKB2FQS5XH1HdA8fwlfyEX1heXo3VG/fUDzdU9xU+gCxgbR
mj51AmjB9kkm+BQwVsLcdS3sgJaomShIkIHx7YodnUsBaA18N4uTzgKjR42a
62IBtKaGa0lrVcDofrB177c4oLWeDQkn1oDxlo2e580HtPZuX4b4NWB0vNzq
3aEDtP44BmdBCjBaDz6KN9UF2jc+WYe6ZmDUtFpW55aBIIX3teAnCjBSH6pU
bT0Hgh68HBEyz4Bhf1SQ7ZEtCHrFdWA9B4BB0HXiTTJB8JxAbM6lEmDQfb+z
+W0DQX8KraUjE31fZOUeP1YHgqGiezZ7RgDDeK6iev42CCbuMQktMwXGUZVf
/EqLIFink3T58CQwtuc1FnsIgRB3qOxZXV6gD3QVTS4VgxDvYS6e+vtAf29x
2A3nASH+79OFygpA7/3w6X5MOwhRT2ZPSRkCvWvw5ttzHSAkcVSMYIsEesvL
Q5c5ZEDoMBNZp5EVfQ9pZkpdwUDIIZBJTGUR6Nd6n8eeR/c/vfCua88Q0PfI
rbHVfwJhYlwCTysBW09e7eQeDRBeF+Y9UMwAm6T9ljuMkkAkbygkIywYrBPF
E7Im9oOoPney8FQhWF13Ge6N3ADR1dZD3LwmYKmveGLHy3HkIUGWzZHIXN9m
4k4hz/9ZWGkCS73xcHowE4j+pjUOjx1Avl2p9oYD2cOosIwGlrrLvH88J0B0
jfSvgfkYWOo8rTljEwGif04nRt31AUut460TjW+Q/S8ERxUii5xYeLiKfJPu
E/wZLDUXPxVPnkWu2kF3tkLOWRe/VQCi65xFsoePIHPKs+fuRi5tfj3EBpbH
OoOb0q8jtzx63oUj5/xz0bwLefBGWZM/clDs0DUuEN3YbJJcMIOskJHm3I5s
Nex6uRMsNeKebad7Iru/sDvLjuyqIbPHDDn0gZkjIGu8/Bi5GbnAWcXwMViq
/3qnO3scuV5fXh19b6t3WG6sXEPukxc/KIH8cMjuLhvyFz5MhkAOtHd+nYK8
scAlegfZYlwkUep//R42lw==
          "]]}}, {{}, 
        {GrayLevel[0], PointSize[Large], 
         PointBox[{{8., 68.93231469526532}, {8., 81.91168472067831}, {8., 
          93.03654847288784}}]}, {}}},
      AspectRatio->1,
      Axes->True,
      AxesOrigin->{0, 0},
      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]}}, {{0, 
          FormBox["0", 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]}}, {}, {}},
      ImageSize->Scaled[0.585],
      PlotLabel->FormBox["\"Temperature Profiles\"", TraditionalForm],
      PlotRange->{{0., 10.}, {0., 140.}},
      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:eJwVx3k41GkAwHFXhjFG426mTB5aIrJypEy/l0jrmEQkVxKryDieaEyH3CpD
yKZstvC4qik2yzTW+zM7jkqyRU2PlmpzNQoRxjj23T++z/f5+MSFG0bKyckp
oJj/H/MOX6EFb06GhmZJlDq/F9gp92O8kupc6BBmvL2BJcbYP1ZEtZVHQ5fc
z7EUfx0se4uzyJcaBA+KfJyrS6lYCV3z262AQzCM4qfkPJ2KVZXJVIVnvGHs
/N3uxbJM7FFgnI6fci88d++vnX9yyzGRQx3NgMmFea6imuYGEjYgX3MivuUS
LI+xmln+HoUNW9S3FWfSYA2v1W0ykYJ9KvaE4ux5+DB+HM8qS8MmZHMd5ZaT
sGk+Rj6j7j42pWKrqnSaCAVV3S+romXYN+Cz2NdLg+027ixW9BVsbqrTbYiC
Q9HbaD3TggVsIWFdk6BSCru5/dmZc1uwJdOQZJsaPuyR3JxypGpiyw1pfueX
1GHf3fOjDsWB2JrFe40s2Qx8Jb/hGsWAAeSmzfO1m3lwQPCs0uFCKlCw9jLq
Lb0Pxeyhma2Ft4FiAefja+t4OFibrdKs+BYoje7TYBjJ4LBimy6PlAiUtyx4
kBz74AeLy5EKYYqAcPHRw63kE/BfWk/cZu3nQEUUYMwjuMMRZU5PoH8dIBoR
hBdbpuHYruV3QUpkoHb20PNGYSicyDEtOtCzHZBac8Y9Z3Wh5KnAfn+6FJD1
u7bZvCuEX4/tiij4vgtohEymZn7xh9PJao0sdQWwnnn4NtXsG5wR7eO7yjPB
ehmmoy8og7PYi6o+RjWg/Hqj8NnRCTjXIFDKr78OND04RsbvS+G8b9KQy/wE
0FzYGmOpbg4XVlX0HtvVAa3rHo8lXseh9MYcc8HDA2i7zNEPnFyBsigCbzTX
G2iPaVyJuukIVwyftz/g5gKdLClve/JRuDpkdoSuIAW61idHyi0/4XKR5Mn6
J9NA95WXMS4j4/JM80rNkk6gl/YH+3LlDVzBK8fK1YAN9I2L+lWOPsEVGXqB
mVUdQF/EGLXt+gdXMsnjEurdwYaEOHUiyRdfZ3R9ydPEElC1tHYXRbBwZeKb
iUlaBaAKGRfaSyW48uqw1YgJGdBYduzbKZ44YYR9qdZDBjaSKzg7Yn1wlcaU
iAN7lsHG1mtXYxl8XDWjYHbPjhKwKVTKCzLswIk7C01JNWZg06pF+mJDE04c
9DozucERGNSPnWIYEHA1jov4A2sZ0D21Q2xvreJqK5cjx796A/pSV3j/0mGc
lLI7j7lbjPzOdVP/IeSfrO8UhCLPmPzM80b2Fz9MLQB0GVWyGL4fOdFY3v4L
MiuJ9nQncm0zSzcJ0Jd1csJv6uMkjs7fSdIiQF85fnfaQYw89dKpl4/MznfQ
6kdeiVe9r4fMjU/70oeTzqqpc+3nkZttte50I//g1uVnDeirxDZ71RbkYD5B
wEFu6Dv/thT5ST6/2gy5s7Hz92vIb0wdc6TIgyUa3KvIn9rPGtgA+tq6I3ec
cpHXZs1PjiEHfBTVsXHSOTvfgDwj5NgO9YzTyHs/8+nbkNNr/EPikQ+mDUT/
gnwvZpxyAvnUA6dtucjtTCtJBDLH9WoVGfm1VUpHGHLOYNnAFLJEU/hbMHJJ
QrjwHvLad7WUI8iVyiPRwv8AbHRZVQ==
          "]]}, 
        {GrayLevel[0.5], Thickness[0.01], Dashing[None], 
         LineBox[CompressedData["
1:eJwVx38403kAwHGF/JqN/WoMq6hOUqmoLrvvRyldGI/f/ZAzdaGSnvzYqEs1
P7pSJs5EyY+4iOc4EyWfT26jus6jJ7qouzwys2cpG2Us3Of+eD/v5xV0kr/8
iIGBwWIc7/8Tgfw59sFlKdDFK7dP9VxFHN8b3VBYnQO9c3a/cB7rJgRuFUc7
bsXDwOYAUCwNIbJW7pAF2x6AkTRqknueE1HIoWpvRoTA+PT6d9GOPkRVid6s
MzUQCj7Wa2tfyYnm/ScZoUt6YM56q+TdFZWEbNtdtgMvF0r87T9vavUg+hfV
xCa2XoK1IwFhM6FexDvX2o7rIjaUhhfHQf4TYuS6H3yd9QV2MCTypMwIQqWf
kt9a9wHKDh0xqtc0Ep9M3c2MkszhU9UFVMpTEVoQpOvtYcPePePKd88YxNSn
Lp9/rRF8qXDw3Zk0S0yfMpY+rJyBfxeEZK8/d5OY/SYyZXNNG3zzUdq8xTyC
+Np4PvTsrCUcclF0tciLiAXXIUqmXgOHMxwHo+VCYDDhcpV+vwEqtno3STRu
YPFGf8ceST1Uaph1D3IrgeG1tOFXGxOh2oaa8DzsPTAa3U3hOurhx9jLcTGT
LWDJymlfkmcvnHhtWbLX8gwwyWj+zZkcC7UPRPEZQg9gKotwajDZC6d+EefH
u6UDc0eTzozWCfilpdRcMaoEFukhfzV1HoI6Vahn9S4uILVnj/lNMqGe0u3U
v9YHkFndaze/FcO5jDUTapUCUCI/nBONh8GFoY7mvO9KgRUv/LbtGi1atCot
zGMAW08wWA9L0OKL3bODifbAurRY/GeUChn2J1fThUWA6pvm6DQkQcYBd129
BtYA6rTzsXWWLmhJT1UwWToIaEW+D9T+MciUmztf7y4HdO8pTkDcHDLrpwbP
iycBXUm5fPSGJ7LYl1I3UCYDjMyZhvUpUYj09tL9kXIbwNwYp7i1bgSRs1bd
9ssKBMyX/k5IT0ZWFqfD74VOg6XnWwQ/VxYjK0WZ0K6UCVhO+X2mUU+RdaPz
8+rPEsCScUfdu/9B1NR4w19JamBz6qSlOSkY0XxLekjiUWBLo23PP5yA6JbX
SDvrvge2ndyfHkvUiD6QT3+6UgTYCR6C20I/xCjuM6LbJQM7ckXaphNBiOkd
xUp5MQjs2gvyTnDbEFOju7GC5wzsD800HFguR0uztY/afYqA/bzrBV2jFLHI
qWf5iXrgUKs8znUwQSxxvFYqGAYcP3qk+815xJozeCTeLgWc2W5+32w4soku
FVYRFOy3u+z7QrBrMnkeK7A1q39sCMRuymtX8gBHb6vW8fdgP6keE7ZhJySz
n23FnuoVHNsKOF8Z2fwbLGTD92elvx8HnLmYuoltr7HnztC9M7AFV7fR+pBN
jGlG2H477NzE8+O92DRR5BYd9n13WvkTbOfLbbzNgDNv3rHFrBU7tKhWswy7
sffsgAS7virp9Avsrqau3wuwW+8Y7kjDflNIyc3D/uNO2FAq4CwY7yv3ysEe
uON3xRU7Ylh2V4BsDhtXnGrQYp+QW15MwrYqX21UiH2hJiwyEZtdVu5Kx753
bMw6Ftut5L3uOPZj3gb1YWzP4nulq7FfbRDKf8D2KfqWLsJWUzvLDmIHFVyJ
nsRe+Gwh3IcdKa4UKf8Dlt1b1A==
          "]]}, 
        {GrayLevel[0], Thickness[0.01], Dashing[None], 
         LineBox[CompressedData["
1:eJwV0Hs01GkAxnEkzIQxNzSDH2asLi6tHRXTnN9bW00bJktSue0OrSmFttLQ
yjWio3LLUCR11kZmSyXKel9zRkolu023jZLThZiSSy5D9t0/nvM9n3+fgDip
ww49PT0DPMn/Jf2ls9xQ+wToWSUlp8R+5O6NPyuLfz8K/bPzvRv6tKT826ro
lopdMMS4TtAgLiWznNaoAzkhMLrmziLJkB1ZTDBGyrduhgdu7vJuYiSQF07r
KKqD/jCzdiI8aqOMvLY9jh1k1AmLgnZLXW/JSLXXRa6dJA9e6Gwv0z9yj3ys
Xy2Lb8yBN4rOrIzOGSNfuda0FGZyoWqqace7hI/km0Jf+CzrC+yYa8vf8Pot
OaAba6twG4L/yGx6CnNDyU8mnhTD/VT4NONRXn2mmBwBAZNdnVzYU1UaEnzJ
jhz7dFv8ko5gH2FE/1u/j5zYO//6rfNT8G0iIzElZTs5vSgsQVDdBAeiTwlF
lqPkzJW0oORpM6gN4aY4reKQc669tCO6z3AYLiNc+CqgN7z0OOuGEo7yOjbt
C3AABh5+vE5FHRzP0M9N+t4GzDuR1PfEIx5OCrWfaY4QGL5bTxPxdFBH42Qz
bH8FRk4TPqaruuCsxiawpJ4DjFOvXV5sLkN6yd66on4KMFFv5SuNNyIDo1dn
hzepAJVnrEptHEbz1IbPpltLwIJDmx/Uq8LR/GO5A3yvSmDanN3vO2qJjNdW
a0+JtcDcut1F0J2PKAIW7eabdEALG0rJ1G5BCxjrmmcUh4GFJLiSs2QEmY4v
r2y95g0sdCTb+tZpRLPwWfRb/2VAP1Oafy9iAFnkxlsRK84Chk8Sj9+rQPSL
4n97bbIAY2JxjJvZUsS4G5QsTPYDzBKfm4N+kYj5Yvlzt1dswFo7RmzaOYvY
lBXiWPYwYL2nHYsuW4Usna0CJfIIwD4ypXRPiEBWa9aH17aNAkuPnW8r3N4g
64M9oMA8Flg+8uMjnTlaWPokUkj/CKzSGuS550sRR8P6gVHRBKz5BRqTiLuI
O0kVixf+AazVonee7T3Ililv/SLQAwv3xplRTQOR3RqDFINxIeAwmcKCqFhE
BE2H5mxrBxyV6HCrYhDZC0dRLa0ecGOXyysTfZG9OnVkSQwL2JhXJX23JwA5
iOouxxWYAZvmopN7RE3IoWmJEXRxAbbhU8oQhzbkyDepsCj5AGy/uqZPXrmO
HDPZj/7yMgJ2Ne93i+yMkePrjl6qoBMQvqwwz/KviOc28NB8bSEgptulmulg
xDvIqLnf+hS7e52tZjM255R7xlLsz86/KP2xHdOi5VcBoeMMTko3YHsEB+s7
Y8ce4HasxA7QF+S5A2KGnS0ts8bO9x1osQLEbGTtsNczxJPT8DMG2PLjXkwN
tlWN+MQX7Lz4NG0XNnHym/37sG94Ms/dwXYPidRqAPGV2rKC0ogt+dRT3Yh9
pSv5uQL7uMX9KBn27frbV4uwi8cWd/6J/aKYlncSu/xpaEItIObmbzu3+ih2
Xfnq1BzsrX3qi3LsB87UtF3Ye9rMMvZjP6bQX3Zjp1dvCYvH7v7w8ZAC+1JM
P12GPaR03teA3SpZNhiFPXYi4+Fh7CfLEtt+wp6Jvx6zHXuQoTobiniJhj82
i5Ox58YXJG7DNvVQRHH+A1mnR8Y=
          "]]}}, {{}, 
        {GrayLevel[0], PointSize[Large], 
         PointBox[{{8., -3.56086}, {8., -6.34801}, {8., -9.73642}}]}, {}}},
      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]}}, {{-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.}, {-10., 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.394963883423554*^9, {3.394963924388444*^9, 3.394963945464834*^9}, 
   3.394964035841164*^9, 3.394964066630714*^9, 3.3949641101625347`*^9, {
   3.394964153732695*^9, 3.394964183568836*^9}, {3.394964219661434*^9, 
   3.3949642570899363`*^9}, {3.394964848354051*^9, 3.394964876941865*^9}, 
   3.394967179169877*^9}]
}, Open  ]]
}, Open  ]]
},
CellGrouping->Manual,
WindowSize->{641, 727},
WindowMargins->{{-1, Automatic}, {Automatic, 0}},
PrintingCopies->1,
PrintingPageRange->{1, Automatic},
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, 152, 4, 30, "Text"],
Cell[745, 29, 256, 5, 30, "Text"],
Cell[1004, 36, 494, 10, 45, "Text"]
}, Open  ]],
Cell[CellGroupData[{
Cell[1535, 51, 339, 8, 45, "Text"],
Cell[1877, 61, 183, 5, 29, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[2097, 71, 294, 7, 28, "Text"],
Cell[2394, 80, 304, 9, 66, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[2735, 94, 157, 3, 28, "Text"],
Cell[2895, 99, 130, 4, 29, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[3062, 108, 264, 5, 28, "Text"],
Cell[3329, 115, 412, 15, 31, "Input"],
Cell[3744, 132, 412, 15, 31, "Input"],
Cell[4159, 149, 412, 15, 31, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[4608, 169, 397, 9, 28, "Text"],
Cell[CellGroupData[{
Cell[5030, 182, 360, 8, 29, "Input"],
Cell[5393, 192, 962, 15, 22, "Print"]
}, Open  ]],
Cell[CellGroupData[{
Cell[6392, 212, 359, 8, 29, "Input"],
Cell[6754, 222, 955, 15, 22, "Print"]
}, Open  ]],
Cell[CellGroupData[{
Cell[7746, 242, 361, 8, 29, "Input"],
Cell[8110, 252, 959, 15, 22, "Print"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[9118, 273, 177, 3, 28, "Text"],
Cell[9298, 278, 2408, 71, 178, "Input"],
Cell[11709, 351, 2408, 71, 178, "Input"],
Cell[14120, 424, 2408, 71, 178, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[16565, 500, 384, 8, 45, "Text"],
Cell[CellGroupData[{
Cell[16974, 512, 434, 9, 29, "Input"],
Cell[17411, 523, 1003, 16, 22, "Print"]
}, Open  ]],
Cell[CellGroupData[{
Cell[18451, 544, 357, 8, 29, "Input"],
Cell[18811, 554, 994, 16, 22, "Print"]
}, Open  ]],
Cell[CellGroupData[{
Cell[19842, 575, 358, 8, 29, "Input"],
Cell[20203, 585, 998, 16, 22, "Print"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[21250, 607, 954, 16, 96, "Text"],
Cell[22207, 625, 160, 3, 29, "Input"],
Cell[CellGroupData[{
Cell[22392, 632, 1274, 32, 102, "Input"],
Cell[23669, 666, 1359, 22, 29, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[25077, 694, 261, 5, 28, "Text"],
Cell[25341, 701, 3545, 77, 210, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[28923, 783, 413, 9, 28, "Text"],
Cell[29339, 794, 3547, 75, 228, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[32923, 874, 339, 8, 45, "Text"],
Cell[33265, 884, 925, 16, 48, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[34227, 905, 222, 4, 28, "Text"],
Cell[34452, 911, 3578, 81, 266, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[38067, 997, 544, 10, 28, "Text"],
Cell[38614, 1009, 3590, 77, 212, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[42241, 1091, 424, 9, 28, "Text"],
Cell[42668, 1102, 878, 15, 48, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[43583, 1122, 504, 10, 28, "Text"],
Cell[CellGroupData[{
Cell[44112, 1136, 536, 11, 48, "Input"],
Cell[44651, 1149, 26905, 475, 757, "Output"]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)