some analysis work

1 files changed, 289 insertions, 220 deletions

diff --git a/tempering_analysis.nb b/tempering_analysis.nb
index 1b6c9aa..044b448 100644
--- a/tempering_analysis.nb
+++ b/tempering_analysis.nb
@@ -10,10 +10,10 @@
 NotebookFileLineBreakTest
 NotebookFileLineBreakTest
 NotebookDataPosition[       158,          7]
-NotebookDataLength[     29108,        698]
-NotebookOptionsPosition[     27171,        661]
-NotebookOutlinePosition[     27505,        676]
-CellTagsIndexPosition[     27462,        673]
+NotebookDataLength[     33282,        767]
+NotebookOptionsPosition[     31337,        730]
+NotebookOutlinePosition[     31675,        745]
+CellTagsIndexPosition[     31632,        742]
 WindowFrame->Normal*)
 
 (* Beginning of Notebook Content *)
@@ -64,16 +64,7 @@ Cell[BoxData[
    "]"}]}]], "Input",
  CellChangeTimes->{{3.778948155736003*^9, 3.778948220094614*^9}, {
   3.778948251774682*^9, 3.778948445776267*^9}},
- CellLabel->"In[1]:=",ExpressionUUID->"6ca216e4-4a7b-4e5f-a718-f4e70c77d6a2"],
-
-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[15]:=",ExpressionUUID->"3b7ef50e-b362-4da3-b734-818bb17a7c2f"],
+ CellLabel->"In[12]:=",ExpressionUUID->"6ca216e4-4a7b-4e5f-a718-f4e70c77d6a2"],
 
 Cell[BoxData[
  RowBox[{
@@ -93,7 +84,7 @@ Cell[BoxData[
     RowBox[{"InverseErfc", "[", "A", "]"}], "2"]}]}]], "Input",
  CellChangeTimes->{{3.7789451133983593`*^9, 3.778945185413146*^9}, {
   3.778945446992125*^9, 3.778945447478743*^9}},
- CellLabel->"In[3]:=",ExpressionUUID->"f756d73a-4921-45d7-acd1-44fe21438a02"],
+ CellLabel->"In[14]:=",ExpressionUUID->"f756d73a-4921-45d7-acd1-44fe21438a02"],
 
 Cell[BoxData[
  RowBox[{
@@ -113,7 +104,7 @@ Cell[BoxData[
    "]"}]}]], "Input",
  CellChangeTimes->{{3.778948523412203*^9, 3.778948610089134*^9}, {
   3.7789486512266293`*^9, 3.778948654153482*^9}},
- CellLabel->"In[4]:=",ExpressionUUID->"c91e85b2-f90e-40d7-bf19-7a03fd9499bb"],
+ CellLabel->"In[15]:=",ExpressionUUID->"c91e85b2-f90e-40d7-bf19-7a03fd9499bb"],
 
 Cell[BoxData[
  RowBox[{
@@ -139,7 +130,7 @@ Cell[BoxData[
        RowBox[{"dat", "[", "\"\<As\>\"", "]"}]}]}], "}"}]}], "]"}]}]], "Input",\
 
  CellChangeTimes->{{3.7789537137281647`*^9, 3.7789537203664207`*^9}},
