Work.

1 files changed, 758 insertions, 198 deletions

diff --git a/caselle_compare.nb b/caselle_compare.nb
index d2c852f..9172109 100644
--- a/caselle_compare.nb
+++ b/caselle_compare.nb
@@ -10,10 +10,10 @@
 NotebookFileLineBreakTest
 NotebookFileLineBreakTest
 NotebookDataPosition[       158,          7]
-NotebookDataLength[    105372,       2426]
-NotebookOptionsPosition[    100295,       2338]
-NotebookOutlinePosition[    100694,       2354]
-CellTagsIndexPosition[    100651,       2351]
+NotebookDataLength[    129143,       2986]
+NotebookOptionsPosition[    122955,       2880]
+NotebookOutlinePosition[    123353,       2896]
+CellTagsIndexPosition[    123310,       2893]
 WindowFrame->Normal*)
 
 (* Beginning of Notebook Content *)
@@ -267,7 +267,7 @@ Cell[BoxData[{
       SuperscriptBox["\[Theta]", "8"]}]}], ")"}]}]}]}], "Input",
  CellChangeTimes->{{3.8321454371387777`*^9, 3.832145559020133*^9}, 
    3.832146444152248*^9, {3.832146620327344*^9, 3.832146620479251*^9}},
- CellLabel->"In[6]:=",ExpressionUUID->"04b6d9c8-662b-43f2-955f-87d9d6c7b622"],
+ CellLabel->"In[46]:=",ExpressionUUID->"04b6d9c8-662b-43f2-955f-87d9d6c7b622"],
 
 Cell[BoxData[
  RowBox[{
@@ -289,7 +289,7 @@ Cell[BoxData[
   3.83214577961791*^9}, {3.832145890778178*^9, 3.832145905218177*^9}, {
   3.832145940706725*^9, 3.8321459464268084`*^9}, {3.832145996299863*^9, 
   3.832146000067889*^9}, {3.832239324678667*^9, 3.832239324966337*^9}},
- CellLabel->"In[8]:=",ExpressionUUID->"bd6b698e-6871-43fa-9d98-7385cc594c9c"],
+ CellLabel->"In[48]:=",ExpressionUUID->"bd6b698e-6871-43fa-9d98-7385cc594c9c"],
 
 Cell[BoxData[
  RowBox[{
@@ -317,12 +317,15 @@ Cell[BoxData[
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Range", "[", 
-       RowBox[{"0", ",", "11"}], "]"}], ",", "Gls"}], "}"}], "]"}]}], 
-  ";"}]], "Input",
+       RowBox[{"0", ",", 
+        RowBox[{
+         RowBox[{"Length", "[", "Gls", "]"}], "-", "1"}]}], "]"}], ",", 
+      "Gls"}], "}"}], "]"}]}], ";"}]], "Input",
  CellChangeTimes->{{3.832146079469277*^9, 3.832146080085175*^9}, {
    3.832146220615906*^9, 3.832146230663995*^9}, 3.836571905035768*^9, {
-   3.836572227321196*^9, 3.836572228850659*^9}},
- CellLabel->"In[11]:=",ExpressionUUID->"ec8b9383-61d6-44c2-b776-3d651fc19660"],
+   3.836572227321196*^9, 3.836572228850659*^9}, {3.843802149837371*^9, 
+   3.843802151956236*^9}},
+ CellLabel->"In[59]:=",ExpressionUUID->"ec8b9383-61d6-44c2-b776-3d651fc19660"],
 
 Cell[CellGroupData[{
 
@@ -367,16 +370,14 @@ Cell[BoxData[
    RowBox[{"Solve", "[", 
     RowBox[{
      RowBox[{
-      RowBox[{"coeffCas", "[", "1", "]"}], "==", 
-      RowBox[{
-       RowBox[{"{", 
+      RowBox[{"datCas", "[", 
+       RowBox[{"[", 
+        RowBox[{"{", 
+         RowBox[{"1", ",", "2"}], "}"}], "]"}], "]"}], "==", 
+      RowBox[{"datNum", "[", 
+       RowBox[{"[", 
         RowBox[{
-         RowBox[{"1", "!"}], ",", 
-         RowBox[{"2", "!"}]}], "}"}], 
-       RowBox[{"datNum", "[", 
-        RowBox[{"[", 
-         RowBox[{
-          RowBox[{"2", ";;", "3"}], ",", "2"}], "]"}], "]"}]}]}], ",", 
+         RowBox[{"2", ";;", "3"}], ",", "2"}], "]"}], "]"}]}], ",", 
      RowBox[{"{", 
       RowBox[{"m0", ",", "h0"}], "}"}]}], "]"}], "[", 
    RowBox[{"[", "1", "]"}], "]"}]}]], "Input",
