path: root/tempering_analysis.nb

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

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

(* CreatedBy='Mathematica 12.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     37301,        869]
NotebookOptionsPosition[     35385,        834]
NotebookOutlinePosition[     35723,        849]
CellTagsIndexPosition[     35680,        846]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"getPData", "[", "dir_", "]"}], ":=", 
  RowBox[{
   RowBox[{
    RowBox[{"With", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"data", "=", 
         RowBox[{"Import", "[", "#", "]"}]}], ",", 
        RowBox[{"metadata", "=", 
         RowBox[{"ToExpression", "/@", 
          RowBox[{"Rest", "[", 
           RowBox[{"StringSplit", "[", 
            RowBox[{
             RowBox[{"FileBaseName", "[", "#", "]"}], ",", "\"\<_\>\""}], 
            "]"}], "]"}]}]}]}], "}"}], ",", 
      RowBox[{"<|", 
       RowBox[{
        RowBox[{"\"\<n\>\"", "\[Rule]", 
         RowBox[{"metadata", "[", 
          RowBox[{"[", "1", "]"}], "]"}]}], ",", 
        RowBox[{"\"\<\[Beta]0\>\"", "\[Rule]", 
         RowBox[{"metadata", "[", 
          RowBox[{"[", "2", "]"}], "]"}]}], ",", 
        RowBox[{"\"\<\[Beta]f\>\"", "\[Rule]", 
         RowBox[{"metadata", "[", 
          RowBox[{"[", "3", "]"}], "]"}]}], ",", 
        RowBox[{"\"\<n\[Beta]\>\"", "\[Rule]", 
         RowBox[{"metadata", "[", 
          RowBox[{"[", "4", "]"}], "]"}]}], ",", 
        RowBox[{"\"\<As\>\"", "\[Rule]", 
         RowBox[{
          RowBox[{"Rest", "[", "data", "]"}], "/", 
          RowBox[{"First", "[", 
           RowBox[{"First", "[", "data", "]"}], "]"}]}]}]}], "|>"}]}], "]"}], 
    "&"}], "/@", 
   RowBox[{"FileNames", "[", 
    RowBox[{"FileNameJoin", "[", 
     RowBox[{"{", 
      RowBox[{"dir", ",", "\"\<probs_*.dat\>\""}], "}"}], "]"}], 
    "]"}]}]}]], "Input",
 CellChangeTimes->{{3.778948155736003*^9, 3.778948220094614*^9}, {
   3.778948251774682*^9, 3.778948445776267*^9}, 3.7849915996103563`*^9, {
   3.784991674882406*^9, 3.784991676833317*^9}},
 CellLabel->"In[1]:=",ExpressionUUID->"6ca216e4-4a7b-4e5f-a718-f4e70c77d6a2"],

Cell[BoxData[
 RowBox[{
  RowBox[{"getEData", "[", "dir_", "]"}], ":=", 
  RowBox[{
   RowBox[{
    RowBox[{"With", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"data", "=", 
         RowBox[{"Import", "[", "#", "]"}]}], ",", 
        RowBox[{"metadata", "=", 
         RowBox[{"ToExpression", "/@", 
          RowBox[{"Rest", "[", 
           RowBox[{"StringSplit", "[", 
            RowBox[{
             RowBox[{"FileBaseName", "[", "#", "]"}], ",", "\"\<_\>\""}], 
            "]"}], "]"}]}]}]}], "}"}], ",", 
      RowBox[{"<|", 
       RowBox[{
        RowBox[{"\"\<n\>\"", "\[Rule]", 
         RowBox[{"metadata", "[", 
          RowBox[{"[", "1", "]"}], "]"}]}], ",", 
        RowBox[{"\"\<\[Beta]0\>\"", "\[Rule]", 
         RowBox[{"metadata", "[", 
          RowBox[{"[", "2", "]"}], "]"}]}], ",", 
        RowBox[{"\"\<\[Beta]f\>\"", "\[Rule]", 
         RowBox[{"metadata", "[", 
          RowBox[{"[", "3", "]"}], "]"}]}], ",", 
        RowBox[{"\"\<n\[Beta]\>\"", "\[Rule]", 
         RowBox[{"metadata", "[", 
          RowBox[{"[", "4", "]"}], "]"}]}], ",", 
        RowBox[{"\"\<Es\>\"", "\[Rule]", 
         RowBox[{
          RowBox[{"First", "@", 
           RowBox[{"Rest", "[", "data", "]"}]}], "/", 
          RowBox[{"First", "[", 
           RowBox[{"First", "[", "data", "]"}], "]"}]}]}]}], "|>"}]}], "]"}], 
    "&"}], "/@", 
   RowBox[{"FileNames", "[", 
    RowBox[{"FileNameJoin", "[", 
     RowBox[{"{", 
      RowBox[{"dir", ",", "\"\<energies_*.dat\>\""}], "}"}], "]"}], 
    "]"}]}]}]], "Input",
 CellChangeTimes->{{3.778948155736003*^9, 3.778948220094614*^9}, {
  3.778948251774682*^9, 3.778948445776267*^9}, {3.7849915996103563`*^9, 
  3.784991614545053*^9}, {3.784991671691276*^9, 3.784991674145091*^9}, {
  3.784991874442992*^9, 3.784991874617395*^9}, {3.784992011338725*^9, 
  3.784992012353281*^9}},
 CellLabel->"In[14]:=",ExpressionUUID->"1ffc3385-0105-4d64-94b0-d5fc2c51fc67"],