- CellLabel->"In[5]:=",ExpressionUUID->"1b068f0a-fa61-487d-893c-5b00f279e758"],
+ CellLabel->"In[16]:=",ExpressionUUID->"1b068f0a-fa61-487d-893c-5b00f279e758"],
 
 Cell[BoxData[
  RowBox[{
@@ -189,7 +180,7 @@ Cell[BoxData[
   3.778948688346896*^9}, {3.7789536224795513`*^9, 3.778953623110043*^9}, {
   3.778953673879937*^9, 3.778953685238799*^9}, {3.778953742582839*^9, 
   3.778953745606895*^9}},
- CellLabel->"In[6]:=",ExpressionUUID->"83b164e2-d1bf-4be5-bf92-d3645cb950ca"],
+ CellLabel->"In[17]:=",ExpressionUUID->"83b164e2-d1bf-4be5-bf92-d3645cb950ca"],
 
 Cell[BoxData[
  RowBox[{
@@ -244,7 +235,7 @@ Cell[BoxData[
   3.779700716133662*^9, 3.7797008200306168`*^9}, {3.779700857824486*^9, 
   3.7797008584386597`*^9}, {3.7797008994884863`*^9, 3.779700901055118*^9}, {
   3.779700939144706*^9, 3.779700939215414*^9}},
- CellLabel->"In[7]:=",ExpressionUUID->"cfeba3ce-dc48-40c8-b8b8-f34cf59a253d"],
+ CellLabel->"In[18]:=",ExpressionUUID->"cfeba3ce-dc48-40c8-b8b8-f34cf59a253d"],
 
 Cell[BoxData[
  RowBox[{
@@ -266,7 +257,16 @@ Cell[BoxData[
           RowBox[{"All", ",", "2"}], "]"}], "]"}], "]"}]}], "}"}], "]"}]}], 
    "]"}]}]], "Input",
  CellChangeTimes->{{3.779700960584425*^9, 3.7797010023114433`*^9}},
- CellLabel->"In[8]:=",ExpressionUUID->"a1830ef0-911d-442d-b5d2-b62c202fce89"],
+ CellLabel->"In[19]:=",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[37]:=",ExpressionUUID->"3b7ef50e-b362-4da3-b734-818bb17a7c2f"],
 
 Cell[CellGroupData[{
 
@@ -275,9 +275,7 @@ Cell[BoxData[
   RowBox[{
    RowBox[{"d", "[", 
     RowBox[{
-     RowBox[{"Select", "[", 
-      RowBox[{
-       RowBox[{"#n\[Beta]", "\[Equal]", "64"}], "&"}], "]"}], ",", 
+     RowBox[{"SortBy", "[", "\"\<n\>\"", "]"}], ",", 
      RowBox[{
       RowBox[{"With", "[", 
        RowBox[{
@@ -293,95 +291,122 @@ Cell[BoxData[
    RowBox[{"Joined", "\[Rule]", "True"}], ",", 
    RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]], "Input",
  CellChangeTimes->CompressedData["
-1:eJwdxVkogwEAB/A5krlTJtLYis2WI02KUZMkVxkxtQcWGqYpV59cL2SON0tR
-yssaNiXKPvHg2DLTjBqtkCixpmwefA9zfP/v4dePp9LKO0NZLBaXhsuF9tZj
-gU9mSxpXYm/tHYEPAk2jeNMh1uErBzWPnUob54TedMRPxoS6vxhvf7OleKzD
-dY83iKcHrJ7de8avTuoFawQHXiy1L3xiSdhY6im9ecbCw1X75izMefoWYLZ8
-kcCVkt5RrP/9FZ7R83fdItyToijFZHtDGW6g8mvw2k4V89e5dh1HiTRGXNKY
-t4sVg6o9vBzQWbAjfJrElLvPitXcCWY9GenCCS2GaxyRUxFipf/x1DMf6t/T
-cVuoPQMXFHUP46Xg8ghutsQ5cJrRw3wz5y7xC32yy7rYMqx5nKLis32yZ9IQ
-xLlbYnYC/luNwReNRdGF9PyXe+apFRMPJw5wM7E2IHH66eUt1S78NjR5i7s+
-ZHf4H8z/9Bs=
+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[24]:=",ExpressionUUID->"73943d9b-70f3-4562-ab84-3881106d7e0c"],
+ CellLabel->"In[46]:=",ExpressionUUID->"73943d9b-70f3-4562-ab84-3881106d7e0c"],
 
 Cell[BoxData[
  GraphicsBox[{{}, {{{}, {}, 
-     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.009166666666666668],
+     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.007333333333333334],
+       AbsoluteThickness[1.6], LineBox[CompressedData["
+1:eJw1yX8s1HEcx3E/MstRuHK5k993x3EnSYzyfZuNobpiLW2t/Go3UUdIUoro
+bKoty2ZSGelO6+iIsNSZULj8TMf5I+t33fy4kB/tVPP+vLfX3nvs6RQvjjpp
+ZGBgQP3b/79+asrRbNTrVnkyte5PVFLkmkwRnYfWUmr5o9uvHtxA6yhxi+5g
+8OtS9CJ1Iir+fEH+XfQKtZaeYC6tqkbrqXLTwbn7wlq0IVR2/EhXJNWjjcH5
++2xVv3Uj2gRCslVT1M9mtCmUPpb/pnW1oTdCS3e+3kb5Ak2DEu0H6zJJB9oC
+suRpQSPCTvRmyMoalhbmdq37iiW0sneZSYp60FagLSyhGex/g7YGpjGDX23a
+h6ZDgOWzXL/j/egtIAwyFNOnVOitoD9mmsGMGkDbQO/1yn2H/AbRDCiu2St8
+qSHeBvdUEZ3Lp4bQtqC0lbFaF4iZILMyn489M4xmgYwShl/QENtBzPNpUVrY
+CHo73DzM5bs8IbaH4KX8nW2sUbQDiJYiOWNiNDhC4Di7ljeEVjrCHx1LEid4
+h90J6hY3rGZmopVO4H92wk0xQLozmAyKk3WMMezOkBjXSStLRYMLtDdMeqe0
+k+4Cp70lFW/nSHeFA7HzfiF73mN3hRi7zyZhl9HAhoRJC+a3TtLZUJ9tZcTT
+k86BcBGPa+Kvxs4B+x2JffQcNHAhuO5i/3Qz6VyYKdJGSOdIdwN2T3gkJ3Ac
+uxuYa7p7Ui+hwR0s3UU+9DrS3SHiV3ac4AvpPGjx6RbFOkxg54HUpji09iga
+PMA3oFV34xrpHjDzsdqWaifdE85VqHYvaEn3hK8B8HDVVoOdD9GMxqu9YWgl
+H3I2USl3MkgXwBG1nD1cQ7oAQi2afPP6SPeCphXFcsMs6V4wVs4ueOoySf0F
+QLMyxA==
+       "]]}, 
+     {RGBColor[0.880722, 0.611041, 0.142051], PointSize[0.007333333333333334],
        AbsoluteThickness[1.6], LineBox[CompressedData["
-1:eJw1yX8s1HEcx/ETF83k5yHhHOe4H+gPrl/q+/4uZiyUa/0Y0oysbjL6pd3W
-roZOIY21uq06TtouKrKS/Pj6kUlEWBx1lx8RaZOZTItq9/68t9fee+zJS86I
-S93AYrGof/v/TTdClbRaC7yVmZTJU9Rg1vHrlpdz0POURB4uOaa5hV6kVG8i
-KtxoNXqZesxyKSgylqFXqXfVmTHnhI/Qa1TxXHxCdF8V2gwK6dhDP5Nr0OZQ
-Vq/cN1RUh2ZDktz31HDLS7QlOLqGhzR4vEZvgnbNh7Rr+c1oa+i5Q28LW2LQ
-NtDaOLdy8mkb2hbo9uXyqNMdJivtYKYpcl21uxNtD7/P5s5eEHWhHeCPQTLt
-F9GNdoQqqcWL2aUetBM8VOmlZ/L60Rzg0DH8zsRBtDNM8DM50gMf0S6QU3E+
-sU8+gnYFRa3YXtChR28Bg75rIbphFO0GiqXy9kLFGHor1MasqJ24n9Du0P1A
-6pxTS+wBhqwWdv6ez2hPqDdeWZe1EXOhuJFJ2AwGk8ELDs6nWgW/QjNeoGZu
-msvFRuw8GDfqZ7bfRzM8CGFEFiqLL9i94erC14aWXWjGG3R1ig10Ouk+sChW
-2rHvke4Dzb7OmsRe0vlA9w5pmlZJ54M297uLg3Acuy8sJzlcCj6MZnxBln6E
-7laSLoA43f6dsTrSBXDjRBSbNUC6H3i8l7LdV0n3g5WxJ2mlvAns/iDLrnHa
-EYlm/CF7I2exPYN0IXRpW6m2UtKFMBqleu5YT7oI8nUjau0Y6SII8wxvMqyR
-LoZqZuDXXt4kdjGU3BWo+2k0SODHRHK8OoV0CVSNJnwbzCU9AGJSspLNKkkP
-AKvbjdyCDtID4dmkTYB2kvRAuCgdDo03m8IeBJV5FUm2XDQTBMKjb2XToVPU
-X9S5NVs=
+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.880722, 0.611041, 0.142051], PointSize[0.009166666666666668],
+     {RGBColor[0.560181, 0.691569, 0.194885], PointSize[0.007333333333333334],
        AbsoluteThickness[1.6], LineBox[CompressedData["
-1:eJw1yX8s1HEcx/E7ThzLb+dHPzgO5/w4LX8ws+/bRmGtH7ejov7x41r0Q4XR
-Vl0WhkK1hNZ2t0w/Fl2IkR9Hd3RuhUl+XXLnKpviosyPdap5f97ba+899mQn
-nxekmdBoNOrf/v+tG6fo5pPC5qx0ast6KvGx4dC4/Dr6O8Uar32jDa9AL1Hi
-KpnkGa0avUI9yRBtOAul6HVqW/thw6OROrSRajAYLMskz9F0aL6x4mR58iXa
-FLjuitqiy81oM4ADa1Ghila0OSQsaPg29q/RTLhyocCHUdmFtoKNuYzfnfQe
-9HZ44CK1og/0om0gpkiiXitRbFlsC9Epqam3BX1oO5hoiko+5/4WbQ/3wnLd
-GnUqtAP8yI9sMWtQox0hMzp40arkHdoJVJv5cYOSQTQLQjezhJw7w2hnaLk5
-15M4Mop2gUaFXYcjfRLtCnuYX/7U359Cu0HwMWWhNkOD3gFdZ4fmlhw+oXeC
-UWcbr6on3gXiyL44Zdg0ejdUi/PUAjmxO4xqRBHxUZ+3DB6g7ui0iulGyz3A
-OLauYrjNYGeDLqTZGCJAy9lQmnv0OKOEdE84YlNe2N9NuifkD49kM5dJ9wJZ
-39WP/d5a7F7gnGJ9dyEeDRxYXSr+KS0gnQMsGHI0viLdG6bFeYOpetK9oWIs
-27zKWofdB5Zlvx6+D0fLfcC78dtUkoh0X1hripU1lZPuC6y8+banbaRzgXlt
-XhqhJZ0LTgf9fFctZrH7QULvAHuKj5b7wdcTrq1VQtJ5kCXcn3Qml3QeSLiJ
-iy9qSPeHW5UXeZou0v0h7UNpjlZHegBYGNodTBh67AFwqUbPm+GgIRCKmcy9
-pvtID4QVs/lTp0WkB8Fsf2yJsYD0ILDJdy0qqyOdD8pMw2akknQ+NEzkUOk6
-PfUX0Ek+tw==
+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.560181, 0.691569, 0.194885], PointSize[0.009166666666666668],
+     {RGBColor[0.922526, 0.385626, 0.209179], PointSize[0.007333333333333334],
        AbsoluteThickness[1.6], LineBox[CompressedData["
-1:eJw1yQ0s1HEcx3EPM0+JRLg4d8cd5zhjlZrW/6th5rHSNbPGJuqOPGx1pZKx
-47pkhdtCSY2tVBhNHnZbd5eHQiQaIeT5LJlMxES1+/6+22ffvfZmxqWeStDT
-0dGh/u3/194XaiMlf7utLpXSeoYak0ufUM0S9CLl2Tg/dCGsEL1C3Widmki0
-e4Beo4K8k9tFehXoTWoipHjysegZepvKEl3RlfJq0Loguns2vFdRj9YH312q
-pNdLDWgDsG2hC/MCm9GGID9XwB1oUKCNIfbyoKzlgBJtCox3Iysp0Wq0GdCD
-FyzFjFa0OchurnrZL7RpnWUBPuvCTUZzB3oPHHsY5/xU+h5tCYKmyNC1E13o
-vVBZTS1W0j6grcDvuubtn/EetDW4Lxey419+RO+DUSo5LehVH9oGVjkRvqFF
-n9C2wORH9tSI+9F2cC+qPZ3bP4CmgXxOZptTOozeD7k53aUxJmNoe2i0zj9t
-rCF2gNHOECo9exxNB/XvXFq8xQTaEURWkuKEYjQwYOmHe/kRs29aqxggvL9x
-cD4EDUzYybhlEyAlnQlb/ApOrJJ0Fuhe2vn+aJ10FkQF1E378SaxO4GA1tlf
-H4NWOYHXNcWWUQHpziCYO/rZQE26MxxvelEj/kk6GzJX7TVBjlPY2RBs2jHo
-FYYGDmzOVym2rpLOgYzooXxxJekuQMllcdJe0l1ALTkclblOuiuIL1bF8enT
-2F1BKEmsO+OPBi5EaLKECUmkc4FeIpgWF5HuBoll4fpFTaS7QZl3WrD5V9J5
-IIpfZi3skM4DfolDNZc1g90daOzZ0Tf+aJU7ONf6dJsISfeA8pHsX3m3SfeA
-Q8MnF1afk86H0NrdYNRJOh/ilX315zWke0JwQJGyy3AWuyeUKwLvCDiz1F84
-ATis
+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.009166666666666668],
+     {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.880722, 0.611041, 0.142051], PointSize[0.009166666666666668],
+     {RGBColor[0.560181, 0.691569, 0.194885], PointSize[0.007333333333333334],
        AbsoluteThickness[1.6]}, 
-     {RGBColor[0.560181, 0.691569, 0.194885], PointSize[0.009166666666666668],
+     {RGBColor[0.922526, 0.385626, 0.209179], PointSize[0.007333333333333334],
        AbsoluteThickness[1.6]}}, {
-     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.009166666666666668],
+     {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.880722, 0.611041, 0.142051], PointSize[0.009166666666666668],
+     {RGBColor[0.560181, 0.691569, 0.194885], PointSize[0.007333333333333334],
        AbsoluteThickness[1.6]}, 
-     {RGBColor[0.560181, 0.691569, 0.194885], PointSize[0.009166666666666668],
-       AbsoluteThickness[1.6]}, {}, {}, {}, {}, {}, {}, {}}, {
-     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.009166666666666668],
+     {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.009166666666666668],
+     {RGBColor[0.880722, 0.611041, 0.142051], PointSize[0.007333333333333334],
        AbsoluteThickness[1.6]}, 
-     {RGBColor[0.560181, 0.691569, 0.194885], PointSize[0.009166666666666668],
-       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},
@@ -393,7 +418,6 @@ ATis
   GridLines->{None, None},
   GridLinesStyle->Directive[
     GrayLevel[0.5, 0.4]],
-  ImageSize->{983., Automatic},
   Method->{
    "OptimizePlotMarkers" -> True, 
     "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
@@ -405,7 +429,7 @@ ATis
          Part[#, 1]], 
         (Identity[#]& )[
          Part[#, 2]]}& )}},
-  PlotRange->{{0.052083333333333336`, 15.875}, {0, 0.6609294624557724}},
+  PlotRange->{{0.052083333333333336`, 15.875}, {0, 0.6579102930533526}},
   PlotRangeClipping->True,
   PlotRangePadding->{{
      Scaled[0.02], 
@@ -427,8 +451,11 @@ ATis
    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}},
- CellLabel->"Out[24]=",ExpressionUUID->"2a29681a-d1fe-4e5f-9729-68c4489a550c"]
+   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}},
+ CellLabel->"Out[46]=",ExpressionUUID->"520bf791-28f4-41af-9c2f-24068e5eb363"]
 }, Open  ]],
 
 Cell[CellGroupData[{
@@ -438,9 +465,7 @@ Cell[BoxData[
   RowBox[{
    RowBox[{"d", "[", 
     RowBox[{
-     RowBox[{"Select", "[", 
-      RowBox[{
-       RowBox[{"#n", "\[Equal]", "20"}], "&"}], "]"}], ",", 
+     RowBox[{"SortBy", "[", "\"\<n\>\"", "]"}], ",", 
      RowBox[{
       RowBox[{"With", "[", 
        RowBox[{
@@ -457,74 +482,119 @@ Cell[BoxData[
    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}},
- CellLabel->"In[28]:=",ExpressionUUID->"60949ea9-18bb-4b5b-868e-6a12e7d013f6"],
+  3.780087398738736*^9}, {3.7800877026590443`*^9, 3.780087726448501*^9}, {
+  3.78498770229777*^9, 3.7849877073353243`*^9}, {3.78498924891437*^9, 
+  3.7849892518156357`*^9}},
+ CellLabel->"In[44]:=",ExpressionUUID->"60949ea9-18bb-4b5b-868e-6a12e7d013f6"],
 
 Cell[BoxData[
  GraphicsBox[{{}, {{{}, {}, 
-     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.009166666666666668],
+     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.007333333333333334],
+       AbsoluteThickness[1.6], LineBox[CompressedData["
+1:eJw10XtIU3EUB3DNTNRi5nO6zdfmY/O1RC2LugdF6WGIK4lgE7EWpmWmm6NE
+HWZJDhMU35kV7Z+StFkDK/PWfM18wzT/EPMZ2nC9SLOxHnh+Fy6XD9/fPef8
+fj+/jCsi6Q4rKyvq3/v/u/18oBrs5haNnTJq20vUPm1j5gOnCrSRijNms+cr
+a9HfqImY3pZsUSv6J5XyKa/tAKjRW1T8SuIxRlc72kKV88+I1UdfoK2hoe5r
+5tB6F9oGmFWpAZXmHrQttP2wdZAIe9F24BykL0rl9aPt4fxi59rO+UG0I5x4
+vTq46jaC3gMWo2bj7LQezYBiztOOldXxbSud4P6sOls1SbwX8mq+1yesj6Kd
+weDoYtSbx9Au8OtPx/4k5hDaFZo7vnTf3SLr3eBtjY3R3ULsDl4uNQkb1thf
+6QHpqunrkUKSM+H5k7IYYdZ7tCcsCOqavB8Re8GCSDPcs4X7U7LAcyEn5Wru
+AJoNPrrSCVsHknPgQrKjbNOX2Bsaew2J5oekvw9ctldYPmrwPMEXxt1EsXIt
+rqd9obQ700MsRYMfmCwjh28E6zD3g6j4aO3FXVgP/OFWT725uQnr0f4wXVe4
+eXOT5FzoV7TIx1g4L82FU8yWtTQBjTkPIg2qsOU7JOcBa2N336GXfZgHQKuu
+/FmhHu+bDoAu2cEZ6TDOA4Ega3q8HMEh/QNh8lxFetQU1oMgmJpclla9wXp0
+EFwLLP1cMETqB0N+9WmuzwyeNx0MCk0ab5ScF/AhfeBIgauQ7J8PlZKq26Zq
+nAcEMHMpsbZMTuYXwO+Z42y5mJxfCCi56mQJh/wfAvfqe+bGtGS+UFDkZ81l
+MMh8ofAqtMhpKZbcTxiUdBdHq0xkf2EwvMjwZBtIHg7mTIl6No70D4e6nIJO
+Nusd5hGQbCrRJPHRdARIc8XtupN66i9WdTxG
+       "]]}, 
+     {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:eJw11HtM01cUB/AKoqwBZBQYrVQohUJ59gFlZo7fUYnoJjiRanAoOhyyhMIS
-QAjhoWxmWwfTITBlEswMDt2gzOoEgQoVBLX0QUEQHQpThqNmILKB1Top557k
-l/5xcz/3e869KeeTzPhPbWg0Wsqbb/E3XLxYRmrHBrOJ/TqRolnrIVX3zPPk
-qLaQOpi6WBPU924fSywuxymtZrFMlPQOa2Pt2Wqq+tRizVBHxnNuv+9zluqz
-rj+nitn0vdKsemqZ1Zun5k5XbVYuNFJL55mpJ2OsrdemLlJWPtVCKSvWRuW+
-9Ttl3a6hQaR+/fhxRQuu20DbF3aDkqNtmG857FugMWNzO/B8O7hqu10656pG
-fyVEmW7uSl64jnntIX0y6VLSoW706FD93+6xsNwe9BxgfkX9DVlrL/Wj1XOE
-87kZ3NJVt9FbBRaTeF3NBg3udwb5ld0Zt+734frb4Js+ZPT5RYueC9ipx1Yb
-n+rwfBfgVDozW90NmJcBpo/2CEpPGKg0q+cKrlmFDrNDBvTcYLpt7z/7OvvR
-c4db11eVdMYY0XMH12Mqp5F5I3rvQH5t2abtowPoeYAisKZr2+MB9Jgwu98t
-49nwEHosiAxzNoZLJtFjgcuOgpvFUjtY8lbDRtP0+pcN0+h5gjlvm0Cd/gd6
-bHBhtDzu3km8NXCeFSPPvDKE978G/m07U/76ZT/O0wvqvrKzGCJIHm+IkSb4
-nmwYxPzesGWuxWO2fRA9DgSaeflp0UacNwfkTeWyzZUGfB8cKHIc8PHar0fP
-B4o7sisf9OrR84HZle5/v/DWo8eFuOTBjj8L9Ohx4dEphX1eoR7754LMeCzW
-MZx4vjDN/YudeFeHni9wH11+MCzVoecHBXJGmWqZDj0/oD+tqihqJ/ftB6LU
-eedxiw49HmiVEtvmFSQfD6ZeZF+eSiL5/CF4rmyhJ4J4/nCRMh4MjSOeP+Qz
-M5jL24kXAKzvqIhmPskXABZVQlZGMfH4wP+6oU7yinh8aM1RpuqOEI8PtvZO
-VRx70m8gePEaNeGJxAuELrPMYeEQ6TcIytP5kY1FxAsCuuFVMBuIFwSKuAMX
-xHuIFwzBuskv451Iv8HQntmlKmWTfCHgbz9SMexP7iMEZtTP6aY8ch8hcMY5
-8aeUKdJvKJhlW7wzfyD5QuFcKUMvDCJeGNAmaPqfGcQLA2/FZ7Uf2OrwPYZB
-zFWTo1uydmk/UwCWUV7WvXHd0nvaKoAOz197RHLs97AAHrrHMtR38HylAGZC
-5Mk2Wux/QgC/NQ2pMtdhHqYQvrnQJTvRqUVPCO9F8drOleL/wWEhxERvam1a
-q0VPCC0BfIGvhnhCUOQ3869JcF5MEbSUvHt/wEA8EeRkpyRdKiGe6M18dp5W
-iYkngnrVSO/nN7ToiUCebLznMUjyiSF+jpPs8SHxxBByoFutohNPDEcjxAl3
-a4knhl1p6W7f9vehJwZD9BOlvkZD/Q9wPq4r
+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.880722, 0.611041, 0.142051], PointSize[0.009166666666666668],
+     {RGBColor[0.922526, 0.385626, 0.209179], PointSize[0.007333333333333334],
        AbsoluteThickness[1.6], LineBox[CompressedData["
-1:eJw90Q0slHEcB/ATksmSq+Q4584d55yTVpHTnl8X0yytNOta4rALa05UitmN
-arI1sZizvLVpLiNmctZVOk4vwnHI2fRmnWRRWNx2i1Xz+/dsz5599t3v+395
-2MmZsfJNNBqN+vv++248k5T41Fx3ed9VasNmKroldNlHWYKep2byXjDyt95D
-L1PzH34ZNw80oFcpi5U/kRjVjLZSwrHiqczSdvQ65d120+Kl6kTbACNvii7T
-a9G2EBcQrCpP6UbbQ7HTktCrpAftAOEVMVKOVo92BLu5K2OJn16ineDhj+ns
-lfuv0c6QfmF0RNrVj94GCzWBSV6NAxsucIEhU+F3g2kIvR2SNGUJjLphtCuY
-juRkyUtH0HR4owo3xSuM6B1wVyyZrPIcRe+Exo7aB8fixtC7oP/4YlhICrEb
-9Gmesc6UjaN3g7ktOXei9x3aHSw1P/v5Z6fRDJC9CgmO59rChj3g81pT6Ee1
-BXNPsPhffFv4aBLNhEiJm/PRVrI/L0ioVNS5OxCzgDnLnW6twfOAN4hiIzLE
-DWidNxjcyrSKVJKzITdGWVmdg/ehY8NprUF2uQUNHLgjr4r66kjmOcCnqtl0
-iQFzH7CyFsfjltA6H1DervwycInkXLC6zCbSaKSfC+G5+9eK0kg/D+qb3g9H
-KMg8D4aUCwdi95HcF1qzn/PoPfj/dL6wnmrfZZ4l/X6wfEgvEUeTfj8YLKiQ
-bkkjOR9spIZvqSfJPB+KdUyu8AYa/EFf23y41kjW94fsazyTJJOsLwAoigxz
-zSK5AIwl9JwnBaQ/ADzSBeeu3yJ5AOyZ2buy2kH6hdCZKhA9/p8LQczyszuo
-If2B8FQ9qPatJ3kgMNoTmgs4pF8E1YrC9SgZ2b8I5Bn5J3q4JA8CVe9vRaea
-zAeB+bwzUxNjoP4Ab0s0ug==
+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.009166666666666668],
+     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.007333333333333334],
        AbsoluteThickness[1.6]}, 
-     {RGBColor[0.880722, 0.611041, 0.142051], PointSize[0.009166666666666668],
+     {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.009166666666666668],
+     {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.009166666666666668],
-       AbsoluteThickness[1.6]}, {}, {}, {}}, {
-     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.009166666666666668],
+     {RGBColor[0.880722, 0.611041, 0.142051], PointSize[0.007333333333333334],
        AbsoluteThickness[1.6]}, 
-     {RGBColor[0.880722, 0.611041, 0.142051], PointSize[0.009166666666666668],
-       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, 0},
+  AxesOrigin->{0.052083333333333336`, 0},
   DisplayFunction->Identity,
   Frame->{{False, False}, {False, False}},
   FrameLabel->{{None, None}, {None, None}},
@@ -543,7 +613,7 @@ zAeB+bwzUxNjoP4Ab0s0ug==
          Part[#, 1]], 
         (Identity[#]& )[
          Part[#, 2]]}& )}},
-  PlotRange->{{0, 19.9}, {0, 2.666747640711386}},
+  PlotRange->{{0.052083333333333336`, 15.875}, {0, 1.8376777821931807`}},
   PlotRangeClipping->True,
   PlotRangePadding->{{
      Scaled[0.02], 
@@ -554,8 +624,10 @@ zAeB+bwzUxNjoP4Ab0s0ug==
  CellChangeTimes->{
   3.7797004888467417`*^9, {3.780087263365243*^9, 3.780087269048139*^9}, {
    3.780087378630178*^9, 3.780087399174458*^9}, {3.780087705097999*^9, 
-   3.780087726696293*^9}},
- CellLabel->"Out[28]=",ExpressionUUID->"1c27332b-7aa0-4f6c-9d7f-d912fe119acd"]
+   3.780087726696293*^9}, {3.784987697054266*^9, 3.7849877075856333`*^9}, 
+   3.7849880628367*^9, 3.784988270798539*^9, 3.7849886153459663`*^9, 
+   3.78498920646163*^9, 3.784989252209024*^9},
+ CellLabel->"Out[44]=",ExpressionUUID->"4b818f0b-341a-496c-8216-10c283d7a417"]
 }, Open  ]],
 
 Cell[BoxData[
@@ -563,59 +635,55 @@ Cell[BoxData[
   RowBox[{"a", "=", 
    RowBox[{
    "Import", "[", 
-    "\"\<~/doc/research/hadamard/code/energies_24_0.200000_20.000000_100.dat\>\
+    "\"\<~/doc/research/hadamard/code/energies_32_0.250000_16.000000_64.dat\>\
 \"", "]"}]}], ";"}]], "Input",
- CellChangeTimes->{{3.779699855009756*^9, 3.779699866112481*^9}},
- CellLabel->"In[8]:=",ExpressionUUID->"bb64ee8f-2ab2-44da-86d3-071e144df4ac"],
+ CellChangeTimes->{{3.779699855009756*^9, 3.779699866112481*^9}, {
+  3.784988180321373*^9, 3.784988193494982*^9}, {3.784988626865487*^9, 
+  3.784988627055429*^9}, {3.784989217849276*^9, 3.7849892179355593`*^9}},
+ CellLabel->"In[40]:=",ExpressionUUID->"bb64ee8f-2ab2-44da-86d3-071e144df4ac"],
 
 Cell[CellGroupData[{
 
 Cell[BoxData[
  RowBox[{"ListPlot", "[", "a", "]"}]], "Input",
  CellChangeTimes->{{3.7796998674408073`*^9, 3.779699868848329*^9}},
