(* Content-type: application/vnd.wolfram.mathematica *)

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

(* CreatedBy='Mathematica 11.3' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     23414,        531]
NotebookOptionsPosition[     21793,        497]
NotebookOutlinePosition[     22191,        513]
CellTagsIndexPosition[     22148,        510]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"e", "[", 
     RowBox[{"x_", ",", "t_"}], "]"}], ":=", 
    RowBox[{
     RowBox[{"8", 
      RowBox[{"t", "/", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"3", "x"}], "-", "1"}], ")"}]}]}], "-", 
     RowBox[{"3", "/", 
      RowBox[{"x", "^", "2"}]}]}]}], ";"}], "\[IndentingNewLine]"}]], "Input",\

 CellChangeTimes->{{3.8432018123667707`*^9, 3.8432018833807793`*^9}, {
   3.843202321214196*^9, 3.843202352184506*^9}, {3.843202666879415*^9, 
   3.843202668083734*^9}, 3.843202711214531*^9},
 CellLabel->"In[28]:=",ExpressionUUID->"74691cb4-c4c5-4779-9e0d-fac0102f15ab"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"f", "[", 
     RowBox[{"x_", ",", "y_", ",", "t_"}], "]"}], ":=", 
    RowBox[{
     RowBox[{"8", "*", 
      RowBox[{"t", "/", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"3", "x"}], "-", "1"}], ")"}]}]}], "-", 
     RowBox[{"8", 
      RowBox[{"t", "/", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"3", "y"}], "-", "1"}], ")"}]}]}], "-", 
     RowBox[{"3", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"1", "/", 
         RowBox[{"x", "^", "2"}]}], "-", 
        RowBox[{"1", "/", 
         RowBox[{"y", "^", "2"}]}]}], ")"}]}]}]}], ";"}], 
  "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{{3.8432027068975058`*^9, 3.843202714480535*^9}, {
  3.843207766911086*^9, 3.843207777700469*^9}, {3.844024849941265*^9, 
  3.844024850844411*^9}},
 CellLabel->"In[29]:=",ExpressionUUID->"c268b59b-fbac-4796-a943-73f95b260c7a"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"g", "[", 
    RowBox[{"x_", ",", "y_", ",", "t_"}], "]"}], ":=", 
   RowBox[{
    RowBox[{"Log", "[", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"3", "x"}], "-", "1"}], ")"}], "/", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"3", "y"}], "-", "1"}], ")"}]}], "]"}], "+", 
    RowBox[{"9", "/", 
     RowBox[{"(", 
      RowBox[{"4", "*", "t", "*", "x"}], ")"}]}], "-", 
    RowBox[{"9", "/", 
     RowBox[{"(", 
      RowBox[{"4", "*", "t", "*", "y"}], ")"}]}], "-", 
    RowBox[{"3", "*", 
     RowBox[{"x", "/", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"3", "x"}], "-", "1"}], ")"}]}]}], "+", 
    RowBox[{"3", "*", 
     RowBox[{"y", "/", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"3", "y"}], "-", "1"}], ")"}]}]}]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.843201906320074*^9, 3.843202002184842*^9}, 
   3.843202669111134*^9, {3.843207610162158*^9, 3.8432076155356407`*^9}, {
   3.844024915430705*^9, 3.8440249177996283`*^9}, {3.844025282174511*^9, 
   3.844025288983574*^9}},
 CellLabel->"In[30]:=",ExpressionUUID->"baca2101-d087-454e-8544-a66806c12aff"],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"Solve", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"0", "==", 
        RowBox[{"f", "[", 
         RowBox[{"x", ",", "y", ",", "0.95"}], "]"}]}], ",", 
       RowBox[{"0", "\[Equal]", 
        RowBox[{"g", "[", 
         RowBox[{"x", ",", "y", ",", "0.95"}], "]"}]}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"x", ",", "y"}], "}"}], ",", 
     RowBox[{"x", ">", 
      RowBox[{"1", "/", "3"}]}], ",", 
     RowBox[{"y", ">", 
      RowBox[{"1", "/", "3"}]}], ",", "Reals"}], "]"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.843202257393837*^9, 3.843202288150055*^9}, {
  3.843202355126652*^9, 3.8432023644987497`*^9}, {3.843202451778936*^9, 
  3.843202455765604*^9}, {3.8432028423278437`*^9, 3.84320284726693*^9}, {
  3.8432028795426493`*^9, 3.8432028814177437`*^9}, {3.8432029298760366`*^9, 
  3.843202930154657*^9}, {3.8432054610241413`*^9, 3.843205479806954*^9}, {
  3.8432055241938868`*^9, 3.843205535829206*^9}, {3.843205639739245*^9, 
  3.843205658678726*^9}, {3.843206675849365*^9, 3.843206680878385*^9}, {
  3.8432078179093447`*^9, 3.843207818034964*^9}, {3.843208095619355*^9, 
  3.8432081251899643`*^9}, {3.8432082324354677`*^9, 3.84320823468462*^9}, {
  3.843208273438744*^9, 3.843208277756687*^9}, {3.843209901434976*^9, 
  3.843209940787924*^9}, {3.8432109228777113`*^9, 3.843210968646092*^9}},
 CellLabel->"In[31]:=",ExpressionUUID->"19b2792f-e0c1-42e3-90a8-3b01487a523f"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"solx", "=", 
   RowBox[{"Solve", "[", 
    RowBox[{
     RowBox[{"0", "\[Equal]", 
      RowBox[{"f", "[", 
       RowBox[{"x", ",", "y", ",", "0.91"}], "]"}]}], ",", "x"}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.843210355727895*^9, 3.843210369745371*^9}, {
   3.843210442494248*^9, 3.843210455866572*^9}, {3.84321049717526*^9, 
   3.84321052450104*^9}, 3.843210707960021*^9, {3.843210769409047*^9, 
   3.843210778561599*^9}, 3.844024920182876*^9, {3.844024954003654*^9, 
   3.844024959052896*^9}, 3.8739694989193172`*^9, 3.873969920143662*^9},
 CellLabel->"In[39]:=",ExpressionUUID->"c928bb18-46d0-4770-9233-e6826484c2e4"],

