deutron.nb

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

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

(* CreatedBy='Mathematica 13.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[    107736,       2535]
NotebookOptionsPosition[    101858,       2433]
NotebookOutlinePosition[    102286,       2450]
CellTagsIndexPosition[    102243,       2447]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", "Input", " ", "*)"}], " ", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"kProton", " ", "=", " ", 
     SqrtBox[
      RowBox[{"2", "*", "2.225", "*", 
       FractionBox[
        RowBox[{"939.565", "*", "938.272"}], 
        RowBox[{"938.272", "+", "939.565"}]]}]]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"radCh", " ", "=", " ", "2.1421"}], ";"}]}]}]], "Input",
 CellChangeTimes->{{3.866683779231073*^9, 3.866683807571497*^9}, {
   3.8666838393818808`*^9, 3.8666838923039827`*^9}, {3.8666847300789223`*^9, 
   3.8666848011687326`*^9}, {3.866684835937563*^9, 3.8666848506595507`*^9}, {
   3.866685313816524*^9, 3.866685351704863*^9}, {3.866689549776765*^9, 
   3.866689552977665*^9}, {3.866691101712513*^9, 3.866691105512731*^9}, {
   3.866691177671338*^9, 3.866691183519827*^9}, {3.866691222726439*^9, 
   3.8666912265676317`*^9}, {3.8666912841399517`*^9, 3.866691327721827*^9}, {
   3.8666913674365187`*^9, 3.866691379354693*^9}, {3.866692223786235*^9, 
   3.866692243517498*^9}, {3.86669234088623*^9, 3.866692356098913*^9}, {
   3.8666925244941597`*^9, 3.8666925247684307`*^9}, {3.866789171218506*^9, 
   3.866789284836467*^9}, {3.866789484746875*^9, 3.8667895047278633`*^9}, 
   3.866791544874084*^9, {3.867203881456174*^9, 3.867203891918221*^9}, {
   3.86859415931714*^9, 3.868594159600313*^9}, {3.868594205588904*^9, 
   3.8685942118497667`*^9}},
 CellLabel->"In[1]:=",ExpressionUUID->"15f55500-85b0-49bd-b8c0-7ea00f0e2245"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
   "\:041d\:0430\:0439\:0434\:0451\:043c", " ", 
    "\:0440\:0430\:0434\:0438\:0443\:0441"}], " ", "*)"}], " ", 
  "\[IndentingNewLine]", 
  RowBox[{"Solve", "[", 
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{
       SuperscriptBox["radii", "2"], "/", "4"}], " ", "+", " ", "0.64", " ", 
      "-", "0.125"}], "==", 
     SuperscriptBox["radCh", "2"]}], ",", "radii", ",", 
    RowBox[{"Assumptions", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{"radii", ">", "0"}], "}"}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.8672039064993887`*^9, 3.867204004483068*^9}, {
  3.8672040350457573`*^9, 3.867204051044759*^9}, {3.868594165302346*^9, 
  3.8685941958762693`*^9}},
 CellLabel->"In[3]:=",ExpressionUUID->"f1c72a75-7ee7-46f9-a8f7-4cd49fd00c7a"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"{", 
   RowBox[{"radii", "\[Rule]", "4.036628499131423`"}], "}"}], "}"}]], "Output",\

 CellChangeTimes->{3.86738855414594*^9, 3.867475734034484*^9, 
  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, 
  3.8698124214068537`*^9},
 CellLabel->"Out[3]=",ExpressionUUID->"242a245a-39df-462d-8951-540e89edf802"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
   "\:0418\:0441\:043f\:043e\:043b\:044c\:0437\:0443\:0435\:043c", " ", 
    "\:0442\:0430\:043a\:043e\:0435", " ", 
    "\:043f\:0440\:0435\:0434\:0441\:0442\:0430\:0432\:043b\:0435\:043d\:0438\
\:0435", " ", "\:0432\:043e\:043b\:043d\:043e\:0432\:043e\:0439", " ", 
    "\:0444\:0443\:043d\:043a\:0446\:0438\:0438"}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"fl1", "[", 
    RowBox[{"k_", ",", "q_", ",", "r_"}], "]"}], " ", ":=", " ", 
   RowBox[{
    RowBox[{"Exp", "[", 
     RowBox[{
      RowBox[{"-", "k"}], "*", "r"}], "]"}], " ", "-", " ", 
    RowBox[{"Exp", "[", 
     RowBox[{
      RowBox[{"-", "q"}], "*", "r"}], "]"}]}]}]}]], "Input",
 CellChangeTimes->{{3.866787316082055*^9, 3.866787385489438*^9}, 
   3.866789623085449*^9, {3.8667897202033997`*^9, 3.866789728514553*^9}, {
   3.868594244912619*^9, 3.86859426617033*^9}, {3.868599535129788*^9, 
   3.868599540788164*^9}},
 CellLabel->"In[4]:=",ExpressionUUID->"40ae1148-daf9-4e15-8762-688f07195dad"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
   "\:041d\:043e\:0440\:043c\:0438\:0440\:0443\:0435\:043c", " ", 
    "\:0432\:043e\:043b\:043d\:043e\:0432\:0443\:044e", " ", 
    "\:0444\:0443\:043d\:043a\:0446\:0438\:044e"}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Norm1", "[", 
    RowBox[{"k_", ",", "q_"}], "]"}], " ", ":=", " ", 
   RowBox[{"Integrate", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"fl1", "[", 
       RowBox[{"k", ",", "q", ",", "r"}], "]"}], "^", "2"}], ",", 
     RowBox[{"{", 
      RowBox[{"r", ",", "0", ",", "Infinity"}], "}"}], ",", 
     RowBox[{"Assumptions", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"k", ">", "0"}], ",", 
        RowBox[{"q", ">", "0"}]}], "}"}]}]}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.866787389630982*^9, 3.866787451078424*^9}, {
  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[5]:=",ExpressionUUID->"7a2bcc4d-5447-48d1-ac24-841589704be5"],

Cell[CellGroupData[{

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

Cell[BoxData[
 RowBox[{
  FractionBox["1", "2"], " ", 
  RowBox[{"(", 
   RowBox[{
    FractionBox["1", "k"], "+", 
    FractionBox["1", "q"], "-", 
    FractionBox["4", 
     RowBox[{"k", "+", "q"}]]}], ")"}]}]], "Output",
 CellChangeTimes->{{3.868603908548182*^9, 3.8686039184646873`*^9}, 
   3.869565521238338*^9, 3.869625957956542*^9, 3.869628024261294*^9, 
   3.869812428243662*^9},
 CellLabel->"Out[6]=",ExpressionUUID->"c3aed629-a6fc-486e-8683-c6b9ccef00bb"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"PsiD", "[", 
   RowBox[{"k_", ",", "q_", ",", "r_"}], "]"}], ":=", " ", 
  FractionBox[
   RowBox[{"fl1", "[", 
    RowBox[{"k", ",", "q", ",", "r"}], "]"}], 
   SqrtBox[
    RowBox[{"Norm1", "[", 
     RowBox[{"k", ",", "q"}], "]"}]]]}]], "Input",
 CellChangeTimes->{{3.866790998002054*^9, 3.8667910749768953`*^9}, {
   3.866791375415957*^9, 3.866791421868765*^9}, {3.866791549994028*^9, 
   3.866791550606978*^9}, {3.8672042627726517`*^9, 3.867204264059602*^9}, 
   3.868423974748228*^9, {3.868510732445642*^9, 3.86851075423566*^9}, {
   3.868591785188673*^9, 3.868591787116836*^9}, {3.868591945385953*^9, 
   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[7]:=",ExpressionUUID->"ed029352-16a6-45e1-b42a-177020211247"],

Cell[CellGroupData[{

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

Cell[BoxData[
 FractionBox[
  RowBox[{
   SqrtBox["2"], " ", 
   RowBox[{"(", 
    RowBox[{
     SuperscriptBox["\[ExponentialE]", 
      RowBox[{
       RowBox[{"-", "k"}], " ", "r"}]], "-", 
     SuperscriptBox["\[ExponentialE]", 
      RowBox[{
       RowBox[{"-", "q"}], " ", "r"}]]}], ")"}]}], 
  SqrtBox[
   RowBox[{
    FractionBox["1", "k"], "+", 
    FractionBox["1", "q"], "-", 
    FractionBox["4", 
     RowBox[{"k", "+", "q"}]]}]]]], "Output",
 CellChangeTimes->{3.868598458698038*^9, 3.8685995629995413`*^9, 
  3.869565525592538*^9, 3.869625962223899*^9, 3.8696280292131033`*^9, 
  3.869812432237939*^9},
 CellLabel->"Out[8]=",ExpressionUUID->"e4db3020-f165-46b8-bef6-9b39b4225155"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{
    "\:041f\:0440\:043e\:0432\:0435\:0440\:0438\:043c", " ", 
     "\:043d\:043e\:0440\:043c\:0438\:0440\:043e\:0432\:043a\:0443"}], ",", 
    " ", 
    RowBox[{
     RowBox[{
     "\:0432\:0435\:0437\:0434\:0435", " ", "\:0433\:0434\:0435", " ", 
      "\:0443", " ", "\:043d\:0430\:0441", " ", 
      "\:0438\:043c\:043f\:0443\:043b\:044c\:0441\:044b", " ", 
      "\:0434\:043e\:043c\:043d\:043e\:0436\:0430\:0435\:043c", " ", 
      "\:043d\:0430", " ", 
      "\:043a\:043e\:044d\:0444\:0444\:0438\:0446\:0438\:0435\:043d\:0442", 
      " ", 
      RowBox[{"1", "/", "hc"}]}], "==", 
     RowBox[{
      RowBox[{"1", "/", "197.327"}], " ", "\:0434\:043b\:044f", " ", 
      "\:0440\:0430\:0437\:043c\:0435\:0440\:043d\:043e\:0441\:0442\:0438"}]}]\
}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{"Integrate", "[", 
   RowBox[{
    SuperscriptBox[
     RowBox[{"(", 
      FractionBox[
       RowBox[{
        SqrtBox["2"], " ", 
        RowBox[{"(", 
         RowBox[{
          SuperscriptBox["\[ExponentialE]", 
           RowBox[{
            RowBox[{"-", "k"}], " ", "r"}]], "-", 
          SuperscriptBox["\[ExponentialE]", 
           RowBox[{
            RowBox[{"-", "q"}], " ", "r"}]]}], ")"}]}], 
       SqrtBox[
        RowBox[{
         FractionBox["1", "k"], "+", 
         FractionBox["1", "q"], "-", 
         FractionBox["4", 
          RowBox[{"k", "+", "q"}]]}]]], ")"}], "2"], ",", 
    RowBox[{"{", 
     RowBox[{"r", ",", "0", ",", "Infinity"}], "}"}], ",", 
    RowBox[{"Assumptions", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"k", ">", "0"}], ",", 
       RowBox[{"q", ">", "0"}]}], "}"}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.86738924390305*^9, 3.867389265162883*^9}, 
   3.868592069790248*^9, 3.868592175306278*^9, {3.868592227741321*^9, 
   3.868592228962371*^9}, {3.8685922770324097`*^9, 3.8685922831167917`*^9}, {
   3.8685929499135513`*^9, 3.868592952952462*^9}, {3.868594341632305*^9, 
   3.868594381115343*^9}, {3.868594497762939*^9, 3.868594504907331*^9}, {
   3.868594547173224*^9, 3.868594552292612*^9}, {3.868598275497266*^9, 
   3.868598275773817*^9}, 3.8685994674679527`*^9, {3.868599511457546*^9, 
   3.868599518423562*^9}},
 CellLabel->"In[9]:=",ExpressionUUID->"cd919aeb-3101-4893-8360-ca90a34b5f8a"],

Cell[BoxData["1"], "Output",
 CellChangeTimes->{
  3.867389267153965*^9, 3.8674757368662777`*^9, 3.868362822870886*^9, 
   3.868364534825906*^9, 3.868364571308099*^9, 3.8683646862164707`*^9, 
   3.868418724949291*^9, 3.868424003660665*^9, 3.8684240538961067`*^9, 
   3.868429726397519*^9, 3.868503338909972*^9, {3.868510760161477*^9, 
   3.8685107836324863`*^9}, 3.868591740324675*^9, 3.8685919503322763`*^9, 
   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, 3.8698124352622843`*^9},
 CellLabel->"Out[9]=",ExpressionUUID->"86a70908-b758-4cb1-8f2d-c9ea5b1964b6"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
   "\:041f\:043e\:0438\:0441\:043a", " ", "q", " ", "\:043a\:0430\:043a", " ",
     "\:0441\:0440\:0435\:0434\:043d\:0435\:0435", " ", 
    SuperscriptBox["r", "2"]}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"Integrate", "[", 
       RowBox[{
        RowBox[{
         SuperscriptBox[
          RowBox[{"PsiD", "[", 
           RowBox[{
            FractionBox["kProton", "197.327"], ",", 
            FractionBox["q", "197.327"], ",", "rad"}], "]"}], "2"], "*", 
         SuperscriptBox["rad", "2"]}], ",", 
        RowBox[{"{", 
         RowBox[{"rad", ",", "0", ",", "Infinity"}], "}"}]}], "]"}], "-", 
      SuperscriptBox["4.04", "2"]}], ",", 
     RowBox[{"{", 
      RowBox[{"q", ",", "225", ",", "230"}], "}"}]}], "]"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"Plot", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"Integrate", "[", 
      RowBox[{
       RowBox[{
        RowBox[{
         SuperscriptBox[
          RowBox[{"(", 
           FractionBox[
            RowBox[{
             SqrtBox["2"], " ", 
             RowBox[{"(", 
              RowBox[{
               SuperscriptBox["\[ExponentialE]", 
                RowBox[{
                 RowBox[{"-", "k"}], " ", "r"}]], "-", 
               SuperscriptBox["\[ExponentialE]", 
                RowBox[{"-", 
                 FractionBox[
                  RowBox[{"q", " ", "r"}], "197.327"]}]]}], ")"}]}], 
            SqrtBox[
             RowBox[{
              FractionBox["1", "k"], "+", 
              FractionBox["197.327", "q"], "-", 
              FractionBox["4", 
               RowBox[{"k", "+", 
                FractionBox["q", "197.327"]}]]}]]], ")"}], "2"], "*", 
         SuperscriptBox["r", "2"]}], "/.", 
        RowBox[{"k", "->", " ", 
         FractionBox["kProton", "197.327"]}]}], " ", ",", 
       RowBox[{"{", 
        RowBox[{"r", ",", "0", ",", "Infinity"}], "}"}]}], "]"}], "-", 
     SuperscriptBox["4.04", "2"]}], ",", 
    RowBox[{"{", 
     RowBox[{"q", ",", "225", ",", "230"}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{
  3.8685946768723087`*^9, {3.868594803071945*^9, 3.868594825178236*^9}, {
   3.8685995856506357`*^9, 3.868599618767168*^9}, {3.8685997043552017`*^9, 
   3.868599788228649*^9}, {3.8685998237866592`*^9, 3.868599944132773*^9}},
 CellLabel->"In[54]:=",ExpressionUUID->"d1262ffb-87ed-4322-8f05-1b89384f8939"],

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    TagBox[
     {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[CompressedData["
1:eJwd0nkw3GcYB/BtjTREjJK4j90444hbxZFXZKRIpClCXW3dpnZtJcYV18Zi
3Ueomw2ydv0coSuo442SwxGt2JWhwRBJHKkRolPi6Pv2j2ee+czzxzPPM19K
AN01+HMSieSCCveAZxYkEplqR5JLCbCfgUCKo1cggc2eFCsVQgBTyWSF/+cD
/PVJCJSBuJ0Jdm/o6ZJRCISdc8lB2F8ffLbSC4FjPfPoGfa6PiuvBoJzSS8O
CynIDvGO84EQ7JpGHKirI3s0l7I2BgB1yaIyWwvZezHy281+wF65Gt1zFpll
KGIk0g+2iyO9X+ojO9NSNxT7wBmm9oKlEXJIYritWi+QHtO996Mp8k6ETrHZ
b2DOHYqMWyDf0alKt+gBT5kBdG0r5POzon623eCmhJu8ky0yxVnF4EYX2J0W
+LfaIS8fLXG8HoKLu/UrM5eQr7gqy9A7AXeQEbtzGTnv9Ev3TD5YG9/kujsj
m8lMuhT8CjYMW2TjXZCPf8zX43aA58I9jdfXkTMp/LHmdjBG8bSxcUeWfOfq
8PAByHpMKvTxRNbSzNC71wY0m63TB7yRrVr6qW0tgLE21P3ODznH2stlgQDd
4XPZ7v7I0i/Slld4QPKuaFh7EPIeX1N9vxFoSd9mCEKRXw15vjHiAAUBw9o6
HHm6ryAONIDnfwtNwiKQz+VppdLrALvF0jQtEjk6eViOxwbfPNAXvIlCLpbV
N4mrBlA8okMiFjmsuHd1qRyIc+2MY24juw9mcki/gKh1Ss9kEnLZ8PFM4yLA
Tr+etMVA5uaYTznlgsYgqWWPNOQ/lTVPlKcBeoVKBoOFXCCkSc4mAI/o86ZV
2cgkemp8JQ0olLxN3svDpsp9MlOC1ScWxdSKkO023vrb0uD+MNUsowTZqEKG
cpgIlacdqjrLkOW9Qzm66dA31+CkoBJZh1jNMcqFMRkcVYta5M3oYz+FFEG3
GSrt+zrkqpTg9q0SdN6w7+h95PHOfw9GyuEiV/eGFg/f+zg5MK8axlnPZzk2
I68o+Rey2PC9lLp/axv+V+T+d/Q6SL7QGTXbgfcFWvmYN8CoLz0K/+nEeSxN
EN26D+u1Wdc8enAeVv46mGqEd0u+cEroQw4y33Qu4UHlqiXPLOyEKBNyNQ82
1PWHlGGLX6buNPAgv+0Wk489wH5Sy+dBwcgCfI+9YS77cYoHZQ+7LPz6kRXX
yZXSTbAsJEzjwgD23PJSfhMk0y6ZXsXOYSp0lTZB7i1Ve29s1UcO2bVNsCdF
8EM09ug207StCc5W2FW2Yj+5VpM60QSV/pCXVoPY8zZnThKwQbhNNsC2mljb
OUVA/VcThtbYFa4FI8oEtFllunhiUxV//1mfgL4iH1j52DmNjwauEHBZbLy0
GnsT6BW5EZAq1cghsDtoucE+BExU8Rt+ih2daSkRTsBjGpZTQuyiGMbCTQLm
6cosvca2/zTYEU9AWeONzQ/YWaf20u4QsOarkaMj7Nihs15ZBPwP/RYtzQ==

       "]]},
     Annotation[#, "Charting`Private`Tag$426866#1"]& ]}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{225.0000000000005, 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  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->{{225, 230}, {-0.053789888713922096`, 0.12440886406701779`}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{3.868594718906636*^9, 3.8685948636051292`*^9, 
  3.868595801986965*^9, 3.868599678253955*^9, 3.86859979678286*^9, 
  3.8685999656321163`*^9, 3.869628052988927*^9},
 CellLabel->"Out[54]=",ExpressionUUID->"01ca38e8-4fdc-4873-bac0-259fd0f0d105"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"FindRoot", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"Integrate", "[", 
     RowBox[{
      RowBox[{
       SuperscriptBox[
        RowBox[{"PsiD", "[", 
         RowBox[{
          FractionBox["kProton", "197.327"], ",", 
          FractionBox["q", "197.327"], ",", "r"}], "]"}], "2"], "*", 
       SuperscriptBox["r", "2"]}], ",", 
      RowBox[{"{", 
       RowBox[{"r", ",", "0", ",", "Infinity"}], "}"}], ",", 
      RowBox[{"Assumptions", "\[Rule]", 
       RowBox[{
        RowBox[{"Re", "[", "q", "]"}], ">", "0"}]}]}], "]"}], "\[Equal]", 
    SuperscriptBox["4.04", "2"]}], ",", 
   RowBox[{"{", 
    RowBox[{"q", ",", "228", ",", "229"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{
  3.8685948815172367`*^9, {3.86859588791683*^9, 3.8685958896563807`*^9}, {
   3.86859605938684*^9, 3.8685960772373877`*^9}},
 CellLabel->"In[55]:=",ExpressionUUID->"cab76d12-2dad-42c6-94b6-db201a8c18ec"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"q", "\[Rule]", "228.46513877063487`"}], "}"}]], "Output",
 CellChangeTimes->{
  3.868594890181797*^9, {3.868595866815669*^9, 3.868595896125098*^9}, 
   3.868596082875332*^9, 3.869565537017137*^9, 3.86962805962097*^9},
 CellLabel->"Out[55]=",ExpressionUUID->"50db03d0-7166-4dd1-b1ad-8af23d2f754f"]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{"Plot", "[", 
   RowBox[{
    RowBox[{"PsiD", "[", 
     RowBox[{
      FractionBox["kProton", "197.327"], ",", 
      FractionBox["228.465", "197.327"], ",", "rad"}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"rad", ",", "0", ",", "20"}], "}"}]}], "]"}], "*)"}]], "Input",
 CellChangeTimes->{
  3.866791446129416*^9, 3.866791571950573*^9, {3.867204290835143*^9, 
   3.867204291028899*^9}, {3.867204355843154*^9, 3.867204356378105*^9}, 
   3.867388638320077*^9, {3.868366587067737*^9, 3.8683665895065203`*^9}, {
   3.8684236246644173`*^9, 3.8684236247974253`*^9}, {3.8685917891913433`*^9, 
   3.868591790296854*^9}, {3.868591869164423*^9, 3.868591875614442*^9}, 
   3.868592070839036*^9, 3.868592298847458*^9, {3.868592922751924*^9, 
   3.868592930081027*^9}, {3.868596211643444*^9, 3.868596220693158*^9}, {
   3.868597808816976*^9, 3.8685978218062897`*^9}, 3.8685978683675537`*^9, {
   3.868597913253585*^9, 3.868597915875341*^9}, 3.868597990681909*^9, {
   3.868598789768701*^9, 
   3.8685987918538237`*^9}},ExpressionUUID->"363c9fee-04cf-4e5c-90ea-\
4a0da118130e"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Plot", "[", 
   RowBox[{
    RowBox[{
     RowBox[{
      FractionBox[
       RowBox[{
        SqrtBox["2"], " ", 
        RowBox[{"(", 
         RowBox[{
          SuperscriptBox["\[ExponentialE]", 
           RowBox[{
            RowBox[{"-", "k"}], " ", "r"}]], "-", 
          SuperscriptBox["\[ExponentialE]", 
           RowBox[{
            RowBox[{"-", "q"}], " ", "r"}]]}], ")"}]}], 
       SqrtBox[
        RowBox[{
         FractionBox["1", "k"], "+", 
         FractionBox["1", "q"], "-", 
         FractionBox["4", 
          RowBox[{"k", "+", "q"}]]}]]], "/.", 
      RowBox[{"k", "->", " ", 
       FractionBox["kProton", "197.327"]}]}], " ", "/.", 
     RowBox[{"q", "->", 
      FractionBox["228.465", "197.327"]}]}], " ", ",", 
    RowBox[{"{", 
     RowBox[{"r", ",", "0", ",", "30"}], "}"}]}], "]"}], 
  "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{{3.868598005232416*^9, 3.86859802860509*^9}, {
   3.868598079627*^9, 3.868598114995576*^9}, {3.868598206253788*^9, 
   3.868598214091436*^9}, 3.86859833210946*^9, {3.868598489203175*^9, 
   3.868598593828712*^9}, {3.868598640088954*^9, 3.868598643934705*^9}, 
   3.868598773730584*^9},
 CellLabel->"In[56]:=",ExpressionUUID->"ec406fd6-54f4-4e5f-8e73-4494b58edd50"],

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    TagBox[
     {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[CompressedData["
1:eJwVlnk01esXxg+iUroNR9wkkSQzZSrH8yoNlKmuqUHmyjUkDTdlOmUoQ6SQ
XGVo4pqFJEkkJdXJPJyvQkXm4Zj5nd8f79rrs96117v3s553ry1p73HIiZdG
o73nnv9HDru6Qjo0VbfUWz1xeUmSbl2i6CW9JiZWREctWSfnDyL8yV606RbM
lF3eJP8RiQumaY18Tcn45NxnPlVwH6sGRG5NNeYiIVVs77rfmaBKjD3GGstR
pRmrxHZ8iXhhEUcv/gr0CwXZ8C99hb/cKasRtUo8YDJV6w6X4b3EmV1D4VXg
3ODjmPS+RgEzSrhvdw1SbhoO+E5UYK90b/twIwvBZ+Amdvc9Uk2NrXLcv6Ip
YZtD4K4P4PHJZZ3hrwNxU5F60PMBJfXeVYNq9Xjb6yTCq/oRqsHLcvrDG1GY
4v1G7uEnrO9VuNazuw12lb8ENekseK+Nmn/S0gZF9fNrflmy0Lhr/J9Tnu2o
PJ/YnhTPwq17L91+JrKRezwtf4f4V1yTeZw+cqUDTMWUkoMidZhq1Lwz1Pgd
sqcWzW7/Vo8itSW2j7d2Qpcmmisl1oAL4U1yNpc7ob3s8k/xvxowsutS2YcN
XRj85XFn/5sG9GYU/37k1I0PaZt9XsU04smS0ILjhd0ICov6fbq6Ec4ORwPo
S3/gDL9+oNZMIzpFZ9cy//uBq96fRw4ea0LbVcau46M/wTkgLJBAb0Y8W0iI
vucXSgsDLHNIM6y02Y3vY34hOpHHe+HvZtQP+Lpp7eiB5gPXcZOyZtRal8Wt
CejFc428kEW2LShT3j1YvaIfAhdPur7waYXvjTXF/rb96AhOkz+X2Aqd7s5r
mrn9IPwuT66XtuJ5/LV1D80HEDzgPsKabUUu/9s9/gmD2E8KJdeebYOnbcxK
zYFBxLFszPrC2xB3PK+fIzGEhsV89UpP22BRdLUow3MIG5+la9DYbWC5SRuL
CQ+DLKt6K6DXDtV7YXkcxjBmH26UPGvVjsh3Y6Is52G8n9tNs/Joh/Gmys6Q
omEU6Ivxxt9rx/smp0ucoyOQTeoclh5ox1aBWvaXayNwn3c+d4uHjZBtGvoZ
GSNYbnLJOZjOxt4IgT8cF0ZgML0rX0ibjTe7H6d+SR7FeyWTZ8+usKESkdUz
8W4U4UW0a36hbCQ2FSptGBxF340k2967bPzj9q7IZecY6Fs1Pzx+xoZSbE8N
X90YRAJWunT8YiPh2/Aquekx5Dtbv/kyxsZShWkL043j0NtT+cuVRqGrTPBb
gus4GiRLXqaupRDfKz+uzs9BcPoFjgYoLFHfvuOYPAcTR/kUovdTOO+n48c0
48DjxqxSkhkFU7rR0k//cvD+t9bpBnsKL23MjcffcDC0pd5vwYWC3NPj0WK9
3PyURCn2WQr8uu7ipzQmcF+8OPNVAIWzIRfsI45NQOd2un1FCIUOlu/jfOYE
NhmJaDFvUig+eVOV59MEWoulXiskUJDNjb2wZXwCT09elxdLpnBn9v4LI7FJ
cEr5N799TOFMVPae+JOT2LzUV9Q8lwK7tehGWfgkri+gYEchhQMyrz/9yJsE
T0S2XMcLCkVnqulCLdz7wxy9vWUUZF58sd5GmwJtdKLSrYJCNH9LorXMFLQM
7H0t31GgmX7v9Ds4hfaoTJVFNRRau0bcauKmcER+LPwpi8J+5ZnckdIpQJn+
PKWeQsElvknR7ikI2HmwHZoobKpYxsCyaVgWNRT1tlCIXEFnOqlO40+hZTza
7RTmrdZXhVpOo6+63M2EovB3ivTyXB8uR2TlKX2j0NSvYNaUMg165dXAxu8U
9mqpx8xXT0NOajLasItCPpPRKj00jSTxL+Eh3RQkP+7ZeGDtDOwtuiRu/aAQ
IWLs5Kkzg+Ye8fGTPynM2lmkxdrPQLkksEDwF4XT/9kMvgyZgWSu/rLLXG7g
OG/vypyBe9vTuOdc1tfzuCRYPwOBO8q8VVzOCb1YqjLDzY/w70nh8oYGPz5L
yVnUfL0yYMrlsI0h+332cblUN6uW+96US2R4itssRFol3olz2flZHKs6ehbf
Xr38qsutr27hgcjQ81nUt8y7qHHr1zN8emxtxyz6FUU2T3VSyLydk6QjMAdz
WnhsNLf/9dTzH/YKc3BkMix4ufrQRj3jqs3mcJtXnf5//boE5AxVLs7hRH9d
91muvu/WfZ+JSZhDZVa3mXsrhQyl+Iy513OYvlHYrtdM4dauQyccf85BPVJG
aKCBwgULwVUfls9D4FaivVcdBV1f73NxlvOIntRlLvtEQeqWmszClXkoZakY
SnD9IPCot9EpeR4r71duW1FNofbjkZ3b+ufh983Uyq+cwgnxnby1zAXo5Lut
Kcnj6qk6mrf9yQJmDRI3FGRx/b0n3enexwV4bV3LezOdwj8F4Y7bl9NIwgcW
neL6nbDnqngsaWQ+O/bSokgKar+aNz47TSOHTMdP8IRSkB55dunUFRoJshUM
7Ajk/l8Bd/lPSTQiG6Quu+8yhS+K7IiEPhpJ++DbXuZIoUKz+KfJAo2YaDyp
abHh+lUvhvCt5iGqbw6LtFhx54G58chpTR7y9oKQ9o2DFOx9Ss01mTxkXFCx
uliNwkhN4nqWCC8RmPyckTbNRleD9/lAOV5yxoYK0h9lo6HDolaLwUtmhbKV
q3+zUTy2wv++PS9RmWYp3G5j4+p6v07XDF5STvUbaL9kg+56Im3xbj5iU+U5
QrvMhtzilLemlnyk7vNDgWxPNkjyj+93/+Yjgg3jkQdOseHa5CamcJuPKErm
WBmYs1GhfyXMtJuPVP3UzTmlxMY58bvud4MXkShXKc3L7e34WstSk6/hJ+En
eky6VdpxS3XvC2PzJWQgJ1Go83krmuNUyoZXCpF9y+/v7mM1ITTV7xPvn6uI
sHCzTv7Leojc7rwe/e8aMmFgvWVslLsX/GWgy+IIkwrjK1fEEz/DSzvpRDVT
lNjdFD439boG+oPT3u6d68jAljQpD9FqBOiY2a3TEidbBypsx4YrsSPzbr7v
KwlyOyRw7+aJcrhODvu/E5IkZ2zv6KuHlWFKfiVb0VyKTLU0DuuFlkAn/b8H
MTGbyOcVQwGjmUUYXm1oa1EuTfyDfiztDH+GrjBbi7OLZYhqRN7FQ/dyYTB4
0aBbYwu5GrPzsK1WNvoclnjl+8kSSq5yenQsnbvfiPAGl28lunfkyn1kn2D1
R5mI7EE5In2zIfst/SHWh2SZamgoEHpUOtKUkxEnr71yzUlFIotUyVdZ9/F+
nZ76aJQSWeR1xImdcQ871pdHKjcoE3mxO+oSP2Oxgi/5nB1ThejQqyyMm25j
Q4CEjrWMKhGZ2B8XNhmFpZuPFRyqUiU+Zvahs6tuQonxZ9bF42pENWg851+J
UJxlr35jSNtGBP1lj3ieCsbgaJi5Qdw24uRh3e/jfA37lN/GGyhuJ7WMyD2f
vwagdHTnzETtdiLifFrA/aMv3BXUWr44qZPcNzKyNpQ3GAOOddZLNUjN43Mu
2xQvYo+0qJFQsga5vVI7MEL5HGT9b/xTrqlJBEV/ZP5+4Ak5J+1Xm1s0SXXN
4OarZ9zhVHrqWr6nFqld+eXV+HUXvHpQ5zgsrE0c/otx/kA5Y+zJwoRyoTb5
vu7yYHuFAxw2py+NN9pB7hru6P9N7DAsF62RO7SD2K/Im6Iqj6ModebbeMhO
wpz+t39oyBpH9DeRnRt0iNdstPLIUwuIuQU/VyrTIdf3VN2ViD4MYUHXHCNL
BomN3H6gYLUZUnf7eA9ZM0h0nu3YwSVmUPOJ2B19jEEK/+jj9MyZwngou77J
jkEM1j6a1OgxRVDD+JSDK4MoPWcx15eZYjLZb5d3AINUfKQfXe1hipad0V8f
pTPIosLiIvEUE5w6n5pgkMkg/u6zL1g3TcDJfObUl80gj45tvBp9xQSrpZom
VAsYZJeHLdGyMMGBJeLrX5YxiFvLUKaQoAlK6h45fq1nEJ3qRQrfvYxhuKJI
8UITg+wTGv/T1M4YTfuqOaKtDHIwjO9otbExxop/h9h0MIhvPlXcvtUYikkq
GT29DLJEH6rhHUYoadG7ENbPIJpWmRN/1RrBkH4YykMMYnXds2VLiRGcg89/
OTfO1eMgu3kg1ghjZUHxIpMMkqW+uOd7oBGY07EOxdMMktjFnPrmZYSV258q
HJ9jEIcD9xb12Rkh0a14fGGBQXK9Evl4TI3wP8cUPHo=
       "]]},
     Annotation[#, "Charting`Private`Tag$517545#1"]& ]}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  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->{{0, 30}, {0., 0.4985826035829891}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{{3.8685985735871162`*^9, 3.868598595359455*^9}, {
   3.8685987585892277`*^9, 3.8685987934634542`*^9}, 3.869628071502441*^9},
 CellLabel->"Out[56]=",ExpressionUUID->"2bd6fbb1-9564-433a-b163-c5529f0a682d"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "\:0424\:0443\:0440\:044c\:0435", "-", 
    "\:043f\:0440\:0435\:043e\:0431\:0440\:0430\:0437\:043e\:0432\:0430\:043d\
\:0438\:0435"}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"PsiP", "[", 
     RowBox[{"k_", ",", "q_", ",", "p_"}], "]"}], ":=", 
    RowBox[{"Simplify", "[", 
     RowBox[{
      RowBox[{
       FractionBox[
        SqrtBox[
         RowBox[{"4", " ", "\[Pi]"}]], "p"], 
       RowBox[{"Integrate", "[", 
        RowBox[{
         RowBox[{
          FractionBox[
           RowBox[{
            SqrtBox["2"], " ", 
            RowBox[{"(", 
             RowBox[{
              SuperscriptBox["\[ExponentialE]", 
               RowBox[{
                RowBox[{"-", "k"}], " ", "r"}]], "-", 
              SuperscriptBox["\[ExponentialE]", 
               RowBox[{
                RowBox[{"-", "q"}], " ", "r"}]]}], ")"}]}], 
           SqrtBox[
            RowBox[{
             FractionBox["1", "k"], "+", 
             FractionBox["1", "q"], "-", 
             FractionBox["4", 
              RowBox[{"k", "+", "q"}]]}]]], "*", 
          RowBox[{"Sin", "[", 
           RowBox[{"p", " ", "r"}], "]"}]}], ",", 
         RowBox[{"{", 
          RowBox[{"r", ",", "0", ",", "Infinity"}], "}"}]}], "]"}]}], ",", 
      RowBox[{"Assumptions", "\[Rule]", 
       RowBox[{
        RowBox[{"k", ">", "0"}], " ", "&&", " ", 
        RowBox[{"q", ">", "0"}]}]}]}], "]"}]}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"PsiP", "[", 
     RowBox[{"k_", ",", "q_", ",", "p_"}], "]"}], ":=", 
    RowBox[{"Simplify", "[", 
     RowBox[{
      RowBox[{
       SqrtBox[
        FractionBox["2", "\[Pi]"]], 
       RowBox[{"Integrate", "[", 
        RowBox[{
         RowBox[{
          FractionBox[
           RowBox[{
            SqrtBox["2"], " ", 
            RowBox[{"(", 
             RowBox[{
              SuperscriptBox["\[ExponentialE]", 
               RowBox[{
                RowBox[{"-", "k"}], " ", "r"}]], "-", 
              SuperscriptBox["\[ExponentialE]", 
               RowBox[{
                RowBox[{"-", "q"}], " ", "r"}]]}], ")"}]}], 
           SqrtBox[
            RowBox[{
             FractionBox["1", "k"], "+", 
             FractionBox["1", "q"], "-", 
             FractionBox["4", 
              RowBox[{"k", "+", "q"}]]}]]], "*", 
          RowBox[{"Sin", "[", 
           RowBox[{"p", " ", "r"}], "]"}]}], ",", 
         RowBox[{"{", 
          RowBox[{"r", ",", "0", ",", "Infinity"}], "}"}]}], "]"}]}], ",", 
      RowBox[{"Assumptions", "\[Rule]", 
       RowBox[{
        RowBox[{"k", ">", "0"}], " ", "&&", " ", 
        RowBox[{"q", ">", "0"}]}]}]}], "]"}]}], "\[IndentingNewLine]", 
   RowBox[{"Print", "[", 
    RowBox[{"PsiP", "[", 
     RowBox[{"k", ",", "q", ",", "p"}], "]"}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.842429619361777*^9, 3.842429677633604*^9}, {
   3.842429720856079*^9, 3.842429753812727*^9}, {3.842430187216261*^9, 
   3.842430258057073*^9}, {3.842430358025407*^9, 3.842430363366953*^9}, {
   3.842431304234517*^9, 3.8424313168119774`*^9}, 3.8424313659165783`*^9, {
   3.84243142345732*^9, 3.8424314893812175`*^9}, {3.842432125684045*^9, 
   3.842432157645306*^9}, {3.8424321906882353`*^9, 3.8424322290137844`*^9}, {
   3.842433503697682*^9, 3.8424335156930647`*^9}, {3.842433964243157*^9, 
   3.8424339654450226`*^9}, {3.842441765367482*^9, 3.842441789073285*^9}, {
   3.8424418280029216`*^9, 3.8424418302671328`*^9}, {3.842441902445393*^9, 
   3.842441904619342*^9}, {3.842442407989282*^9, 3.8424424240173473`*^9}, {
   3.8424446932159224`*^9, 3.84244470052919*^9}, {3.8424483586356497`*^9, 
   3.842448371032787*^9}, 3.8424486628033075`*^9, {3.8424492882931213`*^9, 
   3.842449296532782*^9}, {3.842458930056344*^9, 3.842458959414255*^9}, {
   3.8429757852823305`*^9, 3.842975785602333*^9}, {3.8429758182078094`*^9, 
   3.8429758427135096`*^9}, {3.84297588039443*^9, 3.8429758824232917`*^9}, {
   3.8429759476812525`*^9, 3.8429759822454915`*^9}, {3.8429761768185596`*^9, 
   3.8429761802928667`*^9}, {3.868592046564624*^9, 3.868592100595749*^9}, {
   3.868592455157613*^9, 3.868592462999556*^9}, {3.8685926812186604`*^9, 
   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[10]:=",ExpressionUUID->"589fa93e-41d0-46c0-b8d3-ab8df683f53d"],

Cell[BoxData[
 TemplateBox[{
   FractionBox[
    RowBox[{"2", " ", "p", " ", 
      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[{
     RowBox[{
       RowBox[{"Im", "[", "p", "]"}], "<", "k"}], "&&", 
     RowBox[{
       RowBox[{"Im", "[", "p", "]"}], "<", "q"}]}]},
  "ConditionalExpression"]], "Print",
 CellChangeTimes->{3.868604232135943*^9, 3.868612603799876*^9, 
  3.868612825817466*^9, 3.869565741857703*^9, 3.869626276514852*^9, 
  3.869628183693252*^9, 3.869812555675345*^9},
 CellLabel->
  "During evaluation of \
In[10]:=",ExpressionUUID->"337280b6-e158-4dad-b3cf-8ba39c1d5cf3"]
}, Open  ]],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{
   "\:041f\:0440\:043e\:0432\:0435\:0440\:0438\:043c", " ", 
    "\:043d\:043e\:0440\:043c\:0438\:0440\:043e\:0432\:043a\:0443"}], 
   "*)"}]}]], "Input",
 CellChangeTimes->{{3.868595166146923*^9, 
  3.8685951681358547`*^9}},ExpressionUUID->"58e80c3d-daed-4f30-890a-\
2e8e52e60c0c"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"Integrate", "[", 
    RowBox[{
     RowBox[{"4", " ", "\[Pi]", " ", 
      SuperscriptBox["p", "2"], " ", 
      SuperscriptBox[
       RowBox[{"PsiP", "[", 
        RowBox[{"k", ",", "q", ",", "p"}], "]"}], "2"], " ", 
      FractionBox["1", 
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{"2", " ", "\[Pi]"}], ")"}], "3"]]}], ",", 
     RowBox[{"{", 
      RowBox[{"p", ",", "0", ",", "Infinity"}], "}"}], ",", 
     RowBox[{"Assumptions", "\[Rule]", 
      RowBox[{
       RowBox[{"k", ">", "0"}], " ", "&&", " ", 
       RowBox[{"q", ">", "0"}]}]}]}], "]"}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{"Integrate", "[", 
   RowBox[{
    SuperscriptBox[
     RowBox[{"(", 
      FractionBox[
       RowBox[{"2", " ", "p", " ", 
        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"}]]}]]}]], ")"}], "2"], ",", 
    RowBox[{"{", 
     RowBox[{"p", ",", "0", ",", "Infinity"}], "}"}], ",", 
    RowBox[{"Assumptions", "\[Rule]", 
     RowBox[{
      RowBox[{"k", ">", "0"}], " ", "&&", " ", 
      RowBox[{"q", ">", "0"}]}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.8424297759514723`*^9, 3.8424298904774303`*^9}, {
   3.8424302132641673`*^9, 3.842430221226556*^9}, 3.8424308283157177`*^9, {
   3.842431023031686*^9, 3.842431045030918*^9}, {3.842441926294345*^9, 
   3.8424419379889317`*^9}, 3.8424421645252914`*^9, {3.8424447346828985`*^9, 
   3.8424447531082444`*^9}, {3.8424451876234922`*^9, 3.842445210127918*^9}, {
   3.842445241924023*^9, 3.842445249688614*^9}, {3.8424483961307383`*^9, 
   3.8424484004242554`*^9}, {3.8424589935407553`*^9, 
   3.8424590247120314`*^9}, {3.8429767149416723`*^9, 3.842976731081705*^9}, {
   3.842977308206172*^9, 3.842977317290962*^9}, {3.842977460771119*^9, 
   3.842977479967112*^9}, {3.8429776043458705`*^9, 3.842977607769457*^9}, {
   3.868595004463093*^9, 3.8685950206623774`*^9}, {3.868595172460483*^9, 
   3.868595185484193*^9}, {3.8685952966840687`*^9, 3.868595300860848*^9}, {
   3.868595532269759*^9, 3.8685955362275753`*^9}, {3.8686126439376163`*^9, 
   3.8686126745347157`*^9}, 3.868612839045279*^9, {3.869628388967555*^9, 
   3.8696283893478117`*^9}},
 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, 
  3.869812556427718*^9},
 CellLabel->"Out[12]=",ExpressionUUID->"6beae27f-a003-413e-be5f-41b65ba6a545"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{
    "\:0421\:0442\:0440\:043e\:0438\:043c", " ", 
     "\:0424\:0443\:0440\:044c\:0435", " ", 
     RowBox[{
     "\:043e\:0431\:0440\:0430\:0437", ".", " ", 
      "\:041a\:043e\:044d\:0444\:0444\:0438\:0446\:0438\:0435\:043d\:0442"}], 
     " ", "197.327", " ", "\:043d\:0435", " ", 
     "\:043d\:0443\:0436\:0435\:043d"}], ",", " ", 
    RowBox[{
    "\:043f\:043e\:0442\:043e\:043c\:0443", " ", "\:0447\:0442\:043e", " ", 
     "\:043c\:044b", " ", "\:043d\:0435", " ", 
     "\:043f\:0435\:0440\:0435\:0432\:043e\:0434\:0438\:043c", " ", "\:0432", 
     " ", "\:0424\:0435\:0440\:043c\:0438"}], ",", " ", 
    RowBox[{
    "\:0430", " ", "\:0440\:0430\:0431\:043e\:0442\:0430\:0435\:043c", " ", 
     "\:0432", " ", 
     "\:0438\:043c\:043f\:0443\:043b\:044c\:0441\:043d\:043e\:043c", " ", 
     "\:043f\:0440\:0435\:0434\:0441\:0442\:0430\:0432\:043b\:0435\:043d\:0438\
\:0438"}]}], "  ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{"Plot", "[", 
   RowBox[{
    RowBox[{
     RowBox[{
      SuperscriptBox[
       RowBox[{"(", 
        FractionBox[
         RowBox[{"2", " ", "p", " ", 
          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"}]]}]]}]], ")"}], "2"], "/.", 
      RowBox[{"k", "->", " ", "kProton"}]}], " ", "/.", 
     RowBox[{"q", "->", "228.465"}]}], ",", 
    RowBox[{"{", 
     RowBox[{"p", ",", "0", ",", "300"}], "}"}]}], "]"}]}]], "Input",
 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}, {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: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$94900#1"]& ]}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  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->{{0, 300}, {0., 0.011185286571355386`}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 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, 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[
 RowBox[{"Clear", "[", "x", "]"}]], "Input",
 CellChangeTimes->{{3.869629313118025*^9, 3.869629319664421*^9}},
 CellLabel->"In[68]:=",ExpressionUUID->"1e43fbb6-0494-4504-bdc4-df71c6e9d5e3"],

Cell[BoxData[
 RowBox[{"Clear", "[", "y", "]"}]], "Input",
 CellChangeTimes->{{3.869629323185244*^9, 3.869629326216897*^9}},
 CellLabel->"In[69]:=",ExpressionUUID->"b34abcdd-b1ba-4557-8728-df3df5b28b02"],

Cell[BoxData[
 RowBox[{
  RowBox[{"PsiP1", "[", 
   RowBox[{"k_", ",", "q_", ",", "px_", ",", "py_"}], "]"}], ":=", 
  FractionBox[
   RowBox[{"2", " ", 
    RowBox[{"(", 
     RowBox[{
      FractionBox["1", 
       RowBox[{
        SuperscriptBox["k", "2"], "+", 
        SuperscriptBox[
         RowBox[{"(", 
          RowBox[{"px", "+", "py"}], ")"}], "2"]}]], "-", 
      FractionBox["1", 
       RowBox[{
        SuperscriptBox[
         RowBox[{"(", 
          RowBox[{"px", "+", "py"}], ")"}], "2"], "+", 
        SuperscriptBox["q", "2"]}]]}], ")"}]}], 
   RowBox[{
    SqrtBox["\[Pi]"], " ", 
    SqrtBox[
     RowBox[{
      FractionBox["1", "k"], "+", 
      FractionBox["1", "q"], "-", 
      FractionBox["4", 
       RowBox[{"k", "+", "q"}]]}]]}]]}]], "Input",
 CellChangeTimes->{
  3.869629300173646*^9, {3.869629377151785*^9, 3.869629475904372*^9}},
 CellLabel->"In[13]:=",ExpressionUUID->"e24edc1d-a3ff-4c8b-bb21-62fa5fe16f06"],

Cell[BoxData[
 RowBox[{"Clear", "[", "k", "]"}]], "Input",
 CellChangeTimes->{{3.869630365766315*^9, 3.8696303681353207`*^9}},
 CellLabel->"In[78]:=",ExpressionUUID->"dfd70dbc-5fed-4e5c-a039-a3b3f88789e3"],

Cell[BoxData[
 RowBox[{
  RowBox[{"PsiP2", "[", 
   RowBox[{"k_", ",", "q_", ",", "px_", ",", "py_"}], "]"}], ":=", 
  RowBox[{"Simplify", "[", 
   SuperscriptBox[
    RowBox[{"(", 
     FractionBox[
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{
         FractionBox["1", 
          RowBox[{
           SuperscriptBox["k", "2"], "+", 
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{"px", "+", "py"}], ")"}], "2"]}]], "-", 
         FractionBox["1", 
          RowBox[{
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{"px", "+", "py"}], ")"}], "2"], "+", 
           SuperscriptBox["q", "2"]}]]}], ")"}]}], 
      RowBox[{
       SqrtBox["\[Pi]"], " ", 
       SqrtBox[
        RowBox[{
         FractionBox["1", "k"], "+", 
         FractionBox["1", "q"], "-", 
         FractionBox["4", 
          RowBox[{"k", "+", "q"}]]}]]}]], ")"}], "2"], "]"}]}]], "Input",
 CellChangeTimes->{{3.8696303068554907`*^9, 3.869630331365758*^9}, {
   3.869630391073823*^9, 3.869630393423407*^9}, {3.869630598960945*^9, 
   3.869630604762298*^9}, 3.869814938352353*^9},
 CellLabel->"In[21]:=",ExpressionUUID->"f53dd42a-ae11-462a-bf9f-c07d88f298c5"],

Cell[CellGroupData[{

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

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

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{
    "\:041f\:0430\:0440\:0430\:043b\:043b\:0435\:043b\:044c\:043d\:044b\:0439\
", " ", "\:0438\:043c\:043f\:0443\:043b\:044c\:0441"}], " ", "-", " ", 
    RowBox[{
    "\:0438\:043d\:0442\:0435\:0433\:0440\:0438\:0440\:043e\:0432\:0430\:043d\
\:0438\:0435", " ", "\:043f\:043e", " ", 
     SuperscriptBox["PsiP", "2"], "*", "py", "*", "d", 
     RowBox[{"(", "py", ")"}]}]}], " ", "*)"}]}]], "Input",
 CellChangeTimes->{
  3.868613231554058*^9, {3.868613430567326*^9, 3.868613521155683*^9}, {
   3.868613778472742*^9, 3.8686137787214823`*^9}, {3.869566081566175*^9, 
   3.869566081836238*^9}},ExpressionUUID->"dffe0abc-3dc5-4360-a36a-\
ddf70a8a91e3"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"PsiPar", "[", 
   RowBox[{"k_", ",", "q_", ",", "px_"}], "]"}], "=", 
  RowBox[{"Simplify", "[", 
   RowBox[{"Integrate", "[", 
    RowBox[{
     RowBox[{
      FractionBox[
       RowBox[{"4", " ", "k", " ", "q", " ", 
        SuperscriptBox[
         RowBox[{"(", 
          RowBox[{"k", "+", "q"}], ")"}], "3"]}], 
       RowBox[{"\[Pi]", " ", 
        SuperscriptBox[
         RowBox[{"(", 
          RowBox[{
           SuperscriptBox["k", "2"], "+", 
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{"px", "+", "py"}], ")"}], "2"]}], ")"}], "2"], " ", 
        SuperscriptBox[
         RowBox[{"(", 
          RowBox[{
           SuperscriptBox["px", "2"], "+", 
           RowBox[{"2", " ", "px", " ", "py"}], "+", 
           SuperscriptBox["py", "2"], "+", 
           SuperscriptBox["q", "2"]}], ")"}], "2"]}]], "*", "py"}], ",", 
     RowBox[{"{", 
      RowBox[{"py", ",", "0", ",", "Infinity"}], "}"}], ",", 
     RowBox[{"Assumptions", "\[Rule]", 
      RowBox[{
       RowBox[{"k", ">", "0"}], " ", "&&", " ", 
       RowBox[{"q", ">", "0"}]}]}]}], "]"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.869567112352491*^9, 3.8695671238156633`*^9}, 
   3.869626815028455*^9, {3.869627933495741*^9, 3.8696279413023233`*^9}, {
   3.8696284459881773`*^9, 3.869628460102339*^9}, {3.86962952577512*^9, 
   3.869629595424942*^9}, 3.869629948310508*^9, {3.869630088217606*^9, 
   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}, {3.86981492794865*^9, 3.869814957021761*^9}},
 CellLabel->"In[23]:=",ExpressionUUID->"a4ef5eb4-7778-43b2-8376-e6e4431de66f"],

Cell[BoxData[
 TemplateBox[{
   FractionBox[
    RowBox[{
      RowBox[{
        RowBox[{"-", 
          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[{
              RowBox[{"-", "5"}], " ", 
              SuperscriptBox["k", "2"]}], "+", 
            SuperscriptBox["q", "2"]}], ")"}], " ", 
        RowBox[{"ArcTan", "[", 
          FractionBox["px", "k"], "]"}]}], "+", 
      RowBox[{"2", " ", 
        SuperscriptBox["k", "3"], " ", "px", " ", 
        RowBox[{"(", 
          RowBox[{
            SuperscriptBox["k", "2"], "-", 
            RowBox[{"5", " ", 
              SuperscriptBox["q", "2"]}]}], ")"}], " ", 
        RowBox[{"ArcTan", "[", 
          FractionBox["px", "q"], "]"}]}], "-", 
      RowBox[{"4", " ", 
        SuperscriptBox["k", "3"], " ", 
        SuperscriptBox["q", "3"], " ", 
        RowBox[{"Log", "[", 
          RowBox[{
            SuperscriptBox["k", "2"], "+", 
            SuperscriptBox["px", "2"]}], "]"}]}], "+", 
      RowBox[{"4", " ", 
        SuperscriptBox["k", "3"], " ", 
        SuperscriptBox["q", "3"], " ", 
        RowBox[{"Log", "[", 
          RowBox[{
            SuperscriptBox["px", "2"], "+", 
            SuperscriptBox["q", "2"]}], "]"}]}]}], 
    RowBox[{
      SuperscriptBox["k", "2"], " ", "\[Pi]", " ", 
      SuperscriptBox[
       RowBox[{"(", 
         RowBox[{"k", "-", "q"}], ")"}], "3"], " ", 
      SuperscriptBox["q", "2"]}]], 
   InterpretationBox[
    DynamicModuleBox[{Typeset`open = False}, 
     TemplateBox[{"Expression", 
       StyleBox[
        TagBox[
         TooltipBox["\"condition\"", 
          TagBox[
           RowBox[{
             RowBox[{"(", 
               RowBox[{
                 RowBox[{"k", "\[NotEqual]", 
                   RowBox[{"Im", "[", "px", "]"}]}], "&&", 
                 RowBox[{"(", 
                   RowBox[{
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{
                    SuperscriptBox["k", "2"], "\[GreaterEqual]", 
                    SuperscriptBox[
                    RowBox[{"Im", "[", "px", "]"}], "2"]}], "&&", 
                    RowBox[{
                    SuperscriptBox["q", "2"], "\[GreaterEqual]", 
                    SuperscriptBox[
                    RowBox[{"Im", "[", "px", "]"}], "2"]}]}], ")"}], "||", 
                    RowBox[{
                    RowBox[{"Re", "[", "px", "]"}], "\[GreaterEqual]", 
                    "0"}]}], ")"}], "&&", 
                    RowBox[{"q", "\[NotEqual]", 
                    RowBox[{"Im", "[", "px", "]"}]}], "&&", 
                    RowBox[{
                    RowBox[{"k", "+", 
                    RowBox[{"Im", "[", "px", "]"}]}], "\[NotEqual]", "0"}], 
                    "&&", 
                    RowBox[{
                    RowBox[{"q", "+", 
                    RowBox[{"Im", "[", "px", "]"}]}], "\[NotEqual]", "0"}]}], 
                    ")"}], "||", 
                    RowBox[{
                    RowBox[{"Re", "[", "px", "]"}], ">", "0"}]}], ")"}]}], 
               ")"}], "||", 
             RowBox[{"(", 
               RowBox[{
                 RowBox[{
                   RowBox[{"k", "+", 
                    RowBox[{"Im", "[", "px", "]"}]}], "\[NotEqual]", "0"}], 
                 "&&", 
                 RowBox[{
                   RowBox[{"Re", "[", "px", "]"}], ">", "0"}]}], ")"}]}], 
           Short[#, 7]& ]], Annotation[#, 
          Short[
           Or[
            And[$CellContext`k != Im[$CellContext`px], 
             Or[
              And[
               Or[
                
                And[$CellContext`k^2 >= 
                 Im[$CellContext`px]^2, $CellContext`q^2 >= 
                 Im[$CellContext`px]^2], Re[$CellContext`px] >= 
                0], $CellContext`q != 
               Im[$CellContext`px], $CellContext`k + Im[$CellContext`px] != 
               0, $CellContext`q + Im[$CellContext`px] != 0], 
              Re[$CellContext`px] > 0]], 
            
            And[$CellContext`k + Im[$CellContext`px] != 0, 
             Re[$CellContext`px] > 0]], 7], "Tooltip"]& ], 
        "IconizedCustomName", StripOnInput -> False], 
       GridBox[{{
          RowBox[{
            TagBox["\"Head: \"", "IconizedLabel"], "\[InvisibleSpace]", 
            TagBox["Or", "IconizedItem"]}]}, {
          RowBox[{
            TagBox["\"Byte count: \"", "IconizedLabel"], "\[InvisibleSpace]", 
            
            TagBox["2000", "IconizedItem"]}]}}, 
        GridBoxAlignment -> {"Columns" -> {{Left}}}, DefaultBaseStyle -> 
        "Column", 
        GridBoxItemSize -> {
         "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}], 
       Dynamic[Typeset`open]}, "IconizedObject"]], 
    Or[
     And[$CellContext`k != Im[$CellContext`px], 
      Or[
       And[
        Or[
         And[$CellContext`k^2 >= Im[$CellContext`px]^2, $CellContext`q^2 >= 
          Im[$CellContext`px]^2], Re[$CellContext`px] >= 0], $CellContext`q != 
        Im[$CellContext`px], $CellContext`k + Im[$CellContext`px] != 
        0, $CellContext`q + Im[$CellContext`px] != 0], Re[$CellContext`px] > 
       0]], 
     And[$CellContext`k + Im[$CellContext`px] != 0, Re[$CellContext`px] > 0]],
     SelectWithContents -> True, Selectable -> False]},
  "ConditionalExpression"]], "Output",
 CellChangeTimes->{
  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, 3.86981497908293*^9},
 CellLabel->"Out[23]=",ExpressionUUID->"bd741071-db04-43f7-98ef-6e9e741872b6"]
}, Open  ]],

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

