(* 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[ 71296, 2316] NotebookOptionsPosition[ 64203, 2093] NotebookOutlinePosition[ 64544, 2108] CellTagsIndexPosition[ 64501, 2105] WindowFrame->Normal ContainsDynamic->False*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[TextData[{ StyleBox["Calculation of High-Pressure \nChemical Equilibrium", "Subsubtitle", FontSize->18, FontVariations->{"CompatibilityType"->0}], StyleBox["\n", FontWeight->"Plain", FontVariations->{"CompatibilityType"->0}], StyleBox["Author's Data", FontSize->14, FontWeight->"Bold"], StyleBox[": Housam Binous", FontSize->14], StyleBox["\n", TextAlignment->Center, FontFamily->"MS Shell Dlg", FontSize->8.5, Background->RGBColor[0.605478, 0.996109, 0.605478]], StyleBox["Department of Chemical Engineering\nNational Institute of Applied \ Sciences and Technology\nTunis, TUNISIA\nEmail: binoushousam@yahoo.com ", FontSize->14, FontWeight->"Plain"] }], "Title", TextAlignment->Center, Background->RGBColor[0.605478, 0.996109, 0.605478]], Cell[TextData[{ StyleBox["Acknowledgement :", FontWeight->"Bold"], "\nThis problem was presented in the excellent book by ", StyleBox["S. I. Sandler, Chemical and Engineering Thermodynamics, 3rd \ Edition, Wiley 1999 (Illustration 9.1-8 page 654).", FontWeight->"Bold"] }], "Subsubtitle", Background->RGBColor[0.773449, 0.996109, 0.996109]], Cell[BoxData[ RowBox[{"Off", "[", RowBox[{"General", "::", "\"\\""}], "]"}]], "Input"], Cell[CellGroupData[{ Cell[TextData[{ "Equilibrium constant versus temperature for the reaction:\n\n\t\t\t", Cell[BoxData[ FormBox[ RowBox[{ RowBox[{ RowBox[{ FractionBox["1", "2"], SubscriptBox["N", "2"]}], "+", RowBox[{ FractionBox["3", "2"], SubscriptBox["H", "2"]}]}], "=", SubscriptBox[ StyleBox["NH", FontSlant->"Italic"], "3"]}], TraditionalForm]]] }], "Subsubtitle", Background->RGBColor[1, 1, 0.658824]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Ka298", "=", RowBox[{"Exp", "[", RowBox[{"16.5", " ", RowBox[{"1000", "/", RowBox[{"(", RowBox[{"8.314", " ", "298.15"}], ")"}]}]}], "]"}]}]], "Input"], Cell[BoxData["777.7424108820484`"], "Output", CellChangeTimes->{3.4020637126000752`*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"\[CapitalDelta]H", "[", "T_", "]"}], ":=", RowBox[{ RowBox[{"-", "46100"}], "+", RowBox[{ RowBox[{"(", RowBox[{"T", "-", "298.15"}], ")"}], " ", "\[CapitalDelta]a"}], "+", RowBox[{ RowBox[{"\[CapitalDelta]b", "/", "2"}], " ", RowBox[{"(", RowBox[{ RowBox[{"T", "^", "2"}], "-", RowBox[{"298.15", "^", "2"}]}], ")"}]}], "+", RowBox[{ RowBox[{"\[CapitalDelta]c", "/", "3"}], " ", RowBox[{"(", RowBox[{ RowBox[{"T", "^", "3"}], "-", RowBox[{"298.15", "^", "3"}]}], ")"}]}], "+", RowBox[{ RowBox[{"\[CapitalDelta]d", "/", "4"}], " ", RowBox[{"(", RowBox[{ RowBox[{"T", "^", "4"}], "-", RowBox[{"298.15", "^", "4"}]}], ")"}]}]}]}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"\[CapitalDelta]a", "=", RowBox[{"24.619", "-", RowBox[{"0.5", " ", "27.318"}], "-", RowBox[{"1.5", " ", "26.879"}]}]}]], "Input"], Cell[BoxData[ RowBox[{"-", "29.3585`"}]], "Output", CellChangeTimes->{3.402063714122264*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"\[CapitalDelta]b", "=", RowBox[{ RowBox[{"(", RowBox[{"3.75", "-", RowBox[{"0.5", " ", RowBox[{"(", "0.623", ")"}]}], "-", RowBox[{"1.5", " ", RowBox[{"(", "0.435", ")"}]}]}], ")"}], " ", RowBox[{"10", "^", RowBox[{"-", "2"}]}]}]}]], "Input"], Cell[BoxData["0.02786`"], "Output", CellChangeTimes->{3.4020637148933725`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"\[CapitalDelta]c", "=", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "0.138"}], "+", RowBox[{"0.5", " ", RowBox[{"(", "0.095", ")"}]}], "+", RowBox[{"1.5", " ", RowBox[{"(", "0.033", ")"}]}]}], ")"}], " ", RowBox[{"10", "^", RowBox[{"-", "5"}]}]}]}]], "Input"], Cell[BoxData[ RowBox[{"-", "4.100000000000001`*^-7"}]], "Output", CellChangeTimes->{3.402063715654467*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"\[CapitalDelta]d", "=", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"0.5", " ", RowBox[{"(", "2.871", ")"}]}], "+", RowBox[{"1.5", " ", RowBox[{"(", "0.87", ")"}]}]}], ")"}], " ", RowBox[{"10", "^", RowBox[{"-", "9"}]}]}]}]], "Input"], Cell[BoxData["2.7405`*^-9"], "Output", CellChangeTimes->{3.402063716445605*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"Ka", "[", "T_", "]"}], ":=", RowBox[{"Ka298", " ", RowBox[{"Exp", "[", RowBox[{ RowBox[{ RowBox[{"\[CapitalDelta]a", "/", "8.314"}], " ", RowBox[{"Log", "[", RowBox[{"T", "/", "298.15"}], "]"}]}], "+", RowBox[{ RowBox[{"\[CapitalDelta]b", "/", RowBox[{"(", RowBox[{"2", " ", "8.314"}], ")"}]}], " ", RowBox[{"(", RowBox[{"T", "-", "298.15"}], ")"}]}], "+", RowBox[{ RowBox[{"\[CapitalDelta]c", "/", RowBox[{"(", RowBox[{"6", " ", "8.314"}], ")"}]}], " ", RowBox[{"(", RowBox[{ RowBox[{"T", "^", "2"}], "-", RowBox[{"298.15", "^", "2"}]}], ")"}]}], "+", RowBox[{ RowBox[{"\[CapitalDelta]d", "/", RowBox[{"(", RowBox[{"12", " ", "8.314"}], ")"}]}], " ", RowBox[{"(", RowBox[{ RowBox[{"T", "^", "3"}], "-", RowBox[{"298.15", "^", "3"}]}], ")"}]}], "+", RowBox[{ RowBox[{"1", "/", "8.314"}], " ", RowBox[{"(", RowBox[{"46100", "+", RowBox[{ RowBox[{"(", "298.15", ")"}], " ", "\[CapitalDelta]a"}], "+", RowBox[{ RowBox[{"\[CapitalDelta]b", "/", "2"}], " ", RowBox[{"(", RowBox[{"298.15", "^", "2"}], ")"}]}], "+", RowBox[{ RowBox[{"\[CapitalDelta]c", "/", "3"}], " ", RowBox[{"(", RowBox[{"298.15", "^", "3"}], ")"}]}], "+", RowBox[{ RowBox[{"\[CapitalDelta]d", "/", "4"}], " ", RowBox[{"(", RowBox[{"298.15", "^", "4"}], ")"}]}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"1", "/", "T"}], "-", RowBox[{"1", "/", "298.15"}]}], ")"}]}]}], "]"}]}]}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Ka", "[", "500", "]"}]], "Input"], Cell[BoxData["0.328260263916996`"], "Output", CellChangeTimes->{3.402063718388398*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Ka", "[", "600", "]"}]], "Input"], Cell[BoxData["0.04347004619110203`"], "Output", CellChangeTimes->{3.402063719229608*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Ka", "[", "700", "]"}]], "Input"], Cell[BoxData["0.009901212453597258`"], "Output", CellChangeTimes->{3.40206372003076*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Ka", "[", "800", "]"}]], "Input"], Cell[BoxData["0.0032004900349465112`"], "Output", CellChangeTimes->{3.4020637208719697`*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Reduced Pressure and Temperature", "Subsubtitle", Background->RGBColor[1, 1, 0.658824]], Cell[BoxData[ RowBox[{ RowBox[{"Tre", "[", "i_", "]"}], ":=", RowBox[{"T", "/", RowBox[{"Tc", "[", "i", "]"}]}]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"Pre", "[", "i_", "]"}], ":=", RowBox[{"P", "/", RowBox[{"Pc", "[", "i", "]"}]}]}]], "Input"] }, Open ]], Cell[CellGroupData[{ Cell[TextData[StyleBox["Fugacity coefficients using the Peng Robinson \ Equation of State. Binary interaction parameters for the Peng Robinson EOS \ were taken from HYSYS 3.2.", "Subsubtitle", FontSize->18, FontVariations->{"CompatibilityType"->0}]], "Subsubtitle", Background->RGBColor[1, 1, 0.658824]], Cell[BoxData[ RowBox[{ RowBox[{"ap", "[", "i_", "]"}], ":=", RowBox[{"0.45724", " ", RowBox[{"a", "[", "i", "]"}], " ", RowBox[{"R", "^", "2"}], " ", RowBox[{ RowBox[{ RowBox[{"Tc", "[", "i", "]"}], "^", "2"}], "/", RowBox[{"Pc", "[", "i", "]"}]}], " "}]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"Ap", "[", "i_", "]"}], ":=", RowBox[{ RowBox[{"ap", "[", "i", "]"}], " ", RowBox[{"P", "/", RowBox[{ RowBox[{"(", RowBox[{"R", " ", "T"}], ")"}], "^", "2"}]}]}]}], ";"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"bp", "[", "i_", "]"}], ":=", RowBox[{"0.07780", " ", "R", " ", RowBox[{ RowBox[{"Tc", "[", "i", "]"}], "/", RowBox[{"Pc", "[", "i", "]"}]}]}]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"Bp", "[", "i_", "]"}], ":=", RowBox[{"P", " ", RowBox[{ RowBox[{"bp", "[", "i", "]"}], "/", RowBox[{"(", RowBox[{"R", " ", "T"}], ")"}]}]}]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{ SubscriptBox["\[Phi]", "v"], "[", "i_", "]"}], ":=", RowBox[{"Exp", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ SubscriptBox["Z", "v"], "-", "1"}], ")"}], " ", RowBox[{ RowBox[{"Bp", "[", "i", "]"}], "/", SubscriptBox["B", "v"]}]}], "-", RowBox[{"Log", "[", RowBox[{ SubscriptBox["Z", "v"], "-", SubscriptBox["B", "v"]}], "]"}], "-", RowBox[{ RowBox[{ SubscriptBox["A", "v"], "/", RowBox[{"(", RowBox[{"2", " ", SqrtBox["2"], SubscriptBox["B", "v"]}], ")"}]}], " ", RowBox[{"(", RowBox[{ RowBox[{"2", " ", RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"j", "=", "1"}], "3"], RowBox[{ RowBox[{"y", "[", "j", "]"}], " ", RowBox[{ RowBox[{"Ab", "[", RowBox[{"i", ",", "j"}], "]"}], "/", SubscriptBox["A", "v"]}]}]}]}], "-", RowBox[{ RowBox[{"Bp", "[", "i", "]"}], "/", SubscriptBox["B", "v"]}]}], ")"}], " ", RowBox[{"Log", "[", RowBox[{ RowBox[{"(", RowBox[{ SubscriptBox["Z", "v"], "+", RowBox[{ RowBox[{"(", RowBox[{"1", "+", SqrtBox["2"]}], ")"}], SubscriptBox["B", "v"]}]}], ")"}], "/", RowBox[{"(", RowBox[{ SubscriptBox["Z", "v"], "+", RowBox[{ RowBox[{"(", RowBox[{"1", "-", SqrtBox["2"]}], ")"}], SubscriptBox["B", "v"]}]}], ")"}]}], "]"}]}]}], "]"}]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"a", "[", "i_", "]"}], ":=", RowBox[{ RowBox[{"(", RowBox[{"1", "+", RowBox[{ RowBox[{"m", "[", "i", "]"}], " ", RowBox[{"(", RowBox[{"1", "-", RowBox[{ RowBox[{"Tre", "[", "i", "]"}], "^", "0.5"}]}], ")"}]}]}], ")"}], "^", "2"}]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"m", "[", "i_", "]"}], ":=", RowBox[{"0.37464", "+", RowBox[{"1.54226", " ", RowBox[{"\[Omega]", "[", "i", "]"}]}], "-", RowBox[{"0.26992", " ", RowBox[{ RowBox[{"\[Omega]", "[", "i", "]"}], "^", "2"}]}]}]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{ SubscriptBox["A", "v"], "=", RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "3"], RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"j", "=", "1"}], "3"], RowBox[{ RowBox[{"y", "[", "i", "]"}], " ", RowBox[{"y", "[", "j", "]"}], " ", RowBox[{"Ab", "[", RowBox[{"i", ",", "j"}], "]"}]}]}]}]}], ";"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"Ab", "[", RowBox[{"i_", ",", "j_"}], "]"}], "=", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"Ap", "[", "i", "]"}], " ", RowBox[{"Ap", "[", "j", "]"}]}], ")"}], "^", "0.5"}], RowBox[{"(", RowBox[{"1", "-", RowBox[{"k", "[", RowBox[{"i", ",", "j"}], "]"}]}], ")"}]}]}], ";"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{ SubscriptBox["B", "v"], "=", RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "3"], RowBox[{ RowBox[{"y", "[", "i", "]"}], " ", RowBox[{"Bp", "[", "i", "]"}]}]}]}], ";"}]], "Input"], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"k", "[", RowBox[{"1", ",", "1"}], "]"}], "=", "0"}], ";", RowBox[{ RowBox[{"k", "[", RowBox[{"2", ",", "2"}], "]"}], "=", "0"}], ";", RowBox[{ RowBox[{"k", "[", RowBox[{"3", ",", "3"}], "]"}], "=", "0"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"k", "[", RowBox[{"2", ",", "1"}], "]"}], "=", RowBox[{"k", "[", RowBox[{"1", ",", "2"}], "]"}]}], ";", RowBox[{ RowBox[{"k", "[", RowBox[{"3", ",", "1"}], "]"}], "=", RowBox[{"k", "[", RowBox[{"1", ",", "3"}], "]"}]}], ";", RowBox[{ RowBox[{"k", "[", RowBox[{"2", ",", "3"}], "]"}], "=", RowBox[{"k", "[", RowBox[{"3", ",", "2"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"k", "[", RowBox[{"1", ",", "2"}], "]"}], "=", RowBox[{"-", "0.036"}]}], ";", RowBox[{ RowBox[{"k", "[", RowBox[{"1", ",", "3"}], "]"}], "=", "0"}], ";", RowBox[{ RowBox[{"k", "[", RowBox[{"3", ",", "2"}], "]"}], "=", "0.222"}], ";"}]}], "Input"] }, Open ]], Cell[CellGroupData[{ Cell["Gas constant", "Subsubtitle", Background->RGBColor[1, 1, 0.658824]], Cell[BoxData[ RowBox[{ RowBox[{"R", "=", "1.987"}], ";"}]], "Input"] }, Open ]], Cell[CellGroupData[{ Cell["\<\ Critical Pressure and Temperature and acentric factor for Hydrogen, Nitrogen \ and Ammonia\ \>", "Subsubtitle", Background->RGBColor[1, 1, 0.658824]], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"Pc", "[", "1", "]"}], "=", "190.8"}], ";", RowBox[{ RowBox[{"Tc", "[", "1", "]"}], "=", "59.7"}], ";", RowBox[{ RowBox[{"\[Omega]", "[", "1", "]"}], "=", "0.0"}], ";"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"Pc", "[", "2", "]"}], "=", "492.9"}], ";", RowBox[{ RowBox[{"Tc", "[", "2", "]"}], "=", "227.3"}], ";", RowBox[{ RowBox[{"\[Omega]", "[", "2", "]"}], "=", "0.0206"}], ";"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"Pc", "[", "3", "]"}], "=", "1653.7"}], ";", RowBox[{ RowBox[{"Tc", "[", "3", "]"}], "=", "730.2"}], ";", RowBox[{ RowBox[{"\[Omega]", "[", "3", "]"}], "=", "0.2582"}], ";"}]], "Input"] }, Open ]], Cell[CellGroupData[{ Cell[TextData[{ "The correction factor ", Cell[BoxData[ FormBox[ SubscriptBox["K", "v"], TraditionalForm]]], " and the extent of reaction X versus P for T=500K or 900R" }], "Subsubtitle", Background->RGBColor[1, 1, 0.658824]], Cell[BoxData[ RowBox[{ RowBox[{"T", "=", "900"}], ";"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"P", "=", "0.001"}], ";", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "<", "17"}], ",", RowBox[{"{", RowBox[{ RowBox[{"sol900", "=", RowBox[{"FindRoot", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{"Ka", "[", "500", "]"}], RowBox[{ RowBox[{"(", RowBox[{"0.5", " ", RowBox[{"(", RowBox[{"1", "-", "X"}], ")"}]}], ")"}], "^", "0.5"}], " ", RowBox[{ RowBox[{"(", RowBox[{"1.5", " ", RowBox[{"(", RowBox[{"1", "-", "X"}], ")"}]}], ")"}], "^", "1.5"}], " ", RowBox[{ RowBox[{ SubscriptBox["\[Phi]", "v"], "[", "1", "]"}], "^", "1.5"}], " ", RowBox[{ RowBox[{ SubscriptBox["\[Phi]", "v"], "[", "2", "]"}], "^", "0.5"}]}], "\[Equal]", " ", RowBox[{"X", " ", RowBox[{ RowBox[{"(", RowBox[{"2", "-", "X"}], ")"}], "/", " ", RowBox[{"(", RowBox[{"P", "/", "14.69"}], ")"}]}], " ", RowBox[{ SubscriptBox["\[Phi]", "v"], "[", "3", "]"}]}]}], ",", RowBox[{ RowBox[{ RowBox[{"y", "[", "3", "]"}], RowBox[{"(", RowBox[{"2", "-", "X"}], ")"}]}], "\[Equal]", " ", "X"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"y", "[", "2", "]"}], RowBox[{"(", RowBox[{"2", "-", "X"}], ")"}]}], "\[Equal]", " ", RowBox[{"0.5", RowBox[{"(", RowBox[{"1", "-", "X"}], ")"}]}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"y", "[", "1", "]"}], RowBox[{"(", RowBox[{"2", "-", "X"}], ")"}]}], "\[Equal]", " ", RowBox[{"1.5", " ", RowBox[{"(", RowBox[{"1", "-", "X"}], ")"}]}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["Z", "v"], "^", "3"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", SubscriptBox["B", "v"]}], ")"}], RowBox[{ SubscriptBox["Z", "v"], "^", "2"}]}], "+", RowBox[{ SubscriptBox["Z", "v"], " ", RowBox[{"(", RowBox[{ SubscriptBox["A", "v"], "-", RowBox[{"3", RowBox[{ SubscriptBox["B", "v"], "^", "2"}]}], "-", RowBox[{"2", SubscriptBox["B", "v"]}]}], ")"}]}], "-", RowBox[{ SubscriptBox["A", "v"], " ", SubscriptBox["B", "v"]}], "+", RowBox[{ SubscriptBox["B", "v"], "^", "2"}], "+", RowBox[{ SubscriptBox["B", "v"], "^", "3"}]}], "\[Equal]", "0"}]}], "}"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{ RowBox[{"y", "[", "2", "]"}], ",", "0.2"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"y", "[", "1", "]"}], ",", "0.1"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"y", "[", "3", "]"}], ",", "0.4"}], "}"}], ",", RowBox[{"{", RowBox[{"X", ",", "0.95"}], "}"}], ",", RowBox[{"{", RowBox[{ SubscriptBox["Z", "v"], ",", "0.95"}], "}"}], ",", RowBox[{"MaxIterations", "\[Rule]", " ", "10000"}]}], "]"}]}], ",", RowBox[{ RowBox[{"Pp", "[", "i", "]"}], "=", "P"}], ",", RowBox[{ RowBox[{"Kv", "[", "i", "]"}], "=", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["\[Phi]", "v"], "[", "3", "]"}], "/", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["\[Phi]", "v"], "[", "1", "]"}], "^", "1.5"}], " ", RowBox[{ RowBox[{ SubscriptBox["\[Phi]", "v"], "[", "2", "]"}], "^", "0.5"}]}], ")"}]}], "/.", "sol900"}]}], ",", RowBox[{ RowBox[{"X", "[", "i", "]"}], "=", RowBox[{"X", "/.", "sol900"}]}], ",", RowBox[{"Print", "[", RowBox[{"i", ",", "\"\< \>\"", ",", "P"}], "]"}], ",", RowBox[{"P", "=", RowBox[{"P", "+", "1000"}]}], ",", RowBox[{"i", "=", RowBox[{"i", "+", "1"}]}]}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.402063765816597*^9, 3.402063780547779*^9}}], Cell[CellGroupData[{ Cell[BoxData[ InterpretationBox[ RowBox[{ "1", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "0.001`"}], SequenceForm[1, " ", 0.001], Editable->False]], "Print", CellChangeTimes->{3.4020637840728483`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "2", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1000.001`"}], SequenceForm[2, " ", 1000.001], Editable->False]], "Print", CellChangeTimes->{3.4020637842931647`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "3", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "2000.001`"}], SequenceForm[3, " ", 2000.001], Editable->False]], "Print", CellChangeTimes->{3.4020637845234957`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "4", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "3000.001`"}], SequenceForm[4, " ", 3000.001], Editable->False]], "Print", CellChangeTimes->{3.402063784753827*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "5", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "4000.001`"}], SequenceForm[5, " ", 4000.001], Editable->False]], "Print", CellChangeTimes->{3.402063784974144*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "6", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "5000.001`"}], SequenceForm[6, " ", 5000.001], Editable->False]], "Print", CellChangeTimes->{3.402063785224504*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "7", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "6000.001`"}], SequenceForm[7, " ", 6000.001], Editable->False]], "Print", CellChangeTimes->{3.4020637854348063`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "8", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "7000.001`"}], SequenceForm[8, " ", 7000.001], Editable->False]], "Print", CellChangeTimes->{3.4020637856651373`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "9", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "8000.001`"}], SequenceForm[9, " ", 8000.001], Editable->False]], "Print", CellChangeTimes->{3.4020637859054832`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "10", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "9000.001`"}], SequenceForm[10, " ", 9000.001], Editable->False]], "Print", CellChangeTimes->{3.4020637861358147`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "11", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "10000.001`"}], SequenceForm[11, " ", 10000.001], Editable->False]], "Print", CellChangeTimes->{3.4020637863361025`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "12", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "11000.001`"}], SequenceForm[12, " ", 11000.001], Editable->False]], "Print", CellChangeTimes->{3.4020637865564194`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "13", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "12000.001`"}], SequenceForm[13, " ", 12000.001], Editable->False]], "Print", CellChangeTimes->{3.4020637867466927`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "14", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "13000.001`"}], SequenceForm[14, " ", 13000.001], Editable->False]], "Print", CellChangeTimes->{3.4020637869469805`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "15", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "14000.001`"}], SequenceForm[15, " ", 14000.001], Editable->False]], "Print", CellChangeTimes->{3.402063787147269*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "16", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "15000.001`"}], SequenceForm[16, " ", 15000.001], Editable->False]], "Print", CellChangeTimes->{3.402063787347557*^9}] }, Open ]] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"tbl900Kv", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Pp", "[", "i", "]"}], ",", RowBox[{"Chop", "[", RowBox[{"Kv", "[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "16"}], "}"}]}], "]"}]}], ";"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"pltkv900", "=", RowBox[{"ListPlot", "[", RowBox[{"tbl900Kv", ",", RowBox[{"Joined", "\[Rule]", " ", "True"}], ",", " ", RowBox[{"PlotStyle", "\[Rule]", " ", RowBox[{"Hue", "[", "0.1", "]"}]}], ",", RowBox[{"AspectRatio", "\[Rule]", "1"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.402063793015707*^9, 3.402063805784067*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"tbl900X", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Pp", "[", "i", "]"}], ",", RowBox[{"Chop", "[", RowBox[{"X", "[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "16"}], "}"}]}], "]"}]}], ";"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"pltX900", "=", RowBox[{"ListPlot", "[", RowBox[{"tbl900X", ",", RowBox[{"Joined", "\[Rule]", " ", "True"}], ",", " ", RowBox[{"PlotStyle", "\[Rule]", " ", RowBox[{"Hue", "[", "0.1", "]"}]}], ",", RowBox[{"AspectRatio", "\[Rule]", "1"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.402063812593859*^9, 3.4020638187827587`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[TextData[{ "The correction factor ", Cell[BoxData[ FormBox[ SubscriptBox["K", "v"], TraditionalForm]]], " and the extent of reaction X versus P for T=600K or 1080R" }], "Subsubtitle", Background->RGBColor[1, 1, 0.658824]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"T", "=", RowBox[{"600", " ", "1.8"}]}]], "Input"], Cell[BoxData["1080.`"], "Output", CellChangeTimes->{3.4020638249916863`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"P", "=", "0.001"}], ";", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "<", "17"}], ",", RowBox[{"{", RowBox[{ RowBox[{"sol1080", "=", RowBox[{"FindRoot", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{"Ka", "[", "600", "]"}], RowBox[{ RowBox[{"(", RowBox[{"0.5", " ", RowBox[{"(", RowBox[{"1", "-", "X"}], ")"}]}], ")"}], "^", "0.5"}], " ", RowBox[{ RowBox[{"(", RowBox[{"1.5", " ", RowBox[{"(", RowBox[{"1", "-", "X"}], ")"}]}], ")"}], "^", "1.5"}], " ", RowBox[{ RowBox[{ SubscriptBox["\[Phi]", "v"], "[", "1", "]"}], "^", "1.5"}], " ", RowBox[{ RowBox[{ SubscriptBox["\[Phi]", "v"], "[", "2", "]"}], "^", "0.5"}]}], "\[Equal]", " ", RowBox[{"X", " ", RowBox[{ RowBox[{"(", RowBox[{"2", "-", "X"}], ")"}], "/", RowBox[{"(", RowBox[{"P", " ", "/", "14.69"}], ")"}]}], " ", RowBox[{ SubscriptBox["\[Phi]", "v"], "[", "3", "]"}]}]}], ",", RowBox[{ RowBox[{ RowBox[{"y", "[", "3", "]"}], RowBox[{"(", RowBox[{"2", "-", "X"}], ")"}]}], "\[Equal]", " ", "X"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"y", "[", "2", "]"}], RowBox[{"(", RowBox[{"2", "-", "X"}], ")"}]}], "\[Equal]", " ", RowBox[{"0.5", RowBox[{"(", RowBox[{"1", "-", "X"}], ")"}]}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"y", "[", "1", "]"}], RowBox[{"(", RowBox[{"2", "-", "X"}], ")"}]}], "\[Equal]", " ", RowBox[{"1.5", " ", RowBox[{"(", RowBox[{"1", "-", "X"}], ")"}]}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["Z", "v"], "^", "3"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", SubscriptBox["B", "v"]}], ")"}], RowBox[{ SubscriptBox["Z", "v"], "^", "2"}]}], "+", RowBox[{ SubscriptBox["Z", "v"], " ", RowBox[{"(", RowBox[{ SubscriptBox["A", "v"], "-", RowBox[{"3", RowBox[{ SubscriptBox["B", "v"], "^", "2"}]}], "-", RowBox[{"2", SubscriptBox["B", "v"]}]}], ")"}]}], "-", RowBox[{ SubscriptBox["A", "v"], " ", SubscriptBox["B", "v"]}], "+", RowBox[{ SubscriptBox["B", "v"], "^", "2"}], "+", RowBox[{ SubscriptBox["B", "v"], "^", "3"}]}], "\[Equal]", "0"}]}], "}"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{ RowBox[{"y", "[", "2", "]"}], ",", "0.2"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"y", "[", "1", "]"}], ",", "0.1"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"y", "[", "3", "]"}], ",", "0.4"}], "}"}], ",", RowBox[{"{", RowBox[{"X", ",", "0.6"}], "}"}], ",", RowBox[{"{", RowBox[{ SubscriptBox["Z", "v"], ",", "0.95"}], "}"}], ",", RowBox[{"MaxIterations", "\[Rule]", " ", "10000"}]}], "]"}]}], ",", RowBox[{ RowBox[{"Pp", "[", "i", "]"}], "=", "P"}], ",", RowBox[{ RowBox[{"Kv", "[", "i", "]"}], "=", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["\[Phi]", "v"], "[", "3", "]"}], "/", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["\[Phi]", "v"], "[", "1", "]"}], "^", "1.5"}], " ", RowBox[{ RowBox[{ SubscriptBox["\[Phi]", "v"], "[", "2", "]"}], "^", "0.5"}]}], ")"}]}], "/.", "sol1080"}]}], ",", RowBox[{ RowBox[{"X", "[", "i", "]"}], "=", RowBox[{"X", "/.", "sol1080"}]}], ",", RowBox[{"Print", "[", RowBox[{"i", ",", "\"\< \>\"", ",", "P"}], "]"}], ",", RowBox[{"P", "=", RowBox[{"P", "+", "1000"}]}], ",", RowBox[{"i", "=", RowBox[{"i", "+", "1"}]}]}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.4020638294180512`*^9, 3.402063847443971*^9}}], Cell[CellGroupData[{ Cell[BoxData[ InterpretationBox[ RowBox[{ "1", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "0.001`"}], SequenceForm[1, " ", 0.001], Editable->False]], "Print", CellChangeTimes->{3.4020638508488674`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "2", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1000.001`"}], SequenceForm[2, " ", 1000.001], Editable->False]], "Print", CellChangeTimes->{3.402063851049155*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "3", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "2000.001`"}], SequenceForm[3, " ", 2000.001], Editable->False]], "Print", CellChangeTimes->{3.4020638512794867`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "4", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "3000.001`"}], SequenceForm[4, " ", 3000.001], Editable->False]], "Print", CellChangeTimes->{3.4020638514998035`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "5", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "4000.001`"}], SequenceForm[5, " ", 4000.001], Editable->False]], "Print", CellChangeTimes->{3.40206385172012*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "6", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "5000.001`"}], SequenceForm[6, " ", 5000.001], Editable->False]], "Print", CellChangeTimes->{3.402063851950451*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "7", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "6000.001`"}], SequenceForm[7, " ", 6000.001], Editable->False]], "Print", CellChangeTimes->{3.402063852170768*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "8", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "7000.001`"}], SequenceForm[8, " ", 7000.001], Editable->False]], "Print", CellChangeTimes->{3.402063852401099*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "9", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "8000.001`"}], SequenceForm[9, " ", 8000.001], Editable->False]], "Print", CellChangeTimes->{3.402063852621416*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "10", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "9000.001`"}], SequenceForm[10, " ", 9000.001], Editable->False]], "Print", CellChangeTimes->{3.402063852821704*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "11", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "10000.001`"}], SequenceForm[11, " ", 10000.001], Editable->False]], "Print", CellChangeTimes->{3.4020638530219917`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "12", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "11000.001`"}], SequenceForm[12, " ", 11000.001], Editable->False]], "Print", CellChangeTimes->{3.4020638532122655`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "13", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "12000.001`"}], SequenceForm[13, " ", 12000.001], Editable->False]], "Print", CellChangeTimes->{3.402063853412554*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "14", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "13000.001`"}], SequenceForm[14, " ", 13000.001], Editable->False]], "Print", CellChangeTimes->{3.402063853602827*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "15", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "14000.001`"}], SequenceForm[15, " ", 14000.001], Editable->False]], "Print", CellChangeTimes->{3.402063853803115*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "16", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "15000.001`"}], SequenceForm[16, " ", 15000.001], Editable->False]], "Print", CellChangeTimes->{3.402063853993389*^9}] }, Open ]] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"tbl1080Kv", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Pp", "[", "i", "]"}], ",", RowBox[{"Chop", "[", RowBox[{"Kv", "[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "16"}], "}"}]}], "]"}]}], ";"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"pltkv1080", "=", RowBox[{"ListPlot", "[", RowBox[{"tbl1080Kv", ",", RowBox[{"Joined", "\[Rule]", " ", "True"}], ",", " ", RowBox[{"PlotStyle", "\[Rule]", " ", RowBox[{"Hue", "[", "0.3", "]"}]}], ",", RowBox[{"AspectRatio", "\[Rule]", "1"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.4020638594812803`*^9, 3.4020638707875376`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"tbl1080X", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Pp", "[", "i", "]"}], ",", RowBox[{"Chop", "[", RowBox[{"X", "[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "16"}], "}"}]}], "]"}]}], ";"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"pltX1080", "=", RowBox[{"ListPlot", "[", RowBox[{"tbl1080X", ",", RowBox[{"Joined", "\[Rule]", " ", "True"}], ",", " ", RowBox[{"PlotStyle", "\[Rule]", " ", RowBox[{"Hue", "[", "0.3", "]"}]}], ",", RowBox[{"AspectRatio", "\[Rule]", "1"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.402063878558712*^9, 3.402063886990837*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[TextData[{ "The correction factor ", Cell[BoxData[ FormBox[ SubscriptBox["K", "v"], TraditionalForm]]], " and the extent of reaction X versus P for T=700K or 1260R" }], "Subsubtitle", Background->RGBColor[1, 1, 0.658824]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"T", "=", RowBox[{"700", " ", "1.8"}]}]], "Input"], Cell[BoxData["1260.`"], "Output", CellChangeTimes->{3.4020638928893185`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"P", "=", "0.001"}], ";", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "<", "17"}], ",", RowBox[{"{", RowBox[{ RowBox[{"sol1260", "=", RowBox[{"FindRoot", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{"Ka", "[", "700", "]"}], RowBox[{ RowBox[{"(", RowBox[{"0.5", " ", RowBox[{"(", RowBox[{"1", "-", "X"}], ")"}]}], ")"}], "^", "0.5"}], " ", RowBox[{ RowBox[{"(", RowBox[{"1.5", " ", RowBox[{"(", RowBox[{"1", "-", "X"}], ")"}]}], ")"}], "^", "1.5"}], " ", RowBox[{ RowBox[{ SubscriptBox["\[Phi]", "v"], "[", "1", "]"}], "^", "1.5"}], " ", RowBox[{ RowBox[{ SubscriptBox["\[Phi]", "v"], "[", "2", "]"}], "^", "0.5"}]}], "\[Equal]", " ", RowBox[{"X", " ", RowBox[{ RowBox[{"(", RowBox[{"2", "-", "X"}], ")"}], "/", RowBox[{"(", RowBox[{"P", "/", "14.69"}], ")"}]}], " ", RowBox[{ SubscriptBox["\[Phi]", "v"], "[", "3", "]"}]}]}], ",", RowBox[{ RowBox[{ RowBox[{"y", "[", "3", "]"}], RowBox[{"(", RowBox[{"2", "-", "X"}], ")"}]}], "\[Equal]", " ", "X"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"y", "[", "2", "]"}], RowBox[{"(", RowBox[{"2", "-", "X"}], ")"}]}], "\[Equal]", " ", RowBox[{"0.5", RowBox[{"(", RowBox[{"1", "-", "X"}], ")"}]}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"y", "[", "1", "]"}], RowBox[{"(", RowBox[{"2", "-", "X"}], ")"}]}], "\[Equal]", " ", RowBox[{"1.5", " ", RowBox[{"(", RowBox[{"1", "-", "X"}], ")"}]}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["Z", "v"], "^", "3"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", SubscriptBox["B", "v"]}], ")"}], RowBox[{ SubscriptBox["Z", "v"], "^", "2"}]}], "+", RowBox[{ SubscriptBox["Z", "v"], " ", RowBox[{"(", RowBox[{ SubscriptBox["A", "v"], "-", RowBox[{"3", RowBox[{ SubscriptBox["B", "v"], "^", "2"}]}], "-", RowBox[{"2", SubscriptBox["B", "v"]}]}], ")"}]}], "-", RowBox[{ SubscriptBox["A", "v"], " ", SubscriptBox["B", "v"]}], "+", RowBox[{ SubscriptBox["B", "v"], "^", "2"}], "+", RowBox[{ SubscriptBox["B", "v"], "^", "3"}]}], "\[Equal]", "0"}]}], "}"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{ RowBox[{"y", "[", "2", "]"}], ",", "0.2"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"y", "[", "1", "]"}], ",", "0.1"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"y", "[", "3", "]"}], ",", "0.4"}], "}"}], ",", RowBox[{"{", RowBox[{"X", ",", "0.6"}], "}"}], ",", RowBox[{"{", RowBox[{ SubscriptBox["Z", "v"], ",", "0.95"}], "}"}], ",", RowBox[{"MaxIterations", "\[Rule]", " ", "10000"}]}], "]"}]}], ",", RowBox[{ RowBox[{"Pp", "[", "i", "]"}], "=", "P"}], ",", RowBox[{ RowBox[{"Kv", "[", "i", "]"}], "=", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["\[Phi]", "v"], "[", "3", "]"}], "/", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["\[Phi]", "v"], "[", "1", "]"}], "^", "1.5"}], " ", RowBox[{ RowBox[{ SubscriptBox["\[Phi]", "v"], "[", "2", "]"}], "^", "0.5"}]}], ")"}]}], "/.", "sol1260"}]}], ",", RowBox[{ RowBox[{"X", "[", "i", "]"}], "=", RowBox[{"X", "/.", "sol1260"}]}], ",", RowBox[{"Print", "[", RowBox[{"i", ",", "\"\< \>\"", ",", "P"}], "]"}], ",", RowBox[{"P", "=", RowBox[{"P", "+", "1000"}]}], ",", RowBox[{"i", "=", RowBox[{"i", "+", "1"}]}]}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.4020638961840563`*^9, 3.4020639048765554`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ InterpretationBox[ RowBox[{ "1", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "0.001`"}], SequenceForm[1, " ", 0.001], Editable->False]], "Print", CellChangeTimes->{3.4020639071097665`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "2", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1000.001`"}], SequenceForm[2, " ", 1000.001], Editable->False]], "Print", CellChangeTimes->{3.4020639073100543`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "3", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "2000.001`"}], SequenceForm[3, " ", 2000.001], Editable->False]], "Print", CellChangeTimes->{3.4020639075403852`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "4", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "3000.001`"}], SequenceForm[4, " ", 3000.001], Editable->False]], "Print", CellChangeTimes->{3.4020639077707167`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "5", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "4000.001`"}], SequenceForm[5, " ", 4000.001], Editable->False]], "Print", CellChangeTimes->{3.4020639079910336`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "6", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "5000.001`"}], SequenceForm[6, " ", 5000.001], Editable->False]], "Print", CellChangeTimes->{3.4020639082113504`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "7", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "6000.001`"}], SequenceForm[7, " ", 6000.001], Editable->False]], "Print", CellChangeTimes->{3.4020639084316673`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "8", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "7000.001`"}], SequenceForm[8, " ", 7000.001], Editable->False]], "Print", CellChangeTimes->{3.4020639086619987`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "9", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "8000.001`"}], SequenceForm[9, " ", 8000.001], Editable->False]], "Print", CellChangeTimes->{3.402063908882315*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "10", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "9000.001`"}], SequenceForm[10, " ", 9000.001], Editable->False]], "Print", CellChangeTimes->{3.4020639090926175`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "11", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "10000.001`"}], SequenceForm[11, " ", 10000.001], Editable->False]], "Print", CellChangeTimes->{3.4020639092828913`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "12", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "11000.001`"}], SequenceForm[12, " ", 11000.001], Editable->False]], "Print", CellChangeTimes->{3.402063909483179*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "13", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "12000.001`"}], SequenceForm[13, " ", 12000.001], Editable->False]], "Print", CellChangeTimes->{3.402063909673453*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "14", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "13000.001`"}], SequenceForm[14, " ", 13000.001], Editable->False]], "Print", CellChangeTimes->{3.402063909903784*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "15", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "14000.001`"}], SequenceForm[15, " ", 14000.001], Editable->False]], "Print", CellChangeTimes->{3.4020639101241007`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "16", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "15000.001`"}], SequenceForm[16, " ", 15000.001], Editable->False]], "Print", CellChangeTimes->{3.4020639103143744`*^9}] }, Open ]] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"tbl1260Kv", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Pp", "[", "i", "]"}], ",", RowBox[{"Chop", "[", RowBox[{"Kv", "[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "16"}], "}"}]}], "]"}]}], ";"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"pltkv1260", "=", RowBox[{"ListPlot", "[", RowBox[{"tbl1260Kv", ",", RowBox[{"Joined", "\[Rule]", " ", "True"}], ",", " ", RowBox[{"PlotStyle", "\[Rule]", " ", RowBox[{"Hue", "[", "0.5", "]"}]}], ",", RowBox[{"AspectRatio", "\[Rule]", "1"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.402063915822294*^9, 3.4020639269182496`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"tbl1260X", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Pp", "[", "i", "]"}], ",", RowBox[{"Chop", "[", RowBox[{"X", "[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "16"}], "}"}]}], "]"}]}], ";"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"pltX1260", "=", RowBox[{"ListPlot", "[", RowBox[{"tbl1260X", ",", RowBox[{"Joined", "\[Rule]", " ", "True"}], ",", " ", RowBox[{"PlotStyle", "\[Rule]", " ", RowBox[{"Hue", "[", "0.5", "]"}]}], ",", RowBox[{"AspectRatio", "\[Rule]", "1"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.4020639327967024`*^9, 3.402063941529259*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[TextData[{ "The correction factor ", Cell[BoxData[ FormBox[ SubscriptBox["K", "v"], TraditionalForm]]], " and the extent of reaction X versus P for T=800K or 1440R" }], "Subsubtitle", Background->RGBColor[1, 1, 0.658824]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"T", "=", RowBox[{"800", " ", "1.8"}]}]], "Input"], Cell[BoxData["1440.`"], "Output", CellChangeTimes->{3.4020639448239965`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"P", "=", "0.001"}], ";", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "<", "17"}], ",", RowBox[{"{", RowBox[{ RowBox[{"sol1440", "=", RowBox[{"FindRoot", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{"Ka", "[", "800", "]"}], RowBox[{ RowBox[{"(", RowBox[{"0.5", " ", RowBox[{"(", RowBox[{"1", "-", "X"}], ")"}]}], ")"}], "^", "0.5"}], " ", RowBox[{ RowBox[{"(", RowBox[{"1.5", " ", RowBox[{"(", RowBox[{"1", "-", "X"}], ")"}]}], ")"}], "^", "1.5"}], " ", RowBox[{ RowBox[{ SubscriptBox["\[Phi]", "v"], "[", "1", "]"}], "^", "1.5"}], " ", RowBox[{ RowBox[{ SubscriptBox["\[Phi]", "v"], "[", "2", "]"}], "^", "0.5"}]}], "\[Equal]", " ", RowBox[{"X", " ", RowBox[{ RowBox[{"(", RowBox[{"2", "-", "X"}], ")"}], "/", RowBox[{"(", RowBox[{"P", "/", "14.69"}], ")"}]}], " ", RowBox[{ SubscriptBox["\[Phi]", "v"], "[", "3", "]"}]}]}], ",", RowBox[{ RowBox[{ RowBox[{"y", "[", "3", "]"}], RowBox[{"(", RowBox[{"2", "-", "X"}], ")"}]}], "\[Equal]", " ", "X"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"y", "[", "2", "]"}], RowBox[{"(", RowBox[{"2", "-", "X"}], ")"}]}], "\[Equal]", " ", RowBox[{"0.5", RowBox[{"(", RowBox[{"1", "-", "X"}], ")"}]}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"y", "[", "1", "]"}], RowBox[{"(", RowBox[{"2", "-", "X"}], ")"}]}], "\[Equal]", " ", RowBox[{"1.5", " ", RowBox[{"(", RowBox[{"1", "-", "X"}], ")"}]}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["Z", "v"], "^", "3"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", SubscriptBox["B", "v"]}], ")"}], RowBox[{ SubscriptBox["Z", "v"], "^", "2"}]}], "+", RowBox[{ SubscriptBox["Z", "v"], " ", RowBox[{"(", RowBox[{ SubscriptBox["A", "v"], "-", RowBox[{"3", RowBox[{ SubscriptBox["B", "v"], "^", "2"}]}], "-", RowBox[{"2", SubscriptBox["B", "v"]}]}], ")"}]}], "-", RowBox[{ SubscriptBox["A", "v"], " ", SubscriptBox["B", "v"]}], "+", RowBox[{ SubscriptBox["B", "v"], "^", "2"}], "+", RowBox[{ SubscriptBox["B", "v"], "^", "3"}]}], "\[Equal]", "0"}]}], "}"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{ RowBox[{"y", "[", "2", "]"}], ",", "0.2"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"y", "[", "1", "]"}], ",", "0.1"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"y", "[", "3", "]"}], ",", "0.4"}], "}"}], ",", RowBox[{"{", RowBox[{"X", ",", "0.6"}], "}"}], ",", RowBox[{"{", RowBox[{ SubscriptBox["Z", "v"], ",", "0.95"}], "}"}], ",", RowBox[{"MaxIterations", "\[Rule]", " ", "10000"}]}], "]"}]}], ",", RowBox[{ RowBox[{"Pp", "[", "i", "]"}], "=", "P"}], ",", RowBox[{ RowBox[{"Kv", "[", "i", "]"}], "=", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["\[Phi]", "v"], "[", "3", "]"}], "/", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["\[Phi]", "v"], "[", "1", "]"}], "^", "1.5"}], " ", RowBox[{ RowBox[{ SubscriptBox["\[Phi]", "v"], "[", "2", "]"}], "^", "0.5"}]}], ")"}]}], "/.", "sol1440"}]}], ",", RowBox[{ RowBox[{"X", "[", "i", "]"}], "=", RowBox[{"X", "/.", "sol1440"}]}], ",", RowBox[{"Print", "[", RowBox[{"i", ",", "\"\< \>\"", ",", "P"}], "]"}], ",", RowBox[{"P", "=", RowBox[{"P", "+", "1000"}]}], ",", RowBox[{"i", "=", RowBox[{"i", "+", "1"}]}]}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.402063950512176*^9, 3.402063963641054*^9}}], Cell[CellGroupData[{ Cell[BoxData[ InterpretationBox[ RowBox[{ "1", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "0.001`"}], SequenceForm[1, " ", 0.001], Editable->False]], "Print", CellChangeTimes->{3.402063964952941*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "2", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1000.001`"}], SequenceForm[2, " ", 1000.001], Editable->False]], "Print", CellChangeTimes->{3.4020639651231856`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "3", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "2000.001`"}], SequenceForm[3, " ", 2000.001], Editable->False]], "Print", CellChangeTimes->{3.402063965323474*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "4", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "3000.001`"}], SequenceForm[4, " ", 3000.001], Editable->False]], "Print", CellChangeTimes->{3.4020639655437903`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "5", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "4000.001`"}], SequenceForm[5, " ", 4000.001], Editable->False]], "Print", CellChangeTimes->{3.4020639657240496`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "6", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "5000.001`"}], SequenceForm[6, " ", 5000.001], Editable->False]], "Print", CellChangeTimes->{3.402063965914323*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "7", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "6000.001`"}], SequenceForm[7, " ", 6000.001], Editable->False]], "Print", CellChangeTimes->{3.4020639661346397`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "8", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "7000.001`"}], SequenceForm[8, " ", 7000.001], Editable->False]], "Print", CellChangeTimes->{3.4020639663148994`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "9", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "8000.001`"}], SequenceForm[9, " ", 8000.001], Editable->False]], "Print", CellChangeTimes->{3.4020639665151873`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "10", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "9000.001`"}], SequenceForm[10, " ", 9000.001], Editable->False]], "Print", CellChangeTimes->{3.4020639667254896`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "11", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "10000.001`"}], SequenceForm[11, " ", 10000.001], Editable->False]], "Print", CellChangeTimes->{3.402063966905749*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "12", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "11000.001`"}], SequenceForm[12, " ", 11000.001], Editable->False]], "Print", CellChangeTimes->{3.4020639670960226`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "13", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "12000.001`"}], SequenceForm[13, " ", 12000.001], Editable->False]], "Print", CellChangeTimes->{3.402063967306325*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "14", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "13000.001`"}], SequenceForm[14, " ", 13000.001], Editable->False]], "Print", CellChangeTimes->{3.4020639674865837`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "15", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "14000.001`"}], SequenceForm[15, " ", 14000.001], Editable->False]], "Print", CellChangeTimes->{3.402063967676858*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "16", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "15000.001`"}], SequenceForm[16, " ", 15000.001], Editable->False]], "Print", CellChangeTimes->{3.402063967877146*^9}] }, Open ]] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"tbl1440Kv", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Pp", "[", "i", "]"}], ",", RowBox[{"Chop", "[", RowBox[{"Kv", "[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "16"}], "}"}]}], "]"}]}], ";"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"pltkv1440", "=", RowBox[{"ListPlot", "[", RowBox[{"tbl1440Kv", ",", RowBox[{"Joined", "\[Rule]", " ", "True"}], ",", " ", RowBox[{"PlotStyle", "\[Rule]", " ", RowBox[{"Hue", "[", "0.7", "]"}]}], ",", RowBox[{"AspectRatio", "\[Rule]", "1"}], ",", RowBox[{"DisplayFunction", "\[Rule]", "Identity"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.4020639740360017`*^9, 3.4020639790932736`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"tbl1440X", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Pp", "[", "i", "]"}], ",", RowBox[{"Chop", "[", RowBox[{"X", "[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "16"}], "}"}]}], "]"}]}], ";"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"pltX1440", "=", RowBox[{"ListPlot", "[", RowBox[{"tbl1440X", ",", RowBox[{"Joined", "\[Rule]", " ", "True"}], ",", " ", RowBox[{"PlotStyle", "\[Rule]", " ", RowBox[{"Hue", "[", "0.7", "]"}]}], ",", RowBox[{"AspectRatio", "\[Rule]", "1"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.402063985512504*^9, 3.402063997159251*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[TextData[{ "Plots of the correction factor ", Cell[BoxData[ FormBox[ SubscriptBox["K", "v"], TraditionalForm]]], " and the extent of reaction X versus P for T from 500K to 800K" }], "Subsubtitle", Background->RGBColor[1, 1, 0.658824]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Show", "[", RowBox[{"pltX900", ",", "pltX1080", ",", "pltX1260", ",", "pltX1440", ",", RowBox[{"PlotRange", "\[Rule]", " ", "All"}]}], "]"}]], "Input", CellChangeTimes->{{3.4020640062823696`*^9, 3.4020640204527454`*^9}}], Cell[BoxData[ GraphicsBox[{{{}, {}, {Hue[0.1], LineBox[{{746.5452996587518, 0.6288652115806752}, {1000.001, 0.8423633175968963}, {2000.001, 0.9096860508318435}, {3000.001, 0.9434600901453007}, {4000.001, 0.9635412617409325}, {5000.001, 0.9751272611180043}, {6000.001, 0.9816229662026373}, {7000.001, 0.9854705899928589}, {8000.001, 0.9879530738977143}, {9000.001, 1.0089724234633777`}, {10000.001, 1.0080215378620545`}, {11000.001, 1.0072610965529845`}, {12000.001, 1.0066371105515266`}, {13000.001, 1.0061144793017534`}, {14000.001, 1.005669364758984}, {15000.001, 1.0052849850555339`}}]}}, {{}, {}, {Hue[0.3], LineBox[{{0.001, 1.922024696597683*^-6}, {1000.001, 0.5628402760587687}, {2000.001, 0.69577388543499}, {3000.001, 0.7670998654087623}, {4000.001, 0.814157570981663}, {5000.001, 0.8480723664988894}, {6000.001, 0.8735739664862774}, {7000.001, 0.8931656077619449}, {8000.001, 0.9084095533486498}, {9000.001, 0.9203951264146435}, {10000.001, 0.9299235458562076}, {11000.001, 0.9375924355297982}, {12000.001, 0.9438465315534347}, {13000.001, 0.9490153860306911}, {14000.001, 0.9533430407391547}, {15000.001, 0.9570109057319238}}]}}, {{}, {}, {Hue[0.5], LineBox[{{0.001, 4.3778230684483764`*^-7}, {1000.001, 0.2766849960769658}, {2000.001, 0.4156361118042908}, {3000.001, 0.5055045734394165}, {4000.001, 0.5707569875761376}, {5000.001, 0.6213132467086458}, {6000.001, 0.6620941740547528}, {7000.001, 0.6958766975251396}, {8000.001, 0.7243795580232582}, {9000.001, 0.7487452956915669}, {10000.001, 0.7697785209090604}, {11000.001, 0.7880732200777746}, {12000.001, 0.8040848533152096}, {13000.001, 0.8181731449248372}, {14000.001, 0.830628577294057}, {15000.001, 0.8416895118019874}}]}}, {{}, {}, {Hue[0.7], LineBox[{{0.001, 1.4150979313800966`*^-7}, {1000.001, 0.11960398768285088`}, {2000.001, 0.20880560619471442`}, {3000.001, 0.27935291623012987`}, {4000.001, 0.33736210164098857`}, {5000.001, 0.3863862966505641}, {6000.001, 0.4286579211958242}, {7000.001, 0.4656651897110545}, {8000.001, 0.4984476030279549}, {9000.001, 0.5277593488826187}, {10000.001, 0.554165237632287}, {11000.001, 0.5780998230988788}, {12000.001, 0.5999053182773986}, {13000.001, 0.6198567360113179}, {14000.001, 0.6381790371171789}, {15000.001, 0.6550591153725122}}]}}}, AspectRatio->1, Axes->True, AxesOrigin->{0, 0.65}, PlotRange->All, PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Scaled[0.02]}]], "Output", CellChangeTimes->{{3.4020640074440403`*^9, 3.402064021243883*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Show", "[", RowBox[{ "pltkv900", ",", "pltkv1080", ",", "pltkv1260", ",", "pltkv1440", ",", RowBox[{"PlotRange", "\[Rule]", " ", "All"}]}], "]"}]], "Input", CellChangeTimes->{{3.402064044597464*^9, 3.4020640467405453`*^9}}], Cell[BoxData[ GraphicsBox[{{{}, {}, {Hue[0.1], LineBox[{{0.001, 0.9999999561912054}, {1000.001, 0.739710340255964}, { 2000.001, 0.47743602399416657`}, {3000.001, 0.27928048622011487`}, { 4000.001, 0.15454659612475247`}, {5000.001, 0.08984718658788614}, { 6000.001, 0.05883926687135821}, {7000.001, 0.042904629985059936`}, { 8000.001, 0.033707300769598235`}, {9000.001, 0.021033689366250777`}, { 10000.001, 0.01867934630628976}, {11000.001, 0.016835982108261}, { 12000.001, 0.015345354519875277`}, {13000.001, 0.014109024597384277`}, { 14000.001, 0.013062588166411168`}, {15000.001, 0.012162121982344024`}}]}}, {{}, {}, {Hue[0.3], LineBox[{{0.001, 0.9999999658016204}, {1000.001, 0.9081983570665705}, { 2000.001, 0.7841377731837917}, {3000.001, 0.66141077941578}, {4000.001, 0.5500530687143101}, {5000.001, 0.45412532842469233`}, {6000.001, 0.37463830716941837`}, {7000.001, 0.31066750460841497`}, {8000.001, 0.2601591905720575}, {9000.001, 0.2206342453569687}, {10000.001, 0.1897022250258442}, {11000.001, 0.1653301806551417}, {12000.001, 0.14591375646519192`}, {13000.001, 0.13023935640150408`}, {14000.001, 0.11740796487007615`}, {15000.001, 0.10675824373249315`}}]}}, {{}, {}, {Hue[0.5], LineBox[{{0.001, 0.9999999726601807}, {1000.001, 0.9607128024628265}, { 2000.001, 0.9080655897064472}, {3000.001, 0.8501875735282057}, { 4000.001, 0.7910386939774348}, {5000.001, 0.7328963800769028}, { 6000.001, 0.6771515004669837}, {7000.001, 0.6246475246496477}, { 8000.001, 0.5758557517142711}, {9000.001, 0.5309817638491868}, { 10000.001, 0.4900395170671679}, {11000.001, 0.45290793102121807`}, { 12000.001, 0.41937596441202085`}, {13000.001, 0.3891784996865422}, { 14000.001, 0.3620240872781138}, {15000.001, 0.33761538920020306`}}]}}, {{}, {}, {Hue[0.7], LineBox[{{0.001, 0.9999999777120031}, {1000.001, 0.9753896906705986}, { 2000.001, 0.9473889470148663}, {3000.001, 0.9173566818388598}, { 4000.001, 0.8862071444051042}, {5000.001, 0.8545854141204814}, { 6000.001, 0.8229599950148034}, {7000.001, 0.7916764028418622}, { 8000.001, 0.760990495591278}, {9000.001, 0.7310904813926967}, { 10000.001, 0.7021122204023061}, {11000.001, 0.6741503574770941}, { 12000.001, 0.6472667464952884}, {13000.001, 0.6214970384423774}, { 14000.001, 0.5968559690936162}, {15000.001, 0.5733416840157393}}]}}}, AspectRatio->1, Axes->True, AxesOrigin->{0, 0}, PlotRange->All, PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Scaled[0.02]}]], "Output", CellChangeTimes->{3.4020640479122305`*^9}] }, Open ]] }, Open ]] }, Open ]] }, WindowSize->{1016, 651}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, FrontEndVersion->"6.0 for Microsoft Windows (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, 780, 24, 244, "Title"], Cell[1373, 49, 348, 8, 85, "Subsubtitle"], Cell[1724, 59, 99, 2, 31, "Input"], Cell[CellGroupData[{ Cell[1848, 65, 455, 16, 88, "Subsubtitle"], Cell[CellGroupData[{ Cell[2328, 85, 200, 6, 31, "Input"], Cell[2531, 93, 89, 1, 30, "Output"] }, Open ]], Cell[2635, 97, 787, 25, 31, "Input"], Cell[CellGroupData[{ Cell[3447, 126, 164, 4, 31, "Input"], Cell[3614, 132, 94, 2, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[3745, 139, 312, 10, 31, "Input"], Cell[4060, 151, 79, 1, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[4176, 157, 334, 11, 31, "Input"], Cell[4513, 170, 108, 2, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[4658, 177, 298, 10, 31, "Input"], Cell[4959, 189, 80, 1, 30, "Output"] }, Open ]], Cell[5054, 193, 1789, 53, 92, "Input"], Cell[CellGroupData[{ Cell[6868, 250, 57, 1, 31, "Input"], Cell[6928, 253, 87, 1, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[7052, 259, 57, 1, 31, "Input"], Cell[7112, 262, 89, 1, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[7238, 268, 57, 1, 31, "Input"], Cell[7298, 271, 89, 1, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[7424, 277, 57, 1, 31, "Input"], Cell[7484, 280, 93, 1, 30, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[7626, 287, 94, 1, 44, "Subsubtitle"], Cell[7723, 290, 134, 4, 31, "Input"], Cell[7860, 296, 134, 4, 31, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[8031, 305, 306, 5, 69, "Subsubtitle"], Cell[8340, 312, 300, 9, 31, "Input"], Cell[8643, 323, 260, 9, 31, "Input"], Cell[8906, 334, 205, 6, 31, "Input"], Cell[9114, 342, 209, 7, 31, "Input"], Cell[9326, 351, 1647, 56, 78, "Input"], Cell[10976, 409, 336, 12, 31, "Input"], Cell[11315, 423, 278, 8, 31, "Input"], Cell[11596, 433, 424, 14, 56, "Input"], Cell[12023, 449, 402, 14, 31, "Input"], Cell[12428, 465, 263, 9, 54, "Input"], Cell[12694, 476, 1081, 38, 92, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[13812, 519, 74, 1, 44, "Subsubtitle"], Cell[13889, 522, 71, 2, 31, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[13997, 529, 160, 4, 44, "Subsubtitle"], Cell[14160, 535, 238, 7, 31, "Input"], Cell[14401, 544, 242, 7, 31, "Input"], Cell[14646, 553, 243, 7, 31, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[14926, 565, 234, 7, 44, "Subsubtitle"], Cell[15163, 574, 69, 2, 31, "Input"], Cell[CellGroupData[{ Cell[15257, 580, 4944, 138, 172, "Input"], Cell[CellGroupData[{ Cell[20226, 722, 223, 6, 23, "Print"], Cell[20452, 730, 229, 6, 23, "Print"], Cell[20684, 738, 229, 6, 23, "Print"], Cell[20916, 746, 227, 6, 23, "Print"], Cell[21146, 754, 227, 6, 23, "Print"], Cell[21376, 762, 227, 6, 23, "Print"], Cell[21606, 770, 229, 6, 23, "Print"], Cell[21838, 778, 229, 6, 23, "Print"], Cell[22070, 786, 229, 6, 23, "Print"], Cell[22302, 794, 231, 6, 23, "Print"], Cell[22536, 802, 236, 7, 23, "Print"], Cell[22775, 811, 236, 7, 23, "Print"], Cell[23014, 820, 236, 7, 23, "Print"], Cell[23253, 829, 236, 7, 23, "Print"], Cell[23492, 838, 234, 7, 23, "Print"], Cell[23729, 847, 234, 7, 23, "Print"] }, Open ]] }, Open ]], Cell[23990, 858, 354, 11, 31, "Input"], Cell[24347, 871, 401, 9, 31, "Input"], Cell[24751, 882, 352, 11, 31, "Input"], Cell[25106, 895, 401, 9, 31, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[25544, 909, 235, 7, 44, "Subsubtitle"], Cell[CellGroupData[{ Cell[25804, 920, 76, 2, 31, "Input"], Cell[25883, 924, 77, 1, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[25997, 930, 4948, 138, 172, "Input"], Cell[CellGroupData[{ Cell[30970, 1072, 223, 6, 23, "Print"], Cell[31196, 1080, 227, 6, 23, "Print"], Cell[31426, 1088, 229, 6, 23, "Print"], Cell[31658, 1096, 229, 6, 23, "Print"], Cell[31890, 1104, 226, 6, 23, "Print"], Cell[32119, 1112, 227, 6, 23, "Print"], Cell[32349, 1120, 227, 6, 23, "Print"], Cell[32579, 1128, 227, 6, 23, "Print"], Cell[32809, 1136, 227, 6, 23, "Print"], Cell[33039, 1144, 229, 6, 23, "Print"], Cell[33271, 1152, 236, 7, 23, "Print"], Cell[33510, 1161, 236, 7, 23, "Print"], Cell[33749, 1170, 234, 7, 23, "Print"], Cell[33986, 1179, 234, 7, 23, "Print"], Cell[34223, 1188, 234, 7, 23, "Print"], Cell[34460, 1197, 234, 7, 23, "Print"] }, Open ]] }, Open ]], Cell[34721, 1208, 355, 11, 31, "Input"], Cell[35079, 1221, 407, 9, 31, "Input"], Cell[35489, 1232, 353, 11, 31, "Input"], Cell[35845, 1245, 401, 9, 31, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[36283, 1259, 235, 7, 44, "Subsubtitle"], Cell[CellGroupData[{ Cell[36543, 1270, 76, 2, 31, "Input"], Cell[36622, 1274, 77, 1, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[36736, 1280, 4945, 138, 172, "Input"], Cell[CellGroupData[{ Cell[41706, 1422, 223, 6, 23, "Print"], Cell[41932, 1430, 229, 6, 23, "Print"], Cell[42164, 1438, 229, 6, 23, "Print"], Cell[42396, 1446, 229, 6, 23, "Print"], Cell[42628, 1454, 229, 6, 23, "Print"], Cell[42860, 1462, 229, 6, 23, "Print"], Cell[43092, 1470, 229, 6, 23, "Print"], Cell[43324, 1478, 229, 6, 23, "Print"], Cell[43556, 1486, 227, 6, 23, "Print"], Cell[43786, 1494, 231, 6, 23, "Print"], Cell[44020, 1502, 236, 7, 23, "Print"], Cell[44259, 1511, 234, 7, 23, "Print"], Cell[44496, 1520, 234, 7, 23, "Print"], Cell[44733, 1529, 234, 7, 23, "Print"], Cell[44970, 1538, 236, 7, 23, "Print"], Cell[45209, 1547, 236, 7, 23, "Print"] }, Open ]] }, Open ]], Cell[45472, 1558, 355, 11, 31, "Input"], Cell[45830, 1571, 405, 9, 31, "Input"], Cell[46238, 1582, 353, 11, 31, "Input"], Cell[46594, 1595, 403, 9, 31, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[47034, 1609, 235, 7, 44, "Subsubtitle"], Cell[CellGroupData[{ Cell[47294, 1620, 76, 2, 31, "Input"], Cell[47373, 1624, 77, 1, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[47487, 1630, 4941, 138, 172, "Input"], Cell[CellGroupData[{ Cell[52453, 1772, 221, 6, 23, "Print"], Cell[52677, 1780, 229, 6, 23, "Print"], Cell[52909, 1788, 227, 6, 23, "Print"], Cell[53139, 1796, 229, 6, 23, "Print"], Cell[53371, 1804, 229, 6, 23, "Print"], Cell[53603, 1812, 227, 6, 23, "Print"], Cell[53833, 1820, 229, 6, 23, "Print"], Cell[54065, 1828, 229, 6, 23, "Print"], Cell[54297, 1836, 229, 6, 23, "Print"], Cell[54529, 1844, 231, 6, 23, "Print"], Cell[54763, 1852, 234, 7, 23, "Print"], Cell[55000, 1861, 236, 7, 23, "Print"], Cell[55239, 1870, 234, 7, 23, "Print"], Cell[55476, 1879, 236, 7, 23, "Print"], Cell[55715, 1888, 234, 7, 23, "Print"], Cell[55952, 1897, 234, 7, 23, "Print"] }, Open ]] }, Open ]], Cell[56213, 1908, 355, 11, 31, "Input"], Cell[56571, 1921, 473, 11, 31, "Input"], Cell[57047, 1934, 353, 11, 31, "Input"], Cell[57403, 1947, 401, 9, 31, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[57841, 1961, 248, 7, 44, "Subsubtitle"], Cell[CellGroupData[{ Cell[58114, 1972, 252, 4, 31, "Input"], Cell[58369, 1978, 2763, 50, 366, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[61169, 2033, 257, 5, 31, "Input"], Cell[61429, 2040, 2734, 48, 374, "Output"] }, Open ]] }, Open ]] }, Open ]] } ] *) (* End of internal cache information *)