Cell[BoxData[
 RowBox[{
  RowBox[{"getData", "[", "dir_", "]"}], ":=", 
  RowBox[{"Dataset", "[", 
   RowBox[{
    RowBox[{"Merge", "[", "First", "]"}], "/@", 
    RowBox[{"GatherBy", "[", 
     RowBox[{
      RowBox[{"Join", "[", 
       RowBox[{
        RowBox[{"getPData", "[", "dir", "]"}], ",", 
        RowBox[{"getEData", "[", "dir", "]"}]}], "]"}], ",", " ", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"#n", ",", "#\[Beta]0", ",", "#\[Beta]f", ",", "#n\[Beta]"}], 
        "}"}], "&"}]}], "]"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.784991617461961*^9, 3.784991728472725*^9}},
 CellLabel->"In[3]:=",ExpressionUUID->"10f87413-28fb-4b72-bd1b-eaeef3d02b15"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"c", "[", 
    RowBox[{"\[Beta]i_", ",", "\[Beta]j_"}], "]"}], "[", "A_", "]"}], ":=", 
  RowBox[{"2", 
   FractionBox[
    RowBox[{"(", 
     RowBox[{
      SuperscriptBox["\[Beta]i", "2"], "+", 
      SuperscriptBox["\[Beta]j", "2"]}], ")"}], 
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{"\[Beta]i", "-", "\[Beta]j"}], ")"}], "2"]], 
   SuperscriptBox[
    RowBox[{"InverseErfc", "[", "A", "]"}], "2"]}]}]], "Input",
 CellChangeTimes->{{3.7789451133983593`*^9, 3.778945185413146*^9}, {
  3.778945446992125*^9, 3.778945447478743*^9}},
 CellLabel->"In[4]:=",ExpressionUUID->"f756d73a-4921-45d7-acd1-44fe21438a02"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Ts", "[", "dat_", "]"}], ":=", 
  RowBox[{"Range", "[", 
   RowBox[{
    RowBox[{"dat", "[", "\"\<\[Beta]0\>\"", "]"}], ",", 
    RowBox[{"dat", "[", "\"\<\[Beta]f\>\"", "]"}], ",", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"dat", "[", "\"\<\[Beta]f\>\"", "]"}], "-", 
       RowBox[{"dat", "[", "\"\<\[Beta]0\>\"", "]"}]}], ")"}], "/", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"dat", "[", "\"\<n\[Beta]\>\"", "]"}], "-", "1"}], ")"}]}]}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.778948523412203*^9, 3.778948610089134*^9}, {
  3.7789486512266293`*^9, 3.778948654153482*^9}},
 CellLabel->"In[5]:=",ExpressionUUID->"c91e85b2-f90e-40d7-bf19-7a03fd9499bb"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Cstmp", "[", "dat_", "]"}], ":=", 
  RowBox[{"With", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"ts", "=", 
      RowBox[{"Ts", "[", "dat", "]"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"{", 
         RowBox[{"#", ",", 
          RowBox[{"#", "+", 
           RowBox[{"ts", "[", 
            RowBox[{"[", "2", "]"}], "]"}], "-", 
           RowBox[{"ts", "[", 
            RowBox[{"[", "1", "]"}], "]"}]}]}], "}"}], "&"}], "/@", "ts"}], 
      ",", 
      RowBox[{"Last", "/@", 
       RowBox[{"dat", "[", "\"\<As\>\"", "]"}]}]}], "}"}]}], "]"}]}]], "Input",\

 CellChangeTimes->{{3.7789537137281647`*^9, 3.7789537203664207`*^9}},
 CellLabel->"In[6]:=",ExpressionUUID->"1b068f0a-fa61-487d-893c-5b00f279e758"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Cs", "[", "dat_", "]"}], ":=", 
  RowBox[{"With", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"ts", "=", 
      RowBox[{"Ts", "[", "dat", "]"}]}], "}"}], ",", 
    RowBox[{"MapThread", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{"Mean", "[", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"#1", "[", 
             RowBox[{"[", "1", "]"}], "]"}], ",", 
            RowBox[{"#1", "[", 
             RowBox[{"[", "2", "]"}], "]"}]}], "}"}], "]"}], ",", 
         RowBox[{
          RowBox[{"c", "[", 
           RowBox[{
            RowBox[{"#1", "[", 
             RowBox[{"[", "1", "]"}], "]"}], ",", 
            RowBox[{"#1", "[", 
             RowBox[{"[", "2", "]"}], "]"}]}], "]"}], "[", "#2", "]"}]}], 
        "}"}], "&"}], ",", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{"{", 
           RowBox[{"#", ",", 
            RowBox[{"#", "+", 
             RowBox[{"ts", "[", 
              RowBox[{"[", "2", "]"}], "]"}], "-", 
             RowBox[{"ts", "[", 
              RowBox[{"[", "1", "]"}], "]"}]}]}], "}"}], "&"}], "/@", 
         RowBox[{"Most", "[", "ts", "]"}]}], ",", 
        RowBox[{"Last", "/@", 
         RowBox[{"dat", "[", "\"\<As\>\"", "]"}]}]}], "}"}]}], "]"}]}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.7789479000625353`*^9, 3.7789480157738333`*^9}, {
  3.778948503834503*^9, 3.778948518745183*^9}, {3.778948666921554*^9, 
  3.778948688346896*^9}, {3.7789536224795513`*^9, 3.778953623110043*^9}, {
  3.778953673879937*^9, 3.778953685238799*^9}, {3.778953742582839*^9, 
  3.778953745606895*^9}},
 CellLabel->"In[7]:=",ExpressionUUID->"83b164e2-d1bf-4be5-bf92-d3645cb950ca"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Z0", "[", "n_", "]"}], ":=", 
  RowBox[{
   SuperscriptBox["2", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"n", "-", "1"}], ")"}], 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"n", "/", "4"}], "+", "1"}], ")"}]}]], 
   RowBox[{"Product", "[", 
    RowBox[{
     RowBox[{
      SuperscriptBox["\[Pi]", 
       RowBox[{"k", "/", "2"}]], 
      RowBox[{"Gamma", "[", 
       RowBox[{"k", "/", "2"}], "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"k", ",", "2", ",", "n"}], "}"}]}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.784991222937964*^9, 3.784991266472834*^9}, {
  3.784991403570047*^9, 3.784991404177*^9}, {3.7849915386807756`*^9, 
  3.784991539440803*^9}},
 CellLabel->"In[8]:=",ExpressionUUID->"f81a2199-caac-407b-95bf-2211d7b265a5"],