Cell[CellGroupData[{

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

Cell[BoxData[
 TemplateBox[{
   RowBox[{"0.027166391551979432`", "\[VeryThinSpace]", "-", 
     RowBox[{"0.0007851266342468511`", " ", "px"}], "+", 
     RowBox[{"0.000023643212121192028`", " ", "px", " ", 
       RowBox[{"ArcTan", "[", 
         RowBox[{"0.004377038058345917`", " ", "px"}], "]"}]}], "+", 
     RowBox[{"0.00047618392705276363`", " ", "px", " ", 
       RowBox[{"ArcTan", "[", 
         RowBox[{"0.02187868155499205`", " ", "px"}], "]"}]}], "+", 
     RowBox[{"0.002178093701515095`", " ", 
       RowBox[{"Log", "[", 
         RowBox[{"2089.0926720349003`", "\[VeryThinSpace]", "+", 
           SuperscriptBox["px", "2"]}], "]"}]}], "-", 
     RowBox[{"0.002178093701515095`", " ", 
       RowBox[{"Log", "[", 
         RowBox[{"52196.256225000005`", "\[VeryThinSpace]", "+", 
           SuperscriptBox["px", "2"]}], "]"}]}]}], 
   InterpretationBox[
    DynamicModuleBox[{Typeset`open = False}, 
     TemplateBox[{"Expression", 
       StyleBox[
        TagBox[
         TooltipBox["\"condition\"", 
          TagBox[
           RowBox[{
             RowBox[{
               RowBox[{"-", "45.70659331031903`"}], "<", 
               RowBox[{"Im", "[", "px", "]"}], "<", "45.70659331031903`"}], 
             "||", 
             RowBox[{"(", 
               RowBox[{
                 RowBox[{
                   RowBox[{
                    RowBox[{"228.465`", "\[VeryThinSpace]"}], "+", 
                    RowBox[{"Im", "[", "px", "]"}]}], "\[NotEqual]", "0"}], 
                 "&&", 
                 RowBox[{
                   RowBox[{
                    RowBox[{"45.70659331031903`", "\[VeryThinSpace]"}], "+", 
                    RowBox[{"Im", "[", "px", "]"}]}], "\[NotEqual]", "0"}], 
                 "&&", 
                 RowBox[{
                   RowBox[{"Im", "[", "px", "]"}], "\[NotEqual]", 
                   "45.70659331031903`"}], "&&", 
                 RowBox[{
                   RowBox[{"Im", "[", "px", "]"}], "\[NotEqual]", 
                   "228.465`"}], "&&", 
                 RowBox[{
                   RowBox[{"Re", "[", "px", "]"}], "\[GreaterEqual]", "0"}]}],
                ")"}], "||", 
             RowBox[{
               RowBox[{"Re", "[", "px", "]"}], ">", "0"}]}], Short[#, 7]& ]], 
         Annotation[#, 
          Short[
           Or[
            Inequality[-45.70659331031903, Less, 
             Im[$CellContext`px], Less, 45.70659331031903], 
            And[
            228.465 + Im[$CellContext`px] != 0, 
             45.70659331031903 + Im[$CellContext`px] != 0, 
             Im[$CellContext`px] != 45.70659331031903, Im[$CellContext`px] != 
             228.465, Re[$CellContext`px] >= 0], Re[$CellContext`px] > 0], 7],
           "Tooltip"]& ], "IconizedCustomName", StripOnInput -> False], 
       GridBox[{{
          RowBox[{
            TagBox["\"Head: \"", "IconizedLabel"], "\[InvisibleSpace]", 
            TagBox["Or", "IconizedItem"]}]}, {
          RowBox[{
            TagBox["\"Byte count: \"", "IconizedLabel"], "\[InvisibleSpace]", 
            
            TagBox["1168", "IconizedItem"]}]}}, 
        GridBoxAlignment -> {"Columns" -> {{Left}}}, DefaultBaseStyle -> 
        "Column", 
        GridBoxItemSize -> {
         "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}], 
       Dynamic[Typeset`open]}, "IconizedObject"]], 
    Or[
     Inequality[-45.70659331031903, Less, 
      Im[$CellContext`px], Less, 45.70659331031903], 
     And[228.465 + Im[$CellContext`px] != 0, 
      45.70659331031903 + Im[$CellContext`px] != 0, Im[$CellContext`px] != 
      45.70659331031903, Im[$CellContext`px] != 228.465, Re[$CellContext`px] >= 
      0], Re[$CellContext`px] > 0], SelectWithContents -> True, Selectable -> 
    False]},
  "ConditionalExpression"]], "Print",
 CellChangeTimes->{3.869630787819367*^9, 3.869815274817771*^9},
 CellLabel->
  "During evaluation of \
In[25]:=",ExpressionUUID->"ad5eddd5-5fad-4d49-aef0-c2ad504aea37"]
}, Open  ]],

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"PsiPar", "[", 
      RowBox[{"k", ",", "q", ",", "px_"}], "]"}], "/.", 
     RowBox[{"k", "->", " ", "kProton"}]}], " ", "/.", 
    RowBox[{"q", "->", "228.465"}]}], ",", 
   RowBox[{"{", 
    RowBox[{"px", ",", "0", ",", "50"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.868614314625906*^9, 3.868614370031674*^9}, 
   3.868615202683779*^9, {3.868673640124902*^9, 3.868673659152135*^9}, 
   3.868675910704361*^9, {3.868676129813981*^9, 3.86867615006883*^9}, 
   3.868677386578944*^9, {3.869567104276265*^9, 3.869567107186659*^9}, {
   3.869630793399475*^9, 3.869630807789693*^9}, {3.869630838289406*^9, 
   3.869630873887835*^9}},
 CellLabel->"In[92]:=",ExpressionUUID->"9f6413c0-e44e-4126-b089-740d915bb48d"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"paral", "[", "px", "]"}], ",", 
   RowBox[{"{", 
    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: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$124026#1"]& ]}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  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->{{0, 100}, {0., 0.020156672528549236`}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{
  3.868614371860661*^9, 3.868615484963463*^9, {3.868673641721901*^9, 
   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}, {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[
 RowBox[{"(*", " ", 
  RowBox[{
   RowBox[{
   "\:041f\:043e\:043f\:0435\:0440\:0435\:0447\:043d\:044b\:0439", " ", 
    "\:0438\:043c\:043f\:0443\:043b\:044c\:0441"}], " ", "-", " ", 
   RowBox[{
   "\:0438\:043d\:0442\:0435\:0433\:0440\:0438\:0440\:043e\:0432\:0430\:043d\
\:0438\:0435", " ", "\:043f\:043e", " ", 
    SuperscriptBox["PsiP", "2"], "*", "d", 
    RowBox[{"(", "px", ")"}]}]}], " ", "*)"}]], "Input",
 CellChangeTimes->{{3.868677479644743*^9, 
  3.868677489346261*^9}},ExpressionUUID->"cc447812-1c6e-4129-943e-\
0459fce59f1e"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"PsiPerp", "[", 
   RowBox[{"k_", ",", "q_", ",", "py_"}], "]"}], "=", 
  RowBox[{"Simplify", "[", 
   RowBox[{"Integrate", "[", 
    RowBox[{
     FractionBox[
      RowBox[{"4", " ", "k", " ", "q", " ", 
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{"k", "+", "q"}], ")"}], "3"]}], 
      RowBox[{"\[Pi]", " ", 
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{
          SuperscriptBox["k", "2"], "+", 
          SuperscriptBox[
           RowBox[{"(", 
            RowBox[{"px", "+", "py"}], ")"}], "2"]}], ")"}], "2"], " ", 
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{
          SuperscriptBox["px", "2"], "+", 
          RowBox[{"2", " ", "px", " ", "py"}], "+", 
          SuperscriptBox["py", "2"], "+", 
          SuperscriptBox["q", "2"]}], ")"}], "2"]}]], ",", 
     RowBox[{"{", 
      RowBox[{"px", ",", "0", ",", "Infinity"}], "}"}], ",", 
     RowBox[{"Assumptions", "\[Rule]", 
      RowBox[{
       RowBox[{"k", ">", "0"}], " ", "&&", " ", 
       RowBox[{"q", ">", "0"}]}]}]}], "]"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.8696310639392776`*^9, 3.8696310748684053`*^9}, 
   3.86963113563931*^9, 3.869815311381949*^9},
 CellLabel->"In[29]:=",ExpressionUUID->"b86a339e-b949-4516-8eb5-8294678723fc"],