Cell[BoxData[
 TemplateBox[{
  "Solve", "ratnz", 
   "\"Solve was unable to solve the system with inexact coefficients. The \
answer was obtained by solving a corresponding exact system and numericizing \
the result.\"", 2, 39, 5, 30711370880127001338, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{{3.8432104980045443`*^9, 3.8432105250695467`*^9}, 
   3.843210709141021*^9, {3.84321077005908*^9, 3.843210779102561*^9}, 
   3.8432109786095123`*^9, {3.844024924843418*^9, 3.8440249818362513`*^9}, 
   3.844025292263525*^9, 3.844025325303727*^9, 3.87396942524559*^9, 
   3.87396950042561*^9, 3.8739697729640627`*^9, 3.873969920841055*^9},
 CellLabel->
  "During evaluation of \
In[39]:=",ExpressionUUID->"81dcec80-06d8-48b7-8989-bb73521703ce"]
}, Open  ]],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{{3.873969711326819*^9, 3.873969732987576*^9}},
 CellLabel->"In[33]:=",ExpressionUUID->"0529d9c1-0939-421a-8b78-97076827b2da"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"g", "[", 
     RowBox[{"x", ",", "y", ",", "0.91"}], "]"}], "/.", 
    RowBox[{"solx", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}], ",", 
   RowBox[{"{", 
    RowBox[{"y", ",", "2.2151", ",", "2.21511"}], "}"}], ",", 
   RowBox[{"WorkingPrecision", "\[Rule]", "6"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.8432103861923943`*^9, 3.8432104248005247`*^9}, {
  3.843210474804533*^9, 3.843210475162269*^9}, {3.843210712623311*^9, 
  3.84321083358747*^9}, {3.84402496667848*^9, 3.844025035288475*^9}, {
  3.8440251074170856`*^9, 3.844025136719861*^9}, {3.844025174613268*^9, 
  3.844025229730997*^9}, {3.844025302276883*^9, 3.844025316034048*^9}, {
  3.8440253536613407`*^9, 3.84402536477566*^9}, {3.873969503201144*^9, 
  3.87396959516333*^9}, {3.873969922851433*^9, 3.873970007070918*^9}},
 CellLabel->"In[51]:=",ExpressionUUID->"2200c577-8180-41c9-8244-7b58a8b08f75"],

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    TagBox[
     {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[CompressedData["
1:eJw1zn0s1HEcB/BLLNWuEhVhDetiUZL8Ufnmdrl5WCYbJZznh1hXymMeOsl5
uIvjzlMe7sL9MsnyUJ7qw+LQZmNKWq1IdZWHEaWtLsruc398991r78/7s49Z
yGWvcC0ajcb+/9b/pytuvvyODU7rvHZOSHI7TS+gK7+fFhLv1EWNba/YCcn8
hhJ/NOcSXUhM6R8C0R2RfQKSZpEYgY72Pywgjp6119DDzpvyyK4VOxGa5pqV
TUwN3jegvSZsskgIe6MC/erA9E2y1h/4Ad1fH8Yjw90Tf9B1Fy1SyeIq3zBb
bdG0USI56/PQHj03eeAq+bj/pwd6UbkrhgwzY2PQvSMGweTzuC4fPSX+5UFO
pO+VoXkjUgZ8uxr6BH3leIsPDE4OvkIb6lWHgaosdBntxGrgAjcieXuO2p6P
WXGg19hghe61f5cEW4k+Gx3k9SYNjrmUBqGnrZkZ0F7tlYY+RBvMhITNOeVo
Jz0mH8LnXreiHUpFOXDfOHgUPTCnJwCzPu1ZtF171G0YTzfXyVWb29NaAFP8
RDO0zN2wEIw3fT2B/h3LLIKmFwXn0ft4mWJIb2uOQ1s/6JFAOn0hH31PoiyG
8ZbQ++jFUEYpBIwuK9D1855lYK61OolmlcSXA+GdUqGbGK13IC+hcnee2vzY
txVwqNP2KNq2zqQKtOKmzqDZhS7VoM/aEoVmVnGlkMKJuYXOGvKRgcHSl2q1
eXelszLQHqvoxpyaEd6FIwPPJzA3dj5YA/l2+kvoIRXUAGtWvE2gdoBlYC1Y
mthaqU2zbP5RC27BbBY6zDW7DvqXxIFogTdDDlyVKhn7SuiSg8PE1+uYs+d6
5PC3+WUKmjJSyEEY2ZiG8xFxI3JoGPPLwFxp9UkOyvquHI0ldAqaMqlczb5n
OymI5xTlafYt7KFA2yBaiI5wM6fA/IaRSHPPmgMFM746hZq+zUkKmu2/o3nU
BSYFTjNDYnT4I3cKdBVtErRi2pOCEamsGPuMHT4UlFwXlqCzHP0o4HgnlWru
jQ6i4B++HbYq
       "]]},
     Annotation[#, "Charting`Private`Tag$41005#1"]& ]}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{2.215100000000005, 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{
   "DefaultBoundaryStyle" -> Automatic, 
    "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "DefaultMeshStyle" -> 
    AbsolutePointSize[6], "ScalingFunctions" -> None, 
    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{{2.2151, 2.21511}, {-6.668755752059496*^-7, 
   2.1116136483811943`*^-7}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{{3.8432104261965837`*^9, 3.8432104306709557`*^9}, {
   3.843210476402565*^9, 3.843210503995599*^9}, {3.843210713643531*^9, 
   3.843210834400278*^9}, 3.8432109789339943`*^9, 3.844024926188417*^9, {
   3.8440249905079823`*^9, 3.844025036684132*^9}, {3.844025114025807*^9, 
   3.844025137516275*^9}, {3.844025180773144*^9, 3.844025230528997*^9}, {
   3.844025292587978*^9, 3.844025325954924*^9}, {3.844025359013541*^9, 
   3.844025365690918*^9}, 3.873969425409718*^9, {3.873969503708806*^9, 
   3.87396959554205*^9}, 3.8739697730404778`*^9, {3.873969923528902*^9, 
   3.873970007608473*^9}},
 CellLabel->"Out[51]=",ExpressionUUID->"392da884-076c-4be4-94af-8e6b85d857bf"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"solx", "[", 
    RowBox[{"[", "2", "]"}], "]"}], "/.", 
   RowBox[{"{", 
    RowBox[{"y", "->", "2.2151"}], "}"}]}], "\[IndentingNewLine]"}]], "Input",\

 CellChangeTimes->{{3.8739697385902023`*^9, 3.873969767905293*^9}, {
  3.87397001768839*^9, 3.8739700206093903`*^9}},
 CellLabel->"In[52]:=",ExpressionUUID->"c5c905b6-6cb4-4f90-bee7-90eff7fb098b"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"x", "\[Rule]", "0.6161480659001872`"}], "}"}]], "Output",
 CellChangeTimes->{{3.873969755403537*^9, 3.873969773068412*^9}, 
   3.8739700217525578`*^9},
 CellLabel->"Out[52]=",ExpressionUUID->"3d395bbb-f150-4743-83d9-914bd790167c"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"e", "[", 
      RowBox[{"x", ",", "0.91"}], "]"}], ",", 
     RowBox[{"e", "[", 
      RowBox[{"0.616148", ",", "0.91"}], "]"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", "0.2", ",", "3"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.8432110213107033`*^9, 3.843211058089129*^9}, {
  3.844025378580698*^9, 3.844025379106061*^9}, {3.8739698578946457`*^9, 
  3.873969894832502*^9}, {3.8739700291980886`*^9, 3.873970049328444*^9}},
 CellLabel->"In[54]:=",ExpressionUUID->"212d17cf-b08d-4e7e-b257-d240f363e26f"],

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    TagBox[
     {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[CompressedData["
1:eJwVjX0slHEAx2+oVSwXJi8LaRWXxcoyRd+i64UkrNC8XIbQ0mWJHZlRMlde
ehuuF/NSl+rOZXT0gs0qdnSprZrieZ7f89zd5CIlWTM9/fHd94/Pvt/P2pTT
0WkWAoHAm8//LuiZFsfTDGResmBF/wTifr2vDDMzkAavXGHImECENiW3RkSw
y/569GiPCanOZ6sctxAUj1Tl1VWZIJNdfKAI5Ln4u3mDxARlkPJry34CQYx/
SJOlCRZ9ZnF3OkF136xcF2OEk+diMrIJ9KL6MbW3Eb6lQll/LoEkQdcRJjAi
QbxVNXyBwMOpYEahNqBzIN+RNPI8SuKe4miATlThl9FKINXrXnT+4EDkijCz
hvc15wlUQxyEh14WzfUSOH16OCC7xCFrxNJoPc7vnYfM4RYc3EYrn/q7svCL
vKMas2PhH9TwrsuTxbL2Y7aRcwThtzXfdopYHA4orMIXgvzkD+4HAlkI0nyz
6+4TjBDn8sRYFsp8unchlMC0Z1Mjk8SizUpmEeRDsNgS9PxEOs9XlbhNOhD4
ZCRPn8nlfSa5WGFkUDbZHFd2jYVk+eboZVcZ3IrozLFWsDA12OxQnmfQrnp9
ubqR/6uB5HEmA0o60VuvYaG3cyqhQxhsn/X1Vr9lEacsPzgwT2NqoXv+ozUH
veXPm5n5NAaP1mpG7TgI1u9teZJF4546N3PcmcMfBVW8O5FG0nG/z4aNHKTC
wTGfUBrD/S3a2VAObYOxTR52NFrXlErnwzhMz9jmFC2lUXZO4rUQxUHb/mrM
5S+FYC/XWisJ/6d942/DUXgkr86zL+RARZUaup5RKGdP+a4u5aC8W9E0q6GQ
GhxudKngffvir2QqKbhOLYn1rOVQfCSZ+X2Dgp9LRIeDksM24bqAkyUU/gE4
/VDr
       "]], LineBox[CompressedData["
1:eJwV1nk8FV8bAHBbrmRPWSpZS/IrIaJ4RopSQimKQmSJJNmSbKFNobJLKoWK
0F6o1CtkC8karnvvcCPLXWZc23v8NZ/vZ2aeOeeZ85znqJw8e/CUAB8f3zJ+
Pr7Fa/Ir6/aDdSNgIKVm6Bs7AG3/pXzIbh4BnXcefzQSBuCbkPTe0yJMkB/d
+K2+cAAKyqU8KeFMwNKWa6lODoCSgXLwXDQTkhUNA59zBuDu+81xrCtMKFX5
7G/EG4Dozwce9qei97McjZwEB8GxKfHP2zImKL+clKxYMQgiTBEH7xEmDCeG
Og4ZD4K3Cv+eOse/4L0LG3FIGASt5MkNNwxHwTGBJzWmRIXp2QhiTPEf1Jre
jL99YghWOglvI/vGIU5k/xN3bxo0aG3uy6aPw2hPuFOSHw1ieQ6xMDYO32hK
a98F0GAss6ghYRZ54UwT3wUanPR7UnpdfAI6d+a/9b9Gg23P9u5W0J6A9rpK
tYoiGtA1k/0NfCZA/+fAOcVhGoD6mi8BQxOQ8St9mOpEB09bu3IdfAKE7Kz1
MBc6JEbEP5pgTgA7ILQg6yQdutvH4gOmJoBmmiNh6UOH4PgqqwD+SdgVv4R5
KYQOzxgu7WfXTkLp+Ft/j2Q6yBU9xP2dJsG132O5YjUdxv/TlDjTNgl87/+2
H1vDAPLwzgOjHZOwOgHr6V7LAIFI51u+3ZPQadiY46jKANmmZInTg5PQa215
2Ho9A4z9pyW8xiehIn3u9aotDIgvqZd0WzYFoyorllnvYsCqzWekD5tPgbpi
ZpOIDwPUHa7YtVtMwTfnUbmdvgzYFPUg5ZDVFFTI5m2/cIYBWPMv6YN2UxDk
qr9y8BwDPM+ayti4TEHD3jD99HAGlL2QXL734hS8w/9N37vBAEudclmTl1Og
7XS8b/VzBvgX98uIvJ2C52M7e01KGJCmJS7d9mEKdjhHXnQuZQBD3VvcpxrF
J/bl3X7FgAQFJeHUn8hDuYXUCgbUCFwnRsenoNbbt1+mkQEWv9y6c7RZ8G3w
4s7Kvyj+waROLx0WRK8Zo94dQ/GbKzp09VnQqRt24PQ4il8n11a7nQW9ElsS
JVgoflXjD7YVC7CGsO87eCh+gVHlfh8WtNSoem0RwWF3uHTe7GMW7B/UvRmn
ioMtS30soYgF7NY1NnNqOBzz22YsXcwC+6z7e4M0cAg44dKu8YoFTfqPRp01
ccjaWSxiU82CMMsdIgqbcPgnuvfcgz8sIEscYkyMcMjIjt5pIccGzUcxtaYH
cHgoezepRZENIZwH+dE2OBTfLOg9psQGv+pH3Z9tcaiOago5q8GGQkUXxx2H
cBh1X/0sQ48N+c3r3651xMFM+93yvzZssM/zpGW44cCsGKclXWXDsBln1uk8
DstvrJw4lciGkz21OTFBOOw4ajKzPZkNrVHxgY+DcUjiXJceTmcDpvb0Gx6K
g/6mdaZYARsG0mTWOUTgcOm+c/rE/9jgGSQ2PhiHg1Rs3V47QQ54sE6LlN7F
wch2wn49hQPOvs1Vuak4nFSSc50T5cDVV7E/b6Th8OqDR0iRDAfVe7iUWwYO
R1kLD/lVOfBGWnIFmY3m72EwW4ZxwECnzqjrIQ5bLR+VyERyIMjcrOvsCxye
dvl0dsRwgHb+J1WzFAdlPx2B7HgOVLyMCh9AFk2pPKx2kwNsGTvHfeU4/On+
PauXw4HEaTem5GscEvzF9h3+wIGBkU3oisMMX1uQQhUHIk2TSwnkgDuZuX1f
OKDzyfBVxkf0f9+tmzpVh8ajZ5PVUYHDZgGzzJBODjAOGc2ZfcKhMzUYT+dy
IFmbTv7+ioO15g5pZx4HIhr3vfP4hv7XB4HtyvMcOELe1JhAft6fdKtwCRd2
7Q8RF6zBIWbD063vZbmQ9yXgP4VaHLQq/1zu1uUC3/2FaEoDDvdtHhfnGnCh
/d3FpChkWarv75PGXMA2/z7GRZ4Xntb6a8YFK7+ujQONOLTZLm+dseWCrVth
S34zDhE0S+U1Z7lQuy5+YqoVh/ZAdkZeIBeOTGtfs23DQZv/gbR6CBeEHNRW
FCP3rOEJbLzEhWHDnqXu7TgYOzynb7uB3u+ICK/+hQNZJ/XMvoALnRKH5Ww7
UT04Vqp1PuWCd5ex8T3kQoZPjlMJF9RVK51GkB0Ev948+ZoLHabO9yK7cHiz
Izgg4CsXeNvtNt3rxuH8i66tif1cCPPWUH7Wi0O9SUKJ5BAX4myXRw0jqzbo
rr/D4IJHVyyu3odD6/AN+awxLjhOp/zOQtZRNZ0pmOECe4+I5cU/OFwtZZ7f
uMAFmb57QmXIA6bpoyUCBPwR4/9JR05ymuh7vZSAbeVNGfv7UT3effj5qxwB
Yqo1XIkBHCzUbIwsVxEQ5iVqYoKcWzZTVq9EwOzkgfTTyNZN9o9+aqD7705f
r0Z+TqEk9OsRsMOvPcJrEAehtJfzJw0JeEIKBd9CdlZ3DWUYE1B+3ibhFfIy
s/feY2YE5JS+5c4hnw733TdjQ4D2Gk/lBCoa7/GNsfWHCKAOFaYVIh8x+/su
04EAZ0aAYD2yuYjvum0uBGSdevt16RAOH/5qHae4ExDH/S6viazTzLzT4UkA
32iM127k1Wmn+YP9CSj0l+JEIN++oLVtVyABorTArenIS48z/ZeHEEC7lRVY
isxVO91bdomAYw9tGP3IfhSt5TExBHg2eypxkanMkb228QQc+XjkyDIaDs1l
Pm//JRIgW9ZUq4tskbrhX2UyAb3uvUK7kSvCRtRv3iXg3/om8yPIus5FTs4Z
BBy4cSrBE7kQfG5vzCGgJdO9MRhZSW1DHe8+AR59VxXikFOFRxbqHhHgeiLV
NwVZjFlokFlAQOejPd/uIcc2ep/xfkZAyUs39SJkslQz3/AFATyZ+4kvkf3v
DncLvyTAoqZqrgKZFloo3fGGgCa7yLD/ITs5ee95/IGAb++z5hqQW001o4Kq
CMiLbbnZhrxHdfi1eTUB+8t6NnQhVy0pHJWpQfn+Edrah6w/4qVGrSOA/O1y
dRD5acP6Y2WNKN8nDuynISuX4snRPwnQ9ZlTYiCn3Sn4bvOLAHbsOn4cWTzU
a16pi4AT72NZi/fjjq3f+q8X/e+cei4dmWeC+1YOEHDQ9M2yxXgBKgUPE2kE
5OvP6S5+jyHk1eU0jMaffuTM4nich9dJbRxF60/wwsdO5LYfDAveOFpP6YpK
i/OxevHkUh2LAGFyInVxvp9ve77KIFB97CtTW8yHQci6v14zKJ8gVbuYr+dH
GSqGC2i9C1VfXsynqskTR2FBEg7Yhjss5jtT2TPplzAJiSdGzXORpYTW1eSL
krByon7vbeQEnD57XoKEp/tqfeKRA0tOnZZZScJVRQbPG3k4RePBoAIJzA/e
Z44iuwTTf5euIWEgqml2D/L+Had222iQkJcrEKyBXL1WI0JpAwkNJza7yCBv
E6SXj2mT0CE16z2P1qdGvYdyoj4JshPaf9qQ5x08ZmrNSfAU73QLQnYwFI1x
sCShM6xotTNy6cpSYYYVCXH6N+TNkU/+mpEWPEiC8N+iXxLINQfvaJq4ksDg
9wm6h+pvra5RyQ93EnZJtltEIYdJ9+sd8yLBeUn2CVdkrRYtCPUn4TnTyVUZ
+aZ19eGyCBKUN5/7mI7qnaHt3QPRJPBFx38/j2wqJuHWdJmEsGMb19ggT9Qf
9ft7nQRVMaxAENl+z2SsRiYJd/WW7vZA+81qc+XSzNfo/syAghDav4JVa7Zq
vkf5l/882472tyZ+v49vKkgYvue05TFy9Oe3NW1fSTBd/j3FHHnIxKZPvJUE
Y4aDSRjaP59uixSN/UfCzpa7xV96cBCUV0+WmiJhT2IofgXZmahbcZ9DQtXW
f+HWyBJvVqhUzJIgtjK2tAPt34F6xYZc0Wk4Mn7PfRDt70abej1Or58GvisB
hzp/o/77U2UprjUN/qmpb1KRY4O9it03TUPVSf6Yg8gNFVMc563T0KniY13f
gYPrvqVXbMynoZ2V2FKO+k2Ct0GRvgt6ntfT74X6U/vD5LH5tGnI15yTPtWE
wwWLjpSLWdNQkt+yWQpZibnKgLw3DVf/2139AfVD7y0FkZP5aHyDTy0kkGc+
VUrSyqchpCqxs+QHWs99zC11TdMwdXYqrAP13wC53SF3hHnwJnOXW3U1qm+p
3E4nUR7snBXfYo9MiBDG6hI8qHV29KV/wSGSV7jwagUP2kk7FWHkxD6xax3q
PBg2Yt3H0Pmg4FF7toI5D1rZJTZZ6DzxZ7PHl/tR6P3lSk+GXqJ+suey2HOS
B1Y19ckMdN5h76iu1JvlwTuLFxquyFlb+M9+XOBB74og264HaP0oRv2sF56B
dd8Z3vV56Ptj4WkjsjMQVHi8pCAXh5Lbgcrrt8yAfrRV64EsHCT73PQf+syA
vKt4kV0KDj/PYc4Z3TPAnDA2M45E+aC26gQ1z0L4eZPhg+h8qeJ1pPj3pzlo
utZ/1kgJ9cchxyGJvHlQfDP6mo7O126FOrlb7ywAHyvxu0E5A6rKzZfKb+fD
xIynxFai839dlpfc0ad8mEfHJ8czZgwoWMfFpdX4sXmaZFOqAAM6N+1KfXWH
H/N/HIKJfadDb/esULO4ADYxx4T+eDroCwbbtUUJYFazYXmle+gg+pp+Lntc
ACMFuy/mitHhpegX+x5vQexTBH9P0g8avOh1kl3SLYgdz/6QLZ9Eg8iDbmlP
7IWwuPxg+3lrGvywcxld1iiEce9fCMiToIFf4+q1TiZLsLxVPi/k24dgYHSX
rOrbJdiu0i+6jqlDUHBVm9GqKYxlb9vHinAagtFisUTLDGFso7v/qii1IRAK
dWDrSVMwi7LhHVk4FQIq3IxEYylYrtQm3RulVPh95biUQBwFS2+p1kt/QQXT
Q0fx6XgKFuf2afejEiosY9qmjlyjYDxBbtbH51R4IodN1qZQsE6dy6JjhVTo
CVQquvKAgs1tPrVw8AEVLDb0yC/5QsEkL8lZbbxNhWJ2x/hsNQULN/FzMkqh
guzn1hr2Nwomu/7wJctkKgweqQ+i1VIw++DhKfdbVIiIe9/ytYWC6Vw3fHzv
OhXK+9Ovxg5QMJpP/HrZy1RQeHbHJZxKwfJvlOxUi6VCdEiSQSCNgr0Yi/TX
jaGCtfgVmtswBXvOUiBto6gwbByCmU2g8as6tt26SIUDwoFyRlMUzDig+3xu
OBVe/zzzT4dNwchGMY2SC1SI9T6Vo0xSsD/z8e8bQ6kwoud2Xp5HwbDxppy+
ECrYLjhbSc1SsCythuSxYCq8qXdUEZmnYN0fIlPngqiwJs2eXFigYCKF+FNx
5P8DcJPRaQ==
       "]]},
     Annotation[#, "Charting`Private`Tag$42379#1"]& ], 
    TagBox[
     {RGBColor[0.880722, 0.611041, 0.142051], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[{{0.20000005714285715`, 0.6781557463716945}, {
       0.22748015425357115`, 0.6781557463716945}, {0.25496025136428513`, 
       0.6781557463716945}, {0.31454339983928425`, 0.6781557463716945}, {
       0.33244047619047623`, 0.6781557463716945}}], LineBox[CompressedData["
1:eJxTTMoPSmViYGAwAGIQ7fAt2f5b8lV7j0iOdcWbn9pnXFq2Z5HebTi/Q8BJ
h+XHPTh/8Yw3DlEOD+F83x1JpRO1HsP57//u+nmd+ymc/8BE8OPRM8/gfJnf
gp8mbX4B5y/fJJDGXvUKztea8FGz2/wNnP/zT833t1Lv4HyxaDaLH3ffw/n2
KrIHCx5/QNivq8GXe/kjnO9usEnEdvMnON+1SnDBn6Wf4fxXe94/6e/4Aueb
ui9eJ1T3Fc6veeKuIJv/Dc7Pqsr2/u3/Hc7/F57y+4TzDzjfUu9OSpb6Tzi/
QNy1bDLbL0R4eTTzrPmB4F8sdIiZces3nP/70SWDkvN/4HzF9LC11/f/hfPf
PY54zLfgH5yfuMJgnunk/3D+vk3OnBLWDA4w/slZ6eKRqxD85WrfngsqM8L5
N/Rcpm6ZjODfufWH5TwvE5xvwlwaeLkewefa+rRw9nsEfzPXwZDbGcxw/vo7
0SKstxD8uqDEactCWOD804Hxb7jPIvg5Z2Xko21Z4fwHb1xElLYj+Ms7dJ5d
0mCD89+s5elxn4Hgs5SHfzEWZIfzC/YkWnI1IfhumrclWA8i+LLTQn78/4/g
AwC/wGVR
       "]]},
     Annotation[#, "Charting`Private`Tag$42379#2"]& ], {}}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0.20000000000000046`, 0.5159202649015313},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{
   "DefaultBoundaryStyle" -> Automatic, 
    "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "DefaultMeshStyle" -> 
    AbsolutePointSize[6], "ScalingFunctions" -> None, 
    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{{0.2, 3}, {0.5159202649015302, 0.8364998055196928}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{
  3.843210979426372*^9, {3.843211023352861*^9, 3.843211029786645*^9}, 
   3.843211065751042*^9, 3.844024926991662*^9, 3.844025293086327*^9, 
   3.8440253265188513`*^9, 3.8440253801372423`*^9, 3.873969425495941*^9, 
   3.873969773149023*^9, {3.873969866407733*^9, 3.873969895481106*^9}, {
   3.8739700431339407`*^9, 3.873970049827943*^9}},
 CellLabel->"Out[54]=",ExpressionUUID->"606beb17-67e3-4c51-8ce4-a037917823d8"]
}, Open  ]]
},
WindowSize->{808.5, 621.},
WindowMargins->{{Automatic, 277.5}, {61.5, Automatic}},
FrontEndVersion->"13.1 for Linux x86 (64-bit) (June 16, 2022)",
StyleDefinitions->"Default.nb",
ExpressionUUID->"3637d3cc-3c46-4434-89ab-ed58079b5c28"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 644, 18, 54, "Input",ExpressionUUID->"74691cb4-c4c5-4779-9e0d-fac0102f15ab"],
Cell[1205, 40, 916, 28, 54, "Input",ExpressionUUID->"c268b59b-fbac-4796-a943-73f95b260c7a"],
Cell[2124, 70, 1162, 34, 48, "Input",ExpressionUUID->"baca2101-d087-454e-8544-a66806c12aff"],
Cell[3289, 106, 1517, 30, 71, "Input",ExpressionUUID->"19b2792f-e0c1-42e3-90a8-3b01487a523f"],
Cell[CellGroupData[{
Cell[4831, 140, 670, 14, 46, "Input",ExpressionUUID->"c928bb18-46d0-4770-9233-e6826484c2e4"],
Cell[5504, 156, 765, 14, 44, "Message",ExpressionUUID->"81dcec80-06d8-48b7-8989-bb73521703ce"]
}, Open  ]],
Cell[6284, 173, 225, 3, 71, "Input",ExpressionUUID->"0529d9c1-0939-421a-8b78-97076827b2da"],
Cell[CellGroupData[{
Cell[6534, 180, 939, 18, 46, "Input",ExpressionUUID->"2200c577-8180-41c9-8244-7b58a8b08f75"],
Cell[7476, 200, 3413, 74, 208, "Output",ExpressionUUID->"392da884-076c-4be4-94af-8e6b85d857bf"]
}, Open  ]],
Cell[CellGroupData[{
Cell[10926, 279, 397, 10, 51, "Input",ExpressionUUID->"c5c905b6-6cb4-4f90-bee7-90eff7fb098b"],
Cell[11326, 291, 270, 5, 33, "Output",ExpressionUUID->"3d395bbb-f150-4743-83d9-914bd790167c"]
}, Open  ]],
Cell[CellGroupData[{
Cell[11633, 301, 614, 14, 46, "Input",ExpressionUUID->"212d17cf-b08d-4e7e-b257-d240f363e26f"],
Cell[12250, 317, 9527, 177, 261, "Output",ExpressionUUID->"606beb17-67e3-4c51-8ce4-a037917823d8"]
}, Open  ]]
}
]
*)