@@ -385,15 +386,16 @@ Cell[BoxData[
   3.832146653175795*^9}, {3.832151934783225*^9, 3.8321519353348494`*^9}, {
   3.832152037225222*^9, 3.8321520374089203`*^9}, {3.832239305719432*^9, 
   3.832239316046335*^9}, {3.836984166958787*^9, 3.83698416798269*^9}, {
-  3.836984311313539*^9, 3.836984316617137*^9}},
- CellLabel->"In[12]:=",ExpressionUUID->"7efd85be-fd7d-435a-a870-214d8ec4025d"],
+  3.836984311313539*^9, 3.836984316617137*^9}, {3.8438022146151648`*^9, 
+  3.843802259206324*^9}},
+ CellLabel->"In[72]:=",ExpressionUUID->"7efd85be-fd7d-435a-a870-214d8ec4025d"],
 
 Cell[BoxData[
  TemplateBox[{
   "Solve", "ratnz", 
    "\"Solve was unable to solve the system with inexact coefficients. The \
 answer was obtained by solving a corresponding exact system and numericizing \
-the result.\"", 2, 12, 3, 31608996298011231244, "Local"},
+the result.\"", 2, 72, 17, 31653097672855845075, "Local"},
   "MessageTemplate"]], "Message", "MSG",
  CellChangeTimes->{
   3.832146653505703*^9, 3.8321519356060266`*^9, 3.8321520375652533`*^9, {
@@ -401,25 +403,18 @@ the result.\"", 2, 12, 3, 31608996298011231244, "Local"},
    3.832239326773054*^9}, {3.836570087200712*^9, 3.836570104160351*^9}, {
    3.8369749921742907`*^9, 3.836975002927493*^9}, 3.836984168350904*^9, 
    3.8369843170148077`*^9, 3.8370700000102262`*^9, {3.837072474483852*^9, 
-   3.837072478715604*^9}},
+   3.837072478715604*^9}, {3.8438021401706963`*^9, 3.843802153865232*^9}, {
+   3.843802229529984*^9, 3.843802270916144*^9}},
  CellLabel->
   "During evaluation of \
-In[12]:=",ExpressionUUID->"2e43731b-67a2-4964-9c03-694bb26aa975"],
+In[72]:=",ExpressionUUID->"4bba24d3-fa02-4ad0-86b9-3bdb7e46a000"],
 
 Cell[BoxData[
  RowBox[{"{", 
   RowBox[{
-   RowBox[{"m0", "\[Rule]", 
-    InterpretationBox[
-     TemplateBox[{"-1.00576267284062", "0", "7"},
-      "AroundSeq"],
-     Around[-1.0057626728406195`, 7.407086999594733*^-15]]}], ",", 
+   RowBox[{"m0", "\[Rule]", "566.4709429312807`"}], ",", 
    RowBox[{"h0", "\[Rule]", 
-    InterpretationBox[
-     TemplateBox[{"0.2591817983", "42", "11"},
-      "AroundSeq"],
-     Around[0.2591817983420855, 1.0588943216503216`*^-11]]}]}], 
-  "}"}]], "Output",
+    RowBox[{"-", "5.262843279675711`"}]}]}], "}"}]], "Output",
  CellChangeTimes->{{3.832145993018977*^9, 3.832146075354514*^9}, {
    3.832146456457759*^9, 3.83214647211082*^9}, 3.832146653511766*^9, 
    3.832151935612179*^9, 3.8321520375710917`*^9, {3.8321544206488857`*^9, 
@@ -427,8 +422,9 @@ Cell[BoxData[
    3.8322393267777767`*^9}, {3.836570087516891*^9, 3.8365701041934032`*^9}, {
    3.836974992241055*^9, 3.836975002932447*^9}, 3.836984168357218*^9, 
    3.836984317019814*^9, 3.837070000054738*^9, {3.837072474543331*^9, 
-   3.837072478720648*^9}},
- CellLabel->"Out[12]=",ExpressionUUID->"325e2b46-9d64-4724-8091-36718f8f11a2"]
+   3.837072478720648*^9}, {3.843802140320195*^9, 3.843802153880229*^9}, {
+   3.843802229536105*^9, 3.843802270921445*^9}},
+ CellLabel->"Out[72]=",ExpressionUUID->"e2f0010a-e1b6-4537-a8ae-b59af3108330"]
 }, Open  ]],
 
 Cell[CellGroupData[{
@@ -488,164 +484,710 @@ Cell[BoxData[
  CellLabel->"Out[16]=",ExpressionUUID->"c7c2e58e-e6d9-436b-9cc9-4238f4205a5a"]
 }, Open  ]],
 
-Cell[BoxData[
+Cell[CellGroupData[{
+
+Cell[BoxData[{
  RowBox[{
-  RowBox[{
-   RowBox[{
-    RowBox[{"invDerivativeList", "[", "n_", "]"}], "[", "f_", "]"}], "[", 
-   "x_", "]"}], ":=", 
-  RowBox[{"Module", "[", "\[IndentingNewLine]", 
+  RowBox[{"totCas", "=", "10"}], ";"}], "\[IndentingNewLine]", 
+ RowBox[{
+  RowBox[{"ds", "=", 
    RowBox[{
-    RowBox[{"{", 
-     RowBox[{"xp", ",", "dfs", ",", "fp", ",", "Pns"}], "}"}], ",", 
-    "\[IndentingNewLine]", 
     RowBox[{
-     RowBox[{"dfs", "=", 
+     RowBox[{"invDerivativeList", "[", 
+      RowBox[{"totCas", "+", "1"}], "]"}], "[", 
+     RowBox[{"\[Theta]", "\[Function]", 
+      FractionBox[
+       RowBox[{"hCas", "[", "\[Theta]", "]"}], 
+       SuperscriptBox[
+        RowBox[{"(", 
+         RowBox[{"t", "[", "\[Theta]", "]"}], ")"}], 
+        RowBox[{"\[CapitalDelta]", "[", "2", "]"}]]]}], "]"}], "[", 
+    "\[Theta]0Cas", "]"}]}], ";"}], "\[IndentingNewLine]", 
+ RowBox[{
+  RowBox[{"dMs", "=", 
+   RowBox[{
+    RowBox[{"m0", " ", 
+     RowBox[{"NestList", "[", 
       RowBox[{
-       RowBox[{"Rest", "[", 
-        RowBox[{"NestList", "[", 
-         RowBox[{
+       RowBox[{
+        RowBox[{"D", "[", 
+         RowBox[{"#", ",", "\[Theta]"}], "]"}], "&"}], ",", 
+       FractionBox["\[Theta]", 
+        SuperscriptBox[
+         RowBox[{"(", 
           RowBox[{
-           RowBox[{"D", "[", 
-            RowBox[{"#", ",", "xp"}], "]"}], "&"}], ",", 
-          RowBox[{"f", "[", "xp", "]"}], ",", "n"}], "]"}], "]"}], "/.", 
-       RowBox[{"xp", "->", "x"}]}]}], ";", "\[IndentingNewLine]", 
-     RowBox[{"Pns", "=", 
-      RowBox[{
-       RowBox[{"FoldList", "[", 
+           SuperscriptBox["\[Theta]", "2"], "-", "1"}], ")"}], 
+         RowBox[{"1", "/", "8"}]]], ",", "totCas"}], "]"}]}], "/.", 
+    RowBox[{"\[Theta]", "->", "\[Theta]0Cas"}]}]}], 
+  ";"}], "\[IndentingNewLine]", 
+ RowBox[{"datCas", "=", 
+  RowBox[{"Table", "[", 
+   RowBox[{
+    RowBox[{
+     RowBox[{"-", 
+      RowBox[{"Sum", "[", 
+       RowBox[{
         RowBox[{
-         RowBox[{
-          RowBox[{"{", 
-           RowBox[{"Pm", ",", "m"}], "}"}], "\[Function]", " ", 
+         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}},
+ CellLabel->
+  "In[141]:=",ExpressionUUID->"9e250256-8cd3-402a-bc21-3ee002bb1f82"],
+
+Cell[BoxData[
+ InterpretationBox[
+  TagBox[
+   FrameBox[GridBox[{
+      {
+       ItemBox[
+        TagBox[
+         RowBox[{"{", 
           RowBox[{
            RowBox[{
-            RowBox[{
-             RowBox[{"fp", "'"}], "[", "xp", "]"}], " ", 
-            RowBox[{"D", "[", 
-             RowBox[{"Pm", ",", "xp"}], "]"}]}], "-", 
+            RowBox[{"-", "1.350058396849819`"}], " ", "m0"}], ",", 
            RowBox[{
+            FractionBox["1", "2"], " ", 
             RowBox[{"(", 
-             RowBox[{
-              RowBox[{"2", " ", "m"}], "-", "1"}], ")"}], 
-            RowBox[{
-             RowBox[{"fp", "''"}], "[", "xp", "]"}], " ", "Pm"}]}]}], ",", 
-         "1", ",", 
-         RowBox[{"Range", "[", 
-          RowBox[{"n", "-", "1"}], "]"}]}], "]"}], "/.", 
-       RowBox[{
-        RowBox[{
-         RowBox[{
-          RowBox[{"Derivative", "[", "m_", "]"}], "[", "fp", "]"}], "[", "xp",
-          "]"}], ":>", 
-        RowBox[{"dfs", "[", 
-         RowBox[{"[", "m", "]"}], "]"}]}]}]}], ";", "\[IndentingNewLine]", 
-     RowBox[{"MapIndexed", "[", 
-      RowBox[{
-       RowBox[{
+             RowBox[{"0.`", "\[VeryThinSpace]", "+", 
+              FractionBox[
+               RowBox[{"0.9048826940124446`", " ", "m0"}], 
+               RowBox[{"0.`", "\[VeryThinSpace]", "-", 
+                RowBox[{"35.8650482175897`", " ", "h0"}]}]]}], ")"}]}], ",", 
+           RowBox[{
+            FractionBox["1", "6"], " ", 
+            RowBox[{"(", 
+             RowBox[{"0.`", "\[VeryThinSpace]", "+", 
+              FractionBox[
+               RowBox[{"0.9048826940124446`", " ", 
+                RowBox[{"(", 
+                 TemplateBox[{"1"},
+                  "OutputSizeLimit`Skeleton"], ")"}], " ", "m0"}], 
+               SuperscriptBox[
+                RowBox[{"(", 
+                 TemplateBox[{"1"},
+                  "OutputSizeLimit`Skeleton"], ")"}], "3"]], "-", 
+              FractionBox[
+               RowBox[{"24.93588598416555`", " ", "m0"}], 
+               SuperscriptBox[
+                RowBox[{"(", 
+                 TemplateBox[{"1"},
+                  "OutputSizeLimit`Skeleton"], ")"}], "2"]]}], ")"}]}], ",", 
+           TemplateBox[{"5"},
+            "OutputSizeLimit`Skeleton"], ",", 
+           FractionBox[
+            TemplateBox[{"1"},
+             "OutputSizeLimit`Skeleton"], "362880"], ",", 
+           FractionBox[
+            RowBox[{"0.`", "\[VeryThinSpace]", "-", 
+             FractionBox[
+              TemplateBox[{"1"},
+               "OutputSizeLimit`Skeleton"], 
+              SuperscriptBox[
+               TemplateBox[{"1"},
+                "OutputSizeLimit`Skeleton"], "10"]], "+", 
+             TemplateBox[{"10"},
+              "OutputSizeLimit`Skeleton"], "+", 
+             FractionBox[
+              RowBox[{"0.9048826940124446`", " ", 
+               RowBox[{"(", 
+                RowBox[{
+                 TemplateBox[{"1"},
+                  "OutputSizeLimit`Skeleton"], "+", 
+                 TemplateBox[{"1"},
+                  "OutputSizeLimit`Skeleton"]}], ")"}], " ", "m0"}], 
+              SuperscriptBox[
+               RowBox[{"(", 
+                RowBox[{
+                 TemplateBox[{"3"},
+                  "OutputSizeLimit`Skeleton"], "\[VeryThinSpace]", "-", 
+                 TemplateBox[{"1"},
+                  "OutputSizeLimit`Skeleton"]}], ")"}], "17"]]}], "3628800"], 
+           ",", 
+           FractionBox[
+            RowBox[{"0.`", "\[VeryThinSpace]", "+", 
+             TemplateBox[{"14"},
+              "OutputSizeLimit`Skeleton"], "+", 
+             FractionBox[
+              RowBox[{"0.9048826940124446`", " ", 
+               RowBox[{"(", 
+                TemplateBox[{"1"},
+                 "OutputSizeLimit`Skeleton"], ")"}], " ", "m0"}], 
+              SuperscriptBox[
+               RowBox[{"(", 
+                RowBox[{"0.`", "\[VeryThinSpace]", "-", 
+                 RowBox[{"35.8650482175897`", " ", 
+                  TemplateBox[{"2"},
+                   "OutputSizeLimit`Skeleton"]}]}], ")"}], "19"]]}], 
+            "39916800"]}], "}"}],
+         Short[#, 5]& ],
+        BaseStyle->{Deployed -> False},
+        StripOnInput->False]},
+      {GridBox[{
+         {
+          PaneBox[
+           TagBox[
+            TooltipBox[
+             StyleBox[
+              StyleBox[
+               DynamicBox[ToBoxes[
+                 FEPrivate`FrontEndResource[
+                 "FEStrings", "sizeBriefExplanation"], StandardForm],
+                ImageSizeCache->{34.69921875, {2., 5.}}],
+               StripOnInput->False,
+               DynamicUpdating->True,
+               LineSpacing->{1, 2},
+               LineIndent->0,
+               LinebreakAdjustments->{1., 100, 0, 0, 0}], "OSLText",
+              StripOnInput->False],
+             StyleBox[
+              DynamicBox[
+               ToBoxes[
+                FEPrivate`FrontEndResource["FEStrings", "sizeExplanation"], 
+                StandardForm]], DynamicUpdating -> True, LineIndent -> 0, 
+              LinebreakAdjustments -> {1., 100, 0, 0, 0}, 
+              LineSpacing -> {1, 2}, StripOnInput -> False]],
+            Annotation[#, 
+             Style[
+              Dynamic[
+               FEPrivate`FrontEndResource["FEStrings", "sizeExplanation"]], 
+              DynamicUpdating -> True, LineIndent -> 0, 
+              LinebreakAdjustments -> {1., 100, 0, 0, 0}, 
+              LineSpacing -> {1, 2}], "Tooltip"]& ],
+           Alignment->Center,
+           BaselinePosition->Baseline,
+           ImageSize->{Automatic, {25, Full}}], 
+          ButtonBox[
+           PaneSelectorBox[{False->
+            StyleBox[
+             StyleBox[
+              DynamicBox[ToBoxes[
+                FEPrivate`FrontEndResource["FEStrings", "sizeShowLess"], 
+                StandardForm],
+               ImageSizeCache->{32.33203125, {1., 5.}}],
+              StripOnInput->False,
+              DynamicUpdating->True,
+              LineSpacing->{1, 2},
+              LineIndent->0,
+              LinebreakAdjustments->{1., 100, 0, 0, 0}], "OSLControl",
+             StripOnInput->False], True->
+            StyleBox[
+             StyleBox[
+              DynamicBox[ToBoxes[
+                FEPrivate`FrontEndResource["FEStrings", "sizeShowLess"], 
+                StandardForm]],
+              StripOnInput->False,
+              DynamicUpdating->True,
+              LineSpacing->{1, 2},
+              LineIndent->0,
+              LinebreakAdjustments->{1., 100, 0, 0, 0}], "OSLControlActive",
+             StripOnInput->False]}, Dynamic[
+             CurrentValue["MouseOver"]],
+            Alignment->Center,
+            FrameMargins->0,
+            ImageSize->{Automatic, {25, Full}}],
+           Appearance->None,
+           BaselinePosition->Baseline,
+           
+           ButtonFunction:>OutputSizeLimit`ButtonFunction[
+            OutputSizeLimit`Defer, 144, 31653097672855845075, 5/2],
+           Enabled->True,
+           Evaluator->Automatic,
+           Method->"Queued"], 
+          ButtonBox[
+           PaneSelectorBox[{False->
+            StyleBox[
+             StyleBox[
+              DynamicBox[ToBoxes[
+                FEPrivate`FrontEndResource["FEStrings", "sizeShowMore"], 
+                StandardForm],
+               ImageSizeCache->{35.8125, {1., 5.}}],
+              StripOnInput->False,
+              DynamicUpdating->True,
+              LineSpacing->{1, 2},
+              LineIndent->0,
+              LinebreakAdjustments->{1., 100, 0, 0, 0}], "OSLControl",
+             StripOnInput->False], True->
+            StyleBox[
+             StyleBox[
+              DynamicBox[ToBoxes[
+                FEPrivate`FrontEndResource["FEStrings", "sizeShowMore"], 
+                StandardForm]],
+              StripOnInput->False,
+              DynamicUpdating->True,
+              LineSpacing->{1, 2},
+              LineIndent->0,
+              LinebreakAdjustments->{1., 100, 0, 0, 0}], "OSLControlActive",
+             StripOnInput->False]}, Dynamic[
+             CurrentValue["MouseOver"]],
+            Alignment->Center,
+            FrameMargins->0,
+            ImageSize->{Automatic, {25, Full}}],
+           Appearance->None,
+           BaselinePosition->Baseline,
+           
+           ButtonFunction:>OutputSizeLimit`ButtonFunction[
+            OutputSizeLimit`Defer, 144, 31653097672855845075, 5 2],
+           Enabled->True,
+           Evaluator->Automatic,
+           Method->"Queued"], 
+          ButtonBox[
+           PaneSelectorBox[{False->
+            StyleBox[
+             StyleBox[
+              DynamicBox[ToBoxes[
+                FEPrivate`FrontEndResource["FEStrings", "sizeShowAll"], 
+                StandardForm],
+               ImageSizeCache->{26.54296875, {1., 5.}}],
+              StripOnInput->False,
+              DynamicUpdating->True,
+              LineSpacing->{1, 2},
+              LineIndent->0,
+              LinebreakAdjustments->{1., 100, 0, 0, 0}], "OSLControl",
+             StripOnInput->False], True->
+            StyleBox[
+             StyleBox[
+              DynamicBox[ToBoxes[
+                FEPrivate`FrontEndResource["FEStrings", "sizeShowAll"], 
+                StandardForm],
+               ImageSizeCache->{26.54296875, {1., 5.}}],
+              StripOnInput->False,
+              DynamicUpdating->True,
+              LineSpacing->{1, 2},
+              LineIndent->0,
+              LinebreakAdjustments->{1., 100, 0, 0, 0}], "OSLControlActive",
+             StripOnInput->False]}, Dynamic[
+             CurrentValue["MouseOver"]],
+            Alignment->Center,
+            FrameMargins->0,
+            ImageSize->{Automatic, {25, Full}}],
+           Appearance->None,
+           BaselinePosition->Baseline,
+           
+           ButtonFunction:>OutputSizeLimit`ButtonFunction[
+            OutputSizeLimit`Defer, 144, 31653097672855845075, Infinity],
+           Enabled->True,
+           Evaluator->Automatic,
+           Method->"Queued"], 
+          ButtonBox[
+           PaneSelectorBox[{False->
+            StyleBox[
+             StyleBox[
+              DynamicBox[ToBoxes[
+                FEPrivate`FrontEndResource["FEStrings", "sizeChangeLimit"], 
+                StandardForm],
+               ImageSizeCache->{45.6328125, {1., 5.}}],
+              StripOnInput->False,
+              DynamicUpdating->True,
+              LineSpacing->{1, 2},
+              LineIndent->0,
+              LinebreakAdjustments->{1., 100, 0, 0, 0}], "OSLControl",
+             StripOnInput->False], True->
+            StyleBox[
+             StyleBox[
+              DynamicBox[ToBoxes[
+                FEPrivate`FrontEndResource["FEStrings", "sizeChangeLimit"], 
+                StandardForm],
+               ImageSizeCache->{45.6328125, {1., 5.}}],
+              StripOnInput->False,
+              DynamicUpdating->True,
+              LineSpacing->{1, 2},
+              LineIndent->0,
+              LinebreakAdjustments->{1., 100, 0, 0, 0}], "OSLControlActive",
+             StripOnInput->False]}, Dynamic[
+             CurrentValue["MouseOver"]],
+            Alignment->Center,
+            FrameMargins->0,
+            ImageSize->{Automatic, {25, Full}}],
+           Appearance->None,
+           BaselinePosition->Baseline,
+           ButtonFunction:>FrontEndExecute[{
+              FrontEnd`SetOptions[
+              FrontEnd`$FrontEnd, 
+               FrontEnd`PreferencesSettings -> {"Page" -> "Advanced"}], 
+              FrontEnd`FrontEndToken["PreferencesDialog"]}],
+           Evaluator->None,
+           Method->"Preemptive"]}
+        },
+        AutoDelete->False,
+        FrameStyle->GrayLevel[0.85],
+        GridBoxDividers->{"Columns" -> {False, {True}}},
+        GridBoxItemSize->{"Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}},
+        GridBoxSpacings->{"Columns" -> {{2}}}]}
+     },
+     DefaultBaseStyle->"Column",
+     GridBoxAlignment->{"Columns" -> {{Left}}, "Rows" -> {{Baseline}}},
+     GridBoxDividers->{"Columns" -> {{False}}, "Rows" -> {{False}}},
+     GridBoxItemSize->{"Columns" -> {{Automatic}}, "Rows" -> {{1.}}},
+     GridBoxSpacings->{"Columns" -> {
+         Offset[0.27999999999999997`], {
+          Offset[0.5599999999999999]}, 
+         Offset[0.27999999999999997`]}, "Rows" -> {
+         Offset[0.2], 
+         Offset[1.2], {
+          Offset[0.4]}, 
+         Offset[0.2]}}],
+    BaseStyle->"OutputSizeLimit",
+    FrameMargins->{{12, 12}, {0, 15}},
+    FrameStyle->GrayLevel[0.85],
+    RoundingRadius->5,
+    StripOnInput->False],
+   Deploy,
+   DefaultBaseStyle->"Deploy"],
+  If[31653097672855845075 === $SessionID, 
+   Out[144], Message[
+     MessageName[Syntax, "noinfoker"]]; Missing["NotAvailable"]; 
+   Null]]], "Output",
+ CellChangeTimes->{
+  3.83698713111462*^9, 3.836987321633688*^9, 3.8369882151463833`*^9, 
+   3.83698833312978*^9, 3.837070043781701*^9, 3.837070402550199*^9, 
+   3.837072865049255*^9, 3.84380209279955*^9, 3.843802127647296*^9, {
+   3.8438022397282753`*^9, 3.843802269116067*^9}, 3.8438023757788677`*^9, {
+   3.84380243891853*^9, 3.843802516429851*^9}, 3.8438026852881813`*^9, 
+   3.843802807050653*^9, 3.843802949596843*^9, 3.843803004201252*^9, 
+   3.843803045143682*^9},
+ CellLabel->
+  "Out[144]=",ExpressionUUID->"242f95ed-9db7-4d20-804e-da19445983d7"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+ RowBox[{"m0rule", "=", 
+  RowBox[{
+   RowBox[{"Solve", "[", 
+    RowBox[{
+     RowBox[{
+      RowBox[{"datCas", "[", 
+       RowBox[{"[", 
         RowBox[{"{", 
-         RowBox[{"Pn", ",", "i"}], "}"}], "\[Function]", 
-        RowBox[{"Pn", "/", 
-         SuperscriptBox[
-          RowBox[{"dfs", "[", 
-           RowBox[{"[", "1", "]"}], "]"}], 
+         RowBox[{"1", ",", "2"}], "}"}], "]"}], "]"}], "==", 
+      RowBox[{
+       RowBox[{"{", 
+        RowBox[{"1", ",", "1"}], "}"}], 
+       RowBox[{"N", "[", 
+        RowBox[{"datNum", "[", 
+         RowBox[{"[", 
           RowBox[{
-           RowBox[{"2", 
-            RowBox[{"i", "[", 
-             RowBox[{"[", "1", "]"}], "]"}]}], "-", "1"}]]}]}], ",", "Pns"}], 
-      "]"}]}]}], "\[IndentingNewLine]", "]"}]}]], "Input",
- CellChangeTimes->{{3.83698298732085*^9, 3.836983017577533*^9}, {
-  3.8369830656105223`*^9, 3.836983192220652*^9}, {3.8369852957954073`*^9, 
-  3.836985418901877*^9}, {3.836985491119132*^9, 3.836985540335959*^9}, {
-  3.836985603761368*^9, 3.8369856038652277`*^9}, {3.836986669957045*^9, 
-  3.836986682421062*^9}},
- CellLabel->"In[17]:=",ExpressionUUID->"275ca4b6-77bf-4727-923f-67c9c37ad46e"],
+           RowBox[{"2", ";;", "3"}], ",", "2"}], "]"}], "]"}], "]"}]}]}], ",", 
+     RowBox[{"{", 
+      RowBox[{"m0", ",", "h0"}], "}"}]}], "]"}], "[", 
+   RowBox[{"[", "1", "]"}], "]"}]}]], "Input",
+ CellChangeTimes->{{3.843802521611643*^9, 3.843802522562986*^9}, {
+  3.843802635013496*^9, 3.843802651093405*^9}, {3.843802690006391*^9, 
+  3.843802701983755*^9}, {3.843802816200614*^9, 3.843802837257374*^9}, {
+  3.8438028734585047`*^9, 3.8438028900667257`*^9}, {3.843802924195238*^9, 
+  3.843802955770994*^9}, {3.843803008140307*^9, 3.8438030090357723`*^9}},
+ CellLabel->
+  "In[145]:=",ExpressionUUID->"6610225e-a69c-4a81-b8a1-799f1847532b"],
+
+Cell[BoxData[
+ TemplateBox[{
+  "Solve", "ratnz", 
+   "\"Solve was unable to solve the system with inexact coefficients. The \
+answer was obtained by solving a corresponding exact system and numericizing \
+the result.\"", 2, 145, 49, 31653097672855845075, "Local"},
+  "MessageTemplate"]], "Message", "MSG",
+ CellChangeTimes->{{3.8438025179787197`*^9, 3.843802522855481*^9}, {
+   3.843802638378536*^9, 3.843802651383214*^9}, {3.843802686773093*^9, 
+   3.843802702490966*^9}, {3.843802817095837*^9, 3.843802837439028*^9}, 
+   3.843802890250163*^9, {3.843802926887959*^9, 3.8438029560400467`*^9}, {
+   3.843803005409687*^9, 3.8438030093603773`*^9}, 3.843803045987146*^9},
+ CellLabel->
+  "During evaluation of \
+In[145]:=",ExpressionUUID->"127b9b13-a14c-4d18-9b15-ca5cab6104b0"],
+
+Cell[BoxData[
+ RowBox[{"{", 
+  RowBox[{
+   RowBox[{"m0", "\[Rule]", "1.0057626728406195`"}], ",", 
+   RowBox[{"h0", "\[Rule]", "0.25918179834063165`"}]}], "}"}]], "Output",
+ CellChangeTimes->{{3.843802517987623*^9, 3.8438025228598337`*^9}, {
+   3.843802638382826*^9, 3.8438026513875303`*^9}, {3.843802686777157*^9, 
+   3.843802702496963*^9}, {3.843802817105462*^9, 3.8438028374470243`*^9}, 
+   3.843802890260683*^9, {3.843802926893958*^9, 3.843802956046287*^9}, {
+   3.8438030054153147`*^9, 3.843803009364552*^9}, 3.84380304599266*^9},
+ CellLabel->
+  "Out[145]=",ExpressionUUID->"6f1d88a3-a882-48c5-84c6-fac1de7856da"]
+}, Open  ]],
 
 Cell[CellGroupData[{
 
+Cell[BoxData[
+ RowBox[{"datCas", "/.", "m0rule"}]], "Input",
+ CellChangeTimes->{{3.843802718839593*^9, 3.843802722902604*^9}},
+ CellLabel->
+  "In[146]:=",ExpressionUUID->"1b205c4d-ffa7-4114-8899-bfcb09a7888c"],
+
+Cell[BoxData[
+ RowBox[{"{", 
+  RowBox[{
+   RowBox[{"-", "1.3578383417065958`"}], ",", 
+   RowBox[{"-", "0.0489532897202746`"}], ",", "0.03583699905594043`", ",", 
+   RowBox[{"-", "0.05314325466674752`"}], ",", "0.10995313306417173`", ",", 
+   RowBox[{"-", "0.27695711938473166`"}], ",", "0.7945288687358903`", ",", 
+   RowBox[{"-", "2.4997021825853816`"}], ",", "8.42446931003965`", ",", 
+   RowBox[{"-", "29.943251124061156`"}], ",", "111.02920850534656`"}], 
+  "}"}]], "Output",
+ CellChangeTimes->{
+  3.843802723612238*^9, {3.843802934969544*^9, 3.843802956788475*^9}, 
+   3.8438030099091682`*^9, 3.843803046518009*^9},
+ CellLabel->
+  "Out[146]=",ExpressionUUID->"9d75172f-bdf7-4272-a390-90c81021094e"]
+}, Open  ]],
+
 Cell[BoxData[{
  RowBox[{
-  RowBox[{"totCas", "=", "17"}], ";"}], "\[IndentingNewLine]", 
+  RowBox[{"totCas", "=", "10"}], ";"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"ds", "=", 
    RowBox[{
     RowBox[{
-     RowBox[{"invDerivativeList", "[", "totCas", "]"}], "[", 
+     RowBox[{"invDerivativeList", "[", 
+      RowBox[{"totCas", "+", "1"}], "]"}], "[", 
      RowBox[{"\[Theta]", "\[Function]", 
       FractionBox[
        RowBox[{"hCas", "[", "\[Theta]", "]"}], 
        SuperscriptBox[
         RowBox[{"(", 
-         RowBox[{"t", "[", "\[Theta]", "]"}], ")"}], 
-        RowBox[{"\[CapitalDelta]", "[", "2", "]"}]]]}], "]"}], "[", 
-    "\[Theta]0Cas", "]"}]}], ";"}], "\[IndentingNewLine]", 
+         RowBox[{"-", 
+          RowBox[{"t", "[", "\[Theta]", "]"}]}], ")"}], 
+        RowBox[{"\[CapitalDelta]", "[", "2", "]"}]]]}], "]"}], "[", "0", 
+    "]"}]}], ";"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"dMs", "=", 
    RowBox[{
     RowBox[{"m0", " ", 
-     RowBox[{"Rest", "[", 
-      RowBox[{"NestList", "[", 
+     RowBox[{"NestList", "[", 
+      RowBox[{
        RowBox[{
+        RowBox[{"D", "[", 
+         RowBox[{"#", ",", "\[Theta]"}], "]"}], "&"}], ",", 
+       FractionBox["\[Theta]", 
+        SuperscriptBox[
+         RowBox[{"(", 
+          RowBox[{"1", "-", 
+           SuperscriptBox["\[Theta]", "2"]}], ")"}], 
+         RowBox[{"1", "/", "8"}]]], ",", "totCas"}], "]"}]}], "/.", 
+    RowBox[{"\[Theta]", "->", "0"}]}]}], ";"}], "\[IndentingNewLine]", 
+ RowBox[{
+  RowBox[{"datCasH", "=", 
+   RowBox[{"Table", "[", 
+    RowBox[{
+     RowBox[{
+      RowBox[{"-", 
+       RowBox[{"Sum", "[", 
         RowBox[{
-         RowBox[{"D", "[", 
-          RowBox[{"#", ",", "\[Theta]"}], "]"}], "&"}], ",", 
-        FractionBox["\[Theta]", 
-         SuperscriptBox[
-          RowBox[{"(", 
-           RowBox[{
-            SuperscriptBox["\[Theta]", "2"], "-", "1"}], ")"}], 
-          RowBox[{"1", "/", "8"}]]], ",", "totCas"}], "]"}], "]"}]}], "/.", 
-    RowBox[{"\[Theta]", "->", "\[Theta]0Cas"}]}]}], 
-  ";"}], "\[IndentingNewLine]", 
- RowBox[{"datCas", "=", 
-  RowBox[{"Table", "[", 
+         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.843803094421678*^9, 3.843803146878536*^9}},
+ CellLabel->
+  "In[162]:=",ExpressionUUID->"34fe8bc4-abfb-4504-a7f1-045cd2a55ee6"],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+ RowBox[{"datCasH", "/.", "m0rule"}]], "Input",
+ CellChangeTimes->{{3.8438031226868477`*^9, 3.843803125309888*^9}},
+ CellLabel->
+  "In[166]:=",ExpressionUUID->"40ec85df-7487-4ffe-afab-49583a090315"],
+
+Cell[BoxData[
+ RowBox[{"{", 
+  RowBox[{"0", ",", 
+   RowBox[{"-", "1.9402648628874553`"}], ",", "0.`", ",", 
+   "9.712877740891042`", ",", "0.`", ",", 
+   RowBox[{"-", "122.94298442688967`"}], ",", "0.`", ",", 
+   "2087.0403258684037`", ",", "0.`", ",", 
+   RowBox[{"-", "41080.701036180886`"}], ",", "0.`"}], "}"}]], "Output",
+ CellChangeTimes->{{3.843803125487954*^9, 3.843803147861979*^9}},
+ CellLabel->
+  "Out[166]=",ExpressionUUID->"f1acdb4d-a707-469d-887b-abc0e24714e7"]
+}, Open  ]],
+
+Cell[BoxData[{
+ RowBox[{
+  RowBox[{"totCas", "=", "10"}], ";"}], "\[IndentingNewLine]", 
+ RowBox[{
+  RowBox[{"ds", "=", 
    RowBox[{
     RowBox[{
-     RowBox[{"Sum", "[", 
+     RowBox[{"invDerivativeList", "[", 
+      RowBox[{"totCas", "+", "1"}], "]"}], "[", 
+     RowBox[{"\[Theta]", "\[Function]", 
+      FractionBox[
+       RowBox[{"t", "[", "\[Theta]", "]"}], 
+       SuperscriptBox[
+        RowBox[{"(", 
+         RowBox[{"hCas", "[", "\[Theta]", "]"}], ")"}], 
+        RowBox[{"1", "/", 
+         RowBox[{"\[CapitalDelta]", "[", "2", "]"}]}]]]}], "]"}], "[", "1", 
+    "]"}]}], ";"}], "\[IndentingNewLine]", 
+ RowBox[{
+  RowBox[{"dMs", "=", 
+   RowBox[{
+    RowBox[{"m0", " ", 
+     RowBox[{"NestList", "[", 
       RowBox[{
        RowBox[{
-        RowBox[{"dMs", "[", 
-         RowBox[{"[", "k", "]"}], "]"}], 
-        RowBox[{"BellY", "[", 
-         RowBox[{"n", ",", "k", ",", 
-          RowBox[{"ds", "[", 
-           RowBox[{"[", 
-            RowBox[{";;", 
-             RowBox[{"n", "-", "k", "+", "1"}]}], "]"}], "]"}]}], "]"}]}], 
-       ",", 
-       RowBox[{"{", 
-        RowBox[{"k", ",", "1", ",", "n"}], "}"}]}], "]"}], "/", 
+        RowBox[{"D", "[", 
+         RowBox[{"#", ",", "\[Theta]"}], "]"}], "&"}], ",", 
+       FractionBox["\[Theta]", 
+        SuperscriptBox[
+         RowBox[{"hCas", "[", "\[Theta]", "]"}], 
+         RowBox[{"1", "/", "15"}]]], ",", "totCas"}], "]"}]}], "/.", 
+    RowBox[{"\[Theta]", "->", "1"}]}]}], ";"}], "\[IndentingNewLine]", 
+ RowBox[{
+  RowBox[{"datCasM", "=", 
+   RowBox[{"Table", "[", 
+    RowBox[{
      RowBox[{
-      RowBox[{"(", 
-       RowBox[{"n", "+", "1"}], ")"}], "!"}]}], ",", 
-    RowBox[{"{", 
-     RowBox[{"n", ",", "1", ",", "totCas"}], "}"}]}], "]"}]}]}], "Input",
- CellChangeTimes->CompressedData["
-1:eJxTTMoPSmViYGCQAGIQHXVsxuTiaW8cpepuTAPRt/4wzgbRDhmJB0C0AfeT
-IyDaefWB42B5s2LWciBd8vg9F4iWm+0rAKIn5SyRBNFNksv1QPSRyuv6IDry
-Dp8DiOYzEQbTj2QUA0D0ylfKYPqeomIImJ6UGQqiy/K0UkC0yGE9MN2VltcE
-ons2t7aAaOWou2dAtHDoXzBdsUlNsQJI+3A4KINolltbdEC0jKGiLYhOmRju
-D6J7Ns4NBNE2EZ7nQHS3uDeYFhCceBtEm5TnPATR3usq0yqBdF1aP5iuqduT
-B6L7Zx0E0yq1U65UAek9VV3XQPSOV4nhC4G0jsR8MJ1QPqcKRGc0zQfTXDKC
-p0D0uX4ZMF0eMa1nEZB+eng+mAYA2XPFWw==
-  "],
- CellLabel->"In[18]:=",ExpressionUUID->"9e250256-8cd3-402a-bc21-3ee002bb1f82"],
+      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.843803094421678*^9, 3.843803146878536*^9}, {3.843803198639724*^9, 
+   3.84380324222451*^9}, {3.8438033353723392`*^9, 3.843803345882052*^9}, {
+   3.8438034178358307`*^9, 3.843803431043694*^9}, {3.84380358438281*^9, 
+   3.843803587446344*^9}, 3.843804356868655*^9},
+ CellLabel->
+  "In[191]:=",ExpressionUUID->"b7de9e01-4fea-49ea-97a1-b154bc61886b"],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+ RowBox[{"datCasM", "/.", "m0rule"}]], "Input",
+ CellChangeTimes->{{3.843803230256878*^9, 3.8438032333117857`*^9}},
+ CellLabel->
+  "In[195]:=",ExpressionUUID->"40438766-386e-4ef0-95c5-bb36288917ef"],
 
 Cell[BoxData[
  RowBox[{"{", 
   RowBox[{
-   RowBox[{"-", "0.04895328972000001`"}], ",", "0.0358369990555384`", ",", 
-   RowBox[{"-", "0.05314325466585322`"}], ",", "0.10995313306170523`", ",", 
-   RowBox[{"-", "0.27695711937696654`"}], ",", "0.7945288687091524`", ",", 
-   RowBox[{"-", "2.499702182487203`"}], ",", "8.424469309661326`", ",", 
-   RowBox[{"-", "29.943251122547245`"}], ",", "111.02920849909947`", ",", 
-   RowBox[{"-", "426.1281183859369`"}], ",", "1682.8930954771286`", ",", 
-   RowBox[{"-", "6808.202478556653`"}], ",", "28115.088082925722`", ",", 
-   RowBox[{"-", "118184.25755733704`"}], ",", "504553.6128947737`", ",", 
-   RowBox[{"-", "2.183603153706728`*^6"}]}], "}"}]], "Output",
- CellChangeTimes->{3.83698713111462*^9, 3.836987321633688*^9, 
-  3.8369882151463833`*^9, 3.83698833312978*^9, 3.837070043781701*^9, 
-  3.837070402550199*^9, 3.837072865049255*^9},
- CellLabel->"Out[21]=",ExpressionUUID->"1690ea78-e2c3-4e83-a9e5-8be735e0bca3"]
+   RowBox[{"-", "1.2803840363806758`"}], ",", 
+   RowBox[{"-", "0.08312876211833689`"}], ",", "0.013329954661202004`", ",", 
+   RowBox[{"-", "0.0019834846485221125`"}], ",", "0.00005580926915319621`", 
+   ",", "0.00011134671023551035`", ",", 
+   RowBox[{"-", "0.00005191912354628053`"}], ",", "0.0000129868232839969`", 
+   ",", 
+   RowBox[{"-", "7.456745384784115`*^-7"}], ",", 
+   RowBox[{"-", "1.0845467513314156`*^-6"}], ",", "6.211171366695239`*^-7"}], 
+  "}"}]], "Output",
+ CellChangeTimes->{{3.843803233497788*^9, 3.843803242894781*^9}, 
+   3.84380334832369*^9, 3.843803432165793*^9, 3.843803588716037*^9},
+ CellLabel->
+  "Out[195]=",ExpressionUUID->"462d7199-fbaa-4949-8313-776399b35c4a"]
 }, Open  ]],
 
 Cell[CellGroupData[{
@@ -2336,9 +2878,9 @@ Cell[BoxData[
  CellLabel->"Out[91]=",ExpressionUUID->"b48c9c3b-23cc-4864-8ea2-6e98855eb77d"]
 }, Open  ]]
 },