- CellLabel->"In[9]:=",ExpressionUUID->"27200d6e-a8c7-4ded-bf66-0ba360a8d0c7"],
+ CellLabel->"In[41]:=",ExpressionUUID->"27200d6e-a8c7-4ded-bf66-0ba360a8d0c7"],
 
 Cell[BoxData[
  GraphicsBox[{{}, {{
-     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.009166666666666668],
-       AbsoluteThickness[1.6], PointBox[{{1., 1.*^7}, {1., 1.*^7}}]}, 
-     {RGBColor[0.880722, 0.611041, 0.142051], PointSize[0.009166666666666668],
-       AbsoluteThickness[1.6], PointBox[CompressedData["
-1:eJw1y38w13cAx/FvrUkWKm7JaX5m5Ed+fP1qystX+NIX3y9fvnw1oqVrd8vO
-OhX+4MZNxbo1scba6vyaHcPXTCzlaCXbMGFHJu3K3He5r/zKj7Tt+r4+d+/7
-3OP9fL8sk1Mjj64XCATv/Xf+/7/6NH7LF6c8/bNmOrQXyLPOzrOvpnWRteZe
-vdpHb0FPxz1JywJtjIyKzvhM02dam8Bo9ESxzJc2w6GiqfsD8bQFLjZfM/n8
-DG0F/+ynf5ddom2wDgOJdvW0LcSNu2rRTdvhXHHITp0JejfSriSPFi3Sjrg+
-8qZaT39Wa2cMRo0PJ1nSLrD3FfVeF9KuSIrw7rIMpt1gcNc1qyaWdkdR2Bej
-8cdpIZRiqYX4DO0BYanM9lQ+7Yl01QfdM8W0F1Tdwbd/Lqe9EXgzrGOtgfZB
-f07V45J2ei8EzY2/Vd6j30G/2UqN1xDti/LWo0uBE/Q+vLyprhtU0/uR8yDZ
-ZmWe9oOr0+Hlupc0IO/TMdbdNPfK2UDq93vSNVu1FvjjxU7zG4Wm7P5IsGrV
-77JiF8Fju6Spdje7CGPJ5Qv73NgDkHIio6rQhz0AJQ0t1YVgP4C1HoUkIpj9
-AAy8TBTDYeyBcLyiFLnK2QMxNtOqc0PJHgQ9Pc2RqiPsQfjk0ccO8gz2YHxU
-XzGe+w17MAIK7psKa9nFONvn4y1qYRfjsm7c9P529hCkmST13OlkD0HFbNcx
-g7vsoXBwNvcY6WEPhUmTuPRQL/tBxIy9tun47+wHEffHG5WSQXYJCttz654P
-s0vg7OM2lj/CHgbVZO/A4gP2MCz/MlQmHmcPhzD8VMfZCfZwfFnfmFn9F3sE
-4nozXOsfs0fgV72Sa5cn2aX4c4Pw8MkprSHFjoW8HR5qvpei6XxJqfofrW9J
-YWPkNH5hmnsZ3t7Yfdtaw70Mcw+XhitnuJchMdFo3GqWexkePTUwuDDHfSRW
-f5r2npnnPhI1q51Pgha5j8SnitO9Rc+5j4TRpe7coSXuozDwoeHKlhXuo/BZ
-ivvV8FXuo2C7N7u84AX3UXDRKNf617iX47TI/Lsgwbx2L4eNsi7CYL3W2XJM
-nj9XnLpB61tyqFzc27bpai2IRqt7p1ypz3003r96rPmkMffR2BzTZvaVKffR
-SE9osP3anPsYbN+cUBJozX0M7iSFfavZxX0M6oZKml634z4GoVWqgR/tuVfg
-h4K0t/QduFcgty1z0c2RewWeqcvG7Zy4V8BmY9I6Q2fuYzFnq9j6hEYsIJhP
-aNnDfSwe5kgN8124j4VFyjbVu67cxyHZzzNL4jbf8S8dmScE
+     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
+      0.011000000000000001`], AbsoluteThickness[1.6], 
+      PointBox[{{1., 100000.}, {1., 100000.}}]}, 
+     {RGBColor[0.880722, 0.611041, 0.142051], PointSize[
+      0.011000000000000001`], AbsoluteThickness[1.6], PointBox[CompressedData["
+1:eJw1yQlIE3AUx/FZWuKBYmFGopk6p2jMI2vT6W9ublM3l1JDPOY8yprHjNSS
+lEaaOrAyDyhRIalwWi4rbR2KShBaZFiLYRkaYcdGrRBKoxP//wePx+d9/fI1
+6QfWMRgM/Nv/d21s8emhd548Dw7hkwe0xoHvU0F2E2t0hOekWmL0ZxO746K5
+pvDU71jizTjTaVDFVQiJvTCAsHclK0nE3sh8kCG95C0n3o7FGdvw8HIa8Q74
+m41WFnsfcQAUQtmSWqwgZkJwtNuDF5dBzEJv4ojGsiWTOASDy3OV7k+ziEMR
+r7Dp7pflEO+Eint98PU3JTEbC/UJnxChIg6Hoo6t9lmgjsC2OcYjZUMecSRa
+GseumHzziaNQ4KRobhui3gV290vmOU4BcTRGS5kd0R6FxLshsrCmb7gdJN6D
+pjCeqZWvJuZA52SSFl2j5kIbvP54rUsxcQyUHu9tj7OpY1Feaz7CvkrNwyuZ
+bdO9j9RxUEdxRFkBJcTxeCjwahMpqQF3e8lceCuxFvjD+1GhmKSdj6n+pkT7
+L7Tz8ZnfNN/jVUp6Am6ft2TM8om1CZDH1Lhai2gX4EJE+WG9jnYB2seXvq7o
+aRfC4MCr75qmXQiZxLJy4gPtibDrVbHq7MtIT4R09ecLqy8xQwRv20ozk0O7
+COZacZ4qjXYxcrP7q0yHaBcjzCFA2XOSdgmONbYX6ztol2BDV6DnZT3tSdg4
+dKtzYJT2JET2jQzrn9GeDLcaF+fKt7QnI1epbT6+THsKmL9GjX0OGtJTYJid
+XyzxJGZIUVl1M6svkHYpXLh+Y6ejaJfh7pvo0Gwh7TJUt/jozqbRngpDEDd1
+fy7tqXA05DjrSmmXw6KYGc+ppl0OnTMYWxto3wvrqkOAa6tm4i9L0ywS
        "]]}}, {
-     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.009166666666666668],
-       AbsoluteThickness[1.6]}, 
-     {RGBColor[0.880722, 0.611041, 0.142051], PointSize[0.009166666666666668],
-       AbsoluteThickness[1.6]}, {}, {}, {}}, {
-     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.009166666666666668],
-       AbsoluteThickness[1.6]}, 
-     {RGBColor[0.880722, 0.611041, 0.142051], PointSize[0.009166666666666668],
-       AbsoluteThickness[1.6]}, {}, {}, {}}}, {{}, {}}},
+     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
+      0.011000000000000001`], AbsoluteThickness[1.6]}, 
+     {RGBColor[0.880722, 0.611041, 0.142051], PointSize[
+      0.011000000000000001`], AbsoluteThickness[1.6]}, {}, {}, {}}, {
+     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
+      0.011000000000000001`], AbsoluteThickness[1.6]}, 
+     {RGBColor[0.880722, 0.611041, 0.142051], PointSize[
+      0.011000000000000001`], AbsoluteThickness[1.6]}, {}, {}, {}}}, {{}, {}}},
   AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
   Axes->{True, True},
   AxesLabel->{None, None},
