From f1af831cb8689343a456b4c830c6e5012407d0fa Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Tue, 4 Feb 2020 14:19:04 -0500
Subject: plot fiddling

---
 tempering_analysis.nb | 324 ++++++++++++++++++++++++++++++++++++++++++++------
 1 file changed, 285 insertions(+), 39 deletions(-)

diff --git a/tempering_analysis.nb b/tempering_analysis.nb
index 214f601..4164178 100644
--- a/tempering_analysis.nb
+++ b/tempering_analysis.nb
@@ -10,10 +10,10 @@
 NotebookFileLineBreakTest
 NotebookFileLineBreakTest
 NotebookDataPosition[       158,          7]
-NotebookDataLength[    166178,       3208]
-NotebookOptionsPosition[    162114,       3141]
-NotebookOutlinePosition[    162449,       3156]
-CellTagsIndexPosition[    162406,       3153]
+NotebookDataLength[    181078,       3454]
+NotebookOptionsPosition[    176974,       3387]
+NotebookOutlinePosition[    177309,       3402]
+CellTagsIndexPosition[    177266,       3399]
 WindowFrame->Normal*)
 
 (* Beginning of Notebook Content *)
@@ -1142,27 +1142,31 @@ Cell[CellGroupData[{
 Cell[BoxData[
  RowBox[{"Plot", "[", 
   RowBox[{
-   RowBox[{"{", 
-    RowBox[{
-     RowBox[{"s8", "[", "x", "]"}], ",", 
-     RowBox[{"fit8", "[", "x", "]"}], ",", 
-     RowBox[{"s12", "[", "x", "]"}], ",", 
-     RowBox[{"fit12", "[", "x", "]"}]}], "}"}], ",", 
+   RowBox[{"Evaluate", "[", 
+    RowBox[{"{", 
+     RowBox[{
+      RowBox[{"s8", "[", "x", "]"}], ",", 
+      RowBox[{"fit8", "[", "x", "]"}], ",", 
+      RowBox[{"s12", "[", "x", "]"}], ",", 
+      RowBox[{"fit12", "[", "x", "]"}]}], "}"}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"x", ",", "0", ",", "25"}], "}"}], ",", 
    RowBox[{"AxesLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{"\[Beta]", ",", 
-      RowBox[{
-       RowBox[{"s", "[", "0", "]"}], "-", 
-       RowBox[{"s", "[", "\[Beta]", "]"}]}]}], "}"}]}]}], "]"}]], "Input",
+      RowBox[{"(", 
+       RowBox[{
+        RowBox[{"s", "[", "0", "]"}], "-", 
+        RowBox[{"s", "[", "\[Beta]", "]"}]}], ")"}]}], "}"}]}]}], 
+  "]"}]], "Input",
  CellChangeTimes->{{3.7898191382164*^9, 3.78981915322369*^9}, {
   3.789819203381546*^9, 3.7898192055038*^9}, {3.789819326102756*^9, 
   3.789819347009884*^9}, {3.789819591195232*^9, 3.789819600655766*^9}, {
   3.7898196738325663`*^9, 3.789819677929298*^9}, {3.7898197611734333`*^9, 
   3.78981977649758*^9}, {3.7898322181303864`*^9, 3.789832224367557*^9}, {
-  3.789832474001747*^9, 3.789832488208096*^9}},
- CellLabel->"In[65]:=",ExpressionUUID->"511fd1ec-10b6-4e80-beb1-748a7f18c95e"],
+  3.789832474001747*^9, 3.789832488208096*^9}, {3.7898326615380707`*^9, 
+  3.789832728584391*^9}},
+ CellLabel->"In[73]:=",ExpressionUUID->"511fd1ec-10b6-4e80-beb1-748a7f18c95e"],
 
 Cell[BoxData[
  GraphicsBox[{{{}, {}, 
@@ -1241,7 +1245,7 @@ E5LnnMKip27G7bxEGiS5XggCChYRNbqUa/Yt7dMuWkzDIg49wtBqAw3Ox724
 TSxyGFJH0y006C3owLhvY1GYMe/h8FYaCJT8vvD4FxYRovgGRNpo4FR+/Fnf
 LhZxnT0RXb/viirvzb29/Xu/aAXbdhr8D3MPDy0=
        "]]},
-     Annotation[#, "Charting`Private`Tag$16381#1"]& ], 
+     Annotation[#, "Charting`Private`Tag$22607#1"]& ], 
     TagBox[
      {RGBColor[0.880722, 0.611041, 0.142051], AbsoluteThickness[1.6], Opacity[
       1.], LineBox[CompressedData["
@@ -1292,7 +1296,7 @@ TkBG/MjrbDaG4qjZWdy1BLCSZf5izGDoJ2vKK2fdoY8TXuguYKjBeXhwoI6A
 tpx3JM8lDFlGxIYG1xMgmb9q+/gPhpIG7SUkGwhwLdr39AMHQyE9h7Iq1v2i
 xGdhbQ1DLsyzUraNBPwPLyX/JQ==
        "]]},
-     Annotation[#, "Charting`Private`Tag$16381#2"]& ], 
+     Annotation[#, "Charting`Private`Tag$22607#2"]& ], 
     TagBox[
      {RGBColor[0.560181, 0.691569, 0.194885], AbsoluteThickness[1.6], Opacity[
       1.], LineBox[CompressedData["
@@ -1374,7 +1378,7 @@ qisyz9TTQD/ilbXTj+2+F+qKbzXRIOXeWPHTBRKOkrxUDzbTYOYB/96uJRIW
 Fr370LLt4OSoN0prJFzbO3STs5UGH581EOx+k1D7XXtz6bZ5Xv41SP4vCfXu
 mTCat9HAIvfEi/ZNEjoH+J3Y9ZEGbwpc17a2SPjU4eSl7G3/DzlGjN8=
        "]]},
-     Annotation[#, "Charting`Private`Tag$16381#3"]& ], 
+     Annotation[#, "Charting`Private`Tag$22607#3"]& ], 
     TagBox[
      {RGBColor[0.922526, 0.385626, 0.209179], AbsoluteThickness[1.6], Opacity[
       1.], LineBox[CompressedData["
@@ -1425,7 +1429,7 @@ NBI2WouDSWCm7e1pCpornqB6NeAQHzJUksBe73fBfZXUiMNY+M7N7bMU9P1f
 17P16+wTF5qnwaGgcv2Guu1NODAS6wgOS+t77JWeLF5nsfRV07hf633Vl3gs
 m3GwztFM+8SlIIstZzYTGDjkfXDmrK2t52RJYto6/wf6AQe0
        "]]},
-     Annotation[#, "Charting`Private`Tag$16381#4"]& ]}, {}},
+     Annotation[#, "Charting`Private`Tag$22607#4"]& ]}, {}},
   AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
   Axes->{True, True},
   AxesLabel->{
@@ -1480,8 +1484,250 @@ m3GwztFM+8SlIIstZzYTGDjkfXDmrK2t52RJYto6/wf6AQe0
    3.7898200483861856`*^9, {3.789820594030262*^9, 3.789820707787868*^9}, {
    3.7898224631465263`*^9, 3.7898224693046827`*^9}, 3.7898317351153708`*^9, {
    3.789831769123528*^9, 3.7898317873044233`*^9}, 3.7898320868229*^9, {
-   3.789832219182062*^9, 3.7898322246789417`*^9}, 3.789832488485786*^9},
- CellLabel->"Out[65]=",ExpressionUUID->"8db7b6a9-cb19-4300-b41a-8949de3e6820"]
+   3.789832219182062*^9, 3.7898322246789417`*^9}, 3.789832488485786*^9, {
+   3.7898326687007637`*^9, 3.7898327288395452`*^9}},
+ CellLabel->"Out[73]=",ImageCache->GraphicsData["CompressedBitmap", "\<\
+eJztfAd0XdWZ7g3uQHrykslk3mQmhQCBYOOKTYzBYAKYTkILYIq7jSsY9967
+jWyr2Jbci2y5SHLHVu/FVruSrspV772Xb/7/P3sfrkBXUt7krTUra85alu89
+95z9/b3ss88eP2He5E9mTJg35aMJvxo3Z8KsyVM+mvurZz6dQ6d6fcdi+c5k
+i8Vi+5WFP4M+qj9yTOE/Fsu/yn+9MHHiRKjDuOD5f6Lf7zJ/fsHxrj6Od6HJ
+GoL6uHPG3X92vK5vh+tKZ4xC0dyBQHurfH/W+aVlnz2D4hmD0Fpf3O2lq15H
+yeTBaCpN7O7Sip2TUPrRYNRnX+/u0mqvJSj9cAiqo/Z3d2md7z65tPLy0u4u
+bQq/LJeWHvvkW5eKqC295XyzNVauK3Z5qbshW4ty5NLC9Y9/+9J+HS5tb24y
+Ll3EOmgz9DXO+fUi3WmPo2jeI2itK+zR9eVfPI/imYPQWBjVo+srN05AyaTB
+qE3z6XB9f/nvXsTFxWHRokX4/PPPsWfPHjQ2NhpCr6uTsWpqauT73LlzMWvW
+rF747LPPyDyLi4s7wLBoe8mQvVB9YLmhr/CtHSD7yX/fQ0NDA2bPno3c3Fz5
++eDBg/Dx8TGk3dqK6OhotLS0yPeoqCgsWbIEZWVlfZ1wWO930ND5uWkd4IiS
+8PBwzJs3j/9Z7pK/vbF161a5KiIigj/fpYjOzMwUOUC8zzjHhChwde7b6E0R
+VwW9yPX5Tpjth0uXLsHT09O8vrq6GlOnTkV7e7ulr8LRTFZWVhr22VuZaktm
+kgxesH4wnW2XAZ7p+fj91Pj+/v6abaWk/t/io72m0oBa8hDaGsq/DdVXS8y8
+RYmsvxqXLWfSpElsOU5RiCyUThmJwvl/ZAvuDEVbR15envysrKOP+jksLExs
+k4+YmBi5tK2tTdueY7DVqq1Y/CqKpw9EjfWUxRGxv1J2fHw8Fi9ejC+++AJ7
+9+5lRpS6++D06dNYuHAhWz3WrFmD1NRUQ0VOsKp2z0Xpx4NRHrqmA1ZnyUDf
+U+ftYgQjn/d7fE9j4Hm5J99lVA/v6W0Gvbw1D6Kt2XDrpx1ukwDZS13cVl1u
+2MPih0xNPe0cQ8u/bMYYSkGPoDrxkKWn91Su/QjFUx9FWcCCHt9T67mOMg2R
+d3p8D+/pjYZrJ4SlnJ0PUIZskh/Hdn9bszXGENvaB9BUHN/T29rrqg0BLnqY
+pdHdbdqoyue9gOLZg1B6c56lp/dU75hNMX4o8o4/8e17RKt91PV1Z/YYIvjy
+d2bqf6o7gN5oirlpiGDj/ajLvNTT29rKi8yYUh66uqe3cZAon/1nFM8fgYKz
+47u7TfNWvWseSicOR7b7b9HWWPlH9aOfn594L/9jL7969aplqPrt5MmTcn7y
+5MmS4Di33aV+KywsxNq1ayUxsufn5+dbBslv90rA5vtmzJiBmTNnymebzSb0
+fPnll5JzdMRISEgYpKhWOVTB8N9+ZiLl6Llv3z65jQMOxTaLvo/PT5s2zeG+
+XhIH+SD2zHOKPbnrB90LrN7fyzCGHfeRVi8bAusrjHMa5qO0tFSooRA5TI3E
+iSotLU2+aVlt2rQJt27dkltCQkKwYcMGy2B1/eHDh+Hr6+twfV8ztvORlJQk
+MqyqqhqobnFI+soo+uL8+fNwc3PjvIaCggLMmTOH84C+hbQj8ZtGM29xwkgP
+JNOSqiLl6gdRFrT4UXXp0aNHceHCBflmmNA9juKQg3I4pk+fbhYwlJuEvfLy
+cosex1EkugZYtmwZrFarfNNMOdYfvdTV/J2SkIm3c+dOtgVtKiwa5pPloNO/
+I9nf7475PmhvakTpJyNQvHQ0RZTRdLLd8jD/dhdqa2tVAc0aZM9ITEwc5kCZ
+oyQyMjKYFB2xmcH09HT5pqntzDZIUCI8VjQLgB2GilAtO8daSd/CvsG36ePI
+kSPsF/oWlgnbBonXuMUJI10YRl9lGJWrJ6D00ycoxvwGLZUZDyuRcQXBSmaU
+VatWISgoyDL8HymWPjh16pSEFYbYvHkzD2kZoiAc6js5lDx0rcD65yjI3wyZ
+DBBIDke6lOqMgx86yEOySD8lhboT28U97JRFqu94PKxweIxdu3bJt8fkb8/c
+49NPP3XqHpoHZv/y5cvyTRucQ8EpYyksfYvyDPk2WEmRAzJ3PJplItfyk+7V
+3hQfZFRbmwajyPedzhgeri5lwZK/amJ6KWBWWUBAgAwbGhqKjRs3Guq+iy2Q
+xjp06FCXTBscfNdZ/Svf2cnd3d3lMwc/jpH19fUqngzAnTt3UFJSIr9z3uF7
+s7Oz/z5pUICor6UadxhK176IbI/7KNlWWB5Sv/EwKlXJP+oPVKjsDW9vb4lN
+PKaDaIqKirB+/Xr5iRMuUa48obc0cCtXrpQETdLRt7Bg2Go0ioeHh7ql84pd
+O1FTUxNcXV3FyjlZkimY97FgGWr+/PnSCFI+M+/rjKsfOheU9pPKNR+ibNZY
+ZLv9hvrvM/LjcHU992zMHczI2+VQdWf3Grna5Q8ouTq1w1BsVydOnOjZUH3M
+NqBk96uwH3gQ7S31RtWo/SonJwerV68WB6JMjm3btrGR3tX90GhtQdm0MShf
+/aaEyDrbRceh+0lQ56KJXYAPu90uyiLFWL7jfPz+Sgo1exehbOJI2D0eQPHV
+KQ5S6Cse5pgYd+zYwQavRv2pcZIGHct/Bqib2CSvXbsmv3HnzIZBNKmbfqHY
+Ck8OEHI2n34Bm7zHYOu5QUTgtnMDyac2nX4Orr5LEJYcRH7YKn00X1u5a6qI
+oCbpsGVAJyx1du7uf/B1/ywY/9NpuVf+9kFsbKyEF3ZIjq1GbY6vrbrjhEgr
+mUtanhUXwo5hz4XPyLTG3SVW1Zc+PgV3/9VyefmsZ5B7ZBgKzr36z6zk/zUk
+bUh9uRfCggUL5GRzc7PMgVGu6iJC9lE4NrInr2sbsdn7WTYluYaDfeyhDyQc
+NZclW/p3cv8/i1T/14q4SIRYBM9EcLvEBydcKjm6SrJGEmfjaEdsWggOHvtY
+TCdwzWB4nn4Yt85/SCGr7Q+d3Dryv3FuVBccUPEqts8dCRfUVFx3xcHdio/S
+ZW8hb+oo7PQZJE6w8eQruBx5BW3t7UY58lh3YwxA/aVDxqTs5U3wPfYgtvkM
+lrE2nHwHQfERPBlvjDW8u7HuRltZAUo/Goqa/Ssoko+A7fQ4HLyygMYbxLUE
+Np2chdQc42HIMOfD3aP+Vm2eRlXWE6iM2Co+3ZAXjPwyG/Zc/FBI3OQ9Cvsu
+uKCmvtEgcYjDmNLI6YTVGOovTNYHnKZS8AEU+X9gNtVxtmvE9DMy4tpjr8E/
+LFjKNz4e7Y7IvmhvakDZ1NFE7HSUBS8TQptKk/ToLa2N8AndSOwPln8bjy9F
+fkmljDHI+eia8NqDa4TwpuwYaUFKrs80p32KK21wufimQfiJF+ATeJM1LwM8
+0t3QRJntjgxdd8YFpTfnG1GzIl2P3t7ehhtx7iSaIUL4uqPLYC8yCP9j96Pz
+3GbFsrclqTeXpMhkZcmNWX3VzyVVGdjj+4bQvubYeHjfvEk+Z9D+kPPRv6tu
+bww4J7Q3hl0i2ufK6BTx9eht7a24Eb+PRifayZ7XHF6OnKIqGeMP3Y8uKp3x
+FCqXv4OWyixk77+P3OPj/ibtaUT7a0L76qOv43J4jHpkBTzgfPTvqdtF3kR7
+c0oUykNWGJadH6ZHb2trgX/UFvGXzWdGYqe3BxqajEmE+7sbvR/aKktROvEx
+VG+fhbaGMuR4PoKCsy/JY9oByqgyC8Ox/fxThvCPTkN2YYnhP793GF/85weK
+qJr9K4Xm1rwMVJFgpeC3nrpbDdjQVImDV6eruPEMrkYFGAPe55zg7yuCWwuy
+jJjhtgztrQ3IOzEGOYeHoK2pSsUV1mo7Au8cI10OpX9DsN9/F+nYmAv4rXMI
+TX31ni+kj28tsssTe6a+LHip9rG8snRsOzNeqN925hNU1VUa1P+mu6GJ+jyb
+Qb3rEvmh+OpkMcam4njTh5soBOy5MF/G33JmLGz5ccb4v3Y+/g/V+NV7F8r4
+rTlpaK3Nh93zYeSfeYG02Wwmj1M3j5KlDBcvDUw8yNIyxv+P7sdvLbSLaKp3
+zJEfqhO9RDyVMbuUSfXBzbgYrD3+gtB/Kmgma9sY/1fOx/+RGr/Wc60YTkta
+nKzWKLzwpsSxptJEc/wYax4WexjW43r5ORRVJBnj/3t34/cncy9B2eTHZb6S
+hIKWqizYDz6EfMo8ZE3a526nF2PWzk3YenYEdl4YittZ3tDPmf/NOcqPFRf1
+592NcBMs81aoST6mjGiZsuO+iE8rxqQNXtInMyeXohejubVerv+lA4T41U8U
++TwFXT5vPEXJcTLbxDQZJvQb1Gdd0WOHJ+ZjwiofbDwxRcb2uvE6ymsyZex/
+dU6+RmmKviHk1x4zpk7bGiuQe+xx5Bx6FC3Vdo1yJTwT76+4iG2ntpGQhmG3
+70hY867ILb/oDmWARPyqrTPFXFuyU+RHfoxo33+/PN0idZhO7+WXgPdW+OLo
+1UvYf9UwrYCEbRy4DdX/3DneTxUer1cpmzQSlSv+xtOrhmJSjovsOJ+xG3xP
+SYBTy4ZD4YIZmpAOn7CZgukdMlmCFx8/cw75f9S5ej/jQVK9v5c5d1kW8IVA
+Vid6Kvb6o7a+GXN3foVP1l6CvbASwUm7BW7/tfGm3owx/38sT/m5Vgi13Swc
+Tget+QZqO8WiAp+XxQepnDIpzimqFmrnEdVMPeudbcDFbzRySo3S2hD8P2Z1
+S0f5/kyrlNILT3tVrnpfqJehKOrlHh1BieyPaC5PNSlmn2B17jgRLY5cUH4H
++y49Je6dlm8sy/qJI8VO13Voir69CqYjlT9X5xquHDUKJ28Xc8lBU2mCxJ28
+439Ca12RJvKQf6IQeTUiS8aorM0x7T3JLk+q8GNHIv8762h+pmjhOFi1aao4
+YnNShCaxPuemqD3/9LP82Ph7SuYtrW1YtDcAH672J1utFszahhJ43nhN6vbk
+HF/DJX/Uczqdr8fpKNF/UaRx0cIhkP/RZ01xbaq3uFbBudfQ1lyrKc4vqcHH
+ZKtfuNxCc4vh+PVN5UTxq1TWPIqMwkCD4h/2SP3Ol/V0JPYXiqrmpEgRbdWm
+KRx3NLHVCZ5CbJHvu+xl31Mj3IjKNiLd5STTWGrqC+F2eZwEWC4m+fhBR1r/
+cYuDdDyr83E1rPbsXoteulEV52LQ7Pceh+fvqtPsTluPRQrdSZmmOoqrrNhF
+MeHAtRfR1FIrw6vHrP/IFUb/qs24rU1kLO1PXKB+RMEpsiJyi9BdePEtqhOr
+LZrwyppGTNt0FbO2XUddQ7MO0WTCOiVbHKnuLNT/UptkVbksUSmbMprDkiN6
+VdwewyyprKYSW6GTVyYXiMj2nIk1Dcg/epFApxd8ZUB/zzn0vyloTp5c01R8
+/jLaayq0smqSj0ppmX/6GQqKBRrVzSdeUCOTCkwDa2yuplg4VoxM1yDfdQ78
+f9VtTdFfiWVXrvqAeyANXGe7IC0Qx7aWqsx71en6xhbM2X4DUzZeobjW+LWw
+/YTj4OQvLd0B/7se6tJhY43sjtkc9zVwvf0reXDDk8JNZckG8AA2SeF4y9EI
+s+9ivZwInCD2SdHLQL7XAVnKrv/QHKnnujX7V3HZovEaCyKlX8rxGsjJ8V5F
+8ZHLSYL3VXS2mfXzy+KEyxu3130b6xtc/kobNMVl7hPEC0/u7K1Mipcb5xwe
+KiI21k+1q8ant4S3hXsCJNzll9Zq9LOh0yQ7m5ze4xz9P7TltjRTeTbDbPmJ
+b4t+zNhSkyfWLDVtwAJeNKZI6C85YcIqPyzeF2hGW64IRMdU13QD/58Kgitd
+E/7Uzg7w3PaVfjXHjPW8RFjDc1HKwueKUSv66K13sMf/CZlb4eNu5/C/1hDN
+TdISm2WwAzwPWX1nPyXH38msFa84HqDg2b62H4/SHiaD6nBiNBFdwv/GQfiy
+bIvhD2+Ublzrn2e1uBK37/89P0QT/fdX8FyKcTCbvOEKSirqqY5this5Nidn
+3b8McA7/Ww3f2oJqlwUGPHVkjvAt1TkooI6ShV9yfQaXBgq+H1Lt5fhgpR+W
+uQVJnRCWsk84zy2NtnSH/TuNTXVczb7Fyt9Wclg3sdtb6nnZkWDnHhuFhvxQ
+i54H8QvJEKl7nL+NmoZiyez+0Qu7g73PhCVn81hhhBWSPMUzMxPVZ10liQ+W
+aFoRsUGa6b6KZZapi3esQF8Ky8D58DniambT298BW4LK/RqR3Lv26GZjRTiV
+3W3lxSZia10xzzmqtPEimitsFj1Z1dTcKgJ+nwR9PeaCCDjWdvRbYN9g1Jhp
+6o2GG6ekgy+f8xxaspJ1hcDUVCcclNlOtqyq2250qkX/XEGZkg1r4lp/7PV/
+hjrKN0yD6ucc9EEF2pwQJs+/OU1R1uir2GwuT6Hu4hXF5ng0ldw2RcAdxqT1
+l7Fk/2fCYn55vKU7uD9o8VH3wtlQei8/T/HerxltNRil6ps9uDxsDXl7nV5j
+ac0upwx1XKbTzoUt6g5RL++gvIuq9cYjJbZeCl/mshCOlsVXPhEuOSdSmjKX
+dCbYSrDy8JsyU5VdWCSD9nUOp9dVcXzQc2wVS99ihk37YLXUZfhJaGJInh8g
+p9VJ63qMj8hzkfsKQjTmOLtYHKFXi3K33nDrrHTQrMXG4Iuq+LyXp96UW/5W
+7IfLLRKpDMDTky5+Y7H+5DhMXn9JkjEfvZ1DDtRazElDxaI3DJm6LUN7Y52C
+vEeSYOHFt41IQN0el/5swjrifLZnr0SigLgcw2p6OeCJGz6qNddYhxp348WM
+ii9eQ4s91WRMRJl5CXknjJVtbKsV9mvSM54Kmi4dBTs+t216NrqLxzB6sRXJ
+TmTIbshrTSnkmYxxw80zpty+MuR537HK/NOldWW45e7BKKtqMBjrHGqIgydw
+eSZBhv6n0lRFzLtlnXN5yHIkej4gEOf9/ozmygxcDEoX2XESoa653TnKUG1y
+FLN5roN7cZ49qPc9yOc0EKcMb99xApJw8PcoC1yIVGuSVIPMz+5TMaiua2p1
+DjTckR1easjsrHwPrbnpip17UFqdJg2oN0mMJcf1UVHAYhzxuS4oUzdeRVB8
+rjyganYONcLkqU2edbHB86rPOp99XBVoPV2JXSb8ZNlOm87N3XLaxZlYtfu4
+IK7zDENeSY2lyUAzCB3Qk6kbs9HhJYq8wtf49m1qR2rBFGTJ6wJixQteQfPt
+YDMgcHfOrSMXAeyOTaVJnLnvYn/l6dSjE/D5OjfJJ4cvJXINYWkwKO6nFNjD
+qZueEj1Ki5gctvGWD8pmjjXzbltJvp5Lics4LiIOTdmr6ugBnAiNh08kahZ5
++P5XsHLDFgktfiE26jNbLXWOxDtt4w1i+8lKK14PyC0o/+PFnryE7sSJE50S
+/yc1bHttFWq91kkDxDbCc75kH+ZbDsTb6eCJwoC9JKKXop9doSxoibQpTH+s
+2xPYvmkpZm+9LK0/VU2WGkf6u5440Yuf2YT4WbBe7Ks08PVcz2h1IWV804F4
+Ipnf0qKRHGca9vo/KdVyZV1ub0U1PxiqjNkljRVTneI2CB5bP8WCbedwKTSD
+ip1WS1XPqdblHK9K47V0DhR32uGPcbSXgHMy7SMRgBhpSY3VxOeWxmDHhSHw
+uPoCahtL9coPjqg11pMyj8XEZ7jfj7M738bKLa4481UqRx5LRY9Mpq8C4jWF
+wcHBHaT8DZKN1zLukezCkUPCPZv47vloLTRWllpzL0vAOnDtJapYi5S/cq1K
+qTY3EEU8sU8VChMdsucp7Nq6Akf8Y1BSWW8p60iwk+kfTTAvleGVpBxtnBBs
+vCZ0N9oqilFzYBWZ9TDKT8Ol/KfKVNp0Lql5iqC4MkXHMq5tON1neQ0xLMP1
+YRzb/hEOHD+HxIxSCbbFBq1dzfnoCUAWNC9pdWYJ4xQsxX2zPeOkU3t0CzLT
+fCXMcfWdaD9vrgqhah116ReQ5f26hWPeb2jkQJencXDPWlwPTeZ5iUKDws7S
+wZ81o7Y7Zj8qKcFrOY75vSL+fTFyvtT7ugeVsjZgOWz7B4pQUt0ehM+u93HW
+5wy1x1X5ztGe12hp8bKaQdAmj0LSrr9i78mB5JtjpJ38uucfIG9s1RB/aaff
+kgqMEaP2joS3+0LcCo3iiS1LjnPI8ZpoayyqNk423hilfBe15k9wPzpQZkkK
+KxL1ZEpLTS5yAzch9cBjgpRFiLdcnsOFY9uQmGYX/85yjvaSiRZjPPxhNIqg
+0UuHwtNrkMx/VdXlaTT2+Pq8CFgvzkG6x8OCmO52P67seRNXznvCnl9qyXBA
+k8ruVYXRkm2VhbFiyoRz5/PBOLl/KL66vYFyYvF3Tfk1oiTpPBJOvItMd8Pd
+rK4P4bLrBHx15RRKymssac45ek2htRbnovbQeql/GC1t1mBc3E1oMWtQXV+g
+J6P5WVpOpCfuHHlNZMdoSa4DcdljEoKun0FFZa0lxTnaGwqtraoM9edcUTrz
+KeOFnClD8NWGwbgRtABl1RnfV7e2UIlnC9qL214vmrafTC563eNjBF87ifLK
+6kTnYH/VEmppQmOIL0qWvm4s/floCKKWEpjPe8gsDOJcpx6zDEBTdR6sN3cj
+3nO8yR9L8yuP9xFy5QhKSssst51DvmUqvt142WjnDJSQfTBsOkn05p6ncSfl
+EM9Sqmfqd6OppgjWgL2I83rVhGQjCXB7AwEXdsOenclxKMY56rtaquVFqPLe
+iaKZTxjvLn8yBBGrhiPixlK2yh8qwJb6MqSHeiH2yN9gc3/A9IMw13G4dXIl
+ku9EUiJvjXQO+J46x4VzY9Q1FKx9i/RDkKSf1NmDqaz5C2yZfvx80wDtTVc2
+wH7bDzGnZlD+fVRr847rUNz0morIW6dQWVUV5hxzgoOhFh9aisLpIw1tfjwE
+cStGIvYSuV6t/UcKj12vOCMCsRdWIG7/0wZeH9hIsMFuryLIZysyUu/I0rEg
+56AfaX22tqA28hJy179F+hwiwNnThiB654tIjTtM4mpQT897o6GyAEkBHog+
+9C7p8QHNady+EQg4NAPRzGll+a1vuv0kBd9WVYqis5uRO2+MKVPrvKGI9fwA
+udk32Vo1UhulBnvSdUSfXYh4jyc1j5nu9yHM7QUEnl6BxNhbsgj/unMeJzvY
+bENGHOzuM1EwbYTpKgmLR+HO2fmoKLOqh9G9US88HkDkkQmkzYEa2er2MIL2
+v42QizspnyVINL3sHHmaRiYzqor2RdaWt8lmh5q2m7DqSST7Lkd1VY5CNqad
+irPjSN2bEXXwNVLoAxYNf8d1GAI9P0L4ZVfkZKVK9vZ1Dj9Tx4fGWhRf90DW
+6pfJmgz4gokEv/YZWK+sQ31dkXok3gdtLQ3IvH0VUWeXIHb/sxxxTfh41+EI
+9PoYEVc8kGdPE/hzzuFnK/i22koU+O9ExornTdvKJ/jEteNg9V+Nmqpc9cS4
+D5oba5AefQFR3p8jzmMsu60Jf1txH0Hc2zOSJax5O4efp+Bbq0qQe3ETwT8n
+6tbCT1oxFinnl6K6PF09Cu6DpvoKpEX5IJLgYz2e6QCf4DoYQQffQ8iFbUhP
+DKeGqtlywjn85/KXuviacuRdIu5XjzeFz2SkLBmN5OOzUWLnd2naTQoaqpEe
+44uIMwsRvf8FqoZ/b1KQ6vYQQj1eR/CZlUgkZ62prrIccU7BQu1F9TXIv74X
+tvWvonDSUO1ztnmPIcmNfC7+LLed/6IIaG1pRHbSTURfXI9Iz78Q6MOKgL5E
+zO8Q6fY0Ao98iqgbh5BvT2cXOPhNV1+qYxQlpuKIU0in4ipn+giNbJ8+FElb
+Xobt2jbU1xTpt1W4JSvKvo34a3sQfmQi7riPJI605Q9G8IF3EHx2PZKib6C2
+ttri7pz15ZoAcvlKaxAyDs5A1mdPGAT0kZBqXTIGKSx/WyDL39y/pq6qECnh
+3ojw/gJRpIB00/++5j/oyAxyAndkpsZxX73HOSGrFSGNpXbkXNiA9DXjUTjR
+1EHWp8OQvP11ZNzYifrqwl9qSVC8yLdFIe7qlySJT3DbfZSSRF+kuv4B4R4v
+Iujk54i5eQyFuayDtp3OaVinaGhvbiRtnEb63o+RNfdxLYwiEkbq4ieQTFmq
+MOEi24J6dNkHDbVlsMX5I+rCakR4volkt0EmIcmujyBs/+sIOb0Et4NOo5ja
+JooIW5wTskkR0lCSATuZV9q6l5A/eagmJJdqpJR1zyHNZynKc6JFKwYhxqRp
+ZUkWEoOPI5w0E3ngZbZMTUui6yCEHnhDaIkPOoWivAy6oX29c1q2Kg45NJfe
+uQSb16ewLXqKTcPMSpmzhpN23kDGlS2oKbP9u6KFBkZpnhUJgYcRcfozspIX
+2TU7yuVVhJxagNhbh5GXlcgvWK0yaDG6rB/0eH0Qr7jhPp2ftHNXdvbsWXa5
+zmY+dirxttRXIT/IE2kuHzjqmXedyJg7knqVt5BJ9lZXmaceZSqW8lOREHxM
+ifcVSnuPaJbSqEGKdH8OwWT3kVf2wZYYgvr6OsPNVXP4w79rARG/7b9//35Z
+mK0nR6ir73RGx0X7UFkucq7uQtp2aqRmPab5KqGAaps/EslfvkN87UJdebZ6
+SKrNJhvJ4WcQcW4lwj3fRoLbMIOvftJOxLiORrDXBIScW4c7YefJjO0sDssC
+R+acTjnopxa8AIIv0Qdv0ODt7d0ZO66KnfoiG7L9NyNt619gnzXCkZ2MecTO
+rjeRcXUbqout6qlnX9TXliPjNtVD/tsRdmQSYtzHcoIwLc+NLM/jJQQdm0Nq
+ckVaQjDZU81cR056MvPTVybc2CpZXfz70qVLWX3CjsxXeWgmiomJS1uRuu2v
+rJNeiomsWSOQsvkVpJ9fgbKsMHlAqZloo1IzPyMet28dQhjZWsSB18l/B2ud
+cKEe5zoSoQffIndeSi50FPb022hpbrbM+Hs56SNvxbIR8gu6y5cv5yVMnc5c
+eKnrG0qzkHNtN6y7/4bM+Y+zOixaNTnThiFl9bOwHpmF/FhvtDTWqEeypJrq
+EqTHXxHVhB+dgmiPcew2mqsMt/vI0sYgxJParLOrJUjlZiazX1gm/R2W1kt+
+CgwMNA38m4wcVYw0U+jOC/ZCqvtEpC96kksexUhfFFP5k/7FaKS4EJPXd6Cq
+iPdvaP+tkisn4pL8dCSFkdecX4OwQx8i1n2M1CKaIW7gol2fEtcJZYYCTyI3
+I4HKhxbVyfR07VNCQoJEN36dmg/e8cLJ2qeTirf21maUJFyG7cQCpK55AbnT
+hlu06WXPHE5Z5HmknpiL/HgfKkFq1UNkYoxMr8iehMSQExIOwg69j1i3Mawd
+i6OmYt1GI+Tgu5JOYm96ItMaSfGu1vKOwVhXM2X6oSK/Sslc8dYHbHfNZL6d
+sXRGK5ICYW1BMuxXdiD1y/eQ8flozswdNLbgcSTvfgs2v3XiVFQl/L4DY4lI
+CDmJCJ9VojEODly8GIz1F8e64zacctNrCKEIEXWJrDzmCsqL+ZlL++sGb53l
+yguKp9bmepTE+8J2/HPJ3faZI5TY+yJv8lCkLhtLMegjZH21C1WFPHHRfr+D
+SZUW2pASSbW931aEHp2GSI8X6PKB5I6G5NOosolyfwahRHwYMREfcAT21CgS
+a63lRefUXdKeW5aDvMD9SPOYhPRlTyNvylDTKHKmEXUrqd04MBn2AFfUlKSr
+h9H9qOArRcadW4i97oqwU58hzPMtxLuP4i7TJC3BbSgiPF5GyFEq/Hy3IjH8
+HArsyWLt45yTds1BudU58ci+vA2pe96HbeEYrgI1dfbpw2BdRdQdnGJQV5yu
+Hlz3k2hZaE9BYpgPIi5uJMFNFcEluT2qqOsvaSzObRTCuAI6QbGfNJsU6Sel
+YWtryxPOCbylCGT7qbCFINN3A6xkfQaBQ7RyRXwrqFn0+ISU+yWq8m5zSH9Y
+0dhQV0WKipCIHu6zgqzvI4l/Ka5/NGlkt4pze5xo/AtCjs9DJDWESREXUZCd
+Qr7RNMKBRmkhQk3RtaIiIwxZ/puQ6vI+0pc8ifxJJmWkZKQtfhJWl3eRcWEN
+SpOvo6254RFFWS31fplJweLEYWeXI8TrIzKxZzlTdpDebbcR8qgq5Nh0hF/Y
+iDvBJ8X06muru3ipK9KksQ01Obdhv76LCsFPkLb8GeROH6Zp5AKbM3rKJqpZ
+yfFyww6hvjx7kLa/6nJkJocinmrFcDL7kMMTScPPi4b7KBrJTYTuUK8PEHZm
+MWKuuyE19jq5lJ3DQBdvKcUrGhvL7cgLoaLw8Gykrn8RWbNHcjOuaWQTTF3+
+NFL2TUCG/waUWK+Tu9cZD+H7U1dai1xbDMWXU4j03URymoHwA68hnhI1+YkW
+JU8+GUqmtuDYTESQKG8HnUAWxc/aqrIu3o9JUmRStEZx0mVknF9Nhvgu0heP
+4dBikNnPEOXcx2DdMB5WrxnIvrkHlbmxLAP19L6/PEItzrchNe4qYq66Iez0
+YpLbBKX3QVqk1D0gxu1JhB18m5xmLrG1jdg7jZz0aNTVVnbxGkaaUpyovTgV
+uUEHkHp4FlI3vIRMKrC5hdAEF0ziuP0nWLe8irRj85ATehA1nG2pIBqqCG6k
+/jzHFk8ufhaRftsRenI+BaG3hTqiso+q7sibEE0x3ZFgikM/dU5otg6+1IaX
+pwUZ4YfKgdRVf6ZibQR3BCah+ZOJ0IWjYd36BtJPfI7cEC+T0OGKULbUbFLl
+7WBKnxc3I/TYLPHnGEqhFLyNLrkXWcCfTBJ+7Jy6AkVdS0M1ShIvI/PiOlj3
+TkAaBRo7JfEO1E0aSmIcLWJMperLHuCm9N6iFgT0lh5GKCS/td0J7EIuZQq5
+ub4CRQn+yGDkfR9ShngW2Z8Ol5pPI1MUFOdli0s9MAVZV7ZKkGluqLY47nwr
+QatWD9xYg9KU68ikoGX1mAjrmueROWckm68auL/U+XYCE69jSz+zBHkRxzj4
+s9C72BawSauVqqDKrEjk3KJS/8gcEs5rVNWJv5A4tP0Nhe2zx4n8l4j8qUJ+
+cdIVNNWVWUY7h5C/31E4jeW5KIz2hs1nJax7PiDzeZbCx2PCjcbJofjBcThl
+2+tIOzqbsoQLytKD0ErVVxcvGsvffiY/lISyo5ET6I40ShHW7X9F2pKnKDYN
+VzhKamQcqUvHImXHW3TdfNhv7UNZRihjDewGSr9V2Uwtcqn1BrKuUcQ+9ClS
+tpLoKOflThv2bSjKOCnb/2KwRRG+JOUGmurLu4i58venCookLbaQfXU7Qc00
+oUhi3JQLVO7UYaS5J2Dd+LIYWYbfBhTG+aC+LItt4T+7gfql8k/hKvUmsm98
+SdSSQex4UwTFNk1BSTPGJp05b5SEAU726acWIpcKp3JbMFtuF5FP/v5aobG6
+quyxyAs9LMbLpmFd8xwy5o/iDK0tg6cd2DqEP7JC9oeM8ytREH3K0oWPyt8H
+FVRDZb44Klsvq8u69XWxNvunI5gxo0juJUap73cMPF/vVW3ks94G7XkJ4nC5
+IQe6o+N/yMbbRt/x/R7PHfHSGV5WwD0KX3/u3LlOe19HjP+XyRyevFHbN/IG
+dYJFvVGn/c63sbqdW+HVPrwrmD6YFD8/v56N3rP5Dp5VUyuiZKc37tzS09O/
+nu+QX577ewfuI80st388DcEbcV27ds0g8tW/SwQuLi6sPROGNa9F8krHkbrp
+tnlvDSacD76UR9K7T6mRun/5x9iikjfPZN4cpiUdh5HPU7r+rN/VYp4YwhCs
+cY73SsvKMtYq/EiRxBBs8WxlaiMy89U+3sRIbYnVYfcIjcHXp6SkdMDgcXj/
+LUcMfhdL24LaocycaOFrlXcZGEbRaeGVb2oND0+V8SZkLGjeDocPvoXOm28o
+dOLA5psovA+TUohF7/blZGy96I5P8dZ9fDjZA80c/ubNm7h48WJnHPDB25Hq
+IygoSHYu5YOF4uXlZW6y4SQ8OKLw+2H8Ta+6cTZ0RUUFm9Bd6uIudlozKeTJ
+PAcN36V0yruusCOxlpggPjhmfb2dlVNH04TzToO7d8tbMKbIjh8//o3N2O8R
+e3c8OoNOTk7G9u3bLXptoZON3uQza3jmzJkdDFephEMIwzMZPKR2M5Yv7wms
+BddFOJLvvF8tWcM9Sh3OhlUS0EQ72WhOvnMeIITOiObow3PieodHTlH8+caN
+G0x4j4nmZXz0uyba2bBshmQU2r+c72fWS/jW++XLGIM7GtH169cFgp2I9zjl
+N7LZ5ija9dSIeB9QvYGnpryzYZlJlicpX7uykw205Dujsg86itsIbt/Fli1b
+zMWlTjKAfOaUzTGbl0z8SJHLwU7v6HiXIreTLZjkM59XQbkjCffK/rBqN0c5
+Okkfcp6VcuvWrR8rdPYcDlacNR0d+RvbJ4kO3N3dDe6N1tlcQqj2oVXa+fqc
+siKLxmKJ8NiO1/G4nAZMe1BsyedOUpblO/8Fz2qzMw==\
+\>"],ExpressionUUID->"124892ed-6668-4020-8374-64fc4d7d8b68"]
 }, Open  ]],
 
 Cell[BoxData[
@@ -3177,38 +3423,38 @@ Cell[19807, 520, 30051, 577, 246, "Output",ExpressionUUID->"7b89c7a1-9d18-49a7-b
 Cell[49873, 1100, 732, 18, 55, "Input",ExpressionUUID->"09f155dc-a270-465f-83b4-5452a9d3c6a6"],
 Cell[50608, 1120, 839, 17, 55, "Input",ExpressionUUID->"faf27222-05ad-48db-b551-d6176fd5d567"],
 Cell[CellGroupData[{
-Cell[51472, 1141, 1029, 23, 55, "Input",ExpressionUUID->"511fd1ec-10b6-4e80-beb1-748a7f18c95e"],
-Cell[52504, 1166, 17744, 317, 238, "Output",ExpressionUUID->"8db7b6a9-cb19-4300-b41a-8949de3e6820"]
+Cell[51472, 1141, 1155, 27, 55, "Input",ExpressionUUID->"511fd1ec-10b6-4e80-beb1-748a7f18c95e"],
+Cell[52630, 1170, 32478, 559, 238, 17054, 307, "CachedBoxData", "BoxData", "Output",ExpressionUUID->"124892ed-6668-4020-8374-64fc4d7d8b68"]
 }, Open  ]],
-Cell[70263, 1486, 478, 14, 55, "Input",ExpressionUUID->"5a1efa6e-dfe2-4d60-ae14-3a6c7a1d321c"],
+Cell[85123, 1732, 478, 14, 55, "Input",ExpressionUUID->"5a1efa6e-dfe2-4d60-ae14-3a6c7a1d321c"],
 Cell[CellGroupData[{
-Cell[70766, 1504, 1034, 21, 31, "Input",ExpressionUUID->"90d2bb91-ba31-44f8-b0ca-47994c701e88"],
-Cell[71803, 1527, 1212, 16, 35, "Output",ExpressionUUID->"5f4c18c9-8e83-4549-bb7d-63032de7adc1"]
+Cell[85626, 1750, 1034, 21, 31, "Input",ExpressionUUID->"90d2bb91-ba31-44f8-b0ca-47994c701e88"],
+Cell[86663, 1773, 1212, 16, 35, "Output",ExpressionUUID->"5f4c18c9-8e83-4549-bb7d-63032de7adc1"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[73052, 1548, 601, 15, 31, "Input",ExpressionUUID->"24a1235a-0106-41db-84ea-024b7cb0822c"],
-Cell[73656, 1565, 438, 6, 35, "Output",ExpressionUUID->"f44bda96-870c-49ac-bc9e-91bcd0597cef"]
+Cell[87912, 1794, 601, 15, 31, "Input",ExpressionUUID->"24a1235a-0106-41db-84ea-024b7cb0822c"],
+Cell[88516, 1811, 438, 6, 35, "Output",ExpressionUUID->"f44bda96-870c-49ac-bc9e-91bcd0597cef"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[74131, 1576, 1167, 25, 55, "Input",ExpressionUUID->"57f24fa6-c961-417d-8040-d9f006e7bc42"],
-Cell[75301, 1603, 76969, 1284, 238, "Output",ExpressionUUID->"765eacb6-1851-4d58-84e2-e2ed47100b2f"]
+Cell[88991, 1822, 1167, 25, 55, "Input",ExpressionUUID->"57f24fa6-c961-417d-8040-d9f006e7bc42"],
+Cell[90161, 1849, 76969, 1284, 238, "Output",ExpressionUUID->"765eacb6-1851-4d58-84e2-e2ed47100b2f"]
 }, Open  ]],
-Cell[152285, 2890, 1553, 40, 147, "Input",ExpressionUUID->"8f6a5438-be42-44f3-87d6-da2bc37c171f"],
+Cell[167145, 3136, 1553, 40, 147, "Input",ExpressionUUID->"8f6a5438-be42-44f3-87d6-da2bc37c171f"],
 Cell[CellGroupData[{
-Cell[153863, 2934, 632, 15, 55, "Input",ExpressionUUID->"c4cbaaa3-d693-4c1c-b0a0-8724f52d4d94"],
-Cell[154498, 2951, 1038, 28, 46, "Output",ExpressionUUID->"fa68dd36-a911-4b03-a5e1-774b6ced48f8"]
+Cell[168723, 3180, 632, 15, 55, "Input",ExpressionUUID->"c4cbaaa3-d693-4c1c-b0a0-8724f52d4d94"],
+Cell[169358, 3197, 1038, 28, 46, "Output",ExpressionUUID->"fa68dd36-a911-4b03-a5e1-774b6ced48f8"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[155573, 2984, 1030, 27, 78, "Input",ExpressionUUID->"ea17c51f-a1d7-41fd-a34c-228c2d4ab784"],
-Cell[156606, 3013, 3907, 81, 240, "Output",ExpressionUUID->"734656d0-32f9-40e2-a69d-89571abd59bf"]
+Cell[170433, 3230, 1030, 27, 78, "Input",ExpressionUUID->"ea17c51f-a1d7-41fd-a34c-228c2d4ab784"],
+Cell[171466, 3259, 3907, 81, 240, "Output",ExpressionUUID->"734656d0-32f9-40e2-a69d-89571abd59bf"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[160550, 3099, 347, 6, 31, "Input",ExpressionUUID->"249cb54f-1c35-414f-8168-27f2195185e6"],
-Cell[160900, 3107, 316, 4, 35, "Output",ExpressionUUID->"9b814a6a-324d-4cb2-98b4-7e298a5398a1"]
+Cell[175410, 3345, 347, 6, 31, "Input",ExpressionUUID->"249cb54f-1c35-414f-8168-27f2195185e6"],
+Cell[175760, 3353, 316, 4, 35, "Output",ExpressionUUID->"9b814a6a-324d-4cb2-98b4-7e298a5398a1"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[161253, 3116, 589, 15, 78, "Input",ExpressionUUID->"e66c9322-9282-4dd8-b7bf-4e72fb9ae72c"],
-Cell[161845, 3133, 253, 5, 35, "Output",ExpressionUUID->"a4c4ccfb-7bb8-41fb-a9fe-42202dac8ba9"]
+Cell[176113, 3362, 589, 15, 78, "Input",ExpressionUUID->"e66c9322-9282-4dd8-b7bf-4e72fb9ae72c"],
+Cell[176705, 3379, 253, 5, 35, "Output",ExpressionUUID->"a4c4ccfb-7bb8-41fb-a9fe-42202dac8ba9"]
 }, Open  ]]
 }
 ]
-- 
cgit v1.2.3-54-g00ecf