MC_update.nb

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

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

(* CreatedBy='Mathematica 12.3' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     57006,       1382]
NotebookOptionsPosition[     55694,       1358]
NotebookOutlinePosition[     56088,       1374]
CellTagsIndexPosition[     56045,       1371]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"ic", "[", "n_", "]"}], ":=", 
  RowBox[{"(", " ", 
   RowBox[{"Array", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"#2", "-", 
          RowBox[{"1", "/", "2"}]}], ")"}], ",", 
        RowBox[{"(", 
         RowBox[{"#1", "-", 
          RowBox[{"1", "/", "2"}]}], ")"}]}], "}"}], "&"}], ",", 
     RowBox[{"{", 
      RowBox[{"n", ",", "n"}], "}"}]}], "]"}], ")"}]}]], "Input",
 CellChangeTimes->{{3.855353027382476*^9, 3.8553530431988707`*^9}, {
   3.855353116475442*^9, 3.855353164172847*^9}, {3.855353211613056*^9, 
   3.855353236540526*^9}, {3.855353331348312*^9, 3.8553533364209948`*^9}, {
   3.855353418613307*^9, 3.855353420697453*^9}, 3.8553554074092493`*^9, {
   3.855359382478222*^9, 3.8553593848285*^9}, {3.855360667085102*^9, 
   3.855360695386457*^9}, {3.85536351483589*^9, 3.855363519501548*^9}, {
   3.8554091659742823`*^9, 3.8554091679333887`*^9}, {3.8554098307471943`*^9, 
   3.855409833245586*^9}, {3.855417350456767*^9, 3.8554173634312677`*^9}, {
   3.855417525183695*^9, 3.8554175257223454`*^9}},
 CellLabel->"In[4]:=",ExpressionUUID->"6e53ff02-1f23-4cb5-a63b-5ef983525da5"],