-  AxesOrigin->{0., -4.400895677366526*^9},
+  AxesOrigin->{0., 0},
   DisplayFunction->Identity,
   Frame->{{False, False}, {False, False}},
   FrameLabel->{{None, None}, {None, None}},
@@ -634,7 +702,7 @@ aNnDfSwe5kgN8124j4VFyjbVu67cxyHZzzNL4jbf8S8dmScE
          Part[#, 1]], 
         (Identity[#]& )[
          Part[#, 2]]}& )}},
-  PlotRange->{{0., 100.}, {-5.689886659390474*^9, -4.468737307999365*^9}},
+  PlotRange->{{0., 64.}, {-1.655204112573*^7, 857321.899835}},
   PlotRangeClipping->True,
   PlotRangePadding->{{
      Scaled[0.02], 
@@ -642,8 +710,9 @@ aNnDfSwe5kgN8124j4VFyjbVu67cxyHZzzNL4jbf8S8dmScE
      Scaled[0.05], 
      Scaled[0.05]}},
   Ticks->{Automatic, Automatic}]], "Output",
- CellChangeTimes->{3.7796998691208344`*^9},
- CellLabel->"Out[9]=",ExpressionUUID->"161e209e-ab11-4d64-99a1-d2cc4b5ffda5"]
+ CellChangeTimes->{3.7796998691208344`*^9, 3.784988194503076*^9, 
+  3.7849886277986927`*^9, 3.784989220788183*^9},
+ CellLabel->"Out[41]=",ExpressionUUID->"96361a65-d311-41a5-b251-e9f9bd1db01c"]
 }, Open  ]],
 
 Cell[CellGroupData[{
@@ -659,8 +728,8 @@ Cell[BoxData["0.25`"], "Output",
  CellLabel->"Out[12]=",ExpressionUUID->"5dbe4cc3-7c89-4c96-870e-625a9bf27876"]
 }, Open  ]]
 },