Cell[BoxData[
 RowBox[{
  RowBox[{"\[CapitalDelta]Ss", "[", "dat_", "]"}], ":=", 
  RowBox[{"With", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{"ts", "=", 
       RowBox[{"1", "/", 
        RowBox[{"Ts", "[", "dat", "]"}]}]}], ",", 
      RowBox[{"cs", "=", 
       RowBox[{"Cs", "[", "dat", "]"}]}]}], "}"}], ",", 
    RowBox[{
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Mean", "[", 
         RowBox[{"1", "/", "#1"}], "]"}], ",", 
        FractionBox[
         RowBox[{"#2", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"#1", "[", 
             RowBox[{"[", "1", "]"}], "]"}], "-", 
            RowBox[{"#1", "[", 
             RowBox[{"[", "2", "]"}], "]"}]}], ")"}]}], 
         RowBox[{"Mean", "[", "#1", "]"}]]}], "}"}], "&"}], "@@@", 
     RowBox[{"Thread", "[", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Rest", "@", 
         RowBox[{"Rest", "[", 
          RowBox[{"FoldList", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"{", 
              RowBox[{
               RowBox[{"#1", "[", 
                RowBox[{"[", "2", "]"}], "]"}], ",", "#2"}], "}"}], "&"}], 
            ",", 
            RowBox[{"{", 
             RowBox[{"0", ",", "0"}], "}"}], ",", "ts"}], "]"}], "]"}]}], ",", 
        RowBox[{
         RowBox[{"Cs", "[", "dat", "]"}], "[", 
         RowBox[{"[", 
          RowBox[{"All", ",", "2"}], "]"}], "]"}]}], "}"}], "]"}]}]}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.779561994673626*^9, 3.7795620342479887`*^9}, {
  3.779700498628436*^9, 3.779700513292344*^9}, {3.7797005747090063`*^9, 
  3.779700575188902*^9}, {3.779700627710291*^9, 3.7797006535654173`*^9}, {
  3.779700716133662*^9, 3.7797008200306168`*^9}, {3.779700857824486*^9, 
  3.7797008584386597`*^9}, {3.7797008994884863`*^9, 3.779700901055118*^9}, {
  3.779700939144706*^9, 3.779700939215414*^9}, {3.784991004834445*^9, 
  3.7849910388807697`*^9}, {3.78499127505939*^9, 3.78499130611341*^9}, {
  3.784991358538706*^9, 3.784991380665337*^9}},
 CellLabel->"In[9]:=",ExpressionUUID->"cfeba3ce-dc48-40c8-b8b8-f34cf59a253d"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Ss", "[", "dat_", "]"}], ":=", 
  RowBox[{"With", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"ss", "=", 
      RowBox[{"\[CapitalDelta]Ss", "[", "dat", "]"}]}], "}"}], ",", 
    RowBox[{"Thread", "[", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"ss", "[", 
        RowBox[{"[", 
         RowBox[{"All", ",", "1"}], "]"}], "]"}], ",", 
       RowBox[{"Accumulate", "[", 
        RowBox[{"ss", "[", 
         RowBox[{"[", 
          RowBox[{"All", ",", "2"}], "]"}], "]"}], "]"}]}], "}"}], "]"}]}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.779700960584425*^9, 3.7797010023114433`*^9}, {
  3.7849910515540743`*^9, 3.784991075832299*^9}, {3.784991369146078*^9, 
  3.7849913764007483`*^9}, {3.784991547938244*^9, 3.784991556568716*^9}, {
  3.7849917900181932`*^9, 3.784991803720964*^9}, {3.784991863633103*^9, 
  3.784991884800845*^9}, {3.784992068826951*^9, 3.7849920762732773`*^9}, {
  3.784992122226707*^9, 3.7849921237610617`*^9}, {3.7849921543708973`*^9, 
  3.784992167089695*^9}, {3.784992239922765*^9, 3.784992242761147*^9}},
 CellLabel->"In[26]:=",ExpressionUUID->"a1830ef0-911d-442d-b5d2-b62c202fce89"],

