data added

4ac61a29 · himyss · ffa3296f · 4ac61a29
Commit 4ac61a29 authored Aug 18, 2022 by himyss
Show whitespace changes
Inline Side-by-side

Showing with 741 additions and 495 deletions

deutron.nb deutron.nb +741 -495

No files found.
--- a/deutron.nb
+++ b/deutron.nb
@@ -10,10 +10,10 @@
 NotebookFileLineBreakTest
 NotebookFileLineBreakTest
 NotebookDataPosition[       158,          7]
-NotebookDataLength[     95200,       2289]
-NotebookOptionsPosition[     89553,       2191]
-NotebookOutlinePosition[     89981,       2208]
-CellTagsIndexPosition[     89938,       2205]
+NotebookDataLength[    107736,       2535]
+NotebookOptionsPosition[    101858,       2433]
+NotebookOutlinePosition[    102286,       2450]
+CellTagsIndexPosition[    102243,       2447]
 WindowFrame->Normal*)

 (* Beginning of Notebook Content *)
@@ -46,7 +46,7 @@ Cell[BoxData[
   3.866791544874084*^9, {3.867203881456174*^9, 3.867203891918221*^9}, {
   3.86859415931714*^9, 3.868594159600313*^9}, {3.868594205588904*^9, 
   3.8685942118497667`*^9}},
- CellLabel->"In[45]:=",ExpressionUUID->"15f55500-85b0-49bd-b8c0-7ea00f0e2245"],
+ CellLabel->"In[1]:=",ExpressionUUID->"15f55500-85b0-49bd-b8c0-7ea00f0e2245"],

 Cell[CellGroupData[{

@@ -71,7 +71,7 @@ Cell[BoxData[
 CellChangeTimes->{{3.8672039064993887`*^9, 3.867204004483068*^9}, {
  3.8672040350457573`*^9, 3.867204051044759*^9}, {3.868594165302346*^9, 
  3.8685941958762693`*^9}},
- CellLabel->"In[47]:=",ExpressionUUID->"f1c72a75-7ee7-46f9-a8f7-4cd49fd00c7a"],
+ CellLabel->"In[3]:=",ExpressionUUID->"f1c72a75-7ee7-46f9-a8f7-4cd49fd00c7a"],

 Cell[BoxData[
 RowBox[{"{", 
@@ -82,8 +82,9 @@ Cell[BoxData[
  3.868362805408889*^9, 3.868364543237402*^9, 3.868418704643059*^9, 
  3.8684296801198797`*^9, 3.868503326763413*^9, 3.86851076684522*^9, 
  3.868591730457405*^9, 3.8685945085451393`*^9, 3.868595755218032*^9, 
-  3.86859955362208*^9, 3.8695655152900333`*^9, 3.869628019580538*^9},
- CellLabel->"Out[47]=",ExpressionUUID->"787a36b5-5ebb-450b-97f8-6a66d36d74d3"]
+  3.86859955362208*^9, 3.8695655152900333`*^9, 3.869628019580538*^9, 
+  3.8698124214068537`*^9},
+ CellLabel->"Out[3]=",ExpressionUUID->"242a245a-39df-462d-8951-540e89edf802"]
 }, Open  ]],

 Cell[BoxData[
@@ -110,7 +111,7 @@ Cell[BoxData[
   3.866789623085449*^9, {3.8667897202033997`*^9, 3.866789728514553*^9}, {
   3.868594244912619*^9, 3.86859426617033*^9}, {3.868599535129788*^9, 
   3.868599540788164*^9}},
- CellLabel->"In[48]:=",ExpressionUUID->"40ae1148-daf9-4e15-8762-688f07195dad"],
+ CellLabel->"In[4]:=",ExpressionUUID->"40ae1148-daf9-4e15-8762-688f07195dad"],

 Cell[BoxData[
 RowBox[{
@@ -139,7 +140,7 @@ Cell[BoxData[
  3.86678748396309*^9, 3.866787512998473*^9}, {3.868365116728319*^9, 
  3.868365119682788*^9}, {3.868594293307797*^9, 3.868594304549119*^9}, {
  3.868599541924264*^9, 3.8685995475331373`*^9}},
- CellLabel->"In[49]:=",ExpressionUUID->"7a2bcc4d-5447-48d1-ac24-841589704be5"],
+ CellLabel->"In[5]:=",ExpressionUUID->"7a2bcc4d-5447-48d1-ac24-841589704be5"],

 Cell[CellGroupData[{

@@ -148,7 +149,7 @@ Cell[BoxData[
  RowBox[{"Norm1", "[", 
   RowBox[{"k", ",", "q"}], "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.86860390507694*^9, 3.868603917141555*^9}},
- CellLabel->"In[50]:=",ExpressionUUID->"b71993bf-44b9-4b50-90b9-025a29b6f35d"],
+ CellLabel->"In[6]:=",ExpressionUUID->"b71993bf-44b9-4b50-90b9-025a29b6f35d"],

 Cell[BoxData[
 RowBox[{
@@ -160,8 +161,9 @@ Cell[BoxData[
    FractionBox["4", 
     RowBox[{"k", "+", "q"}]]}], ")"}]}]], "Output",
 CellChangeTimes->{{3.868603908548182*^9, 3.8686039184646873`*^9}, 
-   3.869565521238338*^9, 3.869625957956542*^9, 3.869628024261294*^9},
- CellLabel->"Out[50]=",ExpressionUUID->"36448c3c-ed58-4b87-bca0-c6f7044cd75e"]
+   3.869565521238338*^9, 3.869625957956542*^9, 3.869628024261294*^9, 
+   3.869812428243662*^9},
+ CellLabel->"Out[6]=",ExpressionUUID->"c3aed629-a6fc-486e-8683-c6b9ccef00bb"]
 }, Open  ]],

 Cell[BoxData[
@@ -182,7 +184,7 @@ Cell[BoxData[
   3.868591945528491*^9}, {3.8685920653683558`*^9, 3.868592065938489*^9}, {
   3.868592200352151*^9, 3.868592220304749*^9}, {3.868592933204508*^9, 
   3.868592937488298*^9}, {3.868594321890065*^9, 3.868594328109784*^9}},
- CellLabel->"In[51]:=",ExpressionUUID->"ed029352-16a6-45e1-b42a-177020211247"],
+ CellLabel->"In[7]:=",ExpressionUUID->"ed029352-16a6-45e1-b42a-177020211247"],

 Cell[CellGroupData[{

@@ -191,7 +193,7 @@ Cell[BoxData[
  RowBox[{"PsiD", "[", 
   RowBox[{"k", ",", "q", ",", "r"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.8685984536355467`*^9, 3.868598456868614*^9}},
- CellLabel->"In[52]:=",ExpressionUUID->"1508f981-ab16-4888-9c9d-70043439577a"],
+ CellLabel->"In[8]:=",ExpressionUUID->"1508f981-ab16-4888-9c9d-70043439577a"],

 Cell[BoxData[
 FractionBox[
@@ -212,8 +214,9 @@ Cell[BoxData[
    FractionBox["4", 
     RowBox[{"k", "+", "q"}]]}]]]], "Output",
 CellChangeTimes->{3.868598458698038*^9, 3.8685995629995413`*^9, 
-  3.869565525592538*^9, 3.869625962223899*^9, 3.8696280292131033`*^9},
- CellLabel->"Out[52]=",ExpressionUUID->"d7839f75-9aed-4efb-95bf-f0aa05801eb8"]
+  3.869565525592538*^9, 3.869625962223899*^9, 3.8696280292131033`*^9, 
+  3.869812432237939*^9},
+ CellLabel->"Out[8]=",ExpressionUUID->"e4db3020-f165-46b8-bef6-9b39b4225155"]
 }, Open  ]],

 Cell[CellGroupData[{
@@ -276,7 +279,7 @@ Cell[BoxData[
   3.868594547173224*^9, 3.868594552292612*^9}, {3.868598275497266*^9, 
   3.868598275773817*^9}, 3.8685994674679527`*^9, {3.868599511457546*^9, 
   3.868599518423562*^9}},
- CellLabel->"In[53]:=",ExpressionUUID->"cd919aeb-3101-4893-8360-ca90a34b5f8a"],
+ CellLabel->"In[9]:=",ExpressionUUID->"cd919aeb-3101-4893-8360-ca90a34b5f8a"],

 Cell[BoxData["1"], "Output",
 CellChangeTimes->{
@@ -288,8 +291,8 @@ Cell[BoxData["1"], "Output",
   3.868592286450618*^9, 3.8685929550022717`*^9, {3.868594532156452*^9, 
   3.868594555136402*^9}, 3.868595763986228*^9, 3.868598253224373*^9, 
   3.86859952158809*^9, 3.868599567371044*^9, 3.8695655289855413`*^9, 
-   3.869628033543208*^9},
- CellLabel->"Out[53]=",ExpressionUUID->"b0e60b84-1f1f-4268-863d-06c8900b3b0b"]
+   3.869628033543208*^9, 3.8698124352622843`*^9},
+ CellLabel->"Out[9]=",ExpressionUUID->"86a70908-b758-4cb1-8f2d-c9ea5b1964b6"]
 }, Open  ]],

 Cell[CellGroupData[{
@@ -769,7 +772,7 @@ Cell[BoxData[
   3.868592681493209*^9}, {3.8685953160084543`*^9, 3.8685953222658443`*^9}, {
   3.8685994009874973`*^9, 3.8685994014156713`*^9}, {3.868612439195859*^9, 
   3.8686124712931223`*^9}, {3.868612700424255*^9, 3.868612713038343*^9}},