-WindowSize->{1277, 1076},
-WindowMargins->{{Automatic, 2}, {Automatic, 2}},
+WindowSize->{1916, 1196},
+WindowMargins->{{Automatic, -1916}, {Automatic, 0}},
 FrontEndVersion->"12.0 for Linux x86 (64-bit) (April 8, 2019)",
 StyleDefinitions->"Default.nb"
 ]
@@ -675,30 +744,30 @@ CellTagsIndex->{}
 *)
 (*NotebookFileOutline
 Notebook[{
-Cell[558, 20, 1782, 46, 124, "Input",ExpressionUUID->"6ca216e4-4a7b-4e5f-a718-f4e70c77d6a2"],
-Cell[2343, 68, 337, 7, 31, "Input",ExpressionUUID->"3b7ef50e-b362-4da3-b734-818bb17a7c2f"],
-Cell[2683, 77, 668, 18, 59, "Input",ExpressionUUID->"f756d73a-4921-45d7-acd1-44fe21438a02"],
-Cell[3354, 97, 722, 18, 31, "Input",ExpressionUUID->"c91e85b2-f90e-40d7-bf19-7a03fd9499bb"],
-Cell[4079, 117, 804, 24, 31, "Input",ExpressionUUID->"1b068f0a-fa61-487d-893c-5b00f279e758"],
-Cell[4886, 143, 1753, 48, 55, "Input",ExpressionUUID->"83b164e2-d1bf-4be5-bf92-d3645cb950ca"],
-Cell[6642, 193, 1968, 53, 83, "Input",ExpressionUUID->"cfeba3ce-dc48-40c8-b8b8-f34cf59a253d"],
-Cell[8613, 248, 712, 20, 31, "Input",ExpressionUUID->"a1830ef0-911d-442d-b5d2-b62c202fce89"],
+Cell[558, 20, 1783, 46, 101, "Input",ExpressionUUID->"6ca216e4-4a7b-4e5f-a718-f4e70c77d6a2"],
+Cell[2344, 68, 669, 18, 59, "Input",ExpressionUUID->"f756d73a-4921-45d7-acd1-44fe21438a02"],
+Cell[3016, 88, 723, 18, 31, "Input",ExpressionUUID->"c91e85b2-f90e-40d7-bf19-7a03fd9499bb"],
+Cell[3742, 108, 805, 24, 31, "Input",ExpressionUUID->"1b068f0a-fa61-487d-893c-5b00f279e758"],
+Cell[4550, 134, 1754, 48, 31, "Input",ExpressionUUID->"83b164e2-d1bf-4be5-bf92-d3645cb950ca"],
+Cell[6307, 184, 1969, 53, 52, "Input",ExpressionUUID->"cfeba3ce-dc48-40c8-b8b8-f34cf59a253d"],
+Cell[8279, 239, 713, 20, 31, "Input",ExpressionUUID->"a1830ef0-911d-442d-b5d2-b62c202fce89"],
+Cell[8995, 261, 337, 7, 31, "Input",ExpressionUUID->"3b7ef50e-b362-4da3-b734-818bb17a7c2f"],
 Cell[CellGroupData[{
-Cell[9350, 272, 1317, 33, 39, "Input",ExpressionUUID->"73943d9b-70f3-4562-ab84-3881106d7e0c"],
-Cell[10670, 307, 6446, 123, 630, "Output",ExpressionUUID->"2a29681a-d1fe-4e5f-9729-68c4489a550c"]
+Cell[9357, 272, 1299, 31, 39, "Input",ExpressionUUID->"73943d9b-70f3-4562-ab84-3881106d7e0c"],
+Cell[10659, 305, 8153, 152, 245, "Output",ExpressionUUID->"520bf791-28f4-41af-9c2f-24068e5eb363"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[17153, 435, 979, 25, 39, "Input",ExpressionUUID->"60949ea9-18bb-4b5b-868e-6a12e7d013f6"],
-Cell[18135, 462, 4725, 95, 244, "Output",ExpressionUUID->"1c27332b-7aa0-4f6c-9d7f-d912fe119acd"]
+Cell[18849, 462, 1028, 25, 39, "Input",ExpressionUUID->"60949ea9-18bb-4b5b-868e-6a12e7d013f6"],
+Cell[19880, 489, 7361, 140, 245, "Output",ExpressionUUID->"4b818f0b-341a-496c-8216-10c283d7a417"]
 }, Open  ]],
-Cell[22875, 560, 330, 8, 31, "Input",ExpressionUUID->"bb64ee8f-2ab2-44da-86d3-071e144df4ac"],
+Cell[27256, 632, 476, 10, 31, "Input",ExpressionUUID->"bb64ee8f-2ab2-44da-86d3-071e144df4ac"],
 Cell[CellGroupData[{
-Cell[23230, 572, 207, 3, 31, "Input",ExpressionUUID->"27200d6e-a8c7-4ded-bf66-0ba360a8d0c7"],
-Cell[23440, 577, 3276, 68, 214, "Output",ExpressionUUID->"161e209e-ab11-4d64-99a1-d2cc4b5ffda5"]
+Cell[27757, 646, 208, 3, 31, "Input",ExpressionUUID->"27200d6e-a8c7-4ded-bf66-0ba360a8d0c7"],
+Cell[27968, 651, 2914, 63, 209, "Output",ExpressionUUID->"96361a65-d311-41a5-b251-e9f9bd1db01c"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[26753, 650, 222, 4, 31, "Input",ExpressionUUID->"85d1832b-b409-4be7-b8d2-26df59eff8fb"],
-Cell[26978, 656, 177, 2, 35, "Output",ExpressionUUID->"5dbe4cc3-7c89-4c96-870e-625a9bf27876"]
+Cell[30919, 719, 222, 4, 31, "Input",ExpressionUUID->"85d1832b-b409-4be7-b8d2-26df59eff8fb"],
+Cell[31144, 725, 177, 2, 35, "Output",ExpressionUUID->"5dbe4cc3-7c89-4c96-870e-625a9bf27876"]
 }, Open  ]]
 }
 ]