Cell[BoxData["\[IndentingNewLine]"], "Input",
 CellChangeTimes->{{3.8553553614059477`*^9, 3.855355383573225*^9}, 
   3.8553561955986643`*^9, {3.855356531203327*^9, 3.855356556208063*^9}, {
   3.855356609718405*^9, 3.855356616884458*^9}, {3.855356657204712*^9, 
   3.8553566593587303`*^9}, {3.855359222323327*^9, 3.855359229352557*^9}, {
   3.855359345639184*^9, 3.8553593499972754`*^9}, 3.855359461323807*^9, {
   3.855359672843483*^9, 3.8553596729530487`*^9}, {3.855359725647287*^9, 
   3.855359773408867*^9}, {3.855359805872115*^9, 3.855359825013032*^9}, {
   3.855360199153718*^9, 3.855360199346581*^9}, {3.855360335057404*^9, 
   3.855360346707334*^9}, {3.8553609083522243`*^9, 3.855361011569314*^9}, {
   3.855415621288535*^9, 
   3.85541562230807*^9}},ExpressionUUID->"3d9feee8-3117-4c04-aca1-\
98aed147b460"],

Cell[BoxData[
 RowBox[{"Clear", "[", "\"\<Global`*\>\"", "]"}]], "Input",
 CellChangeTimes->{{3.8553620414877*^9, 3.855362041488234*^9}, 
   3.855362071863247*^9, {3.855362190923656*^9, 3.855362274420238*^9}, {
   3.855362337465797*^9, 3.855362384667659*^9}, {3.855363599915917*^9, 
   3.8553636036864147`*^9}, {3.855363743358302*^9, 3.8553637950277367`*^9}, {
   3.855409997939392*^9, 3.855410007136908*^9}, {3.8554100636977463`*^9, 
   3.855410067695533*^9}, 3.855410617300241*^9},
 CellLabel->
  "(Debug) In[201]:=",ExpressionUUID->"f14827b7-9577-4d0d-992a-2c0c30588c04"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"DIM", "=", "2"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"boxInvTrans", "[", 
    RowBox[{"boxmat_", ",", "pt_"}], "]"}], ":=", 
   RowBox[{"Module", "[", 
    RowBox[{
     RowBox[{"{", "invmat", "}"}], ",", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"invmat", "=", 
       RowBox[{"Inverse", "[", "boxmat", "]"}]}], ";", "\[IndentingNewLine]", 
      
      RowBox[{"invmat", ".", "pt"}]}]}], "\[IndentingNewLine]", "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"boxTrans", "[", 
    RowBox[{"boxmat_", ",", "pt_"}], "]"}], ":=", 
   RowBox[{"boxmat", ".", "pt"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"putInBox", "[", 
    RowBox[{"boxmat_", ",", "pt_"}], "]"}], ":=", 
   RowBox[{"Module", "[", 
    RowBox[{
     RowBox[{"{", "temp", "}"}], ",", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"temp", "=", 
       RowBox[{"boxInvTrans", "[", 
        RowBox[{"boxmat", ",", "pt"}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"Do", "[", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{
         RowBox[{"If", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"Between", "[", 
             RowBox[{
              RowBox[{"temp", "[", 
               RowBox[{"[", "dd", "]"}], "]"}], ",", 
              RowBox[{"{", 
               RowBox[{"0", ",", "1"}], "}"}]}], "]"}], "==", "False"}], ",", 
           
           RowBox[{
            RowBox[{"temp", "[", 
             RowBox[{"[", "dd", "]"}], "]"}], "-=", 
            RowBox[{"Sign", "[", 
             RowBox[{"temp", "[", 
              RowBox[{"[", "dd", "]"}], "]"}], "]"}]}]}], "]"}], ";"}], 
        "\[IndentingNewLine]", ",", 
        RowBox[{"{", 
         RowBox[{"dd", ",", "1", ",", "DIM"}], "}"}]}], "]"}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"boxTrans", "[", 
       RowBox[{"boxmat", ",", "temp"}], "]"}]}]}], "\[IndentingNewLine]", 
    "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"minDist", "[", 
    RowBox[{"boxmat_", ",", "ptA_", ",", "ptB_"}], "]"}], ":=", 
   RowBox[{"Module", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"vA", ",", "vB", ",", "disp"}], "}"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"vA", "=", 
       RowBox[{"boxInvTrans", "[", 
        RowBox[{"boxmat", ",", "ptA"}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"vB", "=", 
       RowBox[{"boxInvTrans", "[", 
        RowBox[{"boxmat", ",", "ptB"}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"disp", "=", 
       RowBox[{"vA", "-", "vB"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"Do", "[", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{
         RowBox[{"If", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"Abs", "[", 
             RowBox[{"disp", "[", 
              RowBox[{"[", "dd", "]"}], "]"}], "]"}], ">", "0.5"}], ",", 
           RowBox[{
            RowBox[{"disp", "[", 
             RowBox[{"[", "dd", "]"}], "]"}], "-=", 
            RowBox[{"Sign", "[", 
             RowBox[{"disp", "[", 
              RowBox[{"[", "dd", "]"}], "]"}], "]"}]}]}], "]"}], ";"}], 
        "\[IndentingNewLine]", ",", 
        RowBox[{"{", 
         RowBox[{"dd", ",", "1", ",", "DIM"}], "}"}]}], "]"}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"boxTrans", "[", 
       RowBox[{"boxmat", ",", "disp"}], "]"}]}]}], "\[IndentingNewLine]", 
    "]"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.855358388530954*^9, 3.855358391767325*^9}, {
  3.8554425287286987`*^9, 3.855442530504261*^9}, {3.855443627316308*^9, 
  3.855443731411318*^9}, {3.855482186172924*^9, 3.8554821942363462`*^9}, {
  3.8567373641625643`*^9, 3.856737398935611*^9}},
 CellLabel->"In[5]:=",ExpressionUUID->"29a3d917-1bbf-4528-a5b1-7b5ccc64b2b9"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"shear", "[", 
     RowBox[{"gamma_", ",", "r_", ",", "n_", ",", "boxmat_"}], "]"}], ":=", 
    RowBox[{"Module", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"shearMax", ",", "rTemp"}], "}"}], ",", "\[IndentingNewLine]", 
      
      RowBox[{
       RowBox[{"shearMax", "=", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"1", ",", "gamma"}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"0", ",", "1"}], "}"}]}], "}"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"rTemp", "=", 
        RowBox[{"Table", "[", 
         RowBox[{
          RowBox[{"putInBox", "[", 
           RowBox[{
            RowBox[{"shearMax", ".", "boxmat"}], ",", 
            RowBox[{"shearMax", ".", 
             RowBox[{"r", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "j"}], "]"}], "]"}]}]}], "]"}], ",", 
          RowBox[{"{", 
           RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"{", 
        RowBox[{"rTemp", ",", 
         RowBox[{"shearMax", ".", "boxmat"}]}], "}"}]}]}], 
     "\[IndentingNewLine]", "]"}]}], ";"}], "\[IndentingNewLine]", 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"eChange", "[", 
     RowBox[{"i_", ",", "j_", ",", "n_", ",", "disList_"}], "]"}], ":=", 
    RowBox[{"(", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{
       RowBox[{"1", "/", "2"}], "*", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Sqrt", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"disList", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "j"}], "]"}], "]"}], "[", 
             RowBox[{"[", "1", "]"}], "]"}], ".", 
            RowBox[{
             RowBox[{"disList", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "j"}], "]"}], "]"}], "[", 
             RowBox[{"[", "1", "]"}], "]"}]}], "]"}], "-", "1"}], ")"}], "^", 
        "2"}]}], "+", 
      RowBox[{
       RowBox[{"1", "/", "2"}], "*", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Sqrt", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"disList", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "j"}], "]"}], "]"}], "[", 
             RowBox[{"[", "2", "]"}], "]"}], ".", 
            RowBox[{
             RowBox[{"disList", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "j"}], "]"}], "]"}], "[", 
             RowBox[{"[", "2", "]"}], "]"}]}], "]"}], "-", "1"}], ")"}], "^", 
        "2"}]}], "+", 
      RowBox[{
       RowBox[{"1", "/", "2"}], "*", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Sqrt", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"disList", "[", 
              RowBox[{"[", 
               RowBox[{
                RowBox[{"Mod", "[", 
                 RowBox[{
                  RowBox[{"i", "-", "1"}], ",", "n", ",", "1"}], "]"}], ",", 
                "j"}], "]"}], "]"}], "[", 
             RowBox[{"[", "2", "]"}], "]"}], ".", 
            RowBox[{
             RowBox[{"disList", "[", 
              RowBox[{"[", 
               RowBox[{
                RowBox[{"Mod", "[", 
                 RowBox[{
                  RowBox[{"i", "-", "1"}], ",", "n", ",", "1"}], "]"}], ",", 
                "j"}], "]"}], "]"}], "[", 
             RowBox[{"[", "2", "]"}], "]"}]}], "]"}], "-", "1"}], ")"}], "^", 
        "2"}]}], "+", 
      RowBox[{
       RowBox[{"1", "/", "2"}], "*", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Sqrt", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"disList", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", 
                RowBox[{"Mod", "[", 
                 RowBox[{
                  RowBox[{"j", "-", "1"}], ",", "n", ",", "1"}], "]"}]}], 
               "]"}], "]"}], "[", 
             RowBox[{"[", "1", "]"}], "]"}], ".", 
            RowBox[{
             RowBox[{"disList", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", 
                RowBox[{"Mod", "[", 
                 RowBox[{
                  RowBox[{"j", "-", "1"}], ",", "n", ",", "1"}], "]"}]}], 
               "]"}], "]"}], "[", 
             RowBox[{"[", "1", "]"}], "]"}]}], "]"}], "-", "1"}], ")"}], "^", 
        "2"}]}]}], ")"}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"disTrial", "[", 
    RowBox[{"i_", ",", "j_", ",", "r_", ",", "n_", ",", "boxMax_"}], "]"}], ":=", 
   RowBox[{"Module", "[", 
    RowBox[{
     RowBox[{"{", "}"}], ",", "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"minDist", "[", 
        RowBox[{"boxMax", ",", 
         RowBox[{"r", "[", 
          RowBox[{"[", 
           RowBox[{"i", ",", "j"}], "]"}], "]"}], ",", 
         RowBox[{"r", "[", 
          RowBox[{"[", 
           RowBox[{"i", ",", 
            RowBox[{"Mod", "[", 
             RowBox[{
              RowBox[{"j", "+", "1"}], ",", "n", ",", "1"}], "]"}]}], "]"}], 
          "]"}]}], "]"}], ",", "\[IndentingNewLine]", 
       RowBox[{"minDist", "[", 
        RowBox[{"boxMax", ",", 
         RowBox[{"r", "[", 
          RowBox[{"[", 
           RowBox[{"i", ",", "j"}], "]"}], "]"}], ",", 
         RowBox[{"r", "[", 
          RowBox[{"[", 
           RowBox[{
            RowBox[{"Mod", "[", 
             RowBox[{
              RowBox[{"i", "+", "1"}], ",", "n", ",", "1"}], "]"}], ",", 
            "j"}], "]"}], "]"}]}], "]"}], ",", "\[IndentingNewLine]", 
       RowBox[{"minDist", "[", 
        RowBox[{"boxMax", ",", 
         RowBox[{"r", "[", 
          RowBox[{"[", 
           RowBox[{"i", ",", 
            RowBox[{"Mod", "[", 
             RowBox[{
              RowBox[{"j", "-", "1"}], ",", "n", ",", "1"}], "]"}]}], "]"}], 
          "]"}], ",", 
         RowBox[{"r", "[", 
          RowBox[{"[", 
           RowBox[{"i", ",", "j"}], "]"}], "]"}]}], "]"}], ",", 
       "\[IndentingNewLine]", 
       RowBox[{"minDist", "[", 
        RowBox[{"boxMax", ",", 
         RowBox[{"r", "[", 
          RowBox[{"[", 
           RowBox[{
            RowBox[{"Mod", "[", 
             RowBox[{
              RowBox[{"i", "-", "1"}], ",", "n", ",", "1"}], "]"}], ",", 
            "j"}], "]"}], "]"}], ",", 
         RowBox[{"r", "[", 
          RowBox[{"[", 
           RowBox[{"i", ",", "j"}], "]"}], "]"}]}], "]"}]}], "}"}]}], 
    "\[IndentingNewLine]", "]"}]}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"eEach", "[", 
     RowBox[{"i_", ",", "j_", ",", "disList_"}], "]"}], ":=", 
    RowBox[{"(", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{
       RowBox[{"1", "/", "2"}], "*", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Sqrt", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"disList", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "j"}], "]"}], "]"}], "[", 
             RowBox[{"[", "1", "]"}], "]"}], ".", 
            RowBox[{
             RowBox[{"disList", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "j"}], "]"}], "]"}], "[", 
             RowBox[{"[", "1", "]"}], "]"}]}], "]"}], "-", "1"}], ")"}], "^", 
        "2"}]}], "+", 
      RowBox[{
       RowBox[{"1", "/", "2"}], "*", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Sqrt", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"disList", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "j"}], "]"}], "]"}], "[", 
             RowBox[{"[", "2", "]"}], "]"}], ".", 
            RowBox[{
             RowBox[{"disList", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "j"}], "]"}], "]"}], "[", 
             RowBox[{"[", "2", "]"}], "]"}]}], "]"}], "-", "1"}], ")"}], "^", 
        "2"}]}]}], ")"}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"eTotal", "[", 
    RowBox[{"disList_", ",", "n_"}], "]"}], ":=", 
   RowBox[{"Module", "[", 
    RowBox[{
     RowBox[{"{", "particleEne", "}"}], ",", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"particleEne", "=", 
       RowBox[{"Table", "[", 
        RowBox[{
         RowBox[{"eEach", "[", 
          RowBox[{"i", ",", "j", ",", "disList"}], "]"}], ",", 
         RowBox[{"{", 
          RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"Total", "[", 
       RowBox[{"particleEne", ",", "2"}], "]"}]}]}], "\[IndentingNewLine]", 
    "]"}]}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"listMinDis", "[", 
    RowBox[{"r_", ",", "n_", ",", "boxMax_"}], "]"}], ":=", 
   RowBox[{"Module", "[", 
    RowBox[{
     RowBox[{"{", "}"}], ",", "\[IndentingNewLine]", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{"minDist", "[", 
          RowBox[{"boxMax", ",", 
           RowBox[{"r", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "j"}], "]"}], "]"}], ",", 
           RowBox[{"r", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", 
              RowBox[{"Mod", "[", 
               RowBox[{
                RowBox[{"j", "+", "1"}], ",", "n", ",", "1"}], "]"}]}], "]"}],
             "]"}]}], "]"}], ",", 
         RowBox[{"minDist", "[", 
          RowBox[{"boxMax", ",", 
           RowBox[{"r", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "j"}], "]"}], "]"}], ",", 
           RowBox[{"r", "[", 
            RowBox[{"[", 
             RowBox[{
              RowBox[{"Mod", "[", 
               RowBox[{
                RowBox[{"i", "+", "1"}], ",", "n", ",", "1"}], "]"}], ",", 
              "j"}], "]"}], "]"}]}], "]"}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], "]"}]}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"listAveEnergy", "[", 
    RowBox[{
    "beta_", ",", "delta_", ",", "r_", ",", "n_", ",", "dgamma_", ",", 
     "numMC_", ",", "numEnsemble_", ",", "boxMax_"}], "]"}], ":=", 
   RowBox[{"Module", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
      "rShear", ",", "disListShear", ",", "sMax", ",", "boxMaxShear"}], "}"}],
      ",", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{
       RowBox[{"{", 
        RowBox[{"rShear", ",", "boxMaxShear"}], "}"}], "=", 
       RowBox[{"shear", "[", 
        RowBox[{"dgamma", ",", "r", ",", "n", ",", "boxMax"}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"disListShear", "=", 
       RowBox[{"listMinDis", "[", 
        RowBox[{"rShear", ",", "n", ",", "boxMaxShear"}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"Do", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"{", 
          RowBox[{"rShear", ",", "disListShear"}], "}"}], "=", 
         RowBox[{"mcstep", "[", 
          RowBox[{
          "beta", ",", "delta", ",", "rShear", ",", "disListShear", ",", "n", 
           ",", "boxMaxShear"}], "]"}]}], ",", "numMC"}], "]"}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"Mean", "[", 
       RowBox[{"Table", "[", 
        RowBox[{
         RowBox[{
          RowBox[{
           RowBox[{"{", 
            RowBox[{"rShear", ",", "disListShear"}], "}"}], "=", 
           RowBox[{"mcstep", "[", 
            RowBox[{
            "beta", ",", "delta", ",", "rShear", ",", "disListShear", ",", 
             "n", ",", "boxMaxShear"}], "]"}]}], ";", 
          RowBox[{"eTotal", "[", 
           RowBox[{"disListShear", ",", "n"}], "]"}]}], ",", "numEnsemble"}], 
        "]"}], "]"}]}]}], "\[IndentingNewLine]", "]"}]}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"listAveEnergy", "[", 
   RowBox[{
   "beta_", ",", "delta_", ",", "r_", ",", "n_", ",", "dgamma_", ",", 
    "numMC_", ",", "numEnsemble_", ",", "boxMax_"}], "]"}], ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"rShear", ",", "disListShear", ",", "sMax", ",", "boxMaxShear"}],
      "}"}], ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{"{", 
       RowBox[{"rShear", ",", "boxMaxShear"}], "}"}], "=", 
      RowBox[{"[", 
       RowBox[{"dgamma", ",", "r", ",", "n", ",", "boxMax"}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"disListShear", "=", 
      RowBox[{"listMinDis", "[", 
       RowBox[{"rShear", ",", "n", ",", "boxMaxShear"}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"Do", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"{", 
         RowBox[{"rShear", ",", "disListShear"}], "}"}], "=", 
        RowBox[{"mcstep", "[", 
         RowBox[{
         "beta", ",", "delta", ",", "rShear", ",", "disListShear", ",", "n", 
          ",", "boxMaxShear"}], "]"}]}], ",", "numMC"}], "]"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"eTotal", "[", 
      RowBox[{"disListShear", ",", "n"}], "]"}]}]}], "\[IndentingNewLine]", 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"listAveEnergySpecial", "[", 
    RowBox[{
    "beta_", ",", "delta_", ",", "r_", ",", "n_", ",", "dgamma_", ",", 
     "numMC_", ",", "numEnsemble_", ",", "boxMax_"}], "]"}], ":=", 
   RowBox[{"Module", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
      "rShear", ",", "disListShear", ",", "sMax", ",", "boxMaxShear"}], "}"}],
      ",", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{
       RowBox[{"{", 
        RowBox[{"rShear", ",", "boxMaxShear"}], "}"}], "=", 
       RowBox[{"shear", "[", 
        RowBox[{"dgamma", ",", "r", ",", "n", ",", "boxMax"}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"disListShear", "=", 
       RowBox[{"listMinDis", "[", 
        RowBox[{"rShear", ",", "n", ",", "boxMaxShear"}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"eTotal", "[", 
       RowBox[{"disListShear", ",", "n"}], "]"}]}]}], "\[IndentingNewLine]", 
    "]"}]}], "\[IndentingNewLine]"}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.855358394641817*^9, 3.855358412688663*^9}, {
   3.855359082572556*^9, 3.85535909595269*^9}, {3.855359136668852*^9, 
   3.855359206445155*^9}, {3.8553607603227043`*^9, 3.8553607619502783`*^9}, {
   3.8553630890193853`*^9, 3.85536310706649*^9}, {3.855363144453569*^9, 
   3.855363229887609*^9}, {3.855363401725505*^9, 3.855363493515493*^9}, {
   3.8553635308287563`*^9, 3.8553635760365458`*^9}, {3.855363638557206*^9, 
   3.855363645217325*^9}, {3.855409218467823*^9, 3.855409244889195*^9}, {
   3.855409385768423*^9, 3.855409411309531*^9}, {3.855409508890375*^9, 
   3.8554095447583847`*^9}, 3.855410170180971*^9, {3.855410224311266*^9, 
   3.855410354123152*^9}, 3.855410405688656*^9, {3.8554104588924503`*^9, 
   3.855410498826261*^9}, {3.855415981439232*^9, 3.8554160292910547`*^9}, {
   3.8554174370523*^9, 3.855417446265965*^9}, {3.855417567935526*^9, 
   3.8554175949577312`*^9}, {3.855418522716997*^9, 3.855418524294128*^9}, {
   3.8554192319763308`*^9, 3.855419241418633*^9}, {3.85541945236158*^9, 
   3.8554196371560926`*^9}, {3.855419838096898*^9, 3.8554198484765882`*^9}, {
   3.8554222061770363`*^9, 3.855422270152869*^9}, {3.855423684378697*^9, 
   3.855423698018071*^9}, {3.855423824900527*^9, 3.855423860330923*^9}, {
   3.855424034824624*^9, 3.8554240743545017`*^9}, {3.8554241221318893`*^9, 
   3.855424166839778*^9}, {3.855439322993415*^9, 3.855439414260227*^9}, {
   3.855439742622282*^9, 3.855439764907737*^9}, 3.85544015090539*^9, {
   3.855444896826777*^9, 3.855444904738309*^9}, {3.856707143754594*^9, 
   3.856707201927875*^9}, {3.8567186476587057`*^9, 3.856718714060413*^9}, {
   3.8567187643963203`*^9, 3.8567187756677628`*^9}, {3.8567188062319183`*^9, 
   3.85671883889991*^9}, {3.856718956076198*^9, 3.8567189796166887`*^9}, {
   3.8567190560260973`*^9, 3.856719088634119*^9}, {3.856719259901279*^9, 
   3.856719428583879*^9}, {3.856719464197022*^9, 3.85671957680366*^9}, {
   3.856721008508321*^9, 3.856721017836043*^9}, {3.856721047969849*^9, 
   3.856721251071269*^9}, {3.8567213132588053`*^9, 3.8567213291132927`*^9}, {
   3.856736337676385*^9, 3.8567364030514193`*^9}, {3.856737695524588*^9, 
   3.856737785769346*^9}, 3.856737901919362*^9, {3.856738286445033*^9, 
   3.856738290871867*^9}, {3.856738330894532*^9, 3.856738389541449*^9}, {
   3.8567384313191977`*^9, 3.8567385253008223`*^9}, 3.856738728543413*^9, {
   3.8567400574509983`*^9, 3.8567400644722643`*^9}, {3.85674013342439*^9, 
   3.856740141618972*^9}, {3.8567415081706047`*^9, 3.856741508299527*^9}, {
   3.856802278851857*^9, 3.856802297111394*^9}, {3.8568027485941143`*^9, 
   3.856802760450788*^9}, {3.856802856508369*^9, 3.85680287950303*^9}, {
   3.856802923172058*^9, 3.856802934131516*^9}, {3.856802989820303*^9, 
   3.8568029983509417`*^9}, {3.856890550305275*^9, 3.856890563602857*^9}, {
   3.856906366838524*^9, 3.856906367718025*^9}, {3.8569064802922897`*^9, 
   3.856906497848835*^9}, {3.856907423689673*^9, 3.8569074244296637`*^9}, {
   3.8577322306936283`*^9, 3.857732236910336*^9}, 3.857922757290866*^9, {
   3.858014130138616*^9, 3.85801413237744*^9}, {3.858167358779636*^9, 
   3.8581673620008593`*^9}},ExpressionUUID->"8f7d00d6-4336-491c-9451-\
3e91ee437232"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"mcstep", "[", 
     RowBox[{
     "beta_", ",", "delta_", ",", "r_", ",", "disList_", ",", "n_", ",", 
      "boxMax_"}], "]"}], ":=", 
    RowBox[{"Module", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
       "i", ",", "j", ",", "dx", ",", "dy", ",", "de1", ",", "de2", ",", 
        "deval", ",", "success", ",", "rTemp", ",", "rTrial", ",", 
        "disListTrial", ",", "disListTemp"}], "}"}], ",", 
      "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"rTemp", "=", "r"}], ";", "\[IndentingNewLine]", 
       RowBox[{"rTrial", "=", "r"}], ";", "\[IndentingNewLine]", 
       RowBox[{"disListTemp", "=", "disList"}], ";", "\[IndentingNewLine]", 
       RowBox[{"i", "=", 
        RowBox[{"RandomInteger", "[", 
         RowBox[{"{", 
          RowBox[{"1", ",", "n"}], "}"}], "]"}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"j", "=", 
        RowBox[{"RandomInteger", "[", 
         RowBox[{"{", 
          RowBox[{"1", ",", "n"}], "}"}], "]"}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"dx", " ", "=", " ", 
        RowBox[{"RandomReal", "[", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"-", "delta"}], ",", "delta"}], "}"}], "]"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"dy", " ", "=", " ", 
        RowBox[{"RandomReal", "[", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"-", "delta"}], ",", "delta"}], "}"}], "]"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"de1", " ", "=", " ", 
        RowBox[{"eChange", "[", 
         RowBox[{"i", ",", "j", ",", "n", ",", "disListTemp"}], "]"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"(*", 
        RowBox[{
         RowBox[{"Print", "[", "de1", "]"}], ";"}], "*)"}], 
       "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{
         RowBox[{"rTrial", "[", 
          RowBox[{"[", 
           RowBox[{"i", ",", "j"}], "]"}], "]"}], "[", 
         RowBox[{"[", "1", "]"}], "]"}], "+=", "dx"}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{
         RowBox[{"rTrial", "[", 
          RowBox[{"[", 
           RowBox[{"i", ",", "j"}], "]"}], "]"}], "[", 
         RowBox[{"[", "2", "]"}], "]"}], "+=", "dy"}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"(*", 
        RowBox[{
         RowBox[{"Print", "[", 
          RowBox[{"rTrial", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "j"}], "]"}], "]"}], "]"}], ";"}], "*)"}], 
       "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"rTrial", "[", 
         RowBox[{"[", 
          RowBox[{"i", ",", "j"}], "]"}], "]"}], "=", 
        RowBox[{"putInBox", "[", 
         RowBox[{"boxMax", ",", 
          RowBox[{"rTrial", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "j"}], "]"}], "]"}]}], "]"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"disListTrial", "=", 
        RowBox[{"disTrial", "[", 
         RowBox[{"i", ",", "j", ",", "rTrial", ",", "n", ",", "boxMax"}], 
         "]"}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"(*", 
        RowBox[{
         RowBox[{"Print", "[", "disListTrial", "]"}], ";"}], "*)"}], 
       "\[IndentingNewLine]", 
       RowBox[{"de2", " ", "=", 
        RowBox[{
         RowBox[{
          RowBox[{"1", "/", "2"}], "*", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"Sqrt", "[", 
              RowBox[{
               RowBox[{"disListTrial", "[", 
                RowBox[{"[", "1", "]"}], "]"}], ".", 
               RowBox[{"disListTrial", "[", 
                RowBox[{"[", "1", "]"}], "]"}]}], "]"}], "-", "1"}], ")"}], 
           "^", "2"}]}], "+", 
         RowBox[{
          RowBox[{"1", "/", "2"}], "*", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"Sqrt", "[", 
              RowBox[{
               RowBox[{"disListTrial", "[", 
                RowBox[{"[", "2", "]"}], "]"}], ".", 
               RowBox[{"disListTrial", "[", 
                RowBox[{"[", "2", "]"}], "]"}]}], "]"}], "-", "1"}], ")"}], 
           "^", "2"}]}], "+", 
         RowBox[{
          RowBox[{"1", "/", "2"}], "*", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"Sqrt", "[", 
              RowBox[{
               RowBox[{"disListTrial", "[", 
                RowBox[{"[", "3", "]"}], "]"}], ".", 
               RowBox[{"disListTrial", "[", 
                RowBox[{"[", "3", "]"}], "]"}]}], "]"}], "-", "1"}], ")"}], 
           "^", "2"}]}], "+", 
         RowBox[{
          RowBox[{"1", "/", "2"}], "*", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"Sqrt", "[", 
              RowBox[{
               RowBox[{"disListTrial", "[", 
                RowBox[{"[", "4", "]"}], "]"}], ".", 
               RowBox[{"disListTrial", "[", 
                RowBox[{"[", "4", "]"}], "]"}]}], "]"}], "-", "1"}], ")"}], 
           "^", "2"}]}]}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"deval", " ", "=", " ", 
        RowBox[{"de2", "-", "de1"}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"success", "=", " ", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{"deval", "<=", "0"}], ",", "True", ",", 
          RowBox[{
           RowBox[{"RandomReal", "[", "]"}], "<", " ", 
           RowBox[{"Exp", "[", 
            RowBox[{
             RowBox[{"-", "beta"}], "*", "deval"}], "]"}]}]}], "]"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"rTemp", "[", 
         RowBox[{"[", 
          RowBox[{"i", ",", "j"}], "]"}], "]"}], "=", 
        RowBox[{"If", "[", 
         RowBox[{"success", ",", 
          RowBox[{"rTrial", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "j"}], "]"}], "]"}], ",", 
          RowBox[{"rTemp", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "j"}], "]"}], "]"}]}], "]"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"disListTemp", "[", 
         RowBox[{"[", 
          RowBox[{"i", ",", "j"}], "]"}], "]"}], "=", 
        RowBox[{"If", "[", 
         RowBox[{"success", ",", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"disListTrial", "[", 
             RowBox[{"[", "1", "]"}], "]"}], ",", 
            RowBox[{"disListTrial", "[", 
             RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", 
          RowBox[{"disListTemp", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "j"}], "]"}], "]"}]}], "]"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{
         RowBox[{"disListTemp", "[", 
          RowBox[{"[", 
           RowBox[{"i", ",", 
            RowBox[{"Mod", "[", 
             RowBox[{
              RowBox[{"j", "-", "1"}], ",", "n", ",", "1"}], "]"}]}], "]"}], 
          "]"}], "[", 
         RowBox[{"[", "1", "]"}], "]"}], "=", 
        RowBox[{"If", "[", 
         RowBox[{"success", ",", 
          RowBox[{"disListTrial", "[", 
           RowBox[{"[", "3", "]"}], "]"}], ",", 
          RowBox[{
           RowBox[{"disListTemp", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", 
              RowBox[{"Mod", "[", 
               RowBox[{
                RowBox[{"j", "-", "1"}], ",", "n", ",", "1"}], "]"}]}], "]"}],
             "]"}], "[", 
           RowBox[{"[", "1", "]"}], "]"}]}], "]"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{
         RowBox[{"disListTemp", "[", 
          RowBox[{"[", 
           RowBox[{
            RowBox[{"Mod", "[", 
             RowBox[{
              RowBox[{"i", "-", "1"}], ",", "n", ",", "1"}], "]"}], ",", 
            "j"}], "]"}], "]"}], "[", 
         RowBox[{"[", "2", "]"}], "]"}], "=", 
        RowBox[{"If", "[", 
         RowBox[{"success", ",", 
          RowBox[{"disListTrial", "[", 
           RowBox[{"[", "4", "]"}], "]"}], ",", 
          RowBox[{
           RowBox[{"disListTemp", "[", 
            RowBox[{"[", 
             RowBox[{
              RowBox[{"Mod", "[", 
               RowBox[{
                RowBox[{"i", "-", "1"}], ",", "n", ",", "1"}], "]"}], ",", 
              "j"}], "]"}], "]"}], "[", 
           RowBox[{"[", "2", "]"}], "]"}]}], "]"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"{", 
        RowBox[{"rTemp", ",", "disListTemp"}], "}"}]}]}], 
     "\[IndentingNewLine]", "\[IndentingNewLine]", "]"}]}], ";"}], 
  "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{{3.855416062549124*^9, 3.855416235472672*^9}, 
   3.855416900881032*^9, {3.8554182262426043`*^9, 3.855418265652061*^9}, {
   3.855418347604722*^9, 3.855418351664425*^9}, {3.8554188552425013`*^9, 
   3.8554188569434643`*^9}, {3.8554192617871847`*^9, 3.855419302031432*^9}, {
   3.855422296414673*^9, 3.855422335055786*^9}, {3.855422879981237*^9, 
   3.855422887225459*^9}, {3.855424373954006*^9, 3.855424374116781*^9}, {
   3.855439142356697*^9, 3.855439283706429*^9}, {3.855440119181386*^9, 
   3.8554401398151493`*^9}, {3.8554403058839197`*^9, 3.855440314946044*^9}, {
   3.855440545483156*^9, 3.855440556256096*^9}, {3.855443127946228*^9, 
   3.855443137856559*^9}, {3.855443257744772*^9, 3.855443257869236*^9}, {
   3.855443328066093*^9, 3.855443329517392*^9}, {3.8554437536628847`*^9, 
   3.855443762648179*^9}, {3.8567076474461517`*^9, 3.856707647575624*^9}, {
   3.856707751342915*^9, 3.856707774226264*^9}, {3.856708574410293*^9, 
   3.856708577579596*^9}, {3.8567086495003633`*^9, 3.856708650576357*^9}, {
   3.856719607425082*^9, 3.856719611797989*^9}, {3.856720230916451*^9, 
   3.856720235702379*^9}, {3.85672026807537*^9, 3.856720277869816*^9}, {
   3.856720317607992*^9, 3.8567204041559877`*^9}, {3.856720525123188*^9, 
   3.856720601927525*^9}, {3.856720633704954*^9, 3.856720769684989*^9}, {
   3.8567208138160458`*^9, 3.856720845384816*^9}, {3.856720875615919*^9, 
   3.8567209005628967`*^9}, {3.85673611364301*^9, 3.856736117818503*^9}, {
   3.856736537727722*^9, 3.8567366022510347`*^9}, {3.856803062187553*^9, 
   3.85680307440697*^9}, {3.856804733293663*^9, 3.856804753456214*^9}, {
   3.856805010626753*^9, 3.856805018085993*^9}, 3.8568050484462433`*^9, {
   3.8568052987538967`*^9, 3.856805305243767*^9}, {3.8575624628006763`*^9, 
   3.857562466978649*^9}, {3.85756269789878*^9, 3.857562700796177*^9}, {
   3.857563769795948*^9, 3.857563776672401*^9}, {3.857564437217531*^9, 
   3.857564479343443*^9}, {3.8575744698281918`*^9, 3.85757447388426*^9}, {
   3.8577322684253387`*^9, 3.857732269454955*^9}, 3.8577325220297947`*^9},
 CellLabel->"In[19]:=",ExpressionUUID->"557bb1dc-20b5-499c-9529-a8a50b411111"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"Thermalizing", " ", "the", " ", "initial", " ", "state"}], "*)"}],
   "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"modulus", "[", 
      RowBox[{"beta_", ",", "delta_", ",", "n_", ",", "numEnsemble_"}], "]"}],
      ":=", 
     RowBox[{"Module", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
        "numMC", ",", "rIni", ",", "dgamma", ",", "disList", ",", "boxMax", 
         ",", " ", "enList", ",", "derive1", ",", "derive2", ",", "g"}], 
        "}"}], ",", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"numMC", "=", 
         RowBox[{"100", "*", 
          RowBox[{"n", "^", "2"}]}]}], ";", "\[IndentingNewLine]", 
        RowBox[{"dgamma", "=", 
         RowBox[{"10", "^", 
          RowBox[{"(", 
           RowBox[{"-", "4"}], ")"}]}]}], ";", "\[IndentingNewLine]", 
        RowBox[{"boxMax", "=", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"n", ",", "0"}], "}"}], ",", 
           RowBox[{"{", 
            RowBox[{"0", ",", "n"}], "}"}]}], "}"}]}], ";", 
        "\[IndentingNewLine]", 
        RowBox[{"rIni", "=", 
         RowBox[{"ic", "[", "n", "]"}]}], ";", "\[IndentingNewLine]", 
        RowBox[{"disList", "=", " ", 
         RowBox[{"listMinDis", "[", 
          RowBox[{"rIni", ",", "n", ",", "boxMax"}], "]"}]}], ";", 
        "\[IndentingNewLine]", 
        RowBox[{"enList", "=", 
         RowBox[{"Table", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"Do", "[", 
             RowBox[{
              RowBox[{
               RowBox[{"{", 
                RowBox[{"rIni", ",", "disList"}], "}"}], "=", 
               RowBox[{"mcstep", "[", 
                RowBox[{
                "beta", ",", "delta", ",", "rIni", ",", "disList", ",", "n", 
                 ",", "boxMax"}], "]"}]}], ",", "numMC"}], "]"}], ";", 
            RowBox[{"{", 
             RowBox[{
              RowBox[{"listAveEnergy", "[", 
               RowBox[{"beta", ",", "delta", ",", "rIni", ",", "n", ",", 
                RowBox[{"-", "dgamma"}], ",", "numMC", ",", "numEnsemble", 
                ",", "boxMax"}], "]"}], ",", 
              RowBox[{"eTotal", "[", 
               RowBox[{"disList", ",", "n"}], "]"}], ",", 
              RowBox[{"listAveEnergy", "[", 
               RowBox[{
               "beta", ",", "delta", ",", "rIni", ",", "n", ",", "dgamma", 
                ",", "numMC", ",", "numEnsemble", ",", "boxMax"}], "]"}]}], 
             "}"}]}], ",", "numEnsemble"}], "]"}]}], ";", 
        "\[IndentingNewLine]", 
        RowBox[{"derive1", "=", 
         RowBox[{"Table", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"(", 
             RowBox[{
              RowBox[{
               RowBox[{"enList", "[", 
                RowBox[{"[", "k", "]"}], "]"}], "[", 
               RowBox[{"[", "3", "]"}], "]"}], "-", 
              RowBox[{
               RowBox[{"enList", "[", 
                RowBox[{"[", "k", "]"}], "]"}], "[", 
               RowBox[{"[", "1", "]"}], "]"}]}], ")"}], "/", 
            RowBox[{"(", 
             RowBox[{"2", "*", "dgamma"}], ")"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"k", ",", "1", ",", "numEnsemble"}], "}"}]}], "]"}]}], 
        ";", "\[IndentingNewLine]", 
        RowBox[{"derive2", "=", 
         RowBox[{"Table", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"(", 
             RowBox[{
              RowBox[{
               RowBox[{"enList", "[", 
                RowBox[{"[", "k", "]"}], "]"}], "[", 
               RowBox[{"[", "3", "]"}], "]"}], "+", 
              RowBox[{
               RowBox[{"enList", "[", 
                RowBox[{"[", "k", "]"}], "]"}], "[", 
               RowBox[{"[", "1", "]"}], "]"}], "-", 
              RowBox[{"2", "*", 
               RowBox[{
                RowBox[{"enList", "[", 
                 RowBox[{"[", "k", "]"}], "]"}], "[", 
                RowBox[{"[", "2", "]"}], "]"}]}]}], ")"}], "/", 
            RowBox[{"(", 
             RowBox[{"dgamma", "^", "2"}], ")"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"k", ",", "1", ",", "numEnsemble"}], "}"}]}], "]"}]}], 
        ";", "\[IndentingNewLine]", 
        RowBox[{"(*", 
         RowBox[{
          RowBox[{"Print", "[", "\"\<initialized\>\"", "]"}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{"Print", "[", "enList", "]"}], ";"}], "*)"}], 
        "\[IndentingNewLine]", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"g", "=", 
           RowBox[{
            RowBox[{"1", "/", 
             RowBox[{"(", 
              RowBox[{"n", "*", "n"}], ")"}]}], "*", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"Mean", "[", "derive2", "]"}], "-", 
              RowBox[{"beta", "*", 
               RowBox[{"Variance", "[", "derive1", "]"}]}]}], ")"}]}]}], ",", 
          
          RowBox[{"Mean", "[", "derive2", "]"}], ",", 
          RowBox[{"beta", "*", 
           RowBox[{"Variance", "[", "derive1", "]"}]}]}], "}"}]}]}], 
      "\[IndentingNewLine]", "]"}]}], ";"}], 
   "\[IndentingNewLine]"}]}]], "Input",
 CellChangeTimes->{{3.8554169275027933`*^9, 3.855417115599173*^9}, {
   3.855417156154703*^9, 3.855417221671495*^9}, {3.8554176186355762`*^9, 
   3.855417665899128*^9}, 3.855417704955613*^9, 3.855418288021029*^9, {
   3.85541836854547*^9, 3.855418397612687*^9}, {3.855418443329171*^9, 
   3.855418465477211*^9}, {3.8554185112247753`*^9, 3.855418585256837*^9}, {
   3.855418628665122*^9, 3.855418642897135*^9}, {3.855418677964911*^9, 
   3.85541877112856*^9}, {3.8554188741242943`*^9, 3.8554188893755083`*^9}, {
   3.855419076164253*^9, 3.855419181349498*^9}, {3.855419313850689*^9, 
   3.8554193489649487`*^9}, {3.855419398760476*^9, 3.855419409520954*^9}, {
   3.855419887782578*^9, 3.855419893848568*^9}, {3.855419933184188*^9, 
   3.855419989088318*^9}, {3.855420110895649*^9, 3.855420163976356*^9}, 
   3.85542024902664*^9, {3.8554202795170507`*^9, 3.855420333253284*^9}, {
   3.8554203735088987`*^9, 3.8554204957212687`*^9}, {3.8554206998299723`*^9, 
   3.85542070033788*^9}, {3.855420747925807*^9, 3.855420766673534*^9}, {
   3.855420828819137*^9, 3.855420909482993*^9}, {3.855421032425424*^9, 
   3.855421212555324*^9}, {3.855422367237555*^9, 3.855422409042226*^9}, {
   3.855422494752393*^9, 3.855422508143042*^9}, {3.85543650873053*^9, 
   3.855436512901725*^9}, {3.855436567882957*^9, 3.855436663870748*^9}, {
   3.855436710544033*^9, 3.8554367455326157`*^9}, {3.85543800207829*^9, 
   3.855438115945752*^9}, {3.855438159495386*^9, 3.855438165319077*^9}, {
   3.85543827877501*^9, 3.855438430249423*^9}, {3.855444079631896*^9, 
   3.855444080241045*^9}, {3.855444167588789*^9, 3.855444170665806*^9}, {
   3.855444205287142*^9, 3.855444216140774*^9}, {3.855444572789743*^9, 
   3.8554445890254908`*^9}, {3.8554448024184647`*^9, 3.855444813277007*^9}, {
   3.8554449837645197`*^9, 3.855445017734879*^9}, 3.8554452440396433`*^9, {
   3.855445387377528*^9, 3.855445497210376*^9}, {3.855445601438259*^9, 
   3.855445613498904*^9}, {3.855445730410639*^9, 3.855445808338111*^9}, {
   3.8554458583397303`*^9, 3.8554459033344097`*^9}, {3.855445936401595*^9, 
   3.855445998970924*^9}, {3.855446084959408*^9, 3.8554461111830807`*^9}, {
   3.8554461799540462`*^9, 3.855446187198318*^9}, {3.855678012941785*^9, 
   3.8556780154041758`*^9}, {3.85568084354602*^9, 3.855680847989657*^9}, {
   3.855680899798092*^9, 3.855680902357273*^9}, {3.855681499853311*^9, 
   3.855681535777308*^9}, 3.855740774603015*^9, {3.856716474651723*^9, 
   3.8567166278432198`*^9}, {3.856718156953952*^9, 3.856718160924007*^9}, {
   3.8567338650975847`*^9, 3.856733869326062*^9}, {3.856735284217683*^9, 
   3.856735314949236*^9}, {3.856739387797792*^9, 3.856739490089493*^9}, {
   3.856739715020706*^9, 3.85673994235324*^9}, 3.856740036377955*^9, {
   3.8567401083531713`*^9, 3.856740111011306*^9}, {3.85674016675718*^9, 
   3.856740198565027*^9}, {3.856740283352112*^9, 3.8567403126229267`*^9}, {
   3.856740362783722*^9, 3.856740538876555*^9}, {3.856740860957964*^9, 
   3.856740865987665*^9}, {3.856740984489504*^9, 3.85674101163253*^9}, 
   3.856741738241521*^9, {3.856890618908861*^9, 3.8568906189817*^9}, {
   3.857919158649033*^9, 3.857919164331808*^9}, {3.8579238871697483`*^9, 
   3.857923887997295*^9}, 3.857946153921307*^9, 3.8579478205749474`*^9, 
   3.8579487185787783`*^9, 3.857949407400175*^9, {3.857949458984645*^9, 
   3.857949481734419*^9}, {3.8579969192866497`*^9, 3.857996921061884*^9}, 
   3.858012812513818*^9, {3.8580135250500727`*^9, 3.858013537765419*^9}, 
   3.858105336196967*^9},
 CellLabel->
  "In[130]:=",ExpressionUUID->"cffde510-1720-4a7f-b240-80202c3a3c96"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "Without", " ", "thermalizing", " ", "the", " ", "initial", " ", "state"}],
    "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"modulusSpecial", "[", 
     RowBox[{"beta_", ",", "delta_", ",", "n_", ",", "numEnsemble_"}], "]"}], 
    ":=", 
    RowBox[{"Module", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
       "numMC", ",", "rIni", ",", "dgamma", ",", "disList", ",", "boxMax", 
        ",", " ", "enList", ",", "derive1", ",", "derive2", ",", "g"}], "}"}],
       ",", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"numMC", "=", 
        RowBox[{"100", "*", 
         RowBox[{"n", "^", "2"}]}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"dgamma", "=", 
        RowBox[{"10", "^", 
         RowBox[{"(", 
          RowBox[{"-", "6"}], ")"}]}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"boxMax", "=", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"n", ",", "0"}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"0", ",", "n"}], "}"}]}], "}"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"rIni", "=", 
        RowBox[{"ic", "[", "n", "]"}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"disList", "=", " ", 
        RowBox[{"listMinDis", "[", 
         RowBox[{"rIni", ",", "n", ",", "boxMax"}], "]"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"enList", "=", 
        RowBox[{"Table", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"Do", "[", 
            RowBox[{
             RowBox[{
              RowBox[{"{", 
               RowBox[{"rIni", ",", "disList"}], "}"}], "=", 
              RowBox[{"mcstep", "[", 
               RowBox[{
               "beta", ",", "delta", ",", "rIni", ",", "disList", ",", "n", 
                ",", "boxMax"}], "]"}]}], ",", "numMC"}], "]"}], ";", 
           RowBox[{"{", 
            RowBox[{
             RowBox[{"listAveEnergySpecial", "[", 
              RowBox[{"beta", ",", "delta", ",", "rIni", ",", "n", ",", 
               RowBox[{"-", "dgamma"}], ",", "numMC", ",", "numEnsemble", ",",
                "boxMax"}], "]"}], ",", 
             RowBox[{"eTotal", "[", 
              RowBox[{"disList", ",", "n"}], "]"}], ",", 
             RowBox[{"listAveEnergySpecial", "[", 
              RowBox[{
              "beta", ",", "delta", ",", "rIni", ",", "n", ",", "dgamma", ",",
                "numMC", ",", "numEnsemble", ",", "boxMax"}], "]"}]}], 
            "}"}]}], ",", "numEnsemble"}], "]"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"derive1", "=", 
        RowBox[{"Table", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{
              RowBox[{"enList", "[", 
               RowBox[{"[", "k", "]"}], "]"}], "[", 
              RowBox[{"[", "3", "]"}], "]"}], "-", 
             RowBox[{
              RowBox[{"enList", "[", 
               RowBox[{"[", "k", "]"}], "]"}], "[", 
              RowBox[{"[", "1", "]"}], "]"}]}], ")"}], "/", 
           RowBox[{"(", 
            RowBox[{"2", "*", "dgamma"}], ")"}]}], ",", 
          RowBox[{"{", 
           RowBox[{"k", ",", "1", ",", "numEnsemble"}], "}"}]}], "]"}]}], ";",
        "\[IndentingNewLine]", 
       RowBox[{"derive2", "=", 
        RowBox[{"Table", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{
              RowBox[{"enList", "[", 
               RowBox[{"[", "k", "]"}], "]"}], "[", 
              RowBox[{"[", "3", "]"}], "]"}], "+", 
             RowBox[{
              RowBox[{"enList", "[", 
               RowBox[{"[", "k", "]"}], "]"}], "[", 
              RowBox[{"[", "1", "]"}], "]"}], "-", 
             RowBox[{"2", "*", 
              RowBox[{
               RowBox[{"enList", "[", 
                RowBox[{"[", "k", "]"}], "]"}], "[", 
               RowBox[{"[", "2", "]"}], "]"}]}]}], ")"}], "/", 
           RowBox[{"(", 
            RowBox[{"dgamma", "^", "2"}], ")"}]}], ",", 
          RowBox[{"{", 
           RowBox[{"k", ",", "1", ",", "numEnsemble"}], "}"}]}], "]"}]}], ";",
        "\[IndentingNewLine]", 
       RowBox[{"(*", 
        RowBox[{
         RowBox[{"Print", "[", "\"\<initialized\>\"", "]"}], ";", 
         "\[IndentingNewLine]", 
         RowBox[{"Print", "[", "enList", "]"}], ";"}], "*)"}], 
       "\[IndentingNewLine]", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"g", "=", 
          RowBox[{
           RowBox[{"1", "/", 
            RowBox[{"(", 
             RowBox[{"n", "*", "n"}], ")"}]}], "*", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"Mean", "[", "derive2", "]"}], "-", 
             RowBox[{"beta", "*", 
              RowBox[{"Variance", "[", "derive1", "]"}]}]}], ")"}]}]}], ",", 
         RowBox[{"Mean", "[", "derive2", "]"}], ",", 
         RowBox[{"beta", "*", 
          RowBox[{"Variance", "[", "derive1", "]"}]}]}], "}"}]}]}], 
     "\[IndentingNewLine]", "]"}]}], ";"}]}]], "Input",
 CellChangeTimes->{{3.855421595379044*^9, 3.855421622913954*^9}, {
   3.8554217276064577`*^9, 3.855421728458218*^9}, {3.855446985279882*^9, 
   3.855447105720294*^9}, {3.8554472451440763`*^9, 3.855447253164412*^9}, {
   3.855447324072555*^9, 3.8554473415558567`*^9}, {3.8554476351230373`*^9, 
   3.855447647589738*^9}, 3.855447802058497*^9, {3.8556814497413673`*^9, 
   3.8556814933117933`*^9}, {3.8557396963291597`*^9, 3.855739701428588*^9}, {
   3.856907436331402*^9, 3.8569074642349453`*^9}, 3.857925039580636*^9, 
   3.857946150820897*^9, {3.85794661823838*^9, 3.85794661831098*^9}, 
   3.857946715242516*^9, 3.857946878253088*^9, 3.85799694333389*^9, 
   3.8580010263817453`*^9, 3.8580034270021257`*^9, {3.858114031326931*^9, 
   3.858114040271694*^9}, {3.8581671719978867`*^9, 3.8581671720989*^9}},
 CellLabel->
  "In[228]:=",ExpressionUUID->"da090306-c072-43f6-af23-db60b6ef2112"],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.857908257696047*^9, 
  3.857908261811071*^9}},ExpressionUUID->"334b6491-5bd5-421e-b511-\
41a5b6eecb5f"],