- CellLabel->"In[57]:=",ExpressionUUID->"589fa93e-41d0-46c0-b8d3-ab8df683f53d"],
+ CellLabel->"In[10]:=",ExpressionUUID->"589fa93e-41d0-46c0-b8d3-ab8df683f53d"],

 Cell[BoxData[
 TemplateBox[{
@@ -801,10 +804,10 @@ Cell[BoxData[
  "ConditionalExpression"]], "Print",
 CellChangeTimes->{3.868604232135943*^9, 3.868612603799876*^9, 
  3.868612825817466*^9, 3.869565741857703*^9, 3.869626276514852*^9, 
-  3.869628183693252*^9},
+  3.869628183693252*^9, 3.869812555675345*^9},
 CellLabel->
  "During evaluation of \
-In[57]:=",ExpressionUUID->"2946d86e-f35b-4b75-b95e-a5caa077a51c"]
+In[10]:=",ExpressionUUID->"337280b6-e158-4dad-b3cf-8ba39c1d5cf3"]
 }, Open  ]],

 Cell[BoxData[
@@ -885,14 +888,15 @@ Cell[BoxData[
   3.868595532269759*^9, 3.8685955362275753`*^9}, {3.8686126439376163`*^9, 
   3.8686126745347157`*^9}, 3.868612839045279*^9, {3.869628388967555*^9, 
   3.8696283893478117`*^9}},
- CellLabel->"In[65]:=",ExpressionUUID->"c95e48ba-da31-4186-9fab-e23d54171786"],
+ CellLabel->"In[12]:=",ExpressionUUID->"c95e48ba-da31-4186-9fab-e23d54171786"],

 Cell[BoxData["1"], "Output",
 CellChangeTimes->{3.868595287863023*^9, 3.86859563951444*^9, 
  3.868596926165069*^9, 3.8685971501862993`*^9, 3.868612678323316*^9, 
  3.8686128422832127`*^9, 3.869565750199133*^9, 3.869626296911549*^9, 
-  3.869628191962613*^9, 3.869628395777919*^9, 3.869628872721054*^9},
- CellLabel->"Out[65]=",ExpressionUUID->"605d794b-f8c9-4677-b077-aada9b92bbd5"]
+  3.869628191962613*^9, 3.869628395777919*^9, 3.869628872721054*^9, 
+  3.869812556427718*^9},
+ CellLabel->"Out[12]=",ExpressionUUID->"6beae27f-a003-413e-be5f-41b65ba6a545"]
 }, Open  ]],

 Cell[CellGroupData[{
@@ -924,6 +928,7 @@ Cell[BoxData[
   RowBox[{
    RowBox[{
     RowBox[{
+      SuperscriptBox[
       RowBox[{"(", 
        FractionBox[
         RowBox[{"2", " ", "p", " ", 
@@ -944,7 +949,7 @@ Cell[BoxData[
            FractionBox["1", "k"], "+", 
            FractionBox["1", "q"], "-", 
            FractionBox["4", 
-            RowBox[{"k", "+", "q"}]]}]]}]], ")"}], "/.", 
+             RowBox[{"k", "+", "q"}]]}]]}]], ")"}], "2"], "/.", 
      RowBox[{"k", "->", " ", "kProton"}]}], " ", "/.", 
     RowBox[{"q", "->", "228.465"}]}], ",", 
    RowBox[{"{", 
@@ -952,97 +957,108 @@ Cell[BoxData[
 CellChangeTimes->{{3.868604821326149*^9, 3.86860483549179*^9}, {
   3.8686048918037643`*^9, 3.8686049478840723`*^9}, {3.868610320523445*^9, 
   3.8686103884198303`*^9}, 3.86861285794902*^9, {3.869628970187204*^9, 
-   3.869628984962192*^9}},
- CellLabel->"In[67]:=",ExpressionUUID->"af705252-b6bb-4c46-94eb-0f4dbbf1db66"],
+   3.869628984962192*^9}, {3.8698130811522417`*^9, 3.869813081599278*^9}},
+ CellLabel->"In[20]:=",ExpressionUUID->"af705252-b6bb-4c46-94eb-0f4dbbf1db66"],

 Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    TagBox[
     {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[CompressedData["
-1:eJwV1Xk4Vev3AHCJQpKGqxRdKkPpIqJjyFrOKKTomlOEyhROKI0yxDFmiBQy
-VErm2T6mJIX4RoNriJQhEpIMZfjt3x/72c/n2Wvv533Xu9ba0qfcTRx5eXh4
-xsjr/+9FZgu/vud06px6a7JrLvKkDqfhRIdxZBnsvKVe1DelDbe6/C9vj3wJ
-XYmZd0KmDCFuIK6FL7ID0ssCxQWnjkP8SPLMbMQwfK0czfKddIXNrjvCfkTM
-QiYLsr1Gr4B5joPpeMRqrFtKkYprD4fZliBiLGIzCp675RjvlQQOIERZP70Z
-OQvdJ/O+JUF7fmSJmtUW/P12RCnmVDLkxibkX5MVxzexAdKFR1LAwfLJI9Ga
-rdgVWZmWvjEV2gdex+yflET2mSajm4Np0LHIUdPZvQNtZz+I83s8hKJNbWtv
-hsmjq0nGGVGRbKj03DagxpVH5R8yOf9pZ0N9qyMxOCqPo3MszYsu2fAh5Pdp
-hv5uXNR9H6D9KhvmFnbW8gnuQRmLcfv313NAZ9CHHXBTAcOuZS/LDudCU4lE
-h98NRcw9Zxyf35QP7RvO5CgXKCI1qPP16b586HYvCPj0SRHnzVc8XjWdD2O7
-mcqISljfnSe9UrIA1qW4c5aXlPB1tr7QuFsB7KH1S0K1Mg4cNR0PWFMIHgr7
-31RpqeDqjt5rN9SKoPcAa1nymApGV+d+XckqgsN0K8VrzipI8UgX97Aogj02
-fuE6iSqYPJnUzH+5CAYjX+vVzKig2OAPc9maIrCacnhWW6CKUn6FW7RpxUCr
-uF34XE4N196JpTVolkDhi8f9O0ENS7W8/w04VAJS7VzRQDM1/Li+1X2PRQks
-jfafowepYcS5fT3q3iVQLqGo8OKzGp7/SlE5lFsCe/0aMhqS1XGN+r4US4lS
-EGPNxjVupOD5KJ7Htl9LYeqVwdyufyhYECAQ7/mrFFr0U639mBS8W87qd+It
-g0CjQzsO+FJwZnT+hahEGUyZ3st70EtBtuDiEQ2jMmhxgKbrTzSQZTylwptb
-Bo+HYv/prtNApf3FRVEV5Ptnv95S79HA3kTj6/wvykDTNdr8u4gm3idGD2R1
-k/HnvwxYeWkiNeftn87V5RDoH7ysBlrocnTmQb91Odjy9ZyKttBC286m/D+O
-5aB1U7lhzFMLHcOuiC+7l8MUpzM844EW7ovj0a8LIOOjFbauF9JG/2398tpZ
-5aCd+r/9Y++0kSWcqi7ysxx+Vm12TnfWwQ5z1+BCrwqYM1ZONPTXwRsJap9T
-LlfAwqDeq5lEHew14k319a8APpFLsgaNOqjr9u7m76gKEDvR8/mnHGDAywPt
-ClkVQFlKtWYMAuZgxh+r7grQjqkInVgA7PfONFXrrwCpjSM9hkKIp41Uti8P
-VYBM5grhWWlErSMndOymKkDxfyrOBsaIB6c14uwECIC/b8tO5yG+eJ/uaLeP
-AKqTtE4jF1H5RvVtTXUC6EU5pikvEXdfo9D4tAjQYzUE6vUhZs+vtXKiE3DU
-ffZzkogu/lJxN75iRsCpGsv7DDdd/EovUnl3gQAHwaHSrb66eFYsVbbqMgGn
-j7FbJwJ1ceS+/KOE6wS4DIcuJibpYsGwaph0MAFe6yqtx5t10fGwaV3obQKC
-Tm7fkrCHioHv/lx9lkeAjkNzo9MBKh6qq72ypoiAmbMXL2vTqRimtraSUUrA
-WXZ7b78NFRkwdySskgD9oJsP90ZT8VShzF82rwjgDVU1X06mYmPoeVuFZgK4
-kZ8E2rOo6Fz+PGa0hYC9dzRdL9RT0a51cEnrLQEiT8dV6mapWJwgvf+fjwQ0
-5N0biOOjYYlx5UxgHwHXi/Xiz6ynYYPVHnprPwGTVenzwgo03Lsv9I3aEAFv
-35jVmp+goZzw8BercQLC369k73GhYen04j7dSTK/Xfk7Fy/QsDqn+Kn4FAGl
-X4SC06Np6CTu6Zn2i4A7M9WHv9fTcCRI011kgcz3H5flmjYapvfX0kIWCRDk
-ES+I6aXhe3O7/uklAi4Jnt9EmaOhcmT6i0cruKAqItUgxE/H2xvnVMZ4uTC2
-oeXCx/V0vMPXVrqLjws2EnLd/gp0TP/5zsh9FRfEpN+Fm1Lo+ISjw7mxmgut
-Mjd05Bl0PBb8WDZYgAuo1J3WcoKO3vf//eAqxIU51ZBjqS50/DZgLXJ4DRcK
-KGr85y/SUaLi+7SUMBd2UKPObomhoxtP9JYHa7nQxdTe+i2FjtK0b5eOiXAh
-1mCkueopHQmO3YcZ0oZH46/eKqdjdbHHwch1XOAzpSnZv6CjSkFQnbgoF6os
-Jz+ptZPPHRcvJJL2OZEcI9BHR9nUXR5r13NB0V6f3v2NjlErKVXepIfPzP7K
-maPjjll7z3bSqa4PMv34Gehwuz5l5wYuWHgaWx7bwMCpe47mTqTX+ywJyf7N
-wDFrak4G6aZLTyvnFBjYlMUobCftf93iXDOFgXVvD12cIa0VyC+VwmBgEWfL
-+nUbuTAdUtjmYcLAxE7nSAnS2REnA2gnGegiObLib9JWN1aad7kwMOPOfJUY
-aQHvzD2eFxmo+u/n2ZWkS88aLK0OYqBAwu7aIfL79scn2lKiGfhyjiFXQ1r0
-aOxDtRQy3iLlYATpatoB39dZDHzdxtloTNr1QLehfRm5n9IjFUKktypcl/r9
-nIG2QXJAkPt/tX3n9K03DPy7j1JykrTPhpcv5T4y8Ln/PYU/ZD53rXK5Vz3C
-QC26eWE46bZ5EXfTGQa28ztYbiJ97XshdYyXiU8LnijHkOezt99MLGAdE9NV
-efRWkQ55lVKVL8/EgCMT3m3keatXUqNZaky8lDD8VJ70QN6QQ68uE7+szrf1
-IesDEhSFha2ZqN1VE/iLrKex0Pa+9DNMrFIa2iBH+u41nyINLyY2Rty1OSLI
-hRnHaqszEUwcfV7hGEjW5wPLU4pLiUx0dvTXjSbr1+TwKt7bj5i4nBwndpuf
-Czn7jZ7U1TCR7+L9rzdWcsFhZe+s5BQT9bNXt/wg+2XD7I3m4iUm5uxmbSon
-+6l2VOa+wRoWdoXoNHuT/bat3Y3pu4uFg6ct6jrmCWhPXYx7Z8bCsPqGMwem
-CfCLSzvrYs/C6931NQTZv4ohDG1eDxa263ycV/tBAMc94osSh4XZB6wNN5H9
-jzqSKmFcFq5emB/3HCZgfN8z/h2vWKgR8aSnZJCAJBnHzvJ3LKxNa9k69YWA
-WeEcv6HvLPS74l597BMBud3a/9OV0sMUFz/xqv/I9fked50L1MNKk85kVXK+
-9Xp/flh6Sw+vdFFz5RsISGWf7fNK0sN67t3nm+sJ2OXqZfKjSA/Ds0tlBmrI
-eXcynDL6WQ8FJc/3q5eR85VRydeDh3A4ZW/f6wcE2K2XSK5ZOITfRUU/LJLz
-e5dI+oerAvooGHfLjOlLwJCQvKj2JtLfFoVCfQhw5tsfUK6gjxyzMJWVngSw
-Zw2cCqz00epwwLkqRwICP15WyyjXR90Ts6ssjQjIfNLdctPLAH1SBqcdJQiQ
-svNPo94wwNM/nmf0biHn35bd3ksRBtiXvvaTyV8EhAb7SPhkGuBbFldWUYQA
-V8f1Lo5dBtjYe24qf7kClKVZAnQ0RCGH5JsT5P+QuFNIXSF8GMvmPdZGZlRA
-azCn7HKGEcq6uu3YKFkBoe97hGVaj+J/H9ueXRUshzcv2q+I0kywOGX/TdZo
-KQwIf1LO6ziGZjbJRX3ZJRDso5imaWyKPFdoW32Ci2FWPc7Ud9QMtwwNTBww
-KYI1Hzr0tK9YYHlsVp2AUiHkJhyf+zZliVFOcqPMvQWwPX/umxHHGisuzVAv
-7MgHyW/v76Xw2WB+U3jIskcujF5pikq7fwI/dlqp2cRlQ46PkV2cuC1ejP/j
-q96bBf/gfxat723x0rOXLZMnnsDGErep/jQ7pGz4rF/fmQm63J2jLrancDdH
-07ol6BE8y9dbQdlmjy+LbCd3HX4Imca+xMdBe8wfWzGRuu4BvA9oaptPd8C1
-3xeHTv5OB2c7ubAJd0dM+XaaoPangbHiT3s5jdMoUc1fGNmXCtfD++6+nz6N
-P43ndF4v3IdN8R3nbWrP4KkPgsWWqikwqTyZZX3tLDbPbkmwy0iC7dwLK5Io
-ThjWbOT+XeUeFE0K/S7hdUZPXe/kVm4iRJVVH2XWOGNkFNq5nb8D7buvHcoP
-c0EdgcsLdjIJcPyfE6m39F3RzMom6SlPPFCkcsSK+N1Q2U+07epUHFyuXq2U
-98YNNb3W1831xgL/uP9D6chz+Dt9/z3V7zEw80qw2s3aHTdneNe3zURDUuPy
-76sbPNAmvr718cZosJsVYvp5eKDTM+d3qw1vQefvbafcP3hgi/EqYgdPFPxi
-HsxJVPNEG0sxysHECCCExAQGYzxxuCtsW4ZhOKg8mFaVnPNE+rGPQVOzoWB0
-80vkLTM2Zt21aKxv5gBLP9JQ3YKN9rcD7+Y0cABFNIR6LNm4rqorK+EZB/bF
-RwbJ2bBRj34v2L2MA5seaVyptmfjxYNhcooZHOh+EeU05sHGh1/PFw5d4sBZ
-fm26Xhgbr36JCqQqcMCucWjFeDgbW57ER5jIcsAqIromNpKNmbezmuylOWD4
-17BmXzQbg3yT+zibOaAkE7PP5w4bV2Qqugys5MAM/ev2Bw/Z+EXJuKq9JwQm
-BGJ7DmWyMa84svdXRwh8fX3w7sRjNuZs8FLc+jYEOv+N/Usrm42BxZ8sTjeG
-QKWDjnB7IRsjuydahUtDoER+pPFCMRvfEuLHNQpCIHcsNliylI1EN/WvM9kh
-kOo1wutUwcZ7kSjUkB4CiZS4WhEuG4sorozZ5BCIWdC5VlzJxtB/GU93J4ZA
-aO2IllU1GwXyL2jZxIVAQGDc/HING4+p5v2KjgqB/wPx8vnf
+1:eJwVlHc81f8Xx0uUlZC0EK4V2aNrdM/BvddMRmaUlazCFQ0y7/0ohUhWSCpK
+SUSRkSKzMhoSkhJFRgj1bfzu76/34/l4vN6P9/uc1+scKc8g24McK1as4Fu5
+YsX/z3sOv39MlfZTevjkT0dv6zI407K/zyb5AZA/qffpqSVSzr+Li5BIboV8
+0xDIUy2hTCvc4kos6QG5dWuqRcLqKemj6c85k/uA1GWZcV+lg6JXEHUxOmQI
+dvtV8goefUWZVa5UOVUyAgcOLRZqhA5RmvhupunljkKeXl2tMe8oJeNr3uJS
+0ji0k6bWjipNUPxa01yqoidAi7f9tULKLMXgekIDI2QKhnPHzcqCFyinDSlC
+WeqzENTTUR7wcJkScpfcFF7yHTpZ3LcNA/9QXLZphtmT5kEhx16Ex3glUJOV
+5bVyF6Dk0d6qoVWcsDFQ+uz3pCUwhaCfK6V4YcWA2K7u1T9haI+5zvZTa2HC
+bOPMnehfQBHM4wonBKFOnt8uMOQPnC/cmcbwEoGijNVcFhN/oTT7tu7R46LQ
+FhczkK6+AkWVV9UIFm+C8u8nk0JLVuLGbmPxD9biILf9+NhHHw68PLbFbqFz
+G+S4h4EtaRVu6hPg0FSTAmZX0He1XE6s1u+oSp4ggWOpt/100hpkCs/Of/ux
+HZ6NetxxM+fGgrXFa6YzlQDFDqx5vpoHdwtuBhu6Miieda65Fc2LYZwNC9ze
+avDb10rcL2Qt1mYb/t6iqgUhBRbhb5UFUCZ25zoOdW0Y6zPtMpkQwOFfNpIT
+ZB3wvXoFL6gLouhUrl+hCxlshHPnQkqE8IbH++GZUQPgv7Pl8k1uYfT6sENN
+TJQCLWbZFiM+whh1+qK64G4AvdiMa9ak9Tihnr9q0xiCzOx5B9VcEXwSJxis
+k2AMS89ZD78lbcSt34cXy3abgTfwkoUWNuIFXtsmPZI59N5NrtJ22YRj4v1Z
+b/8zhzsXMu9GyW3G9ib9k1r3LcHb+WaR4KMtuGGPiqubvDX0jj5L05oVx/Ex
+7xU5nHuh788Zbcp2adS+BbfEa1whWcfyqLWrNAo7fdkcZOUG1CCBe54p0ph7
+lRLfN+oGFR/S1BJ+SOPl3leeQyIHILnpklL3YxLy3d95rajRHWhnSqU8XWRx
+5JMHbqZ6wT2RnrXEWQXc9KXxU2GWL9SFbB3VrlVAZj99ZHrCF5pfHHz4eUIB
+g/xb7SJ2+cGb0798aObbsTIpc8Dqkx8s/yY1cvIoorLq8eBy9QCgfA5nxBNK
+J6S/2667zB0VIn1xcSqIJ/D+mHfyyEQrKTVXa+vgdWvbCJa1h2D9ztN/onb
+aeDUGpuIUd1jsJvqohLlr4FeO9YKzXsdA0W3mHOUbA1s8yi/PfzgGHxOfmb6
+aFED54TvVI7tPw4uc96PG8s18Wr39FO1myfAuOZiRZO8Nr6U83GSVIqEiqc3
+Rkigjbu8LGPMbSNBsrdWkOmgjRWOtWT7E5Hwd2LkCJWljaaqwlW8rZFQLaai
+9PQj+36aQt6E+ynYEdNytSVPBwcWi+4oJEeBqMlSevt6MvY23vAQ84mBuTaL
+ZRllMnbRFJreHI+B5+YF+2LoZJzKCPE8fjYGmFZm0jtPkHE49ldN8F223v5S
+2bX3bN7blqP+k633ho7om7pYS8h5vGbFwo2xC8oDT3TxtppPzlxGLDB9v5zX
+GdTFkoGquYXiWNALTHWcEtDDPP22izfa2frQT6MuR/XweuaXsTa+OGDGJfzT
+Bn0sesZh9i4xDtw5Bz1TnfSxc/xBEiUnDvQJtZZvIfroDVOXiZtxMHem/9zV
+a/r409tl7nErW5+qtEWI1wATtvwY/MQRDwYFXVrfXhngnpWvQjpD42G+fqN/
+oT8FU7a9ZK1AJizbqGVbxlGwNWu+h8uCCb8/m7YtZlNwxzwX95w9EzgFTspZ
+tFNQx9bfMDqACaL7Bz/OywM2PrT0P3yRCeS/BftonwGbBYiuvFEmGKTVJM78
+BqxZp2kpMs0EyfVfBy15EbkyGVsYS0yQLV7JvySF2K3m3/OJhwUqXRr+FjaI
+in09TwqUWQDbLsotlCFG7ck6nhjMAiM/KUp7LSJD5+TmiOMsoN4rtc9vRSyj
+LDTbxbDA1KSFaTqMuMEsOLQ+hQXWQUsfcwUMMdD98JPQUhZ4PnK+TDtsiB6s
+rYknRlngzTN2f8sJQ0wg6fd/n2CBjx3jxQzTELtSE3bYfWdBwHjin+xcQ5xc
+s1+25w8Ljq6r2zfdaYgfuR1s6zYQwDogsSlT0QidZfPqDxgSQPHubPfbaYQf
+1hiGDtIJWPQ9HmFANcIUpu9ZsCTAl9H7fsTNCA3O7jKrcSDAnEVc35FqhIKu
+Cy9O+RPAkajp+C/PCLtMrALUggioTf7A3VtihD6/S/U6QgnYkaUXeKzZCJ95
+9A3diyRA4Na0xpMlIzwW66DsmURAS9ml0XROYzzhbuZkk0pAdKVpxiEhY7T+
+kX5e4SIBs/WFP/mVjJHZ252VkkvAy26HRsf9xjh8w3GrcQkB516vYigGGONB
+z7KmsFICqO/ukv4cM8Y/lAP7Uu8ScP8Tb0JhqjFquV7ujL9PQNZiw+6pZmNU
+CRd/6/SYAOv/Av496jHGdTxrSW+aCOBZsbk87b0xWvDx0na1EHCSJ1SEvGyM
+9j128e0dBGgKSLbwclFRMdjp3fQzAr4JPz82JETFI2ffnvr7ggA3MfmBOCUq
+Xhc1lRzqJUBU6tU5ezIVo2ZtXG6/IuCFbCxFgUbFPwSvo88bAlB14Mrz/VS0
+npZbd6mfgGXN03YFAVTMjZrmEh0goJyszRV6nIpvOxV1IgcJkDZK8d2URsUA
+aeGHvMMEvKMbbJnMp2KBg62Z5gcCLlh87ay/RcVITbUzJiMEWFpnnDpfTUU/
+/phMk48EcNobq3o9pSJ9aFOa5icC6p1nP2j3UpGS/zKDb5SA8P15adzDVOzN
+TmzqYrOKlzl1YJKKyU836cR8JmD80NKP0mUqCos784uPEVAQeK04houGhf9k
+vYvY7BRi42wnTMPvKU77xMYJEAr/yyu3jYaEaNu6aDZ3nLxVt6xEQzMfs8sv
+2BwX7XSkk0zDbOM78nxfCNBncknm02h44XlzixabF05X9ATb0rDzq0OOBZtv
+Jx2INz5AQ8WK7Y172OwSu8rxXQANM+eLmgzZzB1WrBhynIZv2jj/SbH5vq/F
+3zUsGsZa7k+eZb/n5TrTk59Kw+RtbUmlbBa0vnBdO5+GodN7NjuzucF454ln
+JTTkXf7Pcp5dT+DOAUuvBzT0z+igRrJ5i1K05K8mGta0l4stsPvRJkFaON9N
+wx/N+bMubA4Xbm2VH6JhVGTY4F12/2RWB1xq+ErDFh3B3z/Y/e75KRBkv0jD
+8xIQoMjmqKkKo28cdIS+VpoV258dIw6i8evouCB8sNiD7d/ptvz6uwp07Luf
+fMSR7bdOnVGqiTYdtys78Om9J2C0bMz7vSEdz0YU5/EPEQCZKvz8++i4LrRX
+MO4dO5+JvcOFh+iYuBitIMvOV05U+D3do3Rk7N5/tbqPPe8HG1wOJdHRs7lG
+tZydz2vOnip/s+m4azHzmOhLAmx3r+a4WERHv1DflMAeAkq1rG4+eUTHpC3L
+OZPPCfBe9X5JfI6OqkFjwmT2vAgvxXZW/qWjvkZfp1wzAY0Tspct+EwwzZas
+tOYJAVt7D9NPyJig6DjFJK+egN6CP+mvHEzwKvP7a6lKAmLSr/gGeJngPIRf
+Kyhn5+80zYAj2ARHZX0mhcoIOBOU9En1jAlGjguaP7vJnh+KuMbZWhMM2fxm
+y9V8AqbVH3NJt5mgp1X8qepLBOTKHuyvfmWCdVJpso+zCFjiL40ZmzLBoICH
+JbfTCLgzYNBlKGmKw35F9BkW+38nXAOXmaaoxOLxIbP32/uwj9fvnzfFB6+f
+aTz3Yeed4Tt8NNcUwyg8hk5ebH8Dj9p+v2eKN4pkg2xc2f4dOEee+GiKo08C
+FHt3s/crrY5zEM1QquzBVXFVAjyExPIe/TbDk7s2TLpOskBGoPDNKW5zTFPX
+jo0xoIxXgVBAxFz9FE2POL4kQX+nFrx1UrmGHH/gbVWPwsYSxZ+5S7mGO7W
+aJHZwgLmUIT21Wpz7M+ziBIsYEHxzYHnxFELTDtXNEKxYoGkR9wVo1gLrBKS
+kjhlyoKsTdvD/iZZ4NJgxUilEQsSE8LFwostULIt3URoJwsCDwoFHHxngXxh
+DxKjt7FATcqEm4qWeO3NrROGM0x4mFVhtJJ/NwbxF4z1JDLhRcKZBxFXrdCh
+m2l9uioeEl8P8su+sMbLPD6hIBAH3U97IwWNbZH4qZHy2jcGRvk/qJX12eEf
+hvnajW6nICFc5YqejT263eaTkHY8CUs66fYnJhxwldh7x8m9x4DvTZ+pQaQT
+BhpF/nQrOQp3Ml2XJ+ec0Wb8NTn1CAMk7i5PWp3Zh3JxRs1kWjCIT76+lM/p
+hpLudTJHfgTCRGRHypXL+1GcM/e3LeEHpeFWHumb3fGFCCNnW4oPKONbpxev
+3bFB3n4mRNwb1lcdnhu54oGmJeclXrV4gGEtaSLA3RPFXHRv2VUegMd3TVeS
+t3ph5JPBiYV9rlBsc+Lh0GcvPFKamS0l7wyv4zt6fhZ6Y86Nbw3X4xzA30P+
+7EzQQSxImlQX+s8ObFTmveR1fXAmoojDptMGos8N57xe8EGVM8tdWT/3gEhG
+X6hb4yFsaInIyomxhFm12ZJ9Ub7Ym982v+OKGUjUHluZS/bDXwkV1YYSJnBv
+lvdXFYc//kpeK/y4nAopDxqs6Y/8ke9B7OHmO0bQuz3K7O7ZAHz+pJL4Uo/g
+qry/4Lx5IJa7qv/oM6QAWbJU9B7XYZzsPnn7YqA+RDSsUS3rPowVHgLrK0rJ
+wDUdd10q+Qi691cUBXvpwGIbT8PhfUH40jCV4qakBbnt/36dEg7Gv72Gw+91
+NcBjiZceExyMI7v1549nq0H/r62eQW+CkXM61dSYWwV+0HeVZmuH4OcfrNZc
+dyV4yCvK/TktBCWWdX69Wr0dNK4taIovh+C43g3TgIdyYEV8Sj7vwECa8jru
+sE4ZMDFPttRxYmArKawdm2QABXR5B50ZOJ9kqCRYKwPqGckseTcGptmsbq8p
+kQGRIt3IBi8GBv+mGJATZWDgaYrft2AGhi97R34wkwFfLgOq6VkGiu8L8u/t
+JIFH+9jK6XMM1C2W/CjQTAKXpNRHF5IZWFZ3r8mqjgSWG8b1hlMZuPbM1/q3
+t0mgKpumHp7FwJDwW+4bkkmwSP0ice06A23qBedVbUgww31h0KyYgRkdG2iF
+ZiT48mxXzswNBp6U5h7faESC/r0XNujfZqDG1F6RtZokqPOm8PdWMFCuV+KL
+xAYSVCl8bT9WyUBO+qakW2tJcOfbhQTx+wzc9XVvmP5qEhQc/crhV8Ouz8iW
+5+CSNGST0xsFahlIN54r/zsjDWm/KVGVdWw9d3rmpS/SkNj4Vd+lgYGw62a1
+/og0xDPTf/57xMCHrtfWD/dLw/8ABEuWvQ==
       "]]},
-     Annotation[#, "Charting`Private`Tag$584362#1"]& ]}, {}},
+     Annotation[#, "Charting`Private`Tag$94900#1"]& ]}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
@@ -1074,7 +1090,7 @@ aO2IllU1GwXyL2jZxIVAQGDc/HING4+p5v2KjgqB/wPx8vnf
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
-  PlotRange->{{0, 300}, {0., 0.10576051518102296`}},
+  PlotRange->{{0, 300}, {0., 0.011185286571355386`}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
@@ -1085,8 +1101,205 @@ aO2IllU1GwXyL2jZxIVAQGDc/HING4+p5v2KjgqB/wPx8vnf
 CellChangeTimes->{{3.8686048309051447`*^9, 3.868604837191936*^9}, {
   3.8686103236765223`*^9, 3.868610390025558*^9}, 3.8686128611831083`*^9, 
   3.869565753444449*^9, 3.8696263002800617`*^9, 3.869628196004326*^9, 
-   3.869628399732347*^9, 3.869628986502192*^9},
- CellLabel->"Out[67]=",ExpressionUUID->"e98c5906-814e-490d-b5ad-0e56295f99ab"]
+   3.869628399732347*^9, 3.869628986502192*^9, 3.869813083825735*^9},
+ CellLabel->"Out[20]=",ExpressionUUID->"e690fd3a-0037-42ec-aaf7-a005d8ba954a"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+ RowBox[{"Plot", "[", 
+  RowBox[{
+   RowBox[{"{", 
+    RowBox[{
+     RowBox[{
+      RowBox[{
+       RowBox[{"(", 
+        FractionBox[
+         RowBox[{"2", "  ", 
+          RowBox[{"(", 
+           RowBox[{
+            FractionBox["1", 
+             RowBox[{
+              SuperscriptBox["k", "2"], "+", 
+              SuperscriptBox["p", "2"]}]], "-", 
+            FractionBox["1", 
+             RowBox[{
+              SuperscriptBox["p", "2"], "+", 
+              SuperscriptBox["q", "2"]}]]}], ")"}]}], 
+         RowBox[{
+          SqrtBox["\[Pi]"], " ", 
+          SqrtBox[
+           RowBox[{
+            FractionBox["1", "k"], "+", 
+            FractionBox["1", "q"], "-", 
+            FractionBox["4", 
+             RowBox[{"k", "+", "q"}]]}]]}]], ")"}], "/.", 
+       RowBox[{"k", "->", " ", "kProton"}]}], " ", "/.", 
+      RowBox[{"q", "->", "228.465"}]}], ",", 
+     FractionBox["0.0047", 
+      RowBox[{
+       SuperscriptBox[
+        RowBox[{"(", 
+         RowBox[{"p", "/", "44"}], ")"}], "2"], "+", "1"}]]}], "}"}], ",", 
+   RowBox[{"{", 
+    RowBox[{"p", ",", "0", ",", "200"}], "}"}], ",", 
+   RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]], "Input",
+ CellChangeTimes->{{3.869817253296451*^9, 3.869817289793811*^9}, {
+   3.869817352406551*^9, 3.869817388220334*^9}, {3.869817422212582*^9, 
+   3.86981746543668*^9}, 3.869817584991996*^9, {3.869817615929224*^9, 
+   3.869817632940728*^9}},
+ CellLabel->"In[48]:=",ExpressionUUID->"c924a6e4-1400-4616-9786-41ee8ec594ab"],
+
+Cell[BoxData[
+ GraphicsBox[{{{}, {}, 
+    TagBox[
+     {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
+      1.], LineBox[CompressedData["
+1:eJwV1Xkw1YsXAHBki3oSIXlxrxRNz5LlRcs5yZJ937Iv1/1a435RSZJdyFYk
+S5QSXpQikkRlS5t2S0oko6zdbsnN7/v748yZz5yZM2fOnJlD8zloy+Dj4eFJ
+oeL/mZ7FvVi68fmeOoWC2tYqEm4Jq8iHEHWwU+WMQQvl9sO1FcpEG+zvvSN4
+k3Kezpjb6lNdcNFEhl1JuVM9nClIPIXlh01iBZRj1i03LZe9gB6dRscTlPNz
+DfcsZb6BF0Fqvf6U68Qy7v+MHoS3q0wDjSl3Z/TvZzNHQLRw+7+KlGcTV00c
+rv4Igf5yGotXSDh2v06FU/YJxMv3OvdRFuGzC44sGAd/00//FVLOxx+1C5kT
+oLvDRM2HsuLxwrnwxEnokOgYUaJ8rXWX1mz0FKQ4X7w7VklCj15C81fmDAQq
+SK+0p1yzSyJWu24W5GqiIvkpzw36nY+onoOAmrfi1y6TECsrPDpftgAfx1KH
+fl4i4UGzy4rtxd9hOm5ANJ+yiEvNpvACNrTYa4eoUi4osGTOZHJg82enbusK
+Et7rlKWqpv2EWmWX+wMXSdj0eq4qJPEXRF15Pe9J+Zrkma9T0b/hT9qKedcL
+JPTmDIZ/YS6DNac7SbaMhNSz+T15Gjy4JEbTjT1PwjYCDmjW8WDD4t9bhktJ
+YAnlHGVV86KfS2FjRgkJfwy17s6WrcDA2r/u1p8joXzdsGX2Rn6M3vzGa7aQ
+BIPxpBG1Yn58RUgZqVBOT3zDc7BAAEeVPIcyCkiQuR+971umEG6UtHLecJqE
+llzFlxmiwhiW+kpVM48ED58+v21pwsgwNbYxyiWhgndjclDiSgyKox33yCZB
+A9u7J6NF0d8+pcIxgwTTViHLz0wxPP8nxYWZQNWXfdhSs2KY8eXOadl4EjIU
+PJ8o7lqDjbLdft1xJET43I60f7AGS+0E+iRiqXk/hz1oeCmOJ7wC2QGHSRj9
+Nux9mC2Be4rd2tcFk/D6VHpDtq4khnCNPMhAan9quiurjkliZ7NV4yOChPrw
+09feCazDkVT/C8EMEk6wTZf1pKTwXQnfrYMeJChwbxVzdWRQkr/W4pgVCZKl
+jNl1R2VwRTB8TbEgQRgkDFTbZLCkRDAi04yEmbjQKQ/j9biXe8cszZiEu/xK
+evecZPG2HW+kAZDgvir3TfxhOfQ5pZb16B8SHq/WWcnvroAlMs93aQuR0G+2
+9uyTGAXMlF8ye8pPwtu06c2FxQr4tcrdgsFHwphApYHqkAJGbUjQTOKyYIkr
+HefsSsMzChETpxdYsHX614//XOgIui8f/XjPgtQnrWMOjptw5lHHjqF6Fuhn
+Gd6rtFLG84G599hWLIgzvVnVH6qMRrNskyZzFtwVUMzjZiqjY0tS7yETFujF
+8Pjb9SmjUO5g2zd9FmgGtIgu71fBav3U0QZtFmzap+HkqL8VvSYdDW/LskCQ
+IzfNr70Nc+zmL14dDQeG2uJqjdXq+L3BkS0eFA4WKow5KWMtnLB9qypiEQb1
+L49F9uzYgfLD2rwNbqGwc+3WOv7ju9CAt6YqYTQImL9zSy45I8b4jI5ExgbA
+wc4JzToLfZxVy+9scWXCyQOCStxgA1ybFnnHQI0Bvl1H/TUIIzxpJNrab+kL
+7ULrbYJc9uPe1KK/JUq9QUM8Ln8ywhQTk0SUP417gtfC7m2WweYY6SJd/l3S
+HcpLH8pxMi1xn2m50rz9AbhxSm+w39Ya09pXCFyfdIKI27utS8xt8NCqymm3
+UQdo6700n+Zqi+IW03/pfrGDAGPd1gBnO+Rjxqz8ArZQfEWjISjQHkVFVZVe
+XLCGnIXWE2GHHHBS0mE8PdQSPldHRhEsRxyT4mGMi5jDA6Pcs/HxTnhqctOi
+/U0T4LUaGM2Kc8Z2qazauBZj4NnbyFt2xAWj7i8l3fllCKnSHIv+lAPYNMD8
+M7PHAMIO+HZwElzx5AJn+GyXPnw8oeUHZ9wwU1yEqHDYC5J5hQNp5e744916
+5YeiCJfXbD3SVeSBWkdYbcfGd0M7S95wS40nMvd5e6o07YRShTK53fZeSJx2
+UhUZ0oXWmEUe4pkXlgaa+59X3QE0M6ElRSdv1DNm50XF68DNpyrp715647nf
+P3s2SmhDRucfdQV3H6w+V212tkIT8sjJQK8hHxznOL10C9hO/YONYnQXX/zw
+vKvH0UwDqlQEeiM++GKFemRrzjZ1EPzOu77dww8X22OKxQNVoVL9emzphB/O
+04+2JxzZBio+1ZbcEAYOv+k0HW3cCoMmzfUOUwx83JnDKVFXgf6L0f8A6Y+a
+Nkfkj1dugZ83bHbR5/0x3Va329t2Myw3sRqSQ5nUPbQLm8oqgZm6kWEMm4n7
+3KM+171XhC53tYAgksALkwZTUn10WOtLtJREEjgW/mpNVBcdPIjy1c8OEXjW
+lZ3/uoMOP1gS9VoxBCqNszuKmumwKZWzyE0kMK0lOUOnkg5x9W0Z2fkEWj2f
+i74aTwc9IetrDc0EhrTaHSrRo0PyqjS+Ly0ERm+JuL9Fhw7PxTvsZe8SmKQp
+5X1Dgw6EnNav2A4C26VFsp4o06FAQ0Z//yMCs69OXJGUpsN31w8vBoYIfNfX
+HTqyQAP0Xr959QiBL5KDimNnaJDub3sYPhJYftNaUn6KBvTwB3IV41Q/jsxH
+31Ea2CRfYYRMEzjEVh5dekaD4vSPt87PEjg+znC43EeDL9myIv3zBOYKOonZ
+dNMgriijVodD4HuZXIfqNhr0lT3kCfhFoEDR4gfnFhpIX/5jW/SbwMZQ++vC
+t2jgU/PvpcdcAj+Yq3U219Og9loYZ3mZQM8NG+SDa2nwPzht2lg=
+       "]]},
+     Annotation[#, "Charting`Private`Tag$128858#1"]& ], 
+    TagBox[
+     {RGBColor[0.880722, 0.611041, 0.142051], AbsoluteThickness[1.6], Opacity[
+      1.], LineBox[CompressedData["
+1:eJwV1Wk0lQsXB/AjY8lrTLRcyXx0TBGR2juZpzJkyJDZUxLneZxLQkQiQ6KO
+vCiSFDdKTiWZSaUrmRropiJ5Sx0piuSe98Nee/3W/8t/7S97Q2CkS8gKGo12
+gjf/38qnlsrPKz7drt525ZITUnBbhL4+gqgFY+MxPQee22JrLmkSLaCpyBmy
+4DnfaNxHLKcb6PWb8rfwfF+PGSZEPIFJTvh+DZ7j1yzfWS4dACe/OE9Jntl5
+ltt/ZT+D9aFJ/nNAQa14VsePuBGI4q9MHub5QVa/zfew1yAzpXjvBs/c1NWT
+sVVv4KWlkkQ6zwkdtfT50ndwp+N/CV48r1rhepBVMAGjaV2C6jyzca5mNnsS
+zjGWK6e3U6BytHCGmToFRo+mg6/zfL3JzJAb9xHqP342O8TzQ9OUhk9hX8DY
+oBNGtlFQbSaduLmWC7mc2PCTPM+MBF+IrpoBaSat3pDnxHUib7+WzsIx+eKS
+ODMKOhu8+DcVf4PF4R6Q5XmVV7Uqs+A7pBwQ4ru2lYKCAqewL9nzYK2p9q7H
+lIJ/jErTdTJ+QIurv6Abz6rDM1cjUn9CmYub7TMTXj+Zs58+xi0CORBi8mQL
+BY9OjzA/hC3DrR/tfalGFKSfYz/M16fhsKna2vebKWAQsNeglobszl/HzXkm
+hU8fIav48Gjgt7FPBhT8tjRs5pbyo9qZ+SvC+hSUrXnllKsogBHc9l3mehRY
+TBx/rVssgMJJuRqxuhRkpj6jRRYIIk7QIoe0KZDriNs5nS2MFV8PJLlrUdCY
+pzKYJSqCA3x2/QfoFPgFPg5mZIigKDtpX5wmBZf4FNPCU1eix8FL+pnqFOhj
+24OpOFG0t1H4h1KmwK5J2Ol9mDg6nl3RECjPy5cDv8tyxTHSIkVGUo6CLKV9
+vSpmEtiyYWa2QZaC6MC7LLdOCZSd4LNZlOb1fR/VyRmUxIsB3HSz/1DwdvpV
+QOx3aexz569u5KNgOCeTk2sig2yRGlKaxrufrsnKqwkyWLjI9yHkNwl1zDPX
+XwiuwWjTnZMLCyQkf7dbNpWVRck0h8yJWRKUlm4XLxnJ4fME1z9ax0mQOR/C
+XXNEDm+KC2t3viVBBKQtdFrk8Bk3Q6FjjIQvSYc++lnL427475UboyQ0C6iZ
+tnqsQ1n5p0POgyT4rs57dixWAQWYP8vL2kn4W8xopYCvEs5ROTGFRST020ud
+641Xwou3gwYYhSQ8z/isXlishApbWRvusUkYF6y00BlVwsleqYu9p0n4tbQ2
+ydN7A5rezmF3nCBB6/PPub+8lPGGESeYS5KQ3ts0vsddFTOSq9dHWpFgfsqy
+tXKXJq5uZAk/f8uEJLv6q/2HNHFGV0qE/ZoJzYIq+UvZmrhCSEDCeZQJpvG0
+UNfHmthdk23ZPMQEg/2Noss2dIyQE+KwHjBBdae+h7u5Fj7sT1/rc40JQvMK
+nwU2M7C13alSmcWEEN0FMX0xPeyIuSv751IUONJDZmStDbEqvdhvqi0S6gYT
+WA+3bMH3kawe26wI2CqlVStw1Ax1lyv4hJXCIWwxr6TCE5H+aaffxV4CIu9P
+GtQ6mmOylLpIT2sonNwrpLZ00AIPH7m5lHktGIK6j4TqE1Z4fG3C4M+JQGgT
+lncO97LB7gz6p1+OAaAvmcSeirbDU+923WqO2wf+s9sYTgcdUEciRjTvhQ+U
+ne9SmM92wv6zKQMX7PbCzRzTkX6X3bima3Ot7WpPiL67bXeJgzOGrMpdxdno
+Di2PKr5meLvgXPRxAWlPN9hvbdK039MVX4ilnTzY5wLFV/Q54QfcMNPC9pmN
+qzOcnm1KjorZg+FMqx1d/LvhfRXrT4J0R92mR4Ol/Y7QaZV37tgxD7Q/rlYg
+c9Ue+Ha9fHsqyRN9p1xP+w7ZAm3HLb7Sw15YLrbA32xgA+lr5x37T+xFxaNV
+s9EFVhC1N6h9PsUbtS4rT7/YYglvkg2D4awPxg6aNP9+tRNk8gtfZpT5okJz
+SFRwpTlcltA63F3khzUJHP6zJTugjVxvqVG9D2vrT8jcOINwXqlUYZubP7pY
+WuR4BAI0xS/QiD5/tI5R33GlaxtssBf+peIRgNecMo9Wa5hB/RN65ovBAHy8
+kFCxsd0Usu7/1lPyDUSP1jhxmrcJ5FNTB/xHA/GRqsbTQZUtvH+gKK7sFYQh
+UwGyEcLGcJUu+Ch6LAi1HDIdRD5sBqFvfPJtfsGo4GwtODVtCJV6NxLPTwbj
+gRHjJ2NcA6AHVjktRYQgO9OPpbzSAEZsG+r2fAzBpoZfu4q2b4L+8jhtoEJx
+tDxx+q8offhx09lM+Wso1r3KWL/9nh4s3yE5aYfCkFWl0X1mjR7Y61lZxn8P
+Q1Xlw2q1wzrQ7au7P5wi8IcG4eZdrQ1SQURjCYtA1Z5wQdvL2uBHlIn1xRC4
+UbR4ybhMG+ZI6TrDeAJD2/rb5Au0QTV9fmEplcD8jUPZn45pQ1JdS1Yum0DC
+smS4d682mArvvs5pIJBOJViniGpD2uqMFR8aCew4oidRKKQNTyXb3dY1E9jU
+My51nU8bCAXDn4ntBE4N3ep8M8+AAn05c5seAlcW9/W6jjPgm/fYwMtRAmc2
+vVEhmxiAAfLqYq8JVOGM+xTdYUBmqEssvCHwsF3ew66bDFBmdipcmiBwOdTb
+TLGKAc5pV0IiPhM4ybgX/4rNgOLMN7cvcAks/6ZmI5HHgA+561b1fyUwwP4P
+C4tsBiQVZdUYzRPo7kZ/XJvCgMelXbT9PwnUbI31mExkwNrLv12KFgkU1OmR
+VIpjQGC1ccXfSwT2+DQveLEYUHM9an55mUCfTVzJs1EM+Bd55sNZ
+       "]]},
+     Annotation[#, "Charting`Private`Tag$128858#2"]& ]}, {}},
+  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
+  Axes->{True, True},
+  AxesLabel->{None, None},
+  AxesOrigin->{0, 0.00013489789812678603`},
+  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]],
+  ImagePadding->All,
+  Method->{
+   "DefaultBoundaryStyle" -> Automatic, 
+    "DefaultGraphicsInteraction" -> {
+     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
+      "Effects" -> {
+       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
+        "Droplines" -> {
+         "freeformCursorMode" -> True, 
+          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "DefaultMeshStyle" -> 
+    AbsolutePointSize[6], "ScalingFunctions" -> None, 
+    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
+        (Identity[#]& )[
+         Part[#, 1]], 
+        (Identity[#]& )[
+         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
+        (Identity[#]& )[
+         Part[#, 1]], 
+        (Identity[#]& )[
+         Part[#, 2]]}& )}},
+  PlotRange->{All, All},
+  PlotRangeClipping->True,
+  PlotRangePadding->{{
+     Scaled[0.02], 
+     Scaled[0.02]}, {
+     Scaled[0.05], 
+     Scaled[0.05]}},
+  Ticks->{Automatic, Automatic}]], "Output",
+ CellChangeTimes->{{3.869817259717098*^9, 3.869817291231227*^9}, {
+  3.869817410754409*^9, 3.869817466754157*^9}, {3.869817591183922*^9, 
+  3.869817634450905*^9}},
+ CellLabel->"Out[48]=",ExpressionUUID->"dd57294b-38b9-4e1e-a78f-a790c7d844ec"]
 }, Open  ]],

 Cell[BoxData[
@@ -1129,7 +1342,7 @@ Cell[BoxData[
       RowBox[{"k", "+", "q"}]]}]]}]]}]], "Input",
 CellChangeTimes->{
  3.869629300173646*^9, {3.869629377151785*^9, 3.869629475904372*^9}},
- CellLabel->"In[72]:=",ExpressionUUID->"e24edc1d-a3ff-4c8b-bb21-62fa5fe16f06"],
+ CellLabel->"In[13]:=",ExpressionUUID->"e24edc1d-a3ff-4c8b-bb21-62fa5fe16f06"],

 Cell[BoxData[
 RowBox[{"Clear", "[", "k", "]"}]], "Input",
@@ -1145,8 +1358,6 @@ Cell[BoxData[
    RowBox[{"(", 
     FractionBox[
      RowBox[{"2", " ", 
-       RowBox[{"(", 
-        RowBox[{"px", "+", "py"}], ")"}], 
       RowBox[{"(", 
        RowBox[{
         FractionBox["1", 
@@ -1171,8 +1382,8 @@ Cell[BoxData[
          RowBox[{"k", "+", "q"}]]}]]}]], ")"}], "2"], "]"}]}]], "Input",
 CellChangeTimes->{{3.8696303068554907`*^9, 3.869630331365758*^9}, {
   3.869630391073823*^9, 3.869630393423407*^9}, {3.869630598960945*^9, 
-  3.869630604762298*^9}},
- CellLabel->"In[84]:=",ExpressionUUID->"f53dd42a-ae11-462a-bf9f-c07d88f298c5"],
+   3.869630604762298*^9}, 3.869814938352353*^9},
+ CellLabel->"In[21]:=",ExpressionUUID->"f53dd42a-ae11-462a-bf9f-c07d88f298c5"],

 Cell[CellGroupData[{

@@ -1181,14 +1392,11 @@ Cell[BoxData[
  RowBox[{"PsiP2", "[", 
   RowBox[{"k", ",", "q", ",", "px", ",", "py"}], "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.8696304109161053`*^9, 3.869630430062149*^9}},
- CellLabel->"In[85]:=",ExpressionUUID->"417870a7-4839-4f81-8586-d217e9ba44f9"],
+ CellLabel->"In[22]:=",ExpressionUUID->"417870a7-4839-4f81-8586-d217e9ba44f9"],

 Cell[BoxData[
 FractionBox[
-  RowBox[{"4", " ", "k", " ", 
-   SuperscriptBox[
-    RowBox[{"(", 
-     RowBox[{"px", "+", "py"}], ")"}], "2"], " ", "q", " ", 
+  RowBox[{"4", " ", "k", " ", "q", " ", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{"k", "+", "q"}], ")"}], "3"]}], 
@@ -1207,10 +1415,11 @@ Cell[BoxData[
      RowBox[{"2", " ", "px", " ", "py"}], "+", 
      SuperscriptBox["py", "2"], "+", 
      SuperscriptBox["q", "2"]}], ")"}], "2"]}]]], "Print",
- CellChangeTimes->{3.869630432023151*^9, 3.869630611331956*^9},
+ CellChangeTimes->{3.869630432023151*^9, 3.869630611331956*^9, 
+  3.8698149474004307`*^9},
 CellLabel->
  "During evaluation of \
-In[85]:=",ExpressionUUID->"bf7a1e53-2a8c-4713-ad3f-6433be47849b"]
+In[22]:=",ExpressionUUID->"5de6100a-eac1-4d61-85e4-48f2515ee2f3"]
 }, Open  ]],

 Cell[BoxData[
@@ -1242,10 +1451,7 @@ Cell[BoxData[
    RowBox[{
     RowBox[{
      FractionBox[
-       RowBox[{"4", " ", "k", " ", 
-        SuperscriptBox[
-         RowBox[{"(", 
-          RowBox[{"px", "+", "py"}], ")"}], "2"], " ", "q", " ", 
+       RowBox[{"4", " ", "k", " ", "q", " ", 
        SuperscriptBox[
         RowBox[{"(", 
          RowBox[{"k", "+", "q"}], ")"}], "3"]}], 
@@ -1277,76 +1483,67 @@ Cell[BoxData[
   3.86963009471539*^9}, {3.869630131374378*^9, 3.869630131770316*^9}, 
   3.86963019598332*^9, 3.869630440713748*^9, {3.8696304731362457`*^9, 
   3.8696304770559893`*^9}, 3.8696306183960943`*^9, {3.869630881142564*^9, 
-   3.8696308850055428`*^9}},
- CellLabel->"In[91]:=",ExpressionUUID->"a4ef5eb4-7778-43b2-8376-e6e4431de66f"],
+   3.8696308850055428`*^9}, {3.86981492794865*^9, 3.869814957021761*^9}},
+ CellLabel->"In[23]:=",ExpressionUUID->"a4ef5eb4-7778-43b2-8376-e6e4431de66f"],

 Cell[BoxData[
 TemplateBox[{
   FractionBox[
    RowBox[{
-      SuperscriptBox[
-       RowBox[{"(", 
-         RowBox[{"k", "+", "q"}], ")"}], "3"], " ", 
-      RowBox[{"(", 
      RowBox[{
        RowBox[{"-", 
-            RowBox[{"(", 
-              RowBox[{
-                RowBox[{"(", 
-                  RowBox[{"k", "-", "q"}], ")"}], " ", 
-                RowBox[{"(", 
-                  RowBox[{
-                    RowBox[{
-                    SuperscriptBox["k", "2"], " ", "\[Pi]", " ", "px"}], "+", 
-                    
+          SuperscriptBox["k", "5"]}], " ", "\[Pi]", " ", "px"}], "+", 
+      RowBox[{"2", " ", 
+        SuperscriptBox["k", "5"], " ", "q"}], "+", 
+      RowBox[{"5", " ", 
+        SuperscriptBox["k", "3"], " ", "\[Pi]", " ", "px", " ", 
+        SuperscriptBox["q", "2"]}], "-", 
+      RowBox[{"5", " ", 
+        SuperscriptBox["k", "2"], " ", "\[Pi]", " ", "px", " ", 
+        SuperscriptBox["q", "3"]}], "-", 
      RowBox[{"2", " ", "k", " ", 
+        SuperscriptBox["q", "5"]}], "+", 
+      RowBox[{"\[Pi]", " ", "px", " ", 
+        SuperscriptBox["q", "5"]}], "-", 
+      RowBox[{"2", " ", "px", " ", 
+        SuperscriptBox["q", "3"], " ", 
        RowBox[{"(", 
          RowBox[{
-                    RowBox[{"2", " ", "k"}], "-", 
-                    RowBox[{"\[Pi]", " ", "px"}]}], ")"}], " ", "q"}], "+", 
-                    RowBox[{
-                    RowBox[{"(", 
            RowBox[{
-                    RowBox[{"4", " ", "k"}], "+", 
-                    RowBox[{"\[Pi]", " ", "px"}]}], ")"}], " ", 
-                    SuperscriptBox["q", "2"]}]}], ")"}]}], ")"}]}], "-", 
-          RowBox[{"2", " ", "px", " ", "q", " ", 
-            RowBox[{"(", 
-              RowBox[{
-                RowBox[{"3", " ", 
+              RowBox[{"-", "5"}], " ", 
              SuperscriptBox["k", "2"]}], "+", 
            SuperscriptBox["q", "2"]}], ")"}], " ", 
        RowBox[{"ArcTan", "[", 
          FractionBox["px", "k"], "]"}]}], "+", 
-          RowBox[{"2", " ", "k", " ", "px", " ", 
+      RowBox[{"2", " ", 
+        SuperscriptBox["k", "3"], " ", "px", " ", 
        RowBox[{"(", 
          RowBox[{
-                SuperscriptBox["k", "2"], "+", 
-                RowBox[{"3", " ", 
+            SuperscriptBox["k", "2"], "-", 
+            RowBox[{"5", " ", 
              SuperscriptBox["q", "2"]}]}], ")"}], " ", 
        RowBox[{"ArcTan", "[", 
-              FractionBox["px", "q"], "]"}]}], "+", 
-          RowBox[{"2", " ", "k", " ", "q", " ", 
-            RowBox[{"(", 
-              RowBox[{
-                SuperscriptBox["k", "2"], "+", 
-                SuperscriptBox["q", "2"]}], ")"}], " ", 
-            RowBox[{"(", 
-              RowBox[{
+          FractionBox["px", "q"], "]"}]}], "-", 
+      RowBox[{"4", " ", 
+        SuperscriptBox["k", "3"], " ", 
+        SuperscriptBox["q", "3"], " ", 
        RowBox[{"Log", "[", 
          RowBox[{
            SuperscriptBox["k", "2"], "+", 
-                    SuperscriptBox["px", "2"]}], "]"}], "-", 
+            SuperscriptBox["px", "2"]}], "]"}]}], "+", 
+      RowBox[{"4", " ", 
+        SuperscriptBox["k", "3"], " ", 
+        SuperscriptBox["q", "3"], " ", 
        RowBox[{"Log", "[", 
          RowBox[{
            SuperscriptBox["px", "2"], "+", 
-                    SuperscriptBox["q", "2"]}], "]"}]}], ")"}]}]}], ")"}]}], 
-    RowBox[{"\[Pi]", " ", 
+            SuperscriptBox["q", "2"]}], "]"}]}]}], 
+    RowBox[{
+      SuperscriptBox["k", "2"], " ", "\[Pi]", " ", 
      SuperscriptBox[
       RowBox[{"(", 
-         RowBox[{
-           SuperscriptBox["k", "2"], "-", 
-           SuperscriptBox["q", "2"]}], ")"}], "3"]}]], 
+         RowBox[{"k", "-", "q"}], ")"}], "3"], " ", 
+      SuperscriptBox["q", "2"]}]], 
   InterpretationBox[
    DynamicModuleBox[{Typeset`open = False}, 
     TemplateBox[{"Expression", 
@@ -1448,8 +1645,8 @@ Cell[BoxData[
  3.869567239297608*^9, 3.869626935103006*^9, 3.869629942396744*^9, 
   3.869630078096586*^9, 3.869630189648635*^9, 3.869630294281685*^9, {
   3.869630480764154*^9, 3.869630505310848*^9}, 3.8696306349033403`*^9, 
-   3.869630903121903*^9},
- CellLabel->"Out[91]=",ExpressionUUID->"63a40424-d528-4ecd-bf57-728421a50da7"]
+   3.869630903121903*^9, 3.86981497908293*^9},
+ CellLabel->"Out[23]=",ExpressionUUID->"bd741071-db04-43f7-98ef-6e9e741872b6"]
 }, Open  ]],

 Cell[BoxData[
@@ -1459,9 +1656,8 @@ Cell[BoxData[
   RowBox[{"PsiPar", "[", 
    RowBox[{"kProton", ",", "228.465", ",", "px"}], "]"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.8696279016627913`*^9, 3.869627911089171*^9}, {
-   3.869630713019814*^9, 3.8696307684281597`*^9}, 
-   3.869630947432231*^9},ExpressionUUID->"d6114d76-496a-4526-9235-\
-2e2b26389c18"],
+   3.869630713019814*^9, 3.8696307684281597`*^9}, 3.869630947432231*^9},
+ CellLabel->"In[24]:=",ExpressionUUID->"d6114d76-496a-4526-9235-2e2b26389c18"],

 Cell[CellGroupData[{

@@ -1469,24 +1665,23 @@ Cell[BoxData[
 RowBox[{"Print", "[", 
  RowBox[{"paral", "[", "px", "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.8696307762609797`*^9, 3.8696307859832897`*^9}},
- CellLabel->"In[88]:=",ExpressionUUID->"ea14538f-b536-4fe7-a74f-f5e25d13093b"],
+ CellLabel->"In[25]:=",ExpressionUUID->"ea14538f-b536-4fe7-a74f-f5e25d13093b"],

 Cell[BoxData[
 TemplateBox[{
-   RowBox[{
-     RowBox[{"-", "109.13809733550005`"}], "-", 
-     RowBox[{"1.0000000000000004`", " ", "px"}], "-", 
-     RowBox[{"0.7563855522068201`", " ", "px", " ", 
+   RowBox[{"0.027166391551979432`", "\[VeryThinSpace]", "-", 
+     RowBox[{"0.0007851266342468511`", " ", "px"}], "+", 
+     RowBox[{"0.000023643212121192028`", " ", "px", " ", 
       RowBox[{"ArcTan", "[", 
         RowBox[{"0.004377038058345917`", " ", "px"}], "]"}]}], "+", 
-     RowBox[{"1.3930053245744016`", " ", "px", " ", 
+     RowBox[{"0.00047618392705276363`", " ", "px", " ", 
       RowBox[{"ArcTan", "[", 
-         RowBox[{"0.02187868155499205`", " ", "px"}], "]"}]}], "-", 
-     RowBox[{"59.119288258590586`", " ", 
+         RowBox[{"0.02187868155499205`", " ", "px"}], "]"}]}], "+", 
+     RowBox[{"0.002178093701515095`", " ", 
       RowBox[{"Log", "[", 
         RowBox[{"2089.0926720349003`", "\[VeryThinSpace]", "+", 
-           SuperscriptBox["px", "2"]}], "]"}]}], "+", 
-     RowBox[{"59.119288258590586`", " ", 
+           SuperscriptBox["px", "2"]}], "]"}]}], "-", 
+     RowBox[{"0.002178093701515095`", " ", 
       RowBox[{"Log", "[", 
         RowBox[{"52196.256225000005`", "\[VeryThinSpace]", "+", 
           SuperscriptBox["px", "2"]}], "]"}]}]}], 
@@ -1558,10 +1753,10 @@ Cell[BoxData[
      0], Re[$CellContext`px] > 0], SelectWithContents -> True, Selectable -> 
    False]},
  "ConditionalExpression"]], "Print",
- CellChangeTimes->{3.869630787819367*^9},
+ CellChangeTimes->{3.869630787819367*^9, 3.869815274817771*^9},
 CellLabel->
  "During evaluation of \
-In[88]:=",ExpressionUUID->"5d9c5e6a-5231-4a6e-a4f7-21f99a3e5607"]
+In[25]:=",ExpressionUUID->"ad5eddd5-5fad-4d49-aef0-c2ad504aea37"]
 }, Open  ]],

 Cell[BoxData[
@@ -1590,45 +1785,47 @@ Cell[BoxData[
  RowBox[{
   RowBox[{"paral", "[", "px", "]"}], ",", 
   RowBox[{"{", 
-    RowBox[{"px", ",", "0", ",", "200"}], "}"}]}], "]"}]], "Input",
- CellChangeTimes->{{3.8696309646911507`*^9, 3.869631022431979*^9}},
- CellLabel->"In[97]:=",ExpressionUUID->"37c239c3-4ac2-401f-963b-6c070eeb2f73"],
+    RowBox[{"px", ",", "0", ",", "100"}], "}"}]}], "]"}]], "Input",
+ CellChangeTimes->{{3.8696309646911507`*^9, 3.869631022431979*^9}, {
+  3.869815292808971*^9, 3.8698152980950613`*^9}, {3.8698153467792473`*^9, 
+  3.869815347142014*^9}, {3.8698165268242598`*^9, 3.869816537101454*^9}},
+ CellLabel->"In[35]:=",ExpressionUUID->"37c239c3-4ac2-401f-963b-6c070eeb2f73"],

 Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    TagBox[
     {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[CompressedData["
-1:eJwVzGk4FHgAx/GhwzEzj8hRUjRRyMpEQjQ/nhmGxrgmt3YYzB8ppKQsSqtx
-lMqTrRyxQqko6aAQKzSxIvU8S5tYSU9tpOxUmmZ3X3yfz7vv6ojdvlHKFArF
-87/+l1EgryxbNbB1l8NMz5ybALdVzQzjST1LPB/bq8oWoH1/3QVT0sZSNZns
-XwEBuqwSxYtJP6v27PgzZzsBevIHuXPiUVZSkdOrY6YCPHTIanonnmb1Omqr
-rlUTQHpyJHFKrGCFVtbW+0v94NGiwp8UayBTpKXTyvNDH91WbWGYEWJOFl95
-4egLlwLO/RovU9AKX53LgQ+iNnylM+lWsNxil6vB84anWdQHXTcbaFUzKjVP
-8NEw9NPeh3Z24F1vHFa7w8MWLfP6hRmO0JuzVPwq9YB4/lRpVSCQ0v/JPX+e
-i91dr63rPV3QkZ5X/rPQDbnBi03kO9nYcnC8abqZA1H3wWgmcYXa0MrcXC4b
-7SrLfeKCuFg75KOz/qMLmJqZRW+SPSA5euFwwxlnCD86WfB38sC53DpiEQdU
-lD0wkB3jY6LOUd1w31bcOO4wMujrDW4jU1LLc0Rys5N3Kc8H3n+3NSDeAW3S
-qtmcEF+MDZ+vH821Q4ybfUtMoB/CFOyhyje2KLnIvBkXK4D86bUFCtdNOPmx
-5VBCynZMXUoNED2yxmTt3n0kyR/WVVJ7p6yN6HQ9debw4QB0xoxWXbFnQslr
-eLwgMxChyhNG7QZWoDjfUipPDcKm/IvnQjZaQqIn8xw8GozUe/c3q5taICFY
-1CHLCkFJhn62bqQ5xg7ZRLJOh8KuYUnh61JTaBeeHc6pCINy2tQYW2kdqpeY
-p3YX70A1Z0LUXWyC9iRDzrrLP6Jxza7uamtjlBmVGzgJhKhwzv8SPs9AS9pX
-CnksRK+s6P3nP1dj9TaVb2sCwuHQyVZrajFCY79Z3h9D4eAmH6l9nG2I/K7v
-VkZhEThFsRwvPrYShXvexAqfR8D7YcHBpGcrkKazSoMRJEJR85ng9Wx9XDJb
-JE1+KcL9exJtxtVlWPxJaXn7jkiM6j5Z1uahhxqr6+llryNxlCobj5vQgVlE
-LV8eHwWLA0ufNJZoY8S9qWH72yjwe/qy3qYvxWDlgR9Ye6Lh+i6q1CJMC59v
-+DgyZqORtM3f1V6oCcWdpJvZu8SY4ejekqQswTYrV07anBimUkmwXEcD3WEb
-YuL2EHzQjGsL3E+HlojcLd1LkOEW2RWSSMcOUkF/nEIw6R/XKoyl45+kpQ02
-aQS2nE0ZCaF0GEtkX+VHCHI7f9tYAzoyG9ryTxT99+MNDSeo0eGg4n3tZhPB
-vg9SpcQSGrJpOcpTdwkMW+yVrp+mYUCzQ6DfShAgvvJi9jgNxMDmS3oHwVip
-OzPtEA2/MJe5cB8RtBYYm1RF0/Ap5OWT4ecEi2zUu7yYNCB8+Vr6KEFXk5fw
-qjkNedG++1ljBE6a4f00YxoYiZ0GF14RKGjlZgO6NPhkX4yKf08w8Lamdfc3
-Kkryxm6fnyGIDxZNj8xRMXVCX31wluAbP3zGfZqKzOL8OlsZwc6+v4Trx6no
-LX9AiflCkOXH7S8boUKv+rtv8TyBblA9VfspFRGXN1f1yQlyBDk6eb9TUXct
-QaZQEDTPur5b0EPFv0WmIp0=
+1:eJwVxXk01AkcAPChou01Vm2byBtHF2qZsfMk+c33S2QY5nD+yptyjZkxyTG5
+NnKMSKYoMu2K2FDJuoq2NOx4xaaU5O3b0KFTm1IopabZ3T8+72MVEecv0qdQ
+KH7/+X/rIu2pSpqGdb9uk2j6VClcXGhrESs5Cf2MGLuailLQpDbW2EiaoOLG
+vKhgVSn00BPEBpIuGH0Vz+w6WAp/KgfZ78W3oZZm6fNrQilcd1FcmhA/hCuC
+Oq9GKIW+IyMJ4+JJcGwqSb0/WgI+akPuc7EOVvF35alWlkA/1emb+UJjVFLa
+c6OVR8C9yPOP07wVaNGZ7zjlUAwihzkqg2qJ12vs69ljh8DPVvRuuddqnA0/
+3LBGWAitQxlJ151t0LLGoZpKL4DNS+2a5mduQHsFx46rnw/iz0crakk6mqqu
+ZX5etB/iel782OTniAsXV936a7sCDm43WKPdxcRK7oBqf0M2RPbujWZInJCb
+9veXONcs0BiaCmTbnHHAI/PD5FgGMJZklb3c44JVby6vb+XthbBpYgN3lyt+
+zIkKy6ClQXXlNfPZQyykmXu/9ppOhvOHXUYG/RGX6EtVOfOSYM9lgl/h64Z6
+7KQysVAOXX21UwWh7mgUzjVp6UgAqdcmtZTcgt7dD3oZZDycOMNok8V44LuJ
+uAfWU7FwZFqdHZ/iiW+zXxsEZ8rgeX1SsiRxK/LabUtkUilc3Xr0eE6OF64L
+rpsZJcWgxxt+XJTFxowD9JsBZSKguLXrVaV5Iz/q7iLQi4IDJrN+g/k+GMnj
+hhvHRUD89sjuWQUHUxfUXT5qEg5j2cwoOOaL5EjWjvz9O2FZyc/DBdV+6DAs
+ZdtZC6HO2C6tt5yLg8+Oud1ZHgqaRAvPded4uPSsollD3QaVllXmRCAfB3bn
+JNN4IaBOn6NIBvg4rsk749YRBFYcwy+rQgSYWtjsqlkXCBdu2xbeGxIgRSv/
+hdrvD8qer3RLoT+a/RZsdDNRACXylzFho/74tI7tcJ/Nh/Tvad9abwtAi4r4
+tmVJXDhru6Bvz6MA3CDIdXEP8QWDGT1TzY5AHLPMc5pR+MBpesu+yheBqK5/
+ZDbQxAbbiHquNjYIQ1IKDtcaesGI96XWoFdB+LB8HOaOe8LgqZ9+AHkwTgoF
+MpmNB3w8L3C1ngrGO1+MzEQv3UH3e2Jb3u4QHJnkRbv3uAGHvtUz/X0INuup
+k1uqEHqFDlKZnMTiJ54tNA+ApZGSjookEkXrrW8pNwPskFRTB1JInAozUGkd
+AT4kftfKTCcxusjG9oUVwOoDs3PaXBJnOO3MYR0Lslq7lMVlJDobbzKFKyxw
+MeQ3t10iUb2zpUHrxIK8xQX64x0k3svyze6yZ8GdJd2BZp0k9qxsfKpYywKJ
+OfPTvm4S7a2YgcuWs0DFWOHOvkFiFSdGHfCegJnQR3eHR0l84FA2xblAAIab
+rqU+JHHL11BvpwYCCqP9U2GMRF1av+WqGgKsE66a1zwj8Z+zztz5pQQI8s6I
+Yt+Q6Dfx/OBTOQEnCscunnxLoudQ563HMgLGi80WDU6RGJsSpHwSSUBWubLR
+aZbEgSFuwkQAATerrlGkn0jcyJsrm+EQYFL31b/8M4k+Rh8cdVsIiDi3sbZf
+S6LMTeyx2JWAxub4WZ2OxM5ked9KJgH/Ago4JY4=
       "]]},
-     Annotation[#, "Charting`Private`Tag$685471#1"]& ]}, {}},
+     Annotation[#, "Charting`Private`Tag$124026#1"]& ]}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
@@ -1660,7 +1857,7 @@ QaZQEDTPur5b0EPFv0WmIp0=
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
-  PlotRange->{{0, 200}, {0., 81.12437652072705}},
+  PlotRange->{{0, 100}, {0., 0.020156672528549236`}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
@@ -1673,8 +1870,9 @@ QaZQEDTPur5b0EPFv0WmIp0=
   3.8686736607060003`*^9}, 3.868675912288126*^9, {3.86867613698822*^9, 
   3.868676151839303*^9}, 3.8686773890716133`*^9, 3.869630809858492*^9, 
   3.8696308592115803`*^9, 3.8696309557226753`*^9, {3.8696309956820498`*^9, 
-   3.869631024177164*^9}},
- CellLabel->"Out[97]=",ExpressionUUID->"030448a6-0fbd-4cf2-a9d9-83eed4c20ddd"]
+   3.869631024177164*^9}, {3.869815285649044*^9, 3.8698152995861893`*^9}, 
+   3.8698153488166027`*^9, {3.869816532210001*^9, 3.869816539061536*^9}},
+ CellLabel->"Out[35]=",ExpressionUUID->"e052c0db-eb44-4443-8d7c-14a4d951c77f"]
 }, Open  ]],

 Cell[BoxData[
@@ -1702,10 +1900,7 @@ Cell[BoxData[
   RowBox[{"Integrate", "[", 
    RowBox[{
     FractionBox[
-      RowBox[{"4", " ", "k", " ", 
-       SuperscriptBox[
-        RowBox[{"(", 
-         RowBox[{"px", "+", "py"}], ")"}], "2"], " ", "q", " ", 
+      RowBox[{"4", " ", "k", " ", "q", " ", 
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{"k", "+", "q"}], ")"}], "3"]}], 
@@ -1731,78 +1926,112 @@ Cell[BoxData[
       RowBox[{"k", ">", "0"}], " ", "&&", " ", 
       RowBox[{"q", ">", "0"}]}]}]}], "]"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.8696310639392776`*^9, 3.8696310748684053`*^9}, 
-   3.86963113563931*^9},
- CellLabel->"In[99]:=",ExpressionUUID->"b86a339e-b949-4516-8eb5-8294678723fc"],
+   3.86963113563931*^9, 3.869815311381949*^9},
+ CellLabel->"In[29]:=",ExpressionUUID->"b86a339e-b949-4516-8eb5-8294678723fc"],

 Cell[BoxData[
 TemplateBox[{
-   FractionBox[
   RowBox[{
-      FractionBox[
+     FractionBox["1", 
      RowBox[{
+        SuperscriptBox["k", "2"], " ", "\[Pi]", " ", 
        RowBox[{"(", 
-           RowBox[{"k", "-", "q"}], ")"}], " ", 
+          RowBox[{
+            SuperscriptBox["k", "2"], "+", 
+            SuperscriptBox["py", "2"]}], ")"}], " ", 
+        SuperscriptBox["q", "2"], " ", 
+        SuperscriptBox[
         RowBox[{"(", 
           RowBox[{
+             SuperscriptBox["k", "2"], "-", 
+             SuperscriptBox["q", "2"]}], ")"}], "3"], " ", 
+        RowBox[{"(", 
+          RowBox[{
+            SuperscriptBox["py", "2"], "+", 
+            SuperscriptBox["q", "2"]}], ")"}]}]], 
+     RowBox[{
+       SuperscriptBox[
+        RowBox[{"(", 
+          RowBox[{"k", "+", "q"}], ")"}], "3"], " ", 
+       RowBox[{"(", 
+         RowBox[{
+           RowBox[{"5", " ", 
+             SuperscriptBox["k", "4"], " ", "\[Pi]", " ", 
+             SuperscriptBox["q", "3"], " ", 
+             RowBox[{"(", 
+               RowBox[{
+                 SuperscriptBox["py", "2"], "+", 
+                 SuperscriptBox["q", "2"]}], ")"}]}], "+", 
+           RowBox[{"2", " ", "k", " ", "py", " ", 
+             SuperscriptBox["q", "5"], " ", 
+             RowBox[{"(", 
+               RowBox[{
+                 SuperscriptBox["py", "2"], "+", 
+                 SuperscriptBox["q", "2"]}], ")"}]}], "-", 
           RowBox[{"\[Pi]", " ", 
             SuperscriptBox["py", "2"], " ", 
-               SuperscriptBox["q", "2"], " ", 
+             SuperscriptBox["q", "5"], " ", 
             RowBox[{"(", 
               RowBox[{
                 SuperscriptBox["py", "2"], "+", 
                 SuperscriptBox["q", "2"]}], ")"}]}], "-", 
-             RowBox[{"2", " ", 
-               SuperscriptBox["k", "3"], " ", "q", " ", 
+           RowBox[{
+             SuperscriptBox["k", "2"], " ", "\[Pi]", " ", 
+             SuperscriptBox["q", "3"], " ", 
             RowBox[{"(", 
               RowBox[{
                 RowBox[{
-                    RowBox[{"-", "py"}], " ", "q"}], "+", 
-                   RowBox[{"\[Pi]", " ", 
+                   RowBox[{"-", "5"}], " ", 
+                   SuperscriptBox["py", "2"]}], "+", 
+                 SuperscriptBox["q", "2"]}], ")"}], " ", 
             RowBox[{"(", 
               RowBox[{
                 SuperscriptBox["py", "2"], "+", 
-                    SuperscriptBox["q", "2"]}], ")"}]}]}], ")"}]}], "+", 
+                 SuperscriptBox["q", "2"]}], ")"}]}], "+", 
           RowBox[{
-               SuperscriptBox["k", "4"], " ", 
+             SuperscriptBox["k", "7"], " ", 
             RowBox[{"(", 
               RowBox[{
-                   RowBox[{"2", " ", "py", " ", "q"}], "+", 
+                 RowBox[{
+                   RowBox[{"-", "2"}], " ", "py", " ", "q"}], "+", 
                 RowBox[{"\[Pi]", " ", 
                   RowBox[{"(", 
                    RowBox[{
                    SuperscriptBox["py", "2"], "+", 
                    SuperscriptBox["q", "2"]}], ")"}]}]}], ")"}]}], "-", 
-             RowBox[{"2", " ", "k", " ", "py", " ", "q", " ", 
-               RowBox[{"(", 
           RowBox[{
-                   RowBox[{"\[Pi]", " ", "py", " ", 
+             SuperscriptBox["k", "3"], " ", "py", " ", 
+             SuperscriptBox["q", "2"], " ", 
             RowBox[{"(", 
               RowBox[{
-                    SuperscriptBox["py", "2"], "+", 
-                    SuperscriptBox["q", "2"]}], ")"}]}], "-", 
-                   RowBox[{"q", " ", 
+                 RowBox[{"2", " ", 
+                   SuperscriptBox["q", "3"]}], "+", 
+                 RowBox[{"5", " ", "\[Pi]", " ", "py", " ", 
                   RowBox[{"(", 
                    RowBox[{
-                    RowBox[{"2", " ", 
-                    SuperscriptBox["py", "2"]}], "+", 
+                    SuperscriptBox["py", "2"], "+", 
                    SuperscriptBox["q", "2"]}], ")"}]}]}], ")"}]}], "+", 
           RowBox[{
-               SuperscriptBox["k", "2"], " ", 
+             SuperscriptBox["k", "5"], " ", 
             RowBox[{"(", 
               RowBox[{
+                 RowBox[{
+                   RowBox[{"-", "2"}], " ", "py", " ", 
+                   RowBox[{"(", 
+                    RowBox[{"py", "-", "q"}], ")"}], " ", "q", " ", 
+                   RowBox[{"(", 
+                    RowBox[{"py", "+", "q"}], ")"}]}], "+", 
                 RowBox[{"\[Pi]", " ", 
-                    SuperscriptBox[
                   RowBox[{"(", 
                    RowBox[{
-                    SuperscriptBox["py", "2"], "+", 
-                    SuperscriptBox["q", "2"]}], ")"}], "2"]}], "+", 
-                   RowBox[{"2", " ", "py", " ", "q", " ", 
+                    SuperscriptBox["py", "2"], "-", 
+                    RowBox[{"5", " ", 
+                    SuperscriptBox["q", "2"]}]}], ")"}], " ", 
                   RowBox[{"(", 
                    RowBox[{
+                    SuperscriptBox["py", "2"], "+", 
+                    SuperscriptBox["q", "2"]}], ")"}]}]}], ")"}]}], "+", 
           RowBox[{"2", " ", 
-                    SuperscriptBox["py", "2"]}], "+", 
-                    SuperscriptBox["q", "2"]}], ")"}]}]}], ")"}]}]}], ")"}]}], 
-       RowBox[{
             RowBox[{"(", 
               RowBox[{
                 SuperscriptBox["k", "2"], "+", 
@@ -1810,27 +2039,28 @@ Cell[BoxData[
             RowBox[{"(", 
               RowBox[{
                 SuperscriptBox["py", "2"], "+", 
-             SuperscriptBox["q", "2"]}], ")"}]}]], "+", 
-      RowBox[{"2", " ", "q", " ", 
+                 SuperscriptBox["q", "2"]}], ")"}], " ", 
             RowBox[{"(", 
               RowBox[{
-            RowBox[{"3", " ", 
-              SuperscriptBox["k", "2"]}], "+", 
-            SuperscriptBox["q", "2"]}], ")"}], " ", 
+                 RowBox[{
+                   RowBox[{"(", 
+                    RowBox[{
+                    RowBox[{
+                    RowBox[{"-", "5"}], " ", 
+                    SuperscriptBox["k", "2"], " ", 
+                    SuperscriptBox["q", "3"]}], "+", 
+                    SuperscriptBox["q", "5"]}], ")"}], " ", 
                   RowBox[{"ArcTan", "[", 
                    FractionBox["py", "k"], "]"}]}], "-", 
-      RowBox[{"2", " ", "k", " ", 
+                 RowBox[{
+                   SuperscriptBox["k", "3"], " ", 
                   RowBox[{"(", 
                    RowBox[{
-            SuperscriptBox["k", "2"], "+", 
-            RowBox[{"3", " ", 
+                    SuperscriptBox["k", "2"], "-", 
+                    RowBox[{"5", " ", 
                    SuperscriptBox["q", "2"]}]}], ")"}], " ", 
                   RowBox[{"ArcTan", "[", 
-          FractionBox["py", "q"], "]"}]}]}], 
-    RowBox[{"\[Pi]", " ", 
-      SuperscriptBox[
-       RowBox[{"(", 
-         RowBox[{"k", "-", "q"}], ")"}], "3"]}]], 
+                    FractionBox["py", "q"], "]"}]}]}], ")"}]}]}], ")"}]}]}], 
   InterpretationBox[
    DynamicModuleBox[{Typeset`open = False}, 
     TemplateBox[{"Expression", 
@@ -1928,8 +2158,9 @@ Cell[BoxData[
     And[$CellContext`k + Im[$CellContext`py] != 0, Re[$CellContext`py] > 0]],
     SelectWithContents -> True, Selectable -> False]},
  "ConditionalExpression"]], "Output",