Cell[BoxData[
 TemplateBox[{
   RowBox[{
     FractionBox["1", 
      RowBox[{
        SuperscriptBox["k", "2"], " ", "\[Pi]", " ", 
        RowBox[{"(", 
          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", "5"], " ", 
             RowBox[{"(", 
               RowBox[{
                 SuperscriptBox["py", "2"], "+", 
                 SuperscriptBox["q", "2"]}], ")"}]}], "-", 
           RowBox[{
             SuperscriptBox["k", "2"], " ", "\[Pi]", " ", 
             SuperscriptBox["q", "3"], " ", 
             RowBox[{"(", 
               RowBox[{
                 RowBox[{
                   RowBox[{"-", "5"}], " ", 
                   SuperscriptBox["py", "2"]}], "+", 
                 SuperscriptBox["q", "2"]}], ")"}], " ", 
             RowBox[{"(", 
               RowBox[{
                 SuperscriptBox["py", "2"], "+", 
                 SuperscriptBox["q", "2"]}], ")"}]}], "+", 
           RowBox[{
             SuperscriptBox["k", "7"], " ", 
             RowBox[{"(", 
               RowBox[{
                 RowBox[{
                   RowBox[{"-", "2"}], " ", "py", " ", "q"}], "+", 
                 RowBox[{"\[Pi]", " ", 
                   RowBox[{"(", 
                    RowBox[{
                    SuperscriptBox["py", "2"], "+", 
                    SuperscriptBox["q", "2"]}], ")"}]}]}], ")"}]}], "-", 
           RowBox[{
             SuperscriptBox["k", "3"], " ", "py", " ", 
             SuperscriptBox["q", "2"], " ", 
             RowBox[{"(", 
               RowBox[{
                 RowBox[{"2", " ", 
                   SuperscriptBox["q", "3"]}], "+", 
                 RowBox[{"5", " ", "\[Pi]", " ", "py", " ", 
                   RowBox[{"(", 
                    RowBox[{
                    SuperscriptBox["py", "2"], "+", 
                    SuperscriptBox["q", "2"]}], ")"}]}]}], ")"}]}], "+", 
           RowBox[{
             SuperscriptBox["k", "5"], " ", 
             RowBox[{"(", 
               RowBox[{
                 RowBox[{
                   RowBox[{"-", "2"}], " ", "py", " ", 
                   RowBox[{"(", 
                    RowBox[{"py", "-", "q"}], ")"}], " ", "q", " ", 
                   RowBox[{"(", 
                    RowBox[{"py", "+", "q"}], ")"}]}], "+", 
                 RowBox[{"\[Pi]", " ", 
                   RowBox[{"(", 
                    RowBox[{
                    SuperscriptBox["py", "2"], "-", 
                    RowBox[{"5", " ", 
                    SuperscriptBox["q", "2"]}]}], ")"}], " ", 
                   RowBox[{"(", 
                    RowBox[{
                    SuperscriptBox["py", "2"], "+", 
                    SuperscriptBox["q", "2"]}], ")"}]}]}], ")"}]}], "+", 
           RowBox[{"2", " ", 
             RowBox[{"(", 
               RowBox[{
                 SuperscriptBox["k", "2"], "+", 
                 SuperscriptBox["py", "2"]}], ")"}], " ", 
             RowBox[{"(", 
               RowBox[{
                 SuperscriptBox["py", "2"], "+", 
                 SuperscriptBox["q", "2"]}], ")"}], " ", 
             RowBox[{"(", 
               RowBox[{
                 RowBox[{
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{
                    RowBox[{"-", "5"}], " ", 
                    SuperscriptBox["k", "2"], " ", 
                    SuperscriptBox["q", "3"]}], "+", 
                    SuperscriptBox["q", "5"]}], ")"}], " ", 
                   RowBox[{"ArcTan", "[", 
                    FractionBox["py", "k"], "]"}]}], "-", 
                 RowBox[{
                   SuperscriptBox["k", "3"], " ", 
                   RowBox[{"(", 
                    RowBox[{
                    SuperscriptBox["k", "2"], "-", 
                    RowBox[{"5", " ", 
                    SuperscriptBox["q", "2"]}]}], ")"}], " ", 
                   RowBox[{"ArcTan", "[", 
                    FractionBox["py", "q"], "]"}]}]}], ")"}]}]}], ")"}]}]}], 
   InterpretationBox[
    DynamicModuleBox[{Typeset`open = False}, 
     TemplateBox[{"Expression", 
       StyleBox[
        TagBox[
         TooltipBox["\"condition\"", 
          TagBox[
           RowBox[{
             RowBox[{"(", 
               RowBox[{
                 RowBox[{"k", "\[NotEqual]", 
                   RowBox[{"Im", "[", "py", "]"}]}], "&&", 
                 RowBox[{"(", 
                   RowBox[{
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{
                    SuperscriptBox["k", "2"], "\[GreaterEqual]", 
                    SuperscriptBox[
                    RowBox[{"Im", "[", "py", "]"}], "2"]}], "&&", 
                    RowBox[{
                    SuperscriptBox["q", "2"], "\[GreaterEqual]", 
                    SuperscriptBox[
                    RowBox[{"Im", "[", "py", "]"}], "2"]}]}], ")"}], "||", 
                    RowBox[{
                    RowBox[{"Re", "[", "py", "]"}], "\[GreaterEqual]", 
                    "0"}]}], ")"}], "&&", 
                    RowBox[{"q", "\[NotEqual]", 
                    RowBox[{"Im", "[", "py", "]"}]}], "&&", 
                    RowBox[{
                    RowBox[{"k", "+", 
                    RowBox[{"Im", "[", "py", "]"}]}], "\[NotEqual]", "0"}], 
                    "&&", 
                    RowBox[{
                    RowBox[{"q", "+", 
                    RowBox[{"Im", "[", "py", "]"}]}], "\[NotEqual]", "0"}]}], 
                    ")"}], "||", 
                    RowBox[{
                    RowBox[{"Re", "[", "py", "]"}], ">", "0"}]}], ")"}]}], 
               ")"}], "||", 
             RowBox[{"(", 
               RowBox[{
                 RowBox[{
                   RowBox[{"k", "+", 
                    RowBox[{"Im", "[", "py", "]"}]}], "\[NotEqual]", "0"}], 
                 "&&", 
                 RowBox[{
                   RowBox[{"Re", "[", "py", "]"}], ">", "0"}]}], ")"}]}], 
           Short[#, 7]& ]], Annotation[#, 
          Short[
           Or[
            And[$CellContext`k != Im[$CellContext`py], 
             Or[
              And[
               Or[
                
                And[$CellContext`k^2 >= 
                 Im[$CellContext`py]^2, $CellContext`q^2 >= 
                 Im[$CellContext`py]^2], Re[$CellContext`py] >= 
                0], $CellContext`q != 
               Im[$CellContext`py], $CellContext`k + Im[$CellContext`py] != 
               0, $CellContext`q + Im[$CellContext`py] != 0], 
              Re[$CellContext`py] > 0]], 
            
            And[$CellContext`k + Im[$CellContext`py] != 0, 
             Re[$CellContext`py] > 0]], 7], "Tooltip"]& ], 
        "IconizedCustomName", StripOnInput -> False], 
       GridBox[{{
          RowBox[{
            TagBox["\"Head: \"", "IconizedLabel"], "\[InvisibleSpace]", 
            TagBox["Or", "IconizedItem"]}]}, {
          RowBox[{
            TagBox["\"Byte count: \"", "IconizedLabel"], "\[InvisibleSpace]", 
            
            TagBox["2000", "IconizedItem"]}]}}, 
        GridBoxAlignment -> {"Columns" -> {{Left}}}, DefaultBaseStyle -> 
        "Column", 
        GridBoxItemSize -> {
         "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}], 
       Dynamic[Typeset`open]}, "IconizedObject"]], 
    Or[
     And[$CellContext`k != Im[$CellContext`py], 
      Or[
       And[
        Or[
         And[$CellContext`k^2 >= Im[$CellContext`py]^2, $CellContext`q^2 >= 
          Im[$CellContext`py]^2], Re[$CellContext`py] >= 0], $CellContext`q != 
        Im[$CellContext`py], $CellContext`k + Im[$CellContext`py] != 
        0, $CellContext`q + Im[$CellContext`py] != 0], Re[$CellContext`py] > 
       0]], 
     And[$CellContext`k + Im[$CellContext`py] != 0, Re[$CellContext`py] > 0]],
     SelectWithContents -> True, Selectable -> False]},
  "ConditionalExpression"]], "Output",
 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[
 RowBox[{
  RowBox[{"perp", "[", "py_", "]"}], ":=", 
  RowBox[{"Simplify", "[", 
   RowBox[{"PsiPerp", "[", 
    RowBox[{"kProton", ",", "228.465", ",", "py"}], "]"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.869631129704208*^9, 3.869631132293288*^9}, {
  3.869631164376565*^9, 3.8696311890788593`*^9}},
 CellLabel->"In[30]:=",ExpressionUUID->"804a318d-428e-46d0-af1d-4cf0f2a8856f"],

Cell[CellGroupData[{

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

Cell[BoxData[
 TemplateBox[{
   RowBox[{
     FractionBox["1", 
      RowBox[{"1.0904281638730356`*^8", "+", 
        RowBox[{"54285.348897034906`", " ", 
          SuperscriptBox["py", "2"]}], "+", 
        SuperscriptBox["py", "4"]}]], 
     RowBox[{"(", 
       RowBox[{"85612.41941896101`", "\[VeryThinSpace]", "-", 
         RowBox[{"1365.5989463371648`", " ", "py"}], "+", 
         RowBox[{"42.62087326844503`", " ", 
           SuperscriptBox["py", "2"]}], "-", 
         RowBox[{"0.027166391551979432`", " ", 
           SuperscriptBox["py", "3"]}], "+", 
         RowBox[{"0.0007851266342468513`", " ", 
           SuperscriptBox["py", "4"]}], "+", 
         RowBox[{
           RowBox[{"(", 
             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[{"-", "51924.436524199664`"}], "-", 
               RowBox[{"25.849810619219497`", " ", 
                 SuperscriptBox["py", "2"]}], "-", 
               RowBox[{"0.0004761839270527637`", " ", 
                 SuperscriptBox["py", "4"]}]}], ")"}], " ", 
           RowBox[{"ArcTan", "[", 
             RowBox[{"0.02187868155499205`", " ", "py"}], "]"}]}]}], ")"}]}], 
   
   InterpretationBox[
    DynamicModuleBox[{Typeset`open = False}, 
     TemplateBox[{"Expression", 
       StyleBox[
        TagBox[
         TooltipBox["\"condition\"", 
          TagBox[
           RowBox[{
             RowBox[{
               RowBox[{"-", "45.70659331031903`"}], "<", 
               RowBox[{"Im", "[", "py", "]"}], "<", "45.70659331031903`"}], 
             "||", 
             RowBox[{"(", 
               RowBox[{
                 RowBox[{
                   RowBox[{
                    RowBox[{"228.465`", "\[VeryThinSpace]"}], "+", 
                    RowBox[{"Im", "[", "py", "]"}]}], "\[NotEqual]", "0"}], 
                 "&&", 
                 RowBox[{
                   RowBox[{
                    RowBox[{"45.70659331031903`", "\[VeryThinSpace]"}], "+", 
                    RowBox[{"Im", "[", "py", "]"}]}], "\[NotEqual]", "0"}], 
                 "&&", 
                 RowBox[{
                   RowBox[{"Im", "[", "py", "]"}], "\[NotEqual]", 
                   "45.70659331031903`"}], "&&", 
                 RowBox[{
                   RowBox[{"Im", "[", "py", "]"}], "\[NotEqual]", 
                   "228.465`"}], "&&", 
                 RowBox[{
                   RowBox[{"Re", "[", "py", "]"}], "\[GreaterEqual]", "0"}]}],
                ")"}], "||", 
             RowBox[{
               RowBox[{"Re", "[", "py", "]"}], ">", "0"}]}], Short[#, 7]& ]], 
         Annotation[#, 
          Short[
           Or[
            Inequality[-45.70659331031903, Less, 
             Im[$CellContext`py], Less, 45.70659331031903], 
            And[
            228.465 + Im[$CellContext`py] != 0, 
             45.70659331031903 + Im[$CellContext`py] != 0, 
             Im[$CellContext`py] != 45.70659331031903, Im[$CellContext`py] != 
             228.465, Re[$CellContext`py] >= 0], Re[$CellContext`py] > 0], 7],
           "Tooltip"]& ], "IconizedCustomName", StripOnInput -> False], 
       GridBox[{{
          RowBox[{
            TagBox["\"Head: \"", "IconizedLabel"], "\[InvisibleSpace]", 
            TagBox["Or", "IconizedItem"]}]}, {
          RowBox[{
            TagBox["\"Byte count: \"", "IconizedLabel"], "\[InvisibleSpace]", 
            
            TagBox["1168", "IconizedItem"]}]}}, 
        GridBoxAlignment -> {"Columns" -> {{Left}}}, DefaultBaseStyle -> 
        "Column", 
        GridBoxItemSize -> {
         "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}], 
       Dynamic[Typeset`open]}, "IconizedObject"]], 
    Or[
     Inequality[-45.70659331031903, Less, 
      Im[$CellContext`py], Less, 45.70659331031903], 
     And[228.465 + Im[$CellContext`py] != 0, 
      45.70659331031903 + Im[$CellContext`py] != 0, Im[$CellContext`py] != 
      45.70659331031903, Im[$CellContext`py] != 228.465, Re[$CellContext`py] >= 
      0], Re[$CellContext`py] > 0], SelectWithContents -> True, Selectable -> 
    False]},
  "ConditionalExpression"]], "Print",
 CellChangeTimes->{3.869631197009533*^9, 3.86981533848575*^9},
 CellLabel->
  "During evaluation of \
In[31]:=",ExpressionUUID->"adbc5979-c07e-4cb9-9aa5-b7ff43f4bc1a"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"perp", "[", "py", "]"}], "py"}], ",", 
   RowBox[{"{", 
    RowBox[{"py", ",", "0", ",", "300"}], "}"}]}], "]"}]], "Input",
 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: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$122611#1"]& ]}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  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->{{0, 300}, {0., 0.007608918611885073}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 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},
WindowMargins->{{0, Automatic}, {0, Automatic}},
Magnification:>1.1 Inherited,
FrontEndVersion->"13.0 for Linux x86 (64-bit) (December 2, 2021)",
StyleDefinitions->"Default.nb",
ExpressionUUID->"d3eb6006-436a-464a-aded-6b72f7177dd4"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
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[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[6436, 168, 942, 18, 59, "Input",ExpressionUUID->"ed029352-16a6-45e1-b42a-177020211247"],
Cell[CellGroupData[{
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[CellGroupData[{
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[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[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[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[CellGroupData[{
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[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[CellGroupData[{
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[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[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[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[59704, 1424, 744, 16, 57, "Input",ExpressionUUID->"dffe0abc-3dc5-4360-a36a-ddf70a8a91e3"],
Cell[CellGroupData[{
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[68711, 1651, 422, 8, 33, "Input",ExpressionUUID->"d6114d76-496a-4526-9235-2e2b26389c18"],
Cell[CellGroupData[{
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[73356, 1761, 802, 17, 33, "Input",ExpressionUUID->"9f6413c0-e44e-4126-b089-740d915bb48d"],
Cell[CellGroupData[{
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[78528, 1877, 557, 13, 33, "Input",ExpressionUUID->"cc447812-1c6e-4129-943e-0459fce59f1e"],
Cell[CellGroupData[{
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[89783, 2165, 397, 8, 33, "Input",ExpressionUUID->"804a318d-428e-46d0-af1d-4cf0f2a8856f"],
Cell[CellGroupData[{
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[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  ]]
}
]
*)