Some work.HEAD master

author: Jaron Kent-Dobias <jaron@kent-dobias.com> 2023-12-11 18:06:08 +0100
committer: Jaron Kent-Dobias <jaron@kent-dobias.com> 2023-12-11 18:06:08 +0100
commit: f9fc90631253da8bea10694fe0d90dd393c546e6 (patch)
tree: a86403e382f66a945d4b7974aeee9354ab5f0db5
parent: 8aeb86c75c525887119655f72a5f2ef5d68b0940 (diff)
download: mma-f9fc90631253da8bea10694fe0d90dd393c546e6.tar.gz
mma-f9fc90631253da8bea10694fe0d90dd393c546e6.tar.bz2
mma-f9fc90631253da8bea10694fe0d90dd393c546e6.zip
2 files changed, 2011 insertions, 127 deletions
diff --git a/caselle_compare.nb b/caselle_compare.nb
index a64a48c..e8c7f52 100644
--- a/caselle_compare.nb
+++ b/caselle_compare.nb
@@ -10,14 +10,16 @@
 NotebookFileLineBreakTest
 NotebookFileLineBreakTest
 NotebookDataPosition[       158,          7]
-NotebookDataLength[    204728,       4608]
-NotebookOptionsPosition[    193637,       4413]
-NotebookOutlinePosition[    194074,       4430]
-CellTagsIndexPosition[    194031,       4427]
+NotebookDataLength[    234735,       5237]
+NotebookOptionsPosition[    222110,       5018]
+NotebookOutlinePosition[    222547,       5035]
+CellTagsIndexPosition[    222504,       5032]
 WindowFrame->Normal*)
 
 (* Beginning of Notebook Content *)
 Notebook[{
+
+Cell[CellGroupData[{
 Cell[BoxData[
  RowBox[{"Get", "[", 
   "\"\<~/doc/research/first_order_singularities/mma/schofield.wl\>\"", 
@@ -28,7 +30,146 @@ Cell[BoxData[
    3.827832304560932*^9, {3.8278486734515533`*^9, 3.827848673995721*^9}, {
    3.8278488730224733`*^9, 3.827848876357965*^9}, {3.827916449928692*^9, 
    3.827916450999913*^9}},
- CellLabel->"In[1]:=",ExpressionUUID->"94557b52-cc86-4010-b22d-6b74ac4aad03"],
+ CellLabel->"In[3]:=",ExpressionUUID->"94557b52-cc86-4010-b22d-6b74ac4aad03"],
+
+Cell[BoxData[
+ TemplateBox[{
+  "\[Eta]", "shdw", 
+   "\"Symbol \\!\\(\\*RowBox[{\\\"\\\\\\\"\[Eta]\\\\\\\"\\\"}]\\) appears in \
+multiple contexts \\!\\(\\*RowBox[{\\\"{\\\", \
+RowBox[{\\\"\\\\\\\"Schofield`\\\\\\\"\\\", \\\",\\\", \
+\\\"\\\\\\\"IsingScalingFunction`\\\\\\\"\\\"}], \\\"}\\\"}]\\); definitions \
+in context \\!\\(\\*RowBox[{\\\"\\\\\\\"Schofield`\\\\\\\"\\\"}]\\) may \
+shadow or be shadowed by other definitions.\"", 2, 3, 1, 31865477087913293403,
+    "Local", "Schofield`\[Eta]"},
+  "MessageTemplate2"]], "Message", "MSG",
+ CellChangeTimes->{3.8762084299477787`*^9},
+ CellLabel->
+  "During evaluation of \
+In[3]:=",ExpressionUUID->"416ee172-1ad9-434e-b73c-e99e451d8e53"],
+
+Cell[BoxData[
+ TemplateBox[{
+  "g", "shdw", 
+   "\"Symbol \\!\\(\\*RowBox[{\\\"\\\\\\\"g\\\\\\\"\\\"}]\\) appears in \
+multiple contexts \\!\\(\\*RowBox[{\\\"{\\\", \
+RowBox[{\\\"\\\\\\\"Schofield`\\\\\\\"\\\", \\\",\\\", \
+\\\"\\\\\\\"IsingScalingFunction`\\\\\\\"\\\"}], \\\"}\\\"}]\\); definitions \
+in context \\!\\(\\*RowBox[{\\\"\\\\\\\"Schofield`\\\\\\\"\\\"}]\\) may \
+shadow or be shadowed by other definitions.\"", 2, 3, 2, 31865477087913293403,
+    "Local", "Schofield`g"},
+  "MessageTemplate2"]], "Message", "MSG",
+ CellChangeTimes->{3.876208429983996*^9},
+ CellLabel->
+  "During evaluation of \
+In[3]:=",ExpressionUUID->"ec136869-a624-43fb-87cd-ab7ce0e4f6b5"],
+
+Cell[BoxData[
+ TemplateBox[{
+  "\[Xi]", "shdw", 
+   "\"Symbol \\!\\(\\*RowBox[{\\\"\\\\\\\"\[Xi]\\\\\\\"\\\"}]\\) appears in \
+multiple contexts \\!\\(\\*RowBox[{\\\"{\\\", \
+RowBox[{\\\"\\\\\\\"Schofield`\\\\\\\"\\\", \\\",\\\", \
+\\\"\\\\\\\"IsingScalingFunction`\\\\\\\"\\\"}], \\\"}\\\"}]\\); definitions \
+in context \\!\\(\\*RowBox[{\\\"\\\\\\\"Schofield`\\\\\\\"\\\"}]\\) may \
+shadow or be shadowed by other definitions.\"", 2, 3, 3, 31865477087913293403,
+    "Local", "Schofield`\[Xi]"},
+  "MessageTemplate2"]], "Message", "MSG",
+ CellChangeTimes->{3.876208429988597*^9},
+ CellLabel->
+  "During evaluation of \
+In[3]:=",ExpressionUUID->"c71f26b5-fddb-4e00-990e-08b74adae32a"],
+
+Cell[BoxData[
+ TemplateBox[{
+  "ruleB", "shdw", 
+   "\"Symbol \\!\\(\\*RowBox[{\\\"\\\\\\\"ruleB\\\\\\\"\\\"}]\\) appears in \
+multiple contexts \\!\\(\\*RowBox[{\\\"{\\\", \
+RowBox[{\\\"\\\\\\\"Schofield`\\\\\\\"\\\", \\\",\\\", \
+\\\"\\\\\\\"IsingScalingFunction`\\\\\\\"\\\"}], \\\"}\\\"}]\\); definitions \
+in context \\!\\(\\*RowBox[{\\\"\\\\\\\"Schofield`\\\\\\\"\\\"}]\\) may \
+shadow or be shadowed by other definitions.\"", 2, 3, 4, 31865477087913293403,
+    "Local", "Schofield`ruleB"},
+  "MessageTemplate2"]], "Message", "MSG",
+ CellChangeTimes->{3.876208430021934*^9},
+ CellLabel->
+  "During evaluation of \
+In[3]:=",ExpressionUUID->"2eb6cbef-108c-4bdb-b0d3-f595ef80164f"],
+
+Cell[BoxData[
+ TemplateBox[{
+  "g$", "shdw", 
+   "\"Symbol \\!\\(\\*RowBox[{\\\"\\\\\\\"g$\\\\\\\"\\\"}]\\) appears in \
+multiple contexts \\!\\(\\*RowBox[{\\\"{\\\", \
+RowBox[{\\\"\\\\\\\"Schofield`\\\\\\\"\\\", \\\",\\\", \
+\\\"\\\\\\\"IsingScalingFunction`\\\\\\\"\\\"}], \\\"}\\\"}]\\); definitions \
+in context \\!\\(\\*RowBox[{\\\"\\\\\\\"Schofield`\\\\\\\"\\\"}]\\) may \
+shadow or be shadowed by other definitions.\"", 2, 3, 5, 31865477087913293403,
+    "Local", "Schofield`g$"},
+  "MessageTemplate2"]], "Message", "MSG",
+ CellChangeTimes->{3.87620843004827*^9},
+ CellLabel->
+  "During evaluation of \
+In[3]:=",ExpressionUUID->"2d64117a-f0ef-4d22-873c-490eeb5b19e8"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+ RowBox[{
+  RowBox[{
+   RowBox[{"g", "[", 
+    RowBox[{"\[Theta]0Cas", ",", "gsCas"}], "]"}], "[", "\[Theta]", "]"}], "//",
+   "Simplify"}]], "Input",
+ CellChangeTimes->{{3.884691107722596*^9, 3.8846911288583593`*^9}, {
+  3.884691221764353*^9, 3.8846912227880697`*^9}},
+ CellLabel->"In[42]:=",ExpressionUUID->"fe6a6d03-8b9a-4d30-9d0f-09497bed842b"],
+
+Cell[BoxData[
+ RowBox[{
+  RowBox[{"0.27286889266851116`", " ", "\[Theta]"}], "-", 
+  RowBox[{"0.2940020307476348`", " ", 
+   SuperscriptBox["\[Theta]", "3"]}], "+", 
+  RowBox[{"0.04000495532616536`", " ", 
+   SuperscriptBox["\[Theta]", "5"]}], "+", 
+  RowBox[{"0.006119426651084101`", " ", 
+   SuperscriptBox["\[Theta]", "7"]}], "+", 
+  RowBox[{"0.0014992967087917594`", " ", 
+   SuperscriptBox["\[Theta]", "9"]}], "+", 
+  RowBox[{"0.0016724305244486051`", " ", 
+   SuperscriptBox["\[Theta]", "11"]}]}]], "Output",
+ CellChangeTimes->{{3.8846911187487593`*^9, 3.884691129095536*^9}, 
+   3.884691222976877*^9},
+ CellLabel->"Out[42]=",ExpressionUUID->"142568cc-1192-404c-88bb-d7eaa3fbf94e"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+ RowBox[{
+  RowBox[{
+   RowBox[{"hCas", "[", "\[Theta]", "]"}], "h0"}], "//", 
+  "Simplify"}]], "Input",
+ CellChangeTimes->{{3.884691136810926*^9, 3.884691148498569*^9}, {
+  3.8846912155321903`*^9, 3.884691226484233*^9}},
+ CellLabel->"In[43]:=",ExpressionUUID->"15c6347f-e08e-4478-999e-a055b83a9e4f"],
+
+Cell[BoxData[
+ RowBox[{
+  RowBox[{"0.27286889266851116`", " ", "\[Theta]"}], "-", 
+  RowBox[{"0.2940020307476348`", " ", 
+   SuperscriptBox["\[Theta]", "3"]}], "+", 
+  RowBox[{"0.04000495532616536`", " ", 
+   SuperscriptBox["\[Theta]", "5"]}], "+", 
+  RowBox[{"0.006119426651084101`", " ", 
+   SuperscriptBox["\[Theta]", "7"]}], "+", 
+  RowBox[{"0.0014992967087917599`", " ", 
+   SuperscriptBox["\[Theta]", "9"]}], "+", 
+  RowBox[{"0.0016724305244486053`", " ", 
+   SuperscriptBox["\[Theta]", "11"]}]}]], "Output",
+ CellChangeTimes->{{3.8846912189151897`*^9, 3.8846912266590776`*^9}},
+ CellLabel->"Out[43]=",ExpressionUUID->"99fdb326-dc62-43d2-a350-3a194a15d647"]
+}, Open  ]],
 
 Cell[CellGroupData[{
 
@@ -95,31 +236,37 @@ Cell[BoxData[{
    3.84380599279574*^9, 3.8438060263082247`*^9}, {3.843806057721389*^9, 
    3.8438061646410103`*^9}, {3.8438062079468822`*^9, 3.843806212028804*^9}, {
    3.843806310535351*^9, 3.843806319271604*^9}, {3.843811727076584*^9, 
-   3.843811750155121*^9}},
- CellLabel->"In[30]:=",ExpressionUUID->"04b6d9c8-662b-43f2-955f-87d9d6c7b622"],
+   3.843811750155121*^9}, {3.876210266924039*^9, 3.876210267043165*^9}, {
+   3.8762103842382107`*^9, 3.876210426032144*^9}, {3.876210485023739*^9, 
+   3.876210485447196*^9}},
+ CellLabel->"In[37]:=",ExpressionUUID->"04b6d9c8-662b-43f2-955f-87d9d6c7b622"],
 
 Cell[BoxData["0.27286889266851116`"], "Output",
  CellChangeTimes->{{3.8438061496682377`*^9, 3.843806165132856*^9}, 
    3.8438062124458447`*^9, {3.8438063135040703`*^9, 3.843806319765049*^9}, 
    3.8438114231817093`*^9, {3.843811741230702*^9, 3.8438117504040937`*^9}, 
    3.8438127470406*^9, 3.8438833063506927`*^9, 3.8439767083823442`*^9, {
-   3.844135505327883*^9, 3.84413553304257*^9}, 3.844143298259859*^9},
- CellLabel->"Out[32]=",ExpressionUUID->"747fa724-da10-4534-a92d-12fdfb1ff525"],
+   3.844135505327883*^9, 3.84413553304257*^9}, 3.844143298259859*^9, 
+   3.876208430215836*^9, 3.876210287689782*^9, {3.876210430601141*^9, 
+   3.876210486050713*^9}, 3.884691139112475*^9, 3.8846912032549753`*^9},
+ CellLabel->"Out[39]=",ExpressionUUID->"2797a478-d270-4869-ba7f-428591635fb9"],
 
 Cell[BoxData["1.007010684856479`"], "Output",
  CellChangeTimes->{{3.8438061496682377`*^9, 3.843806165132856*^9}, 
    3.8438062124458447`*^9, {3.8438063135040703`*^9, 3.843806319765049*^9}, 
    3.8438114231817093`*^9, {3.843811741230702*^9, 3.8438117504040937`*^9}, 
    3.8438127470406*^9, 3.8438833063506927`*^9, 3.8439767083823442`*^9, {
-   3.844135505327883*^9, 3.84413553304257*^9}, 3.844143298261297*^9},
- CellLabel->"Out[33]=",ExpressionUUID->"c1b4bacd-7ed5-4a85-9925-55ead180e80f"]
+   3.844135505327883*^9, 3.84413553304257*^9}, 3.844143298259859*^9, 
+   3.876208430215836*^9, 3.876210287689782*^9, {3.876210430601141*^9, 
+   3.876210486050713*^9}, 3.884691139112475*^9, 3.884691203256315*^9},
+ CellLabel->"Out[40]=",ExpressionUUID->"567108b3-a280-49d8-9f91-9b9feb4211c2"]
 }, Open  ]],
 
 Cell[CellGroupData[{
 
 Cell[BoxData[{
  RowBox[{
-  RowBox[{"totCas", "=", "13"}], ";"}], "\[IndentingNewLine]", 
+  RowBox[{"totCas", "=", "4"}], ";"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"ds", "=", 
    RowBox[{
@@ -131,8 +278,8 @@ Cell[BoxData[{
        RowBox[{"h0", " ", 
         RowBox[{"hCas", "[", "\[Theta]", "]"}]}], 
        SuperscriptBox[
-        RowBox[{"(", 
-         RowBox[{"t", "[", "\[Theta]", "]"}], ")"}], 
+        RowBox[{"RealAbs", "[", 
+         RowBox[{"t", "[", "\[Theta]", "]"}], "]"}], 
         RowBox[{"\[CapitalDelta]", "[", "2", "]"}]]]}], "]"}], "[", 
     "\[Theta]0Cas", "]"}]}], ";"}], "\[IndentingNewLine]", 
  RowBox[{
@@ -201,21 +348,18 @@ Cell[BoxData[{
    3.8438064612671556`*^9, 3.843806483348926*^9}, {3.843811490926489*^9, 
    3.843811490990367*^9}, {3.843811762123721*^9, 3.843811792156164*^9}, {
    3.84381190693526*^9, 3.843811907005992*^9}, {3.843812765889268*^9, 
-   3.8438127660003633`*^9}, {3.8438833187875977`*^9, 3.843883318858028*^9}},
+   3.8438127660003633`*^9}, {3.8438833187875977`*^9, 3.843883318858028*^9}, 
+   3.876208442731238*^9, {3.876209028372431*^9, 3.876209031292045*^9}},
  CellLabel->
-  "In[360]:=",ExpressionUUID->"9e250256-8cd3-402a-bc21-3ee002bb1f82"],
+  "In[188]:=",ExpressionUUID->"9e250256-8cd3-402a-bc21-3ee002bb1f82"],
 
 Cell[BoxData[
  RowBox[{"{", 
   RowBox[{
    RowBox[{"-", "1.3595232308079763`"}], ",", 
    RowBox[{"-", "0.04655549169678349`"}], ",", "0.03237211997439275`", ",", 
-   RowBox[{"-", "0.04559719478179554`"}], ",", "0.0896082623558699`", ",", 
-   RowBox[{"-", "0.2143894700632707`"}], ",", "0.584185898399437`", ",", 
-   RowBox[{"-", "1.7457422693875277`"}], ",", "5.58836650656107`", ",", 
-   RowBox[{"-", "18.86651788362811`"}], ",", "66.44778572365097`", ",", 
-   RowBox[{"-", "242.23334827652096`"}], ",", "908.6585949684338`", ",", 
-   RowBox[{"-", "3491.621379722553`"}]}], "}"}]], "Output",
+   RowBox[{"-", "0.04559719478179543`"}], ",", "0.08960826235586972`"}], 
+  "}"}]], "Output",
  CellChangeTimes->{
   3.83698713111462*^9, 3.836987321633688*^9, 3.8369882151463833`*^9, 
    3.83698833312978*^9, 3.837070043781701*^9, 3.837070402550199*^9, 
@@ -226,9 +370,471 @@ Cell[BoxData[
    3.843803045143682*^9, {3.843806458261032*^9, 3.843806484750998*^9}, {
    3.8438117627002707`*^9, 3.843811778845757*^9}, 3.843811897566523*^9, 
    3.843812875019771*^9, {3.8438833095937767`*^9, 3.843883327729155*^9}, 
-   3.844157549053947*^9},
+   3.844157549053947*^9, {3.8762084391800756`*^9, 3.876208443033372*^9}, 
+   3.876209031472166*^9, 3.8762104487066717`*^9, 3.876210487363555*^9},
+ CellLabel->
+  "Out[191]=",ExpressionUUID->"cd2ce696-c253-4ff8-b993-64c7e05c4dd9"]
+}, Open  ]],
+
+Cell[BoxData[{
+ RowBox[{
+  RowBox[{"totCas", "=", "4"}], ";"}], "\[IndentingNewLine]", 
+ RowBox[{
+  RowBox[{"ds", "=", 
+   RowBox[{
+    RowBox[{
+     RowBox[{"invDerivativeList", "[", 
+      RowBox[{"totCas", "+", "1"}], "]"}], "[", 
+     RowBox[{"\[Theta]", "\[Function]", 
+      FractionBox[
+       RowBox[{"h0", " ", 
+        RowBox[{"hCas", "[", "\[Theta]", "]"}]}], 
+       SuperscriptBox[
+        RowBox[{"RealAbs", "[", 
+         RowBox[{"t", "[", "\[Theta]", "]"}], "]"}], 
+        RowBox[{"\[CapitalDelta]", "[", "2", "]"}]]]}], "]"}], "[", 
+    "\[Theta]", "]"}]}], ";"}], "\[IndentingNewLine]", 
+ RowBox[{
+  RowBox[{"dMs", "=", 
+   RowBox[{"m0", " ", 
+    RowBox[{"NestList", "[", 
+     RowBox[{
+      RowBox[{
+       RowBox[{"D", "[", 
+        RowBox[{"#", ",", "\[Theta]"}], "]"}], "&"}], ",", 
+      FractionBox["\[Theta]", 
+       SuperscriptBox[
+        RowBox[{"RealAbs", "[", 
+         RowBox[{
+          SuperscriptBox["\[Theta]", "2"], "-", "1"}], "]"}], 
+        RowBox[{"1", "/", "8"}]]], ",", "totCas"}], "]"}]}]}], 
+  ";"}], "\[IndentingNewLine]", 
+ RowBox[{
+  RowBox[{"datCas", "=", 
+   RowBox[{"Table", "[", 
+    RowBox[{
+     RowBox[{
+      RowBox[{"-", 
+       RowBox[{"Sum", "[", 
+        RowBox[{
+         RowBox[{
+          RowBox[{"dMs", "[", 
+           RowBox[{"[", 
+            RowBox[{"k", "+", "1"}], "]"}], "]"}], 
+          RowBox[{"BellY", "[", 
+           RowBox[{"n", ",", "k", ",", 
+            RowBox[{"ds", "[", 
+             RowBox[{"[", 
+              RowBox[{";;", 
+               RowBox[{"n", "-", "k", "+", "1"}]}], "]"}], "]"}]}], "]"}]}], 
+         ",", 
+         RowBox[{"{", 
+          RowBox[{"k", ",", "0", ",", "n"}], "}"}]}], "]"}]}], "/", 
+      RowBox[{
+       RowBox[{"(", 
+        RowBox[{"n", "+", "1"}], ")"}], "!"}]}], ",", 
+     RowBox[{"{", 
+      RowBox[{"n", ",", "0", ",", "totCas"}], "}"}]}], "]"}]}], 
+  ";"}]}], "Input",
+ CellChangeTimes->{{3.836976284774213*^9, 3.8369763107653933`*^9}, {
+   3.8369763440621157`*^9, 3.836976643043976*^9}, {3.8369766791263657`*^9, 
+   3.836976702020906*^9}, {3.8369833396004457`*^9, 3.8369833834965153`*^9}, {
+   3.836983426425185*^9, 3.836983501138253*^9}, {3.8369836692218637`*^9, 
+   3.836983678733614*^9}, {3.836983808464398*^9, 3.8369838086001043`*^9}, {
+   3.8369839370347757`*^9, 3.836983937122395*^9}, {3.836983969035384*^9, 
+   3.836983979299178*^9}, {3.836984097325984*^9, 3.8369840974613132`*^9}, {
+   3.836984339450017*^9, 3.8369843561781673`*^9}, {3.836984934917253*^9, 
+   3.836984935916635*^9}, {3.8369856092092857`*^9, 3.836985626001013*^9}, {
+   3.836985701651613*^9, 3.836985833037245*^9}, {3.836985978736498*^9, 
+   3.836985996927923*^9}, {3.836986994292021*^9, 3.836986994346624*^9}, {
+   3.836987100533333*^9, 3.836987147389551*^9}, {3.836988211802526*^9, 
+   3.8369882124812613`*^9}, {3.83698827789044*^9, 3.836988278050117*^9}, {
+   3.836991140640276*^9, 3.836991156327482*^9}, {3.837070011059902*^9, 
+   3.8370700129717007`*^9}, {3.837070292889572*^9, 3.8370702929846687`*^9}, {
+   3.8370709285346727`*^9, 3.8370709288925543`*^9}, {3.837072484698299*^9, 
+   3.837072484992663*^9}, {3.843802126725583*^9, 3.843802127123851*^9}, {
+   3.8438022381502447`*^9, 3.843802268518549*^9}, {3.8438024340742607`*^9, 
+   3.843802515907176*^9}, {3.84380268121432*^9, 3.843802684734229*^9}, 
+   3.843802806265029*^9, {3.843802947563074*^9, 3.843802949058984*^9}, 
+   3.843803003228539*^9, {3.843803040877533*^9, 3.843803044636691*^9}, {
+   3.8438064612671556`*^9, 3.843806483348926*^9}, {3.843811490926489*^9, 
+   3.843811490990367*^9}, {3.843811762123721*^9, 3.843811792156164*^9}, {
+   3.84381190693526*^9, 3.843811907005992*^9}, {3.843812765889268*^9, 
+   3.8438127660003633`*^9}, {3.8438833187875977`*^9, 3.843883318858028*^9}, 
+   3.876208442731238*^9, {3.876208483378458*^9, 3.8762084890345097`*^9}, {
+   3.8762090363083553`*^9, 3.876209068516684*^9}},
+ CellLabel->
+  "In[192]:=",ExpressionUUID->"4b92dc3c-7bc7-4b5c-a5c3-e2aacf85ed52"],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+ RowBox[{"pCasM", "=", 
+  RowBox[{"ParametricPlot", "[", 
+   RowBox[{
+    RowBox[{"{", 
+     RowBox[{
+      FractionBox[
+       RowBox[{"h0", "  ", 
+        RowBox[{"hCas", "[", "\[Theta]", "]"}]}], 
+       SuperscriptBox[
+        RowBox[{"Abs", "[", 
+         RowBox[{"(", 
+          RowBox[{
+           SuperscriptBox["\[Theta]", "2"], "-", "1"}], ")"}], "]"}], 
+        RowBox[{"15", "/", "8"}]]], ",", 
+      RowBox[{"-", " ", 
+       RowBox[{"datCas", "[", 
+        RowBox[{"[", "1", "]"}], "]"}]}]}], "}"}], ",", 
+    RowBox[{"{", 
+     RowBox[{"\[Theta]", ",", "1.001", ",", "\[Theta]0Cas"}], "}"}], ",", 
+    RowBox[{"AspectRatio", "->", "1"}], ",", 
+    RowBox[{"PlotRange", "->", 
+     RowBox[{"{", 
+      RowBox[{
+       RowBox[{"{", 
+        RowBox[{"0", ",", "3.3"}], "}"}], ",", "Automatic"}], "}"}]}], ",", 
+    RowBox[{"PlotStyle", "->", 
+     RowBox[{
+      RowBox[{"ColorData", "[", "97", "]"}], "[", "3", "]"}]}], ",", 
+    RowBox[{"PlotPoints", "->", "40"}]}], "]"}]}]], "Input",
+ CellChangeTimes->{{3.876208490482359*^9, 3.876208512090119*^9}, {
+   3.8762085423630447`*^9, 3.876208589275639*^9}, 3.8762086279652348`*^9, {
+   3.8762086597653513`*^9, 3.876208673349307*^9}, {3.8762087044624357`*^9, 
+   3.8762087621668367`*^9}, {3.8762088101927567`*^9, 3.87620901371556*^9}, {
+   3.876209044101116*^9, 3.876209165118668*^9}, {3.876209218072238*^9, 
+   3.876209218151585*^9}, {3.8762103169644537`*^9, 3.876210332444528*^9}, 
+   3.8762104042786694`*^9, 3.876210437783276*^9},
+ CellLabel->
+  "In[196]:=",ExpressionUUID->"9c63e4ce-bd4c-4c9e-9b17-80e65259d14e"],
+
+Cell[BoxData[
+ GraphicsBox[{{{}, {}, 
+    TagBox[
+     {RGBColor[0.560181, 0.691569, 0.194885], AbsoluteThickness[1.6], Opacity[
+      1.], FaceForm[Opacity[0.3]], LineBox[CompressedData["
+1:eJwtzgk01I8TAPAlIhUKuULIEUXHhkIzzhCFXNkSITnKbUspRMgtUVSuXG0U
+RVIiutBBpPJLpehw9Nevdn33/u//vf+8N2/e5828mVE/EO4aKEwike4J8n91
+qiwsuIKoQ1vPnriZdCHspapJXSioxR6x9WQ6VQhJ2rW3bR1qMCAwZ3tukBDm
+Wq4A54VX0f7z2jY3OyFsH39Xe0S0Cje/6CuI2CKE+6meP6/LV6KByvVCnp6g
+r0PI2R8pRxdmwFkvKSEcWyqTuXJ/GQ5/bQhWFBZCu/4iZiD7Mgon2Roa0UlY
+WmF63uVoKb5WOZUY+YGEPpmtARdcirE0Tj0guZGEQ+dQS/ZIEXZcsxO+VUVC
+ZUd+RWLxeSSWTtE0LpJQdviVve/SQhxqSui4nULC8+SsO6u9zuGcvojH63gS
+5sysaKReL0AnNVV3rUgSylQ4iBuE5+O905MUfx8SFioGCwd9yMPROekiB3cS
+0jiiX0ac89BNpKF4vyMJo8jb3BaV5qBP5bNDS0wF+yb7V+5YkYm5a2Il7ZRI
++FuyX+aa01mU8nWril4u6B86RtzJzkClbU5fHkqQ8FeH9vVprXT03fRtQco8
+H251i4oNXUtFxS9WY8EP+PA7AbKSlFNxbOz22ls3+fD+V19Le34KWieMN2+s
+4oO6hUF8Y/ZpfOTbZ/8ynQ8NLAlW6/0knJPsc3L24IMCRDo2uSYha/B4sqYD
+H35SFGR8ZxJxxtebrLGNDxkZf6LX6SWi7UxT5WVtPnBU382VFCUgs9fkRziT
+B/3p20eMpI+hR//f223VPOBHP5eVSz+KT7Nq57klPHi5oeSItchRjN5YvOhg
+Hg9K37pc+ShOxSztDQGlx3nwZzorNn5bLAYpvptO3c0D/1IbfbvhGAyDqZDn
+9jwgkS7/MjgcgxtmLv82RB7o6zdLvayNRlKi95Ud63jwaeSM607zKMz2Htr6
+WJwHku47Lw//iESLdZTir3wumC2TEqsqikQzan2y0jwXjk6fJj8mItAvZGDF
+wCQXshJ0Y+SvheOmnMZwp0dciF3fyCZ1hmJl944BcioXMncvtxvfGYrFUppN
+dQlcyDGsFgocD8Ek695HBnFcKPsS/ltYMgQpokttIoK4EN0V5hV/+hCKtj/w
+7nLgwkPijUj5QCC2MBcaeslyIc3Ukr83ORC5z/uyXyzlQtP5A/fubQ5Eqq/c
+tLMYF6zI9sfXVwTgRN2Fp4ksDpRtedvRl+GPxhrKnqfGOTAPbeXHL/qhclrL
+4IsmDpDZ1ssOBvhhgV6F0eHrHECZu17b1vthOo+wWVHLgeoa76DKPl9cze+Y
+TLnEgffFZ5VqF/kiefudkYE0Djxx/tStuGcf9jzm1Un4cuDV1O92BfF9qPns
+6hJjCgcM/jQr+bftxfZmt/2hHhy42Tkg/mDlXuyrsI+eduTAkmC9vNd/vFEt
+kv3GdisH8kVthg1HvPB98NYuzxUc4A2rNdGLvDBY5d5TYhkH/j6KvVjv5YX9
+g8eOlC/lwHiq7ILej55Ie6jqskCUA6lVXrQQugfeS2k5KfeXDY/EOmlNdu6I
+GyXpe4fYUMf86Zce7oq5/BlX6SI2WB7SISq3u+K+HmLVsQI2BFuKVJutcsX3
+/pLZ33LYsMz41fjEaxfMpJ2TfZXGBpWkXVUD5i5YG5jc8eYYG9p+yGvF6Drj
+M11+700fNix/xj7rlumIzQq8eh99wf3SlQaFbo4oVLX5q50uGx5evVN5WNUR
+uw+sKzPSYkPNCgmpyVs7MP9ue/JKNcF/Of4X6ycdkDhzjLxMhg1RP/LEnxyw
+x5iUhQcvsFnQ8E7fdarVFn2V3L78+5IFeZQS79AsW1xgJ1K38DkLyiVJG30P
+2GK8ec81lV4WPDN2enpFyhYVi987ufew4H5tfKxVuA2WKeyjTN1hQZbFqfDP
+FtYo/ozmVlPBAt5oo88ufUt8cGGTJFBZEKhGHry4xBK/vGt2uhHDguqAQ4cN
+f1mgR8eulRpRLKh6XOlu0myBG+J63KQOs2Co7KRKrbkF8nVtWsX8WTC7M1Lc
+yA9R+fSHhLBdLBgfsWrfQzNHz3p9p6k1LPjz6uqF84nmONfWYdyvwwLxwTWb
+9TzM8Vu1xFyDFgsW5NSM7hE2xxNNop1H1VkQnLZvqs3HDOl6qlmaiizoPtX9
+6Ka2KfIGAir7xVmQtsHkza9ZE6StWPTa7ScTdl9opVc8N8EC5dVlId+ZcEV0
+fWchzQSzOmtnkyaZ8ER1stcoxAQtLl5a2DLOhKoG4SLKlDH+NZT6z9pRJrx9
+0T19lW6ENwrHxvz6mRCUlnaeSt6MaeZp8kONTDhCff7xjfxm1EFfmbnrTPAR
+LtUMZZMxS80sR5LGhMeNaVEmPWQcVG4jdtUyYbHXIrtd7mQ0nO0LmihjQijl
+gykneRMuPTnyxyWfCX9v9N2RFdqId/3NWzbGMmHkW0n+wIQBRhDc/Q3IhA+x
+eyS+dRng3XjjYE1gQkDzWnmjywaYM8tPKDFnQr/QU/ppDwNcp015nLmVCc4h
+c469/evQuTKhPp4smGePNRd0rcWwJzpd6bpMiDL+/FttQA+3ELouPsuY4HYl
+8WlYlA5OcbrZLZMEDDpfrrxM0cGo2jLmiQkCFkk3Gcxa62B6nFud9VcCYqor
+E7/L6yBEl8mNfCZgN/XfzVad2th+YJ+K0AcCyreONQXLaGNTSaoL9TUBPyyP
+R6sNrMY47RDWp04CJLPXij/P18Bqq9Y9yZcIyFiYqzB6SgN18/te7SklwE6H
+eMk/rIH3/RynN5QI5sepZ/MdNJC6fJ30RDEBnjvj3/aKaGBEajnJ+RwB89rR
+d0MT1HFU9ythc5aAQM/UsJiTqzAzi7W75igBRr7FkfwIFTSbmI9U9iTA0of7
+WUlGAWsitZ6OqxAgHtXUQCHJ4t+Y4cVl0/PQ3Tyj/idDGuuXxF3Pa58Hs7mI
+hUNzi9Fs/xM93fx5sDrURqFEiyH1OynY89A8/L6pKXfr0wKcmfkp+8ZqHrJU
+H6mPrRTCwOGMlJzV8zB5wvMGz4QLdqlqCo6i8/BPjPknzx1MEFO4oOM5xQDK
+iRft0zfooC6kr2b/igEiIFK6svFf+LnpmndDKwMMTMts0XoOrIO7uxwqGNBZ
+XdNCCM3AJuOQgrpMBiiojp+b8PgB6Xd3dvceZcC/jJQMUcNJ4IlttDocxADT
+Ygbp47ovUCClvLHHnQGjE0OPd4V9Alp9hbGDLQMccuuOU2+NQpdBVWSGEQN0
+G/RsvC6NgIz+RZqoLgO8sykl0cZDMOuo4m2jzIDptNbZ4o+v4Bbtbrm5JAMU
+9UNd7gX2gVPHcMsTYQbEFTabHip4DMvF9xzQJuggd+LkJdWPXeC0LEnBapYO
+QTYH5STe3oUMWY2TtK90aDX4hU/Um2GYnfvp/Xs6/BZZr/TLvhbirXVnNwzS
+wdVdRcM0uRiUFI4HremlQ4nY0ngNtWKY1FistVDgxXH388n3iyBxt9n0l2d0
+CJ4ub3jHOA81147lXxI4kLTYU/VIIQBrtY20wEdVK85lB+RD6YeD25lP6HAN
+VQM+fj4Db7eL1A8+EuxjuF6S1DoDS7bQ4m4InKP4yy03OBUe/nPTMVtg7mBA
+Q9Kf03BGfmCZvcDLc071aS1OAr8OLdLDHjrMX02aKSyIg+g7IaHN3XR4IbHl
+xvjqWNjfKPEgT+BNzWOtfXei4YXhT4VwgSs3D4dGfIqAwc97vusLHCL9nSYS
+dhBU019/q35Ih/LY186uJf4gSjE8niLw4oDWr0O9vkD93i7vL7CBprP8uJcH
+nHTxCF4lcNeAHVu20xlErA6o8rvoQB14t89rly3Q5c7/MyZw+aiMVZGDEXDe
+zJTfF5j0/wjPDTpSKvB/ASWdn3s=
+       "]]},
+     Annotation[#, "Charting`Private`Tag$384300#1"]& ]}, {}},
+  AspectRatio->1,
+  Axes->{True, True},
+  AxesLabel->{None, None},
+  AxesOrigin->{0, 1.3595232308079763`},
+  DisplayFunction->Identity,
+  FrameLabel->{{None, None}, {None, None}},
+  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
+  GridLinesStyle->Directive[
+    GrayLevel[0.5, 0.4]],
+  ImagePadding->All,
+  Method->{
+   "DefaultGraphicsInteraction" -> {
+     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
+      "Effects" -> {
+       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
+        "Droplines" -> {
+         "freeformCursorMode" -> True, 
+          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "ScalingFunctions" -> 
+    None},
+  PlotRange->{{0, 3.3}, {1.3595232308079763`, 1.8828734633575026`}},
+  PlotRangeClipping->True,
+  PlotRangePadding->{{0, 0}, {
+     Scaled[0.05], 
+     Scaled[0.05]}},
+  Ticks->{Automatic, Automatic}]], "Output",
+ CellChangeTimes->{
+  3.876209165468998*^9, 3.876209218440271*^9, {3.876210314216745*^9, 
+   3.876210332711735*^9}, {3.8762104336122227`*^9, 3.876210449468999*^9}, 
+   3.876210487907978*^9},
+ CellLabel->
+  "Out[196]=",ExpressionUUID->"79c580a9-d885-45d2-b73d-816f14fe791e"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+ RowBox[{"pCasS", "=", 
+  RowBox[{"ParametricPlot", "[", 
+   RowBox[{
+    RowBox[{"{", 
+     RowBox[{
+      FractionBox[
+       RowBox[{"h0", " ", 
+        RowBox[{"hCas", "[", "\[Theta]", "]"}]}], 
+       SuperscriptBox[
+        RowBox[{"Abs", "[", 
+         RowBox[{"(", 
+          RowBox[{
+           SuperscriptBox["\[Theta]", "2"], "-", "1"}], ")"}], "]"}], 
+        RowBox[{"15", "/", "8"}]]], ",", 
+      RowBox[{"2", " ", 
+       RowBox[{"Abs", "@", 
+        RowBox[{"datCas", "[", 
+         RowBox[{"[", "2", "]"}], "]"}]}]}]}], "}"}], ",", 
+    RowBox[{"{", 
+     RowBox[{"\[Theta]", ",", "1", ",", "\[Theta]0Cas"}], "}"}], ",", 
+    RowBox[{"AspectRatio", "->", "1"}], ",", 
+    RowBox[{"PlotRange", "->", 
+     RowBox[{"{", 
+      RowBox[{
+       RowBox[{"{", 
+        RowBox[{"0", ",", "3.3"}], "}"}], ",", "Automatic"}], "}"}]}], ",", 
+    RowBox[{"PlotStyle", "->", 
+     RowBox[{
+      RowBox[{"ColorData", "[", "97", "]"}], "[", "3", "]"}]}], ",", 
+    RowBox[{"PlotPoints", "->", "40"}]}], "]"}]}]], "Input",
+ CellChangeTimes->{{3.876208490482359*^9, 3.876208512090119*^9}, {
+   3.8762085423630447`*^9, 3.876208589275639*^9}, 3.8762086279652348`*^9, {
+   3.8762086597653513`*^9, 3.876208673349307*^9}, {3.8762087044624357`*^9, 
+   3.8762087621668367`*^9}, {3.8762088101927567`*^9, 
+   3.8762089596587143`*^9}, {3.876209089221918*^9, 3.876209122589672*^9}, {
+   3.876209181039466*^9, 3.876209181767116*^9}, 3.876209220528407*^9, {
+   3.8762103380532513`*^9, 3.8762103385007963`*^9}},
+ CellLabel->
+  "In[198]:=",ExpressionUUID->"9a7cc83f-920d-4c9e-9d42-998a933b225c"],
+
+Cell[BoxData[
+ GraphicsBox[{{{}, {}, 
+    TagBox[
+     {RGBColor[0.560181, 0.691569, 0.194885], AbsoluteThickness[1.6], Opacity[
+      1.], FaceForm[Opacity[0.3]], LineBox[CompressedData["
+1:eJwtlmc4F47bxc2MQigje6/ICkXf+zYyUzaRmb0lRBLZexPJbEmlKN80jKyM
+bNkjJdEgK4T+v+e6nvPmvPm8OOfFua7D5+Bj5ERGQkKSTkpC8n8u/dmtiym0
+UZXzetUVOnFyQqR2zZKWBqpypR1bvGcYSmhMFH9u+ZpPdYZ6QYGtN5dwe7Z/
++UDmQVWvSwpps/IPCKxfb7QqGpGrRvV+3tyjfE7go3HULPTYQlva0qv3c94S
+Xhwb4gujW8Hn11tb7lk2Ez5VafCyUX/HlzriT6kedBDSB3rfydHMoQ/vTX+u
+j70EuRqP2+YTM5ibMB936+kgIVmi7rDhw3F85MTCszgzQihpUKBa5B1GnZof
+5zsdJwkVysvUuWUDGN4QmczP9YkQFjV7QWauF3E88lt75WfCHO/fTzxOH7B+
+lfN8mttXwrE14DTua8ejCW9F3uosECYURw3LJ1tRIGHG5f65H4R215i8LOdm
+XH574Gl1yBKh/X5f/aLhOwxXU+IXc/pNcNe0bWVIakBRDomuPe5VQm+RKK/o
+gTpMO6MZX7C4RhCQla8ddX2NNzw4yN+3bRA2onj3+DhrsVEx1123ZpNwuLfR
+6ikrEa9QfWZffrFNePvGapPe4AUGTPWpx7XuEDSWx2i3X1fjGUnXJ1/m9wjs
+5qaUn42rkE5D5/WPJBJQnUgkzWR8ho4H33lxPyGFzrjIjLbpSvxnrhbnP0wG
+P987TzENPcaIAFbaUWoKsFAItOv9W4Eb/qKtKmqUkPeFjHle/yFqsc8aJUfu
+g2HG/FTR1gf4kWnkb20XFXBj4LCI3X20T5x5XsVFA1UknaZsbPfwT+UVetdA
+WshKENCuXryD+QPaEp0f9wP9vV/i/ENlWG9/p+YDgQ4kdSU5fw6WIjXdt3H7
+J/TwLSbZ8cd8CVrs6fO8ZToI07lr0SoHSlCKv3SLov8gdA5/MJTHYtxZc1Ih
+y2aEMGuidi1FEZ6WPzX6xIoJioQmfvaX3sZZL9nfFMLMILX16LCrRQEKHRXp
+31tlhgSBFxseArfQnzVwLaflEIQsnuSiJctHddt3m235h8EyWkW3aOUmcvQs
+3U3yZwEt0066NfJclIxM73hzjA0ucL+KkObPwY+Nb8K4mdlhf0OCkal+Nn7m
+oTci2WYHVzXnk25RWSj2duSp9ZcjwJXoIpjQmokuFHyCsv0cIDVbXLTIlIm/
+n1wwCG7ihHyh23d6XDOQO1hzTeolF5DGfG1KlErHJqvJI0MVPDBAKRGgdDsN
+uxrsOmxMeUHEeHcfE3Ma8lh4n7jxkRdmm2JoVXlSkZG4s2p0ng9UszTMvgum
+YJQjkatnkg/YQtK/myskY/+/o1/WLvKDld/z0k6jJBy94Bf27js/2BDXb/A/
+TMDVFwSlsySCoKpnp5swEI+M0+VWKzcF4bVo+C4lVTxmF4iHhMsJAcsUmUBz
+fCxOnt25YuQtDCxfbrUsDMcgj3+OVQa9CBR+I8+cl4zBF/kcp54/FQF7Fcqf
+hUtReFGW3zNnWxS4kvyyYyyj8KtcP5VxqRhMnXl9qqQzEuM4pQmLeuLAvMZc
+xVp/Aw2UNC6X3JGA/APjZB/xBoaQ2VPUGR+FkGt+voltEbgvKvjcEwpJiNF6
+9frjbDjOUdHHk3tKQe5jnYS1wHDM05FQc+I/BtnFDDVsDOF41Z5TL2vsGCyP
+Z14zyw7DJ4e2D4k/kIboyKcxwePX8PwC5b/lfTLw7jc99XHRa3hx4l/KLWcZ
+SJQ+9EN58CpmEg7V3xSVhZbg6dF12as4uAnrCwmyUEigIQrmhqCbntkuxy9Z
++H5P4N6CdzBym0h28RDlYIm2fojMLAhPMaV/zzQ+Dp8PDrsrjwWiP0+90Oib
+42BlFLOl7RCIY3l6pNvCCnBEvqJE/VoACo+fleraVQDZ37rMXxkCUKr3rX6o
+uyKobj9Ybrt7+b89XnxIOaIIccMf4P2YP6Z+jOkqeaEEV/2utZcH++NTbxWf
+aqETMFF/sMGAwx871SNoCnJOwPsdmeQo50u4ZZq0shRyEvxWM3in2P2QVyU5
+K1NfBczmGpuVu3zRMwfm9JtVoKk03PZ+uC96FJiujimfgvnlM951Uz5ItfhU
+y0OSAGw/84sYjX1wRSk63es+AepLPpFtdHijX1ZXEfIDXNo8nPan2QtvVy2m
+StAhMBk9jx/T9MLpTxs9XXoIRzVr5MY6PdFhIyl8KAGhRi6C8v2UBx6oLNDm
+pFWF8wdOVzQKuaN5S4N4MJ0ajLkP1z1vdsPpyvK0JV01YLD81Nng5IYnrOcL
+6uLUoOkXc8rXx6547+SvP06U6mBxxyln1MQFudg5/w6SaUDdp7PipOQuGEv7
+cmC/ugYwp8SdH6xyxkblgb95kRrwrrbf9xyrM0ZdFc55TXkapDSsyT1WHFGS
+SBLHQK8JeyEXduoeOeKYywfvHQNNoHJ+2ezh6oheXdJOkKUJeng97tHni7h3
+++/XVg4tkIueuP35mwOmJidwjEhpgxW189NDjx2wturleb/L2uAf3fLb5pID
+WksOXLV9pQ0/H4pVC5A44LJu8PgJLR24MR3s9ErEHoXTY4Q+O+qCpI1yxr9V
+O4yOZCnkeawL4ZPctkcb7bBC6etM5bouLBW4eIjY2GFsq8Wpzlg9EO9Jy5DQ
+sEX+nKb7C5VngI3FOYVp0QbT8ppaG7bPAPXUrXsh6Ta44e15b/G0PiSRnhd+
+9MUa/V3NRBQm9YHmWQlH9a0LGJR6tPYBwznY2nv/q/7MBZTnverEZH0OaP6I
+ZDP8s8Ipb4PhDw/PQeKSGuMFFytcu7LzvkDAABK4qe+snbbEka6y5OVGA3CZ
+Pe1IqmCBXzotfST2GUHh16bfGVvmyNfDznhTygisMkPDyevMUWaUl1fC3AjK
+F/Jv+eqYo0i0hlHzAyMIHzCRp3E3w4GKgeJGfWPgZrsTwCtvhmFhipp2QcZA
+4t8m9WfPFHeDjFYFS4zBvGp0YzDHFA3Oa8vxrhtDbEzx4Gq3CSr1mqxNFJqA
+51FJmQeFJrhvizHiYbsJtJFYz5L6mOBwT/Tz22smkP6oMSmP0QRVDpvPrOma
+wp1RzN9vY4yrdoUUTlumwB+pm/uS0QiXGMevqlmbg2d/VcWNRUNMsDA74JRg
+Dl0NTivtTYZI0tOW/IhoDoM5LH1tQYbIVyqm9JDJAopYL2rpzBmgTRPs2XRY
+AIXrcG1K+Tkk6FCw3kVLsBE0WCBzO4c0bclvGX0todvnkpKS2Dnc31X4uLTI
+EpbxGY5VnMWUmVNqMv8sgVPVJqiiVh/Nm+1FSBqsoGTI+af4uh769BmYX9C2
+BsEETfvpV3oodcLqMoZZAzGfqvZ0hB4qdfxMOvvcGmzKX45xHtRDUUfhiG+8
+NpDgIOBQflwXh118jjLt2IDP638faQq1cV9TiBK7tR2o7Obw3vfSRqcj4m4W
+wXaQIjAt13NKG6lONio/yrYDkrHLocRpLVwdyXua3m0HXR3OzM/EtHAiexuM
+Ve0hUezQi/LB03jdoVX1kbgDxLTS5DI8Oo269XLUctoOYBtOY/U58jRKZrBx
+dDs5wLisz8Vu+dO4fOclo2qJAzR2TSzE3dJAl9jcanb2i//lC1K3uKqO9W43
+ndIPOEJnW+SJ3xbq2PyqJ9tSwhGKQuWMZxXUkfP3w3eKuo7Q8bhN8/2qGjrG
+7/8nGecIOt5U2ZF+ashA+jqUhtIJ8hO8qDuuquLZ9BSvK2TOsBjS/dvMThWF
++vu8rHidoXluVZ7/tCpW+ZVknyM4AzmRddyXQRXTQiT1PUKc4cVsV8PYPcSs
+zbvrPqvOsLJb7pDVAJjMRkVR8c0FGkjIMwbTABM+8HcUUrlC7t8581B7QPN/
+VmFlwq6Q8iqEbI8cMGlHoWza0RUszwwpPtQnoLTk3Wv8s66g2eV4pGZFBalu
+xgX1TLuB7nH1zB/vVXBkyi6y958bCHOL9gYXqeDZjC+aszzu0GloMJmlr4Jm
+u1XzinbuIDpHs1D9WBmDSR3kKz65Q44AsW4w9CT6JOvwV857gDRLpKyooRLq
+zCmxsJB4g1fFJ9MFSSV8qVLYRybkDeeOxoSv0iqhprrBSzJdb8if8O+ablHE
+g5uGm8pZ3mDCE7xJBoro0BAykinmAymNYuV9pxSw/q3et0UJX0g2MJ7o4FZA
+2g6Vui2CLyzsj3hFSqKA5geW9jMZ+cLLwKFloabjeIPRgMr2ii+wXyrbp3jm
+ON5TsNVzb/EFtH2ivekqj+yGm0Ea9n6wmy3rE6Uvj8p2ASsLl/2geWSN4rSs
+PIrzyr2/GecHnWFvvOx35PCZTFgY01M/UPouwrGVIYeyhD9zAf/84AgPuZBS
+uywufS3j+1t0CUofhjyYMZNBSbaOzNV5f1gQ6p26qiqDJmXz27Y7/hDaZn9O
+5agMcvw9Lzd08DKQPa6JlyCTQbbd40xjJy7DIFOW6qtKaZzw6bXWT7oMoRnm
+FUaM0rhXOz09IhsALSHDnQo7xzCJL6C8TysAgsLkx6g/HsNuiaqjAxcCoMzc
+O8Yz4RgahVpYLscEwA5N+wf+DSmcWdgVyxz/j//o7jc0Iom/1Acsg2MCQVJt
+L1y6RhJ7Xmwz7d4KhD8xIbNNmZKo94TSO/5ZIIQz7atNPyeJV8/IUb2ZCITA
+nqOdb9qP4noqa1GcXBB4XLkl+rNNAjPtqGTqZoOAttZC4+M3MVyPcHFJ1woG
+JUufsnPdYtg7t0leeiEYStSY235Vi2GT19BerV8w8PIyqd6+LobVaif3kxQE
+Q/2Lky7VbGLI13je/styMLy9nlsXbyiKY2zf2OkKQqCV3f6QyIQwsui7pw7/
+uQo1Jt+1xluEkdXfZr2PLhTK1K9HPa4Uxtr4pZsDAqFwyo1boyJSGHdOr2X9
+OBcKWZ1SwvskhVHrw5pGwINQqJs9eYwuWgjZXJrKZyyvwZdz2qPMOoLI5y7X
+ytYSBmyjiYsZ8oLIc+Cn+fxkGHjvzDtL8QriTfEl7bqNMJieUOiv/SOArGP7
+NaJFr4NY/cs5zgcC+NU2/FpQ0nVw7rhGuUAngJObJH/XxcPhXZTqktkPPjxs
+sN835k84CLsw6lON86FcEVXSJYoIyE885d7dzod1sUdbHBgj4NY+cs3i+3w4
+U1qVdE4iAhxO/Dpz1pEP054z2FvaRoC3U/uvqk+8GJo0VSr/PgIuyLXcMhzm
+Qc9c9Zo/BTdA5sL2+Y7XPKiWldAeUX4D+ioNFS8W8+CjumsXGWpuwAyNTOKa
+Gw/uUFaKHe+5AXsON86m7XHjyVm51YekkaA0/5AvUpIbk8NoOMAlEiSvDHlx
+l3AihWP1+33Ho+AoLoapxXNib2DbFIlaFHCXtUpH+HGitdIB772zUeD17cq4
+vhonPmKN3CB3iwLFp0mcvHMceI/tSZnS7ShIuHSvIvYYB16+T2TW3hcN/Pr0
+Jzw+sqNPxsUrO2PR4Ji7mUf5jh1NMtYv8H2LBr9qfdsXj9nxR39tos76f3zE
+oVKTaHZk8LVpLWOIAeO01amU4+x4vqT57nWNGODyL0kczWPDv1MDrwafxEC5
+ZGLvfj9WjId/Vb+jY+G/pxHKZcuKJ8V6RK9nxgL5SKyKqj4rKk4oTTKUxMIr
+T7rcOjFWvFrmTaL8JhakGDhFUz6zYNfFa8IvV2KhYkXxmaclC1Y3mZbq28VB
+tW97YIrJYfRuNPw2eyoepoweoK36YQxUibWk0o+HbO4z/Gqyh/H8F/f0Yxfi
+4WjGq/ZTBw/jmPV1vcSQeBjs092h+XAIic+IX3yI8fDFVZPhk84hPGP+Qven
+TAKID+59XjzHjHyVI19fiCWC4W5Oeh8yYyF5P7vRiUSgfKgW1ybDjENig+Qr
+2olA3+NK8omZGVNCtP6quCZCBZPgctsIE5amqnr+u5cI7+fyh02cmHB979HW
+caEkqDtRZbYQx4i1bv6zzoLJwFyzEHAllBE/8re/dZJLBnnF6HkmX0aMj6xr
+d1FLBun7Pklu5oyYIXzyUoBdMpCw0HoUCDMidar5Tk1BMtydfCZG13YQb/9N
+fNR/OAUu2HnONDEcRJLdTBV32lTYbspvvUNxEPsyGxK72VOBZKDxB/cqA7oN
+apEqiKWCrIPQAYo2BvwSIWBDr50K8qd4HF74MKDzHfz4JToV6j3pwkua6fGw
+0bLKIdI0GExqUqEPpUO1vLfSUWFpEGbxS8XLmQ699j8WYIpNg+TpVIUZAzoc
+DgxdLElNA8njH2pIhenQU2ttpbU4Ddazd8yk+w7gM8WxBeGmNOCPdZ4TkTiA
+w5sEDQJ1Onzqowya+0GLVxIFTdoz0yGTo3Qxd5QW1VieX0wvSIf1VpdbVq20
+SLxKqWp1Nx2oO5cdWYppUSdI0nGtJh1o7ar9zE1oMUffSEtlLB2WqHPzN+pp
+8Lkguzk/Xwasm7l/EimjxvRb7gaaTzMgPcdDSyqDGk8Wv0rlqs0ATd1tFvUI
+aryxdCFmozEDuOIfV+faUuMK17ebjwcyoNArXeo7JzVKjWmflf6TAb9nv7dJ
+5lPhXtHwQCwhE/yP7mW8KtyHqmFe7GHdmXByiHXsTxcF0gn1i4WtZsGW0peb
+LA0UyMXbSsm/kwXS2z9WCdUUeORJjFkbRTZ0hk6df5RHgW9jv5qysGRD4nr3
+3KYzBbJ2+b1vUcoG9s92QobkFCi6z8jWPSwbnB5dvN2oRo4We0Hu7bQ5IMcw
+vjYxRoqUl3++HRbJBZ+D9Y5KPaTI1rvWmi+TC5XCXa0FTaQ4VzAUaqucCzRe
+8oTwR6TokHn918+zuRDgabrvWRgpelH9/cIdmAurwXNhVYKk2NY/vt3bnAtP
+A4uV/waSYBXP2SeXL94ErbLbet+c9uAKh4z3s7t5UGko/mvWeA/6Gk4OSlfm
+wQlblaEF1T3wOC1x/9nLPJhQstHk5N6Duc+qmcTOPGC4cCiZaXgXCgI+hA/+
+zgPXA0NZynq7sFJcmW1EyAf6KOn9J1R2wG1bs8NmJB+2Yi13ugnb0HWDM/np
+wQKQ7Q67ryW1DcojCRTWHAUwOCyx2cm1DZoKp3n2CxfA7/4Auz87W7C8RUbw
+VC4AuooSiZHXW0CT48yp5lwAhzjo3r9U3oIsa3JllTcFwGmUL/dMfROq7NJX
+XN1uw+PB1DPBlzbAv+fBn1vvC4FzRrWl13EDLm999AodKAS5s7nqMuYbQDNK
+7WMzVQjn7Q/I0alsQOTQIqPIWiFoji/tl963AW1vuoc6eIqg/zX9cY6CdTBM
+sFW8EFQEjh6FF0271yCF+kZNP3MxRNBvkc7oroJjGRuVU2sxpB2z9KDBVfin
+1uBwt7MYUgJSlE8cX4WizXjnr73FML3qulvBswrkEWec3ceL4dChOzcm11ag
+w+XA6ZDlYvg6KHnds3gFtn18/WuPlECnJefVIzu/YZL5p9wLnxIQfj6g39K8
+DHfiw3SLOEqB/Vtp5a1XyzAwaUsSxlcKpJlnzK48XYbJsa9B1iKlsPekx0aj
+YBmguqSCW64UvpzeHBK+vAzTyb8EHuqWgvdsvQKj0DI4VdLozAWXgmt+23IN
+yxJovZj89X2kFASVYp+X9v0Aje2qHYW8MpBjUdFUJl0At632YzLUdyF0iV/g
+zOs5WCKgf3vwPVD+JOHiUDQLjwvkfkcdfgBcZUt1/MwzMC8T1figqBx0wn/W
+PdcehwoL2RfuphWwn+xw+TLzMCjSijyqoHsM9r817BhvDkBVvZ9MyLMnsDRz
+P0ZSuRfeUtxs0VB9Cg5sbQs6kx3wRmt++CP3Mzhr8qkvcKkFCssthI/tPINl
+9XAHcYFGsOi49XJ+uQoOStp3dyi/Ag5fiX8uo9Ug8TNpdNyxCqguVQjPjz4H
+Ro+bsTu37gOnh6VTxIcX8Olrxw+ZvlyIJ0qYsDTUwGO/oIevPXPhyQS9p1Nz
+DRR9N9b1o86FBkGK7ufva8BxsmDMFnPgcyq1r0lfDeyuSn/VqsqCLtJnErmz
+NTCQ6sImU5gOVrVqV8UpiXBl17SzQiIWIu8k6MTpEkFmyMLM2SsGBu7MZs+c
+JULjWuZl/8poODzlKnHCmAg1rz50RMlHQV/qZ5nvVkRo2GaNuq4WAZ50nRVG
+XkSwavml6TQeCJUKR/WPpREh/bxeOZXbf7+3mlksKYsIdvQPtHT/+IPQKFF/
+4SYRTAWWySsP+8HUO3a/OyVEiHcsFJ+odgatkms83NVEODl/cJZ59CJk/o7E
+a0QirL58rh37zw4GT9I/nHhNBEdZ08RTxWag6ZRhVtBMBPWfLVaRawaQyRWW
+t/OeCJ3qF+8E5WgC7bU+QesP//XpfkWdm64AGz7Rm2/7iEDy/0qczTjE85EI
+/wNkCn8r
+       "]]},
+     Annotation[#, "Charting`Private`Tag$384371#1"]& ]}, {}},
+  AspectRatio->1,
+  Axes->{True, True},
+  AxesLabel->{None, None},
+  AxesOrigin->{0, 0},
+  DisplayFunction->Identity,
+  FrameLabel->{{None, None}, {None, None}},
+  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
+  GridLinesStyle->Directive[
+    GrayLevel[0.5, 0.4]],
+  ImagePadding->All,
+  Method->{
+   "DefaultGraphicsInteraction" -> {
+     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
+      "Effects" -> {
+       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
+        "Droplines" -> {
+         "freeformCursorMode" -> True, 
+          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "ScalingFunctions" -> 
+    None},
+  PlotRange->{{0, 3.3}, {0., 0.09311098339356698}},
+  PlotRangeClipping->True,
+  PlotRangePadding->{{0, 0}, {
+     Scaled[0.05], 
+     Scaled[0.05]}},
+  Ticks->{Automatic, Automatic}]], "Output",
+ CellChangeTimes->{{3.8762085510207253`*^9, 3.87620858950519*^9}, 
+   3.876208628576564*^9, 3.876208673611944*^9, {3.876208726286438*^9, 
+   3.876208762616218*^9}, {3.876208810719953*^9, 3.8762089672639933`*^9}, {
+   3.876209089578055*^9, 3.876209122850472*^9}, 3.876209182088485*^9, 
+   3.876209221479105*^9, {3.876210314516953*^9, 3.8762103388902884`*^9}, 
+   3.876210498192659*^9},
  CellLabel->
-  "Out[363]=",ExpressionUUID->"1fe1a989-6f4d-41a0-9680-c055dfc0798b"]
+  "Out[198]=",ExpressionUUID->"72455824-b956-48f2-805c-31dcb93fc169"]
 }, Open  ]],
 
 Cell[CellGroupData[{
@@ -2385,7 +2991,6 @@ rW6+hjMZrUR2qv0eCs2GZ6hQiiKw8WrrF2rHZuOttEddz8TAmZwN3+XLaXCC
 ivJgcNdeDK+lWgEs+ywogw0/GSLoYewooUER0sKSJJsOpnsqjFCSX0WUBb+n
 YMujfA5G2vROKPVvL8Hzptiqlfxf6YGuwLEX3qpPiHXOIvUa7ZWbR3yTfQ3D
 8pEYURJeyFzHjnlD/Ac6Hgav
-  
   "],ExpressionUUID->"f81ed9c7-d940-4314-963c-c36ea0215ef3"],
 
 Cell[BoxData[
@@ -3721,8 +4326,8 @@ Cell[BoxData[
              0.16666666666666666`, 3.584705786922063}, 
              Offset[{3, 0}, {0.16666666666666666`, 3.584705786922063}]}, {{
              0.16666666666666666`, 3.584705786922063}, 
-             Offset[{-3, 0}, {0.16666666666666666`, 
-              3.584705786922063}]}}], {{{1., 0.}, {0., 1.}}, {0., 0.}}]},
+             Offset[{-3, 0}, {0.16666666666666666`, 3.584705786922063}]}}], {{
+           {1., 0.}, {0., 1.}}, {0., 0.}}]},
         Antialiasing->False]}, 
       {RGBColor[0.368417, 0.506779, 0.709798], 
        StyleBox[{
@@ -4411,10 +5016,10 @@ Cell[BoxData[
  CellLabel->"Out[91]=",ExpressionUUID->"b48c9c3b-23cc-4864-8ea2-6e98855eb77d"]
 }, Open  ]]
 },
-WindowSize->{635.25, 528.},
-WindowMargins->{{2.25, Automatic}, {2.25, Automatic}},
+WindowSize->{1426.5, 740.25},
+WindowMargins->{{0, Automatic}, {Automatic, 0}},
 TaggingRules->{"TryRealOnly" -> False},
-FrontEndVersion->"12.3 for Linux x86 (64-bit) (July 9, 2021)",
+FrontEndVersion->"13.2 for Linux x86 (64-bit) (December 7, 2022)",
 StyleDefinitions->"Default.nb",
 ExpressionUUID->"d257fbda-2be1-4221-9ac8-5153716cfc00"
 ]
@@ -4429,186 +5034,210 @@ CellTagsIndex->{}
 *)
 (*NotebookFileOutline
 Notebook[{
-Cell[558, 20, 596, 10, 24, "Input",ExpressionUUID->"94557b52-cc86-4010-b22d-6b74ac4aad03"],
 Cell[CellGroupData[{
-Cell[1179, 34, 2287, 64, 112, "Input",ExpressionUUID->"04b6d9c8-662b-43f2-955f-87d9d6c7b622"],
-Cell[3469, 100, 489, 6, 25, "Output",ExpressionUUID->"747fa724-da10-4534-a92d-12fdfb1ff525"],
-Cell[3961, 108, 487, 6, 25, "Output",ExpressionUUID->"c1b4bacd-7ed5-4a85-9925-55ead180e80f"]
+Cell[580, 22, 596, 10, 24, "Input",ExpressionUUID->"94557b52-cc86-4010-b22d-6b74ac4aad03"],
+Cell[1179, 34, 689, 14, 20, "Message",ExpressionUUID->"416ee172-1ad9-434e-b73c-e99e451d8e53"],
+Cell[1871, 50, 672, 14, 20, "Message",ExpressionUUID->"ec136869-a624-43fb-87cd-ab7ce0e4f6b5"],
+Cell[2546, 66, 684, 14, 20, "Message",ExpressionUUID->"c71f26b5-fddb-4e00-990e-08b74adae32a"],
+Cell[3233, 82, 684, 14, 20, "Message",ExpressionUUID->"2eb6cbef-108c-4bdb-b0d3-f595ef80164f"],
+Cell[3920, 98, 674, 14, 20, "Message",ExpressionUUID->"2d64117a-f0ef-4d22-873c-490eeb5b19e8"]
+}, Open  ]],
+Cell[CellGroupData[{
+Cell[4631, 117, 362, 8, 22, "Input",ExpressionUUID->"fe6a6d03-8b9a-4d30-9d0f-09497bed842b"],
+Cell[4996, 127, 687, 15, 26, "Output",ExpressionUUID->"142568cc-1192-404c-88bb-d7eaa3fbf94e"]
+}, Open  ]],
+Cell[CellGroupData[{
+Cell[5720, 147, 314, 7, 23, "Input",ExpressionUUID->"15c6347f-e08e-4478-999e-a055b83a9e4f"],
+Cell[6037, 156, 663, 14, 26, "Output",ExpressionUUID->"99fdb326-dc62-43d2-a350-3a194a15d647"]
+}, Open  ]],
+Cell[CellGroupData[{
+Cell[6737, 175, 2435, 66, 112, "Input",ExpressionUUID->"04b6d9c8-662b-43f2-955f-87d9d6c7b622"],
+Cell[9175, 243, 633, 8, 25, "Output",ExpressionUUID->"2797a478-d270-4869-ba7f-428591635fb9"],
+Cell[9811, 253, 629, 8, 25, "Output",ExpressionUUID->"567108b3-a280-49d8-9f91-9b9feb4211c2"]
+}, Open  ]],
+Cell[CellGroupData[{
+Cell[10477, 266, 3932, 87, 115, "Input",ExpressionUUID->"9e250256-8cd3-402a-bc21-3ee002bb1f82"],
+Cell[14412, 355, 1153, 20, 25, "Output",ExpressionUUID->"cd2ce696-c253-4ff8-b993-64c7e05c4dd9"]
+}, Open  ]],
+Cell[15580, 378, 3957, 88, 115, "Input",ExpressionUUID->"4b92dc3c-7bc7-4b5c-a5c3-e2aacf85ed52"],
+Cell[CellGroupData[{
+Cell[19562, 470, 1589, 38, 44, "Input",ExpressionUUID->"9c63e4ce-bd4c-4c9e-9b17-80e65259d14e"],
+Cell[21154, 510, 5371, 102, 268, "Output",ExpressionUUID->"79c580a9-d885-45d2-b73d-816f14fe791e"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[4485, 119, 3850, 86, 109, "Input",ExpressionUUID->"9e250256-8cd3-402a-bc21-3ee002bb1f82"],
-Cell[8338, 207, 1369, 23, 42, "Output",ExpressionUUID->"1fe1a989-6f4d-41a0-9680-c055dfc0798b"]
+Cell[26562, 617, 1597, 39, 44, "Input",ExpressionUUID->"9a7cc83f-920d-4c9e-9d42-998a933b225c"],
+Cell[28162, 658, 10024, 178, 264, "Output",ExpressionUUID->"72455824-b956-48f2-805c-31dcb93fc169"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[9744, 235, 3851, 86, 109, "Input",ExpressionUUID->"34fe8bc4-abfb-4504-a7f1-045cd2a55ee6"],
-Cell[13598, 323, 670, 14, 44, "Output",ExpressionUUID->"e9b15296-174c-4b92-bc6b-e3136dbafca2"]
+Cell[38223, 841, 3851, 86, 132, "Input",ExpressionUUID->"34fe8bc4-abfb-4504-a7f1-045cd2a55ee6"],
+Cell[42077, 929, 670, 14, 44, "Output",ExpressionUUID->"e9b15296-174c-4b92-bc6b-e3136dbafca2"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[14305, 342, 302, 8, 24, "Input",ExpressionUUID->"2e8e8755-9de4-46f0-8d34-70439b2b441b"],
-Cell[14610, 352, 206, 4, 25, "Output",ExpressionUUID->"f78d77d4-a743-42fe-aad7-f383120ebe66"]
+Cell[42784, 948, 302, 8, 24, "Input",ExpressionUUID->"2e8e8755-9de4-46f0-8d34-70439b2b441b"],
+Cell[43089, 958, 206, 4, 26, "Output",ExpressionUUID->"f78d77d4-a743-42fe-aad7-f383120ebe66"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[14853, 361, 551, 15, 24, "Input",ExpressionUUID->"e826229b-9150-4894-b288-5f834e4105f5"],
-Cell[15407, 378, 820, 23, 52, "Output",ExpressionUUID->"ba4067d3-83a2-4fb1-b40d-dda17275ca39"]
+Cell[43332, 967, 551, 15, 25, "Input",ExpressionUUID->"e826229b-9150-4894-b288-5f834e4105f5"],
+Cell[43886, 984, 820, 23, 54, "Output",ExpressionUUID->"ba4067d3-83a2-4fb1-b40d-dda17275ca39"]
 }, Open  ]],
-Cell[16242, 404, 2411, 65, 86, "Input",ExpressionUUID->"aac833ad-2b24-4138-bd6f-72fced45c28b"],
+Cell[44721, 1010, 2411, 65, 100, "Input",ExpressionUUID->"aac833ad-2b24-4138-bd6f-72fced45c28b"],
 Cell[CellGroupData[{
-Cell[18678, 473, 406, 9, 24, "Input",ExpressionUUID->"192d6104-6702-44ac-becd-0eb638099ac5"],
-Cell[19087, 484, 1368, 33, 40, "Output",ExpressionUUID->"d7b8e166-951c-43b4-b661-6bb65d35f1c8"]
+Cell[47157, 1079, 406, 9, 24, "Input",ExpressionUUID->"192d6104-6702-44ac-becd-0eb638099ac5"],
+Cell[47566, 1090, 1368, 33, 43, "Output",ExpressionUUID->"d7b8e166-951c-43b4-b661-6bb65d35f1c8"]
 }, Open  ]],
-Cell[20470, 520, 4890, 100, 86, "Input",ExpressionUUID->"1800b59b-60e8-4839-8aff-d82ae8bbba1a"],
+Cell[48949, 1126, 4890, 100, 100, "Input",ExpressionUUID->"1800b59b-60e8-4839-8aff-d82ae8bbba1a"],
 Cell[CellGroupData[{
-Cell[25385, 624, 874, 18, 24, "Input",ExpressionUUID->"04521673-5a45-41f8-ba39-1383a2daafcd"],
-Cell[26262, 644, 534, 12, 25, "Output",ExpressionUUID->"4344f182-cac0-40ef-9167-1216ed9218e0"]
+Cell[53864, 1230, 874, 18, 24, "Input",ExpressionUUID->"04521673-5a45-41f8-ba39-1383a2daafcd"],
+Cell[54741, 1250, 534, 12, 28, "Output",ExpressionUUID->"4344f182-cac0-40ef-9167-1216ed9218e0"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[26833, 661, 304, 7, 22, "Input",ExpressionUUID->"340432fd-b01b-4746-840a-e312eec4ade0"],
-Cell[27140, 670, 12065, 270, 177, "Output",ExpressionUUID->"0fe99cd0-b508-4c82-8cb5-d58d82b4dd4e"]
+Cell[55312, 1267, 304, 7, 22, "Input",ExpressionUUID->"340432fd-b01b-4746-840a-e312eec4ade0"],
+Cell[55619, 1276, 12065, 270, 179, "Output",ExpressionUUID->"0fe99cd0-b508-4c82-8cb5-d58d82b4dd4e"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[39242, 945, 281, 6, 22, "Input",ExpressionUUID->"33938e6c-adaa-464a-9564-ced71dd22d4e"],
-Cell[39526, 953, 238, 5, 25, "Output",ExpressionUUID->"321f4d39-6591-48a7-9389-9935b198b499"]
+Cell[67721, 1551, 281, 6, 22, "Input",ExpressionUUID->"33938e6c-adaa-464a-9564-ced71dd22d4e"],
+Cell[68005, 1559, 238, 5, 25, "Output",ExpressionUUID->"321f4d39-6591-48a7-9389-9935b198b499"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[39801, 963, 1589, 42, 92, "Input",ExpressionUUID->"5ed4b47f-e272-40ab-8c23-58385be3e0c6"],
-Cell[41393, 1007, 3039, 47, 95, "Message",ExpressionUUID->"06e3abf0-2019-4c5e-99b9-d2a9da7384db"],
-Cell[44435, 1056, 2133, 34, 57, "Message",ExpressionUUID->"465d3232-1f00-42ad-89f9-f7dbe9649552"],
-Cell[46571, 1092, 1416, 24, 39, "Message",ExpressionUUID->"51fecbe8-1bd4-4ff4-b67b-93b66d40b947"],
-Cell[47990, 1118, 628, 12, 22, "Message",ExpressionUUID->"4e40cd38-d9b9-4ca7-8ccd-de1b917cec88"],
-Cell[48621, 1132, 4228, 87, 182, "Output",ExpressionUUID->"02496185-9d40-4431-9ca9-6085c166d125"]
+Cell[68280, 1569, 1589, 42, 94, "Input",ExpressionUUID->"5ed4b47f-e272-40ab-8c23-58385be3e0c6"],
+Cell[69872, 1613, 3039, 47, 107, "Message",ExpressionUUID->"06e3abf0-2019-4c5e-99b9-d2a9da7384db"],
+Cell[72914, 1662, 2133, 34, 57, "Message",ExpressionUUID->"465d3232-1f00-42ad-89f9-f7dbe9649552"],
+Cell[75050, 1698, 1416, 24, 39, "Message",ExpressionUUID->"51fecbe8-1bd4-4ff4-b67b-93b66d40b947"],
+Cell[76469, 1724, 628, 12, 22, "Message",ExpressionUUID->"4e40cd38-d9b9-4ca7-8ccd-de1b917cec88"],
+Cell[77100, 1738, 4228, 87, 182, "Output",ExpressionUUID->"02496185-9d40-4431-9ca9-6085c166d125"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[52886, 1224, 241, 5, 22, "Input",ExpressionUUID->"eb2501cb-e52c-43ae-bbbf-2c91bfb69b81"],
-Cell[53130, 1231, 267, 7, 26, "Output",ExpressionUUID->"ddfad34c-6774-4ef8-b8bc-c5de09b41700"]
+Cell[81365, 1830, 241, 5, 22, "Input",ExpressionUUID->"eb2501cb-e52c-43ae-bbbf-2c91bfb69b81"],
+Cell[81609, 1837, 267, 7, 29, "Output",ExpressionUUID->"ddfad34c-6774-4ef8-b8bc-c5de09b41700"]
 }, Open  ]],
-Cell[53412, 1241, 4519, 90, 106, "Input",ExpressionUUID->"b7de9e01-4fea-49ea-97a1-b154bc61886b"],
+Cell[81891, 1847, 4519, 90, 112, "Input",ExpressionUUID->"b7de9e01-4fea-49ea-97a1-b154bc61886b"],
 Cell[CellGroupData[{
-Cell[57956, 1335, 1242, 35, 61, "Input",ExpressionUUID->"850f45b2-810d-4cc6-a9fa-3b038d25d4c3"],
-Cell[59201, 1372, 14451, 259, 182, "Output",ExpressionUUID->"664ef5af-b4d1-4773-bb48-a1fb5f99d3ea"]
+Cell[86435, 1941, 1242, 35, 42, "Input",ExpressionUUID->"850f45b2-810d-4cc6-a9fa-3b038d25d4c3"],
+Cell[87680, 1978, 14451, 259, 182, "Output",ExpressionUUID->"664ef5af-b4d1-4773-bb48-a1fb5f99d3ea"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[73689, 1636, 326, 7, 22, "Input",ExpressionUUID->"c2ec938b-df20-4457-ac8f-5fb89f93d9f3"],
-Cell[74018, 1645, 188, 4, 39, "Output",ExpressionUUID->"7a8c2bb5-90c4-4aeb-82ef-a923f5a729e3"]
+Cell[102168, 2242, 326, 7, 22, "Input",ExpressionUUID->"c2ec938b-df20-4457-ac8f-5fb89f93d9f3"],
+Cell[102497, 2251, 188, 4, 39, "Output",ExpressionUUID->"7a8c2bb5-90c4-4aeb-82ef-a923f5a729e3"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[74243, 1654, 255, 5, 22, "Input",ExpressionUUID->"f775ff7d-d0b6-4518-90fc-dce2b496637c"],
-Cell[74501, 1661, 194, 3, 39, "Output",ExpressionUUID->"60b8b008-7807-40a6-b3d9-1f2e8b589a73"]
+Cell[102722, 2260, 255, 5, 22, "Input",ExpressionUUID->"f775ff7d-d0b6-4518-90fc-dce2b496637c"],
+Cell[102980, 2267, 194, 3, 39, "Output",ExpressionUUID->"60b8b008-7807-40a6-b3d9-1f2e8b589a73"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[74732, 1669, 306, 7, 22, "Input",ExpressionUUID->"33e05555-1738-48c2-b7aa-e8fd44cecd72"],
-Cell[75041, 1678, 221, 3, 25, "Output",ExpressionUUID->"a53b9d58-d244-4548-8b2d-0153f5f04983"]
+Cell[103211, 2275, 306, 7, 22, "Input",ExpressionUUID->"33e05555-1738-48c2-b7aa-e8fd44cecd72"],
+Cell[103520, 2284, 221, 3, 25, "Output",ExpressionUUID->"a53b9d58-d244-4548-8b2d-0153f5f04983"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[75299, 1686, 4200, 108, 191, "Input",ExpressionUUID->"79242c4a-ef9e-453f-aabd-057c68f68aca"],
-Cell[79502, 1796, 19071, 333, 179, "Output",ExpressionUUID->"9e7d57d6-efa3-474c-ad30-e142da216a92"]
+Cell[103778, 2292, 4200, 108, 136, "Input",ExpressionUUID->"79242c4a-ef9e-453f-aabd-057c68f68aca"],
+Cell[107981, 2402, 19071, 333, 179, "Output",ExpressionUUID->"9e7d57d6-efa3-474c-ad30-e142da216a92"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[98610, 2134, 203, 4, 22, "Input",ExpressionUUID->"e2ebf167-6996-49ef-9ca3-9dd16132c4ea"],
-Cell[98816, 2140, 569, 14, 25, "Output",ExpressionUUID->"c33b6abd-b8d8-415e-b881-93fd6c52f0db"]
+Cell[127089, 2740, 203, 4, 22, "Input",ExpressionUUID->"e2ebf167-6996-49ef-9ca3-9dd16132c4ea"],
+Cell[127295, 2746, 569, 14, 25, "Output",ExpressionUUID->"c33b6abd-b8d8-415e-b881-93fd6c52f0db"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[99422, 2159, 501, 13, 24, "Input",ExpressionUUID->"e153916d-9f1b-43b9-a6c4-f59eefa45d16"],
-Cell[99926, 2174, 222, 4, 25, "Output",ExpressionUUID->"ca5cf147-c698-4d49-a753-1cf32997e5f5"]
+Cell[127901, 2765, 501, 13, 25, "Input",ExpressionUUID->"e153916d-9f1b-43b9-a6c4-f59eefa45d16"],
+Cell[128405, 2780, 222, 4, 25, "Output",ExpressionUUID->"ca5cf147-c698-4d49-a753-1cf32997e5f5"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[100185, 2183, 243, 5, 22, "Input",ExpressionUUID->"c29b0035-e029-4b18-8ef4-ecd12fb49333"],
-Cell[100431, 2190, 1247, 39, 70, "Output",ExpressionUUID->"a77b28a6-eacd-4614-b43a-ed39428597d6"]
+Cell[128664, 2789, 243, 5, 22, "Input",ExpressionUUID->"c29b0035-e029-4b18-8ef4-ecd12fb49333"],
+Cell[128910, 2796, 1247, 39, 70, "Output",ExpressionUUID->"a77b28a6-eacd-4614-b43a-ed39428597d6"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[101715, 2234, 212, 4, 22, "Input",ExpressionUUID->"40438766-386e-4ef0-95c5-bb36288917ef"],
-Cell[101930, 2240, 735, 15, 44, "Output",ExpressionUUID->"462d7199-fbaa-4949-8313-776399b35c4a"]
+Cell[130194, 2840, 212, 4, 22, "Input",ExpressionUUID->"40438766-386e-4ef0-95c5-bb36288917ef"],
+Cell[130409, 2846, 735, 15, 27, "Output",ExpressionUUID->"462d7199-fbaa-4949-8313-776399b35c4a"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[102702, 2260, 4388, 128, 175, "Input",ExpressionUUID->"f81ed9c7-d940-4314-963c-c36ea0215ef3"],
-Cell[107093, 2390, 1199, 20, 60, "Output",ExpressionUUID->"2234af5e-068b-46a3-beaf-99a907ef5184"]
+Cell[131181, 2866, 4385, 127, 145, "Input",ExpressionUUID->"f81ed9c7-d940-4314-963c-c36ea0215ef3"],
+Cell[135569, 2995, 1199, 20, 44, "Output",ExpressionUUID->"2234af5e-068b-46a3-beaf-99a907ef5184"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[108329, 2415, 689, 18, 24, "Input",ExpressionUUID->"922bbf13-0a40-4ffe-b65d-fb54bcae380b"],
-Cell[109021, 2435, 952, 16, 44, "Output",ExpressionUUID->"043e197a-7692-4187-83b6-ce5829ad4a5d"]
+Cell[136805, 3020, 689, 18, 24, "Input",ExpressionUUID->"922bbf13-0a40-4ffe-b65d-fb54bcae380b"],
+Cell[137497, 3040, 952, 16, 27, "Output",ExpressionUUID->"043e197a-7692-4187-83b6-ce5829ad4a5d"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[110010, 2456, 783, 19, 24, "Input",ExpressionUUID->"649e9e6d-3d2c-4639-9c47-2ea1ddcf4b02"],
-Cell[110796, 2477, 874, 15, 44, "Output",ExpressionUUID->"2bf68a8e-7713-4e0f-8e35-9b2bda25554f"]
+Cell[138486, 3061, 783, 19, 24, "Input",ExpressionUUID->"649e9e6d-3d2c-4639-9c47-2ea1ddcf4b02"],
+Cell[139272, 3082, 874, 15, 27, "Output",ExpressionUUID->"2bf68a8e-7713-4e0f-8e35-9b2bda25554f"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[111707, 2497, 756, 18, 24, "Input",ExpressionUUID->"65da12ab-b919-4aad-899a-77b90d54ba33"],
-Cell[112466, 2517, 783, 15, 27, "Output",ExpressionUUID->"bf118231-ac60-4b08-8d5b-c2f02ee72a5b"]
+Cell[140183, 3102, 756, 18, 24, "Input",ExpressionUUID->"65da12ab-b919-4aad-899a-77b90d54ba33"],
+Cell[140942, 3122, 783, 15, 27, "Output",ExpressionUUID->"bf118231-ac60-4b08-8d5b-c2f02ee72a5b"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[113286, 2537, 797, 19, 24, "Input",ExpressionUUID->"2476ffb9-8112-4899-a159-f51883b8e216"],
-Cell[114086, 2558, 687, 14, 27, "Output",ExpressionUUID->"4472fe7c-a13d-48f2-8ce9-7c86fc80c436"]
+Cell[141762, 3142, 797, 19, 24, "Input",ExpressionUUID->"2476ffb9-8112-4899-a159-f51883b8e216"],
+Cell[142562, 3163, 687, 14, 27, "Output",ExpressionUUID->"4472fe7c-a13d-48f2-8ce9-7c86fc80c436"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[114810, 2577, 819, 19, 24, "Input",ExpressionUUID->"e595b072-4f65-493a-afbb-d0fd95bb7c1c"],
-Cell[115632, 2598, 1078, 22, 60, "Output",ExpressionUUID->"fccdbef1-244e-44c6-b3c4-5001cf802c3b"]
+Cell[143286, 3182, 819, 19, 24, "Input",ExpressionUUID->"e595b072-4f65-493a-afbb-d0fd95bb7c1c"],
+Cell[144108, 3203, 1078, 22, 44, "Output",ExpressionUUID->"fccdbef1-244e-44c6-b3c4-5001cf802c3b"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[116747, 2625, 859, 20, 24, "Input",ExpressionUUID->"428ba532-a5a7-47a0-834d-51e3a093d96a"],
-Cell[117609, 2647, 963, 19, 60, "Output",ExpressionUUID->"a74c17d3-52e9-4bc4-8151-87ccfed07273"]
+Cell[145223, 3230, 859, 20, 24, "Input",ExpressionUUID->"428ba532-a5a7-47a0-834d-51e3a093d96a"],
+Cell[146085, 3252, 963, 19, 44, "Output",ExpressionUUID->"a74c17d3-52e9-4bc4-8151-87ccfed07273"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[118609, 2671, 891, 22, 41, "Input",ExpressionUUID->"11720fd4-85fa-4891-a6d6-e316b2c63a60"],
-Cell[119503, 2695, 559, 9, 25, "Output",ExpressionUUID->"bfc7450c-8753-4ad4-a680-b80674ac5a25"]
+Cell[147085, 3276, 891, 22, 24, "Input",ExpressionUUID->"11720fd4-85fa-4891-a6d6-e316b2c63a60"],
+Cell[147979, 3300, 559, 9, 25, "Output",ExpressionUUID->"bfc7450c-8753-4ad4-a680-b80674ac5a25"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[120099, 2709, 937, 22, 41, "Input",ExpressionUUID->"dd088f15-420b-4212-8186-b62091ac9d6c"],
-Cell[121039, 2733, 535, 9, 25, "Output",ExpressionUUID->"86d97ab8-bfe5-478d-a4ba-5bc364b4d54e"]
+Cell[148575, 3314, 937, 22, 24, "Input",ExpressionUUID->"dd088f15-420b-4212-8186-b62091ac9d6c"],
+Cell[149515, 3338, 535, 9, 25, "Output",ExpressionUUID->"86d97ab8-bfe5-478d-a4ba-5bc364b4d54e"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[121611, 2747, 1006, 24, 41, "Input",ExpressionUUID->"e7fdb778-fd9a-4c13-95dc-aefbf5b8a29f"],
-Cell[122620, 2773, 425, 7, 25, "Output",ExpressionUUID->"ac5b89c4-24ae-436b-b19b-5a19cf92f196"]
+Cell[150087, 3352, 1006, 24, 24, "Input",ExpressionUUID->"e7fdb778-fd9a-4c13-95dc-aefbf5b8a29f"],
+Cell[151096, 3378, 425, 7, 25, "Output",ExpressionUUID->"ac5b89c4-24ae-436b-b19b-5a19cf92f196"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[123082, 2785, 1052, 25, 41, "Input",ExpressionUUID->"3f815a4c-10cf-489f-9fc2-62e52069c4c0"],
-Cell[124137, 2812, 446, 8, 25, "Output",ExpressionUUID->"0058d268-597e-4f7a-a395-c0fdf2e2f0fc"]
+Cell[151558, 3390, 1052, 25, 24, "Input",ExpressionUUID->"3f815a4c-10cf-489f-9fc2-62e52069c4c0"],
+Cell[152613, 3417, 446, 8, 25, "Output",ExpressionUUID->"0058d268-597e-4f7a-a395-c0fdf2e2f0fc"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[124620, 2825, 850, 18, 41, "Input",ExpressionUUID->"cd97b893-1f61-415c-8763-743df3cfe10e"],
-Cell[125473, 2845, 524, 9, 25, "Output",ExpressionUUID->"9bbd03d3-a0ba-4157-8a3b-5f539ce94bce"]
+Cell[153096, 3430, 850, 18, 24, "Input",ExpressionUUID->"cd97b893-1f61-415c-8763-743df3cfe10e"],
+Cell[153949, 3450, 524, 9, 25, "Output",ExpressionUUID->"9bbd03d3-a0ba-4157-8a3b-5f539ce94bce"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[126034, 2859, 953, 20, 41, "Input",ExpressionUUID->"7cc8abed-034f-48f9-a04c-7a7cd336cffc"],
-Cell[126990, 2881, 311, 5, 25, "Output",ExpressionUUID->"07f77345-76ab-4e5d-a2e6-e24bb55feac9"]
+Cell[154510, 3464, 953, 20, 24, "Input",ExpressionUUID->"7cc8abed-034f-48f9-a04c-7a7cd336cffc"],
+Cell[155466, 3486, 311, 5, 25, "Output",ExpressionUUID->"07f77345-76ab-4e5d-a2e6-e24bb55feac9"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[127338, 2891, 847, 18, 41, "Input",ExpressionUUID->"ed90fc79-b0f9-464f-acb1-e5aa0d9691f1"],
-Cell[128188, 2911, 308, 5, 25, "Output",ExpressionUUID->"4514c09f-2394-460a-8977-2ada54b3b9d0"]
+Cell[155814, 3496, 847, 18, 24, "Input",ExpressionUUID->"ed90fc79-b0f9-464f-acb1-e5aa0d9691f1"],
+Cell[156664, 3516, 308, 5, 25, "Output",ExpressionUUID->"4514c09f-2394-460a-8977-2ada54b3b9d0"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[128533, 2921, 996, 20, 41, "Input",ExpressionUUID->"6dd628a1-23df-42ea-b634-4bef16413391"],
-Cell[129532, 2943, 549, 9, 25, "Output",ExpressionUUID->"80acb19e-2d21-40df-ae02-eb1e1c8d4e00"]
+Cell[157009, 3526, 996, 20, 24, "Input",ExpressionUUID->"6dd628a1-23df-42ea-b634-4bef16413391"],
+Cell[158008, 3548, 549, 9, 25, "Output",ExpressionUUID->"80acb19e-2d21-40df-ae02-eb1e1c8d4e00"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[130118, 2957, 228, 5, 22, "Input",ExpressionUUID->"feb08355-9ad1-4257-9f97-e394177d7e32"],
-Cell[130349, 2964, 274, 8, 25, "Output",ExpressionUUID->"87e29f54-9d68-4261-96d4-95d24dcdc36a"]
+Cell[158594, 3562, 228, 5, 22, "Input",ExpressionUUID->"feb08355-9ad1-4257-9f97-e394177d7e32"],
+Cell[158825, 3569, 274, 8, 28, "Output",ExpressionUUID->"87e29f54-9d68-4261-96d4-95d24dcdc36a"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[130660, 2977, 1542, 27, 55, "Input",ExpressionUUID->"e61ed118-a102-4290-8440-cba17b87145b"],
-Cell[132205, 3006, 12321, 270, 180, "Output",ExpressionUUID->"2e45389f-59c7-4eaf-9b3a-270214fa013e"]
+Cell[159136, 3582, 1542, 27, 55, "Input",ExpressionUUID->"e61ed118-a102-4290-8440-cba17b87145b"],
+Cell[160681, 3611, 12321, 270, 180, "Output",ExpressionUUID->"2e45389f-59c7-4eaf-9b3a-270214fa013e"]
 }, Open  ]],
-Cell[144541, 3279, 175, 3, 22, "Input",ExpressionUUID->"0f87aa05-2cc9-472a-ab27-39ee15777442"],
-Cell[144719, 3284, 356, 7, 24, "Input",ExpressionUUID->"627c17e0-e2d8-4347-8fe9-812183a42941"],
+Cell[173017, 3884, 175, 3, 22, "Input",ExpressionUUID->"0f87aa05-2cc9-472a-ab27-39ee15777442"],
+Cell[173195, 3889, 356, 7, 24, "Input",ExpressionUUID->"627c17e0-e2d8-4347-8fe9-812183a42941"],
 Cell[CellGroupData[{
-Cell[145100, 3295, 176, 2, 22, "Input",ExpressionUUID->"f32c32be-7fcb-4660-bf43-cbc119acda96"],
-Cell[145279, 3299, 2097, 65, 77, "Output",ExpressionUUID->"01b48fb7-d8c9-4012-92f6-488cf7e7087c"]
+Cell[173576, 3900, 176, 2, 22, "Input",ExpressionUUID->"f32c32be-7fcb-4660-bf43-cbc119acda96"],
+Cell[173755, 3904, 2097, 65, 47, "Output",ExpressionUUID->"01b48fb7-d8c9-4012-92f6-488cf7e7087c"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[147413, 3369, 860, 24, 24, "Input",ExpressionUUID->"bfa6a31f-5eb9-4063-b798-e45ab1e3b53a"],
-Cell[148276, 3395, 1944, 57, 77, "Output",ExpressionUUID->"0ca79efb-1867-4b1f-96fc-4a691d54b72e"]
+Cell[175889, 3974, 860, 24, 24, "Input",ExpressionUUID->"bfa6a31f-5eb9-4063-b798-e45ab1e3b53a"],
+Cell[176752, 4000, 1944, 57, 60, "Output",ExpressionUUID->"0ca79efb-1867-4b1f-96fc-4a691d54b72e"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[150257, 3457, 2622, 77, 93, "Input",ExpressionUUID->"242aacdf-e4c7-423c-a26a-12188ba2476e"],
-Cell[152882, 3536, 21838, 439, 180, "Output",ExpressionUUID->"073d831e-11a6-4203-90ce-0df429c1a3f0"]
+Cell[178733, 4062, 2622, 77, 93, "Input",ExpressionUUID->"242aacdf-e4c7-423c-a26a-12188ba2476e"],
+Cell[181358, 4141, 21835, 439, 180, "Output",ExpressionUUID->"073d831e-11a6-4203-90ce-0df429c1a3f0"]
 }, Open  ]],
-Cell[174735, 3978, 725, 15, 24, "Input",ExpressionUUID->"80c6c675-eade-40e1-8d0f-13a6578b30ff"],
-Cell[175463, 3995, 764, 21, 22, "Input",ExpressionUUID->"d8847f2e-21fe-4fb7-84ae-676f83acdf5a"],
+Cell[203208, 4583, 725, 15, 24, "Input",ExpressionUUID->"80c6c675-eade-40e1-8d0f-13a6578b30ff"],
+Cell[203936, 4600, 764, 21, 25, "Input",ExpressionUUID->"d8847f2e-21fe-4fb7-84ae-676f83acdf5a"],
 Cell[CellGroupData[{
-Cell[176252, 4020, 679, 17, 40, "Input",ExpressionUUID->"15d60f60-2f87-40da-a2c6-66d640111d75"],
-Cell[176934, 4039, 1298, 32, 44, "Output",ExpressionUUID->"6f1d902a-0431-4a99-adc2-97eb3c423317"]
+Cell[204725, 4625, 679, 17, 40, "Input",ExpressionUUID->"15d60f60-2f87-40da-a2c6-66d640111d75"],
+Cell[205407, 4644, 1298, 32, 48, "Output",ExpressionUUID->"6f1d902a-0431-4a99-adc2-97eb3c423317"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[178269, 4076, 736, 18, 57, "Input",ExpressionUUID->"bec4c929-c1a6-4fd0-be18-caf3643ba3a4"],
-Cell[179008, 4096, 14613, 314, 179, "Output",ExpressionUUID->"b48c9c3b-23cc-4864-8ea2-6e98855eb77d"]
+Cell[206742, 4681, 736, 18, 57, "Input",ExpressionUUID->"bec4c929-c1a6-4fd0-be18-caf3643ba3a4"],
+Cell[207481, 4701, 14613, 314, 179, "Output",ExpressionUUID->"b48c9c3b-23cc-4864-8ea2-6e98855eb77d"]
 }, Open  ]]
 }
 ]
diff --git a/caselle_compare_2nb.nb b/caselle_compare_2nb.nb
new file mode 100644
index 0000000..68c1caa
--- /dev/null
+++ b/caselle_compare_2nb.nb
@@ -0,0 +1,1255 @@
+(* Content-type: application/vnd.wolfram.mathematica *)
+
+(*** Wolfram Notebook File ***)
+(* http://www.wolfram.com/nb *)
+
+(* CreatedBy='Mathematica 13.3' *)
+
+(*CacheID: 234*)
+(* Internal cache information:
+NotebookFileLineBreakTest
+NotebookFileLineBreakTest
+NotebookDataPosition[       158,          7]
+NotebookDataLength[     51082,       1245]
+NotebookOptionsPosition[     46863,       1174]
+NotebookOutlinePosition[     47255,       1190]
+CellTagsIndexPosition[     47212,       1187]
+WindowFrame->Normal*)
+
+(* Beginning of Notebook Content *)
+Notebook[{
+Cell[BoxData[
+ RowBox[{"Get", "[", 
+  "\"\<~/doc/research/first_order_singularities/paper/IsingScalingFunction.wl\
+\>\"", "]"}]], "Input",
+ CellChangeTimes->{{3.904294245329715*^9, 3.904294272680094*^9}},
+ CellLabel->"In[1]:=",ExpressionUUID->"8861b958-81f2-4de2-b6f6-2949b5d18820"],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+ RowBox[{"castest", "=", 
+  RowBox[{
+   RowBox[{"Thread", "[", 
+    RowBox[{"{", 
+     RowBox[{
+      RowBox[{"Range", "[", 
+       RowBox[{"2", ",", "7"}], "]"}], ",", 
+      RowBox[{
+       RowBox[{"-", 
+        RowBox[{"Table", "[", 
+         RowBox[{
+          RowBox[{"1", "/", "2"}], ",", 
+          RowBox[{"{", 
+           RowBox[{"i", ",", "6"}], "}"}]}], "]"}]}], 
+       RowBox[{
+        RowBox[{"DScriptMCasD\[Eta]List", "[", 
+         RowBox[{"DataCas", "[", "2", "]"}], "]"}], "[", "5", "]"}]}]}], 
+     "}"}], "]"}], "//", "N"}]}]], "Input",
+ CellChangeTimes->{{3.904294276264942*^9, 3.904294310664956*^9}, {
+   3.904294355643044*^9, 3.904294367146028*^9}, {3.904294574669956*^9, 
+   3.904294581302743*^9}, {3.904294757978594*^9, 3.9042947589373713`*^9}, {
+   3.904294959349716*^9, 3.9042950197269907`*^9}, {3.904295069040153*^9, 
+   3.9042950707038517`*^9}, {3.9042951038167133`*^9, 3.904295139768777*^9}, 
+   3.904295178490686*^9, {3.9042958346868277`*^9, 3.904295843638595*^9}, {
+   3.90429733361961*^9, 3.904297333811123*^9}},
+ CellLabel->"In[2]:=",ExpressionUUID->"c8881b4c-983b-40d9-adac-dce1d52aa135"],
+
+Cell[BoxData[
+ RowBox[{"{", 
+  RowBox[{
+   RowBox[{"{", 
+    RowBox[{"2.`", ",", 
+     RowBox[{"-", "\[Infinity]"}]}], "}"}], ",", 
+   RowBox[{"{", 
+    RowBox[{"3.`", ",", "Indeterminate"}], "}"}], ",", 
+   RowBox[{"{", 
+    RowBox[{"4.`", ",", 
+     RowBox[{"-", "\[Infinity]"}]}], "}"}], ",", 
+   RowBox[{"{", 
+    RowBox[{"5.`", ",", "Indeterminate"}], "}"}], ",", 
+   RowBox[{"{", 
+    RowBox[{"6.`", ",", 
+     RowBox[{"-", "\[Infinity]"}]}], "}"}], ",", 
+   RowBox[{"{", 
+    RowBox[{"7.`", ",", "Indeterminate"}], "}"}]}], "}"}]], "Output",
+ CellChangeTimes->{{3.904295840517601*^9, 3.904295863121759*^9}, {
+   3.9042972009737062`*^9, 3.904297202948735*^9}, 3.904297316268012*^9, 
+   3.904297382202108*^9},
+ CellLabel->"Out[2]=",ExpressionUUID->"9e9a7f87-7c4b-4db8-b644-a242b762d190"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+ RowBox[{"ourtest", "=", 
+  RowBox[{
+   RowBox[{
+    RowBox[{"(", 
+     RowBox[{"DScriptF0D\[Eta]List", "@@", 
+      RowBox[{"PrepareArgument", "[", 
+       RowBox[{"Data", "[", "2", "]"}], "]"}]}], ")"}], "[", 
+    RowBox[{"5", ",", "1"}], "]"}], "//", "N"}]}]], "Input",
+ CellChangeTimes->{{3.904294860388315*^9, 3.9042949515897093`*^9}, {
+  3.904295044511445*^9, 3.9042950455511217`*^9}},
+ CellLabel->"In[6]:=",ExpressionUUID->"9e33e378-57cc-4077-8c23-ee4b0a4ba44d"],
+
+Cell[BoxData[
+ RowBox[{"{", 
+  RowBox[{
+   RowBox[{
+    RowBox[{"-", "1.1974378008271835`"}], "+", 
+    RowBox[{"0.`", " ", "\[ImaginaryI]"}]}], ",", 
+   RowBox[{
+    RowBox[{"-", "0.31907264609656283`"}], "+", 
+    RowBox[{"0.`", " ", "\[ImaginaryI]"}]}], ",", 
+   RowBox[{"0.11091443628541109`", "\[VeryThinSpace]", "+", 
+    RowBox[{"0.`", " ", "\[ImaginaryI]"}]}], ",", 
+   RowBox[{"0.016509650519848356`", "\[VeryThinSpace]", "+", 
+    RowBox[{"0.`", " ", "\[ImaginaryI]"}]}], ",", 
+   RowBox[{
+    RowBox[{"-", "0.00031411211128264577`"}], "+", 
+    RowBox[{"0.`", " ", "\[ImaginaryI]"}]}], ",", 
+   RowBox[{
+    RowBox[{"-", "0.0005146254923647323`"}], "+", 
+    RowBox[{"0.`", " ", "\[ImaginaryI]"}]}]}], "}"}]], "Output",
+ CellChangeTimes->{{3.90429490389577*^9, 3.9042949528420477`*^9}, {
+   3.904294998727845*^9, 3.904295046529663*^9}, {3.9042951169849367`*^9, 
+   3.904295142011281*^9}, 3.904295180276587*^9, 3.904295866984446*^9},
+ CellLabel->"Out[6]=",ExpressionUUID->"3b03069f-e883-42b0-a938-fedc5c785081"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+ RowBox[{"ListLogPlot", "[", 
+  RowBox[{"Abs", "@", 
+   RowBox[{"{", 
+    RowBox[{"\[CapitalPhi]s", ",", "ourtest", ",", "castest"}], "}"}]}], 
+  "]"}]], "Input",
+ CellChangeTimes->{{3.9042950344177027`*^9, 3.904295073567689*^9}},
+ CellLabel->"In[7]:=",ExpressionUUID->"f1660eeb-8ae8-4771-b70c-47e3a8eb7754"],
+
+Cell[BoxData[
+ GraphicsBox[{{{
+     {RGBColor[0.368417, 0.506779, 0.709798], 
+      StyleBox[{{}, {
+         LineBox[{{3., -2.199250858729474}, {3., -2.1992508587114377`}}], 
+         LineBox[{{3., -2.1992508587114377`}, {3., -2.199250858693401}}]}},
+       Antialiasing->False]}, 
+     {RGBColor[0.368417, 0.506779, 0.709798], 
+      StyleBox[{{}, {
+         LineBox[{{4., -4.108835370883289}, {4., -4.108835369665774}}], 
+         LineBox[{{4., -4.108835369665774}, {4., -4.108835368448258}}]}},
+       Antialiasing->False]}, 
+     {RGBColor[0.368417, 0.506779, 0.709798], 
+      StyleBox[{{}, {
+         LineBox[{{5., -8.23957016697712}, {5., -8.239569788186014}}], 
+         LineBox[{{5., -8.239569788186014}, {5., -8.23956940939505}}]}},
+       Antialiasing->False]}, 
+     {RGBColor[0.368417, 0.506779, 0.709798], 
+      StyleBox[{{}, {
+         LineBox[{{6., -7.57318515764723}, {6., -7.573184963114589}}], 
+         LineBox[{{6., -7.573184963114589}, {6., -7.573184768581987}}]}},
+       Antialiasing->False]}, 
+     {RGBColor[0.368417, 0.506779, 0.709798], 
+      StyleBox[{{}, {
+         LineBox[{{7., -8.473827235511466}, {7., -8.473822447718424}}], 
+         LineBox[{{7., -8.473822447718424}, {7., -8.473817659948304}}]}},
+       Antialiasing->False]}, 
+     {RGBColor[0.368417, 0.506779, 0.709798], 
+      StyleBox[{{}, {
+         LineBox[{{8., -10.012900713523315`}, {8., -10.012878401308493`}}], 
+         LineBox[{{8., -10.012878401308493`}, {8., -10.012856089591494`}}]}},
+       Antialiasing->False]}, 
+     {RGBColor[0.368417, 0.506779, 0.709798], 
+      StyleBox[{{}, {
+         LineBox[{{9., -14.970693198120777`}, {9., -14.967523623359499`}}], 
+         LineBox[{{9., -14.967523623359499`}, {9., -14.96436406306913}}]}},
+       Antialiasing->False]}, 
+     {RGBColor[0.368417, 0.506779, 0.709798], 
+      StyleBox[{{}, {
+         LineBox[{{10., -12.356895535264757`}, {10., -12.354572653848617`}}], 
+         LineBox[{{10., -12.354572653848617`}, {10., -12.352255155708255`}}]}},
+       Antialiasing->False]}, 
+     {RGBColor[0.368417, 0.506779, 0.709798], 
+      StyleBox[{{}, {
+         LineBox[{{11., -13.13241371325783}, {11., -13.127375919227873`}}], 
+         LineBox[{{11., -13.127375919227873`}, {11., -13.122363377404328`}}]}},
+       Antialiasing->False]}}, {
+     {RGBColor[0.368417, 0.506779, 0.709798], 
+      StyleBox[{
+        GeometricTransformationBox[
+         LineBox[{}], {{{1., 0.}, {0., 1.}}, {0., 0.}}], 
+        GeometricTransformationBox[
+         LineBox[{{{3., -2.199250858693401}, 
+            Offset[{3, 0}, {3., -2.199250858693401}]}, {{
+            3., -2.199250858693401}, 
+            Offset[{-3, 0}, {3., -2.199250858693401}]}, {{
+            3., -2.199250858729474}, 
+            Offset[{3, 0}, {3., -2.199250858729474}]}, {{
+            3., -2.199250858729474}, 
+            Offset[{-3, 0}, {3., -2.199250858729474}]}}], {{{1., 0.}, {0., 
+          1.}}, {0., 0.}}]},
+       Antialiasing->False]}, 
+     {RGBColor[0.368417, 0.506779, 0.709798], 
+      StyleBox[{
+        GeometricTransformationBox[
+         LineBox[{}], {{{1., 0.}, {0., 1.}}, {0., 0.}}], 
+        GeometricTransformationBox[
+         LineBox[{{{4., -4.108835368448258}, 
+            Offset[{3, 0}, {4., -4.108835368448258}]}, {{
+            4., -4.108835368448258}, 
+            Offset[{-3, 0}, {4., -4.108835368448258}]}, {{
+            4., -4.108835370883289}, 
+            Offset[{3, 0}, {4., -4.108835370883289}]}, {{
+            4., -4.108835370883289}, 
+            Offset[{-3, 0}, {4., -4.108835370883289}]}}], {{{1., 0.}, {0., 
+          1.}}, {0., 0.}}]},
+       Antialiasing->False]}, 
+     {RGBColor[0.368417, 0.506779, 0.709798], 
+      StyleBox[{
+        GeometricTransformationBox[
+         LineBox[{}], {{{1., 0.}, {0., 1.}}, {0., 0.}}], 
+        GeometricTransformationBox[
+         LineBox[{{{5., -8.23956940939505}, 
+            Offset[{3, 0}, {5., -8.23956940939505}]}, {{
+            5., -8.23956940939505}, 
+            Offset[{-3, 0}, {5., -8.23956940939505}]}, {{
+            5., -8.23957016697712}, 
+            Offset[{3, 0}, {5., -8.23957016697712}]}, {{
+            5., -8.23957016697712}, 
+            Offset[{-3, 0}, {5., -8.23957016697712}]}}], {{{1., 0.}, {0., 
+          1.}}, {0., 0.}}]},
+       Antialiasing->False]}, 
+     {RGBColor[0.368417, 0.506779, 0.709798], 
+      StyleBox[{
+        GeometricTransformationBox[
+         LineBox[{}], {{{1., 0.}, {0., 1.}}, {0., 0.}}], 
+        GeometricTransformationBox[
+         LineBox[{{{6., -7.573184768581987}, 
+            Offset[{3, 0}, {6., -7.573184768581987}]}, {{
+            6., -7.573184768581987}, 
+            Offset[{-3, 0}, {6., -7.573184768581987}]}, {{
+            6., -7.57318515764723}, 
+            Offset[{3, 0}, {6., -7.57318515764723}]}, {{
+            6., -7.57318515764723}, 
+            Offset[{-3, 0}, {6., -7.57318515764723}]}}], {{{1., 0.}, {0., 
+          1.}}, {0., 0.}}]},
+       Antialiasing->False]}, 
+     {RGBColor[0.368417, 0.506779, 0.709798], 
+      StyleBox[{
+        GeometricTransformationBox[
+         LineBox[{}], {{{1., 0.}, {0., 1.}}, {0., 0.}}], 
+        GeometricTransformationBox[
+         LineBox[{{{7., -8.473817659948304}, 
+            Offset[{3, 0}, {7., -8.473817659948304}]}, {{
+            7., -8.473817659948304}, 
+            Offset[{-3, 0}, {7., -8.473817659948304}]}, {{
+            7., -8.473827235511466}, 
+            Offset[{3, 0}, {7., -8.473827235511466}]}, {{
+            7., -8.473827235511466}, 
+            Offset[{-3, 0}, {7., -8.473827235511466}]}}], {{{1., 0.}, {0., 
+          1.}}, {0., 0.}}]},
+       Antialiasing->False]}, 
+     {RGBColor[0.368417, 0.506779, 0.709798], 
+      StyleBox[{
+        GeometricTransformationBox[
+         LineBox[{}], {{{1., 0.}, {0., 1.}}, {0., 0.}}], 
+        GeometricTransformationBox[
+         LineBox[{{{8., -10.012856089591494`}, 
+            Offset[{3, 0}, {8., -10.012856089591494`}]}, {{
+            8., -10.012856089591494`}, 
+            Offset[{-3, 0}, {8., -10.012856089591494`}]}, {{
+            8., -10.012900713523315`}, 
+            Offset[{3, 0}, {8., -10.012900713523315`}]}, {{
+            8., -10.012900713523315`}, 
+            Offset[{-3, 0}, {8., -10.012900713523315`}]}}], {{{1., 0.}, {0., 
+          1.}}, {0., 0.}}]},
+       Antialiasing->False]}, 
+     {RGBColor[0.368417, 0.506779, 0.709798], 
+      StyleBox[{
+        GeometricTransformationBox[
+         LineBox[{}], {{{1., 0.}, {0., 1.}}, {0., 0.}}], 
+        GeometricTransformationBox[
+         LineBox[{{{9., -14.96436406306913}, 
+            Offset[{3, 0}, {9., -14.96436406306913}]}, {{
+            9., -14.96436406306913}, 
+            Offset[{-3, 0}, {9., -14.96436406306913}]}, {{
+            9., -14.970693198120777`}, 
+            Offset[{3, 0}, {9., -14.970693198120777`}]}, {{
+            9., -14.970693198120777`}, 
+            Offset[{-3, 0}, {9., -14.970693198120777`}]}}], {{{1., 0.}, {0., 
+          1.}}, {0., 0.}}]},
+       Antialiasing->False]}, 
+     {RGBColor[0.368417, 0.506779, 0.709798], 
+      StyleBox[{
+        GeometricTransformationBox[
+         LineBox[{}], {{{1., 0.}, {0., 1.}}, {0., 0.}}], 
+        GeometricTransformationBox[
+         LineBox[{{{10., -12.352255155708255`}, 
+            Offset[{3, 0}, {10., -12.352255155708255`}]}, {{
+            10., -12.352255155708255`}, 
+            Offset[{-3, 0}, {10., -12.352255155708255`}]}, {{
+            10., -12.356895535264757`}, 
+            Offset[{3, 0}, {10., -12.356895535264757`}]}, {{
+            10., -12.356895535264757`}, 
+            Offset[{-3, 0}, {10., -12.356895535264757`}]}}], {{{1., 0.}, {0., 
+          1.}}, {0., 0.}}]},
+       Antialiasing->False]}, 
+     {RGBColor[0.368417, 0.506779, 0.709798], 
+      StyleBox[{
+        GeometricTransformationBox[
+         LineBox[{}], {{{1., 0.}, {0., 1.}}, {0., 0.}}], 
+        GeometricTransformationBox[
+         LineBox[{{{11., -13.122363377404328`}, 
+            Offset[{3, 0}, {11., -13.122363377404328`}]}, {{
+            11., -13.122363377404328`}, 
+            Offset[{-3, 0}, {11., -13.122363377404328`}]}, {{
+            11., -13.13241371325783}, 
+            Offset[{3, 0}, {11., -13.13241371325783}]}, {{
+            11., -13.13241371325783}, 
+            Offset[{-3, 0}, {11., -13.13241371325783}]}}], {{{1., 0.}, {0., 
+          1.}}, {0., 0.}}]},
+       Antialiasing->False]}}}, 
+   InterpretationBox[{
+     TagBox[{
+       TagBox[
+        {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
+         0.012833333333333334`], AbsoluteThickness[2], 
+         PointBox[{{1., 0.18043092383827275`}, {2., -1.1431595731895612`}, {
+          3., -2.1992508587114377`}, {4., -4.108835369665774}, {
+          5., -8.239569788186014}, {6., -7.573184963114589}, {
+          7., -8.473822447718424}, {8., -10.012878401308493`}, {
+          9., -14.967523623359499`}, {10., -12.354572653848617`}, {
+          11., -13.127375919227873`}}]},
+        Annotation[#, "Charting`Private`Tag#1"]& ], 
+       TagBox[
+        {RGBColor[0.880722, 0.611041, 0.142051], PointSize[
+         0.012833333333333334`], AbsoluteThickness[2], 
+         PointBox[{{1., 0.1801841080998339}, {2., -1.1423364714137663`}, {
+          3., -2.1989962191985946`}, {4., -4.103810189060431}, {
+          5., -8.065760593479036}, {6., -7.57207112112998}}]},
+        Annotation[#, "Charting`Private`Tag#2"]& ], 
+       TagBox[
+        {RGBColor[0.560181, 0.691569, 0.194885], PointSize[
+         0.012833333333333334`], AbsoluteThickness[2], 
+         PointBox[{{2., -0.9218075646421646}, {3., -5.349285388582857}, {
+          4., -4.037348264148434}, {5., -4.875683974574631}, {
+          6., -7.010966114383607}, {7., -7.45246624254697}}]},
+        Annotation[#, "Charting`Private`Tag#3"]& ]},
+      {"WolframDynamicHighlight", <|
+       "Label" -> {"XYLabel"}, "Ball" -> {"IndicatedBall"}|>}], 
+     DynamicBox[(Charting`HighlightActionBox["DynamicHighlight", {}, 
+       Slot["HighlightElements"], 
+       Slot["LayoutOptions"], 
+       Slot["Meta"], 
+       Charting`HighlightActionFunction["DynamicHighlight", {{
+          Annotation[{
+            Directive[
+             PointSize[0.012833333333333334`], 
+             RGBColor[0.368417, 0.506779, 0.709798], 
+             AbsoluteThickness[2]], 
+            
+            Point[{{1., 0.18043092383827275`}, {2., -1.1431595731895612`}, {
+             3., -2.1992508587114377`}, {4., -4.108835369665774}, {
+             5., -8.239569788186014}, {6., -7.573184963114589}, {
+             7., -8.473822447718424}, {8., -10.012878401308493`}, {
+             9., -14.967523623359499`}, {10., -12.354572653848617`}, {
+             11., -13.127375919227873`}}]}, "Charting`Private`Tag#1"], 
+          Annotation[{
+            Directive[
+             PointSize[0.012833333333333334`], 
+             RGBColor[0.880722, 0.611041, 0.142051], 
+             AbsoluteThickness[2]], 
+            
+            Point[{{1., 0.1801841080998339}, {2., -1.1423364714137663`}, {
+             3., -2.1989962191985946`}, {4., -4.103810189060431}, {
+             5., -8.065760593479036}, {6., -7.57207112112998}}]}, 
+           "Charting`Private`Tag#2"], 
+          Annotation[{
+            Directive[
+             PointSize[0.012833333333333334`], 
+             RGBColor[0.560181, 0.691569, 0.194885], 
+             AbsoluteThickness[2]], 
+            
+            Point[{{2., -0.9218075646421646}, {3., -5.349285388582857}, {
+             4., -4.037348264148434}, {5., -4.875683974574631}, {
+             6., -7.010966114383607}, {7., -7.45246624254697}}]}, 
+           "Charting`Private`Tag#3"]}}, <|
+        "HighlightElements" -> <|
+          "Label" -> {"XYLabel"}, "Ball" -> {"IndicatedBall"}|>, 
+         "LayoutOptions" -> <|
+          "PlotRange" -> {{0., 
+             11.}, {-16.156354643745974`, 
+              0.18043092383827276183497595745137112076`15.256310972792258}}, 
+           "Frame" -> {{False, False}, {False, False}}, 
+           "AxesOrigin" -> {0., -16.156354643745974`}, 
+           "ImageSize" -> {360, 360/GoldenRatio}, "Axes" -> {True, True}, 
+           "LabelStyle" -> {}, "AspectRatio" -> GoldenRatio^(-1), 
+           "DefaultStyle" -> {
+             Directive[
+              PointSize[0.012833333333333334`], 
+              RGBColor[0.560181, 0.691569, 0.194885], 
+              AbsoluteThickness[2]], 
+             Directive[
+              PointSize[0.012833333333333334`], 
+              RGBColor[0.880722, 0.611041, 0.142051], 
+              AbsoluteThickness[2]], 
+             Directive[
+              PointSize[0.012833333333333334`], 
+              RGBColor[0.368417, 0.506779, 0.709798], 
+              AbsoluteThickness[2]]}, 
+           "HighlightLabelingFunctions" -> <|"CoordinatesToolOptions" -> ({
+               Identity[
+                Part[#, 1]], 
+               Exp[
+                Part[#, 2]]}& ), 
+             "ScalingFunctions" -> {{Identity, Identity}, {Log, Exp}}|>, 
+           "Primitives" -> {}, "GCFlag" -> False|>, 
+         "Meta" -> <|
+          "DefaultHighlight" -> {"Dynamic", None}, "Index" -> {}, "Function" -> 
+           ListPlot, "GroupHighlight" -> False|>|>]]& )[<|
+       "HighlightElements" -> <|
+         "Label" -> {"XYLabel"}, "Ball" -> {"IndicatedBall"}|>, 
+        "LayoutOptions" -> <|
+         "PlotRange" -> {{0., 
+            11.}, {-16.156354643745974`, 
+             0.18043092383827276183497595745137112076`15.256310972792258}}, 
+          "Frame" -> {{False, False}, {False, False}}, 
+          "AxesOrigin" -> {0., -16.156354643745974`}, 
+          "ImageSize" -> {360, 360/GoldenRatio}, "Axes" -> {True, True}, 
+          "LabelStyle" -> {}, "AspectRatio" -> GoldenRatio^(-1), 
+          "DefaultStyle" -> {
+            Directive[
+             PointSize[0.012833333333333334`], 
+             RGBColor[0.560181, 0.691569, 0.194885], 
+             AbsoluteThickness[2]], 
+            Directive[
+             PointSize[0.012833333333333334`], 
+             RGBColor[0.880722, 0.611041, 0.142051], 
+             AbsoluteThickness[2]], 
+            Directive[
+             PointSize[0.012833333333333334`], 
+             RGBColor[0.368417, 0.506779, 0.709798], 
+             AbsoluteThickness[2]]}, 
+          "HighlightLabelingFunctions" -> <|"CoordinatesToolOptions" -> ({
+              Identity[
+               Part[#, 1]], 
+              Exp[
+               Part[#, 2]]}& ), 
+            "ScalingFunctions" -> {{Identity, Identity}, {Log, Exp}}|>, 
+          "Primitives" -> {}, "GCFlag" -> False|>, 
+        "Meta" -> <|
+         "DefaultHighlight" -> {"Dynamic", None}, "Index" -> {}, "Function" -> 
+          ListPlot, "GroupHighlight" -> False|>|>],
+      ImageSizeCache->{{4.503599627370496*^15, -4.503599627370496*^15}, {
+       4.503599627370496*^15, -4.503599627370496*^15}}]},
+    Annotation[{{
+       Annotation[{
+         Directive[
+          PointSize[0.012833333333333334`], 
+          RGBColor[0.368417, 0.506779, 0.709798], 
+          AbsoluteThickness[2]], 
+         Point[{{1., 0.18043092383827275`}, {2., -1.1431595731895612`}, {
+          3., -2.1992508587114377`}, {4., -4.108835369665774}, {
+          5., -8.239569788186014}, {6., -7.573184963114589}, {
+          7., -8.473822447718424}, {8., -10.012878401308493`}, {
+          9., -14.967523623359499`}, {10., -12.354572653848617`}, {
+          11., -13.127375919227873`}}]}, "Charting`Private`Tag#1"], 
+       Annotation[{
+         Directive[
+          PointSize[0.012833333333333334`], 
+          RGBColor[0.880722, 0.611041, 0.142051], 
+          AbsoluteThickness[2]], 
+         Point[{{1., 0.1801841080998339}, {2., -1.1423364714137663`}, {
+          3., -2.1989962191985946`}, {4., -4.103810189060431}, {
+          5., -8.065760593479036}, {6., -7.57207112112998}}]}, 
+        "Charting`Private`Tag#2"], 
+       Annotation[{
+         Directive[
+          PointSize[0.012833333333333334`], 
+          RGBColor[0.560181, 0.691569, 0.194885], 
+          AbsoluteThickness[2]], 
+         Point[{{2., -0.9218075646421646}, {3., -5.349285388582857}, {
+          4., -4.037348264148434}, {5., -4.875683974574631}, {
+          6., -7.010966114383607}, {7., -7.45246624254697}}]}, 
+        "Charting`Private`Tag#3"]}}, <|
+     "HighlightElements" -> <|
+       "Label" -> {"XYLabel"}, "Ball" -> {"IndicatedBall"}|>, 
+      "LayoutOptions" -> <|
+       "PlotRange" -> {{0., 
+          11.}, {-16.156354643745974`, 
+           0.18043092383827276183497595745137112076`15.256310972792258}}, 
+        "Frame" -> {{False, False}, {False, False}}, 
+        "AxesOrigin" -> {0., -16.156354643745974`}, 
+        "ImageSize" -> {360, 360/GoldenRatio}, "Axes" -> {True, True}, 
+        "LabelStyle" -> {}, "AspectRatio" -> GoldenRatio^(-1), "DefaultStyle" -> {
+          Directive[
+           PointSize[0.012833333333333334`], 
+           RGBColor[0.560181, 0.691569, 0.194885], 
+           AbsoluteThickness[2]], 
+          Directive[
+           PointSize[0.012833333333333334`], 
+           RGBColor[0.880722, 0.611041, 0.142051], 
+           AbsoluteThickness[2]], 
+          Directive[
+           PointSize[0.012833333333333334`], 
+           RGBColor[0.368417, 0.506779, 0.709798], 
+           AbsoluteThickness[2]]}, 
+        "HighlightLabelingFunctions" -> <|"CoordinatesToolOptions" -> ({
+            Identity[
+             Part[#, 1]], 
+            Exp[
+             Part[#, 2]]}& ), 
+          "ScalingFunctions" -> {{Identity, Identity}, {Log, Exp}}|>, 
+        "Primitives" -> {}, "GCFlag" -> False|>, 
+      "Meta" -> <|
+       "DefaultHighlight" -> {"Dynamic", None}, "Index" -> {}, "Function" -> 
+        ListPlot, "GroupHighlight" -> False|>|>, 
+     "DynamicHighlight"]], {{}, {}}},
+  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
+  Axes->{True, True},
+  AxesLabel->{None, None},
+  AxesOrigin->{0., -16.156354643745974`},
+  DisplayFunction->Identity,
+  Frame->{{False, False}, {False, False}},
+  FrameLabel->{{None, None}, {None, None}},
+  FrameTicks->{{
+     Charting`ScaledTicks[{Log, Exp}, {Log, Exp}, "Nice", WorkingPrecision -> 
+      15.954589770191003`, RotateLabel -> 0], 
+     Charting`ScaledFrameTicks[{Log, Exp}]}, {Automatic, Automatic}},
+  GridLines->{None, None},
+  GridLinesStyle->Directive[
+    GrayLevel[0.5, 0.4]],
+  Method->{
+   "AxisPadding" -> Scaled[0.02], "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], "DefaultPlotStyle" -> {
+      Directive[
+       RGBColor[0.368417, 0.506779, 0.709798], 
+       AbsoluteThickness[2]], 
+      Directive[
+       RGBColor[0.880722, 0.611041, 0.142051], 
+       AbsoluteThickness[2]], 
+      Directive[
+       RGBColor[0.560181, 0.691569, 0.194885], 
+       AbsoluteThickness[2]], 
+      Directive[
+       RGBColor[0.922526, 0.385626, 0.209179], 
+       AbsoluteThickness[2]], 
+      Directive[
+       RGBColor[0.528488, 0.470624, 0.701351], 
+       AbsoluteThickness[2]], 
+      Directive[
+       RGBColor[0.772079, 0.431554, 0.102387], 
+       AbsoluteThickness[2]], 
+      Directive[
+       RGBColor[0.363898, 0.618501, 0.782349], 
+       AbsoluteThickness[2]], 
+      Directive[
+       RGBColor[1, 0.75, 0], 
+       AbsoluteThickness[2]], 
+      Directive[
+       RGBColor[0.647624, 0.37816, 0.614037], 
+       AbsoluteThickness[2]], 
+      Directive[
+       RGBColor[0.571589, 0.586483, 0.], 
+       AbsoluteThickness[2]], 
+      Directive[
+       RGBColor[0.915, 0.3325, 0.2125], 
+       AbsoluteThickness[2]], 
+      Directive[
+       RGBColor[0.40082222609352647`, 0.5220066643438841, 0.85], 
+       AbsoluteThickness[2]], 
+      Directive[
+       RGBColor[0.9728288904374106, 0.621644452187053, 0.07336199581899142], 
+       AbsoluteThickness[2]], 
+      Directive[
+       RGBColor[0.736782672705901, 0.358, 0.5030266573755369], 
+       AbsoluteThickness[2]], 
+      Directive[
+       RGBColor[0.28026441037696703`, 0.715, 0.4292089322474965], 
+       AbsoluteThickness[2]]}, "DomainPadding" -> Scaled[0.02], 
+    "PointSizeFunction" -> "SmallPointSize", "RangePadding" -> Scaled[0.05], 
+    "OptimizePlotMarkers" -> True, "IncludeHighlighting" -> "CurrentPoint", 
+    "HighlightStyle" -> Automatic, "OptimizePlotMarkers" -> True, 
+    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
+        Identity[
+         Part[#, 1]], 
+        Exp[
+         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
+        Identity[
+         Part[#, 1]], 
+        Exp[
+         Part[#, 2]]}& )}},
+  PlotRange->{{0., 11.}, {-16.156354643745974`, 0.18043092383827275`}},
+  PlotRangeClipping->True,
+  PlotRangePadding->{{
+     Scaled[0.02], 
+     Scaled[0.02]}, {
+     Scaled[0.02], 
+     Scaled[0.05]}},
+  Ticks->FrontEndValueCache[{Automatic, 
+     Charting`ScaledTicks[{Log, Exp}, {Log, Exp}, "Nice", WorkingPrecision -> 
+      15.954589770191003`, RotateLabel -> 0]}, {
+    Automatic, {{-11.512925464970229`, 
+       FormBox[
+        TemplateBox[{"10", 
+          RowBox[{"-", "5"}]}, "Superscript", SyntaxForm -> SuperscriptBox], 
+        TraditionalForm], {0.01, 0.}}, {-4.605170185988091, 
+       FormBox[
+        TagBox[
+         InterpretationBox[
+          StyleBox["\"0.01\"", ShowStringCharacters -> False], 
+          0.01`15.954589770191003, AutoDelete -> True], NumberForm[#, {
+           DirectedInfinity[1], 2}]& ], TraditionalForm], {0.01, 
+       0.}}, {-18.420680743952367`, 
+       FormBox[
+        TemplateBox[{0, 0}, "Spacer2"], TraditionalForm], {0.005, 
+       0.}}, {-16.11809565095832, 
+       FormBox[
+        TemplateBox[{0, 0}, "Spacer2"], TraditionalForm], {0.005, 
+       0.}}, {-13.815510557964274`, 
+       FormBox[
+        TemplateBox[{0, 0}, "Spacer2"], TraditionalForm], {0.005, 
+       0.}}, {-9.210340371976182, 
+       FormBox[
+        TemplateBox[{0, 0}, "Spacer2"], TraditionalForm], {0.005, 
+       0.}}, {-6.907755278982137, 
+       FormBox[
+        TemplateBox[{0, 0}, "Spacer2"], TraditionalForm], {0.005, 
+       0.}}, {-2.3025850929940455`, 
+       FormBox[
+        TemplateBox[{0, 0}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {0., 
+       FormBox[
+        TemplateBox[{0, 0}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
+      2.302585092994046, 
+       FormBox[
+        TemplateBox[{0, 0}, "Spacer2"], TraditionalForm], {0.005, 
+       0.}}}}]]], "Output",
+ CellChangeTimes->{{3.904295051561068*^9, 3.904295073785268*^9}, {
+   3.904295117050866*^9, 3.904295142078001*^9}, 3.904295180341357*^9, 
+   3.904295867220212*^9},
+ CellLabel->"Out[7]=",ExpressionUUID->"247ec957-b845-4e37-a49c-a627e31cf9d6"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+ RowBox[{"Sqrt", "[", "1.16951", "]"}]], "Input",
+ CellChangeTimes->{{3.904295309316988*^9, 3.904295313676331*^9}, {
+  3.904295367862049*^9, 3.904295370165359*^9}, {3.904295411006793*^9, 
+  3.9042954146702623`*^9}, {3.9042954722160263`*^9, 
+  3.904295475087414*^9}},ExpressionUUID->"915f4570-3c11-443f-ba20-\
+0e101a65f85d"],
+
+Cell[BoxData["1.079078310411251`"], "Output",
+ CellChangeTimes->{3.904295313928844*^9, 3.904295370720072*^9, 
+  3.904295415247345*^9},
+ CellLabel->"Out[51]=",ExpressionUUID->"88b90945-96a7-44d9-ba64-fcd3bc5329f8"]
+}, Open  ]],
+
+Cell[BoxData["1"], "Input",
+ CellChangeTimes->{
+  3.904295156529663*^9},ExpressionUUID->"70f06980-fceb-4ccd-8db7-\
+3677340566ea"],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+ RowBox[{"\[CapitalPhi]s", "[", 
+  RowBox[{"[", "2", "]"}], "]"}]], "Input",
+ CellChangeTimes->{{3.904296019300721*^9, 3.9042960208259706`*^9}},
+ CellLabel->"In[20]:=",ExpressionUUID->"9495ae9c-88b9-4989-ad80-7c6dea597541"],
+
+Cell[BoxData[
+ RowBox[{"-", "0.31881012489061`"}]], "Output",
+ CellChangeTimes->{3.904296021007484*^9},
+ CellLabel->"Out[20]=",ExpressionUUID->"a322351c-5c4e-46e6-aa7b-88fe1c6c027e"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+ RowBox[{
+  RowBox[{
+   RowBox[{"DScriptMCasD\[Eta]List", "[", 
+    RowBox[{"DataCas", "[", "4", "]"}], "]"}], "[", 
+   RowBox[{"1", ",", "1"}], "]"}], "[", 
+  RowBox[{"[", "2", "]"}], "]"}]], "Input",
+ CellChangeTimes->{{3.9042960401946173`*^9, 3.904296055058302*^9}},
+ CellLabel->"In[24]:=",ExpressionUUID->"9452f3f3-3fa2-4aa8-80b7-294e0f586368"],
+
+Cell[BoxData["0.011323023797314785`"], "Output",
+ CellChangeTimes->{{3.904296030351758*^9, 3.9042960552267313`*^9}},
+ CellLabel->"Out[24]=",ExpressionUUID->"d9182823-a681-4a0d-9190-7835b33ee3dd"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+ RowBox[{
+  RowBox[{"DScriptMCasD\[Xi]List", "[", 
+   RowBox[{"DataCas", "[", "5", "]"}], "]"}], "[", 
+  RowBox[{"1", ",", "1"}], "]"}]], "Input",
+ CellChangeTimes->{{3.904296093868104*^9, 3.9042961241237803`*^9}},
+ CellLabel->"In[25]:=",ExpressionUUID->"09467d82-09bb-46f0-b554-8f62ebb9b1cd"],
+
+Cell[BoxData[
+ TemplateBox[{
+  "Power", "infy", 
+   "\"Infinite expression \\!\\(\\*FractionBox[\\\"1\\\", \
+SuperscriptBox[\\\"0\\\", RowBox[{\\\"15\\\", \\\"/\\\", \\\"8\\\"}]]]\\) \
+encountered.\"", 2, 25, 3, 23725455768631783963, "Local"},
+  "MessageTemplate"]], "Message", "MSG",
+ CellChangeTimes->{3.904296124472308*^9},
+ CellLabel->
+  "During evaluation of \
+In[25]:=",ExpressionUUID->"1265041d-082b-46cd-a8f0-dcfab5815f20"],
+
+Cell[BoxData[
+ TemplateBox[{
+  "Power", "infy", 
+   "\"Infinite expression \\!\\(\\*FractionBox[\\\"1\\\", \
+SuperscriptBox[\\\"0\\\", RowBox[{\\\"31\\\", \\\"/\\\", \\\"8\\\"}]]]\\) \
+encountered.\"", 2, 25, 4, 23725455768631783963, "Local"},
+  "MessageTemplate"]], "Message", "MSG",
+ CellChangeTimes->{3.904296124477688*^9},
+ CellLabel->
+  "During evaluation of \
+In[25]:=",ExpressionUUID->"b39558c0-2040-4fe7-afbe-3e5a7612d064"],
+
+Cell[BoxData[
+ TemplateBox[{
+  "Infinity", "indet", 
+   "\"Indeterminate expression \\!\\(\\*RowBox[{\\\"0.`\\\", \\\" \\\", \
+\\\"ComplexInfinity\\\"}]\\) encountered.\"", 2, 25, 5, 23725455768631783963, 
+   "Local"},
+  "MessageTemplate"]], "Message", "MSG",
+ CellChangeTimes->{3.904296124512054*^9},
+ CellLabel->
+  "During evaluation of \
+In[25]:=",ExpressionUUID->"978573a6-966f-4bf6-bbd1-573fa7bc88ab"],
+
+Cell[BoxData[
+ TemplateBox[{
+  "Power", "infy", 
+   "\"Infinite expression \\!\\(\\*FractionBox[\\\"1\\\", \
+SuperscriptBox[\\\"0\\\", RowBox[{\\\"15\\\", \\\"/\\\", \\\"8\\\"}]]]\\) \
+encountered.\"", 2, 25, 6, 23725455768631783963, "Local"},
+  "MessageTemplate"]], "Message", "MSG",
+ CellChangeTimes->{3.904296124516045*^9},
+ CellLabel->
+  "During evaluation of \
+In[25]:=",ExpressionUUID->"98fa9263-6a66-46e6-ac01-4b3ad3eb595f"],
+
+Cell[BoxData[
+ TemplateBox[{
+  "General", "stop", 
+   "\"Further output of \\!\\(\\*StyleBox[RowBox[{\\\"Power\\\", \\\"::\\\", \
+\\\"infy\\\"}], \\\"MessageName\\\"]\\) will be suppressed during this \
+calculation.\"", 2, 25, 7, 23725455768631783963, "Local"},
+  "MessageTemplate"]], "Message", "MSG",
+ CellChangeTimes->{3.904296124519102*^9},
+ CellLabel->
+  "During evaluation of \
+In[25]:=",ExpressionUUID->"1f113347-0bdd-4685-985a-2a875f0c613e"],
+
+Cell[BoxData[
+ TemplateBox[{
+  "Infinity", "indet", 
+   "\"Indeterminate expression \\!\\(\\*RowBox[{\\\"0.`\\\", \\\" \\\", \
+\\\"ComplexInfinity\\\"}]\\) encountered.\"", 2, 25, 8, 23725455768631783963, 
+   "Local"},
+  "MessageTemplate"]], "Message", "MSG",
+ CellChangeTimes->{3.904296124524341*^9},
+ CellLabel->
+  "During evaluation of \
+In[25]:=",ExpressionUUID->"2667a3df-89c1-494d-9667-2c9d5dbe1502"],
+
+Cell[BoxData[
+ TemplateBox[{
+  "Infinity", "indet", 
+   "\"Indeterminate expression \\!\\(\\*RowBox[{\\\"0.`\\\", \\\" \\\", \
+\\\"ComplexInfinity\\\"}]\\) encountered.\"", 2, 25, 9, 23725455768631783963, 
+   "Local"},
+  "MessageTemplate"]], "Message", "MSG",
+ CellChangeTimes->{3.904296124527913*^9},
+ CellLabel->
+  "During evaluation of \
+In[25]:=",ExpressionUUID->"7158adc8-878f-4330-9c68-c5bdb98230b2"],
+
+Cell[BoxData[
+ TemplateBox[{
+  "General", "stop", 
+   "\"Further output of \\!\\(\\*StyleBox[RowBox[{\\\"Infinity\\\", \
+\\\"::\\\", \\\"indet\\\"}], \\\"MessageName\\\"]\\) will be suppressed \
+during this calculation.\"", 2, 25, 10, 23725455768631783963, "Local"},
+  "MessageTemplate"]], "Message", "MSG",
+ CellChangeTimes->{3.9042961245311823`*^9},
+ CellLabel->
+  "During evaluation of \
+In[25]:=",ExpressionUUID->"39dd93a4-41e2-48ee-9e6e-314799940caf"],
+
+Cell[BoxData[
+ RowBox[{"{", 
+  RowBox[{"ComplexInfinity", ",", "Indeterminate"}], "}"}]], "Output",
+ CellChangeTimes->{3.904296124534574*^9},
+ CellLabel->"Out[25]=",ExpressionUUID->"8de7e965-03e3-431b-a4be-9da4e88bb702"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+ RowBox[{"m0Cas", "[", 
+  RowBox[{"DataCas", "[", "4", "]"}], "]"}]], "Input",
+ CellChangeTimes->{{3.904296133996367*^9, 3.904296143572279*^9}},
+ CellLabel->"In[26]:=",ExpressionUUID->"ef5f94b5-2b79-43b9-a8e6-89ab72ea29d7"],
+
+Cell[BoxData[
+ RowBox[{"m0Cas", "[", 
+  RowBox[{"\[LeftAssociation]", 
+   RowBox[{
+    RowBox[{"\<\"\[Theta]0\"\>", "\[Rule]", "1.079078310411251`"}], ",", 
+    RowBox[{"\<\"gs\"\>", "\[Rule]", 
+     RowBox[{"{", 
+      RowBox[{"1", ",", 
+       RowBox[{"-", "0.219388`"}], ",", 
+       RowBox[{"-", "0.041791`"}], ",", 
+       RowBox[{"-", "0.013488`"}]}], "}"}]}], ",", 
+    RowBox[{"\<\"\[Rho]\"\>", "\[Rule]", "2.0077`"}]}], 
+   "\[RightAssociation]"}], "]"}]], "Output",
+ CellChangeTimes->{3.904296144137517*^9},
+ CellLabel->"Out[26]=",ExpressionUUID->"f3f52d82-06af-4cfc-beb9-e8b699409c75"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+ RowBox[{"castest", "=", 
+  RowBox[{
+   RowBox[{"Table", "[", 
+    RowBox[{
+     RowBox[{"{", 
+      RowBox[{"k", ",", 
+       RowBox[{
+        RowBox[{
+         RowBox[{
+          RowBox[{
+           RowBox[{"DScriptMCasD\[Eta]List", "[", 
+            RowBox[{"DataCas", "[", "k", "]"}], "]"}], "[", 
+           RowBox[{"1", ",", "1"}], "]"}], "[", 
+          RowBox[{"[", "1", "]"}], "]"}], "/", "2"}], "+", 
+        RowBox[{"\[CapitalPhi]s", "[", 
+         RowBox[{"[", "2", "]"}], "]"}]}]}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{"k", ",", "2", ",", "5"}], "}"}]}], "]"}], "//", 
+   "N"}]}]], "Input",
+ CellChangeTimes->{{3.904295945120613*^9, 3.904296011066345*^9}, {
+  3.904297170680503*^9, 3.904297171023962*^9}},
+ CellLabel->"In[30]:=",ExpressionUUID->"aae85514-d241-40dd-a75a-048f5ee9648e"],
+
+Cell[BoxData[
+ TemplateBox[{
+  "Power", "infy", 
+   "\"Infinite expression \\!\\(\\*FractionBox[\\\"1\\\", \
+SuperscriptBox[\\\"0\\\", RowBox[{\\\"15\\\", \\\"/\\\", \\\"8\\\"}]]]\\) \
+encountered.\"", 2, 30, 19, 23725455768631783963, "Local"},
+  "MessageTemplate"]], "Message", "MSG",
+ CellChangeTimes->{{3.9042971426455708`*^9, 3.904297171216605*^9}},
+ CellLabel->
+  "During evaluation of \
+In[30]:=",ExpressionUUID->"630d7734-9c48-4e11-a310-b32933f2ee4a"],
+
+Cell[BoxData[
+ TemplateBox[{
+  "Power", "infy", 
+   "\"Infinite expression \\!\\(\\*FractionBox[\\\"1\\\", \
+SuperscriptBox[\\\"0\\\", RowBox[{\\\"31\\\", \\\"/\\\", \\\"8\\\"}]]]\\) \
+encountered.\"", 2, 30, 20, 23725455768631783963, "Local"},
+  "MessageTemplate"]], "Message", "MSG",
+ CellChangeTimes->{{3.9042971426455708`*^9, 3.904297171222148*^9}},
+ CellLabel->
+  "During evaluation of \
+In[30]:=",ExpressionUUID->"21d383f3-2ba5-4877-9527-1d64a782954d"],
+
+Cell[BoxData[
+ TemplateBox[{
+  "Infinity", "indet", 
+   "\"Indeterminate expression \\!\\(\\*RowBox[{\\\"0.`\\\", \\\" \\\", \
+\\\"ComplexInfinity\\\"}]\\) encountered.\"", 2, 30, 21, 23725455768631783963,
+    "Local"},
+  "MessageTemplate"]], "Message", "MSG",
+ CellChangeTimes->{{3.9042971426455708`*^9, 3.904297171225709*^9}},
+ CellLabel->
+  "During evaluation of \
+In[30]:=",ExpressionUUID->"486a6448-019f-42a4-85b2-35779b497211"],
+
+Cell[BoxData[
+ TemplateBox[{
+  "Power", "infy", 
+   "\"Infinite expression \\!\\(\\*FractionBox[\\\"1\\\", \
+SuperscriptBox[\\\"0\\\", RowBox[{\\\"15\\\", \\\"/\\\", \\\"8\\\"}]]]\\) \
+encountered.\"", 2, 30, 22, 23725455768631783963, "Local"},
+  "MessageTemplate"]], "Message", "MSG",
+ CellChangeTimes->{{3.9042971426455708`*^9, 3.904297171228862*^9}},
+ CellLabel->
+  "During evaluation of \
+In[30]:=",ExpressionUUID->"6f6b8276-f2c5-4162-80eb-18e9e87f9949"],
+
+Cell[BoxData[
+ TemplateBox[{
+  "General", "stop", 
+   "\"Further output of \\!\\(\\*StyleBox[RowBox[{\\\"Power\\\", \\\"::\\\", \
+\\\"infy\\\"}], \\\"MessageName\\\"]\\) will be suppressed during this \
+calculation.\"", 2, 30, 23, 23725455768631783963, "Local"},
+  "MessageTemplate"]], "Message", "MSG",
+ CellChangeTimes->{{3.9042971426455708`*^9, 3.9042971712319403`*^9}},
+ CellLabel->
+  "During evaluation of \
+In[30]:=",ExpressionUUID->"b628eb32-e81d-49d3-a1cd-149fbc61e9e5"],
+
+Cell[BoxData[
+ TemplateBox[{
+  "Infinity", "indet", 
+   "\"Indeterminate expression \\!\\(\\*RowBox[{\\\"0.`\\\", \\\" \\\", \
+\\\"ComplexInfinity\\\"}]\\) encountered.\"", 2, 30, 24, 23725455768631783963,
+    "Local"},
+  "MessageTemplate"]], "Message", "MSG",
+ CellChangeTimes->{{3.9042971426455708`*^9, 3.904297171235475*^9}},
+ CellLabel->
+  "During evaluation of \
+In[30]:=",ExpressionUUID->"265c2b61-3342-4898-8309-a0cad5834f2b"],
+
+Cell[BoxData[
+ TemplateBox[{
+  "Infinity", "indet", 
+   "\"Indeterminate expression \\!\\(\\*RowBox[{\\\"0.`\\\", \\\" \\\", \
+\\\"ComplexInfinity\\\"}]\\) encountered.\"", 2, 30, 25, 23725455768631783963,
+    "Local"},
+  "MessageTemplate"]], "Message", "MSG",
+ CellChangeTimes->{{3.9042971426455708`*^9, 3.9042971712386103`*^9}},
+ CellLabel->
+  "During evaluation of \
+In[30]:=",ExpressionUUID->"79976062-d0ab-49e3-a892-a33c3e1c43ff"],
+
+Cell[BoxData[
+ TemplateBox[{
+  "General", "stop", 
+   "\"Further output of \\!\\(\\*StyleBox[RowBox[{\\\"Infinity\\\", \
+\\\"::\\\", \\\"indet\\\"}], \\\"MessageName\\\"]\\) will be suppressed \
+during this calculation.\"", 2, 30, 26, 23725455768631783963, "Local"},
+  "MessageTemplate"]], "Message", "MSG",
+ CellChangeTimes->{{3.9042971426455708`*^9, 3.9042971712415447`*^9}},
+ CellLabel->
+  "During evaluation of \
+In[30]:=",ExpressionUUID->"ee547acd-2ff6-4b62-ab7a-257e268f3fb5"],
+
+Cell[BoxData[
+ RowBox[{"{", 
+  RowBox[{
+   RowBox[{"{", 
+    RowBox[{"2.`", ",", 
+     RowBox[{
+      RowBox[{"-", "0.31881012489061`"}], "+", 
+      RowBox[{"0.5`", " ", 
+       RowBox[{"(", 
+        RowBox[{
+         RowBox[{
+          RowBox[{"-", "0.039788735772973836`"}], " ", 
+          RowBox[{"Log", "[", 
+           SuperscriptBox[
+            RowBox[{"\[Eta]", "[", 
+             RowBox[{"1.0736759287606292`", ",", 
+              RowBox[{"{", 
+               RowBox[{"0.2731881074761689`", ",", 
+                RowBox[{"-", "0.05693458710289341`"}]}], "}"}]}], "]"}], "2"],
+            "]"}]}], "-", 
+         RowBox[{"1.197733383797993`", " ", 
+          SuperscriptBox[
+           RowBox[{"\[Eta]", "[", 
+            RowBox[{"1.0736759287606292`", ",", 
+             RowBox[{"{", 
+              RowBox[{"0.2731881074761689`", ",", 
+               RowBox[{"-", "0.05693458710289341`"}]}], "}"}]}], "]"}], 
+           "2"]}]}], ")"}]}]}]}], "}"}], ",", 
+   RowBox[{"{", 
+    RowBox[{"3.`", ",", 
+     RowBox[{
+      RowBox[{"-", "0.31881012489061`"}], "+", 
+      RowBox[{"0.5`", " ", 
+       RowBox[{"(", 
+        RowBox[{
+         RowBox[{
+          RowBox[{"-", "0.039788735772973836`"}], " ", 
+          RowBox[{"Log", "[", 
+           SuperscriptBox[
+            RowBox[{"\[Eta]", "[", 
+             RowBox[{"1.0767543823918249`", ",", 
+              RowBox[{"{", 
+               RowBox[{"0.272718114610426`", ",", 
+                RowBox[{"-", "0.05881847936860363`"}], ",", 
+                RowBox[{"-", "0.010745911869994617`"}]}], "}"}]}], "]"}], 
+            "2"], "]"}]}], "-", 
+         RowBox[{"1.197733383797993`", " ", 
+          SuperscriptBox[
+           RowBox[{"\[Eta]", "[", 
+            RowBox[{"1.0767543823918249`", ",", 
+             RowBox[{"{", 
+              RowBox[{"0.272718114610426`", ",", 
+               RowBox[{"-", "0.05881847936860363`"}], ",", 
+               RowBox[{"-", "0.010745911869994617`"}]}], "}"}]}], "]"}], 
+           "2"]}]}], ")"}]}]}]}], "}"}], ",", 
+   RowBox[{"{", 
+    RowBox[{"4.`", ",", 
+     RowBox[{
+      RowBox[{"-", "0.31881012489061`"}], "+", 
+      RowBox[{"0.5`", " ", 
+       RowBox[{"(", 
+        RowBox[{
+         RowBox[{
+          RowBox[{"-", "0.039788735772973836`"}], " ", 
+          RowBox[{"Log", "[", 
+           SuperscriptBox[
+            RowBox[{"\[Eta]", "[", 
+             RowBox[{"1.079078310411251`", ",", 
+              RowBox[{"{", 
+               RowBox[{"0.272718114610426`", ",", 
+                RowBox[{"-", "0.05983108172815214`"}], ",", 
+                RowBox[{"-", "0.011397162727684313`"}], ",", 
+                RowBox[{"-", "0.003678421929865426`"}]}], "}"}]}], "]"}], 
+            "2"], "]"}]}], "-", 
+         RowBox[{"1.197733383797993`", " ", 
+          SuperscriptBox[
+           RowBox[{"\[Eta]", "[", 
+            RowBox[{"1.079078310411251`", ",", 
+             RowBox[{"{", 
+              RowBox[{"0.272718114610426`", ",", 
+               RowBox[{"-", "0.05983108172815214`"}], ",", 
+               RowBox[{"-", "0.011397162727684313`"}], ",", 
+               RowBox[{"-", "0.003678421929865426`"}]}], "}"}]}], "]"}], 
+           "2"]}]}], ")"}]}]}]}], "}"}], ",", 
+   RowBox[{"{", 
+    RowBox[{"5.`", ",", 
+     RowBox[{
+      RowBox[{"-", "0.31881012489061`"}], "+", 
+      RowBox[{"0.5`", " ", 
+       RowBox[{"(", 
+        RowBox[{
+         RowBox[{
+          RowBox[{"-", "0.039788735772973836`"}], " ", 
+          RowBox[{"Log", "[", 
+           SuperscriptBox[
+            RowBox[{"\[Eta]", "[", 
+             RowBox[{"1.0814388563390904`", ",", 
+              RowBox[{"{", 
+               RowBox[{"0.27286889266851116`", ",", 
+                RowBox[{"-", "0.06068304017165753`"}], ",", 
+                RowBox[{"-", "0.011882621669035656`"}], ",", 
+                RowBox[{"-", "0.004040915431527982`"}], ",", 
+                RowBox[{"-", "0.001955924222647888`"}]}], "}"}]}], "]"}], 
+            "2"], "]"}]}], "-", 
+         RowBox[{"1.197733383797993`", " ", 
+          SuperscriptBox[
+           RowBox[{"\[Eta]", "[", 
+            RowBox[{"1.0814388563390904`", ",", 
+             RowBox[{"{", 
+              RowBox[{"0.27286889266851116`", ",", 
+               RowBox[{"-", "0.06068304017165753`"}], ",", 
+               RowBox[{"-", "0.011882621669035656`"}], ",", 
+               RowBox[{"-", "0.004040915431527982`"}], ",", 
+               RowBox[{"-", "0.001955924222647888`"}]}], "}"}]}], "]"}], 
+           "2"]}]}], ")"}]}]}]}], "}"}]}], "}"}]], "Output",
+ CellChangeTimes->{{3.904295966671215*^9, 3.9042960112588997`*^9}, {
+  3.904297142674638*^9, 3.904297171244854*^9}},
+ CellLabel->"Out[30]=",ExpressionUUID->"95aa2eb1-67be-41e1-b35f-bdc4394e9da1"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+ RowBox[{"ListPlot", "[", "castest", "]"}]], "Input",
+ CellChangeTimes->{{3.904295977553668*^9, 3.904295979304986*^9}},
+ CellLabel->"In[29]:=",ExpressionUUID->"76f935a1-981b-49a8-9e8b-337039807f78"],
+
+Cell[BoxData[
+ GraphicsBox[{{}, {}, {{}, {}}},
+  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
+  Axes->{True, True},
+  AxesLabel->{None, None},
+  AxesOrigin->{0, 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]],
+  Method->{
+   "AxisPadding" -> Scaled[0.02], "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], "DefaultPlotStyle" -> {
+      Directive[
+       RGBColor[0.368417, 0.506779, 0.709798], 
+       AbsoluteThickness[2]], 
+      Directive[
+       RGBColor[0.880722, 0.611041, 0.142051], 
+       AbsoluteThickness[2]], 
+      Directive[
+       RGBColor[0.560181, 0.691569, 0.194885], 
+       AbsoluteThickness[2]], 
+      Directive[
+       RGBColor[0.922526, 0.385626, 0.209179], 
+       AbsoluteThickness[2]], 
+      Directive[
+       RGBColor[0.528488, 0.470624, 0.701351], 
+       AbsoluteThickness[2]], 
+      Directive[
+       RGBColor[0.772079, 0.431554, 0.102387], 
+       AbsoluteThickness[2]], 
+      Directive[
+       RGBColor[0.363898, 0.618501, 0.782349], 
+       AbsoluteThickness[2]], 
+      Directive[
+       RGBColor[1, 0.75, 0], 
+       AbsoluteThickness[2]], 
+      Directive[
+       RGBColor[0.647624, 0.37816, 0.614037], 
+       AbsoluteThickness[2]], 
+      Directive[
+       RGBColor[0.571589, 0.586483, 0.], 
+       AbsoluteThickness[2]], 
+      Directive[
+       RGBColor[0.915, 0.3325, 0.2125], 
+       AbsoluteThickness[2]], 
+      Directive[
+       RGBColor[0.40082222609352647`, 0.5220066643438841, 0.85], 
+       AbsoluteThickness[2]], 
+      Directive[
+       RGBColor[0.9728288904374106, 0.621644452187053, 0.07336199581899142], 
+       AbsoluteThickness[2]], 
+      Directive[
+       RGBColor[0.736782672705901, 0.358, 0.5030266573755369], 
+       AbsoluteThickness[2]], 
+      Directive[
+       RGBColor[0.28026441037696703`, 0.715, 0.4292089322474965], 
+       AbsoluteThickness[2]]}, "DomainPadding" -> Scaled[0.02], 
+    "PointSizeFunction" -> "SmallPointSize", "RangePadding" -> Scaled[0.05], 
+    "OptimizePlotMarkers" -> True, "IncludeHighlighting" -> "CurrentPoint", 
+    "HighlightStyle" -> Automatic, "OptimizePlotMarkers" -> True, 
+    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
+        Identity[
+         Part[#, 1]], 
+        Identity[
+         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
+        Identity[
+         Part[#, 1]], 
+        Identity[
+         Part[#, 2]]}& )}},
+  PlotRange->{{0, 0.5463762149523378}, {-0.11386917420578682`, 0}},
+  PlotRangeClipping->True,
+  PlotRangePadding->{{
+     Scaled[0.02], 
+     Scaled[0.02]}, {
+     Scaled[0.05], 
+     Scaled[0.02]}},
+  Ticks->{Automatic, Automatic}]], "Output",
+ CellChangeTimes->{{3.904295979514624*^9, 3.904296011858904*^9}, 
+   3.904297165423973*^9},
+ CellLabel->"Out[29]=",ExpressionUUID->"3a315a7e-9c7c-48f2-afdd-5bff99d3cad0"]
+}, Open  ]]
+},
+WindowSize->{630, 1010.25},
+WindowMargins->{{0, Automatic}, {Automatic, 0}},
+FrontEndVersion->"13.3 for Linux x86 (64-bit) (July 24, 2023)",
+StyleDefinitions->"Default.nb",
+ExpressionUUID->"9a646958-7bcf-4952-a2b9-e394adf10ab8"
+]
+(* End of Notebook Content *)
+
+(* Internal cache information *)
+(*CellTagsOutline
+CellTagsIndex->{}
+*)
+(*CellTagsIndex
+CellTagsIndex->{}
+*)
+(*NotebookFileOutline
+Notebook[{
+Cell[558, 20, 282, 5, 22, "Input",ExpressionUUID->"8861b958-81f2-4de2-b6f6-2949b5d18820"],
+Cell[CellGroupData[{
+Cell[865, 29, 1137, 26, 22, "Input",ExpressionUUID->"c8881b4c-983b-40d9-adac-dce1d52aa135"],
+Cell[2005, 57, 792, 21, 25, "Output",ExpressionUUID->"9e9a7f87-7c4b-4db8-b644-a242b762d190"]
+}, Open  ]],
+Cell[CellGroupData[{
+Cell[2834, 83, 483, 11, 22, "Input",ExpressionUUID->"9e33e378-57cc-4077-8c23-ee4b0a4ba44d"],
+Cell[3320, 96, 1021, 22, 44, "Output",ExpressionUUID->"3b03069f-e883-42b0-a938-fedc5c785081"]
+}, Open  ]],
+Cell[CellGroupData[{
+Cell[4378, 123, 322, 7, 22, "Input",ExpressionUUID->"f1660eeb-8ae8-4771-b70c-47e3a8eb7754"],
+Cell[4703, 132, 22681, 507, 178, "Output",ExpressionUUID->"247ec957-b845-4e37-a49c-a627e31cf9d6"]
+}, Open  ]],
+Cell[CellGroupData[{
+Cell[27421, 644, 337, 6, 22, "Input",ExpressionUUID->"915f4570-3c11-443f-ba20-0e101a65f85d"],
+Cell[27761, 652, 213, 3, 25, "Output",ExpressionUUID->"88b90945-96a7-44d9-ba64-fcd3bc5329f8"]
+}, Open  ]],
+Cell[27989, 658, 129, 3, 22, "Input",ExpressionUUID->"70f06980-fceb-4ccd-8db7-3677340566ea"],
+Cell[CellGroupData[{
+Cell[28143, 665, 237, 4, 22, "Input",ExpressionUUID->"9495ae9c-88b9-4989-ad80-7c6dea597541"],
+Cell[28383, 671, 182, 3, 25, "Output",ExpressionUUID->"a322351c-5c4e-46e6-aa7b-88fe1c6c027e"]
+}, Open  ]],
+Cell[CellGroupData[{
+Cell[28602, 679, 362, 8, 22, "Input",ExpressionUUID->"9452f3f3-3fa2-4aa8-80b7-294e0f586368"],
+Cell[28967, 689, 195, 2, 25, "Output",ExpressionUUID->"d9182823-a681-4a0d-9190-7835b33ee3dd"]
+}, Open  ]],
+Cell[CellGroupData[{
+Cell[29199, 696, 307, 6, 22, "Input",ExpressionUUID->"09467d82-09bb-46f0-b554-8f62ebb9b1cd"],
+Cell[29509, 704, 431, 10, 32, "Message",ExpressionUUID->"1265041d-082b-46cd-a8f0-dcfab5815f20"],
+Cell[29943, 716, 431, 10, 32, "Message",ExpressionUUID->"b39558c0-2040-4fe7-afbe-3e5a7612d064"],
+Cell[30377, 728, 406, 10, 22, "Message",ExpressionUUID->"978573a6-966f-4bf6-bbd1-573fa7bc88ab"],
+Cell[30786, 740, 431, 10, 32, "Message",ExpressionUUID->"98fa9263-6a66-46e6-ac01-4b3ad3eb595f"],
+Cell[31220, 752, 449, 10, 22, "Message",ExpressionUUID->"1f113347-0bdd-4685-985a-2a875f0c613e"],
+Cell[31672, 764, 406, 10, 22, "Message",ExpressionUUID->"2667a3df-89c1-494d-9667-2c9d5dbe1502"],
+Cell[32081, 776, 406, 10, 22, "Message",ExpressionUUID->"7158adc8-878f-4330-9c68-c5bdb98230b2"],
+Cell[32490, 788, 456, 10, 22, "Message",ExpressionUUID->"39dd93a4-41e2-48ee-9e6e-314799940caf"],
+Cell[32949, 800, 220, 4, 25, "Output",ExpressionUUID->"8de7e965-03e3-431b-a4be-9da4e88bb702"]
+}, Open  ]],
+Cell[CellGroupData[{
+Cell[33206, 809, 237, 4, 22, "Input",ExpressionUUID->"ef5f94b5-2b79-43b9-a8e6-89ab72ea29d7"],
+Cell[33446, 815, 596, 14, 25, "Output",ExpressionUUID->"f3f52d82-06af-4cfc-beb9-e8b699409c75"]
+}, Open  ]],
+Cell[CellGroupData[{
+Cell[34079, 834, 821, 22, 22, "Input",ExpressionUUID->"aae85514-d241-40dd-a75a-048f5ee9648e"],
+Cell[34903, 858, 458, 10, 32, "Message",ExpressionUUID->"630d7734-9c48-4e11-a310-b32933f2ee4a"],
+Cell[35364, 870, 458, 10, 32, "Message",ExpressionUUID->"21d383f3-2ba5-4877-9527-1d64a782954d"],
+Cell[35825, 882, 433, 10, 22, "Message",ExpressionUUID->"486a6448-019f-42a4-85b2-35779b497211"],
+Cell[36261, 894, 458, 10, 32, "Message",ExpressionUUID->"6f6b8276-f2c5-4162-80eb-18e9e87f9949"],
+Cell[36722, 906, 478, 10, 22, "Message",ExpressionUUID->"b628eb32-e81d-49d3-a1cd-149fbc61e9e5"],
+Cell[37203, 918, 433, 10, 22, "Message",ExpressionUUID->"265c2b61-3342-4898-8309-a0cad5834f2b"],
+Cell[37639, 930, 435, 10, 22, "Message",ExpressionUUID->"79976062-d0ab-49e3-a892-a33c3e1c43ff"],
+Cell[38077, 942, 482, 10, 22, "Message",ExpressionUUID->"ee547acd-2ff6-4b62-ab7a-257e268f3fb5"],
+Cell[38562, 954, 4694, 117, 163, "Output",ExpressionUUID->"95aa2eb1-67be-41e1-b35f-bdc4394e9da1"]
+}, Open  ]],
+Cell[CellGroupData[{
+Cell[43293, 1076, 212, 3, 22, "Input",ExpressionUUID->"76f935a1-981b-49a8-9e8b-337039807f78"],
+Cell[43508, 1081, 3339, 90, 169, "Output",ExpressionUUID->"3a315a7e-9c7c-48f2-afdd-5bff99d3cad0"]
+}, Open  ]]
+}
+]
+*)
+
+(* End of internal cache information *)
+
author	Jaron Kent-Dobias <jaron@kent-dobias.com>	2023-12-11 18:06:08 +0100
committer	Jaron Kent-Dobias <jaron@kent-dobias.com>	2023-12-11 18:06:08 +0100
commit	f9fc90631253da8bea10694fe0d90dd393c546e6 (patch)
tree	a86403e382f66a945d4b7974aeee9354ab5f0db5
parent	8aeb86c75c525887119655f72a5f2ef5d68b0940 (diff)
download	mma-f9fc90631253da8bea10694fe0d90dd393c546e6.tar.gz mma-f9fc90631253da8bea10694fe0d90dd393c546e6.tar.bz2 mma-f9fc90631253da8bea10694fe0d90dd393c546e6.zip