(* 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[ 199437, 4784] NotebookOptionsPosition[ 180755, 4280] NotebookOutlinePosition[ 181096, 4295] CellTagsIndexPosition[ 181053, 4292] WindowFrame->Normal ContainsDynamic->False*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[TextData[{ StyleBox["Rectifying and Stripping Profiles\nfor the Reactive Cascade \n", FontSize->18], StyleBox["Author's Data: ", FontSize->16], StyleBox["Housam BINOUS\nDepartment of Chemical Engineering\nNational \ Institute of Applied Sciences and Technology\nTunis, TUNISIA\nEmail: \ binoushousam@yahoo.com ", FontSize->14, FontWeight->"Plain"] }], "Title", CellChangeTimes->{{3.410094392802664*^9, 3.410094413041766*^9}, { 3.41017622093452*^9, 3.410176243266632*^9}}, TextAlignment->Center, Background->RGBColor[0.531258, 0.996109, 0.765637]], Cell[TextData[{ StyleBox["Reference :", FontWeight->"Bold"], "\n\n", StyleBox["Doherty and Malone, Conceptual Design of Distillation Systems, \ McGraw-Hill, 2001. \n", FontFamily->"Arial"] }], "Subsubtitle", CellChangeTimes->{ 3.400035572909675*^9, {3.408263426288051*^9, 3.4082634490307536`*^9}, { 3.409313583211299*^9, 3.4093137754076643`*^9}, {3.410094867304965*^9, 3.4100948789817553`*^9}, {3.4101718220091763`*^9, 3.4101719154735713`*^9}}, Background->RGBColor[0.773449, 0.996109, 0.996109]], Cell[CellGroupData[{ Cell["Relative volatilities", "Subsubtitle", CellGroupingRules->{GroupTogetherGrouping, 10000.}, CellChangeTimes->{{3.4101720380698557`*^9, 3.4101720550142207`*^9}, 3.4101840250663104`*^9, 3.4101849917463284`*^9}, Background->RGBColor[0.996109, 0.996109, 0.640635]], Cell[BoxData[ RowBox[{ RowBox[{ SubscriptBox["\[Alpha]", "AC"], "=", "5"}], ";", RowBox[{ SubscriptBox["\[Alpha]", "BC"], "=", "3"}], ";", RowBox[{ SubscriptBox["\[Alpha]", "CC"], "=", "1"}], ";"}]], "Input", CellGroupingRules->{GroupTogetherGrouping, 10000.}, CellChangeTimes->{{3.410171934871464*^9, 3.410171939097541*^9}, { 3.410171999985093*^9, 3.4101720550142207`*^9}, 3.4101840250663104`*^9, 3.4101849917463284`*^9}] }, Closed]], Cell[CellGroupData[{ Cell["Molecular Weigths", "Subsubtitle", CellChangeTimes->{{3.4101720380698557`*^9, 3.4101720550142207`*^9}, { 3.41018400046093*^9, 3.4101840250663104`*^9}}, Background->RGBColor[0.996109, 0.996109, 0.640635]], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"M", "[", "1", "]"}], "=", "86.18"}], ";", RowBox[{ RowBox[{"M", "[", "2", "]"}], "=", "100.21"}], ";", RowBox[{ RowBox[{"M", "[", "3", "]"}], "=", "186.39"}], ";"}]], "Input", CellChangeTimes->{{3.410184053837682*^9, 3.4101840947965775`*^9}, { 3.410184174901763*^9, 3.410184186258093*^9}}] }, Closed]], Cell[CellGroupData[{ Cell["Equilibrium constant for reaction", "Subsubtitle", CellChangeTimes->{ 3.410091630290366*^9, {3.410092930199544*^9, 3.4100929425773425`*^9}, { 3.4101729669455137`*^9, 3.4101729684076157`*^9}}, Background->RGBColor[0.996109, 0.996109, 0.640635]], Cell[BoxData[ RowBox[{ RowBox[{"Keq", "=", "2"}], ";"}]], "Input", CellChangeTimes->{{3.4100929450008273`*^9, 3.410092951099597*^9}, { 3.4101729723933473`*^9, 3.4101729831688414`*^9}}] }, Closed]], Cell[CellGroupData[{ Cell["\<\ Rectifying Cascade Trajectory at Damk\[ODoubleDot]hler number Da=0.277\ \>", "Subsubtitle", CellChangeTimes->{ 3.410091630290366*^9, {3.410092930199544*^9, 3.4100929425773425`*^9}, { 3.410093659297936*^9, 3.4100937216576047`*^9}, {3.410093789705453*^9, 3.4100937938814573`*^9}, {3.4100939430860033`*^9, 3.410093947842843*^9}, { 3.410172356417618*^9, 3.410172359772442*^9}, {3.4101759060216975`*^9, 3.4101759104080048`*^9}}, Background->RGBColor[0.996109, 0.996109, 0.640635]], Cell[BoxData[ RowBox[{ RowBox[{"t", "=", "0"}], ";", RowBox[{ RowBox[{"y", "[", RowBox[{"1", ",", RowBox[{"-", "1"}]}], "]"}], "=", "0.4"}], ";", RowBox[{ RowBox[{"y", "[", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "]"}], "=", "0.4"}], ";", RowBox[{ RowBox[{"y", "[", RowBox[{"3", ",", RowBox[{"-", "1"}]}], "]"}], "=", "0.2"}], ";", RowBox[{ RowBox[{"x", "[", RowBox[{"1", ",", RowBox[{"-", "1"}]}], "]"}], "=", "0.4"}], ";", RowBox[{ RowBox[{"x", "[", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "]"}], "=", "0.4"}], ";", RowBox[{ RowBox[{"x", "[", RowBox[{"3", ",", RowBox[{"-", "1"}]}], "]"}], "=", "0.2"}], ";"}]], "Input", CellChangeTimes->{{3.410091642257574*^9, 3.4100917039863358`*^9}, { 3.4100920383070655`*^9, 3.410092045727736*^9}, {3.410092192739128*^9, 3.4100921996590786`*^9}, {3.4100926628350925`*^9, 3.4100926712471886`*^9}, { 3.4101724027141886`*^9, 3.4101724283210096`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"\[Phi]", "=", "0.5"}], ";", RowBox[{"Da", "=", "0.277"}], ";", " "}]], "Input", CellChangeTimes->{{3.41009177143332*^9, 3.410091803018738*^9}, { 3.410172437514229*^9, 3.410172446336915*^9}, {3.410172602381296*^9, 3.410172613166805*^9}, 3.4101740491115923`*^9, {3.410175895015872*^9, 3.4101758951260304`*^9}, {3.4101760247223806`*^9, 3.4101760248425536`*^9}}],\ Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"While", "[", RowBox[{ RowBox[{"t", "<", "60"}], ",", "\[IndentingNewLine]", RowBox[{"{", "\[IndentingNewLine]", RowBox[{ RowBox[{"sol", "=", RowBox[{"Chop", "[", RowBox[{"FindRoot", "[", RowBox[{ RowBox[{"{", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"y", "[", RowBox[{"1", ",", RowBox[{"t", "-", "1"}]}], "]"}], "\[Equal]", " ", RowBox[{ RowBox[{"\[Phi]", " ", FractionBox[ RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "3"], RowBox[{ RowBox[{"y", "[", RowBox[{"i", ",", RowBox[{"t", "-", "1"}]}], "]"}], " ", RowBox[{"M", "[", "i", "]"}]}]}], RowBox[{ RowBox[{"y1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"y2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"y3", " ", RowBox[{"M", "[", "3", "]"}]}]}]], " ", "y1"}], "+", RowBox[{ RowBox[{"(", RowBox[{"1", "-", RowBox[{"\[Phi]", " ", FractionBox[ RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "3"], RowBox[{ RowBox[{"y", "[", RowBox[{"i", ",", RowBox[{"t", "-", "1"}]}], "]"}], " ", RowBox[{"M", "[", "i", "]"}]}]}], RowBox[{ RowBox[{"y1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"y2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"y3", " ", RowBox[{"M", "[", "3", "]"}]}]}]]}]}], ")"}], " ", "x1"}], "-", RowBox[{ FractionBox[ RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "3"], RowBox[{ RowBox[{"y", "[", RowBox[{"i", ",", RowBox[{"t", "-", "1"}]}], "]"}], " ", RowBox[{"M", "[", "i", "]"}]}]}], RowBox[{ RowBox[{"x1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"x2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"x3", " ", RowBox[{"M", "[", "3", "]"}]}]}]], "Da", " ", RowBox[{"(", RowBox[{ RowBox[{"x1", " ", "x2"}], "-", RowBox[{"x3", "/", "Keq"}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "x1"}], ")"}]}]}]}], ",", RowBox[{ RowBox[{"y", "[", RowBox[{"2", ",", RowBox[{"t", "-", "1"}]}], "]"}], "\[Equal]", " ", RowBox[{ RowBox[{"\[Phi]", " ", FractionBox[ RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "3"], RowBox[{ RowBox[{"y", "[", RowBox[{"i", ",", RowBox[{"t", "-", "1"}]}], "]"}], " ", RowBox[{"M", "[", "i", "]"}]}]}], RowBox[{ RowBox[{"y1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"y2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"y3", " ", RowBox[{"M", "[", "3", "]"}]}]}]], " ", "y2"}], "+", RowBox[{ RowBox[{"(", RowBox[{"1", "-", RowBox[{"\[Phi]", FractionBox[ RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "3"], RowBox[{ RowBox[{"y", "[", RowBox[{"i", ",", RowBox[{"t", "-", "1"}]}], "]"}], " ", RowBox[{"M", "[", "i", "]"}]}]}], RowBox[{ RowBox[{"y1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"y2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"y3", " ", RowBox[{"M", "[", "3", "]"}]}]}]]}]}], ")"}], " ", "x2"}], "-", RowBox[{ FractionBox[ RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "3"], RowBox[{ RowBox[{"y", "[", RowBox[{"i", ",", RowBox[{"t", "-", "1"}]}], "]"}], " ", RowBox[{"M", "[", "i", "]"}]}]}], RowBox[{ RowBox[{"x1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"x2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"x3", " ", RowBox[{"M", "[", "3", "]"}]}]}]], "Da", " ", RowBox[{"(", RowBox[{ RowBox[{"x1", " ", "x2"}], "-", RowBox[{"x3", "/", "Keq"}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "x2"}], ")"}]}]}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "[", RowBox[{"3", ",", RowBox[{"t", "-", "1"}]}], "]"}], "\[Equal]", " ", RowBox[{ RowBox[{"\[Phi]", " ", FractionBox[ RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "3"], RowBox[{ RowBox[{"y", "[", RowBox[{"i", ",", RowBox[{"t", "-", "1"}]}], "]"}], " ", RowBox[{"M", "[", "i", "]"}]}]}], RowBox[{ RowBox[{"y1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"y2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"y3", " ", RowBox[{"M", "[", "3", "]"}]}]}]], "y3"}], "+", RowBox[{ RowBox[{"(", RowBox[{"1", "-", RowBox[{"\[Phi]", FractionBox[ RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "3"], RowBox[{ RowBox[{"y", "[", RowBox[{"i", ",", RowBox[{"t", "-", "1"}]}], "]"}], " ", RowBox[{"M", "[", "i", "]"}]}]}], RowBox[{ RowBox[{"y1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"y2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"y3", " ", RowBox[{"M", "[", "3", "]"}]}]}]]}]}], ")"}], " ", "x3"}], "-", RowBox[{ FractionBox[ RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "3"], RowBox[{ RowBox[{"y", "[", RowBox[{"i", ",", RowBox[{"t", "-", "1"}]}], "]"}], " ", RowBox[{"M", "[", "i", "]"}]}]}], RowBox[{ RowBox[{"x1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"x2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"x3", " ", RowBox[{"M", "[", "3", "]"}]}]}]], "Da", " ", RowBox[{"(", RowBox[{ RowBox[{"x1", " ", "x2"}], "-", RowBox[{"x3", "/", "Keq"}]}], ")"}], " ", RowBox[{"(", RowBox[{"1", "+", "x3"}], ")"}]}]}]}], ",", "\[IndentingNewLine]", RowBox[{"y1", "==", RowBox[{"x1", "/", RowBox[{"(", RowBox[{"x1", "+", RowBox[{ RowBox[{ SubscriptBox["\[Alpha]", "BC"], "/", SubscriptBox["\[Alpha]", "AC"]}], " ", "x2"}], "+", RowBox[{ RowBox[{"1", "/", SubscriptBox["\[Alpha]", "AC"]}], " ", "x3"}]}], ")"}]}]}], ",", "\[IndentingNewLine]", RowBox[{"y2", "==", RowBox[{"x2", "/", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["\[Alpha]", "AC"], "/", SubscriptBox["\[Alpha]", "BC"]}], "x1"}], "+", "x2", "+", RowBox[{ RowBox[{"1", "/", SubscriptBox["\[Alpha]", "BC"]}], " ", "x3"}]}], ")"}]}]}], ",", "\[IndentingNewLine]", RowBox[{"y3", "==", RowBox[{"x3", "/", RowBox[{"(", RowBox[{ RowBox[{ SubscriptBox["\[Alpha]", "AC"], " ", "x1"}], "+", RowBox[{ SubscriptBox["\[Alpha]", "BC"], " ", "x2"}], "+", "x3"}], ")"}]}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"x3", ",", RowBox[{"x", "[", RowBox[{"3", ",", RowBox[{"t", "-", "1"}]}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"x1", ",", RowBox[{"x", "[", RowBox[{"1", ",", RowBox[{"t", "-", "1"}]}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"x2", ",", RowBox[{"x", "[", RowBox[{"2", ",", RowBox[{"t", "-", "1"}]}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"y1", ",", RowBox[{"y", "[", RowBox[{"1", ",", RowBox[{"t", "-", "1"}]}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"y2", ",", RowBox[{"y", "[", RowBox[{"2", ",", RowBox[{"t", "-", "1"}]}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"y3", ",", RowBox[{"y", "[", RowBox[{"3", ",", RowBox[{"t", "-", "1"}]}], "]"}]}], "}"}], ",", RowBox[{"MaxIterations", "\[Rule]", "100000"}]}], "]"}], "]"}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"x", "[", RowBox[{"3", ",", "t"}], "]"}], "=", RowBox[{"x3", "/.", "sol"}]}], ",", RowBox[{ RowBox[{"x", "[", RowBox[{"2", ",", "t"}], "]"}], "=", RowBox[{"x2", "/.", "sol"}]}], ",", RowBox[{ RowBox[{"x", "[", RowBox[{"1", ",", "t"}], "]"}], "=", RowBox[{"x1", "/.", "sol"}]}], ",", RowBox[{ RowBox[{"y", "[", RowBox[{"3", ",", "t"}], "]"}], "=", RowBox[{"y3", "/.", "sol"}]}], ",", RowBox[{ RowBox[{"y", "[", RowBox[{"2", ",", "t"}], "]"}], "=", RowBox[{"y2", "/.", "sol"}]}], ",", RowBox[{ RowBox[{"y", "[", RowBox[{"1", ",", "t"}], "]"}], "=", RowBox[{"y1", "/.", "sol"}]}], ",", "\[IndentingNewLine]", RowBox[{"t", "++"}], ",", RowBox[{"Print", "[", "t", "]"}]}], "}"}]}], "]"}], "//", "Timing"}], "\[IndentingNewLine]"}]], "Input", CellChangeTimes->{ 3.4100917219922276`*^9, {3.4100918343938527`*^9, 3.4100919867729635`*^9}, { 3.410092098133091*^9, 3.410092190515931*^9}, {3.4100922333074627`*^9, 3.410092299903222*^9}, {3.410092360229968*^9, 3.4100924322435184`*^9}, { 3.4100924634684176`*^9, 3.410092518737891*^9}, {3.4100925494921136`*^9, 3.410092579335026*^9}, {3.410092618471301*^9, 3.4100926208847713`*^9}, { 3.410092705075832*^9, 3.4100927353894205`*^9}, {3.410092807382942*^9, 3.4100928089051313`*^9}, {3.4100929098803263`*^9, 3.4100929207659793`*^9}, {3.410092969736395*^9, 3.4100929716491456`*^9}, { 3.4100930091731024`*^9, 3.4100930574425106`*^9}, {3.4100930906102033`*^9, 3.4100931742404575`*^9}, {3.41009326282784*^9, 3.4100932847794046`*^9}, { 3.410093450117149*^9, 3.4100934508081427`*^9}, {3.4101726245732064`*^9, 3.4101727053894143`*^9}, {3.410172777883656*^9, 3.4101727890196686`*^9}, { 3.4101728699760785`*^9, 3.410172894451272*^9}, {3.410173078485901*^9, 3.4101731062658463`*^9}, {3.410174053517928*^9, 3.410174064153221*^9}, { 3.4101758345489244`*^9, 3.4101758386848717`*^9}, {3.410176028988515*^9, 3.4101760399142256`*^9}, 3.410183978689624*^9, {3.410184132751154*^9, 3.4101843189689217`*^9}, {3.4101845461355715`*^9, 3.410184592091653*^9}, { 3.4101846495242367`*^9, 3.41018465171739*^9}, {3.4101847762664833`*^9, 3.410184782265109*^9}}], Cell[CellGroupData[{ Cell[BoxData["1"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.410185299639056*^9}], Cell[BoxData["2"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101852997291856`*^9}], Cell[BoxData["3"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.410185299809301*^9}], Cell[BoxData["4"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.410185299839344*^9}], Cell[BoxData["5"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101852998994303`*^9}], Cell[BoxData["6"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101852999595165`*^9}], Cell[BoxData["7"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853000196033`*^9}], Cell[BoxData["8"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853000696754`*^9}], Cell[BoxData["9"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.410185300119747*^9}], Cell[BoxData["10"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.410185300189848*^9}], Cell[BoxData["11"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853002299056`*^9}], Cell[BoxData["12"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853002899923`*^9}], Cell[BoxData["13"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853003500786`*^9}], Cell[BoxData["14"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.410185300390136*^9}], Cell[BoxData["15"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853004502225`*^9}], Cell[BoxData["16"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853005002947`*^9}], Cell[BoxData["17"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853005703955`*^9}], Cell[BoxData["18"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853006104527`*^9}], Cell[BoxData["19"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.410185300670539*^9}], Cell[BoxData["20"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853007406397`*^9}], Cell[BoxData["21"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.410185300780698*^9}], Cell[BoxData["22"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.410185300840784*^9}], Cell[BoxData["23"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.410185300910885*^9}], Cell[BoxData["24"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853009509425`*^9}], Cell[BoxData["25"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.410185301011029*^9}], Cell[BoxData["26"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.410185301061101*^9}], Cell[BoxData["27"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853011312017`*^9}], Cell[BoxData["28"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853011812735`*^9}], Cell[BoxData["29"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853012413597`*^9}], Cell[BoxData["30"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853013114605`*^9}], Cell[BoxData["31"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853013515186`*^9}], Cell[BoxData["32"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853014015903`*^9}], Cell[BoxData["33"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853014616766`*^9}], Cell[BoxData["34"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853015317774`*^9}], Cell[BoxData["35"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853015818496`*^9}], Cell[BoxData["36"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853016419363`*^9}], Cell[BoxData["37"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.410185301712037*^9}], Cell[BoxData["38"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853017520943`*^9}], Cell[BoxData["39"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853018121805`*^9}], Cell[BoxData["40"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853018722672`*^9}], Cell[BoxData["41"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853019323535`*^9}], Cell[BoxData["42"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853019824257`*^9}], Cell[BoxData["43"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.410185302042512*^9}], Cell[BoxData["44"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853021126127`*^9}], Cell[BoxData["45"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.410185302162685*^9}], Cell[BoxData["46"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.410185302212757*^9}], Cell[BoxData["47"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.410185302282858*^9}], Cell[BoxData["48"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853023329296`*^9}], Cell[BoxData["49"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853023830013`*^9}], Cell[BoxData["50"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853024631166`*^9}], Cell[BoxData["51"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.410185302513189*^9}], Cell[BoxData["52"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.410185302573275*^9}], Cell[BoxData["53"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853026533904`*^9}], Cell[BoxData["54"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.410185302693448*^9}], Cell[BoxData["55"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.41018530274352*^9}], Cell[BoxData["56"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.410185302793592*^9}], Cell[BoxData["57"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.410185302843664*^9}], Cell[BoxData["58"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853028837214`*^9}], Cell[BoxData["59"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853029337935`*^9}], Cell[BoxData["60"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101853029838657`*^9}] }, Open ]], Cell[BoxData[ RowBox[{"{", RowBox[{"0.831`", ",", "Null"}], "}"}]], "Output", CellChangeTimes->{{3.4101731015390496`*^9, 3.4101731109025135`*^9}, 3.4101740866555777`*^9, 3.410174233276408*^9, 3.4101747868423967`*^9, 3.4101755031123424`*^9, {3.410175916616933*^9, 3.4101759312980433`*^9}, 3.410176057319253*^9, 3.4101775787469587`*^9, 3.410184340029205*^9, 3.410184661972136*^9, 3.410184799990597*^9, 3.4101853031541104`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"plt1", "=", RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"y", "[", RowBox[{"1", ",", "t"}], "]"}], ",", RowBox[{"y", "[", RowBox[{"2", ",", "t"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"t", ",", RowBox[{"-", "1"}], ",", "59"}], "}"}]}], "]"}], ",", RowBox[{"PlotRange", "\[Rule]", " ", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "1"}], "}"}]}], "}"}]}], ",", RowBox[{"AspectRatio", "\[Rule]", " ", "1"}], ",", RowBox[{"Epilog", "\[Rule]", " ", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "1"}], "}"}]}], "}"}], "]"}]}]}], "]"}]}]], "Input",\ CellChangeTimes->{{3.4100935088315763`*^9, 3.410093509402397*^9}, { 3.41017311487823*^9, 3.4101731172015715`*^9}, {3.4101742461549263`*^9, 3.410174247476827*^9}, {3.41017469207613*^9, 3.410174725203765*^9}}], Cell[BoxData[ GraphicsBox[{{}, {Hue[0.67, 0.6, 0.6], PointBox[{{0.4, 0.4}, {0.4949929607977413, 0.40493887514997967`}, { 0.571779221101841, 0.3718403417446142}, {0.640695401182614, 0.3205436054170047}, {0.7061137160449863, 0.2630089307360585}, { 0.7680020134457367, 0.206096077266995}, {0.8242767908674061, 0.15425625431291942`}, {0.8725480479417976, 0.11038775989849024`}, { 0.9113122625938601, 0.07580674072133686}, {0.9404831958705043, 0.050276476649935314`}, {0.9612073114825679, 0.03245201484311536}, { 0.9752613041681355, 0.02053950986231399}, {0.9844653509918122, 0.012826581500989545`}, {0.9903463560839156, 0.007939889756300202}, { 0.9940418968825027, 0.004887521284135303}, {0.9963388429941256, 0.00299809449749849}, {0.9977565056968245, 0.0018351213447786714`}, { 0.9986276035430713, 0.0011217841925089207`}, {0.9991613749355367, 0.0006851764458319741}, {0.9994878840289279, 0.0004182935206446715}, { 0.9996873979253934, 0.00025528727077556336`}, {0.9998092314178021, 0.00015577492186401236`}, {0.9998835993565666, 0.00009504239931458159}, { 0.999928982847548, 0.000057983931994224556`}, {0.9999566742308854, 0.00003537365237695732}, {0.999973568967401, 0.00002157948884136345}, { 0.9999838760028441, 0.000013164237255051422`}, {0.9999901638387704, 8.030566108085357*^-6}, {0.9999939996703204, 4.8988498657298065`*^-6}, { 0.9999963396511315, 2.988412732476743*^-6}, {0.9999977671036588, 1.8229976093411578`*^-6}, {0.9999986378847354, 1.1120672589666394`*^-6}, {0.9999991690810198, 6.783841074829269*^-7}, { 0.9999994931222744, 4.138281840408347*^-7}, {0.9999996907942534, 2.5244358380909236`*^-7}, {0.999999811378238, 1.5399567737212213`*^-7}, { 0.999999884936926, 9.394045863001603*^-8}, {0.9999999298092139, 5.730556419620167*^-8}, {0.9999999571822143, 3.495754118987411*^-8}, { 0.9999999738802935, 2.1324799591785054`*^-8}, {0.9999999840664565, 1.3008554337424276`*^-8}, {0.9999999902802196, 7.935478331485382*^-9}, { 0.9999999940707398, 4.840800484844209*^-9}, {0.9999999963830326, 2.952985101468699*^-9}, {0.9999999977935776, 1.8013799630108517`*^-9}, { 0.9999999986540382, 1.0988777989269314`*^-9}, {0.9999999991789364, 6.703374308476096*^-10}, {0.9999999994991349, 4.089192369701577*^-10}, { 0.9999999997307184, 2.393555006091168*^-10}, {0.9999999998423796, 1.401035962088464*^-10}, {0.999999999907739, 0.}, {1., 0.}, {1., 0.}, { 1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}}]}, {}}, AspectRatio->1, Axes->True, AxesOrigin->{0, 0}, Epilog->LineBox[{{1, 0}, {0, 1}}], PlotRange->{{0, 1}, {0, 1}}, PlotRangeClipping->True, PlotRangePadding->{Automatic, Automatic}]], "Output", CellChangeTimes->{{3.4100934996784143`*^9, 3.4100935105340242`*^9}, 3.410173119524912*^9, 3.4101740879774785`*^9, {3.4101742345582514`*^9, 3.4101742481077347`*^9}, 3.410174726355421*^9, 3.4101747884346867`*^9, 3.4101755047947617`*^9, 3.410175932539829*^9, 3.4101760590116863`*^9, 3.4101775792977505`*^9, 3.4101843425328045`*^9, 3.4101846631838784`*^9, 3.4101848011923246`*^9, 3.4101853047764435`*^9}] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["\<\ Stripping Cascade Trajectory at Damk\[ODoubleDot]hler number Da=0.277\ \>", "Subsubtitle", CellChangeTimes->{ 3.410091630290366*^9, {3.410092930199544*^9, 3.4100929425773425`*^9}, { 3.410093659297936*^9, 3.4100937216576047`*^9}, {3.4100937524218416`*^9, 3.4100937804120893`*^9}, {3.410093953370792*^9, 3.410093953751339*^9}, { 3.410173188454027*^9, 3.4101731977473907`*^9}}, Background->RGBColor[0.996109, 0.996109, 0.640635]], Cell[BoxData[ RowBox[{ RowBox[{"t", "=", "0"}], ";", RowBox[{ RowBox[{"y", "[", RowBox[{"1", ",", RowBox[{"-", "1"}]}], "]"}], "=", "0.4"}], ";", RowBox[{ RowBox[{"y", "[", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "]"}], "=", "0.4"}], ";", RowBox[{ RowBox[{"y", "[", RowBox[{"3", ",", RowBox[{"-", "1"}]}], "]"}], "=", "0.2"}], ";", RowBox[{ RowBox[{"x", "[", RowBox[{"1", ",", RowBox[{"-", "1"}]}], "]"}], "=", "0.4"}], ";", RowBox[{ RowBox[{"x", "[", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "]"}], "=", "0.4"}], ";", RowBox[{ RowBox[{"x", "[", RowBox[{"3", ",", RowBox[{"-", "1"}]}], "]"}], "=", "0.2"}], ";"}]], "Input", CellChangeTimes->{{3.410091642257574*^9, 3.4100917039863358`*^9}, { 3.4100920383070655`*^9, 3.410092045727736*^9}, {3.410092192739128*^9, 3.4100921996590786`*^9}, {3.4100926628350925`*^9, 3.4100926712471886`*^9}, { 3.4101724027141886`*^9, 3.4101724283210096`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"\[Phi]", "=", "0.5"}], ";", RowBox[{"Da", "=", "0.277"}], ";", " "}]], "Input", CellChangeTimes->{{3.41009177143332*^9, 3.410091803018738*^9}, { 3.410172437514229*^9, 3.410172446336915*^9}, {3.410172602381296*^9, 3.410172613166805*^9}, {3.4101739474754467`*^9, 3.4101739505698957`*^9}, { 3.410174126793293*^9, 3.410174129296893*^9}, 3.4101746810302467`*^9, 3.410175306469584*^9, 3.410175354338416*^9, {3.410175941632904*^9, 3.4101759417130194`*^9}, {3.4101760642492175`*^9, 3.410176064379405*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"While", "[", RowBox[{ RowBox[{"t", "<", "60"}], ",", "\[IndentingNewLine]", RowBox[{"{", "\[IndentingNewLine]", RowBox[{ RowBox[{"sol", "=", RowBox[{"Chop", "[", RowBox[{"FindRoot", "[", RowBox[{ RowBox[{"{", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"x", "[", RowBox[{"1", ",", RowBox[{"t", "-", "1"}]}], "]"}], "\[Equal]", " ", RowBox[{ RowBox[{"\[Phi]", " ", FractionBox[ RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "3"], RowBox[{ RowBox[{"x", "[", RowBox[{"i", ",", RowBox[{"t", "-", "1"}]}], "]"}], " ", RowBox[{"M", "[", "i", "]"}]}]}], RowBox[{ RowBox[{"y1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"y2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"y3", " ", RowBox[{"M", "[", "3", "]"}]}]}]], "y1"}], "+", RowBox[{ RowBox[{"(", RowBox[{"1", "-", RowBox[{"\[Phi]", FractionBox[ RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "3"], RowBox[{ RowBox[{"x", "[", RowBox[{"i", ",", RowBox[{"t", "-", "1"}]}], "]"}], " ", RowBox[{"M", "[", "i", "]"}]}]}], RowBox[{ RowBox[{"y1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"y2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"y3", " ", RowBox[{"M", "[", "3", "]"}]}]}]]}]}], ")"}], " ", "x1"}], "-", RowBox[{"Da", " ", FractionBox[ RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "3"], RowBox[{ RowBox[{"x", "[", RowBox[{"i", ",", RowBox[{"t", "-", "1"}]}], "]"}], " ", RowBox[{"M", "[", "i", "]"}]}]}], RowBox[{ RowBox[{"x1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"x2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"x3", " ", RowBox[{"M", "[", "3", "]"}]}]}]], RowBox[{"(", RowBox[{ RowBox[{"x1", " ", "x2"}], "-", RowBox[{"x3", "/", "Keq"}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "x1"}], ")"}]}]}]}], ",", RowBox[{ RowBox[{"x", "[", RowBox[{"2", ",", RowBox[{"t", "-", "1"}]}], "]"}], "\[Equal]", " ", RowBox[{ RowBox[{"\[Phi]", " ", FractionBox[ RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "3"], RowBox[{ RowBox[{"x", "[", RowBox[{"i", ",", RowBox[{"t", "-", "1"}]}], "]"}], " ", RowBox[{"M", "[", "i", "]"}]}]}], RowBox[{ RowBox[{"y1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"y2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"y3", " ", RowBox[{"M", "[", "3", "]"}]}]}]], "y2"}], "+", RowBox[{ RowBox[{"(", RowBox[{"1", "-", RowBox[{"\[Phi]", FractionBox[ RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "3"], RowBox[{ RowBox[{"x", "[", RowBox[{"i", ",", RowBox[{"t", "-", "1"}]}], "]"}], " ", RowBox[{"M", "[", "i", "]"}]}]}], RowBox[{ RowBox[{"y1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"y2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"y3", " ", RowBox[{"M", "[", "3", "]"}]}]}]]}]}], ")"}], " ", "x2"}], "-", RowBox[{"Da", " ", FractionBox[ RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "3"], RowBox[{ RowBox[{"x", "[", RowBox[{"i", ",", RowBox[{"t", "-", "1"}]}], "]"}], " ", RowBox[{"M", "[", "i", "]"}]}]}], RowBox[{ RowBox[{"x1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"x2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"x3", " ", RowBox[{"M", "[", "3", "]"}]}]}]], RowBox[{"(", RowBox[{ RowBox[{"x1", " ", "x2"}], "-", RowBox[{"x3", "/", "Keq"}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "x2"}], ")"}]}]}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"x1", "+", "x2", "+", "x3"}], "==", "1"}], ",", "\[IndentingNewLine]", RowBox[{"y1", "==", RowBox[{"x1", "/", RowBox[{"(", RowBox[{"x1", "+", RowBox[{ RowBox[{ SubscriptBox["\[Alpha]", "BC"], "/", SubscriptBox["\[Alpha]", "AC"]}], " ", "x2"}], "+", RowBox[{ RowBox[{"1", "/", SubscriptBox["\[Alpha]", "AC"]}], " ", "x3"}]}], ")"}]}]}], ",", "\[IndentingNewLine]", RowBox[{"y2", "==", RowBox[{"x2", "/", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["\[Alpha]", "AC"], "/", SubscriptBox["\[Alpha]", "BC"]}], "x1"}], "+", "x2", "+", RowBox[{ RowBox[{"1", "/", SubscriptBox["\[Alpha]", "BC"]}], " ", "x3"}]}], ")"}]}]}], ",", "\[IndentingNewLine]", RowBox[{"y3", "==", RowBox[{"x3", "/", RowBox[{"(", RowBox[{ RowBox[{ SubscriptBox["\[Alpha]", "AC"], " ", "x1"}], "+", RowBox[{ SubscriptBox["\[Alpha]", "BC"], " ", "x2"}], "+", "x3"}], ")"}]}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"x3", ",", RowBox[{"x", "[", RowBox[{"3", ",", RowBox[{"t", "-", "1"}]}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"x1", ",", RowBox[{"x", "[", RowBox[{"1", ",", RowBox[{"t", "-", "1"}]}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"x2", ",", RowBox[{"x", "[", RowBox[{"2", ",", RowBox[{"t", "-", "1"}]}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"y1", ",", RowBox[{"y", "[", RowBox[{"1", ",", RowBox[{"t", "-", "1"}]}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"y2", ",", RowBox[{"y", "[", RowBox[{"2", ",", RowBox[{"t", "-", "1"}]}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"y3", ",", RowBox[{"y", "[", RowBox[{"3", ",", RowBox[{"t", "-", "1"}]}], "]"}]}], "}"}], ",", RowBox[{"MaxIterations", "\[Rule]", "1000000"}]}], "]"}], "]"}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"x", "[", RowBox[{"3", ",", "t"}], "]"}], "=", RowBox[{"x3", "/.", "sol"}]}], ",", RowBox[{ RowBox[{"x", "[", RowBox[{"2", ",", "t"}], "]"}], "=", RowBox[{"x2", "/.", "sol"}]}], ",", RowBox[{ RowBox[{"x", "[", RowBox[{"1", ",", "t"}], "]"}], "=", RowBox[{"x1", "/.", "sol"}]}], ",", RowBox[{ RowBox[{"y", "[", RowBox[{"3", ",", "t"}], "]"}], "=", RowBox[{"y3", "/.", "sol"}]}], ",", RowBox[{ RowBox[{"y", "[", RowBox[{"2", ",", "t"}], "]"}], "=", RowBox[{"y2", "/.", "sol"}]}], ",", RowBox[{ RowBox[{"y", "[", RowBox[{"1", ",", "t"}], "]"}], "=", RowBox[{"y1", "/.", "sol"}]}], ",", "\[IndentingNewLine]", RowBox[{"t", "++"}], ",", RowBox[{"Print", "[", "t", "]"}]}], "}"}]}], "]"}], "//", "Timing"}], "\[IndentingNewLine]", "\[IndentingNewLine]"}]], "Input", CellChangeTimes->{ 3.4100917219922276`*^9, {3.4100918343938527`*^9, 3.4100919867729635`*^9}, { 3.410092098133091*^9, 3.410092190515931*^9}, {3.4100922333074627`*^9, 3.410092299903222*^9}, {3.410092360229968*^9, 3.4100924322435184`*^9}, { 3.4100924634684176`*^9, 3.410092518737891*^9}, {3.4100925494921136`*^9, 3.410092579335026*^9}, {3.410092618471301*^9, 3.4100926208847713`*^9}, { 3.410092705075832*^9, 3.4100927353894205`*^9}, {3.410092807382942*^9, 3.4100928089051313`*^9}, {3.4100929098803263`*^9, 3.4100929207659793`*^9}, {3.410092969736395*^9, 3.4100929716491456`*^9}, { 3.4100930091731024`*^9, 3.4100930574425106`*^9}, {3.4100930906102033`*^9, 3.4100931742404575`*^9}, {3.41009326282784*^9, 3.4100932847794046`*^9}, { 3.410093450117149*^9, 3.4100934508081427`*^9}, {3.4101731565181055`*^9, 3.410173158531*^9}, {3.4101732281310797`*^9, 3.410173245936683*^9}, 3.410173288197451*^9, {3.4101734446524224`*^9, 3.410173457360696*^9}, { 3.410173575921178*^9, 3.4101735769426465`*^9}, {3.410173748779736*^9, 3.410173756420723*^9}, {3.410174117439843*^9, 3.4101741343441505`*^9}, { 3.410174523293432*^9, 3.4101745339687824`*^9}, {3.4101758444832096`*^9, 3.410175846015413*^9}, {3.410176071319384*^9, 3.4101760755755043`*^9}, { 3.410184352777536*^9, 3.4101844114919634`*^9}, {3.4101846812999277`*^9, 3.4101846977335587`*^9}, {3.410184737981432*^9, 3.4101847431388483`*^9}}], Cell[CellGroupData[{ Cell[BoxData["1"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.410185308631987*^9}], Cell[BoxData["2"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.4101853087221165`*^9}], Cell[BoxData["3"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.410185308802232*^9}], Cell[BoxData["4"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.41018530884229*^9}], Cell[BoxData["5"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.4101853088923616`*^9}], Cell[BoxData["6"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.410185308952448*^9}], Cell[BoxData["7"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.41018530900252*^9}], Cell[BoxData["8"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.4101853090626063`*^9}], Cell[BoxData["9"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.4101853091126785`*^9}], Cell[BoxData["10"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.410185309172765*^9}], Cell[BoxData["11"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.4101853092128224`*^9}], Cell[BoxData["12"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.4101853092729087`*^9}], Cell[BoxData["13"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.410185309322981*^9}], Cell[BoxData["14"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.4101853093930817`*^9}], Cell[BoxData["15"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.410185309433139*^9}], Cell[BoxData["16"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.4101853094932256`*^9}], Cell[BoxData["17"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.4101853095533123`*^9}], Cell[BoxData["18"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.4101853095933695`*^9}], Cell[BoxData["19"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.4101853096534557`*^9}], Cell[BoxData["20"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.4101853097135425`*^9}], Cell[BoxData["21"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.4101853097636147`*^9}], Cell[BoxData["22"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.410185309823701*^9}], Cell[BoxData["23"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.410185309883787*^9}], Cell[BoxData["24"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.4101853099438734`*^9}], Cell[BoxData["25"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.4101853100039597`*^9}], Cell[BoxData["26"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.4101853100540323`*^9}], Cell[BoxData["27"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.410185310124133*^9}], Cell[BoxData["28"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.410185310174205*^9}], Cell[BoxData["29"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.4101853102242765`*^9}], Cell[BoxData["30"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.4101853102843633`*^9}], Cell[BoxData["31"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.410185310354464*^9}], Cell[BoxData["32"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.4101853104045362`*^9}], Cell[BoxData["33"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.410185310454608*^9}], Cell[BoxData["34"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.4101853105247087`*^9}], Cell[BoxData["35"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.4101853105747805`*^9}], Cell[BoxData["36"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.410185310624853*^9}], Cell[BoxData["37"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.4101853106849394`*^9}], Cell[BoxData["38"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.41018531075504*^9}], Cell[BoxData["39"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.410185310805112*^9}], Cell[BoxData["40"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.410185310865198*^9}], Cell[BoxData["41"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.4101853109152703`*^9}], Cell[BoxData["42"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.410185310975357*^9}], Cell[BoxData["43"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.4101853110354433`*^9}], Cell[BoxData["44"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.410185311105544*^9}], Cell[BoxData["45"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.410185311155616*^9}], Cell[BoxData["46"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.410185311215702*^9}], Cell[BoxData["47"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.4101853112657747`*^9}], Cell[BoxData["48"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.410185311325861*^9}], Cell[BoxData["49"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.410185311385947*^9}], Cell[BoxData["50"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.410185311456048*^9}], Cell[BoxData["51"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.4101853114961057`*^9}], Cell[BoxData["52"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.410185311556192*^9}], Cell[BoxData["53"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.410185311606264*^9}], Cell[BoxData["54"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.410185311676365*^9}], Cell[BoxData["55"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.4101853117264366`*^9}], Cell[BoxData["56"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.410185311776509*^9}], Cell[BoxData["57"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.4101853118365955`*^9}], Cell[BoxData["58"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.4101853119066963`*^9}], Cell[BoxData["59"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.410185311956768*^9}], Cell[BoxData["60"], "Print", CellChangeTimes->{3.4101739745243406`*^9, 3.4101741348949423`*^9, 3.410174538385133*^9, 3.410174795965515*^9, 3.4101753078916287`*^9, 3.4101753592855296`*^9, 3.4101755092010975`*^9, 3.410175943585712*^9, 3.410176077007563*^9, 3.4101775833535824`*^9, 3.410184418001323*^9, 3.410184699425992*^9, 3.4101847473549104`*^9, 3.4101853120068398`*^9}] }, Open ]], Cell[BoxData[ RowBox[{"{", RowBox[{"0.7709999999999999`", ",", "Null"}], "}"}]], "Output", CellChangeTimes->{ 3.410092595498267*^9, 3.4100926298777027`*^9, 3.4100926750426464`*^9, { 3.4100927314036894`*^9, 3.410092743901661*^9}, 3.410092813902317*^9, { 3.4100929581297054`*^9, 3.410092977137037*^9}, 3.410093190403699*^9, 3.410093308002798*^9, 3.410093481812725*^9, 3.410093587644904*^9, { 3.4101732505132637`*^9, 3.4101732623903427`*^9}, 3.410173300154645*^9, 3.4101734225806847`*^9, 3.410173465772792*^9, 3.4101735879584866`*^9, 3.4101737734952755`*^9, {3.4101739684856577`*^9, 3.4101739846889567`*^9}, 3.4101741378692193`*^9, 3.4101745413093376`*^9, 3.410174798909749*^9, 3.410175310895949*^9, 3.410175362319893*^9, 3.410175512215432*^9, 3.4101759470006227`*^9, 3.410176080182128*^9, 3.410177586678363*^9, 3.410184421396205*^9, 3.410184703001133*^9, 3.4101847508899937`*^9, 3.410185312056912*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"plt2", "=", RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"x", "[", RowBox[{"1", ",", "t"}], "]"}], ",", RowBox[{"x", "[", RowBox[{"2", ",", "t"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"t", ",", RowBox[{"-", "1"}], ",", "59"}], "}"}]}], "]"}], ",", RowBox[{"PlotRange", "\[Rule]", " ", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "1"}], "}"}]}], "}"}]}], ",", RowBox[{"AspectRatio", "\[Rule]", " ", "1"}], ",", RowBox[{"PlotStyle", "\[Rule]", " ", "Red"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.4100935088315763`*^9, 3.410093509402397*^9}, { 3.410173302958677*^9, 3.410173305141816*^9}, {3.410173992710491*^9, 3.410173995113947*^9}, {3.4101741421053104`*^9, 3.410174143096736*^9}, { 3.4101742191460896`*^9, 3.4101742206482496`*^9}, {3.410174738482859*^9, 3.4101747403455377`*^9}, {3.4101748052789073`*^9, 3.4101748126995773`*^9}}],\ Cell[BoxData[ GraphicsBox[{{}, {RGBColor[1, 0, 0], PointBox[{{0.4, 0.4}, {0.29636181594182454`, 0.4040745016908241}, { 0.21167269559066593`, 0.3909182765295349}, {0.1523136868921553, 0.3621884116897882}, {0.1158784112102993, 0.32523548881175873`}, { 0.095726290812115, 0.288013894380049}, {0.08539515687788057, 0.25565260376295534`}, {0.08036620662609187, 0.23006231046571013`}, { 0.07799627052036644, 0.21104886666571498`}, {0.0768979556017048, 0.1974918178785205}, {0.07639071109827053, 0.18808488354511366`}, { 0.07615441420618946, 0.18167393001779214`}, {0.07604209622801274, 0.17735629714042697`}, {0.07598703506844022, 0.17447110856925221`}, { 0.07595895090137143, 0.17255302299481026`}, {0.07594396867631163, 0.17128218041539595`}, {0.07593560409041895, 0.17044204525056525`}, { 0.07593073565185665, 0.16988745582847625`}, {0.0759278016361807, 0.16952171207359018`}, {0.0759259849103766, 0.1692806613797533}, { 0.07592483745649234, 0.16912185788192674`}, {0.07592410254612789, 0.16901726702021833`}, {0.07592362736915027, 0.1689483938641411}, { 0.07592331818102135, 0.16890304614746393`}, {0.07592311616184715, 0.16887319042825677`}, {0.07592298380921039, 0.16885353521357244`}, { 0.07592289694787095, 0.16884059582630265`}, {0.07592283987847237, 0.16883207777538617`}, {0.07592280235631227, 0.16882647038761356`}, { 0.07592277767503038, 0.16882277910829158`}, {0.07592276143558867, 0.16882034919608688`}, {0.07592275074867631, 0.1688187496287754}, { 0.07592274371500331, 0.16881769666526325`}, {0.07592273908541314, 0.1688170035214051}, {0.07592273603806671, 0.16881654723981201`}, { 0.07592273403215047, 0.16881624687973992`}, {0.07592273271173404, 0.16881604915944862`}, {0.07592273184254669, 0.16881591900466136`}, { 0.07592273127038537, 0.16881583332673197`}, {0.07592273089374649, 0.16881577692691183`}, {0.07592273064581445, 0.16881573980020334`}, { 0.07592273048260666, 0.16881571536054485`}, {0.07592273037517079, 0.1688156992724801}, {0.07592273030444822, 0.16881568868207755`}, { 0.07592273025789321, 0.16881568171065958`}, {0.0759227302272471, 0.16881567712153572`}, {0.07592273020707346, 0.16881567410062126`}, { 0.07592273019379366, 0.16881567211202286`}, {0.07592273018505184, 0.16881567080297438`}, {0.07592273017929731, 0.16881566994125793`}, { 0.07592273017550923, 0.16881566937400985`}, {0.07592273017301561, 0.16881566900060355`}, {0.07592273017137416, 0.16881566875479884`}, { 0.0759227301702936, 0.1688156685929912}, {0.0759227301695823, 0.16881566848647697`}, {0.07592273016911406, 0.16881566841636111`}, { 0.07592273016880585, 0.1688156683702055}, {0.07592273016860294, 0.16881566833982226`}, {0.07592273016846939, 0.16881566831982167`}, { 0.07592273016838147, 0.1688156683066558}, {0.07592273016832359, 0.168815668297989}}]}, {}}, AspectRatio->1, Axes->True, AxesOrigin->{0, 0}, PlotRange->{{0, 1}, {0, 1}}, PlotRangeClipping->True, PlotRangePadding->{Automatic, Automatic}]], "Output", CellChangeTimes->{{3.4100934996784143`*^9, 3.4100935105340242`*^9}, 3.4100935905390654`*^9, 3.410173306593904*^9, 3.410173468186262*^9, 3.4101735894906893`*^9, 3.41017377926357*^9, {3.410173986171088*^9, 3.4101739962756176`*^9}, 3.410174144598896*^9, 3.4101745426913247`*^9, 3.410174814492155*^9, 3.4101753121177053`*^9, 3.4101753630910015`*^9, 3.410175513937909*^9, 3.4101759483125086`*^9, 3.410176082285152*^9, 3.410177587850048*^9, 3.410184422808235*^9, 3.4101847045233216`*^9, 3.41018475220188*^9, 3.410185313308712*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"x", "[", RowBox[{"1", ",", "59"}], "]"}], ",", RowBox[{"x", "[", RowBox[{"2", ",", "59"}], "]"}]}], "}"}]], "Input", CellChangeTimes->{{3.4101734813251553`*^9, 3.4101735363943405`*^9}, { 3.4101737905197554`*^9, 3.4101737958574305`*^9}, {3.410174170205717*^9, 3.4101742096324096`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"0.07592273016832359`", ",", "0.168815668297989`"}], "}"}]], "Output", CellChangeTimes->{{3.4101735008131776`*^9, 3.410173536975176*^9}, 3.410173796588482*^9, 3.410174176685034*^9, 3.410174210313389*^9, 3.4101745473079634`*^9, 3.410174820250435*^9, 3.410175315262227*^9, 3.4101753651339393`*^9, 3.410175515640357*^9, 3.4101759555529203`*^9, 3.4101760841578445`*^9, 3.410177589181963*^9, 3.4101844265436063`*^9, 3.4101847078881598`*^9, 3.4101847539944572`*^9, 3.4101853147608*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"FindRoot", "[", RowBox[{ RowBox[{"{", "\[IndentingNewLine]", RowBox[{ RowBox[{"x1", "\[Equal]", " ", RowBox[{ RowBox[{"\[Phi]", " ", FractionBox[ RowBox[{ RowBox[{"x1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"x2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"x3", " ", RowBox[{"M", "[", "3", "]"}]}]}], RowBox[{ RowBox[{"y1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"y2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"y3", " ", RowBox[{"M", "[", "3", "]"}]}]}]], "y1"}], "+", RowBox[{ RowBox[{"(", RowBox[{"1", "-", RowBox[{"\[Phi]", FractionBox[ RowBox[{ RowBox[{"x1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"x2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"x3", " ", RowBox[{"M", "[", "3", "]"}]}]}], RowBox[{ RowBox[{"y1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"y2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"y3", " ", RowBox[{"M", "[", "3", "]"}]}]}]]}]}], ")"}], " ", "x1"}], "-", RowBox[{"Da", " ", FractionBox[ RowBox[{ RowBox[{"x1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"x2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"x3", " ", RowBox[{"M", "[", "3", "]"}]}]}], RowBox[{ RowBox[{"x1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"x2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"x3", " ", RowBox[{"M", "[", "3", "]"}]}]}]], RowBox[{"(", RowBox[{ RowBox[{"x1", " ", "x2"}], "-", RowBox[{"x3", "/", "Keq"}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "x1"}], ")"}]}]}]}], ",", RowBox[{"x2", "\[Equal]", " ", RowBox[{ RowBox[{"\[Phi]", " ", FractionBox[ RowBox[{ RowBox[{"x1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"x2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"x3", " ", RowBox[{"M", "[", "3", "]"}]}]}], RowBox[{ RowBox[{"y1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"y2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"y3", " ", RowBox[{"M", "[", "3", "]"}]}]}]], "y2"}], "+", RowBox[{ RowBox[{"(", RowBox[{"1", "-", RowBox[{"\[Phi]", FractionBox[ RowBox[{ RowBox[{"x1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"x2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"x3", " ", RowBox[{"M", "[", "3", "]"}]}]}], RowBox[{ RowBox[{"y1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"y2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"y3", " ", RowBox[{"M", "[", "3", "]"}]}]}]]}]}], ")"}], " ", "x2"}], "-", RowBox[{"Da", " ", FractionBox[ RowBox[{ RowBox[{"x1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"x2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"x3", " ", RowBox[{"M", "[", "3", "]"}]}]}], RowBox[{ RowBox[{"x1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"x2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"x3", " ", RowBox[{"M", "[", "3", "]"}]}]}]], RowBox[{"(", RowBox[{ RowBox[{"x1", " ", "x2"}], "-", RowBox[{"x3", "/", "Keq"}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "x2"}], ")"}]}]}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"x1", "+", "x2", "+", "x3"}], "==", "1"}], ",", "\[IndentingNewLine]", RowBox[{"y1", "==", RowBox[{"x1", "/", RowBox[{"(", RowBox[{"x1", "+", RowBox[{ RowBox[{ SubscriptBox["\[Alpha]", "BC"], "/", SubscriptBox["\[Alpha]", "AC"]}], " ", "x2"}], "+", RowBox[{ RowBox[{"1", "/", SubscriptBox["\[Alpha]", "AC"]}], " ", "x3"}]}], ")"}]}]}], ",", "\[IndentingNewLine]", RowBox[{"y2", "==", RowBox[{"x2", "/", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["\[Alpha]", "AC"], "/", SubscriptBox["\[Alpha]", "BC"]}], "x1"}], "+", "x2", "+", RowBox[{ RowBox[{"1", "/", SubscriptBox["\[Alpha]", "BC"]}], " ", "x3"}]}], ")"}]}]}], ",", "\[IndentingNewLine]", RowBox[{"y3", "==", RowBox[{"x3", "/", RowBox[{"(", RowBox[{ RowBox[{ SubscriptBox["\[Alpha]", "AC"], " ", "x1"}], "+", RowBox[{ SubscriptBox["\[Alpha]", "BC"], " ", "x2"}], "+", "x3"}], ")"}]}]}]} ], "}"}], ",", RowBox[{"{", RowBox[{"x3", ",", "0.7"}], "}"}], ",", RowBox[{"{", RowBox[{"x1", ",", "0.1"}], "}"}], ",", RowBox[{"{", RowBox[{"x2", ",", "0.2"}], "}"}], ",", RowBox[{"{", RowBox[{"y1", ",", "0.1"}], "}"}], ",", RowBox[{"{", RowBox[{"y2", ",", "0.3"}], "}"}], ",", RowBox[{"{", RowBox[{"y3", ",", "0.6"}], "}"}], ",", RowBox[{"MaxIterations", "\[Rule]", "1000000"}]}], "]"}]], "Input", CellChangeTimes->{{3.410185074555402*^9, 3.4101850785311184`*^9}, { 3.410185109185197*^9, 3.4101852192134094`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"x3", "\[Rule]", "0.7552616015504943`"}], ",", RowBox[{"x1", "\[Rule]", "0.07592273016821209`"}], ",", RowBox[{"x2", "\[Rule]", "0.16881566828129366`"}], ",", RowBox[{"y1", "\[Rule]", "0.23128526355358356`"}], ",", RowBox[{"y2", "\[Rule]", "0.30856037113449974`"}], ",", RowBox[{"y3", "\[Rule]", "0.46015436531191667`"}]}], "}"}]], "Output", CellChangeTimes->{{3.410185156212819*^9, 3.410185182260274*^9}, 3.410185219974504*^9, 3.41018531626296*^9}] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["\<\ Stripping & Rectifying Profiles at Damk\[ODoubleDot]hler number Da=0.277\ \>", "Subsubtitle", CellChangeTimes->{ 3.410091630290366*^9, {3.410092930199544*^9, 3.4100929425773425`*^9}, { 3.410093659297936*^9, 3.4100937216576047`*^9}, {3.4100937524218416`*^9, 3.4100937804120893`*^9}, {3.410093953370792*^9, 3.410093953751339*^9}, { 3.410173188454027*^9, 3.4101731977473907`*^9}, {3.410174895007931*^9, 3.4101748974714737`*^9}, {3.4101774943756385`*^9, 3.410177501175416*^9}, { 3.410177545008445*^9, 3.4101775472416563`*^9}}, Background->RGBColor[0.996109, 0.996109, 0.640635]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Show", "[", RowBox[{"plt1", ",", "plt2"}], "]"}]], "Input", CellChangeTimes->{{3.4101757384807854`*^9, 3.410175744549512*^9}}], Cell[BoxData[ GraphicsBox[{{{}, {Hue[0.67, 0.6, 0.6], PointBox[{{0.4, 0.4}, {0.4949929607977413, 0.40493887514997967`}, { 0.571779221101841, 0.3718403417446142}, {0.640695401182614, 0.3205436054170047}, {0.7061137160449863, 0.2630089307360585}, { 0.7680020134457367, 0.206096077266995}, {0.8242767908674061, 0.15425625431291942`}, {0.8725480479417976, 0.11038775989849024`}, { 0.9113122625938601, 0.07580674072133686}, {0.9404831958705043, 0.050276476649935314`}, {0.9612073114825679, 0.03245201484311536}, { 0.9752613041681355, 0.02053950986231399}, {0.9844653509918122, 0.012826581500989545`}, {0.9903463560839156, 0.007939889756300202}, { 0.9940418968825027, 0.004887521284135303}, {0.9963388429941256, 0.00299809449749849}, {0.9977565056968245, 0.0018351213447786714`}, { 0.9986276035430713, 0.0011217841925089207`}, {0.9991613749355367, 0.0006851764458319741}, {0.9994878840289279, 0.0004182935206446715}, { 0.9996873979253934, 0.00025528727077556336`}, {0.9998092314178021, 0.00015577492186401236`}, {0.9998835993565666, 0.00009504239931458159}, {0.999928982847548, 0.000057983931994224556`}, {0.9999566742308854, 0.00003537365237695732}, {0.999973568967401, 0.00002157948884136345}, { 0.9999838760028441, 0.000013164237255051422`}, {0.9999901638387704, 8.030566108085357*^-6}, {0.9999939996703204, 4.8988498657298065`*^-6}, { 0.9999963396511315, 2.988412732476743*^-6}, {0.9999977671036588, 1.8229976093411578`*^-6}, {0.9999986378847354, 1.1120672589666394`*^-6}, {0.9999991690810198, 6.783841074829269*^-7}, { 0.9999994931222744, 4.138281840408347*^-7}, {0.9999996907942534, 2.5244358380909236`*^-7}, {0.999999811378238, 1.5399567737212213`*^-7}, {0.999999884936926, 9.394045863001603*^-8}, { 0.9999999298092139, 5.730556419620167*^-8}, {0.9999999571822143, 3.495754118987411*^-8}, {0.9999999738802935, 2.1324799591785054`*^-8}, { 0.9999999840664565, 1.3008554337424276`*^-8}, {0.9999999902802196, 7.935478331485382*^-9}, {0.9999999940707398, 4.840800484844209*^-9}, { 0.9999999963830326, 2.952985101468699*^-9}, {0.9999999977935776, 1.8013799630108517`*^-9}, {0.9999999986540382, 1.0988777989269314`*^-9}, {0.9999999991789364, 6.703374308476096*^-10}, {0.9999999994991349, 4.089192369701577*^-10}, { 0.9999999997307184, 2.393555006091168*^-10}, {0.9999999998423796, 1.401035962088464*^-10}, {0.999999999907739, 0.}, {1., 0.}, {1., 0.}, { 1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, { 1., 0.}}]}, {}}, {{}, {RGBColor[1, 0, 0], PointBox[{{0.4, 0.4}, {0.29636181594182454`, 0.4040745016908241}, { 0.21167269559066593`, 0.3909182765295349}, {0.1523136868921553, 0.3621884116897882}, {0.1158784112102993, 0.32523548881175873`}, { 0.095726290812115, 0.288013894380049}, {0.08539515687788057, 0.25565260376295534`}, {0.08036620662609187, 0.23006231046571013`}, { 0.07799627052036644, 0.21104886666571498`}, {0.0768979556017048, 0.1974918178785205}, {0.07639071109827053, 0.18808488354511366`}, { 0.07615441420618946, 0.18167393001779214`}, {0.07604209622801274, 0.17735629714042697`}, {0.07598703506844022, 0.17447110856925221`}, { 0.07595895090137143, 0.17255302299481026`}, {0.07594396867631163, 0.17128218041539595`}, {0.07593560409041895, 0.17044204525056525`}, { 0.07593073565185665, 0.16988745582847625`}, {0.0759278016361807, 0.16952171207359018`}, {0.0759259849103766, 0.1692806613797533}, { 0.07592483745649234, 0.16912185788192674`}, {0.07592410254612789, 0.16901726702021833`}, {0.07592362736915027, 0.1689483938641411}, { 0.07592331818102135, 0.16890304614746393`}, {0.07592311616184715, 0.16887319042825677`}, {0.07592298380921039, 0.16885353521357244`}, { 0.07592289694787095, 0.16884059582630265`}, {0.07592283987847237, 0.16883207777538617`}, {0.07592280235631227, 0.16882647038761356`}, { 0.07592277767503038, 0.16882277910829158`}, {0.07592276143558867, 0.16882034919608688`}, {0.07592275074867631, 0.1688187496287754}, { 0.07592274371500331, 0.16881769666526325`}, {0.07592273908541314, 0.1688170035214051}, {0.07592273603806671, 0.16881654723981201`}, { 0.07592273403215047, 0.16881624687973992`}, {0.07592273271173404, 0.16881604915944862`}, {0.07592273184254669, 0.16881591900466136`}, { 0.07592273127038537, 0.16881583332673197`}, {0.07592273089374649, 0.16881577692691183`}, {0.07592273064581445, 0.16881573980020334`}, { 0.07592273048260666, 0.16881571536054485`}, {0.07592273037517079, 0.1688156992724801}, {0.07592273030444822, 0.16881568868207755`}, { 0.07592273025789321, 0.16881568171065958`}, {0.0759227302272471, 0.16881567712153572`}, {0.07592273020707346, 0.16881567410062126`}, { 0.07592273019379366, 0.16881567211202286`}, {0.07592273018505184, 0.16881567080297438`}, {0.07592273017929731, 0.16881566994125793`}, { 0.07592273017550923, 0.16881566937400985`}, {0.07592273017301561, 0.16881566900060355`}, {0.07592273017137416, 0.16881566875479884`}, { 0.0759227301702936, 0.1688156685929912}, {0.0759227301695823, 0.16881566848647697`}, {0.07592273016911406, 0.16881566841636111`}, { 0.07592273016880585, 0.1688156683702055}, {0.07592273016860294, 0.16881566833982226`}, {0.07592273016846939, 0.16881566831982167`}, { 0.07592273016838147, 0.1688156683066558}, {0.07592273016832359, 0.168815668297989}}]}, {}}}, AspectRatio->1, Axes->True, AxesOrigin->{0, 0}, Epilog->LineBox[{{1, 0}, {0, 1}}], PlotRange->{{0, 1}, {0, 1}}, PlotRangeClipping->True, PlotRangePadding->{Automatic, Automatic}]], "Output", CellChangeTimes->{3.410175745330635*^9, 3.410176192042976*^9, 3.4101776019002514`*^9, 3.4101848231939616`*^9, 3.410185324835286*^9}] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["\<\ Rectifying Cascade Trajectory at Damk\[ODoubleDot]hler number Da=11.1\ \>", "Subsubtitle", CellChangeTimes->{ 3.410091630290366*^9, {3.410092930199544*^9, 3.4100929425773425`*^9}, { 3.410093659297936*^9, 3.4100937216576047`*^9}, {3.410093789705453*^9, 3.4100937938814573`*^9}, {3.4100939430860033`*^9, 3.410093947842843*^9}, { 3.410172356417618*^9, 3.410172359772442*^9}, {3.410174866637136*^9, 3.410174869040592*^9}}, Background->RGBColor[0.996109, 0.996109, 0.640635]], Cell[BoxData[ RowBox[{ RowBox[{"t", "=", "0"}], ";", RowBox[{ RowBox[{"y", "[", RowBox[{"1", ",", RowBox[{"-", "1"}]}], "]"}], "=", "0.4"}], ";", RowBox[{ RowBox[{"y", "[", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "]"}], "=", "0.4"}], ";", RowBox[{ RowBox[{"y", "[", RowBox[{"3", ",", RowBox[{"-", "1"}]}], "]"}], "=", "0.2"}], ";", RowBox[{ RowBox[{"x", "[", RowBox[{"1", ",", RowBox[{"-", "1"}]}], "]"}], "=", "0.4"}], ";", RowBox[{ RowBox[{"x", "[", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "]"}], "=", "0.4"}], ";", RowBox[{ RowBox[{"x", "[", RowBox[{"3", ",", RowBox[{"-", "1"}]}], "]"}], "=", "0.2"}], ";"}]], "Input", CellChangeTimes->{{3.410091642257574*^9, 3.4100917039863358`*^9}, { 3.4100920383070655`*^9, 3.410092045727736*^9}, {3.410092192739128*^9, 3.4100921996590786`*^9}, {3.4100926628350925`*^9, 3.4100926712471886`*^9}, { 3.4101724027141886`*^9, 3.4101724283210096`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"\[Phi]", "=", "0.5"}], ";", RowBox[{"Da", "=", "11.1"}], ";", " "}]], "Input", CellChangeTimes->{{3.41009177143332*^9, 3.410091803018738*^9}, { 3.410172437514229*^9, 3.410172446336915*^9}, {3.410172602381296*^9, 3.410172613166805*^9}, 3.4101740491115923`*^9, {3.4101748724855456`*^9, 3.410174874698728*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"While", "[", RowBox[{ RowBox[{"t", "<", "60"}], ",", "\[IndentingNewLine]", RowBox[{"{", "\[IndentingNewLine]", RowBox[{ RowBox[{"sol", "=", RowBox[{"Chop", "[", RowBox[{"FindRoot", "[", RowBox[{ RowBox[{"{", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"y", "[", RowBox[{"1", ",", RowBox[{"t", "-", "1"}]}], "]"}], "\[Equal]", " ", RowBox[{ RowBox[{"\[Phi]", " ", FractionBox[ RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "3"], RowBox[{ RowBox[{"y", "[", RowBox[{"i", ",", RowBox[{"t", "-", "1"}]}], "]"}], " ", RowBox[{"M", "[", "i", "]"}]}]}], RowBox[{ RowBox[{"y1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"y2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"y3", " ", RowBox[{"M", "[", "3", "]"}]}]}]], " ", "y1"}], "+", RowBox[{ RowBox[{"(", RowBox[{"1", "-", RowBox[{"\[Phi]", " ", FractionBox[ RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "3"], RowBox[{ RowBox[{"y", "[", RowBox[{"i", ",", RowBox[{"t", "-", "1"}]}], "]"}], " ", RowBox[{"M", "[", "i", "]"}]}]}], RowBox[{ RowBox[{"y1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"y2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"y3", " ", RowBox[{"M", "[", "3", "]"}]}]}]]}]}], ")"}], " ", "x1"}], "-", RowBox[{ FractionBox[ RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "3"], RowBox[{ RowBox[{"y", "[", RowBox[{"i", ",", RowBox[{"t", "-", "1"}]}], "]"}], " ", RowBox[{"M", "[", "i", "]"}]}]}], RowBox[{ RowBox[{"x1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"x2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"x3", " ", RowBox[{"M", "[", "3", "]"}]}]}]], "Da", " ", RowBox[{"(", RowBox[{ RowBox[{"x1", " ", "x2"}], "-", RowBox[{"x3", "/", "Keq"}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "x1"}], ")"}]}]}]}], ",", RowBox[{ RowBox[{"y", "[", RowBox[{"2", ",", RowBox[{"t", "-", "1"}]}], "]"}], "\[Equal]", " ", RowBox[{ RowBox[{"\[Phi]", " ", FractionBox[ RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "3"], RowBox[{ RowBox[{"y", "[", RowBox[{"i", ",", RowBox[{"t", "-", "1"}]}], "]"}], " ", RowBox[{"M", "[", "i", "]"}]}]}], RowBox[{ RowBox[{"y1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"y2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"y3", " ", RowBox[{"M", "[", "3", "]"}]}]}]], " ", "y2"}], "+", RowBox[{ RowBox[{"(", RowBox[{"1", "-", RowBox[{"\[Phi]", FractionBox[ RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "3"], RowBox[{ RowBox[{"y", "[", RowBox[{"i", ",", RowBox[{"t", "-", "1"}]}], "]"}], " ", RowBox[{"M", "[", "i", "]"}]}]}], RowBox[{ RowBox[{"y1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"y2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"y3", " ", RowBox[{"M", "[", "3", "]"}]}]}]]}]}], ")"}], " ", "x2"}], "-", RowBox[{ FractionBox[ RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "3"], RowBox[{ RowBox[{"y", "[", RowBox[{"i", ",", RowBox[{"t", "-", "1"}]}], "]"}], " ", RowBox[{"M", "[", "i", "]"}]}]}], RowBox[{ RowBox[{"x1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"x2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"x3", " ", RowBox[{"M", "[", "3", "]"}]}]}]], "Da", " ", RowBox[{"(", RowBox[{ RowBox[{"x1", " ", "x2"}], "-", RowBox[{"x3", "/", "Keq"}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "x2"}], ")"}]}]}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "[", RowBox[{"3", ",", RowBox[{"t", "-", "1"}]}], "]"}], "\[Equal]", " ", RowBox[{ RowBox[{"\[Phi]", " ", FractionBox[ RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "3"], RowBox[{ RowBox[{"y", "[", RowBox[{"i", ",", RowBox[{"t", "-", "1"}]}], "]"}], " ", RowBox[{"M", "[", "i", "]"}]}]}], RowBox[{ RowBox[{"y1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"y2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"y3", " ", RowBox[{"M", "[", "3", "]"}]}]}]], "y3"}], "+", RowBox[{ RowBox[{"(", RowBox[{"1", "-", RowBox[{"\[Phi]", FractionBox[ RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "3"], RowBox[{ RowBox[{"y", "[", RowBox[{"i", ",", RowBox[{"t", "-", "1"}]}], "]"}], " ", RowBox[{"M", "[", "i", "]"}]}]}], RowBox[{ RowBox[{"y1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"y2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"y3", " ", RowBox[{"M", "[", "3", "]"}]}]}]]}]}], ")"}], " ", "x3"}], "-", RowBox[{ FractionBox[ RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "3"], RowBox[{ RowBox[{"y", "[", RowBox[{"i", ",", RowBox[{"t", "-", "1"}]}], "]"}], " ", RowBox[{"M", "[", "i", "]"}]}]}], RowBox[{ RowBox[{"x1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"x2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"x3", " ", RowBox[{"M", "[", "3", "]"}]}]}]], "Da", " ", RowBox[{"(", RowBox[{ RowBox[{"x1", " ", "x2"}], "-", RowBox[{"x3", "/", "Keq"}]}], ")"}], " ", RowBox[{"(", RowBox[{"1", "+", "x3"}], ")"}]}]}]}], ",", "\[IndentingNewLine]", RowBox[{"y1", "==", RowBox[{"x1", "/", RowBox[{"(", RowBox[{"x1", "+", RowBox[{ RowBox[{ SubscriptBox["\[Alpha]", "BC"], "/", SubscriptBox["\[Alpha]", "AC"]}], " ", "x2"}], "+", RowBox[{ RowBox[{"1", "/", SubscriptBox["\[Alpha]", "AC"]}], " ", "x3"}]}], ")"}]}]}], ",", "\[IndentingNewLine]", RowBox[{"y2", "==", RowBox[{"x2", "/", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["\[Alpha]", "AC"], "/", SubscriptBox["\[Alpha]", "BC"]}], "x1"}], "+", "x2", "+", RowBox[{ RowBox[{"1", "/", SubscriptBox["\[Alpha]", "BC"]}], " ", "x3"}]}], ")"}]}]}], ",", "\[IndentingNewLine]", RowBox[{"y3", "==", RowBox[{"x3", "/", RowBox[{"(", RowBox[{ RowBox[{ SubscriptBox["\[Alpha]", "AC"], " ", "x1"}], "+", RowBox[{ SubscriptBox["\[Alpha]", "BC"], " ", "x2"}], "+", "x3"}], ")"}]}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"x3", ",", RowBox[{"x", "[", RowBox[{"3", ",", RowBox[{"t", "-", "1"}]}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"x1", ",", RowBox[{"x", "[", RowBox[{"1", ",", RowBox[{"t", "-", "1"}]}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"x2", ",", RowBox[{"x", "[", RowBox[{"2", ",", RowBox[{"t", "-", "1"}]}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"y1", ",", RowBox[{"y", "[", RowBox[{"1", ",", RowBox[{"t", "-", "1"}]}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"y2", ",", RowBox[{"y", "[", RowBox[{"2", ",", RowBox[{"t", "-", "1"}]}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"y3", ",", RowBox[{"y", "[", RowBox[{"3", ",", RowBox[{"t", "-", "1"}]}], "]"}]}], "}"}], ",", RowBox[{"MaxIterations", "\[Rule]", "100000"}]}], "]"}], "]"}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"x", "[", RowBox[{"3", ",", "t"}], "]"}], "=", RowBox[{"x3", "/.", "sol"}]}], ",", RowBox[{ RowBox[{"x", "[", RowBox[{"2", ",", "t"}], "]"}], "=", RowBox[{"x2", "/.", "sol"}]}], ",", RowBox[{ RowBox[{"x", "[", RowBox[{"1", ",", "t"}], "]"}], "=", RowBox[{"x1", "/.", "sol"}]}], ",", RowBox[{ RowBox[{"y", "[", RowBox[{"3", ",", "t"}], "]"}], "=", RowBox[{"y3", "/.", "sol"}]}], ",", RowBox[{ RowBox[{"y", "[", RowBox[{"2", ",", "t"}], "]"}], "=", RowBox[{"y2", "/.", "sol"}]}], ",", RowBox[{ RowBox[{"y", "[", RowBox[{"1", ",", "t"}], "]"}], "=", RowBox[{"y1", "/.", "sol"}]}], ",", "\[IndentingNewLine]", RowBox[{"t", "++"}], ",", RowBox[{"Print", "[", "t", "]"}]}], "}"}]}], "]"}], "//", "Timing"}], "\[IndentingNewLine]"}]], "Input", CellChangeTimes->{ 3.4100917219922276`*^9, {3.4100918343938527`*^9, 3.4100919867729635`*^9}, { 3.410092098133091*^9, 3.410092190515931*^9}, {3.4100922333074627`*^9, 3.410092299903222*^9}, {3.410092360229968*^9, 3.4100924322435184`*^9}, { 3.4100924634684176`*^9, 3.410092518737891*^9}, {3.4100925494921136`*^9, 3.410092579335026*^9}, {3.410092618471301*^9, 3.4100926208847713`*^9}, { 3.410092705075832*^9, 3.4100927353894205`*^9}, {3.410092807382942*^9, 3.4100928089051313`*^9}, {3.4100929098803263`*^9, 3.4100929207659793`*^9}, {3.410092969736395*^9, 3.4100929716491456`*^9}, { 3.4100930091731024`*^9, 3.4100930574425106`*^9}, {3.4100930906102033`*^9, 3.4100931742404575`*^9}, {3.41009326282784*^9, 3.4100932847794046`*^9}, { 3.410093450117149*^9, 3.4100934508081427`*^9}, {3.4101726245732064`*^9, 3.4101727053894143`*^9}, {3.410172777883656*^9, 3.4101727890196686`*^9}, { 3.4101728699760785`*^9, 3.410172894451272*^9}, {3.410173078485901*^9, 3.4101731062658463`*^9}, {3.410174053517928*^9, 3.410174064153221*^9}, { 3.4101758345489244`*^9, 3.4101758386848717`*^9}, {3.410176028988515*^9, 3.4101760399142256`*^9}, 3.410183978689624*^9, {3.410184132751154*^9, 3.4101843189689217`*^9}, {3.4101845461355715`*^9, 3.410184592091653*^9}, { 3.4101846495242367`*^9, 3.41018465171739*^9}, {3.4101847762664833`*^9, 3.410184782265109*^9}}], Cell[CellGroupData[{ Cell[BoxData["1"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.41018532903132*^9}], Cell[BoxData["2"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.410185329101421*^9}], Cell[BoxData["3"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.4101853291715217`*^9}], Cell[BoxData["4"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.410185329221594*^9}], Cell[BoxData["5"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.4101853292916946`*^9}], Cell[BoxData["6"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.4101853293417664`*^9}], Cell[BoxData["7"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.410185329391838*^9}], Cell[BoxData["8"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.410185329461939*^9}], Cell[BoxData["9"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.410185329501997*^9}], Cell[BoxData["10"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.4101853295620832`*^9}], Cell[BoxData["11"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.4101853296221695`*^9}], Cell[BoxData["12"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.410185329662227*^9}], Cell[BoxData["13"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.410185329722314*^9}], Cell[BoxData["14"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.4101853297924147`*^9}], Cell[BoxData["15"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.410185329832472*^9}], Cell[BoxData["16"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.410185329882544*^9}], Cell[BoxData["17"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.410185329932616*^9}], Cell[BoxData["18"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.410185330002717*^9}], Cell[BoxData["19"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.4101853300527887`*^9}], Cell[BoxData["20"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.4101853301028605`*^9}], Cell[BoxData["21"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.4101853301729617`*^9}], Cell[BoxData["22"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.4101853302130194`*^9}], Cell[BoxData["23"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.4101853302731056`*^9}], Cell[BoxData["24"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.410185330333192*^9}], Cell[BoxData["25"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.410185330383264*^9}], Cell[BoxData["26"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.4101853304433503`*^9}], Cell[BoxData["27"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.410185330503437*^9}], Cell[BoxData["28"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.4101853305535088`*^9}], Cell[BoxData["29"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.410185330613595*^9}], Cell[BoxData["30"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.4101853306736813`*^9}], Cell[BoxData["31"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.410185330723754*^9}], Cell[BoxData["32"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.41018533078384*^9}], Cell[BoxData["33"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.410185330853941*^9}], Cell[BoxData["34"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.410185330893998*^9}], Cell[BoxData["35"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.410185330954085*^9}], Cell[BoxData["36"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.410185331014171*^9}], Cell[BoxData["37"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.410185331074258*^9}], Cell[BoxData["38"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.4101853311243296`*^9}], Cell[BoxData["39"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.410185331184416*^9}], Cell[BoxData["40"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.4101853312545166`*^9}], Cell[BoxData["41"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.4101853312945747`*^9}], Cell[BoxData["42"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.410185331354661*^9}], Cell[BoxData["43"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.410185331424762*^9}], Cell[BoxData["44"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.410185331464819*^9}], Cell[BoxData["45"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.4101853315249057`*^9}], Cell[BoxData["46"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.410185331584992*^9}], Cell[BoxData["47"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.410185331635064*^9}], Cell[BoxData["48"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.4101853316951504`*^9}], Cell[BoxData["49"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.4101853317552366`*^9}], Cell[BoxData["50"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.410185331805309*^9}], Cell[BoxData["51"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.4101853318653955`*^9}], Cell[BoxData["52"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.410185331925482*^9}], Cell[BoxData["53"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.4101853319755535`*^9}], Cell[BoxData["54"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.41018533203564*^9}], Cell[BoxData["55"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.4101853321057405`*^9}], Cell[BoxData["56"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.4101853321558127`*^9}], Cell[BoxData["57"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.410185332205885*^9}], Cell[BoxData["58"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.4101853322860003`*^9}], Cell[BoxData["59"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.4101853323260574`*^9}], Cell[BoxData["60"], "Print", CellChangeTimes->{3.41017592801332*^9, 3.4101760541747313`*^9, 3.41017757538212*^9, 3.410184336574237*^9, 3.410184658427038*^9, 3.4101847965156*^9, 3.4101848723947086`*^9, 3.4101853323761296`*^9}] }, Open ]], Cell[BoxData[ RowBox[{"{", RowBox[{"0.9010000000000001`", ",", "Null"}], "}"}]], "Output", CellChangeTimes->{{3.4101731015390496`*^9, 3.4101731109025135`*^9}, 3.4101740866555777`*^9, 3.410174233276408*^9, 3.4101747868423967`*^9, 3.4101755031123424`*^9, {3.410175916616933*^9, 3.4101759312980433`*^9}, 3.410176057319253*^9, 3.4101775787469587`*^9, 3.410184340029205*^9, 3.410184661972136*^9, 3.410184799990597*^9, 3.4101848758697057`*^9, 3.410185332426202*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"plt1", "=", RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"y", "[", RowBox[{"1", ",", "t"}], "]"}], ",", RowBox[{"y", "[", RowBox[{"2", ",", "t"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"t", ",", RowBox[{"-", "1"}], ",", "59"}], "}"}]}], "]"}], ",", RowBox[{"PlotRange", "\[Rule]", " ", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "1"}], "}"}]}], "}"}]}], ",", RowBox[{"AspectRatio", "\[Rule]", " ", "1"}], ",", RowBox[{"Epilog", "\[Rule]", " ", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "1"}], "}"}]}], "}"}], "]"}]}]}], "]"}]}]], "Input",\ CellChangeTimes->{{3.4100935088315763`*^9, 3.410093509402397*^9}, { 3.41017311487823*^9, 3.4101731172015715`*^9}, {3.4101742461549263`*^9, 3.410174247476827*^9}, {3.41017469207613*^9, 3.410174725203765*^9}}], Cell[BoxData[ GraphicsBox[{{}, {Hue[0.67, 0.6, 0.6], PointBox[{{0.4, 0.4}, {0.503594279051528, 0.40982616451610426`}, { 0.5612885402024901, 0.3574724874974058}, {0.6257823541156673, 0.29653080369536766`}, {0.6973562737431916, 0.2316487711972538}, { 0.7717784276704847, 0.16768972755127018`}, {0.8418828600892818, 0.11108411710985212`}, {0.8999081304977128, 0.06730337278742694}, { 0.9415711279522833, 0.037868818643307725`}, {0.967900917646925, 0.020262368947436226`}, {0.983041004379155, 0.010529516460701263`}, { 0.9912401241724712, 0.005387624039712777}, {0.9955303439501497, 0.0027349953373423407`}, {0.997733976091146, 0.0013828955634971496`}, { 0.9988549515333836, 0.0006978383188074628}, {0.9994223625980388, 0.00035179063471812527`}, {0.9997088490061523, 0.0001772533424359455}, { 0.9998533117895213, 0.00008928823189032727}, {0.999926111247593, 0.00004497161133667719}, {0.9999627853290821, 0.00002264929544413429}, { 0.9999812575562045, 0.000011406616930786974`}, {0.9999905609941331, 5.7444958794818506`*^-6}, {0.9999952464260604, 2.8929665684137027`*^-6}, {0.9999976060713162, 1.4569112884354696`*^-6}, {0.9999987944072404, 7.33705704405612*^-7}, { 0.9999993928593909, 3.6949641295177665`*^-7}, {0.9999996942422003, 1.8607941320844104`*^-7}, {0.9999998460195405, 9.371008101727996*^-8}, { 0.9999999224550395, 4.7192636490002857`*^-8}, {0.9999999609481603, 2.3766330838558073`*^-8}, {0.9999999803333962, 1.1968782971801074`*^-8}, {0.9999999900958493, 6.027508921973992*^-9}, { 0.999999995012245, 3.0354685056733018`*^-9}, {0.9999999974881543, 1.5286694976795835`*^-9}, {0.9999999987350283, 7.698417648017802*^-10}, { 0.9999999993629571, 3.876942292046757*^-10}, {0.9999999996791836, 1.9524377891995226`*^-10}, {0.9999999998384361, 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}}]}, {}}, AspectRatio->1, Axes->True, AxesOrigin->{0, 0}, Epilog->LineBox[{{1, 0}, {0, 1}}], PlotRange->{{0, 1}, {0, 1}}, PlotRangeClipping->True, PlotRangePadding->{Automatic, Automatic}]], "Output", CellChangeTimes->{{3.4100934996784143`*^9, 3.4100935105340242`*^9}, 3.410173119524912*^9, 3.4101740879774785`*^9, {3.4101742345582514`*^9, 3.4101742481077347`*^9}, 3.410174726355421*^9, 3.410174886045043*^9, 3.410175538913822*^9, 3.4101761258778353`*^9, 3.4101776180434637`*^9, 3.410184877782456*^9, 3.410185333688016*^9}] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["\<\ Stripping Cascade Trajectory at Damk\[ODoubleDot]hler number Da=11.1\ \>", "Subsubtitle", CellChangeTimes->{ 3.410091630290366*^9, {3.410092930199544*^9, 3.4100929425773425`*^9}, { 3.410093659297936*^9, 3.4100937216576047`*^9}, {3.4100937524218416`*^9, 3.4100937804120893`*^9}, {3.410093953370792*^9, 3.410093953751339*^9}, { 3.410173188454027*^9, 3.4101731977473907`*^9}, {3.410174895007931*^9, 3.4101748974714737`*^9}}, Background->RGBColor[0.996109, 0.996109, 0.640635]], Cell[BoxData[ RowBox[{ RowBox[{"t", "=", "0"}], ";", RowBox[{ RowBox[{"y", "[", RowBox[{"1", ",", RowBox[{"-", "1"}]}], "]"}], "=", "0.4"}], ";", RowBox[{ RowBox[{"y", "[", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "]"}], "=", "0.4"}], ";", RowBox[{ RowBox[{"y", "[", RowBox[{"3", ",", RowBox[{"-", "1"}]}], "]"}], "=", "0.2"}], ";", RowBox[{ RowBox[{"x", "[", RowBox[{"1", ",", RowBox[{"-", "1"}]}], "]"}], "=", "0.4"}], ";", RowBox[{ RowBox[{"x", "[", RowBox[{"2", ",", RowBox[{"-", "1"}]}], "]"}], "=", "0.4"}], ";", RowBox[{ RowBox[{"x", "[", RowBox[{"3", ",", RowBox[{"-", "1"}]}], "]"}], "=", "0.2"}], ";"}]], "Input", CellChangeTimes->{{3.410091642257574*^9, 3.4100917039863358`*^9}, { 3.4100920383070655`*^9, 3.410092045727736*^9}, {3.410092192739128*^9, 3.4100921996590786`*^9}, {3.4100926628350925`*^9, 3.4100926712471886`*^9}, { 3.4101724027141886`*^9, 3.4101724283210096`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"\[Phi]", "=", "0.5"}], ";", RowBox[{"Da", "=", "11.1"}], ";", " "}]], "Input", CellChangeTimes->{{3.41009177143332*^9, 3.410091803018738*^9}, { 3.410172437514229*^9, 3.410172446336915*^9}, {3.410172602381296*^9, 3.410172613166805*^9}, {3.4101739474754467`*^9, 3.4101739505698957`*^9}, { 3.410174126793293*^9, 3.410174129296893*^9}, 3.4101746810302467`*^9, { 3.410174903960805*^9, 3.4101749062140446`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"While", "[", RowBox[{ RowBox[{"t", "<", "200"}], ",", "\[IndentingNewLine]", RowBox[{"{", "\[IndentingNewLine]", RowBox[{ RowBox[{"sol", "=", RowBox[{"Chop", "[", RowBox[{"FindRoot", "[", RowBox[{ RowBox[{"{", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"x", "[", RowBox[{"1", ",", RowBox[{"t", "-", "1"}]}], "]"}], "\[Equal]", " ", RowBox[{ RowBox[{"\[Phi]", " ", FractionBox[ RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "3"], RowBox[{ RowBox[{"x", "[", RowBox[{"i", ",", RowBox[{"t", "-", "1"}]}], "]"}], " ", RowBox[{"M", "[", "i", "]"}]}]}], RowBox[{ RowBox[{"y1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"y2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"y3", " ", RowBox[{"M", "[", "3", "]"}]}]}]], "y1"}], "+", RowBox[{ RowBox[{"(", RowBox[{"1", "-", RowBox[{"\[Phi]", FractionBox[ RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "3"], RowBox[{ RowBox[{"x", "[", RowBox[{"i", ",", RowBox[{"t", "-", "1"}]}], "]"}], " ", RowBox[{"M", "[", "i", "]"}]}]}], RowBox[{ RowBox[{"y1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"y2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"y3", " ", RowBox[{"M", "[", "3", "]"}]}]}]]}]}], ")"}], " ", "x1"}], "-", RowBox[{"Da", " ", FractionBox[ RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "3"], RowBox[{ RowBox[{"x", "[", RowBox[{"i", ",", RowBox[{"t", "-", "1"}]}], "]"}], " ", RowBox[{"M", "[", "i", "]"}]}]}], RowBox[{ RowBox[{"x1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"x2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"x3", " ", RowBox[{"M", "[", "3", "]"}]}]}]], RowBox[{"(", RowBox[{ RowBox[{"x1", " ", "x2"}], "-", RowBox[{"x3", "/", "Keq"}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "x1"}], ")"}]}]}]}], ",", RowBox[{ RowBox[{"x", "[", RowBox[{"2", ",", RowBox[{"t", "-", "1"}]}], "]"}], "\[Equal]", " ", RowBox[{ RowBox[{"\[Phi]", " ", FractionBox[ RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "3"], RowBox[{ RowBox[{"x", "[", RowBox[{"i", ",", RowBox[{"t", "-", "1"}]}], "]"}], " ", RowBox[{"M", "[", "i", "]"}]}]}], RowBox[{ RowBox[{"y1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"y2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"y3", " ", RowBox[{"M", "[", "3", "]"}]}]}]], "y2"}], "+", RowBox[{ RowBox[{"(", RowBox[{"1", "-", RowBox[{"\[Phi]", FractionBox[ RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "3"], RowBox[{ RowBox[{"x", "[", RowBox[{"i", ",", RowBox[{"t", "-", "1"}]}], "]"}], " ", RowBox[{"M", "[", "i", "]"}]}]}], RowBox[{ RowBox[{"y1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"y2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"y3", " ", RowBox[{"M", "[", "3", "]"}]}]}]]}]}], ")"}], " ", "x2"}], "-", RowBox[{"Da", " ", FractionBox[ RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "3"], RowBox[{ RowBox[{"x", "[", RowBox[{"i", ",", RowBox[{"t", "-", "1"}]}], "]"}], " ", RowBox[{"M", "[", "i", "]"}]}]}], RowBox[{ RowBox[{"x1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"x2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"x3", " ", RowBox[{"M", "[", "3", "]"}]}]}]], RowBox[{"(", RowBox[{ RowBox[{"x1", " ", "x2"}], "-", RowBox[{"x3", "/", "Keq"}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "x2"}], ")"}]}]}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"x1", "+", "x2", "+", "x3"}], "==", "1"}], ",", "\[IndentingNewLine]", RowBox[{"y1", "==", RowBox[{"x1", "/", RowBox[{"(", RowBox[{"x1", "+", RowBox[{ RowBox[{ SubscriptBox["\[Alpha]", "BC"], "/", SubscriptBox["\[Alpha]", "AC"]}], " ", "x2"}], "+", RowBox[{ RowBox[{"1", "/", SubscriptBox["\[Alpha]", "AC"]}], " ", "x3"}]}], ")"}]}]}], ",", "\[IndentingNewLine]", RowBox[{"y2", "==", RowBox[{"x2", "/", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["\[Alpha]", "AC"], "/", SubscriptBox["\[Alpha]", "BC"]}], "x1"}], "+", "x2", "+", RowBox[{ RowBox[{"1", "/", SubscriptBox["\[Alpha]", "BC"]}], " ", "x3"}]}], ")"}]}]}], ",", "\[IndentingNewLine]", RowBox[{"y3", "==", RowBox[{"x3", "/", RowBox[{"(", RowBox[{ RowBox[{ SubscriptBox["\[Alpha]", "AC"], " ", "x1"}], "+", RowBox[{ SubscriptBox["\[Alpha]", "BC"], " ", "x2"}], "+", "x3"}], ")"}]}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"x3", ",", RowBox[{"x", "[", RowBox[{"3", ",", RowBox[{"t", "-", "1"}]}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"x1", ",", RowBox[{"x", "[", RowBox[{"1", ",", RowBox[{"t", "-", "1"}]}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"x2", ",", RowBox[{"x", "[", RowBox[{"2", ",", RowBox[{"t", "-", "1"}]}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"y1", ",", RowBox[{"y", "[", RowBox[{"1", ",", RowBox[{"t", "-", "1"}]}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"y2", ",", RowBox[{"y", "[", RowBox[{"2", ",", RowBox[{"t", "-", "1"}]}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"y3", ",", RowBox[{"y", "[", RowBox[{"3", ",", RowBox[{"t", "-", "1"}]}], "]"}]}], "}"}], ",", RowBox[{"MaxIterations", "\[Rule]", "1000000"}]}], "]"}], "]"}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"x", "[", RowBox[{"3", ",", "t"}], "]"}], "=", RowBox[{"x3", "/.", "sol"}]}], ",", RowBox[{ RowBox[{"x", "[", RowBox[{"2", ",", "t"}], "]"}], "=", RowBox[{"x2", "/.", "sol"}]}], ",", RowBox[{ RowBox[{"x", "[", RowBox[{"1", ",", "t"}], "]"}], "=", RowBox[{"x1", "/.", "sol"}]}], ",", RowBox[{ RowBox[{"y", "[", RowBox[{"3", ",", "t"}], "]"}], "=", RowBox[{"y3", "/.", "sol"}]}], ",", RowBox[{ RowBox[{"y", "[", RowBox[{"2", ",", "t"}], "]"}], "=", RowBox[{"y2", "/.", "sol"}]}], ",", RowBox[{ RowBox[{"y", "[", RowBox[{"1", ",", "t"}], "]"}], "=", RowBox[{"y1", "/.", "sol"}]}], ",", "\[IndentingNewLine]", RowBox[{"t", "++"}], ",", RowBox[{"Print", "[", "t", "]"}]}], "}"}]}], "]"}], "//", "Timing"}], "\[IndentingNewLine]"}]], "Input", CellChangeTimes->{ 3.4100917219922276`*^9, {3.4100918343938527`*^9, 3.4100919867729635`*^9}, { 3.410092098133091*^9, 3.410092190515931*^9}, {3.4100922333074627`*^9, 3.410092299903222*^9}, {3.410092360229968*^9, 3.4100924322435184`*^9}, { 3.4100924634684176`*^9, 3.410092518737891*^9}, {3.4100925494921136`*^9, 3.410092579335026*^9}, {3.410092618471301*^9, 3.4100926208847713`*^9}, { 3.410092705075832*^9, 3.4100927353894205`*^9}, {3.410092807382942*^9, 3.4100928089051313`*^9}, {3.4100929098803263`*^9, 3.4100929207659793`*^9}, {3.410092969736395*^9, 3.4100929716491456`*^9}, { 3.4100930091731024`*^9, 3.4100930574425106`*^9}, {3.4100930906102033`*^9, 3.4100931742404575`*^9}, {3.41009326282784*^9, 3.4100932847794046`*^9}, { 3.410093450117149*^9, 3.4100934508081427`*^9}, {3.4101731565181055`*^9, 3.410173158531*^9}, {3.4101732281310797`*^9, 3.410173245936683*^9}, 3.410173288197451*^9, {3.4101734446524224`*^9, 3.410173457360696*^9}, { 3.410173575921178*^9, 3.4101735769426465`*^9}, {3.410173748779736*^9, 3.410173756420723*^9}, {3.410174117439843*^9, 3.4101741343441505`*^9}, { 3.410174523293432*^9, 3.4101745339687824`*^9}, {3.4101758444832096`*^9, 3.410175846015413*^9}, {3.410176071319384*^9, 3.4101760755755043`*^9}, { 3.410184352777536*^9, 3.4101844114919634`*^9}, {3.4101846812999277`*^9, 3.4101846977335587`*^9}, {3.410184737981432*^9, 3.4101847431388483`*^9}, 3.4101849065438128`*^9, {3.4101853399069586`*^9, 3.4101853408883696`*^9}}],\ Cell[CellGroupData[{ Cell[BoxData["1"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185341769637*^9}], Cell[BoxData["2"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185341839738*^9}], Cell[BoxData["3"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853419098387`*^9}], Cell[BoxData["4"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853419599104`*^9}], Cell[BoxData["5"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185342030011*^9}], Cell[BoxData["6"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185342080083*^9}], Cell[BoxData["7"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185342130155*^9}], Cell[BoxData["8"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185342190242*^9}], Cell[BoxData["9"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853422403135`*^9}], Cell[BoxData["10"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853423004*^9}], Cell[BoxData["11"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853423705006`*^9}], Cell[BoxData["12"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185342400544*^9}], Cell[BoxData["13"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853424606304`*^9}], Cell[BoxData["14"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853425207167`*^9}], Cell[BoxData["15"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185342580803*^9}], Cell[BoxData["16"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185342620861*^9}], Cell[BoxData["17"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853426809473`*^9}], Cell[BoxData["18"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853427410336`*^9}], Cell[BoxData["19"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853428111343`*^9}], Cell[BoxData["20"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185342851192*^9}], Cell[BoxData["21"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185342901264*^9}], Cell[BoxData["22"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853429613504`*^9}], Cell[BoxData["23"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185343031451*^9}], Cell[BoxData["24"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185343081523*^9}], Cell[BoxData["25"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185343131595*^9}], Cell[BoxData["26"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853432117105`*^9}], Cell[BoxData["27"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185343251768*^9}], Cell[BoxData["28"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853433118544`*^9}], Cell[BoxData["29"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185343361926*^9}], Cell[BoxData["30"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185343432027*^9}], Cell[BoxData["31"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853434820995`*^9}], Cell[BoxData["32"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185343542186*^9}], Cell[BoxData["33"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853435922575`*^9}], Cell[BoxData["34"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853436523438`*^9}], Cell[BoxData["35"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853437124305`*^9}], Cell[BoxData["36"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853437825313`*^9}], Cell[BoxData["37"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853438326035`*^9}], Cell[BoxData["38"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185343882675*^9}], Cell[BoxData["39"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853439627905`*^9}], Cell[BoxData["40"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853440028477`*^9}], Cell[BoxData["41"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853440629344`*^9}], Cell[BoxData["42"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853441130066`*^9}], Cell[BoxData["43"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853441831074`*^9}], Cell[BoxData["44"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185344233179*^9}], Cell[BoxData["45"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853442932653`*^9}], Cell[BoxData["46"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853443433375`*^9}], Cell[BoxData["47"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853444034243`*^9}], Cell[BoxData["48"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853444635105`*^9}], Cell[BoxData["49"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853445336113`*^9}], Cell[BoxData["50"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185344583683*^9}], Cell[BoxData["51"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185344633755*^9}], Cell[BoxData["52"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185344703856*^9}], Cell[BoxData["53"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185344753928*^9}], Cell[BoxData["54"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185344804*^9}], Cell[BoxData["55"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185344864086*^9}], Cell[BoxData["56"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185344934187*^9}], Cell[BoxData["57"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185344984259*^9}], Cell[BoxData["58"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185345044346*^9}], Cell[BoxData["59"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853450944176`*^9}], Cell[BoxData["60"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185345154504*^9}], Cell[BoxData["61"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.41018534521459*^9}], Cell[BoxData["62"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185345274677*^9}], Cell[BoxData["63"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853453347635`*^9}], Cell[BoxData["64"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853453848352`*^9}], Cell[BoxData["65"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185345454936*^9}], Cell[BoxData["66"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853455050077`*^9}], Cell[BoxData["67"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.41018534555508*^9}], Cell[BoxData["68"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853456151667`*^9}], Cell[BoxData["69"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853456852674`*^9}], Cell[BoxData["70"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185345735339*^9}], Cell[BoxData["71"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185345785411*^9}], Cell[BoxData["72"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853458454976`*^9}], Cell[BoxData["73"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853459155984`*^9}], Cell[BoxData["74"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853459656706`*^9}], Cell[BoxData["75"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185346025757*^9}], Cell[BoxData["76"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185346085843*^9}], Cell[BoxData["77"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853461359153`*^9}], Cell[BoxData["78"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853461960015`*^9}], Cell[BoxData["79"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853462661023`*^9}], Cell[BoxData["80"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853463161745`*^9}], Cell[BoxData["81"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853463762608`*^9}], Cell[BoxData["82"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853464263325`*^9}], Cell[BoxData["83"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185346486419*^9}], Cell[BoxData["84"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185346546506*^9}], Cell[BoxData["85"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853466166067`*^9}], Cell[BoxData["86"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853466666784`*^9}], Cell[BoxData["87"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.41018534671675*^9}], Cell[BoxData["88"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185346776837*^9}], Cell[BoxData["89"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853468469377`*^9}], Cell[BoxData["90"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.41018534689701*^9}], Cell[BoxData["91"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185346957096*^9}], Cell[BoxData["92"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185347027197*^9}], Cell[BoxData["93"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185347067254*^9}], Cell[BoxData["94"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185347127341*^9}], Cell[BoxData["95"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853471874275`*^9}], Cell[BoxData["96"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185347247514*^9}], Cell[BoxData["97"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853473076*^9}], Cell[BoxData["98"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853473576717`*^9}], Cell[BoxData["99"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185347437787*^9}], Cell[BoxData["100"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853474778447`*^9}], Cell[BoxData["101"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853475379314`*^9}], Cell[BoxData["102"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853475980177`*^9}], Cell[BoxData["103"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853476681185`*^9}], Cell[BoxData["104"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853477081757`*^9}], Cell[BoxData["105"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853477682624`*^9}], Cell[BoxData["106"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185347828349*^9}], Cell[BoxData["107"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.41018534789845*^9}], Cell[BoxData["108"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185347938507*^9}], Cell[BoxData["109"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853479985933`*^9}], Cell[BoxData["110"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.41018534805868*^9}], Cell[BoxData["111"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185348128781*^9}], Cell[BoxData["112"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853481688385`*^9}], Cell[BoxData["113"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853482289248`*^9}], Cell[BoxData["114"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185348289011*^9}], Cell[BoxData["115"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185348359112*^9}], Cell[BoxData["116"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.41018534839917*^9}], Cell[BoxData["117"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185348459256*^9}], Cell[BoxData["118"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853485193424`*^9}], Cell[BoxData["119"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853485794287`*^9}], Cell[BoxData["120"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185348639515*^9}], Cell[BoxData["121"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853486996017`*^9}], Cell[BoxData["122"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853487697024`*^9}], Cell[BoxData["123"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.41018534880976*^9}], Cell[BoxData["124"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853488698463`*^9}], Cell[BoxData["125"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853489299326`*^9}], Cell[BoxData["126"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853490000334`*^9}], Cell[BoxData["127"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853490400915`*^9}], Cell[BoxData["128"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185349100178*^9}], Cell[BoxData["129"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185349160264*^9}], Cell[BoxData["130"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185349230365*^9}], Cell[BoxData["131"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853492804365`*^9}], Cell[BoxData["132"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853493405232`*^9}], Cell[BoxData["133"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853493905954`*^9}], Cell[BoxData["134"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185349460696*^9}], Cell[BoxData["135"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185349510768*^9}], Cell[BoxData["136"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185349570854*^9}], Cell[BoxData["137"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185349640955*^9}], Cell[BoxData["138"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185349691027*^9}], Cell[BoxData["139"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185349751114*^9}], Cell[BoxData["140"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853498011856`*^9}], Cell[BoxData["141"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185349861272*^9}], Cell[BoxData["142"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185349921358*^9}], Cell[BoxData["143"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185349991459*^9}], Cell[BoxData["144"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185350041531*^9}], Cell[BoxData["145"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185350101618*^9}], Cell[BoxData["146"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185350161704*^9}], Cell[BoxData["147"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853502217903`*^9}], Cell[BoxData["148"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853502718625`*^9}], Cell[BoxData["149"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853503519773`*^9}], Cell[BoxData["150"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853503920355`*^9}], Cell[BoxData["151"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853504521217`*^9}], Cell[BoxData["152"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185350512208*^9}], Cell[BoxData["153"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185350582309*^9}], Cell[BoxData["154"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853506223664`*^9}], Cell[BoxData["155"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185350682453*^9}], Cell[BoxData["156"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853507425394`*^9}], Cell[BoxData["157"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.41018535081264*^9}], Cell[BoxData["158"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185350862712*^9}], Cell[BoxData["159"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185350922798*^9}], Cell[BoxData["160"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185350982885*^9}], Cell[BoxData["161"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185351042971*^9}], Cell[BoxData["162"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185351103058*^9}], Cell[BoxData["163"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853511731586`*^9}], Cell[BoxData["164"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185351213216*^9}], Cell[BoxData["165"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185351273302*^9}], Cell[BoxData["166"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185351333389*^9}], Cell[BoxData["167"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853514034896`*^9}], Cell[BoxData["168"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185351453562*^9}], Cell[BoxData["169"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185351513648*^9}], Cell[BoxData["170"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853515737343`*^9}], Cell[BoxData["171"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185351643835*^9}], Cell[BoxData["172"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853516838927`*^9}], Cell[BoxData["173"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853517439795`*^9}], Cell[BoxData["174"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853518040657`*^9}], Cell[BoxData["175"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853518741665`*^9}], Cell[BoxData["176"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185351924238*^9}], Cell[BoxData["177"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853519843245`*^9}], Cell[BoxData["178"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185352044411*^9}], Cell[BoxData["179"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185352104498*^9}], Cell[BoxData["180"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185352164584*^9}], Cell[BoxData["181"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185352234685*^9}], Cell[BoxData["182"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185352274742*^9}], Cell[BoxData["183"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185352334829*^9}], Cell[BoxData["184"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185352394915*^9}], Cell[BoxData["185"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185352455002*^9}], Cell[BoxData["186"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185352515088*^9}], Cell[BoxData["187"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853525751743`*^9}], Cell[BoxData["188"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185352645275*^9}], Cell[BoxData["189"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853526953473`*^9}], Cell[BoxData["190"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853527454195`*^9}], Cell[BoxData["191"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853528155203`*^9}], Cell[BoxData["192"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185352885621*^9}], Cell[BoxData["193"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853529256783`*^9}], Cell[BoxData["194"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853529857645`*^9}], Cell[BoxData["195"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185353045851*^9}], Cell[BoxData["196"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185353115952*^9}], Cell[BoxData["197"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185353166024*^9}], Cell[BoxData["198"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853532261105`*^9}], Cell[BoxData["199"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.4101853532861967`*^9}], Cell[BoxData["200"], "Print", CellChangeTimes->{3.410184907575296*^9, 3.410185353346283*^9}] }, Open ]], Cell[BoxData[ RowBox[{"{", RowBox[{"2.534`", ",", "Null"}], "}"}]], "Output", CellChangeTimes->{3.4101849108199615`*^9, 3.410185353386341*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"plt2", "=", RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"x", "[", RowBox[{"1", ",", "t"}], "]"}], ",", RowBox[{"x", "[", RowBox[{"2", ",", "t"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"t", ",", RowBox[{"-", "1"}], ",", "199"}], "}"}]}], "]"}], ",", RowBox[{"PlotRange", "\[Rule]", " ", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "1"}], "}"}]}], "}"}]}], ",", RowBox[{"AspectRatio", "\[Rule]", " ", "1"}], ",", RowBox[{"PlotStyle", "\[Rule]", "Red"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.4100935088315763`*^9, 3.410093509402397*^9}, { 3.410173302958677*^9, 3.410173305141816*^9}, {3.410173992710491*^9, 3.410173995113947*^9}, {3.4101741421053104`*^9, 3.410174143096736*^9}, { 3.4101742191460896`*^9, 3.4101742206482496`*^9}, {3.410174738482859*^9, 3.4101747403455377`*^9}, {3.410175567264589*^9, 3.4101755761974335`*^9}, { 3.4101756525872765`*^9, 3.410175653067968*^9}, {3.410176167457624*^9, 3.4101761698210225`*^9}, {3.410185356641021*^9, 3.410185357201827*^9}}], Cell[BoxData[ GraphicsBox[{{}, {RGBColor[1, 0, 0], PointBox[CompressedData[" 1:eJw10gs01dkeB3B55B1iypQucXKIRNGM54+QwjGNETKG8qiphiKPqByixBRT 5HG851J6SKiuPP7HMxEi71ee509FepoZl9y991l3r7XXXp+1fmv9vr+992av k46+ggICAs1o45OTjtcA/P8sULv5OZ3ugk2LQQYry4OgPvnNrRj1TnBoXVYP HR2BhPmA2fTQdjCmGJL/1RkDrbJKejjlOYxZXXnF1R+HkZVwNutQM+Tbn6Sq JCZAd/+yY0l+EzgcqrYPrpiAfOqKGNfyKdy8PKRfYjcJo47KXfPqjdBzg1mU 2DAJxrfawMqsAZxerTR1a06BvPrfBzQj6iH9Z4P9u9hT4PFBvGD9qzp4rf7Q N6NpCk6XsFpvu9SB76aZTpYID2rm+1tXeLVQXxy/T9mIB1TA2Bani7UwbSzi cPEoD5iXWHdKdtSCpIqr6P0EHii6Gepw39aAupzObOIDHlwZuCXCuV8D/krz KV/aeHDM/cmTirAaOLGx49mlaR7EMlsmx+xqYK1Xy4EXyzxos7ZVEmbUgO8a rygrORombTjX9q6qgbOP3bfvU6VBTedJ0mwjF+TqlPT09GiYU4msmPqDC5Kx DjxLMxrExWMdT3txoWq880dqHw2/qUdbPDHkQjKruaf7JxrW+JS++esbLoRp +cQ/dKfhBqO84/cFCtxck7xjfGhgWa0zmRmkIJzTvenICRosN6rdT6qjQDnJ Qig0gIbENc1fGEUUaAa0JnaEoH7yto3mHAqSm0cWL5+lwfPDneLUOApaqr+7 W8am4VD217CpcAo0GgTSD0fTwNBfXfLGjwKxdtm8vEs0pPr2hQR7UaAqJNAa GUfDbbuKCGdXCq65L8VLXKGhyMlzu+8PFGw1lpzdm0CDfKYc77wNBYK1Mk1W f9AwUF7hEGWO8j0rmRK5TsPNAxHjPxlRIHpEViUpiQZJ+SDPEX0KbvgJxXxM pmFIP9dfTJcC2R8TfTRT0Hz1TaZ1WhQMchSCTVJpcPx1s8Eik4K8yvtD29No eOG9qySPQYHHd4zdouk0MHXin97dTIFRQuVUI3JOW6+UtDLyapMcfw4NfzYO nKpUoiDFwuVvoQwaMsu7JO9voODlQrHWZeT++qWVYUUKrMfu7FhGNvC3PWix noKv3bqtPpk0/Jq3+KH7GwqONURH1yBfYR13u65AQdiA5WO5LBoe5bTXB8lT cMeQSj+IvJQWPBe6lgKXqxmQjryh3EwhTY6C3HurOJ3IzOzVU52yFDDcDdlC 2TRc8Go+rorcNUPf2468ZWKPdKwMBQFZEy3OyJ6vF88tr6FgVXCVXxiy3XPz 4mjkExJfFNOw5yN61iO/T/xsXoacpWPsXC6N5o1vsn+OPJC79NQbOUelMHUM +VrNBK2IPLx4MPATcknZcHqvFAXnBD0YwjmoXsx/KgN5T4Hg3FrkbYHvp48i R1c7Casgd/5lOWqIzOZ+ztDC7hzbKIeccrhr1ABZkgNCc5IUCCf3L5ghuw9c lWpDNt0XKmSDzNt0k3qArGCty3BAXm63vJ+KbDwtEuSEfHvoxLkoZOruhIQb cqq2cI8fcuvMp2UP5GehyX7uyKsYi57eyC9KtYRZyJ56f1gcRZZ55+sByH1b j5ccRw4ZTdi7E1la92qVH/LXjpFkDeTwukj/U8hzYs82KCNHh4V3BiKbrD/6 eB3yg4cx74KQu7nvbWSQhUZz20KQD50zKxNDjmrkBZ9BFm/1GhVE3hbEmQtD flpp+3hZggL3FZbRWeRW8+OMf5DvJc0eOoccyAxU/IJ8Ueml73nkzPNroj8g DzcN2kcgj7LvOb9DnigI+JaNvKBtE/sWWTwzqx37QZWhxGvkC8NSJyORY1jr e2nkDGnlZWx53diJKeTrNu/Co5CZ+7ZoTCL3SLi8wbZaeyR3HPkf+xX7C8hT 7eetx5AF/+Oej22RZrt2FLlau+IDduisueArZLlSk13RyP6hbnIjyAM104HY pt2rTYaRS5XzCrEPSP+LPYRsvVe/H9uwsLtvEPc7dkAwBjlepsQS21shn4n9 r50XagaQU5OK9mI/N3axx2aLzvtii+ie4fUjx5n2sbFLAj3jsY87PUvBhpUy Q+zq7d/exa50sPrYh6xUu64au/HlqjLsWC2RNuzaz0vh2KGNKkPY7PjD+7B9 lrOnsRNKT2/CVmUUf8Tukbq60IusoHJ2CdvIakM39qSkhshFnHcs7RH2cl2v FHb8wRscbEmPfHnsEY/bF7DV8u58iy3iEu6HffgXAWVsUdFyN2xW0xM17Dav jbbYPvUzTGypSF9j7H7NLC1s8UIPHexUQZ4OdlRMjyr2yYpWPex+n2xF7A2O nvrYuWnXZbBnY2/uwlYfSRTFfrx053vsy0OhAthhTkFG2G8fbV/sQXb+SdwE +15ezmdsjb9OmmKzrpbOY5skFJlhF3YceIvt+/wpYO+mz0xjN/dxzbErtaWn sAu8Mi2wfYpFxrHNiz13Y2/94vwK2ypPyhJ7R8/CEHZJUwGx7f7hAezrR7ZZ YU/4r/Rhx7ELiEdUD/Zinx6QscYe1p7sxl4v6k98WTitC5vpV0t8qz7wJXak jeQe7FT5k53YPUX2xC+U4jqw+1MvEbueqX6BPT72hPijqjixlt00cXbAsXbs 8gsyNtiJosNt2H3KO4mZbR7EfrcdiTVm5luxZ1L8+fUPrxLLOsYS/0btIk78 PovYp+r1c/IeWQ+Iu+ILiAeba4lPMY8Rq610EMcV7SRuC3lFfM1RmF/v/YZY 23uwhdzXp8/ESeaPiFNZX4kVPZKJy26s3outBaHE6m+liROCfyEuPaZAfD1n D/EnzQ3EZUd2EJdYKBMfilAhNmhQI376pyxxch2T+GWqIPE7Oy1iV42FZuxC Px1iyy1vibN26hH3+owT7y/cyc9X108s1WNAHKvUSSxb9x3xR+cW4oIQQ+IW lwbi3f8YEbuJcYlzfjAhjnKqIG6NMCUO2fqY2Op3M+JT4aXEgmFAfHZXMXG8 nTmxjuM94kURC+LMutvE8nf5Xhdxi3iP6W5i4XMFxFUU360P/018WM+S+Obm P5v5/4Xvwqpc4o73fMdF5hD/Zm5FzPk1m3juEt+lp7KIf27kW/NGJnH+Et+1 LzOIM7ZZ8/Oq8/3FlW+NBA7xQTbf2qJ8u+bxnZOYTnyay3fDFr6VBvle/SyN uPwD3/8DU01O6Q== "]]}, {}}, AspectRatio->1, Axes->True, AxesOrigin->{0, 0}, PlotRange->{{0, 1}, {0, 1}}, PlotRangeClipping->True, PlotRangePadding->{Automatic, Automatic}]], "Output", CellChangeTimes->{{3.4100934996784143`*^9, 3.4100935105340242`*^9}, 3.4100935905390654`*^9, 3.410173306593904*^9, 3.410173468186262*^9, 3.4101735894906893`*^9, 3.41017377926357*^9, {3.410173986171088*^9, 3.4101739962756176`*^9}, 3.410174144598896*^9, 3.4101745426913247`*^9, 3.410174912653304*^9, 3.4101755776094637`*^9, 3.4101756537789907`*^9, 3.4101761712230387`*^9, 3.4101776473856564`*^9, 3.4101849126425824`*^9, 3.4101853580830946`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"x", "[", RowBox[{"1", ",", "199"}], "]"}], ",", RowBox[{"x", "[", RowBox[{"2", ",", "199"}], "]"}]}], "}"}]], "Input", CellChangeTimes->{{3.4101734813251553`*^9, 3.4101735363943405`*^9}, { 3.4101737905197554`*^9, 3.4101737958574305`*^9}, {3.410174170205717*^9, 3.4101742096324096`*^9}, {3.4101755815551376`*^9, 3.4101755858012433`*^9}, { 3.4101756573441167`*^9, 3.4101756611696177`*^9}, {3.410184926632699*^9, 3.4101849326813965`*^9}, {3.410185362499445*^9, 3.4101853660245132`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"0.11735283061787535`", ",", "0.7069727336489029`"}], "}"}]], "Output", CellChangeTimes->{{3.4101735008131776`*^9, 3.410173536975176*^9}, 3.410173796588482*^9, 3.410174176685034*^9, 3.410174210313389*^9, 3.4101745473079634`*^9, 3.4101749252914767`*^9, 3.4101755866124096`*^9, 3.4101756621109715`*^9, 3.4101761733661203`*^9, 3.4101776528535185`*^9, { 3.4101849159473343`*^9, 3.410184933412448*^9}, 3.4101853668757377`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"FindRoot", "[", RowBox[{ RowBox[{"{", "\[IndentingNewLine]", RowBox[{ RowBox[{"x1", "\[Equal]", " ", RowBox[{ RowBox[{"\[Phi]", " ", FractionBox[ RowBox[{ RowBox[{"x1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"x2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"x3", " ", RowBox[{"M", "[", "3", "]"}]}]}], RowBox[{ RowBox[{"y1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"y2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"y3", " ", RowBox[{"M", "[", "3", "]"}]}]}]], "y1"}], "+", RowBox[{ RowBox[{"(", RowBox[{"1", "-", RowBox[{"\[Phi]", FractionBox[ RowBox[{ RowBox[{"x1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"x2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"x3", " ", RowBox[{"M", "[", "3", "]"}]}]}], RowBox[{ RowBox[{"y1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"y2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"y3", " ", RowBox[{"M", "[", "3", "]"}]}]}]]}]}], ")"}], " ", "x1"}], "-", RowBox[{"Da", " ", FractionBox[ RowBox[{ RowBox[{"x1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"x2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"x3", " ", RowBox[{"M", "[", "3", "]"}]}]}], RowBox[{ RowBox[{"x1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"x2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"x3", " ", RowBox[{"M", "[", "3", "]"}]}]}]], RowBox[{"(", RowBox[{ RowBox[{"x1", " ", "x2"}], "-", RowBox[{"x3", "/", "Keq"}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "x1"}], ")"}]}]}]}], ",", RowBox[{"x2", "\[Equal]", " ", RowBox[{ RowBox[{"\[Phi]", " ", FractionBox[ RowBox[{ RowBox[{"x1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"x2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"x3", " ", RowBox[{"M", "[", "3", "]"}]}]}], RowBox[{ RowBox[{"y1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"y2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"y3", " ", RowBox[{"M", "[", "3", "]"}]}]}]], "y2"}], "+", RowBox[{ RowBox[{"(", RowBox[{"1", "-", RowBox[{"\[Phi]", FractionBox[ RowBox[{ RowBox[{"x1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"x2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"x3", " ", RowBox[{"M", "[", "3", "]"}]}]}], RowBox[{ RowBox[{"y1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"y2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"y3", " ", RowBox[{"M", "[", "3", "]"}]}]}]]}]}], ")"}], " ", "x2"}], "-", RowBox[{"Da", " ", FractionBox[ RowBox[{ RowBox[{"x1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"x2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"x3", " ", RowBox[{"M", "[", "3", "]"}]}]}], RowBox[{ RowBox[{"x1", " ", RowBox[{"M", "[", "1", "]"}]}], "+", RowBox[{"x2", " ", RowBox[{"M", "[", "2", "]"}]}], "+", RowBox[{"x3", " ", RowBox[{"M", "[", "3", "]"}]}]}]], RowBox[{"(", RowBox[{ RowBox[{"x1", " ", "x2"}], "-", RowBox[{"x3", "/", "Keq"}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "x2"}], ")"}]}]}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"x1", "+", "x2", "+", "x3"}], "==", "1"}], ",", "\[IndentingNewLine]", RowBox[{"y1", "==", RowBox[{"x1", "/", RowBox[{"(", RowBox[{"x1", "+", RowBox[{ RowBox[{ SubscriptBox["\[Alpha]", "BC"], "/", SubscriptBox["\[Alpha]", "AC"]}], " ", "x2"}], "+", RowBox[{ RowBox[{"1", "/", SubscriptBox["\[Alpha]", "AC"]}], " ", "x3"}]}], ")"}]}]}], ",", "\[IndentingNewLine]", RowBox[{"y2", "==", RowBox[{"x2", "/", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["\[Alpha]", "AC"], "/", SubscriptBox["\[Alpha]", "BC"]}], "x1"}], "+", "x2", "+", RowBox[{ RowBox[{"1", "/", SubscriptBox["\[Alpha]", "BC"]}], " ", "x3"}]}], ")"}]}]}], ",", "\[IndentingNewLine]", RowBox[{"y3", "==", RowBox[{"x3", "/", RowBox[{"(", RowBox[{ RowBox[{ SubscriptBox["\[Alpha]", "AC"], " ", "x1"}], "+", RowBox[{ SubscriptBox["\[Alpha]", "BC"], " ", "x2"}], "+", "x3"}], ")"}]}]}]} ], "}"}], ",", RowBox[{"{", RowBox[{"x3", ",", "0.7"}], "}"}], ",", RowBox[{"{", RowBox[{"x1", ",", "0.1"}], "}"}], ",", RowBox[{"{", RowBox[{"x2", ",", "0.2"}], "}"}], ",", RowBox[{"{", RowBox[{"y1", ",", "0.1"}], "}"}], ",", RowBox[{"{", RowBox[{"y2", ",", "0.3"}], "}"}], ",", RowBox[{"{", RowBox[{"y3", ",", "0.6"}], "}"}], ",", RowBox[{"MaxIterations", "\[Rule]", "1000000"}]}], "]"}]], "Input", CellChangeTimes->{{3.410185074555402*^9, 3.4101850785311184`*^9}, { 3.410185109185197*^9, 3.4101852192134094`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"x3", "\[Rule]", "0.17567443573060715`"}], ",", RowBox[{"x1", "\[Rule]", "0.11735283061529864`"}], ",", RowBox[{"x2", "\[Rule]", "0.7069727336540942`"}], ",", RowBox[{"y1", "\[Rule]", "0.2035003628959231`"}], ",", RowBox[{"y2", "\[Rule]", "0.7355725827922662`"}], ",", RowBox[{"y3", "\[Rule]", "0.06092705431181062`"}]}], "}"}]], "Output", CellChangeTimes->{3.410185245841699*^9, 3.4101853691490064`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Chop", "[", RowBox[{ RowBox[{ RowBox[{"x1", " ", "x2"}], "-", RowBox[{"x3", "/", "Keq"}]}], "/.", "sol"}], "]"}]], "Input", CellChangeTimes->{{3.410175699104165*^9, 3.410175725652339*^9}}], Cell[BoxData[ RowBox[{"-", "0.004871966403254893`"}]], "Output", CellChangeTimes->{{3.4101757023588448`*^9, 3.410175726273232*^9}, 3.410176190741104*^9, 3.4101776573199406`*^9, 3.4101849418946447`*^9, 3.410185376629763*^9}] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["\<\ Stripping & Rectifying Profiles at Damk\[ODoubleDot]hler number Da=11.1\ \>", "Subsubtitle", CellChangeTimes->{ 3.410091630290366*^9, {3.410092930199544*^9, 3.4100929425773425`*^9}, { 3.410093659297936*^9, 3.4100937216576047`*^9}, {3.4100937524218416`*^9, 3.4100937804120893`*^9}, {3.410093953370792*^9, 3.410093953751339*^9}, { 3.410173188454027*^9, 3.4101731977473907`*^9}, {3.410174895007931*^9, 3.4101748974714737`*^9}, {3.4101774943756385`*^9, 3.410177501175416*^9}}, Background->RGBColor[0.996109, 0.996109, 0.640635]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Show", "[", RowBox[{"plt1", ",", "plt2"}], "]"}]], "Input", CellChangeTimes->{{3.4101757384807854`*^9, 3.410175744549512*^9}}], Cell[BoxData[ GraphicsBox[{{{}, {Hue[0.67, 0.6, 0.6], PointBox[{{0.4, 0.4}, {0.503594279051528, 0.40982616451610426`}, { 0.5612885402024901, 0.3574724874974058}, {0.6257823541156673, 0.29653080369536766`}, {0.6973562737431916, 0.2316487711972538}, { 0.7717784276704847, 0.16768972755127018`}, {0.8418828600892818, 0.11108411710985212`}, {0.8999081304977128, 0.06730337278742694}, { 0.9415711279522833, 0.037868818643307725`}, {0.967900917646925, 0.020262368947436226`}, {0.983041004379155, 0.010529516460701263`}, { 0.9912401241724712, 0.005387624039712777}, {0.9955303439501497, 0.0027349953373423407`}, {0.997733976091146, 0.0013828955634971496`}, { 0.9988549515333836, 0.0006978383188074628}, {0.9994223625980388, 0.00035179063471812527`}, {0.9997088490061523, 0.0001772533424359455}, { 0.9998533117895213, 0.00008928823189032727}, {0.999926111247593, 0.00004497161133667719}, {0.9999627853290821, 0.00002264929544413429}, { 0.9999812575562045, 0.000011406616930786974`}, {0.9999905609941331, 5.7444958794818506`*^-6}, {0.9999952464260604, 2.8929665684137027`*^-6}, {0.9999976060713162, 1.4569112884354696`*^-6}, {0.9999987944072404, 7.33705704405612*^-7}, { 0.9999993928593909, 3.6949641295177665`*^-7}, {0.9999996942422003, 1.8607941320844104`*^-7}, {0.9999998460195405, 9.371008101727996*^-8}, { 0.9999999224550395, 4.7192636490002857`*^-8}, {0.9999999609481603, 2.3766330838558073`*^-8}, {0.9999999803333962, 1.1968782971801074`*^-8}, {0.9999999900958493, 6.027508921973992*^-9}, { 0.999999995012245, 3.0354685056733018`*^-9}, {0.9999999974881543, 1.5286694976795835`*^-9}, {0.9999999987350283, 7.698417648017802*^-10}, {0.9999999993629571, 3.876942292046757*^-10}, { 0.9999999996791836, 1.9524377891995226`*^-10}, {0.9999999998384361, 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}, {1., 0.}}]}, {}}, {{}, {RGBColor[1, 0, 0], PointBox[CompressedData[" 1:eJw10gs01dkeB3B55B1iypQucXKIRNGM54+QwjGNETKG8qiphiKPqByixBRT 5HG851J6SKiuPP7HMxEi71ee509FepoZl9y991l3r7XXXp+1fmv9vr+992av k46+ggICAs1o45OTjtcA/P8sULv5OZ3ugk2LQQYry4OgPvnNrRj1TnBoXVYP HR2BhPmA2fTQdjCmGJL/1RkDrbJKejjlOYxZXXnF1R+HkZVwNutQM+Tbn6Sq JCZAd/+yY0l+EzgcqrYPrpiAfOqKGNfyKdy8PKRfYjcJo47KXfPqjdBzg1mU 2DAJxrfawMqsAZxerTR1a06BvPrfBzQj6iH9Z4P9u9hT4PFBvGD9qzp4rf7Q N6NpCk6XsFpvu9SB76aZTpYID2rm+1tXeLVQXxy/T9mIB1TA2Bani7UwbSzi cPEoD5iXWHdKdtSCpIqr6P0EHii6Gepw39aAupzObOIDHlwZuCXCuV8D/krz KV/aeHDM/cmTirAaOLGx49mlaR7EMlsmx+xqYK1Xy4EXyzxos7ZVEmbUgO8a rygrORombTjX9q6qgbOP3bfvU6VBTedJ0mwjF+TqlPT09GiYU4msmPqDC5Kx DjxLMxrExWMdT3txoWq880dqHw2/qUdbPDHkQjKruaf7JxrW+JS++esbLoRp +cQ/dKfhBqO84/cFCtxck7xjfGhgWa0zmRmkIJzTvenICRosN6rdT6qjQDnJ Qig0gIbENc1fGEUUaAa0JnaEoH7yto3mHAqSm0cWL5+lwfPDneLUOApaqr+7 W8am4VD217CpcAo0GgTSD0fTwNBfXfLGjwKxdtm8vEs0pPr2hQR7UaAqJNAa GUfDbbuKCGdXCq65L8VLXKGhyMlzu+8PFGw1lpzdm0CDfKYc77wNBYK1Mk1W f9AwUF7hEGWO8j0rmRK5TsPNAxHjPxlRIHpEViUpiQZJ+SDPEX0KbvgJxXxM pmFIP9dfTJcC2R8TfTRT0Hz1TaZ1WhQMchSCTVJpcPx1s8Eik4K8yvtD29No eOG9qySPQYHHd4zdouk0MHXin97dTIFRQuVUI3JOW6+UtDLyapMcfw4NfzYO nKpUoiDFwuVvoQwaMsu7JO9voODlQrHWZeT++qWVYUUKrMfu7FhGNvC3PWix noKv3bqtPpk0/Jq3+KH7GwqONURH1yBfYR13u65AQdiA5WO5LBoe5bTXB8lT cMeQSj+IvJQWPBe6lgKXqxmQjryh3EwhTY6C3HurOJ3IzOzVU52yFDDcDdlC 2TRc8Go+rorcNUPf2468ZWKPdKwMBQFZEy3OyJ6vF88tr6FgVXCVXxiy3XPz 4mjkExJfFNOw5yN61iO/T/xsXoacpWPsXC6N5o1vsn+OPJC79NQbOUelMHUM +VrNBK2IPLx4MPATcknZcHqvFAXnBD0YwjmoXsx/KgN5T4Hg3FrkbYHvp48i R1c7Casgd/5lOWqIzOZ+ztDC7hzbKIeccrhr1ABZkgNCc5IUCCf3L5ghuw9c lWpDNt0XKmSDzNt0k3qArGCty3BAXm63vJ+KbDwtEuSEfHvoxLkoZOruhIQb cqq2cI8fcuvMp2UP5GehyX7uyKsYi57eyC9KtYRZyJ56f1gcRZZ55+sByH1b j5ccRw4ZTdi7E1la92qVH/LXjpFkDeTwukj/U8hzYs82KCNHh4V3BiKbrD/6 eB3yg4cx74KQu7nvbWSQhUZz20KQD50zKxNDjmrkBZ9BFm/1GhVE3hbEmQtD flpp+3hZggL3FZbRWeRW8+OMf5DvJc0eOoccyAxU/IJ8Ueml73nkzPNroj8g DzcN2kcgj7LvOb9DnigI+JaNvKBtE/sWWTwzqx37QZWhxGvkC8NSJyORY1jr e2nkDGnlZWx53diJKeTrNu/Co5CZ+7ZoTCL3SLi8wbZaeyR3HPkf+xX7C8hT 7eetx5AF/+Oej22RZrt2FLlau+IDduisueArZLlSk13RyP6hbnIjyAM104HY pt2rTYaRS5XzCrEPSP+LPYRsvVe/H9uwsLtvEPc7dkAwBjlepsQS21shn4n9 r50XagaQU5OK9mI/N3axx2aLzvtii+ie4fUjx5n2sbFLAj3jsY87PUvBhpUy Q+zq7d/exa50sPrYh6xUu64au/HlqjLsWC2RNuzaz0vh2KGNKkPY7PjD+7B9 lrOnsRNKT2/CVmUUf8Tukbq60IusoHJ2CdvIakM39qSkhshFnHcs7RH2cl2v FHb8wRscbEmPfHnsEY/bF7DV8u58iy3iEu6HffgXAWVsUdFyN2xW0xM17Dav jbbYPvUzTGypSF9j7H7NLC1s8UIPHexUQZ4OdlRMjyr2yYpWPex+n2xF7A2O nvrYuWnXZbBnY2/uwlYfSRTFfrx053vsy0OhAthhTkFG2G8fbV/sQXb+SdwE +15ezmdsjb9OmmKzrpbOY5skFJlhF3YceIvt+/wpYO+mz0xjN/dxzbErtaWn sAu8Mi2wfYpFxrHNiz13Y2/94vwK2ypPyhJ7R8/CEHZJUwGx7f7hAezrR7ZZ YU/4r/Rhx7ELiEdUD/Zinx6QscYe1p7sxl4v6k98WTitC5vpV0t8qz7wJXak jeQe7FT5k53YPUX2xC+U4jqw+1MvEbueqX6BPT72hPijqjixlt00cXbAsXbs 8gsyNtiJosNt2H3KO4mZbR7EfrcdiTVm5luxZ1L8+fUPrxLLOsYS/0btIk78 PovYp+r1c/IeWQ+Iu+ILiAeba4lPMY8Rq610EMcV7SRuC3lFfM1RmF/v/YZY 23uwhdzXp8/ESeaPiFNZX4kVPZKJy26s3outBaHE6m+liROCfyEuPaZAfD1n D/EnzQ3EZUd2EJdYKBMfilAhNmhQI376pyxxch2T+GWqIPE7Oy1iV42FZuxC Px1iyy1vibN26hH3+owT7y/cyc9X108s1WNAHKvUSSxb9x3xR+cW4oIQQ+IW lwbi3f8YEbuJcYlzfjAhjnKqIG6NMCUO2fqY2Op3M+JT4aXEgmFAfHZXMXG8 nTmxjuM94kURC+LMutvE8nf5Xhdxi3iP6W5i4XMFxFUU360P/018WM+S+Obm P5v5/4Xvwqpc4o73fMdF5hD/Zm5FzPk1m3juEt+lp7KIf27kW/NGJnH+Et+1 LzOIM7ZZ8/Oq8/3FlW+NBA7xQTbf2qJ8u+bxnZOYTnyay3fDFr6VBvle/SyN uPwD3/8DU01O6Q== "]]}, {}}}, AspectRatio->1, Axes->True, AxesOrigin->{0, 0}, Epilog->LineBox[{{1, 0}, {0, 1}}], PlotRange->{{0, 1}, {0, 1}}, PlotRangeClipping->True, PlotRangePadding->{Automatic, Automatic}]], "Output", CellChangeTimes->{3.410175745330635*^9, 3.410176192042976*^9, 3.4101776625074*^9, 3.4101849473224497`*^9, 3.410185378983147*^9}] }, Open ]] }, Closed]] }, Open ]] }, WindowSize->{615, 523}, WindowMargins->{{3, Automatic}, {Automatic, 26}}, 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, 567, 14, 247, "Title"], Cell[1160, 39, 516, 13, 121, "Subsubtitle"], Cell[CellGroupData[{ Cell[1701, 56, 272, 4, 44, "Subsubtitle", CellGroupingRules->{GroupTogetherGrouping, 10000.}], Cell[1976, 62, 449, 11, 31, "Input", CellGroupingRules->{GroupTogetherGrouping, 10000.}] }, Closed]], Cell[CellGroupData[{ Cell[2462, 78, 213, 3, 44, "Subsubtitle"], Cell[2678, 83, 351, 9, 31, "Input"] }, Closed]], Cell[CellGroupData[{ Cell[3066, 97, 256, 4, 44, "Subsubtitle"], Cell[3325, 103, 190, 4, 31, "Input"] }, Closed]], Cell[CellGroupData[{ Cell[3552, 112, 503, 9, 44, "Subsubtitle"], Cell[4058, 123, 1000, 30, 72, "Input"], Cell[5061, 155, 408, 7, 31, "Input"], Cell[CellGroupData[{ Cell[5496, 167, 13795, 327, 663, "Input"], Cell[CellGroupData[{ Cell[19316, 498, 205, 3, 23, "Print"], Cell[19524, 503, 207, 3, 23, "Print"], Cell[19734, 508, 205, 3, 23, "Print"], Cell[19942, 513, 205, 3, 23, "Print"], Cell[20150, 518, 207, 3, 23, "Print"], Cell[20360, 523, 207, 3, 23, "Print"], Cell[20570, 528, 207, 3, 23, "Print"], Cell[20780, 533, 207, 3, 23, "Print"], Cell[20990, 538, 205, 3, 23, "Print"], Cell[21198, 543, 206, 3, 23, "Print"], Cell[21407, 548, 208, 3, 23, "Print"], Cell[21618, 553, 208, 3, 23, "Print"], Cell[21829, 558, 208, 3, 23, "Print"], Cell[22040, 563, 206, 3, 23, "Print"], Cell[22249, 568, 208, 3, 23, "Print"], Cell[22460, 573, 208, 3, 23, "Print"], Cell[22671, 578, 208, 3, 23, "Print"], Cell[22882, 583, 208, 3, 23, "Print"], Cell[23093, 588, 206, 3, 23, "Print"], Cell[23302, 593, 208, 3, 23, "Print"], Cell[23513, 598, 206, 3, 23, "Print"], Cell[23722, 603, 206, 3, 23, "Print"], Cell[23931, 608, 206, 3, 23, "Print"], Cell[24140, 613, 208, 3, 23, "Print"], Cell[24351, 618, 206, 3, 23, "Print"], Cell[24560, 623, 206, 3, 23, "Print"], Cell[24769, 628, 208, 3, 23, "Print"], Cell[24980, 633, 208, 3, 23, "Print"], Cell[25191, 638, 208, 3, 23, "Print"], Cell[25402, 643, 208, 3, 23, "Print"], Cell[25613, 648, 208, 3, 23, "Print"], Cell[25824, 653, 208, 3, 23, "Print"], Cell[26035, 658, 208, 3, 23, "Print"], Cell[26246, 663, 208, 3, 23, "Print"], Cell[26457, 668, 208, 3, 23, "Print"], Cell[26668, 673, 208, 3, 23, "Print"], Cell[26879, 678, 206, 3, 23, "Print"], Cell[27088, 683, 208, 3, 23, "Print"], Cell[27299, 688, 208, 3, 23, "Print"], Cell[27510, 693, 208, 3, 23, "Print"], Cell[27721, 698, 208, 3, 23, "Print"], Cell[27932, 703, 208, 3, 23, "Print"], Cell[28143, 708, 206, 3, 23, "Print"], Cell[28352, 713, 208, 3, 23, "Print"], Cell[28563, 718, 206, 3, 23, "Print"], Cell[28772, 723, 206, 3, 23, "Print"], Cell[28981, 728, 206, 3, 23, "Print"], Cell[29190, 733, 208, 3, 23, "Print"], Cell[29401, 738, 208, 3, 23, "Print"], Cell[29612, 743, 208, 3, 23, "Print"], Cell[29823, 748, 206, 3, 23, "Print"], Cell[30032, 753, 206, 3, 23, "Print"], Cell[30241, 758, 208, 3, 23, "Print"], Cell[30452, 763, 206, 3, 23, "Print"], Cell[30661, 768, 205, 3, 23, "Print"], Cell[30869, 773, 206, 3, 23, "Print"], Cell[31078, 778, 206, 3, 23, "Print"], Cell[31287, 783, 208, 3, 23, "Print"], Cell[31498, 788, 208, 3, 23, "Print"], Cell[31709, 793, 208, 3, 23, "Print"] }, Open ]], Cell[31932, 799, 445, 7, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[32414, 811, 1185, 34, 72, "Input"], Cell[33602, 847, 3264, 51, 374, "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[36915, 904, 450, 8, 44, "Subsubtitle"], Cell[37368, 914, 1000, 30, 72, "Input"], Cell[38371, 946, 554, 9, 31, "Input"], Cell[CellGroupData[{ Cell[38950, 959, 11172, 263, 580, "Input"], Cell[CellGroupData[{ Cell[50147, 1226, 381, 5, 23, "Print"], Cell[50531, 1233, 383, 5, 23, "Print"], Cell[50917, 1240, 381, 5, 23, "Print"], Cell[51301, 1247, 380, 5, 23, "Print"], Cell[51684, 1254, 383, 5, 23, "Print"], Cell[52070, 1261, 381, 5, 23, "Print"], Cell[52454, 1268, 380, 5, 23, "Print"], Cell[52837, 1275, 383, 5, 23, "Print"], Cell[53223, 1282, 383, 5, 23, "Print"], Cell[53609, 1289, 382, 5, 23, "Print"], Cell[53994, 1296, 384, 5, 23, "Print"], Cell[54381, 1303, 384, 5, 23, "Print"], Cell[54768, 1310, 382, 5, 23, "Print"], Cell[55153, 1317, 384, 5, 23, "Print"], Cell[55540, 1324, 382, 5, 23, "Print"], Cell[55925, 1331, 384, 5, 23, "Print"], Cell[56312, 1338, 384, 5, 23, "Print"], Cell[56699, 1345, 384, 5, 23, "Print"], Cell[57086, 1352, 384, 5, 23, "Print"], Cell[57473, 1359, 384, 5, 23, "Print"], Cell[57860, 1366, 384, 5, 23, "Print"], Cell[58247, 1373, 382, 5, 23, "Print"], Cell[58632, 1380, 382, 5, 23, "Print"], Cell[59017, 1387, 384, 5, 23, "Print"], Cell[59404, 1394, 384, 5, 23, "Print"], Cell[59791, 1401, 384, 5, 23, "Print"], Cell[60178, 1408, 382, 5, 23, "Print"], Cell[60563, 1415, 382, 5, 23, "Print"], Cell[60948, 1422, 384, 5, 23, "Print"], Cell[61335, 1429, 384, 5, 23, "Print"], Cell[61722, 1436, 382, 5, 23, "Print"], Cell[62107, 1443, 384, 5, 23, "Print"], Cell[62494, 1450, 382, 5, 23, "Print"], Cell[62879, 1457, 384, 5, 23, "Print"], Cell[63266, 1464, 384, 5, 23, "Print"], Cell[63653, 1471, 382, 5, 23, "Print"], Cell[64038, 1478, 384, 5, 23, "Print"], Cell[64425, 1485, 381, 5, 23, "Print"], Cell[64809, 1492, 382, 5, 23, "Print"], Cell[65194, 1499, 382, 5, 23, "Print"], Cell[65579, 1506, 384, 5, 23, "Print"], Cell[65966, 1513, 382, 5, 23, "Print"], Cell[66351, 1520, 384, 5, 23, "Print"], Cell[66738, 1527, 382, 5, 23, "Print"], Cell[67123, 1534, 382, 5, 23, "Print"], Cell[67508, 1541, 382, 5, 23, "Print"], Cell[67893, 1548, 384, 5, 23, "Print"], Cell[68280, 1555, 382, 5, 23, "Print"], Cell[68665, 1562, 382, 5, 23, "Print"], Cell[69050, 1569, 382, 5, 23, "Print"], Cell[69435, 1576, 384, 5, 23, "Print"], Cell[69822, 1583, 382, 5, 23, "Print"], Cell[70207, 1590, 382, 5, 23, "Print"], Cell[70592, 1597, 382, 5, 23, "Print"], Cell[70977, 1604, 384, 5, 23, "Print"], Cell[71364, 1611, 382, 5, 23, "Print"], Cell[71749, 1618, 384, 5, 23, "Print"], Cell[72136, 1625, 384, 5, 23, "Print"], Cell[72523, 1632, 382, 5, 23, "Print"], Cell[72908, 1639, 384, 5, 23, "Print"] }, Open ]], Cell[73307, 1647, 947, 15, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[74291, 1667, 1135, 28, 52, "Input"], Cell[75431, 1698, 3714, 57, 374, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[79182, 1760, 356, 9, 31, "Input"], Cell[79541, 1771, 542, 9, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[80120, 1785, 6123, 178, 337, "Input"], Cell[86246, 1965, 525, 10, 50, "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[86820, 1981, 605, 10, 44, "Subsubtitle"], Cell[CellGroupData[{ Cell[87450, 1995, 154, 3, 31, "Input"], Cell[87607, 2000, 6040, 91, 374, "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[93696, 2097, 498, 9, 44, "Subsubtitle"], Cell[94197, 2108, 1000, 30, 72, "Input"], Cell[95200, 2140, 357, 7, 31, "Input"], Cell[CellGroupData[{ Cell[95582, 2151, 13795, 327, 663, "Input"], Cell[CellGroupData[{ Cell[109402, 2482, 228, 3, 23, "Print"], Cell[109633, 2487, 229, 3, 23, "Print"], Cell[109865, 2492, 231, 3, 23, "Print"], Cell[110099, 2497, 229, 3, 23, "Print"], Cell[110331, 2502, 231, 3, 23, "Print"], Cell[110565, 2507, 231, 3, 23, "Print"], Cell[110799, 2512, 229, 3, 23, "Print"], Cell[111031, 2517, 229, 3, 23, "Print"], Cell[111263, 2522, 229, 3, 23, "Print"], Cell[111495, 2527, 232, 3, 23, "Print"], Cell[111730, 2532, 232, 3, 23, "Print"], Cell[111965, 2537, 230, 3, 23, "Print"], Cell[112198, 2542, 230, 3, 23, "Print"], Cell[112431, 2547, 232, 3, 23, "Print"], Cell[112666, 2552, 230, 3, 23, "Print"], Cell[112899, 2557, 230, 3, 23, "Print"], Cell[113132, 2562, 230, 3, 23, "Print"], Cell[113365, 2567, 230, 3, 23, "Print"], Cell[113598, 2572, 232, 3, 23, "Print"], Cell[113833, 2577, 232, 3, 23, "Print"], Cell[114068, 2582, 232, 3, 23, "Print"], Cell[114303, 2587, 232, 3, 23, "Print"], Cell[114538, 2592, 232, 3, 23, "Print"], Cell[114773, 2597, 230, 3, 23, "Print"], Cell[115006, 2602, 230, 3, 23, "Print"], Cell[115239, 2607, 232, 3, 23, "Print"], Cell[115474, 2612, 230, 3, 23, "Print"], Cell[115707, 2617, 232, 3, 23, "Print"], Cell[115942, 2622, 230, 3, 23, "Print"], Cell[116175, 2627, 232, 3, 23, "Print"], Cell[116410, 2632, 230, 3, 23, "Print"], Cell[116643, 2637, 229, 3, 23, "Print"], Cell[116875, 2642, 230, 3, 23, "Print"], Cell[117108, 2647, 230, 3, 23, "Print"], Cell[117341, 2652, 230, 3, 23, "Print"], Cell[117574, 2657, 230, 3, 23, "Print"], Cell[117807, 2662, 230, 3, 23, "Print"], Cell[118040, 2667, 232, 3, 23, "Print"], Cell[118275, 2672, 230, 3, 23, "Print"], Cell[118508, 2677, 232, 3, 23, "Print"], Cell[118743, 2682, 232, 3, 23, "Print"], Cell[118978, 2687, 230, 3, 23, "Print"], Cell[119211, 2692, 230, 3, 23, "Print"], Cell[119444, 2697, 230, 3, 23, "Print"], Cell[119677, 2702, 232, 3, 23, "Print"], Cell[119912, 2707, 230, 3, 23, "Print"], Cell[120145, 2712, 230, 3, 23, "Print"], Cell[120378, 2717, 232, 3, 23, "Print"], Cell[120613, 2722, 232, 3, 23, "Print"], Cell[120848, 2727, 230, 3, 23, "Print"], Cell[121081, 2732, 232, 3, 23, "Print"], Cell[121316, 2737, 230, 3, 23, "Print"], Cell[121549, 2742, 232, 3, 23, "Print"], Cell[121784, 2747, 229, 3, 23, "Print"], Cell[122016, 2752, 232, 3, 23, "Print"], Cell[122251, 2757, 232, 3, 23, "Print"], Cell[122486, 2762, 230, 3, 23, "Print"], Cell[122719, 2767, 232, 3, 23, "Print"], Cell[122954, 2772, 232, 3, 23, "Print"], Cell[123189, 2777, 232, 3, 23, "Print"] }, Open ]], Cell[123436, 2783, 484, 8, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[123957, 2796, 1185, 34, 72, "Input"], Cell[125145, 2832, 2702, 44, 374, "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[127896, 2882, 501, 9, 44, "Subsubtitle"], Cell[128400, 2893, 1000, 30, 72, "Input"], Cell[129403, 2925, 457, 8, 31, "Input"], Cell[CellGroupData[{ Cell[129885, 2937, 11228, 264, 560, "Input"], Cell[CellGroupData[{ Cell[141140, 3206, 91, 1, 23, "Print"], Cell[141234, 3209, 91, 1, 23, "Print"], Cell[141328, 3212, 93, 1, 23, "Print"], Cell[141424, 3215, 93, 1, 23, "Print"], Cell[141520, 3218, 91, 1, 23, "Print"], Cell[141614, 3221, 91, 1, 23, "Print"], Cell[141708, 3224, 91, 1, 23, "Print"], Cell[141802, 3227, 91, 1, 23, "Print"], Cell[141896, 3230, 93, 1, 23, "Print"], Cell[141992, 3233, 90, 1, 23, "Print"], Cell[142085, 3236, 94, 1, 23, "Print"], Cell[142182, 3239, 92, 1, 23, "Print"], Cell[142277, 3242, 94, 1, 23, "Print"], Cell[142374, 3245, 94, 1, 23, "Print"], Cell[142471, 3248, 92, 1, 23, "Print"], Cell[142566, 3251, 92, 1, 23, "Print"], Cell[142661, 3254, 94, 1, 23, "Print"], Cell[142758, 3257, 94, 1, 23, "Print"], Cell[142855, 3260, 94, 1, 23, "Print"], Cell[142952, 3263, 92, 1, 23, "Print"], Cell[143047, 3266, 92, 1, 23, "Print"], Cell[143142, 3269, 94, 1, 23, "Print"], Cell[143239, 3272, 92, 1, 23, "Print"], Cell[143334, 3275, 92, 1, 23, "Print"], Cell[143429, 3278, 92, 1, 23, "Print"], Cell[143524, 3281, 94, 1, 23, "Print"], Cell[143621, 3284, 92, 1, 23, "Print"], Cell[143716, 3287, 94, 1, 23, "Print"], Cell[143813, 3290, 92, 1, 23, "Print"], Cell[143908, 3293, 92, 1, 23, "Print"], Cell[144003, 3296, 94, 1, 23, "Print"], Cell[144100, 3299, 92, 1, 23, "Print"], Cell[144195, 3302, 94, 1, 23, "Print"], Cell[144292, 3305, 94, 1, 23, "Print"], Cell[144389, 3308, 94, 1, 23, "Print"], Cell[144486, 3311, 94, 1, 23, "Print"], Cell[144583, 3314, 94, 1, 23, "Print"], Cell[144680, 3317, 92, 1, 23, "Print"], Cell[144775, 3320, 94, 1, 23, "Print"], Cell[144872, 3323, 94, 1, 23, "Print"], Cell[144969, 3326, 94, 1, 23, "Print"], Cell[145066, 3329, 94, 1, 23, "Print"], Cell[145163, 3332, 94, 1, 23, "Print"], Cell[145260, 3335, 92, 1, 23, "Print"], Cell[145355, 3338, 94, 1, 23, "Print"], Cell[145452, 3341, 94, 1, 23, "Print"], Cell[145549, 3344, 94, 1, 23, "Print"], Cell[145646, 3347, 94, 1, 23, "Print"], Cell[145743, 3350, 94, 1, 23, "Print"], Cell[145840, 3353, 92, 1, 23, "Print"], Cell[145935, 3356, 92, 1, 23, "Print"], Cell[146030, 3359, 92, 1, 23, "Print"], Cell[146125, 3362, 92, 1, 23, "Print"], Cell[146220, 3365, 89, 1, 23, "Print"], Cell[146312, 3368, 92, 1, 23, "Print"], Cell[146407, 3371, 92, 1, 23, "Print"], Cell[146502, 3374, 92, 1, 23, "Print"], Cell[146597, 3377, 92, 1, 23, "Print"], Cell[146692, 3380, 94, 1, 23, "Print"], Cell[146789, 3383, 92, 1, 23, "Print"], Cell[146884, 3386, 91, 1, 23, "Print"], Cell[146978, 3389, 92, 1, 23, "Print"], Cell[147073, 3392, 94, 1, 23, "Print"], Cell[147170, 3395, 94, 1, 23, "Print"], Cell[147267, 3398, 92, 1, 23, "Print"], Cell[147362, 3401, 94, 1, 23, "Print"], Cell[147459, 3404, 91, 1, 23, "Print"], Cell[147553, 3407, 94, 1, 23, "Print"], Cell[147650, 3410, 94, 1, 23, "Print"], Cell[147747, 3413, 92, 1, 23, "Print"], Cell[147842, 3416, 92, 1, 23, "Print"], Cell[147937, 3419, 94, 1, 23, "Print"], Cell[148034, 3422, 94, 1, 23, "Print"], Cell[148131, 3425, 94, 1, 23, "Print"], Cell[148228, 3428, 92, 1, 23, "Print"], Cell[148323, 3431, 92, 1, 23, "Print"], Cell[148418, 3434, 94, 1, 23, "Print"], Cell[148515, 3437, 94, 1, 23, "Print"], Cell[148612, 3440, 94, 1, 23, "Print"], Cell[148709, 3443, 94, 1, 23, "Print"], Cell[148806, 3446, 94, 1, 23, "Print"], Cell[148903, 3449, 94, 1, 23, "Print"], Cell[149000, 3452, 92, 1, 23, "Print"], Cell[149095, 3455, 92, 1, 23, "Print"], Cell[149190, 3458, 94, 1, 23, "Print"], Cell[149287, 3461, 94, 1, 23, "Print"], Cell[149384, 3464, 91, 1, 23, "Print"], Cell[149478, 3467, 92, 1, 23, "Print"], Cell[149573, 3470, 94, 1, 23, "Print"], Cell[149670, 3473, 91, 1, 23, "Print"], Cell[149764, 3476, 92, 1, 23, "Print"], Cell[149859, 3479, 92, 1, 23, "Print"], Cell[149954, 3482, 92, 1, 23, "Print"], Cell[150049, 3485, 92, 1, 23, "Print"], Cell[150144, 3488, 94, 1, 23, "Print"], Cell[150241, 3491, 92, 1, 23, "Print"], Cell[150336, 3494, 90, 1, 23, "Print"], Cell[150429, 3497, 94, 1, 23, "Print"], Cell[150526, 3500, 92, 1, 23, "Print"], Cell[150621, 3503, 95, 1, 23, "Print"], Cell[150719, 3506, 95, 1, 23, "Print"], Cell[150817, 3509, 95, 1, 23, "Print"], Cell[150915, 3512, 95, 1, 23, "Print"], Cell[151013, 3515, 95, 1, 23, "Print"], Cell[151111, 3518, 95, 1, 23, "Print"], Cell[151209, 3521, 93, 1, 23, "Print"], Cell[151305, 3524, 92, 1, 23, "Print"], Cell[151400, 3527, 93, 1, 23, "Print"], Cell[151496, 3530, 95, 1, 23, "Print"], Cell[151594, 3533, 92, 1, 23, "Print"], Cell[151689, 3536, 93, 1, 23, "Print"], Cell[151785, 3539, 95, 1, 23, "Print"], Cell[151883, 3542, 95, 1, 23, "Print"], Cell[151981, 3545, 93, 1, 23, "Print"], Cell[152077, 3548, 93, 1, 23, "Print"], Cell[152173, 3551, 92, 1, 23, "Print"], Cell[152268, 3554, 93, 1, 23, "Print"], Cell[152364, 3557, 95, 1, 23, "Print"], Cell[152462, 3560, 95, 1, 23, "Print"], Cell[152560, 3563, 93, 1, 23, "Print"], Cell[152656, 3566, 95, 1, 23, "Print"], Cell[152754, 3569, 95, 1, 23, "Print"], Cell[152852, 3572, 92, 1, 23, "Print"], Cell[152947, 3575, 95, 1, 23, "Print"], Cell[153045, 3578, 95, 1, 23, "Print"], Cell[153143, 3581, 95, 1, 23, "Print"], Cell[153241, 3584, 95, 1, 23, "Print"], Cell[153339, 3587, 93, 1, 23, "Print"], Cell[153435, 3590, 93, 1, 23, "Print"], Cell[153531, 3593, 93, 1, 23, "Print"], Cell[153627, 3596, 95, 1, 23, "Print"], Cell[153725, 3599, 95, 1, 23, "Print"], Cell[153823, 3602, 95, 1, 23, "Print"], Cell[153921, 3605, 93, 1, 23, "Print"], Cell[154017, 3608, 93, 1, 23, "Print"], Cell[154113, 3611, 93, 1, 23, "Print"], Cell[154209, 3614, 93, 1, 23, "Print"], Cell[154305, 3617, 93, 1, 23, "Print"], Cell[154401, 3620, 93, 1, 23, "Print"], Cell[154497, 3623, 95, 1, 23, "Print"], Cell[154595, 3626, 93, 1, 23, "Print"], Cell[154691, 3629, 93, 1, 23, "Print"], Cell[154787, 3632, 93, 1, 23, "Print"], Cell[154883, 3635, 93, 1, 23, "Print"], Cell[154979, 3638, 93, 1, 23, "Print"], Cell[155075, 3641, 93, 1, 23, "Print"], Cell[155171, 3644, 95, 1, 23, "Print"], Cell[155269, 3647, 95, 1, 23, "Print"], Cell[155367, 3650, 95, 1, 23, "Print"], Cell[155465, 3653, 95, 1, 23, "Print"], Cell[155563, 3656, 95, 1, 23, "Print"], Cell[155661, 3659, 93, 1, 23, "Print"], Cell[155757, 3662, 93, 1, 23, "Print"], Cell[155853, 3665, 95, 1, 23, "Print"], Cell[155951, 3668, 93, 1, 23, "Print"], Cell[156047, 3671, 95, 1, 23, "Print"], Cell[156145, 3674, 92, 1, 23, "Print"], Cell[156240, 3677, 93, 1, 23, "Print"], Cell[156336, 3680, 93, 1, 23, "Print"], Cell[156432, 3683, 93, 1, 23, "Print"], Cell[156528, 3686, 93, 1, 23, "Print"], Cell[156624, 3689, 93, 1, 23, "Print"], Cell[156720, 3692, 95, 1, 23, "Print"], Cell[156818, 3695, 93, 1, 23, "Print"], Cell[156914, 3698, 93, 1, 23, "Print"], Cell[157010, 3701, 93, 1, 23, "Print"], Cell[157106, 3704, 95, 1, 23, "Print"], Cell[157204, 3707, 93, 1, 23, "Print"], Cell[157300, 3710, 93, 1, 23, "Print"], Cell[157396, 3713, 95, 1, 23, "Print"], Cell[157494, 3716, 93, 1, 23, "Print"], Cell[157590, 3719, 95, 1, 23, "Print"], Cell[157688, 3722, 95, 1, 23, "Print"], Cell[157786, 3725, 95, 1, 23, "Print"], Cell[157884, 3728, 95, 1, 23, "Print"], Cell[157982, 3731, 93, 1, 23, "Print"], Cell[158078, 3734, 95, 1, 23, "Print"], Cell[158176, 3737, 93, 1, 23, "Print"], Cell[158272, 3740, 93, 1, 23, "Print"], Cell[158368, 3743, 93, 1, 23, "Print"], Cell[158464, 3746, 93, 1, 23, "Print"], Cell[158560, 3749, 93, 1, 23, "Print"], Cell[158656, 3752, 93, 1, 23, "Print"], Cell[158752, 3755, 93, 1, 23, "Print"], Cell[158848, 3758, 93, 1, 23, "Print"], Cell[158944, 3761, 93, 1, 23, "Print"], Cell[159040, 3764, 95, 1, 23, "Print"], Cell[159138, 3767, 93, 1, 23, "Print"], Cell[159234, 3770, 95, 1, 23, "Print"], Cell[159332, 3773, 95, 1, 23, "Print"], Cell[159430, 3776, 95, 1, 23, "Print"], Cell[159528, 3779, 93, 1, 23, "Print"], Cell[159624, 3782, 95, 1, 23, "Print"], Cell[159722, 3785, 95, 1, 23, "Print"], Cell[159820, 3788, 93, 1, 23, "Print"], Cell[159916, 3791, 93, 1, 23, "Print"], Cell[160012, 3794, 93, 1, 23, "Print"], Cell[160108, 3797, 95, 1, 23, "Print"], Cell[160206, 3800, 95, 1, 23, "Print"], Cell[160304, 3803, 93, 1, 23, "Print"] }, Open ]], Cell[160412, 3807, 147, 3, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[160596, 3815, 1277, 30, 52, "Input"], Cell[161876, 3847, 3984, 70, 374, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[165897, 3922, 560, 11, 31, "Input"], Cell[166460, 3935, 479, 8, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[166976, 3948, 6123, 178, 337, "Input"], Cell[173102, 4128, 476, 9, 50, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[173615, 4142, 228, 6, 31, "Input"], Cell[173846, 4150, 233, 4, 30, "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[174128, 4160, 552, 9, 44, "Subsubtitle"], Cell[CellGroupData[{ Cell[174705, 4173, 154, 3, 31, "Input"], Cell[174862, 4178, 5853, 97, 374, "Output"] }, Open ]] }, Closed]] }, Open ]] } ] *) (* End of internal cache information *)