-WindowSize->{891.75, 1060.5},
-WindowMargins->{{2.25, Automatic}, {2.25, Automatic}},
-FrontEndVersion->"12.3 for Linux x86 (64-bit) (May 11, 2021)",
+WindowSize->{635.25, 1060.5},
+WindowMargins->{{642, Automatic}, {2.25, Automatic}},
+FrontEndVersion->"12.3 for Linux x86 (64-bit) (July 9, 2021)",
 StyleDefinitions->"Default.nb",
 ExpressionUUID->"d257fbda-2be1-4221-9ac8-5153716cfc00"
 ]
@@ -2356,77 +2898,95 @@ Notebook[{
 Cell[558, 20, 596, 10, 24, "Input",ExpressionUUID->"94557b52-cc86-4010-b22d-6b74ac4aad03"],
 Cell[1157, 32, 7539, 107, 55, "Input",ExpressionUUID->"872503c3-0779-4b87-805f-13a088f459e3"],
 Cell[8699, 141, 1256, 31, 105, "Input",ExpressionUUID->"02a1f6e5-4144-4552-9463-2a226c1ee186"],
-Cell[9958, 174, 2710, 69, 121, "Input",ExpressionUUID->"d8831e63-87f3-498c-8951-21f347e539ff"],
-Cell[12671, 245, 941, 24, 59, "Input",ExpressionUUID->"04b6d9c8-662b-43f2-955f-87d9d6c7b622"],
-Cell[13615, 271, 915, 20, 40, "Input",ExpressionUUID->"bd6b698e-6871-43fa-9d98-7385cc594c9c"],
-Cell[14533, 293, 752, 17, 38, "Input",ExpressionUUID->"7bddfed2-a761-4c62-bcdc-0460988fbd67"],
-Cell[15288, 312, 494, 12, 22, "Input",ExpressionUUID->"ec8b9383-61d6-44c2-b776-3d651fc19660"],
+Cell[9958, 174, 2710, 69, 185, "Input",ExpressionUUID->"d8831e63-87f3-498c-8951-21f347e539ff"],
+Cell[12671, 245, 942, 24, 59, "Input",ExpressionUUID->"04b6d9c8-662b-43f2-955f-87d9d6c7b622"],
+Cell[13616, 271, 916, 20, 40, "Input",ExpressionUUID->"bd6b698e-6871-43fa-9d98-7385cc594c9c"],
+Cell[14535, 293, 752, 17, 38, "Input",ExpressionUUID->"7bddfed2-a761-4c62-bcdc-0460988fbd67"],
+Cell[15290, 312, 618, 15, 22, "Input",ExpressionUUID->"ec8b9383-61d6-44c2-b776-3d651fc19660"],
+Cell[CellGroupData[{
+Cell[15933, 331, 943, 23, 59, "Input",ExpressionUUID->"cbeb62aa-7e89-4ec5-a1cb-d5012bac06fd"],
+Cell[16879, 356, 334, 5, 25, "Output",ExpressionUUID->"9ddcc7c6-de02-4d31-9401-bbbc5eefe582"]
+}, Open  ]],
+Cell[CellGroupData[{
+Cell[17250, 366, 1041, 24, 22, "Input",ExpressionUUID->"7efd85be-fd7d-435a-a870-214d8ec4025d"],
+Cell[18294, 392, 924, 17, 34, "Message",ExpressionUUID->"4bba24d3-fa02-4ad0-86b9-3bdb7e46a000"],
+Cell[19221, 411, 877, 15, 25, "Output",ExpressionUUID->"e2f0010a-e1b6-4537-a8ae-b59af3108330"]
+}, Open  ]],
+Cell[CellGroupData[{
+Cell[20135, 431, 523, 12, 55, "Input",ExpressionUUID->"28731a52-d222-49b4-aebb-262f9a5d4dd5"],
+Cell[20661, 445, 240, 3, 25, "Output",ExpressionUUID->"65626e49-d3af-4115-b5d3-96f63b95604e"],
+Cell[20904, 450, 257, 4, 25, "Output",ExpressionUUID->"fe83a208-ae15-4c69-a7c9-b98329954997"]
+}, Open  ]],
+Cell[CellGroupData[{
+Cell[21198, 459, 531, 11, 24, "Input",ExpressionUUID->"f42f7737-f59e-4289-8153-b78626e30f0e"],
+Cell[21732, 472, 583, 11, 25, "Output",ExpressionUUID->"c7c2e58e-e6d9-436b-9cc9-4238f4205a5a"]
+}, Open  ]],
 Cell[CellGroupData[{
-Cell[15807, 328, 943, 23, 41, "Input",ExpressionUUID->"cbeb62aa-7e89-4ec5-a1cb-d5012bac06fd"],
-Cell[16753, 353, 334, 5, 25, "Output",ExpressionUUID->"9ddcc7c6-de02-4d31-9401-bbbc5eefe582"]
+Cell[22352, 488, 3523, 81, 109, "Input",ExpressionUUID->"9e250256-8cd3-402a-bc21-3ee002bb1f82"],
+Cell[25878, 571, 12378, 312, 144, "Output",ExpressionUUID->"242f95ed-9db7-4d20-804e-da19445983d7"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[17124, 363, 1038, 25, 22, "Input",ExpressionUUID->"7efd85be-fd7d-435a-a870-214d8ec4025d"],
-Cell[18165, 390, 825, 16, 22, "Message",ExpressionUUID->"2e43731b-67a2-4964-9c03-694bb26aa975"],
-Cell[18993, 408, 1025, 22, 26, "Output",ExpressionUUID->"325e2b46-9d64-4724-8091-36718f8f11a2"]
+Cell[38293, 888, 1058, 27, 22, "Input",ExpressionUUID->"6610225e-a69c-4a81-b8a1-799f1847532b"],
+Cell[39354, 917, 773, 14, 34, "Message",ExpressionUUID->"127b9b13-a14c-4d18-9b15-ca5cab6104b0"],
+Cell[40130, 933, 618, 11, 25, "Output",ExpressionUUID->"6f1d88a3-a882-48c5-84c6-fac1de7856da"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[20055, 435, 523, 12, 55, "Input",ExpressionUUID->"28731a52-d222-49b4-aebb-262f9a5d4dd5"],
-Cell[20581, 449, 240, 3, 25, "Output",ExpressionUUID->"65626e49-d3af-4115-b5d3-96f63b95604e"],
-Cell[20824, 454, 257, 4, 25, "Output",ExpressionUUID->"fe83a208-ae15-4c69-a7c9-b98329954997"]
+Cell[40785, 949, 209, 4, 22, "Input",ExpressionUUID->"1b205c4d-ffa7-4114-8899-bfcb09a7888c"],
+Cell[40997, 955, 704, 14, 42, "Output",ExpressionUUID->"9d75172f-bdf7-4272-a390-90c81021094e"]
 }, Open  ]],
+Cell[41716, 972, 3620, 84, 125, "Input",ExpressionUUID->"34fe8bc4-abfb-4504-a7f1-045cd2a55ee6"],
 Cell[CellGroupData[{
-Cell[21118, 463, 531, 11, 24, "Input",ExpressionUUID->"f42f7737-f59e-4289-8153-b78626e30f0e"],
-Cell[21652, 476, 583, 11, 25, "Output",ExpressionUUID->"c7c2e58e-e6d9-436b-9cc9-4238f4205a5a"]
+Cell[45361, 1060, 212, 4, 22, "Input",ExpressionUUID->"40ec85df-7487-4ffe-afab-49583a090315"],
+Cell[45576, 1066, 476, 10, 25, "Output",ExpressionUUID->"f1acdb4d-a707-469d-887b-abc0e24714e7"]
 }, Open  ]],
-Cell[22250, 490, 2575, 70, 108, "Input",ExpressionUUID->"275ca4b6-77bf-4727-923f-67c9c37ad46e"],
+Cell[46067, 1079, 3789, 85, 122, "Input",ExpressionUUID->"b7de9e01-4fea-49ea-97a1-b154bc61886b"],
 Cell[CellGroupData[{
-Cell[24850, 564, 2358, 66, 109, "Input",ExpressionUUID->"9e250256-8cd3-402a-bc21-3ee002bb1f82"],
-Cell[27211, 632, 960, 15, 44, "Output",ExpressionUUID->"1690ea78-e2c3-4e83-a9e5-8be735e0bca3"]
+Cell[49881, 1168, 212, 4, 22, "Input",ExpressionUUID->"40438766-386e-4ef0-95c5-bb36288917ef"],
+Cell[50096, 1174, 735, 15, 44, "Output",ExpressionUUID->"462d7199-fbaa-4949-8313-776399b35c4a"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[28208, 652, 4388, 128, 156, "Input",ExpressionUUID->"f81ed9c7-d940-4314-963c-c36ea0215ef3"],
-Cell[32599, 782, 1199, 20, 44, "Output",ExpressionUUID->"2234af5e-068b-46a3-beaf-99a907ef5184"]
+Cell[50868, 1194, 4388, 128, 175, "Input",ExpressionUUID->"f81ed9c7-d940-4314-963c-c36ea0215ef3"],
+Cell[55259, 1324, 1199, 20, 60, "Output",ExpressionUUID->"2234af5e-068b-46a3-beaf-99a907ef5184"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[33835, 807, 882, 21, 41, "Input",ExpressionUUID->"dd088f15-420b-4212-8186-b62091ac9d6c"],
-Cell[34720, 830, 435, 6, 25, "Output",ExpressionUUID->"21b5cf8e-7f0a-45f7-ab93-9e8ebe789a07"]
+Cell[56495, 1349, 882, 21, 41, "Input",ExpressionUUID->"dd088f15-420b-4212-8186-b62091ac9d6c"],
+Cell[57380, 1372, 435, 6, 25, "Output",ExpressionUUID->"21b5cf8e-7f0a-45f7-ab93-9e8ebe789a07"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[35192, 841, 1195, 29, 41, "Input",ExpressionUUID->"e7fdb778-fd9a-4c13-95dc-aefbf5b8a29f"],
-Cell[36390, 872, 349, 5, 25, "Output",ExpressionUUID->"a5e80f1a-5f57-4a7a-ad79-17c92282d141"]
+Cell[57852, 1383, 1195, 29, 41, "Input",ExpressionUUID->"e7fdb778-fd9a-4c13-95dc-aefbf5b8a29f"],
+Cell[59050, 1414, 349, 5, 25, "Output",ExpressionUUID->"a5e80f1a-5f57-4a7a-ad79-17c92282d141"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[36776, 882, 228, 5, 22, "Input",ExpressionUUID->"feb08355-9ad1-4257-9f97-e394177d7e32"],
-Cell[37007, 889, 274, 8, 25, "Output",ExpressionUUID->"87e29f54-9d68-4261-96d4-95d24dcdc36a"]
+Cell[59436, 1424, 228, 5, 22, "Input",ExpressionUUID->"feb08355-9ad1-4257-9f97-e394177d7e32"],
+Cell[59667, 1431, 274, 8, 25, "Output",ExpressionUUID->"87e29f54-9d68-4261-96d4-95d24dcdc36a"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[37318, 902, 1542, 27, 55, "Input",ExpressionUUID->"e61ed118-a102-4290-8440-cba17b87145b"],
-Cell[38863, 931, 12321, 270, 180, "Output",ExpressionUUID->"2e45389f-59c7-4eaf-9b3a-270214fa013e"]
+Cell[59978, 1444, 1542, 27, 55, "Input",ExpressionUUID->"e61ed118-a102-4290-8440-cba17b87145b"],
+Cell[61523, 1473, 12321, 270, 180, "Output",ExpressionUUID->"2e45389f-59c7-4eaf-9b3a-270214fa013e"]
 }, Open  ]],
-Cell[51199, 1204, 175, 3, 22, "Input",ExpressionUUID->"0f87aa05-2cc9-472a-ab27-39ee15777442"],
-Cell[51377, 1209, 356, 7, 24, "Input",ExpressionUUID->"627c17e0-e2d8-4347-8fe9-812183a42941"],
+Cell[73859, 1746, 175, 3, 22, "Input",ExpressionUUID->"0f87aa05-2cc9-472a-ab27-39ee15777442"],
+Cell[74037, 1751, 356, 7, 24, "Input",ExpressionUUID->"627c17e0-e2d8-4347-8fe9-812183a42941"],
 Cell[CellGroupData[{
-Cell[51758, 1220, 176, 2, 22, "Input",ExpressionUUID->"f32c32be-7fcb-4660-bf43-cbc119acda96"],
-Cell[51937, 1224, 2097, 65, 61, "Output",ExpressionUUID->"01b48fb7-d8c9-4012-92f6-488cf7e7087c"]
+Cell[74418, 1762, 176, 2, 22, "Input",ExpressionUUID->"f32c32be-7fcb-4660-bf43-cbc119acda96"],
+Cell[74597, 1766, 2097, 65, 77, "Output",ExpressionUUID->"01b48fb7-d8c9-4012-92f6-488cf7e7087c"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[54071, 1294, 860, 24, 24, "Input",ExpressionUUID->"bfa6a31f-5eb9-4063-b798-e45ab1e3b53a"],
-Cell[54934, 1320, 1944, 57, 60, "Output",ExpressionUUID->"0ca79efb-1867-4b1f-96fc-4a691d54b72e"]
+Cell[76731, 1836, 860, 24, 24, "Input",ExpressionUUID->"bfa6a31f-5eb9-4063-b798-e45ab1e3b53a"],
+Cell[77594, 1862, 1944, 57, 77, "Output",ExpressionUUID->"0ca79efb-1867-4b1f-96fc-4a691d54b72e"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[56915, 1382, 2622, 77, 93, "Input",ExpressionUUID->"242aacdf-e4c7-423c-a26a-12188ba2476e"],
-Cell[59540, 1461, 21838, 439, 180, "Output",ExpressionUUID->"073d831e-11a6-4203-90ce-0df429c1a3f0"]
+Cell[79575, 1924, 2622, 77, 93, "Input",ExpressionUUID->"242aacdf-e4c7-423c-a26a-12188ba2476e"],
+Cell[82200, 2003, 21838, 439, 180, "Output",ExpressionUUID->"073d831e-11a6-4203-90ce-0df429c1a3f0"]
 }, Open  ]],
-Cell[81393, 1903, 725, 15, 24, "Input",ExpressionUUID->"80c6c675-eade-40e1-8d0f-13a6578b30ff"],
-Cell[82121, 1920, 764, 21, 22, "Input",ExpressionUUID->"d8847f2e-21fe-4fb7-84ae-676f83acdf5a"],
+Cell[104053, 2445, 725, 15, 24, "Input",ExpressionUUID->"80c6c675-eade-40e1-8d0f-13a6578b30ff"],
+Cell[104781, 2462, 764, 21, 22, "Input",ExpressionUUID->"d8847f2e-21fe-4fb7-84ae-676f83acdf5a"],
 Cell[CellGroupData[{
-Cell[82910, 1945, 679, 17, 40, "Input",ExpressionUUID->"15d60f60-2f87-40da-a2c6-66d640111d75"],
-Cell[83592, 1964, 1298, 32, 44, "Output",ExpressionUUID->"6f1d902a-0431-4a99-adc2-97eb3c423317"]
+Cell[105570, 2487, 679, 17, 40, "Input",ExpressionUUID->"15d60f60-2f87-40da-a2c6-66d640111d75"],
+Cell[106252, 2506, 1298, 32, 44, "Output",ExpressionUUID->"6f1d902a-0431-4a99-adc2-97eb3c423317"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[84927, 2001, 736, 18, 57, "Input",ExpressionUUID->"bec4c929-c1a6-4fd0-be18-caf3643ba3a4"],
-Cell[85666, 2021, 14613, 314, 179, "Output",ExpressionUUID->"b48c9c3b-23cc-4864-8ea2-6e98855eb77d"]
+Cell[107587, 2543, 736, 18, 57, "Input",ExpressionUUID->"bec4c929-c1a6-4fd0-be18-caf3643ba3a4"],
+Cell[108326, 2563, 14613, 314, 179, "Output",ExpressionUUID->"b48c9c3b-23cc-4864-8ea2-6e98855eb77d"]
 }, Open  ]]
 }
 ]