- CellChangeTimes->{3.869631098566168*^9, 3.8696311495711813`*^9},
- CellLabel->"Out[99]=",ExpressionUUID->"e2a7570f-7b25-4460-9f7b-a0b86b653c41"]
+ CellChangeTimes->{3.869631098566168*^9, 3.8696311495711813`*^9, 
+  3.8698129177682133`*^9, 3.8698153302243233`*^9},
+ CellLabel->"Out[29]=",ExpressionUUID->"31429318-32ae-4391-9ae8-abf4945da8e9"]
 }, Open  ]],

 Cell[BoxData[
@@ -1940,8 +2171,7 @@ Cell[BoxData[
    RowBox[{"kProton", ",", "228.465", ",", "py"}], "]"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.869631129704208*^9, 3.869631132293288*^9}, {
  3.869631164376565*^9, 3.8696311890788593`*^9}},
- CellLabel->
-  "In[101]:=",ExpressionUUID->"804a318d-428e-46d0-af1d-4cf0f2a8856f"],
+ CellLabel->"In[30]:=",ExpressionUUID->"804a318d-428e-46d0-af1d-4cf0f2a8856f"],

 Cell[CellGroupData[{

@@ -1949,8 +2179,7 @@ Cell[BoxData[
 RowBox[{"Print", "[", 
  RowBox[{"perp", "[", "py", "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.869631185940134*^9, 3.869631194927557*^9}},
- CellLabel->
-  "In[102]:=",ExpressionUUID->"3c17de4c-76db-4de9-9355-68bece89103c"],
+ CellLabel->"In[31]:=",ExpressionUUID->"3c17de4c-76db-4de9-9355-68bece89103c"],

 Cell[BoxData[
 TemplateBox[{
@@ -1961,30 +2190,31 @@ Cell[BoxData[
          SuperscriptBox["py", "2"]}], "+", 
        SuperscriptBox["py", "4"]}]], 
     RowBox[{"(", 
-       RowBox[{"1.090428163873036`*^8", "+", 
-         RowBox[{"2.962299845908087`*^6", " ", "py"}], "+", 
-         RowBox[{"54285.34889703492`", " ", 
-           SuperscriptBox["py", "2"]}], "+", 
-         RowBox[{"109.13809733550004`", " ", 
+       RowBox[{"85612.41941896101`", "\[VeryThinSpace]", "-", 
+         RowBox[{"1365.5989463371648`", " ", "py"}], "+", 
+         RowBox[{"42.62087326844503`", " ", 
+           SuperscriptBox["py", "2"]}], "-", 
+         RowBox[{"0.027166391551979432`", " ", 
           SuperscriptBox["py", "3"]}], "+", 
-         RowBox[{"1.0000000000000002`", " ", 
+         RowBox[{"0.0007851266342468513`", " ", 
           SuperscriptBox["py", "4"]}], "+", 
         RowBox[{
           RowBox[{"(", 
-             RowBox[{"8.247841088729748`*^7", "+", 
-               RowBox[{"41060.65360222363`", " ", 
-                 SuperscriptBox["py", "2"]}], "+", 
-               RowBox[{"0.75638555220682`", " ", 
+             RowBox[{
+               RowBox[{"-", "2578.122438137212`"}], "-", 
+               RowBox[{"1.283480019045514`", " ", 
+                 SuperscriptBox["py", "2"]}], "-", 
+               RowBox[{"0.000023643212121192028`", " ", 
                 SuperscriptBox["py", "4"]}]}], ")"}], " ", 
           RowBox[{"ArcTan", "[", 
             RowBox[{"0.004377038058345917`", " ", "py"}], "]"}]}], "+", 
         RowBox[{
           RowBox[{"(", 
             RowBox[{
-               RowBox[{"-", "1.5189722383410266`*^8"}], "-", 
-               RowBox[{"75619.78005994874`", " ", 
+               RowBox[{"-", "51924.436524199664`"}], "-", 
+               RowBox[{"25.849810619219497`", " ", 
                 SuperscriptBox["py", "2"]}], "-", 
-               RowBox[{"1.3930053245744014`", " ", 
+               RowBox[{"0.0004761839270527637`", " ", 
                 SuperscriptBox["py", "4"]}]}], ")"}], " ", 
           RowBox[{"ArcTan", "[", 
             RowBox[{"0.02187868155499205`", " ", "py"}], "]"}]}]}], ")"}]}], 
@@ -2057,10 +2287,10 @@ Cell[BoxData[
      0], Re[$CellContext`py] > 0], SelectWithContents -> True, Selectable -> 
    False]},
  "ConditionalExpression"]], "Print",
- CellChangeTimes->{3.869631197009533*^9},
+ CellChangeTimes->{3.869631197009533*^9, 3.86981533848575*^9},
 CellLabel->
  "During evaluation of \
-In[102]:=",ExpressionUUID->"a7a3cfc5-6831-43bc-a688-44c0eec53d7e"]
+In[31]:=",ExpressionUUID->"adbc5979-c07e-4cb9-9aa5-b7ff43f4bc1a"]
 }, Open  ]],

 Cell[CellGroupData[{
@@ -2068,83 +2298,94 @@ Cell[CellGroupData[{
 Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
-   RowBox[{"perp", "[", "py", "]"}], ",", 
+   RowBox[{
+    RowBox[{"perp", "[", "py", "]"}], "py"}], ",", 
   RowBox[{"{", 
    RowBox[{"py", ",", "0", ",", "300"}], "}"}]}], "]"}]], "Input",
- CellChangeTimes->{{3.869631214201191*^9, 3.869631226589404*^9}},
- CellLabel->
-  "In[104]:=",ExpressionUUID->"57f1c25e-794f-4adf-bd78-75639085813a"],
+ CellChangeTimes->{{3.869631214201191*^9, 3.869631226589404*^9}, {
+  3.869812855793929*^9, 3.869812867982595*^9}},
+ CellLabel->"In[32]:=",ExpressionUUID->"57f1c25e-794f-4adf-bd78-75639085813a"],

 Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    TagBox[
     {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[CompressedData["
-1:eJwd1AkwVf0bwHGREpIlecuS7EuuLUXkeZKs5SVRlhYlStYbKkWL5RJu4aZE
-kSRkS0lRJBGSVK8oSxSV7R77Fvn/+p+ZM2c+M9+Z35x5znPWHfLZdYSTg4Oj
-hNx/nw/t5yaG8j4bcvz/Goaomv0tNswS8DxdN7ewQMGVLxfPSDFfw2jMjh9/
-zVa6z30p5z3EZG3L/mtWD+vtYmYLFOl5aP315rSQq+f8OmBifbD3H+JhtUe0
-4Jxu+OylbzFPXMWXHb85pQdWu7kV/SZO7Ls5ORX7E7YEv0yYIT72Ot6x+Fw/
-FF8Ta5skNrjLKKf7DYFKombAGHHkVkOh65rDoPRflxVF7FeoWxWYMwJjmbO7
-+4kd12oH2MmOAX+7g38PsTFTTXFDyjgYbPqa1EFMm1f8LCw6CcnDfkXNxGKe
-MtEjsVOgqcVxv4GYo01iS9OSGVC8es73JXG/uRiVf24WOm4NUY+JPz4Ruh07
-/Rvst1ivyiF+pshv6+k3D4HbHjclE2cmLuG27P8DjT2b+GOIay+eb2NpcuBb
-janioL/vUxq8YTGbA5/HSpe5Ez8YCYo9kbMIp7UGeGyJFZRP/fjmxolD6Ukx
-BsQ3DgbALlkutJYL0pYjXnGdfr3yKxfOb/7Zv4w47J3PiEbKYuSXNksb+kPB
-zBIvi7S93BijJ7T1HbG3ocedFaJLMDTuVWk+8Z48Vzt27FLMldSj3Igbelzy
-91nw4O4Kk3AgRokDS98uWYZph1tzRIlVoh2e3j/Hiz/jGpxK5ym49dJeWNyA
-D48Ozr6MJBaZtT1+aZoPeRq3++0mnjtqJXnMbznuOx7p9H2OzCfNMrBVTQCF
-4j6fyCT+0WL2zrRfAIfZs9FuxEfv3MYETUG07NQO6PhNgQxID9VfFsSwrSjL
-Im7/ciuJky2IivZLo02JbYRTRv1yhNDMnWWUMUsBf/6a1GweYbw4M+Wzk7jG
-PMmy200YhzfuEh6bId/jhcQMa1kRXN4tQdcmHpcQtYm8IIKOUV3LG6cpKHiS
-MF/xVQRbVxu1uRLLDV+xV09ZifkxsbURUxR0Rq/gcp9diZwFMa+FiJMUmQW3
-9opi0+Tux9cnKRA4EM0jILoKI3+HytyYoGDqbXjpYKwYLnTOBBwbo8AVeHWF
-xsVwrc6y7dWjFHwoZBbrOP6DwpuXlYgT5ydcKwxRWI2fa5p0ngyT3iE7U7Bi
-DY4ucIfnDZG+Xk1OR14cF1kL27QPkvkYFN12iBbHd3rSStzE4mvLUtL3SKCn
-ytF8437S9zTEbxiWRCywVjz9g/T2NoIO9lL4PP2nXWAvOe91c2zwMymU8tzJ
-6dVDQVROZ2RN5FqctTKLM/lGgbTuKci9Io3fBiPLsjooaPHdHFJeJY3nSr/z
-+rZTwMyee9Y0KY0hYs72Gm1k3msubJ5wXoeRlFNGcivp56N0DJVlMKB0d5D4
-R9Jv3OFv7SyDXTxLlma/J/voI/Dw0GUZ3NNoska9iYKirngNxoQMZlFuN+Xe
-kr4qWbWpUhbv27SmPqyhYPvcPo9v47JYaerUx1VNztsgnT2uKIdsk/PWO6oo
-8MjMUFjNlENVUwHT6grSR+WtO+Qoj2dVNsXYPyF9pc9B/1h5zNR0PuX2mJw3
-q5ka8UIeVR+pK3g/Iu/r+VjivoICBj96l+JaSHqrCrHxUQW0m1rcNniPgocr
-3y+PiFbC54de+R26SvbdT7xHp0wJ+RKfKa5MoOBV45HS3n4ljH13QqP8CgWf
-Imfdtlso47QlS+Z3NAXTc7IvFi9TQasm1pDYRQoMewPpoRGqeDaAVtnqSYGp
-UaWZ9mNVjBjw117uQYF1Kt/a772qSMtyyNVzp8DFIe2N0fb1uGRyT0aQC9n3
-t3VynIvVcEe/sH2KHQX1xRIt5y/Q0GyTdbizAZm/sHuexgMaXtt1WExQj4I2
-nwehXV00fLTq1aUyHQoGlU00ENUxe2CD0CyN/E9u+UQt/FHHcZoKQ3Qd2fdt
-3ZJQroHfr73rrOaiwFd1Q9NzfS3kfih+QaiKDZ2bTBckbbXw30suIS3lbNhp
-7EgL8dDC3gCuhqulbFDZdz7GMEkL1VKow78fsKGX2WBWMamFp9JtOjzS2OA4
-6lr54oE2Vp0c5Ow8w4ZtT68WVSnq4J+iBIaWGhuKqrO6ZUEHj5iILuNVYoP0
-hzLBMHsddO6gu32RYcOf/m5v43AdZGeOvjnyDxueSNBUq7/p4Eq/kF+KXGxY
-f77mTs3NjRgfpt6yv2UIVplOsepEdFG8gu5TEzQEo7WW03Jquhg8ue74Nv8h
-eGuR5nTeRBerTVU5nnoNQZiVucym07pYOv5EknmQ9HbJBRmdusjcIG7YvZ30
-rlB/LlsPy+6Mz2uvIP1FxoIO6OPYnLl1c8ogjD0X80j3MMSttql1fpkDEH5A
-6p9rKkY4IqRfwMztg9zYA6HbDmzHWflPyoFuP0H8tLPndJgZUuZnG+eVe+Fe
-dtvbCH9LTF2r3pOp+B0aGVElZ+5YofuXnTz1Et1wqbmdX77RGs3i9ta6lXRC
-U/WHs4LbdiFXc/tHGqsNevi7NApabPFKg5vyvxdbgRFIu73Zxg5jCrdvjRpp
-hqmNLLvT/fY41ffJ8FfgR+D71GJmcHYvUr94OaT+fQ/515ynB0YdMCGmyzuR
-7x1IFU4PWEU5YV2OCss8+Q1IDjQn31q8D2+IWVPSVbXQf7b+8u3U/fhSReyY
-GU8N5AVaubBWH8QBacZkaU4VqGHr3sbmg3jcP6zatLgSRIq9Rrtvu+CI/ITP
-1rsVsLVMtv/4wUN4wGXou374M6gsNFukK34Yu9dZtNbUPoV7NqdLO3oP4xO7
-ups9siXQHFr/fibdFUe/hVaojDwCDxfFaMrnCGYt8Q2MfV8ENrSxw4p6bqid
-sbVdcKoQzsV8vdE87oaXpXz7vD/lw8rElhP7XrhjrnemNt+n+zCsMZzjFHIU
-d+0AgT7bbJAqO7koRfcYcrCbTnFI3oOHw7yzxZweqO+eNqwzngGXS8qtTSo8
-cH5Y1DxR7g58UA4xL4w+ji1hNh8vHLwNzmr7065YeGKpv802m7pU0JXOW/WQ
-2wtL5U+MmSrehDPlS9ULmrzwspXD0/H4G8DNvnh3HdMbR/LHAxXcrsNk7bJy
-LycfbGd2SN4xToSUuoXZYGFftFzT8uC4CwtcpnhNzvv64voTc2ahtfHweVb8
-kM8nXzzBmqvk3xgHEyZb8pJ0/ND4Q1l9SeVlKOVdxdMb74ePjSbG+Q/EglbG
-uLbktB+K7tW84sQVDVYR35lX7OlIO3Wk7/uiKDC1YO7YuJeOeXEya+xmIgEF
-9HjbHehosl76ft1wJGgmMsMV99ExK7H357OvkbAyU+9s+WE6yqT/SX/zPBLa
-qi8fG/Slo3cu+7/IoEg4ym1gbBZNx4HfkjMnJxngUvdjETuGjqw/15x2sRng
-GBtXkcCkY/LgCKX+gwE7RH9u/hpHxxdf72WNNzNAXT5eM/A6HXlrRrUyihkw
-afxLKuMuHbulhk7cCGAAxZPQbn6PjheqbHbe82LAr4YtN6gsOhrlehmVHGHA
-590Jovq5dDxweG9Qlx0Dnrka8n8ooqPPpZnZPRsZUKzUV3fyER0tB0VigmkM
-yB9MYEg+puNokSAtU4EBaf59nMee0tHVyTaVQ4wBSbqsFwJldLzaMHBcawUD
-4ucMQx49o6NINZ+x+1IGXHrRp+9YTsf/lPYrpC5EQGgYa2ahgo72IUyhL1MR
-8D8X86yZ
+1:eJwVlHk41PsXx5GuLZKUFvu+3BjLjLF9z8lepCiyZ8uWYobs8x2h5IqytFCh
+ZEsU0qYs7RSVFtxU3EgLkqyD/Ob3x+f5PK/neX/Oc855f85R8At32svHw8PT
+xT3/v+tdFqZGq3sJXcFNumalBkT6Y+9ux6ybkB7Unb3+swSc+Dc5QTbrCQRW
+PinteasIeYN5HfxZ3XDrVJHiyUY9OPXt/PRM5jA4Cf42o6SZg1SYYsavzBmw
+deNflzWCXZXBziPZQrgFH+1t6yDLzwf9K3x2iqIru//8dn0zBdQeo9Ax19C
+OGy6+sYBGz/QzHC7XcUWRr9hFbt0M39YCHaQCWGIouDJL/wq8nvBUeLcBOPy
+KuyOkQnVbw6GmY7Dd0YypXCocTDoS2A4BIAwfdWkFGZNCrrrvg6HrmtZDVT3
+dfjr0u8kO7kIqMk9fY1UXY+pj1sd3dMjIMCtsky8eQP2fbaJUXZlQNfg8xyD
+cRnUExseuNLPBHRxFHdzkcWh6cg37SsjoebJ20zWXVlkzF0ciyciIf3yx6OP
+j8phA6vtR05BJMjTY+HKCXlMOTMQWbUtCroX06mEhiJuKOMpfl50ELJo9lE7
+PBVx6EerNe3pQbAMF6v3O66Iav7MG4bjB6GuP4eSNqWIpeHKzGUQDVkPzmq9
+bFXC1OhmobbuaLBKr1bwc1fBMSkXWsB0DCy0hvtEZarggdr/Pseui4U6jm7R
+kRYV3PhBsFzLOBbkw25IV6mqIuVNoeT+xFhYcGiWmpxQxS2j+tR9C7FQL/lK
+9EiGOjqI6lrWf4+Du4yNg9RGdXz2quoTz/J4eNi5987Qd3XcGaDDvyQbD++O
+cgKttmrgwBY/ZzGneJhdUGrhF9LEDwnmCinX44EYimamHNFCpkOC29p9CWBj
+3mqrf0ML91rS759OSoAdRSJyn4e08Ev7haQXJxPA1634mbnV39joWhbj35wA
+qR1tynz8mxDFTpWqrUyE9gbp7qRD2ihy/fbKqAuJ0CURVE2p1carOVRIq0uE
+9+G1Kf392ljyKui6y4NEGNGwpiDqoGyGwpjWYCKsLAxPX/qjgyuKnVbfUWTB
+scW+q7ulKDhZUz54W5cF4n/zndXRoeD6ff6zsciC9en2jE/eFDzX4bCX5sUC
+TYsBGWiiYEut4JNbuSyoYfwlJPWOgtmf7srcK2aBfrHW5NgoBX8eDKhLqWaB
+6eLB9kIZXfR4s1vC7TELtt0UjvmTqIuFdfISOlMs6BrS8XuXp4uty7WatvOQ
+4CLpvK3mii7any8qooqQ4M0oUvLu00VHo7pEJ3kSIrQMXt4z0UOVP25j/9qQ
+8NHQZklmpx5u/jRnwetIwjZLd20yVA+Zhbu+/XIjQdMr6RiRr4dNjLAU/X0k
+DGU9t22e1sOXd+OyNmeQsOvspxg5MX2Muqez8CCHhPvlE2VsFX301dSMXl1A
+QnHLuuW4Sx+Vfmt931BBgvtEQGtLrT6+CXWyedtKwtM/MePybfq4tR0z256Q
+YCiSIXeoXx9NR22zMztIkFSuZW1eaYCTinPqR3pISKY8vHJB1QA50W2BjR9I
+GXa/Z6XMMAPLnw3W/8jodN50eh+mAGKvq7vNx8hwdRPPEQx1QBhstSvaZyE
+qgNKZ5LPGmCV8ITNiikSjqZtmTZvN0CRbOkM+UUSZnI9VUoGDJDRUnV8mIcN
+gcXhu5bNGeCbvknLZH42WNw+WfdAjYop7w6vMBBhQ92jigEloOKSyfkHdmJs
+kO9qFE91oeJy3RUH6avYcPxjJwzup6LqzkImZzUb/nwfOGB5mIqdGZ9Ec9ay
+IWxm8vylc1Q0C9h4ZmkdG94vE+zgv07FNE9PP4uNbLglra316D8qBhV1zjnI
+sUFVY7O7CoeK0jNtQVIKbDhJ3ZV+eBUNs+29dt9UZAO/edCtIXUaRsdYrtNR
+ZgPTIX7YCmkonD09R6qwod89c23Zbhru9im1K1dlw/agYqu/wml4uM3Us0KN
+DU2R9VGBR2hoj37xyeps+Dvpccnj8zTUMIuapGqwoeBYb5dqAw0jG9LXPeCy
+UP4Ib9pzGjacYhhoarIhtnSJMvyZhiWWdSf3cflLrYSPzTwNV2xIT0/jsnOT
+yvFyCUMsH/nmSXL5YTu9SUDTED/lkp47uKzXbTcatNkQdZllbxa58S989pZ+
+6mqI13ijV6VxWXycYaceYYg6/je8R7j5kQup8UfTDFFgfFBAm8ujgmcqvxYa
+otoLPVd7bj2ea6p6bG8YovWx1AIbbr3tCk0ClR2GaFfHWC/P7YeR9iua0JAh
+GgZbyr1TYkOF8eDekAVDfC4+NxXC7edam5m8ttV07B8xG34vz4aJp3azypvo
+qKDJkf6L60fH1mKPJGs69pSfUBKS4b5/Ptn0fg8de0Ve7x7awIZUhy2KhnF0
+vNetvPYs11+fl+cP5+TQ8dx4/B1Nrv8mThNfR6voSNXyKM7j/o8J57NXL32k
+o9q+knOzomzo7P4pwTNLR6/R4aDfwmyodLOM9lhlhCXNBz62CXDjeY2aSlga
+Yc+sop4ILzefAGhnVxrhma17ZL79IqHiS+6m9/eNkH1u+c/2URJSg7+eoPUZ
+YX/NGs2MbyQYh2XvHhUzRsmxP3oV/Vx95OdB9yhjPP0q+I0Xd/5Spwxtb2Qa
+o9WLGbrpUxL2xByrWlVujOYWxhXz90lYm0BlPu01Rul7b8Zkb3H1yWlLVDDB
+Us08zlIxCT78fX7ZriaYq6966id3H5gcoTweYZigfEdUy/08EibSe4+VXDJB
+3vkTHMmjXH221oZVwqY4nLSUwNjPnd/iFwYjb0xxQ8X0d3t9EtYpKp+xGTNF
+rQV+ar4WCb9LYucvCpjh2Z78ex1KJFwuV7jvZmyGO78Zl42v5uqvRm5/UmSG
+lwNE2Cq/WfD7nlToxVAC2+IU5dK4+3XWkZJvn0ygzgXTFatLWbAwZPt0Op/A
+mGmeycxzLOAXi1e1ayNQunMp2TGDBWu9+/77rQbo6GtHXxPMAvqfYg+rIcDS
+3qv1d2RZAHInVSevIm6r3+DUfigRDu+RXXda0xw7e+M8XdQT4ErmnhSLPVbo
+ndMRt3oiBjbGeYbNptriSPX+xfJHUVBe+b7jSJQd6jbwaGbyMaEzLf1mQokD
+RpF0e/XKA/DP274VKp07cE3jNoVOu1B4+agrUdzCCUNmIjN6TgfC4Ip+ytXu
+nSjm4KFOaPhDWrT2BWNHZ7wmN6VkxfCBGVqec9x3FzxpbrDYMuQBIu+6bU0T
+XTEgL7T/7pPdUHPac/bHhBt+qbslO3lkF8hem/3hkO6BSlumKpPTHUHmx9uz
+hfxeWGQa7CPm6gDfE9uPXyjyxpReyddP/myB6mgH37z1PnhUYz1/nJ81bMIe
+1863PkgVuCdclmwBqxv2Twxc8MVWc4GWvh6EzY1K3/f5+KFmQeb8wwgzaL1m
+y0vf6I+Wdy9fCek2gnLHuDsfhvxxqdYubDieBm9T2l/NXQxA5ZmeyunH+hDq
+q5bxM3wv3qQsXnoOuuCo/dtfzSgQj4n6h3400gb2sU8FbycDcZJnh2vCQ02Q
+PNUd6dUShK/bazfyJavBOGX8sgcZjP6zpZGMPmWQbYzhPUcPQWvOra+ipxSh
+flyY08AXis6xkQfp++Xh+M2mHdbNoejQpqnqsFMGujTILdcy9mEAFibLzK0H
+z03exSe2hmGVy8AVVWMpoMtXr61fvh+tLMS6JXolIaFJQOfqy/0YNJWpF5Mj
+AcvHkksVsg5gu236y+I6cZh+KtS03yMcPT7WDhzqFIVzbUsclkQEhh2yko6o
+FwbfGWHrpIgI1PeND6IoCEIvZ6Nf+LsItLu4uXxN+XKYsjarzqcyUDv4ZcGL
+M8vgjvBawaEcBq5zmxWWfcYLepcm9WVmGVicf/H4YxEecDjyOeuECxM9rPK3
+n9i0SNhszbKnuTLRuYAvABUWCRQzEu5zY6LgAz7PaclFQvdU1mE1LyY+TTH/
+HrmwQEiWGSU2+TPRa+qoesOzBeL9o+MhIxFMFFoKid22b4EIXm5qaZvBxJ2K
+icotNfOEb9sX3rFjTEwqy7qnWjJPuGdmN+dmMTE6MaE35/Q8Yb9m2PhTNhOt
+dxBZCUnzhI5Kjm70GSbO3NxekOU0T0xbfpW9VMrEraWWt2/NcIifgrl9W8qZ
+FODUI0e4RBfn5sV/Kxgoouz5Vf6AIfo3ZW7xuQKE1l0M5/Odg5xN4BY0VXH
+xOygslf0Qg7RoP6tLeY6E023t5yRz+UQNSO5aTI3mLgn8nyh6FEOURz1jS/k
+NhNVKgo8ZhgcIp+e1yLWyI2XbyI6GcghchYI8vpdJhZIX5iY8uAQ/7R8M3Fv
+YqJZ7kqhxR0cIiU1b26pmYmBJh3bhKw5xP8AUvjwbw==
       "]]},
-     Annotation[#, "Charting`Private`Tag$708369#1"]& ]}, {}},
+     Annotation[#, "Charting`Private`Tag$122611#1"]& ]}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
@@ -2176,7 +2417,7 @@ yB9MYEg+puNokSAtU4EBaf59nMee0tHVyTaVQ4wBSbqsFwJldLzaMHBcawUD
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
-  PlotRange->{{0, 300}, {0., 1.}},
+  PlotRange->{{0, 300}, {0., 0.007608918611885073}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
@@ -2184,9 +2425,10 @@ yB9MYEg+puNokSAtU4EBaf59nMee0tHVyTaVQ4wBSbqsFwJldLzaMHBcawUD
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
- CellChangeTimes->{{3.869631219529129*^9, 3.869631228317721*^9}},
- CellLabel->
-  "Out[104]=",ExpressionUUID->"6204eb92-3b0b-4108-b9ed-f12beecf3760"]
+ CellChangeTimes->{{3.869631219529129*^9, 3.869631228317721*^9}, {
+   3.869812862920465*^9, 3.869812869354279*^9}, 3.869812953045562*^9, 
+   3.86981534062702*^9},
+ CellLabel->"Out[32]=",ExpressionUUID->"92946697-5ec6-4805-ba6b-0310b33ef2d5"]
 }, Open  ]]
 },
 WindowSize->{1389.75, 768.75},
@@ -2207,89 +2449,93 @@ CellTagsIndex->{}
 *)
 (*NotebookFileOutline
 Notebook[{
-Cell[558, 20, 1512, 28, 109, "Input",ExpressionUUID->"15f55500-85b0-49bd-b8c0-7ea00f0e2245"],
+Cell[558, 20, 1511, 28, 109, "Input",ExpressionUUID->"15f55500-85b0-49bd-b8c0-7ea00f0e2245"],
+Cell[CellGroupData[{
+Cell[2094, 52, 811, 21, 59, "Input",ExpressionUUID->"f1c72a75-7ee7-46f9-a8f7-4cd49fd00c7a"],
+Cell[2908, 75, 574, 11, 37, "Output",ExpressionUUID->"242a245a-39df-462d-8951-540e89edf802"]
+}, Open  ]],
+Cell[3497, 89, 1047, 24, 57, "Input",ExpressionUUID->"40ae1148-daf9-4e15-8762-688f07195dad"],
+Cell[4547, 115, 1121, 27, 57, "Input",ExpressionUUID->"7a2bcc4d-5447-48d1-ac24-841589704be5"],
 Cell[CellGroupData[{
-Cell[2095, 52, 812, 21, 59, "Input",ExpressionUUID->"f1c72a75-7ee7-46f9-a8f7-4cd49fd00c7a"],
-Cell[2910, 75, 548, 10, 37, "Output",ExpressionUUID->"787a36b5-5ebb-450b-97f8-6a66d36d74d3"]
+Cell[5693, 146, 260, 5, 33, "Input",ExpressionUUID->"b71993bf-44b9-4b50-90b9-025a29b6f35d"],
+Cell[5956, 153, 465, 12, 54, "Output",ExpressionUUID->"c3aed629-a6fc-486e-8683-c6b9ccef00bb"]
 }, Open  ]],
-Cell[3473, 88, 1048, 24, 57, "Input",ExpressionUUID->"40ae1148-daf9-4e15-8762-688f07195dad"],
-Cell[4524, 114, 1122, 27, 57, "Input",ExpressionUUID->"7a2bcc4d-5447-48d1-ac24-841589704be5"],
+Cell[6436, 168, 942, 18, 59, "Input",ExpressionUUID->"ed029352-16a6-45e1-b42a-177020211247"],
 Cell[CellGroupData[{
-Cell[5671, 145, 261, 5, 33, "Input",ExpressionUUID->"b71993bf-44b9-4b50-90b9-025a29b6f35d"],
-Cell[5935, 152, 440, 11, 54, "Output",ExpressionUUID->"36448c3c-ed58-4b87-bca0-c6f7044cd75e"]
+Cell[7403, 190, 273, 5, 57, "Input",ExpressionUUID->"1508f981-ab16-4888-9c9d-70043439577a"],
+Cell[7679, 197, 696, 21, 74, "Output",ExpressionUUID->"e4db3020-f165-46b8-bef6-9b39b4225155"]
 }, Open  ]],
-Cell[6390, 166, 943, 18, 59, "Input",ExpressionUUID->"ed029352-16a6-45e1-b42a-177020211247"],
 Cell[CellGroupData[{
-Cell[7358, 188, 274, 5, 57, "Input",ExpressionUUID->"1508f981-ab16-4888-9c9d-70043439577a"],
-Cell[7635, 195, 672, 20, 74, "Output",ExpressionUUID->"d7839f75-9aed-4efb-95bf-f0aa05801eb8"]
+Cell[8412, 223, 2346, 58, 114, "Input",ExpressionUUID->"cd919aeb-3101-4893-8360-ca90a34b5f8a"],
+Cell[10761, 283, 752, 11, 37, "Output",ExpressionUUID->"86a70908-b758-4cb1-8f2d-c9ea5b1964b6"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[8344, 220, 2347, 58, 114, "Input",ExpressionUUID->"cd919aeb-3101-4893-8360-ca90a34b5f8a"],
-Cell[10694, 280, 729, 11, 37, "Output",ExpressionUUID->"b0e60b84-1f1f-4268-863d-06c8900b3b0b"]
+Cell[11550, 299, 2475, 66, 155, "Input",ExpressionUUID->"d1262ffb-87ed-4322-8f05-1b89384f8939"],
+Cell[14028, 367, 3450, 75, 249, "Output",ExpressionUUID->"01ca38e8-4fdc-4873-bac0-259fd0f0d105"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[11460, 296, 2475, 66, 155, "Input",ExpressionUUID->"d1262ffb-87ed-4322-8f05-1b89384f8939"],
-Cell[13938, 364, 3450, 75, 249, "Output",ExpressionUUID->"01ca38e8-4fdc-4873-bac0-259fd0f0d105"]
+Cell[17515, 447, 921, 24, 50, "Input",ExpressionUUID->"cab76d12-2dad-42c6-94b6-db201a8c18ec"],
+Cell[18439, 473, 336, 6, 37, "Output",ExpressionUUID->"50db03d0-7166-4dd1-b1ad-8af23d2f754f"]
 }, Open  ]],
+Cell[18790, 482, 1109, 22, 37, "Input",ExpressionUUID->"363c9fee-04cf-4e5c-90ea-4a0da118130e"],
 Cell[CellGroupData[{
-Cell[17425, 444, 921, 24, 50, "Input",ExpressionUUID->"cab76d12-2dad-42c6-94b6-db201a8c18ec"],
-Cell[18349, 470, 336, 6, 37, "Output",ExpressionUUID->"50db03d0-7166-4dd1-b1ad-8af23d2f754f"]
+Cell[19924, 508, 1277, 35, 105, "Input",ExpressionUUID->"ec406fd6-54f4-4e5f-8e73-4494b58edd50"],
+Cell[21204, 545, 6098, 118, 260, "Output",ExpressionUUID->"2bd6fbb1-9564-433a-b163-c5529f0a682d"]
 }, Open  ]],
-Cell[18700, 479, 1109, 22, 37, "Input",ExpressionUUID->"363c9fee-04cf-4e5c-90ea-4a0da118130e"],
 Cell[CellGroupData[{
-Cell[19834, 505, 1277, 35, 105, "Input",ExpressionUUID->"ec406fd6-54f4-4e5f-8e73-4494b58edd50"],
-Cell[21114, 542, 6098, 118, 260, "Output",ExpressionUUID->"2bd6fbb1-9564-433a-b163-c5529f0a682d"]
+Cell[27339, 668, 4519, 106, 187, "Input",ExpressionUUID->"589fa93e-41d0-46c0-b8d3-ab8df683f53d"],
+Cell[31861, 776, 1095, 33, 82, "Print",ExpressionUUID->"337280b6-e158-4dad-b3cf-8ba39c1d5cf3"]
 }, Open  ]],
+Cell[32971, 812, 363, 9, 57, "Input",ExpressionUUID->"58e80c3d-daed-4f30-890a-2e8e52e60c0c"],
 Cell[CellGroupData[{
-Cell[27249, 665, 4519, 106, 187, "Input",ExpressionUUID->"589fa93e-41d0-46c0-b8d3-ab8df683f53d"],
-Cell[31771, 773, 1073, 33, 82, "Print",ExpressionUUID->"2946d86e-f35b-4b75-b95e-a5caa077a51c"]
+Cell[33359, 825, 2744, 65, 124, "Input",ExpressionUUID->"c95e48ba-da31-4186-9fab-e23d54171786"],
+Cell[36106, 892, 406, 6, 37, "Output",ExpressionUUID->"6beae27f-a003-413e-be5f-41b65ba6a545"]
 }, Open  ]],
-Cell[32859, 809, 363, 9, 57, "Input",ExpressionUUID->"58e80c3d-daed-4f30-890a-2e8e52e60c0c"],
 Cell[CellGroupData[{
-Cell[33247, 822, 2744, 65, 124, "Input",ExpressionUUID->"c95e48ba-da31-4186-9fab-e23d54171786"],
-Cell[35994, 889, 381, 5, 37, "Output",ExpressionUUID->"605d794b-f8c9-4677-b077-aada9b92bbd5"]
+Cell[36549, 903, 2335, 57, 114, "Input",ExpressionUUID->"af705252-b6bb-4c46-94eb-0f4dbbf1db66"],
+Cell[38887, 962, 7577, 142, 253, "Output",ExpressionUUID->"e690fd3a-0037-42ec-aaf7-a005d8ba954a"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[36412, 899, 2238, 56, 113, "Input",ExpressionUUID->"af705252-b6bb-4c46-94eb-0f4dbbf1db66"],
-Cell[38653, 957, 6896, 131, 256, "Output",ExpressionUUID->"e98c5906-814e-490d-b5ad-0e56295f99ab"]
+Cell[46501, 1109, 1501, 42, 85, "Input",ExpressionUUID->"c924a6e4-1400-4616-9786-41ee8ec594ab"],
+Cell[48005, 1153, 7778, 148, 277, "Output",ExpressionUUID->"dd57294b-38b9-4e1e-a78f-a790c7d844ec"]
 }, Open  ]],
-Cell[45564, 1091, 203, 3, 33, "Input",ExpressionUUID->"1e43fbb6-0494-4504-bdc4-df71c6e9d5e3"],
-Cell[45770, 1096, 203, 3, 33, "Input",ExpressionUUID->"b34abcdd-b1ba-4557-8728-df3df5b28b02"],
-Cell[45976, 1101, 945, 30, 82, "Input",ExpressionUUID->"e24edc1d-a3ff-4c8b-bb21-62fa5fe16f06"],
-Cell[46924, 1133, 205, 3, 33, "Input",ExpressionUUID->"dfd70dbc-5fed-4e5c-a039-a3b3f88789e3"],
-Cell[47132, 1138, 1241, 36, 86, "Input",ExpressionUUID->"f53dd42a-ae11-462a-bf9f-c07d88f298c5"],
+Cell[55798, 1304, 203, 3, 33, "Input",ExpressionUUID->"1e43fbb6-0494-4504-bdc4-df71c6e9d5e3"],
+Cell[56004, 1309, 203, 3, 33, "Input",ExpressionUUID->"b34abcdd-b1ba-4557-8728-df3df5b28b02"],
+Cell[56210, 1314, 945, 30, 82, "Input",ExpressionUUID->"e24edc1d-a3ff-4c8b-bb21-62fa5fe16f06"],
+Cell[57158, 1346, 205, 3, 33, "Input",ExpressionUUID->"dfd70dbc-5fed-4e5c-a039-a3b3f88789e3"],
+Cell[57366, 1351, 1201, 34, 86, "Input",ExpressionUUID->"f53dd42a-ae11-462a-bf9f-c07d88f298c5"],
 Cell[CellGroupData[{
-Cell[48398, 1178, 283, 5, 33, "Input",ExpressionUUID->"417870a7-4839-4f81-8586-d217e9ba44f9"],
-Cell[48684, 1185, 872, 27, 50, "Print",ExpressionUUID->"bf7a1e53-2a8c-4713-ad3f-6433be47849b"]
+Cell[58592, 1389, 283, 5, 33, "Input",ExpressionUUID->"417870a7-4839-4f81-8586-d217e9ba44f9"],
+Cell[58878, 1396, 811, 25, 50, "Print",ExpressionUUID->"5de6100a-eac1-4d61-85e4-48f2515ee2f3"]
 }, Open  ]],
-Cell[49571, 1215, 744, 16, 57, "Input",ExpressionUUID->"dffe0abc-3dc5-4360-a36a-ddf70a8a91e3"],
+Cell[59704, 1424, 744, 16, 57, "Input",ExpressionUUID->"dffe0abc-3dc5-4360-a36a-ddf70a8a91e3"],
 Cell[CellGroupData[{
-Cell[50340, 1235, 1868, 45, 60, "Input",ExpressionUUID->"a4ef5eb4-7778-43b2-8376-e6e4431de66f"],
-Cell[52211, 1282, 6867, 169, 83, "Output",ExpressionUUID->"63a40424-d528-4ecd-bf57-728421a50da7"]
+Cell[60473, 1444, 1810, 42, 60, "Input",ExpressionUUID->"a4ef5eb4-7778-43b2-8376-e6e4431de66f"],
+Cell[62286, 1488, 6410, 160, 110, "Output",ExpressionUUID->"bd741071-db04-43f7-98ef-6e9e741872b6"]
 }, Open  ]],
-Cell[59093, 1454, 404, 9, 33, "Input",ExpressionUUID->"d6114d76-496a-4526-9235-2e2b26389c18"],
+Cell[68711, 1651, 422, 8, 33, "Input",ExpressionUUID->"d6114d76-496a-4526-9235-2e2b26389c18"],
 Cell[CellGroupData[{
-Cell[59522, 1467, 240, 4, 33, "Input",ExpressionUUID->"ea14538f-b536-4fe7-a74f-f5e25d13093b"],
-Cell[59765, 1473, 3901, 90, 43, "Print",ExpressionUUID->"5d9c5e6a-5231-4a6e-a4f7-21f99a3e5607"]
+Cell[69158, 1663, 240, 4, 33, "Input",ExpressionUUID->"ea14538f-b536-4fe7-a74f-f5e25d13093b"],
+Cell[69401, 1669, 3940, 89, 63, "Print",ExpressionUUID->"ad5eddd5-5fad-4d49-aef0-c2ad504aea37"]
 }, Open  ]],
-Cell[63681, 1566, 802, 17, 33, "Input",ExpressionUUID->"9f6413c0-e44e-4126-b089-740d915bb48d"],
+Cell[73356, 1761, 802, 17, 33, "Input",ExpressionUUID->"9f6413c0-e44e-4126-b089-740d915bb48d"],
 Cell[CellGroupData[{
-Cell[64508, 1587, 323, 7, 33, "Input",ExpressionUUID->"37c239c3-4ac2-401f-963b-6c070eeb2f73"],
-Cell[64834, 1596, 3712, 80, 284, "Output",ExpressionUUID->"030448a6-0fbd-4cf2-a9d9-83eed4c20ddd"]
+Cell[74183, 1782, 473, 9, 33, "Input",ExpressionUUID->"37c239c3-4ac2-401f-963b-6c070eeb2f73"],
+Cell[74659, 1793, 3854, 81, 253, "Output",ExpressionUUID->"e052c0db-eb44-4443-8d7c-14a4d951c77f"]
 }, Open  ]],
-Cell[68561, 1679, 557, 13, 33, "Input",ExpressionUUID->"cc447812-1c6e-4129-943e-0459fce59f1e"],
+Cell[78528, 1877, 557, 13, 33, "Input",ExpressionUUID->"cc447812-1c6e-4129-943e-0459fce59f1e"],
 Cell[CellGroupData[{
-Cell[69143, 1696, 1377, 38, 60, "Input",ExpressionUUID->"b86a339e-b949-4516-8eb5-8294678723fc"],
-Cell[70523, 1736, 7831, 195, 122, "Output",ExpressionUUID->"e2a7570f-7b25-4460-9f7b-a0b86b653c41"]
+Cell[79110, 1894, 1299, 35, 60, "Input",ExpressionUUID->"b86a339e-b949-4516-8eb5-8294678723fc"],
+Cell[80412, 1931, 9356, 231, 150, "Output",ExpressionUUID->"31429318-32ae-4391-9ae8-abf4945da8e9"]
 }, Open  ]],
-Cell[78369, 1934, 401, 9, 33, "Input",ExpressionUUID->"804a318d-428e-46d0-af1d-4cf0f2a8856f"],
+Cell[89783, 2165, 397, 8, 33, "Input",ExpressionUUID->"804a318d-428e-46d0-af1d-4cf0f2a8856f"],
 Cell[CellGroupData[{
-Cell[78795, 1947, 239, 5, 33, "Input",ExpressionUUID->"3c17de4c-76db-4de9-9355-68bece89103c"],
-Cell[79037, 1954, 4590, 108, 82, "Print",ExpressionUUID->"a7a3cfc5-6831-43bc-a688-44c0eec53d7e"]
+Cell[90205, 2177, 235, 4, 33, "Input",ExpressionUUID->"3c17de4c-76db-4de9-9355-68bece89103c"],
+Cell[90443, 2183, 4666, 109, 82, "Print",ExpressionUUID->"adbc5979-c07e-4cb9-9aa5-b7ff43f4bc1a"]
 }, Open  ]],
 Cell[CellGroupData[{
-Cell[83664, 2067, 324, 8, 33, "Input",ExpressionUUID->"57f1c25e-794f-4adf-bd78-75639085813a"],
-Cell[83991, 2077, 5546, 111, 282, "Output",ExpressionUUID->"6204eb92-3b0b-4108-b9ed-f12beecf3760"]
+Cell[95146, 2297, 390, 9, 33, "Input",ExpressionUUID->"57f1c25e-794f-4adf-bd78-75639085813a"],
+Cell[95539, 2308, 6303, 122, 257, "Output",ExpressionUUID->"92946697-5ec6-4805-ba6b-0310b33ef2d5"]
 }, Open  ]]
 }
 ]