(* 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[ 142515, 2887] NotebookOptionsPosition[ 139150, 2762] NotebookOutlinePosition[ 139557, 2780] CellTagsIndexPosition[ 139514, 2777] WindowFrame->Normal ContainsDynamic->False*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell["APPENDIX C", "Text", CellChangeTimes->{{3.3931670131553783`*^9, 3.39316701760128*^9}}, TextAlignment->Left, FontSize->16, FontWeight->"Bold"], Cell["Simultaneous determination of k and Tc with n=1.", "Text", CellChangeTimes->{ 3.3932541151780024`*^9, {3.393256437730439*^9, 3.393256443380356*^9}}, TextAlignment->Left, FontSize->16, FontWeight->"Bold"], Cell["\<\ Programmed by: Mark D. Normand and Micha Peleg Last modified: August 1, 2007\ \>", "Text", CellChangeTimes->{{3.393166039594693*^9, 3.39316605992082*^9}, { 3.393241773739553*^9, 3.393241777356629*^9}, 3.39324181747005*^9, { 3.393241865841612*^9, 3.3932418680779676`*^9}, 3.393330050497901*^9, { 3.393339938926812*^9, 3.393339939354993*^9}, {3.393340931847633*^9, 3.3933409388732967`*^9}, 3.393768549884416*^9, 3.393844866871951*^9, { 3.394191726851714*^9, 3.3941917269981422`*^9}, {3.394965254744875*^9, 3.394965256657276*^9}}, TextAlignment->Left, FontSize->14] }, Open ]], Cell[CellGroupData[{ Cell["\<\ Assign the two parameters k and Tc their original values. These are the \ values we are trying to recover. (In this case, n=1 by definition.)\ \>", "Text", CellChangeTimes->{{3.393167057906272*^9, 3.393167104947035*^9}, 3.393242067205596*^9, {3.3932524597927837`*^9, 3.3932524674366493`*^9}, 3.393252764466267*^9, {3.3933300596388693`*^9, 3.3933300633790083`*^9}}, TextAlignment->Left, FontSize->14], Cell[BoxData[ RowBox[{ RowBox[{"k", "=", "0.035"}], ";", " ", RowBox[{"Tc", "=", "110."}], ";"}]], "Input", CellChangeTimes->{{3.393167044544138*^9, 3.393167047360354*^9}, 3.393241927165762*^9}, FontSize->14] }, Open ]], Cell[CellGroupData[{ Cell["\<\ Assign the final time values, tCycle1, tCycle2 and tCycle3, for the three \ survival curves.\ \>", "Text", CellChangeTimes->{{3.393339719425305*^9, 3.39333974191691*^9}, 3.393341354175434*^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.393166362624907*^9, 3.393166454118064*^9}, { 3.393166577468137*^9, 3.393166582658086*^9}, {3.393166725191359*^9, 3.393166781654152*^9}, {3.393166882404248*^9, 3.3931668998421392`*^9}, 3.393252230849867*^9, 3.3933303398865747`*^9}, TextAlignment->Left, FontSize->14], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"T1", "[", "t_", "]"}], ":=", RowBox[{ RowBox[{ RowBox[{"110", "/", RowBox[{"(", RowBox[{"1", "+", RowBox[{"Exp", "[", RowBox[{"9", "*", RowBox[{"(", RowBox[{"1", "-", "t"}], ")"}]}], "]"}]}], ")"}]}], "/", RowBox[{"(", RowBox[{"1", "+", RowBox[{"Exp", "[", RowBox[{"9", "*", RowBox[{"(", RowBox[{"t", "-", "6."}], ")"}]}], "]"}]}], ")"}]}], "+", "20"}]}], ";"}]], "Input", CellChangeTimes->{{3.39333080816358*^9, 3.393330816299164*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"T2", "[", "t_", "]"}], ":=", RowBox[{ RowBox[{ RowBox[{"133", "/", RowBox[{"(", RowBox[{"1", "+", RowBox[{"Exp", "[", RowBox[{"9", "*", RowBox[{"(", RowBox[{"1", "-", "t"}], ")"}]}], "]"}]}], ")"}]}], "/", RowBox[{"(", RowBox[{"1", "+", RowBox[{"Exp", "[", RowBox[{"9", "*", RowBox[{"(", RowBox[{"t", "-", "5."}], ")"}]}], "]"}]}], ")"}]}], "+", "20"}]}], ";"}]], "Input", CellChangeTimes->{{3.393330809518353*^9, 3.393330815059342*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"T3", "[", "t_", "]"}], ":=", RowBox[{ RowBox[{ RowBox[{"160", "/", RowBox[{"(", RowBox[{"1", "+", RowBox[{"Exp", "[", RowBox[{"9", "*", RowBox[{"(", RowBox[{"1", "-", "t"}], ")"}]}], "]"}]}], ")"}]}], "/", RowBox[{"(", RowBox[{"1", "+", RowBox[{"Exp", "[", RowBox[{"9", "*", RowBox[{"(", RowBox[{"t", "-", "4.5"}], ")"}]}], "]"}]}], ")"}]}], "+", "20"}]}], ";"}]], "Input", CellChangeTimes->{3.393330810824532*^9}] }, Open ]], Cell[CellGroupData[{ Cell["\<\ Compute and save the three temperatures at time tCycle3 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.3933303708437433`*^9, 3.393330378065383*^9}, { 3.393330985699769*^9, 3.393330989439679*^9}, {3.393339427103109*^9, 3.393339429086478*^9}, {3.393772884477716*^9, 3.3937728861669083`*^9}}, TextAlignment->Left, FontSize->14], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"tCycle3T1", "=", RowBox[{"T1", "[", "tCycle3", "]"}]}], ";", RowBox[{"Print", "[", RowBox[{"\"\\"", ",", "tCycle3T1"}], "]"}]}]], "Input", CellChangeTimes->{{3.393322100432173*^9, 3.3933221278763723`*^9}, 3.39332554916901*^9, {3.393325607619711*^9, 3.393325612349854*^9}, 3.393326084597148*^9, {3.393330385984529*^9, 3.3933304059236383`*^9}, { 3.393330996182164*^9, 3.393331003269291*^9}, {3.393331046968734*^9, 3.393331048972906*^9}}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"T1[tCycle3]=\"\>", "\[InvisibleSpace]", "20.000001675297746`"}], SequenceForm["T1[tCycle3]=", 20.000001675297746`], 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.393330481761557*^9, 3.393330838455962*^9, 3.3933310811866903`*^9, 3.393339336616654*^9, 3.3933399263090687`*^9, 3.393340379244706*^9, 3.393341226105783*^9, 3.3937690275234203`*^9, 3.3938449477756777`*^9, 3.3941917306380367`*^9, 3.394202733968403*^9, 3.394207933868371*^9, 3.394965368021757*^9, 3.3949668332989883`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"tCycle3T2", "=", RowBox[{"T2", "[", "tCycle3", "]"}]}], ";", RowBox[{"Print", "[", RowBox[{"\"\\"", ",", "tCycle3T2"}], "]"}]}]], "Input", CellChangeTimes->{{3.39332215124137*^9, 3.3933221528608713`*^9}, { 3.393325557500794*^9, 3.393325579892871*^9}, 3.393325619031274*^9, 3.39332609024548*^9, {3.393330392021524*^9, 3.393330410491646*^9}, { 3.393330997885365*^9, 3.393331004684204*^9}, {3.393331050948956*^9, 3.393331053101078*^9}}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"T2[tCycle3]=\"\>", "\[InvisibleSpace]", "20.000000000249976`"}], SequenceForm["T2[tCycle3]=", 20.000000000249976`], 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.393330481869207*^9, 3.3933308385259132`*^9, 3.393331081294511*^9, 3.3933393367146387`*^9, 3.393339926410222*^9, 3.39334037936642*^9, 3.393341226211865*^9, 3.393769027627445*^9, 3.393844947878705*^9, 3.3941917307816553`*^9, 3.394202734060958*^9, 3.3942079339447927`*^9, 3.394965368139865*^9, 3.394966833377267*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"tCycle3T3", "=", RowBox[{"T3", "[", "tCycle3", "]"}]}], ";", RowBox[{"Print", "[", RowBox[{"\"\\"", ",", "tCycle3T3"}], "]"}]}]], "Input", CellChangeTimes->{{3.393322170793523*^9, 3.393322172397147*^9}, { 3.393325566670063*^9, 3.393325582959031*^9}, 3.393325626854475*^9, 3.3933261158418827`*^9, {3.3933303956527243`*^9, 3.393330414259577*^9}, 3.393331055262907*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"T3[tCycle3]=\"\>", "\[InvisibleSpace]", "20.00000000000334`"}], SequenceForm["T3[tCycle3]=", 20.00000000000334], 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.393330481976028*^9, 3.393330838605731*^9, 3.393331081400239*^9, 3.393339336807788*^9, 3.393339926501018*^9, 3.393340379487528*^9, 3.39334122632506*^9, 3.3937690277338057`*^9, 3.393844947967454*^9, 3.394191730920383*^9, 3.394202734145694*^9, 3.3942079340263243`*^9, 3.394965368264694*^9, 3.394966833457037*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Define the three survival functions as function modules.", "Text", CellChangeTimes->{{3.393340195045807*^9, 3.393340195302422*^9}}, TextAlignment->Left, FontSize->14], Cell[BoxData[ RowBox[{ RowBox[{"LS1", "[", RowBox[{ RowBox[{"k_", "?", "NumericQ"}], ",", RowBox[{"Tc_", "?", "NumericQ"}], ",", RowBox[{"tFinal_", "?", "NumericQ"}]}], "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{"b1", ",", "t", ",", "T1", ",", "x"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"T1", "[", "t_", "]"}], ":=", RowBox[{ RowBox[{ RowBox[{"110", "/", RowBox[{"(", RowBox[{"1", "+", RowBox[{"Exp", "[", RowBox[{"9.", "*", RowBox[{"(", RowBox[{"1", "-", "t"}], ")"}]}], "]"}]}], ")"}]}], "/", RowBox[{"(", RowBox[{"1", "+", RowBox[{"Exp", "[", RowBox[{"9", "*", RowBox[{"(", RowBox[{"t", "-", "6."}], ")"}]}], "]"}]}], ")"}]}], "+", "20"}]}], ";", RowBox[{ RowBox[{"b1", "[", "t_", "]"}], ":=", RowBox[{"Log", "[", RowBox[{"1.", "+", RowBox[{"Exp", "[", RowBox[{"k", "*", RowBox[{"(", RowBox[{ RowBox[{"T1", "[", "t", "]"}], "-", "Tc"}], ")"}]}], "]"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"bInt1", "[", "t_", "]"}], ":=", RowBox[{"NIntegrate", "[", RowBox[{ RowBox[{"-", RowBox[{"b1", "[", "x", "]"}]}], ",", RowBox[{"{", RowBox[{"x", ",", "0", ",", "t"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"bInt1", "[", "tFinal", "]"}]}]}], "]"}]}]], "Input", CellChangeTimes->{{3.393339887526169*^9, 3.3933398878273153`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"LS2", "[", RowBox[{ RowBox[{"k_", "?", "NumericQ"}], ",", RowBox[{"Tc_", "?", "NumericQ"}], ",", RowBox[{"tFinal_", "?", "NumericQ"}]}], "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{"b2", ",", "t", ",", "T2", ",", "x"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"T2", "[", "t_", "]"}], ":=", RowBox[{ RowBox[{ RowBox[{"133", "/", RowBox[{"(", RowBox[{"1", "+", RowBox[{"Exp", "[", RowBox[{"9.", "*", RowBox[{"(", RowBox[{"1", "-", "t"}], ")"}]}], "]"}]}], ")"}]}], "/", RowBox[{"(", RowBox[{"1", "+", RowBox[{"Exp", "[", RowBox[{"9", "*", RowBox[{"(", RowBox[{"t", "-", "5."}], ")"}]}], "]"}]}], ")"}]}], "+", "20"}]}], ";", RowBox[{ RowBox[{"b2", "[", "t_", "]"}], ":=", RowBox[{"Log", "[", RowBox[{"1.", "+", RowBox[{"Exp", "[", RowBox[{"k", "*", RowBox[{"(", RowBox[{ RowBox[{"T2", "[", "t", "]"}], "-", "Tc"}], ")"}]}], "]"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"bInt2", "[", "t_", "]"}], ":=", RowBox[{"NIntegrate", "[", RowBox[{ RowBox[{"-", RowBox[{"b2", "[", "x", "]"}]}], ",", RowBox[{"{", RowBox[{"x", ",", "0", ",", "t"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"bInt2", "[", "tFinal", "]"}]}]}], "]"}]}]], "Input", CellChangeTimes->{{3.393339891562159*^9, 3.3933398917635098`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"LS3", "[", RowBox[{ RowBox[{"k_", "?", "NumericQ"}], ",", RowBox[{"Tc_", "?", "NumericQ"}], ",", RowBox[{"tFinal_", "?", "NumericQ"}]}], "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{"b3", ",", "t", ",", "T3", ",", "x"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"T3", "[", "t_", "]"}], ":=", RowBox[{ RowBox[{ RowBox[{"160", "/", RowBox[{"(", RowBox[{"1", "+", RowBox[{"Exp", "[", RowBox[{"9.", "*", RowBox[{"(", RowBox[{"1", "-", "t"}], ")"}]}], "]"}]}], ")"}]}], "/", RowBox[{"(", RowBox[{"1", "+", RowBox[{"Exp", "[", RowBox[{"9", "*", RowBox[{"(", RowBox[{"t", "-", "4.5"}], ")"}]}], "]"}]}], ")"}]}], "+", "20"}]}], ";", RowBox[{ RowBox[{"b3", "[", "t_", "]"}], ":=", RowBox[{"Log", "[", RowBox[{"1.", "+", RowBox[{"Exp", "[", RowBox[{"k", "*", RowBox[{"(", RowBox[{ RowBox[{"T3", "[", "t", "]"}], "-", "Tc"}], ")"}]}], "]"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"bInt3", "[", "t_", "]"}], ":=", RowBox[{"NIntegrate", "[", RowBox[{ RowBox[{"-", RowBox[{"b3", "[", "x", "]"}]}], ",", RowBox[{"{", RowBox[{"x", ",", "0", ",", "t"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"bInt3", "[", "tFinal", "]"}]}]}], "]"}]}]], "Input", CellChangeTimes->{{3.3933398949497232`*^9, 3.393339895148252*^9}}] }, Open ]], Cell[CellGroupData[{ Cell["\<\ Compute and save the three survival values at each of the three final times \ (tCycle1, tCycle2 & tCycle3) using the original values for the parameters k \ and Tc.\ \>", "Text", CellChangeTimes->{ 3.393252270846343*^9, {3.393252527849209*^9, 3.39325252851716*^9}, { 3.39333009959886*^9, 3.393330105871035*^9}, {3.3933301622362003`*^9, 3.393330197332992*^9}}, TextAlignment->Left, FontSize->14], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"lstCycle1", "=", RowBox[{"LS1", "[", RowBox[{"k", ",", "Tc", ",", "tCycle1"}], "]"}]}], ";", RowBox[{"Print", "[", RowBox[{"\"\\"", ",", "lstCycle1"}], "]"}]}]], "Input",\ CellChangeTimes->{{3.393330222143549*^9, 3.393330243018742*^9}}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"LS1[k,Tc,tCycle1]=\"\>", "\[InvisibleSpace]", RowBox[{"-", "5.323457255213671`"}]}], SequenceForm["LS1[k,Tc,tCycle1]=", -5.323457255213671], Editable->False]], "Print", CellChangeTimes->{3.3933305781739807`*^9, 3.393330838764044*^9, 3.393331091953264*^9, 3.393339337178993*^9, 3.393339926724555*^9, 3.393340379717238*^9, 3.393341226585587*^9, 3.393769028117043*^9, 3.393844948182076*^9, 3.3941917311577673`*^9, 3.394202734504816*^9, 3.394207934480723*^9, 3.394965368614175*^9, 3.3949668336188793`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"lstCycle2", "=", RowBox[{"LS2", "[", RowBox[{"k", ",", "Tc", ",", "tCycle2"}], "]"}]}], ";", RowBox[{"Print", "[", RowBox[{"\"\\"", ",", "lstCycle2"}], "]"}]}]], "Input",\ CellChangeTimes->{{3.3933302478016157`*^9, 3.3933302690761337`*^9}}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"LS2[k,Tc,tCycle2]=\"\>", "\[InvisibleSpace]", RowBox[{"-", "6.533848330579998`"}]}], SequenceForm["LS2[k,Tc,tCycle2]=", -6.533848330579998], Editable->False]], "Print", CellChangeTimes->{3.393330585178144*^9, 3.39333083886591*^9, 3.393331092054915*^9, 3.393339337298172*^9, 3.39333992683287*^9, 3.393340379849668*^9, 3.393341226695669*^9, 3.3937690284154043`*^9, 3.393844948286063*^9, 3.394191731279695*^9, 3.39420273459414*^9, 3.394207934744274*^9, 3.394965368711134*^9, 3.394966833708251*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"lstCycle3", "=", RowBox[{"LS3", "[", RowBox[{"k", ",", "Tc", ",", "tCycle3"}], "]"}]}], ";", RowBox[{"Print", "[", RowBox[{"\"\\"", ",", "lstCycle3"}], "]"}]}]], "Input",\ CellChangeTimes->{ 3.3931662602221727`*^9, {3.393330273658127*^9, 3.393330287048992*^9}}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"LS3[k,Tc,tCycle3]=\"\>", "\[InvisibleSpace]", RowBox[{"-", "8.430087595093248`"}]}], SequenceForm["LS3[k,Tc,tCycle3]=", -8.430087595093248], Editable->False]], "Print", CellChangeTimes->{3.393330590432136*^9, 3.393330838959113*^9, 3.393331092155212*^9, 3.39333933757596*^9, 3.393339926934671*^9, 3.393340380150382*^9, 3.393341226819251*^9, 3.393769028524839*^9, 3.3938449484064493`*^9, 3.39419173158877*^9, 3.394202734687237*^9, 3.39420793491243*^9, 3.3949653688137627`*^9, 3.394966833790856*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["\<\ Find k and Tc using LS1[tCycle1] and LS3[tCycle3]. NOTICE THAT IN EACH \ COMBINATION THE PROGRAM FINDS THE ORIGINAL GENERATION PARAMETER VALUES \ k=0.035 and Tc=110.\ \>", "Text", CellChangeTimes->{{3.393166362624907*^9, 3.393166454118064*^9}, { 3.3931671615911922`*^9, 3.393167165924527*^9}, {3.3932422183664618`*^9, 3.3932422232985697`*^9}, {3.393242314952285*^9, 3.393242315424522*^9}, { 3.393242459395752*^9, 3.393242480010466*^9}, {3.393243035099559*^9, 3.393243039247012*^9}, {3.3932525627153597`*^9, 3.393252563470331*^9}, 3.3932534949476748`*^9, {3.393340581694598*^9, 3.393340591005661*^9}}, TextAlignment->Left, FontSize->14], Cell[BoxData[ RowBox[{"Clear", "[", RowBox[{"k", ",", "Tc"}], "]"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"frResult", "=", RowBox[{"FindRoot", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"LS1", "[", RowBox[{"k", ",", "Tc", ",", "tCycle1"}], "]"}], "==", "lstCycle1"}], ",", RowBox[{ RowBox[{"LS3", "[", RowBox[{"k", ",", "Tc", ",", "tCycle3"}], "]"}], "\[Equal]", "lstCycle3"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"k", ",", ".01", ",", ".015"}], "}"}], ",", RowBox[{"{", RowBox[{"Tc", ",", "90.", ",", "95."}], "}"}]}], "}"}], ",", RowBox[{"MaxIterations", "\[Rule]", "2000"}]}], "]"}]}]], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"k", "\[Rule]", "0.03500000000000068`"}], ",", RowBox[{"Tc", "\[Rule]", "110.00000000000229`"}]}], "}"}]], "Output", CellChangeTimes->{3.393166623615265*^9, 3.393243787630814*^9, 3.3932496426723623`*^9, 3.393250869654623*^9, 3.393252669594225*^9, 3.39333059845035*^9, 3.393330839657639*^9, 3.393331092849626*^9, 3.393339338278491*^9, 3.3933399279946938`*^9, 3.3933403809072847`*^9, 3.3933412278842287`*^9, 3.393769029195719*^9, 3.3938449491222878`*^9, 3.3941917323075933`*^9, 3.394202735756583*^9, 3.39420793571881*^9, 3.3949653698083057`*^9, 3.394966834506035*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Find k and Tc using LS2(tCycle2) and LS3(tCycle3).", "Text", CellChangeTimes->{{3.393166362624907*^9, 3.393166454118064*^9}, { 3.393166526203062*^9, 3.3931665315765142`*^9}, 3.393166587298316*^9, 3.393167175604994*^9, {3.393242243623486*^9, 3.393242245954938*^9}}, TextAlignment->Left, FontSize->14], Cell[BoxData[ RowBox[{"Clear", "[", RowBox[{"k", ",", "Tc"}], "]"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"frResult", "=", RowBox[{"FindRoot", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"LS2", "[", RowBox[{"k", ",", "Tc", ",", "tCycle2"}], "]"}], "==", "lstCycle2"}], ",", RowBox[{ RowBox[{"LS3", "[", RowBox[{"k", ",", "Tc", ",", "tCycle3"}], "]"}], "\[Equal]", "lstCycle3"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"k", ",", ".01", ",", ".015"}], "}"}], ",", RowBox[{"{", RowBox[{"Tc", ",", "90.", ",", "95."}], "}"}]}], "}"}], ",", RowBox[{"MaxIterations", "\[Rule]", "2000"}]}], "]"}]}]], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"k", "\[Rule]", "0.0350000000002601`"}], ",", RowBox[{"Tc", "\[Rule]", "110.00000000240489`"}]}], "}"}]], "Output", CellChangeTimes->{3.393166625337246*^9, 3.393243788467682*^9, 3.393249643287869*^9, 3.393250870446741*^9, 3.3932526703281813`*^9, 3.3933306041872807`*^9, 3.393330840786152*^9, 3.39333109355436*^9, 3.393339338953945*^9, 3.393339928725037*^9, 3.393340381715424*^9, 3.393341228685154*^9, 3.393769029870781*^9, 3.3938449497831993`*^9, 3.39419173298339*^9, 3.394202736340352*^9, 3.394207936340229*^9, 3.394965370477409*^9, 3.394966835131833*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Find k and Tc using LS1(tCycle1) and LS2(tCycle2).", "Text", CellChangeTimes->{{3.393166362624907*^9, 3.393166454118064*^9}, { 3.393166577468137*^9, 3.393166582658086*^9}, {3.393167193317575*^9, 3.393167196101915*^9}, {3.393242265526527*^9, 3.393242265924616*^9}}, TextAlignment->Left, FontSize->14], Cell[BoxData[ RowBox[{"Clear", "[", RowBox[{"k", ",", "Tc"}], "]"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"frResult", "=", RowBox[{"FindRoot", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"LS1", "[", RowBox[{"k", ",", "Tc", ",", "tCycle1"}], "]"}], "==", "lstCycle1"}], ",", RowBox[{ RowBox[{"LS2", "[", RowBox[{"k", ",", "Tc", ",", "tCycle2"}], "]"}], "\[Equal]", "lstCycle2"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"k", ",", ".01", ",", ".02"}], "}"}], ",", RowBox[{"{", RowBox[{"Tc", ",", "90.", ",", "95."}], "}"}]}], "}"}], ",", RowBox[{"MaxIterations", "\[Rule]", "2000"}]}], "]"}]}]], "Input", CellChangeTimes->{3.3931666688903503`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"k", "\[Rule]", "0.03499999999998977`"}], ",", RowBox[{"Tc", "\[Rule]", "109.9999999999968`"}]}], "}"}]], "Output", CellChangeTimes->{3.393166629907181*^9, 3.393166676962244*^9, 3.393243789196649*^9, 3.3932496439824247`*^9, 3.393250871183515*^9, 3.393252671328047*^9, 3.39333060835049*^9, 3.393330841984303*^9, 3.393331094364418*^9, 3.3933393396663*^9, 3.393339929459578*^9, 3.3933403824545193`*^9, 3.393341229414462*^9, 3.393769030545883*^9, 3.393844950498653*^9, 3.394191733753697*^9, 3.394202737125207*^9, 3.3942079370193653`*^9, 3.394965371226275*^9, 3.394966835831992*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Plot the three temperature profile functions.", "Text", CellChangeTimes->{{3.393166362624907*^9, 3.393166454118064*^9}, { 3.393166577468137*^9, 3.393166582658086*^9}, {3.393166725191359*^9, 3.393166781654152*^9}, {3.393166882404248*^9, 3.3931668998421392`*^9}, 3.393252230849867*^9, {3.393330330715049*^9, 3.393330332134222*^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.", ",", "200."}], "}"}]}], "}"}]}], ",", 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]", "\"\\""}], ",", RowBox[{"Frame", "\[Rule]", "True"}], ",", RowBox[{"RotateLabel", "\[Rule]", "True"}], ",", RowBox[{"BaseStyle", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"FontFamily", "\[Rule]", "\"\\""}], ",", RowBox[{"FontWeight", "\[Rule]", "\"\\""}], ",", RowBox[{"FontSize", "\[Rule]", "14"}]}], "}"}]}], ",", RowBox[{"FrameLabel", "\[Rule]", RowBox[{"{", RowBox[{"\"\<\>\"", ",", "\"\\""}], "}"}]}], ",", RowBox[{"FrameTicks", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "2", ",", "4", ",", "6", ",", "8", ",", "10"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "50", ",", "100", ",", "150", ",", "200"}], "}"}], ",", RowBox[{"{", "}"}], ",", RowBox[{"{", "}"}]}], "}"}]}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.3932494781597433`*^9, 3.39324949828539*^9}, { 3.393249561255631*^9, 3.393249658020599*^9}, 3.3932508164105663`*^9, 3.393250877249597*^9, {3.3933308823753023`*^9, 3.39333093160261*^9}, 3.3933312008690434`*^9, {3.3933312326698236`*^9, 3.393331236179345*^9}, 3.3933312752089977`*^9, 3.3937685217923098`*^9, {3.394191561909532*^9, 3.39419161012357*^9}, 3.394191697830122*^9}] }, Open ]], Cell[CellGroupData[{ Cell["\<\ ListPlot the temperature at 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.39333063583978*^9, 3.393330647794845*^9}, {3.3933311342538*^9, 3.393331141059312*^9}}, TextAlignment->Left, FontSize->14], Cell[BoxData[ RowBox[{ RowBox[{"tempPtsPlot", "=", RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"tCycle3", ",", "tCycle3T1"}], "}"}], ",", RowBox[{"{", RowBox[{"tCycle3", ",", "tCycle3T2"}], "}"}], ",", RowBox[{"{", RowBox[{"tCycle3", ",", "tCycle3T3"}], "}"}]}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0.", ",", "10."}], "}"}], ",", RowBox[{"{", RowBox[{"0.", ",", "200."}], "}"}]}], "}"}]}], ",", RowBox[{"AspectRatio", "\[Rule]", "1"}], ",", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"Directive", "[", RowBox[{ RowBox[{"PointSize", "[", "Large", "]"}], ",", "Black"}], "]"}]}], ",", RowBox[{"PlotLabel", "\[Rule]", "\"\\""}], ",", RowBox[{"Frame", "\[Rule]", "True"}], ",", RowBox[{"RotateLabel", "\[Rule]", "True"}], ",", RowBox[{"BaseStyle", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"FontFamily", "\[Rule]", "\"\\""}], ",", RowBox[{"FontWeight", "\[Rule]", "\"\\""}], ",", RowBox[{"FontSize", "\[Rule]", "14"}]}], "}"}]}], ",", RowBox[{"FrameLabel", "\[Rule]", RowBox[{"{", RowBox[{ "\"\