Cell[BoxData[
 RowBox[{
  RowBox[{"d", "=", 
   RowBox[{"getData", "[", "\"\<~/doc/research/hadamard/code/\>\"", "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.778948447393362*^9, 3.778948476191824*^9}, {
  3.7795618755287123`*^9, 3.779561880399358*^9}},
 CellLabel->"In[34]:=",ExpressionUUID->"3b7ef50e-b362-4da3-b734-818bb17a7c2f"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{
   RowBox[{"d", "[", 
    RowBox[{
     RowBox[{"SortBy", "[", "\"\<n\>\"", "]"}], ",", 
     RowBox[{
      RowBox[{"With", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"L", "=", 
          SuperscriptBox["#n", "2"]}], "}"}], ",", 
        RowBox[{
         RowBox[{
          RowBox[{"{", 
           RowBox[{"#1", ",", 
            RowBox[{"#2", "/", "L"}]}], "}"}], "&"}], "@@@", 
         RowBox[{"Ss", "[", "#", "]"}]}]}], "]"}], "&"}]}], "]"}], ",", 
   RowBox[{"Joined", "\[Rule]", "True"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]], "Input",
 CellChangeTimes->CompressedData["
1:eJwdxVkogwEAB/A5kjmXQiRsxeYWk5ojn1wZyshRHhAapmmuJtcLOV9kKUop
LXfJyrZ4cJtpRm0oJCuxphjle3B9/+/h149dKxXVOzIYDDYFZ/J0lbtcG3Hk
21uFrQVXcqy1l3bjZX3UCD7Xk2PYUHXkt0e9usPxx3JxqwCvfzFTcU+d8RYv
yR/usHhY9YifDKQFS7haK07Vjb9hvlNP4D712pCajfO21sKx38MXFzNFE3Kc
w2/uxorfX94BNWfTFImbAirSsKamOB0Xk/FCPLeRR/9xIp3HbpGSRZxSEreJ
K9prVXjaPqLGeudBDSZNLYdYHNxHr9C4GjGrXHmBXWKyHA6pf26K6LcVLyG4
2lEXihOSGzvx1Pd0Fy5Te+lx0OIN/eWoKeWdZyPOCj3TseR+gPSOsBGPGuU3
jl2JYrLw36wHPi1Jdk+i5lhu6QdmVtnYpy04DEvtfMM7tag834ifO/rNuOGV
uMKyvmNZfqKNMMvC23HG5Ge2kHpBEJKLo+uuhTjzWVmE/wGo7glY
  "],
 CellLabel->"In[35]:=",ExpressionUUID->"73943d9b-70f3-4562-ab84-3881106d7e0c"],

Cell[BoxData[
 GraphicsBox[{{}, {{{}, {}, 
     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.007333333333333334],
       AbsoluteThickness[1.6], LineBox[CompressedData["
1:eJw1yQ0s1HEcx/HrIswajjzHnftzd9y5ZdxM4/9ttpJs8jAqustDDRlTzUpT
1JRqxYaURCs2auO6nplc6elwDqfcXSbjOlYpWm4YVs33990+++61Nyc9P/4w
k8Fg0P/2/6+fjk6PLPAeKcyh122kF/MtqHZmKfoHHXpoSDNhuIr+TTePDozn
TtegzXSgKtzxnU0Depk2DUtEOVFN6DVaNyE0ln1pRW+AEy+44XaadvRG6KqK
KL92/hHaEvJKs5lBx56hraCAappfSOpE28CNhor2uuButC30ys0lIdWv0Juh
pynmtiKtB20H383Hw54nvV13iT20uULI3cL3aAfwvzC7d2e0Cs2C2LoZfQez
D+0IRpGsovpIP9oJtPcsrXavqNFbwM30cn9CqgbtDC6bWjzlEYNoFyia+1Z0
aozYFazVU7RzxhDaDQo6Wm11X4ndgTU5LjGlDqM9QKNfGXY3EHtChnp1znmH
Fr0V5hcGI60bib3ARctnem0aQXuDmTrgtJqCBjb0WgbXxnWhlWwIS+nsLbT7
iJ0Dic0npZFJaCUHHtJVxWfuk+4DtRx3lmKRdB8wbPv55mnsJ+xc2GMuk4lv
opVcUH/Qno2eIJ2CRn5xXJBgFDsFXfXsoXO5aPAFj2TVTFYb6b5wObNVZzVL
uh/UcMUqIU+H3Q+uxPCk0jw08OBo8sqirZx0HjBTVl//MpHOB4nDxfpigR47
HxIVS0+mM9AggIRdk3f+1JAugFssWYi+n3R/4DSrnaIZBuz+wChzFTiEoiEA
MmuXri9lkR4AsrRMaVQ96UJwNIz6ZPWRLoQxZXmLYo10EVBTRpFc9Bm7CB7n
HmQX7UNDIFRK+vXyS6QHQvzpBx2V3aSLgdo+YFw2kS6GbItp+3mPMfov2T42
7A==
       "]]}, 
     {RGBColor[0.880722, 0.611041, 0.142051], PointSize[0.007333333333333334],
       AbsoluteThickness[1.6], LineBox[CompressedData["
1:eJw10X0s1HEcwHFd5PEWylN5uPNwHMft9KDT1e/jYdwqzVA3sTmkTZFWcSrd
9IRqtagmFu1Maqu7cDNOzX43OQ/FTtF5TqLZiM3q0Kal+Xw/22efvfb+88NO
y4nLYJiYmFDr+/9uzBClHO678Lohh9rwNPVcl7daHnMDPU/x6yrCytofoJco
uUrTsiSoQBuppzIx5565Av2HsthnEFRG1KH/UpHi93emta/Qm4B+Oxucpq5H
b4ZfVT8jHDrUaDMQ1bROnfRtRptDyQen/mJFK9oSUgclzEZJG9oa1qLzqFyp
Fs0E6THNQGdbO3orSAN1V/ILdBsutIUhbmWSvrkbbQdnjMqDuXa9aHtILhLG
pofr0dtgtoJ9dHCmH70duipVSenxn9EOcD2GVTT/ZADtCEL7wqZUxSDaCU6p
GUfsbn5BO4PDkL52Nt6AdoHOzBV3hvUQegd4znW97FAS74TpXrb+XdQw2hVa
0rsNCb3EbiDruT8gPjyCdofi/ItRwj5iD7D6nqJnhY9uGFiwn1tKxzahaRYY
ao0FLJcx7Oz1f7Re21OCptlwjs+a2L1IuifcPXt1RiIZx+4J0fSnxDENGrzg
0Qves3qnCexeoKsuVodmo8EbtIeq2aXdpHtDY7Jsl4r1FbsP5PXE1VifRtM+
ENCYKhdrSeeAzbSjs5fFJHYOMG2cQjJ5aPCFHg/NmjiBdF/YsjfU9MBl0v0g
uEidlVFDuh/kTP6OXOghnQtm2SGPBQukc8FbbDq+xfYbdn9wEVgllvPRtD+s
tMtZ5QmkB4DCLW656jzpAfBw0njbvJR0HnAula2olKTzYPEWQxT+kfRAOD5S
2y+YIz0Qst4YRpctp7AHwQmhrqyBg6aDgBaNJ/eFkc6HVanrD5sU0vkgn2gR
MWVT1D8AFzSG
       "]]}, 
     {RGBColor[0.560181, 0.691569, 0.194885], PointSize[0.007333333333333334],
       AbsoluteThickness[1.6], LineBox[CompressedData["
1:eJw1yXlI03EYx/Et+5XmVmrLPHObOrc5DzJK+qPf4wWG9c/Sf/JYTQwVvKY1
QsR0Li2J0gStxCODOROyzBIXuolJzlaKdGjeuoSl4ZGRK7vY833gw8OLN0+e
I03bwWAw6H/7/233kTb6WkZc6Vza5kU6mnGqe5NZhl6mk0yKDFZ4FXqdVsxk
UirqLvo7zSlqCHefb0Zb6fu83g99JzXo33REP6s56Ws7mgmfPR2O5ysfo+0g
KfkXL6/xKZqCNuXIrnPW5+jdUG7cqCyI1aEdIHPDIk439KIdwXNlNY4tMKDZ
8C1e7J1g6Ufvg7MFEwnUwwGbrzhBdaL2Jp05iHaG+j8DqbU+r9AukCIojnZ1
MaL3Qy8/qoueHkZzYFs7JJtsM6EPwB5dxer22Fu0K1gl2Wu8W6Pog1Cqpppa
F9+h3cClh+rOvTGBdoe+2zVTgZxJtAc4C7fiXYeJPSEnvXD8cuEU2gu6p9aG
r+6dRnvDg9lDWlEN8SHIOGFUbHnMoH1A0/Oaw61BAxc01XMlK/azNuu5cGHZ
LI2IQQMPzudntTNVpPOg/YxSnfWCdD6wmqKZ4xuk82FpLCTiaMAcdl/YYLNM
Q8lovS9cq83raagk3Q9AHtlRrCPdD96YO0Wxq6T7gzP3sHzTex67PzATc/lL
p9EggD6VRvtSSboAssubUzX1pAdAmLW0aqee9AAI/TTYurVAuhDSvhiKO+wX
sAth2uuI4YkEDSKICWWrHaWki0BcRgX/KCJdDNrRuLCoFtLFoJaZWh4NkB4I
snvhKZ1m0gOhoE4msqMWsUvgWElFo5sArZeAyWk9dymS9CBIy3jmUCcnPQgU
K9fHg1SkBwM/hTKPaUgPBmHqe2m2kfQQuHTxp5fOQnoIaLscF+6wzfRfY/Mz
SA==
       "]]}, 
     {RGBColor[0.922526, 0.385626, 0.209179], PointSize[0.007333333333333334],
       AbsoluteThickness[1.6], LineBox[CompressedData["
1:eJw1yX8s1HEcx/GTOJyScK6S7s65cxy2mli/vu/Ntf4yJD+alkb6Nd2h4ZpW
GZHJ76LU1NkaJuW6sxibL3Wjc8MunN90RKXMWNcv0495f97ba+899uTFy48n
bmIwGNS//f8bN0J9OgZ/ihkp1IbnqFfWZ25DXg76C+XmH7lZWVOKXqGSzRed
+GFVaDP1ba11XMGpQf+i5peM2UJDLXqdyouRDpmeN6ItIFfvnnP1twptCY7p
oXqlpBltBcl9XxeyKlvQTAjXxdolsNvRtjBWYp6419SBZsGhkhPL+ied6C3Q
GxibLYt7jXaAcKuu1jDQbvjmNpD1T8ul7G60I+h4q+qPqz3o7VChsihPGdKh
naDfJJfbVevRzhCcpORfSu1Du0BWkbq4OnQAzYbAdFbl1gAD2hUuhOTuPbrf
iOZAIcNXqbYeQe8AffH52nYd8U6QHak3KdJG0bsgIH5J2sgeQ7tB0yPNunMz
8W5QdKnp6zCOdodIg24fd4B4D2Q6SZNYMRMbBi702p7kCwbRNBd6cgyKrtBJ
7Dx4+cE1P7EHTfNgZWzxu+LgFHY+xFUPtxW8QNN8MGksuac409g9oEybqlm+
iaY9IPhaxhX1Z9IFIFhzdZgSzWAXgN1EVWthFBo84cEqM+JpLumeMNvgNzqo
Il0Il4sch7snSRcCW3rO/oD9e+wiaHhjYXM4EE2LIEjrXf44jnQvWAmUFKfl
k+4F+rqEG21q0sWgut/BqhsnXQwacU20mGHC7g1aMIlLhGjaGx6+873DDyHd
B7xE+k5KRroPLJrfuj8rJV0Ct34a51ga0iUQ3UKdZRpJ9wWbwaGq7h+k+0KF
y/xCGWcWux848E7PtAShaT8osKaVEVGk+wPTLzO4Pp10f1C1ZMxk352l/gJL
ADfH
       "]]}}, {
     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.007333333333333334],
       AbsoluteThickness[1.6]}, 
     {RGBColor[0.880722, 0.611041, 0.142051], PointSize[0.007333333333333334],
       AbsoluteThickness[1.6]}, 
     {RGBColor[0.560181, 0.691569, 0.194885], PointSize[0.007333333333333334],
       AbsoluteThickness[1.6]}, 
     {RGBColor[0.922526, 0.385626, 0.209179], PointSize[0.007333333333333334],
       AbsoluteThickness[1.6]}}, {
     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.007333333333333334],
       AbsoluteThickness[1.6]}, 
     {RGBColor[0.880722, 0.611041, 0.142051], PointSize[0.007333333333333334],
       AbsoluteThickness[1.6]}, 
     {RGBColor[0.560181, 0.691569, 0.194885], PointSize[0.007333333333333334],
       AbsoluteThickness[1.6]}, 
     {RGBColor[0.922526, 0.385626, 0.209179], PointSize[0.007333333333333334],
       AbsoluteThickness[
      1.6]}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}}, {
     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.007333333333333334],
       AbsoluteThickness[1.6]}, 
     {RGBColor[0.880722, 0.611041, 0.142051], PointSize[0.007333333333333334],
       AbsoluteThickness[1.6]}, 
     {RGBColor[0.560181, 0.691569, 0.194885], PointSize[0.007333333333333334],
       AbsoluteThickness[1.6]}, 
     {RGBColor[0.922526, 0.385626, 0.209179], PointSize[0.007333333333333334],
       AbsoluteThickness[
      1.6]}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}}}, {{}, {}}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0.052083333333333336`, 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  Method->{
   "OptimizePlotMarkers" -> True, 
    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{{0.052083333333333336`, 15.875}, {0, 0.6579102930533526}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.02], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{
  3.778946374557211*^9, 3.7789464676677723`*^9, 3.778946508497999*^9, 
   3.778946604369417*^9, {3.7789466690966673`*^9, 3.778946701500992*^9}, {
   3.778946731994258*^9, 3.778946749349575*^9}, 3.7789467962636337`*^9, 
   3.778946833931172*^9, 3.778947082771483*^9, {3.778947327915792*^9, 
   3.778947345892826*^9}, 3.778948029994212*^9, 3.7789481053151817`*^9, {
   3.778948692166095*^9, 3.7789487038141737`*^9}, 3.778953585064103*^9, 
   3.778953624385591*^9, 3.778953682031535*^9, 3.778953746826598*^9, {
   3.77956194067076*^9, 3.779561968768979*^9}, 3.779624858106914*^9, {
   3.779624924084155*^9, 3.779624933643355*^9}, {3.779625037660088*^9, 
   3.779625061656746*^9}, 3.779699807199407*^9, {3.779700827690935*^9, 
   3.7797008596365232`*^9}, 3.7797009048423357`*^9, 3.779700941258234*^9, {
   3.779701005894755*^9, 3.779701033694045*^9}, 3.78008725076254*^9, 
   3.780087374491346*^9, 3.780087426664627*^9, {3.780087476814011*^9, 
   3.7800874839642763`*^9}, {3.784987528733613*^9, 3.7849875536676188`*^9}, 
   3.7849876962409067`*^9, 3.7849880617038116`*^9, 3.784988269808079*^9, 
   3.7849886106727123`*^9, 3.784989205039175*^9, {3.784989237650683*^9, 
   3.784989241626642*^9}, {3.7849892794456577`*^9, 3.784989301810355*^9}, 
   3.784991010484539*^9, {3.7849910412866077`*^9, 3.7849910778183727`*^9}, 
   3.784991309318139*^9, 3.784991386947817*^9, {3.784991549707535*^9, 
   3.7849915579184523`*^9}, 3.784991736881611*^9, {3.784991879079419*^9, 
   3.784991887673251*^9}, 3.7849919728655863`*^9, 3.784992017379757*^9, 
   3.78499208028115*^9, 3.784992125434637*^9, {3.784992155783885*^9, 
   3.784992168646801*^9}, 3.7849922441034927`*^9, 3.784992833706853*^9},
 CellLabel->"Out[35]=",ExpressionUUID->"4e0edb5f-c69d-4119-94ee-c0c8715ba2f1"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{
   RowBox[{"d", "[", 
    RowBox[{
     RowBox[{"SortBy", "[", "\"\<n\>\"", "]"}], ",", 
     RowBox[{
      RowBox[{"With", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"L", "=", 
          SuperscriptBox["#n", "2"]}], "}"}], ",", 
        RowBox[{
         RowBox[{
          RowBox[{"{", 
           RowBox[{"#1", ",", 
            RowBox[{"#2", "/", "L"}]}], "}"}], "&"}], "@@@", 
         RowBox[{"Cs", "[", "#", "]"}]}]}], "]"}], "&"}]}], "]"}], ",", 
   RowBox[{"Joined", "\[Rule]", "True"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.779700477940461*^9, 3.77970048742832*^9}, {
  3.780087268580113*^9, 3.7800872686230392`*^9}, {3.780087382286021*^9, 
  3.780087398738736*^9}, {3.7800877026590443`*^9, 3.780087726448501*^9}, {
  3.78498770229777*^9, 3.7849877073353243`*^9}, {3.78498924891437*^9, 
  3.7849892518156357`*^9}},
 CellLabel->"In[36]:=",ExpressionUUID->"60949ea9-18bb-4b5b-868e-6a12e7d013f6"],

Cell[BoxData[
 GraphicsBox[{{}, {{{}, {}, 
     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.007333333333333334],
       AbsoluteThickness[1.6], LineBox[CompressedData["
1:eJw10ntI01EUB3BtyTKztJxuWm1zmzof06hpFvk7Ba3HVHQ5yWGmtZUoPaRM
tHD7JywhfFQqJFoyJ6UWlhpp6EozH5mPVmpmuVpCIJZTtllmDzz3Bz9+fPje
e+45lx/32BmZaoWdnR317/3/XX5GqXzzDQ/22Dlq2SZqIEopyVi8ip6mMrS6
6ficErSZ+tbg5PFDWYm2UGLmlpjMJB36J6Xuqt+a5dWAXqIiGz1naiyNaHu4
yWS10Hc8QdPgkGr+LUfwDO0AxqMpu99rO9F0kAqLKk8c70I7Qm7rn6q9wm60
E3iabAFt9b1oZyiP/ehZwiNeBxnZRtqI/PWyNS6Q7+4mP5g3hHaFhV8K02Pr
AHo9aOf5U+kG4g1g3nVFxvDqQ7tB3/a03uyIQTQD5l0KzBOF/Wh3GBrdrbtz
sgftAeGrDkTbrK/QTAgJq3pk30lyFpilp8turcF+NZ4g32ZvnqnG+TRecC38
bH7UFMk3gmln6H3vGbwPzSZoPXzPvUL3Ar0ZnB9+V4gKyX42TAh4OV+5mAMH
Bhmy8MxmzPUcgMT0S0/VJOcCe2HlMI/RgTkXWgaSb6sSX2LuDbVpFxRNkXi+
3htSOzjK5FycB3gwydCpJSyS80BsodmSZkl9PsTFlsYnuGI9PR+aqpdSWG9I
LoBUo8O4yojWC8BQyL9M209yHxhPU1eMhOL/ofeB7ILQ+cpocr4vyOaUlOsH
Ut8XQuTJlrwkXA9+UOMoOtUehuv1fjAQHTdCL8V+QQhl3c/rme1kfiH87hG3
WRkk94c6pkIrZZH+/MG/WHpRHEPuJwDaDIZuyzDpLwBM16Oy/CiyPxDK99ma
HfpJHgjn67LW5qwm+4OAVes1N5FN6geBdZbe/OkL9gMiWEw9Uj2xh8wnAspR
4lT0mcwfDA/uWovfRZD+gyFhTDrJlvRSfwGOLjJQ
       "]]}, 
     {RGBColor[0.880722, 0.611041, 0.142051], PointSize[0.007333333333333334],
       AbsoluteThickness[1.6], LineBox[CompressedData["
1:eJw10HtIU3EUB/DpMJ2p6dLyNXNOnW5TQaGH5e7R1BaGaJI2smYr/wmx2B9p
VDYECcOcli0LwSgEHylJljIr1Eyd5janTueDHo7It1Mwsbke7PwuXC4fvuec
3/ldtvTq6Vx7Go1G/Xv/f23PJJXAmc6M/1RA2WyiyutPcNdkZeglatcCw/zG
VI1ep14znVzEW8/Rm9Rvs8IgyKlHb1PnzYpcHr0VbaVMlsJh5ZM2tB10Om1J
2y3taDqsX25hh4S9RztAQWabXVBvN9oR8psfT4msH9EMqIqr2Lob3Y/eDfLF
JuuxvmG0K3iwQ1UG/gh6D9DuvXDzEY7bLHeH2e9Vm1UNU2gPsBMu90jTZtFM
SHaTDLxa+ILeC0m3IyYiN4g9wVA56PU1idR7wcviMneHrGn0PmAMLjScEk6g
90P10E5agHYM7Q1jxao/fqZRtA/E0s2NjDg92hdGLnn1xhl1aD+YO+eiKssn
9oeT2Wltxl94XzkLLM1F8UcdiQNAc0PFSYnRog+A5ODq2nbJZ5shEJT6Yguf
wvquQJgXK1IvHBjEnA2GYL2ocY3kbEgS85/pK0h/EPy0Kn12WkkeBHqz1VMz
Q3IOsB84z0weRndxoFLEFNXkkvnBUHBFMH72EOkPhh+6HHWekPSHQKxsWuOQ
PYB5CDR9cFe31GE9hIKsdkMbfX0I81AovebvPKVWY86FvKyedOUWOZ8LD+Wa
UX8+mR8G9O4Vfm2+BvMwSBlKPz4hxv8F4WC/0pHBGiXzwyE1Y7Kv9ggaeNDh
bXxXYsLzunjwtE/WqVvGecCHt+LtVX0lyfmw6PqNd19C7i+AO63ZddoUkgvA
eKtJmVpI5kfAUsjFPmsy2T8CYnwTbrJqyP6RMCctWi7nkv0iISqRFd3oSuZH
wSNJgyNTQfqjYL40cfWMrp/6C+N+O/s=
       "]]}, 
     {RGBColor[0.560181, 0.691569, 0.194885], PointSize[0.007333333333333334],
       AbsoluteThickness[1.6], LineBox[CompressedData["
1:eJw10n8slHEcB/AjzaTSFfKzHDp3jhKaaqvnnX6Y7raaurVYzrEa/VAyuqk/
kBRSSa01ZenabWHFlO1q5ZHMlOOccphKv8iNstv6gejH7vN9tmfPXvs8n8/z
+Xy+jyj5aNx+e4FAwP27/z9tVx9Xt0a17naihrP5Eze+63qbU2MpeYyTd55P
f4LrZCtXL/n6dnunlvyDO9qel5rcUU2e4gy94jNvLtWTZznlwLBVLn1ItkNJ
hJmPzdeT52B4vNB7bPoJeS68iveuPnCzmeyIDJ+Sme+CFrITyprOVVkbWsnO
6L2ar4r80kZeAIOa3yusbCe7QDd8N70m+KXNuYuwdOS1Hj0dZCGqGj/+yf1s
IC9G3rC+0m2OibwEDwazna8UdpNdsdHyoWdrBou7IeHn3FeS9WayOwb3uWc6
lX4gL4WnRL/p2rFJsgccou3K24t+kz2RtOfNMlPWR7IXMuNHE7Z+Y/W9scgy
87jOyuwD3cKw2YlEI9kXJ2pVk+c6OsnLMFpjKWx2YfHlkJXfeJtnpPngh+Ou
ZbXZY102835wntqW+q6cDBFa1WdFfXconxfhXlpOiqmfxf0hu+hWoahg+f6Q
JhsztR70fQSgMWuDUH2Q5QdgxFJdcDCW9o1ADMX4XUhQ0Pt8IN7rDKo/PrRf
rIC2znTaUcniK/DYUayUyKkexJDtjnd45EPz8GJEVD6Nyln7guJBUKcYlA31
dP58EPqnk7aM3GdxCQZioi0XTPQ/8BI82xxX+bSfDCkCdsQWzdOw+lJ4NI9O
R9Wy/oJRXVazq/ga6y8YRVcV83UTLF8G+187m0J8aV5ehipNlyw8iMVDIDzl
u9MayeqHwCwwBxrfU78IxfJpjfbQT5YfCs3zy1ExcrbflQjb/O7wyTjqh1+J
12lHWsJvsf2sQtpQROpMMTufVVD0u0QXKLq5v/0/O38=
       "]]}, 
     {RGBColor[0.922526, 0.385626, 0.209179], PointSize[0.007333333333333334],
       AbsoluteThickness[1.6], LineBox[CompressedData["
1:eJw90g0slHEcB3CiNZk5V+pwuHPOyzmnbNeU2vPVpIgdJ0WthJStVGe20mJ2
YauuQ5rRG+W0Gg239Iodk2VDXudSQnO0mPKSvCyl3f/fsz179tn3//v9/s//
efjx5+SJa0xMTJjV+9/TeOmZz6eHV/oz0xijR5nR50NRyTZq4klGYKcsVbYV
E88w6bqpuPicMuJ55sZVBzN/UQXxEiOTx8Wn3tESrzDOqoNbm5ZriU3hIy/T
B2pfE5uhKqTj6Yq8gXgtZmV/8jjTjcTrEJS0kK22aya2wEj0Xk39lRZiS5wI
/sn5yH5HbIV76cM9bqGtxNZ4oNi+X+vYbnQmC/qciZ3SkA5iG5SznkQejewk
ZoOvDbL2PdlFvAHzsb6J3QnUGzGgezPZq+smtoWFOMzxU1w/8Sa8nLorF0cZ
iDfDr1rRxTv1m5iD9cK5wvyZHmI7tAT2OY3lUNujviDU5XIlneeA8UPJal4E
NReVynyJ6iG1I1QDYwNVLLp/J7xNYLbVcWnuDCZrIkXxiuTgQXFN7h+kfG+0
jgfbIkO4nzfN+cjIjalZDCDW8VFi3lqX8I2shwucBhemKlU0d4FnkTTGqpac
JwQYNk8qr06l/QU4smVHmiqC9ndFxq5Y84ZfNHdFvWQw7+seWi/EkkbxpeI7
zYVI3x0hCTtMvh/cwM9q7JtVk/U6N5gWcMOLlbS/Ow6cWe4Qgda7r/4/Z/cJ
Fmm9B3ymOYsez2juAbXmvMEyhtZ7oijkuu5RJs09UVIdzrI8TutFuD1x82L0
/1yEgNhCWx9nWu+FQWHP4yHfNpJ7IdjLXrPwgr6fGLLQY3pVNN2/GGyF9MJc
Mj1fbyzLpIaaDzT3Ruk0t9w0l86XQKK/NDKaQHMJONns3toUMg8+6FQ3/ygp
oPtb9a38pvH77cxffAM2LA==
       "]]}}, {
     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.007333333333333334],
       AbsoluteThickness[1.6]}, 
     {RGBColor[0.880722, 0.611041, 0.142051], PointSize[0.007333333333333334],
       AbsoluteThickness[1.6]}, 
     {RGBColor[0.560181, 0.691569, 0.194885], PointSize[0.007333333333333334],
       AbsoluteThickness[1.6]}, 
     {RGBColor[0.922526, 0.385626, 0.209179], PointSize[0.007333333333333334],
       AbsoluteThickness[1.6]}}, {
     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.007333333333333334],
       AbsoluteThickness[1.6]}, 
     {RGBColor[0.880722, 0.611041, 0.142051], PointSize[0.007333333333333334],
       AbsoluteThickness[1.6]}, 
     {RGBColor[0.560181, 0.691569, 0.194885], PointSize[0.007333333333333334],
       AbsoluteThickness[1.6]}, 
     {RGBColor[0.922526, 0.385626, 0.209179], PointSize[0.007333333333333334],
       AbsoluteThickness[
      1.6]}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}}, {
     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.007333333333333334],
       AbsoluteThickness[1.6]}, 
     {RGBColor[0.880722, 0.611041, 0.142051], PointSize[0.007333333333333334],
       AbsoluteThickness[1.6]}, 
     {RGBColor[0.560181, 0.691569, 0.194885], PointSize[0.007333333333333334],
       AbsoluteThickness[1.6]}, 
     {RGBColor[0.922526, 0.385626, 0.209179], PointSize[0.007333333333333334],
       AbsoluteThickness[
      1.6]}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}}}, {{}, {}}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0.052083333333333336`, 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  Method->{
   "OptimizePlotMarkers" -> True, 
    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{{0.052083333333333336`, 15.875}, {0, 1.8376777821931807`}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.02], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{
  3.7797004888467417`*^9, {3.780087263365243*^9, 3.780087269048139*^9}, {
   3.780087378630178*^9, 3.780087399174458*^9}, {3.780087705097999*^9, 
   3.780087726696293*^9}, {3.784987697054266*^9, 3.7849877075856333`*^9}, 
   3.7849880628367*^9, 3.784988270798539*^9, 3.7849886153459663`*^9, 
   3.78498920646163*^9, 3.784989252209024*^9, 3.7849913833561907`*^9, 
   3.784992834391604*^9},
 CellLabel->"Out[36]=",ExpressionUUID->"8d2214e2-fecc-4e88-ad31-f8070f2e3de8"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{
   RowBox[{"d", "[", 
    RowBox[{
     RowBox[{"Select", "[", 
      RowBox[{
       RowBox[{"KeyExistsQ", "[", 
        RowBox[{"#", ",", "\"\<Es\>\""}], "]"}], "&"}], "]"}], ",", 
     RowBox[{
      RowBox[{"With", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"L", "=", 
          SuperscriptBox["#n", "2"]}], "}"}], ",", 
        RowBox[{"Thread", "[", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"Ts", "[", "#", "]"}], ",", 
           RowBox[{"#Es", "/", 
            SuperscriptBox["#n", "2"]}]}], "}"}], "]"}]}], "]"}], "&"}]}], 
    "]"}], ",", 
   RowBox[{"Joined", "\[Rule]", "True"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.784992259322311*^9, 3.7849923458651857`*^9}, {
  3.784992651707872*^9, 3.784992651985375*^9}, {3.784992846100277*^9, 
  3.784992888545384*^9}},
 CellLabel->"In[40]:=",ExpressionUUID->"ee8c7b42-c7b8-4a5f-a611-21e18e0fcbb2"],

Cell[BoxData[
 GraphicsBox[{{}, {{{}, {}, 
     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
      0.011000000000000001`], AbsoluteThickness[1.6], LineBox[CompressedData["

1:eJw1yX0s1AEcx/FDGnYpDiEPh7tzx7k7tZA//D7rLIfRhppal5iHabOOGmVy
t86GEu0ILSG3rtWK1LhleUhbam0eWjJbs1aJJFft6KZQzX2/22ffvfYOyDqV
kmPLYrHwb///5o0zlU87U+vOfR3c9Humi5PX7jFEnmc4i8YnhWvk78ydrvOG
3ohFq82M+FWK7bYCsoU5lHRfa2kj/2FEHwslxyfILKyvCj4YN8i2iJs2coal
36zegnInYYlCSd4KvW+IpvQi2QFMk6OqpofsBI9bBq1lhszGiWxfryr7Jaud
8daUoDWJyDvwup6t9E4ku2DJkCyLzCe7YkGXxhdryRyoDENl8hayG2rb/A9o
HpLdcaZ5Jt7vOdkDqfcCB8unyDvRP8Yz2M6RPTG7V72fMZO9wOZXOETZmKz2
xs/RB7+Pscm7UN3smDjiQfZBuC4uy5lL9kWRpDsqQUj2Q/5S8WSYjOyPHqkj
VxFJ5oKdXpVrF2O1hguFy6/0z3LqAXip4ipH46kHIFo3tYGD1AMh9xz4EpNG
PRAd4yv6tXTqQRh7YbN8REk9CMUJmp6GTOo81KV0rj7Loc7DjdbLI1fzqfMR
nBW9Ii2gzof5LKu8RkVdALnbmxjLaeoCNJq7J+1KqAejKFyd51pKPRi7f2xX
vCujLsRsv/90n5q6EM31NfaxF6iLcDJqjtVfQV0EPWOXIaikHoLMxAhebzX1
EKznBmUMXaIeihXZozH7WuqheLywJ9b9CnUxhtsHCkJ11MUQ6bL3meuph6HP
5+hdy1XqYejglDY0NlGXwIcpVR++Rl2CQv3tT97XqUthKZ9wb2yhLkXSvJ8x
qZW6DKKby8l9babBv1tsSh4=
       "]]}}, 
    {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.011000000000000001`],
      AbsoluteThickness[1.6]}, {
     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
      0.011000000000000001`], AbsoluteThickness[1.6]}, {}}, {
     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
      0.011000000000000001`], AbsoluteThickness[1.6]}, {}}}, {{}, {}}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, -0.8513855643374515},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  Method->{
   "OptimizePlotMarkers" -> True, 
    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{{0, 16.}, {-0.9882480184724375, -0.8585888513971874}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{{3.784992312221224*^9, 3.784992346103744*^9}, 
   3.7849926522781267`*^9, {3.7849928355405827`*^9, 3.7849928888523283`*^9}},
 CellLabel->"Out[40]=",ExpressionUUID->"0fd34c5b-9bf7-4c9c-b984-ed2f76da4c8d"]
}, Open  ]]
},
WindowSize->{1916, 1196},
WindowMargins->{{Automatic, -1918}, {Automatic, 2}},
FrontEndVersion->"12.0 for Linux x86 (64-bit) (April 8, 2019)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 1786, 46, 78, "Input",ExpressionUUID->"6ca216e4-4a7b-4e5f-a718-f4e70c77d6a2"],
Cell[2347, 68, 1946, 49, 78, "Input",ExpressionUUID->"1ffc3385-0105-4d64-94b0-d5fc2c51fc67"],
Cell[4296, 119, 678, 17, 31, "Input",ExpressionUUID->"10f87413-28fb-4b72-bd1b-eaeef3d02b15"],
Cell[4977, 138, 668, 18, 59, "Input",ExpressionUUID->"f756d73a-4921-45d7-acd1-44fe21438a02"],
Cell[5648, 158, 722, 18, 31, "Input",ExpressionUUID->"c91e85b2-f90e-40d7-bf19-7a03fd9499bb"],
Cell[6373, 178, 804, 24, 31, "Input",ExpressionUUID->"1b068f0a-fa61-487d-893c-5b00f279e758"],
Cell[7180, 204, 1753, 48, 31, "Input",ExpressionUUID->"83b164e2-d1bf-4be5-bf92-d3645cb950ca"],
Cell[8936, 254, 786, 23, 39, "Input",ExpressionUUID->"f81a2199-caac-407b-95bf-2211d7b265a5"],
Cell[9725, 279, 2112, 55, 52, "Input",ExpressionUUID->"cfeba3ce-dc48-40c8-b8b8-f34cf59a253d"],
Cell[11840, 336, 1157, 26, 31, "Input",ExpressionUUID->"a1830ef0-911d-442d-b5d2-b62c202fce89"],
Cell[13000, 364, 337, 7, 31, "Input",ExpressionUUID->"3b7ef50e-b362-4da3-b734-818bb17a7c2f"],
Cell[CellGroupData[{
Cell[13362, 375, 1299, 31, 39, "Input",ExpressionUUID->"73943d9b-70f3-4562-ab84-3881106d7e0c"],
Cell[14664, 408, 8586, 158, 245, "Output",ExpressionUUID->"4e0edb5f-c69d-4119-94ee-c0c8715ba2f1"]
}, Open  ]],
Cell[CellGroupData[{
Cell[23287, 571, 1028, 25, 39, "Input",ExpressionUUID->"60949ea9-18bb-4b5b-868e-6a12e7d013f6"],
Cell[24318, 598, 7415, 141, 245, "Output",ExpressionUUID->"8d2214e2-fecc-4e88-ad31-f8070f2e3de8"]
}, Open  ]],
Cell[CellGroupData[{
Cell[31770, 744, 947, 26, 39, "Input",ExpressionUUID->"ee8c7b42-c7b8-4a5f-a611-21e18e0fcbb2"],
Cell[32720, 772, 2649, 59, 224, "Output",ExpressionUUID->"0fd34c5b-9bf7-4c9c-b984-ed2f76da4c8d"]
}, Open  ]]
}
]
*)