Cell[BoxData[
 RowBox[{
  RowBox[{"findDelta", "[", "b_", "]"}], ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
     "trialD", ",", "trialDPre", ",", "increaseBy", ",", "successf", ",", 
      "successRate", ",", "nf", ",", "rf", ",", "boxf", ",", "dlf", ",", 
      "rfPre", ",", "dlfPre", ",", "rfTest", ",", "dlfTest", ",", "lowD", ",",
       "highD"}], "}"}], ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"increaseBy", "=", "0.01"}], ";", "\[IndentingNewLine]", 
     RowBox[{"trialD", "=", 
      RowBox[{
       RowBox[{"Sqrt", "[", 
        RowBox[{"1", "/", "b"}], "]"}], "/", "5"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"nf", "=", "3"}], ";", "\[IndentingNewLine]", 
     RowBox[{"rf", "=", 
      RowBox[{"ic", "[", "nf", "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"boxf", "=", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"nf", ",", "0"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"0", ",", "nf"}], "}"}]}], "}"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"dlf", "=", 
      RowBox[{"listMinDis", "[", 
       RowBox[{"rf", ",", "nf", ",", "boxf"}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"Do", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"{", 
         RowBox[{"rf", ",", "dlf"}], "}"}], "=", 
        RowBox[{"mcstep", "[", 
         RowBox[{
         "b", ",", "trialD", ",", "rf", ",", "dlf", ",", "nf", ",", "boxf"}], 
         "]"}]}], ",", "3000"}], "]"}], ";", "\[IndentingNewLine]", 
     RowBox[{"successRate", "=", "1"}], ";", "\[IndentingNewLine]", 
     RowBox[{"trialDPre", "=", "trialD"}], ";", "\[IndentingNewLine]", 
     RowBox[{"While", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"successRate", ">", "0.5"}], "||", 
        RowBox[{"successRate", "<", "0.3"}]}], ",", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"increaseBy", "=", 
         RowBox[{"If", "[", 
          RowBox[{
           RowBox[{"successRate", "<", "0.3"}], ",", 
           RowBox[{
            RowBox[{"trialD", "=", "trialDPre"}], ";", 
            RowBox[{"increaseBy", "*", "0.5"}]}], ",", "increaseBy"}], 
          "]"}]}], ";", "\[IndentingNewLine]", 
        RowBox[{"trialD", "=", 
         RowBox[{"trialD", "*", 
          RowBox[{"(", 
           RowBox[{"1", "+", "increaseBy"}], ")"}]}]}], ";", 
        "\[IndentingNewLine]", 
        RowBox[{"trialDPre", "=", "trialD"}], ";", "\[IndentingNewLine]", 
        RowBox[{"rfTest", "=", "rf"}], ";", "\[IndentingNewLine]", 
        RowBox[{"dlfTest", "=", "dlf"}], ";", "\[IndentingNewLine]", 
        RowBox[{"successf", "=", "0"}], ";", "\[IndentingNewLine]", 
        RowBox[{"Do", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"rfPre", "=", "rfTest"}], ";", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"rfTest", ",", "dlfTest"}], "}"}], "=", 
            RowBox[{"mcstep", "[", 
             RowBox[{
             "b", ",", "trialD", ",", "rfTest", ",", "dlfTest", ",", "nf", 
              ",", "boxf"}], "]"}]}], ";", "\[IndentingNewLine]", 
           RowBox[{"successf", "=", 
            RowBox[{"If", "[", 
             RowBox[{
              RowBox[{"rfTest", "!=", "rfPre"}], ",", 
              RowBox[{"successf", "+", "1"}], ",", "successf"}], "]"}]}]}], 
          "\[IndentingNewLine]", ",", "100"}], "]"}], ";", 
        "\[IndentingNewLine]", 
        RowBox[{"successRate", "=", 
         RowBox[{"successf", "/", "100"}]}], ";"}]}], "\[IndentingNewLine]", 
      "]"}], ";", "\[IndentingNewLine]", "trialD"}]}], "\[IndentingNewLine]", 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.857908268171382*^9, 3.857908459010537*^9}, {
   3.8579085463809443`*^9, 3.857908694059441*^9}, {3.8579087335625153`*^9, 
   3.857908744400371*^9}, {3.857909122310005*^9, 3.857909183615252*^9}, {
   3.857909410282154*^9, 3.857909476769559*^9}, {3.857909756570613*^9, 
   3.85790976744886*^9}, {3.857909918453526*^9, 3.8579100016861258`*^9}, {
   3.8579104752488947`*^9, 3.857910476351824*^9}, {3.8579105552635927`*^9, 
   3.857910601613515*^9}, {3.857910780742284*^9, 3.857911103949154*^9}, {
   3.857911191539918*^9, 3.857911263526147*^9}, {3.857912893360674*^9, 
   3.8579129245627813`*^9}, {3.857913076766694*^9, 3.857913236649433*^9}, {
   3.8579134494927692`*^9, 3.857913488337515*^9}, {3.8579135202878532`*^9, 
   3.857913520392115*^9}, {3.857913616380848*^9, 3.857913623255724*^9}, {
   3.8579140498925047`*^9, 3.857914050816557*^9}, {3.857914094834467*^9, 
   3.8579141027472467`*^9}, {3.857914242170517*^9, 3.857914299700845*^9}, 
   3.857914397917735*^9, {3.857914499931308*^9, 3.857914503252205*^9}, {
   3.857916712405786*^9, 3.8579167128010902`*^9}},
 CellLabel->"In[22]:=",ExpressionUUID->"2b6d9722-345b-4c9a-950d-7c4ec7701bb3"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"beltaList", "=", 
   RowBox[{"Power", "[", 
    RowBox[{"10", ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"6", "-", 
        RowBox[{"i", "*", "0.1"}]}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "0", ",", "20"}], "}"}]}], "]"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"deltaList", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"findDelta", "[", 
      RowBox[{"beltaList", "[", 
       RowBox[{"[", "i", "]"}], "]"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "21"}], "}"}]}], "]"}]}], ";"}]}], "Input",
 CellChangeTimes->{
  3.857916939841792*^9},ExpressionUUID->"6d19614d-1674-487b-bcd7-\
8cd9e3c77942"]
},
WindowSize->{808, 758},
WindowMargins->{{211, Automatic}, {-19, Automatic}},
FrontEndVersion->"12.3 for Mac OS X x86 (64-bit) (July 9, 2021)",
StyleDefinitions->"Default.nb",
ExpressionUUID->"9c9c600c-6471-46ea-bd2e-59a2777f58bb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 1194, 26, 30, "Input",ExpressionUUID->"6e53ff02-1f23-4cb5-a63b-5ef983525da5"],
Cell[1755, 48, 815, 12, 52, "Input",ExpressionUUID->"3d9feee8-3117-4c04-aca1-98aed147b460"],
Cell[2573, 62, 574, 9, 30, "Input",ExpressionUUID->"f14827b7-9577-4d0d-992a-2c0c30588c04"],
Cell[3150, 73, 3859, 102, 472, "Input",ExpressionUUID->"29a3d917-1bbf-4528-a5b1-7b5ccc64b2b9"],
Cell[7012, 177, 17695, 452, 1270, "Input",ExpressionUUID->"8f7d00d6-4336-491c-9451-3e91ee437232"],
Cell[24710, 631, 10600, 258, 724, "Input",ExpressionUUID->"557bb1dc-20b5-499c-9529-a8a50b411111"],
Cell[35313, 891, 8743, 183, 472, "Input",ExpressionUUID->"cffde510-1720-4a7f-b240-80202c3a3c96"],
Cell[44059, 1076, 5929, 144, 481, "Input",ExpressionUUID->"da090306-c072-43f6-af23-db60b6ef2112"],
Cell[49991, 1222, 152, 3, 30, "Input",ExpressionUUID->"334b6491-5bd5-421e-b511-41a5b6eecb5f"],
Cell[50146, 1227, 4806, 104, 598, "Input",ExpressionUUID->"2b6d9722-345b-4c9a-950d-7c4ec7701bb3"],
Cell[54955, 1333, 735, 23, 52, "Input",ExpressionUUID->"6d19614d-1674-487b-bcd7-8cd9e3c77942"]
}
]
*)