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

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

(* CreatedBy='Mathematica 12.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     36218,        839]
NotebookOptionsPosition[     34231,        798]
NotebookOutlinePosition[     34629,        814]
CellTagsIndexPosition[     34586,        811]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"Mp", "=", "938.272"}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"mass", " ", "of", " ", "proton"}], ",", " ", "MeV"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Mn", "=", "939.565"}], ";"}], 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"mass", " ", "of", " ", "neutron"}], ",", " ", "MeV"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Ebind", "=", "7.718"}], " ", ";"}], 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"binding", " ", "energy", " ", 
     SuperscriptBox[
      RowBox[{"of", " "}], "3"], "He"}], ",", " ", "MeV"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Esep", "=", "5.49351"}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"1", "p", " ", "separation", " ", "energy", " ", 
     SuperscriptBox[
      RowBox[{"for", " "}], "3"], "He"}], ",", " ", "MeV"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"p", "=", "197.327"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"mass", "=", "625.411"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"range", "=", "0.912"}], ";"}]}], "Input",
 CellChangeTimes->{{3.8762198970264225`*^9, 3.8762199052580624`*^9}, {
   3.8771827496672716`*^9, 3.8771827824193807`*^9}, {3.878286651057734*^9, 
   3.8782866641181517`*^9}, 3.87828681545533*^9, {3.878287607680107*^9, 
   3.8782876162121983`*^9}, {3.878288500983369*^9, 3.878288503593914*^9}, {
   3.878288562531385*^9, 3.878288589887306*^9}, {3.8782886404025097`*^9, 
   3.8782886526954603`*^9}, {3.879567012355792*^9, 3.879567019384864*^9}, {
   3.879568558358831*^9, 3.879568558418344*^9}, {3.87956872817546*^9, 
   3.879568728576707*^9}, {3.879571401760208*^9, 3.87957144648076*^9}, {
   3.87957152299862*^9, 3.8795715245943327`*^9}, 3.880085368016633*^9, {
   3.880086070054988*^9, 3.880086097201236*^9}, {3.8800861558228188`*^9, 
   3.880086188062545*^9}, 3.880086260893177*^9, {3.880086310836638*^9, 
   3.8800863244282427`*^9}, {3.880086371096949*^9, 3.8800865058217907`*^9}},
 CellLabel->
  "In[350]:=",ExpressionUUID->"ad6df89e-d3d8-490c-ae91-dac6d80d1b5a"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"fIn", "[", 
   RowBox[{"q_", ",", "r_", ",", "ang_"}], "]"}], ":=", 
  RowBox[{"q", " ", "r", " ", 
   RowBox[{"SphericalBesselJ", "[", 
    RowBox[{"ang", ",", 
     RowBox[{"q", " ", "r"}]}], "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fOut", "[", 
   RowBox[{"k_", ",", "r_", ",", "ang_"}], "]"}], ":=", 
  RowBox[{
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      FractionBox["2", "Pi"], "k", " ", "r"}], ")"}], 
    FractionBox["1", "2"]], " ", 
   RowBox[{"BesselK", "[", 
    RowBox[{
     RowBox[{"ang", "+", 
      FractionBox["1", "2"]}], ",", 
     RowBox[{"k", " ", "r"}]}], "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"dfInR", "[", 
   RowBox[{"q_", ",", "r_", ",", "ang_"}], "]"}], ":=", 
  RowBox[{"D", "[", 
   RowBox[{
    RowBox[{"fIn", "[", 
     RowBox[{"q", ",", "r", ",", "ang"}], "]"}], ",", "r"}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"dfIn", "[", 
   RowBox[{"q_", ",", "ang_"}], "]"}], ":=", 
  RowBox[{"Simplify", "[", 
   RowBox[{
    RowBox[{"dfInR", "[", 
     RowBox[{"q", ",", "r", ",", "ang"}], "]"}], "/.", 
    RowBox[{"r", "->", "range"}]}], "]"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"dfOutR", "[", 
   RowBox[{"k_", ",", "r_", ",", "ang_"}], "]"}], ":=", 
  RowBox[{"D", "[", 
   RowBox[{
    RowBox[{
     SuperscriptBox[
      RowBox[{"(", 
       RowBox[{
        FractionBox["2", "Pi"], "k", " ", "r"}], ")"}], 
      FractionBox["1", "2"]], " ", 
     RowBox[{"BesselK", "[", 
      RowBox[{
       RowBox[{"ang", "+", 
        FractionBox["1", "2"]}], ",", 
       RowBox[{"k", " ", "r"}]}], "]"}]}], ",", "r"}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"dfOut", "[", 
   RowBox[{"q_", ",", "ang_"}], "]"}], ":=", 
  RowBox[{"Simplify", "[", 
   RowBox[{
    RowBox[{"dfOutR", "[", 
     RowBox[{"q", ",", "r", ",", "ang"}], "]"}], "/.", 
    RowBox[{"r", "->", "range"}]}], "]"}]}]}], "Input",
 CellChangeTimes->{{3.876220044319666*^9, 3.87622004450217*^9}, {
   3.8765038952022886`*^9, 3.8765039216405926`*^9}, 3.876505041667235*^9, 
   3.8765050912134867`*^9, {3.876505234788604*^9, 3.876505237898678*^9}, {
   3.8765054584182944`*^9, 3.876505476330624*^9}, {3.876505721940834*^9, 
   3.8765057522605066`*^9}, 3.8765058613666496`*^9, 3.8769000872264614`*^9, 
   3.876912136986064*^9, 3.8771829534368153`*^9, {3.8771829917988296`*^9, 
   3.8771830129482393`*^9}, {3.8782869289935923`*^9, 3.878286953801031*^9}, {
   3.87828729637827*^9, 3.878287488606941*^9}, {3.878287533404892*^9, 
   3.8782876443143806`*^9}, {3.878288596072308*^9, 3.8782886055224047`*^9}, {
   3.878288683259207*^9, 3.8782887288588037`*^9}, {3.878288791621035*^9, 
   3.8782888301901093`*^9}},
 CellLabel->
  "In[357]:=",ExpressionUUID->"2fc4ef7c-d275-441b-9d3d-7d7239f74130"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{
    FractionBox[
     RowBox[{"fIn", "[", 
      RowBox[{
       FractionBox[
        SqrtBox[
         RowBox[{"2", " ", "mass", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "Esep"}], "+", "U"}], ")"}]}]], "p"], ",", "range", 
       ",", "0"}], "]"}], 
     RowBox[{"dfIn", "[", 
      RowBox[{
       FractionBox[
        SqrtBox[
         RowBox[{"2", " ", "mass", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "Esep"}], "+", "U"}], ")"}]}]], "p"], ",", "0"}], 
      "]"}]], "-", 
    FractionBox[
     RowBox[{"fOut", "[", 
      RowBox[{
       FractionBox[
        SqrtBox[
         RowBox[{"2", " ", "mass", " ", "Esep"}]], "p"], ",", "range", ",", 
       "0"}], "]"}], 
     RowBox[{"dfOut", "[", 
      RowBox[{
       FractionBox[
        SqrtBox[
         RowBox[{"2", " ", "mass", " ", "Esep"}]], "p"], ",", "0"}], "]"}]]}],
    ",", 
   RowBox[{"{", 
    RowBox[{"U", ",", "0", ",", "200"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.8762201382091026`*^9, 3.8762201688017263`*^9}, {
   3.8762202544209175`*^9, 3.876220291312565*^9}, {3.8762203214092555`*^9, 
   3.8762203787933655`*^9}, {3.876220573416912*^9, 3.8762206019928536`*^9}, {
   3.8765051098111005`*^9, 3.8765051099810953`*^9}, {3.8765051569651012`*^9, 
   3.8765051650170374`*^9}, {3.876505811152088*^9, 3.8765058192839603`*^9}, {
   3.8768987217644353`*^9, 3.876898731622919*^9}, {3.876901070970166*^9, 
   3.8769010836795783`*^9}, {3.876902823818077*^9, 3.876902840384226*^9}, {
   3.8771829609304247`*^9, 3.8771829704376745`*^9}, {3.87828765248904*^9, 
   3.878287675426609*^9}, {3.878288507460178*^9, 3.8782885088455267`*^9}, {
   3.878288638199546*^9, 3.878288679533744*^9}, {3.878288714552636*^9, 
   3.878288714697687*^9}, {3.878288777693362*^9, 3.87828878780177*^9}, {
   3.878288877525551*^9, 3.8782889008012238`*^9}, {3.878288973320475*^9, 
   3.878289061721037*^9}, {3.8795689049227552`*^9, 3.87956893336374*^9}, {
   3.879568964437098*^9, 3.879568981740798*^9}, {3.879569036636216*^9, 
   3.879569049028247*^9}, {3.879571489450306*^9, 3.879571489673706*^9}, {
   3.879571536098033*^9, 3.879571540253443*^9}, 3.880086196183363*^9, 
   3.8800862894777327`*^9, {3.880086334717277*^9, 3.880086336693862*^9}, {
   3.880086380215939*^9, 3.8800863803996*^9}},
 CellLabel->
  "In[363]:=",ExpressionUUID->"c2ad0be7-22ee-4300-8ad8-45a6e113c1b0"],

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    TagBox[
     {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[{{4.081632653061224*^-6, 3.2503281874312813`}, {
       4.081632653061224*^-6, 3.2503281874312813`}}], LineBox[CompressedData["
1:eJwVlHc81QsDxiXJCmfhcJQZjWuUFdEj4nQyfmfZGZWVrp3MpnI7RqJXyuiK
CmkYeTMSsnIVRcaNLpWUXBoaRnp7/3g+z5/P9/njeVT3hHB8hYWEhE7/0v99
G1OPd+CQOEbo7tfi9SXwaJWRuMhuFfg3WbxYeVIC28/saLzmqI2Xmy92EY8l
4Ks7v0p/lR7CjTq2LDEkYb/O96OcrQF6vFTNf4uTREVfwsGHJiYgarUWaocl
YUZef0vkyFa0qVAtRzhS8F/IyLviAgjcfdv/eiaFkLaJzbfst0Pu4Ijfgtcq
CNxENX8csIb2dflbqpLS2Nse56cfYIMzotmkhAfSaFpJZwe5MrH6plbN1DkZ
vCnUYQ35MOHS43W8OVsGUrC2st3PhFSR2tTFXBm4RIcYqscxsZWnVcYtksHM
21b681wmmgsUNoxXykC5M/wla5QJC8u8bEGvDGJSuiLW++/Egw1FT2MkZKFP
Opr1LpIFskcBoglZ1IoNaLfHszAnLDKew5HFdiGduqJEFgpiU8WaeLLgzTz/
x+scC9eY/r0yrrKIemS8rr+SBfr12aL7PrKoPz1T1/yJha7xyaLACFkwhb3G
ckJ34Ueu+7LP52Xh/dl8o8MBO3CHH4havZHFKrJQT1+kHezWUpok3/3K12uO
cE+ww9Euv67+97KgBtvUBaTZoY3UMBbxURYdbx1YJ27bodi10rVzURZ6L7wC
7s7aob4zUv09hQShh8eK1BLssTf1jBvZmoQbE1bMkpP2qL6vbLrBlgQ3UdEp
3TR7ZFwYe2LDIqHKSrDZ/JI93tVEmpwiSAhoyGh2brJH6Eo6c+1uEnoqC0dT
VjigLexIzXQUCQX5rYxvqQ4ovp6RWV9GQq7HxgF+lgOa+ndVqd8mIVsx82xV
vgPKip/PpFaSkJblLRp+ywFVDdkdQbUkxKYtzEz1OCC1/0uMfQcJnMP6D15S
HCH83bm57DUJIp55gY8vOILxZ8+LhTVk/FQS0fjtsiPu/TCz7FYnY35o/4vk
Uke8Ol/pfFWLjA98E+7OOkcI3dhf7aVLxrB979aWYUd4k5v6flqQ8Xe+iNiI
DAFPDyTNe5JRmWb6/CmHQO7v+f9dfpmMPTnSK2qdCNCijMWkr5JBLn6pU+BG
oHm825RRSkZok+B4yB4CqYYnPKwryNg4+/c6qXACZNNv+i3NZBS5xkXvOEtA
bInJSnpNxjmNetrdbgKdN4gemw0UWOmnb7vUSyBmMVpzSJeCT+b7Ak8NENB7
UKAVakAB4Sx1j/cPgYaMHscycwqkBB77Pk4T8NOZWtxNUJA4s1ixXpoN5wlF
/xtRFETWmhN5dmzgqneLQTsFd9zahYcINm7u2aU40kXB13niDpXPhu61h5op
TymIMdunmLKbjd3Zp1OXRig4XC94HRvMRnVf8VfDLxQkNQxEu6Sz8cLphM28
BhUXm8MKKH1sGBSlm5gIqHi+Z4HrOMjGgK9pX086FYzlJ0WTh9lY2xT2Kuw8
FflW2UHC4794Tr7Pe1RERWFLg+HHL2zI1Uk6/XufirI2yc7H8hy4HEpeHP5G
xf3OK59Ou3Mg9YS0SyiEhgqFnrAELw4+mI+Z74iiochv/kPoXg4q2p8UZibQ
IBB2nHEO4qCfC4p9Cg080+/vNeM4CCzflGFQSsNkCWu86SIHIsFLNhMTNFBP
zwzMDXGw7+NkfmiQHEQH6M7/jnDQsTK+VzJSDnMa1v2jYxx4HttiVBkvhxeN
2X3t7zjI0+VPr02TQ8l3yyf/meNAt8eYfKVcDhYB5zr16Vxs9SEiqAtyCLTd
ci/QhQs7nbCkzGx5TGvTW8s9uFCITFr/V6E8wiXmuua8uXCYPXabdkse8Y/u
Dp8O5KJ2wUXzdas80rkmiyWxXHBmtBQOzMqj2svYbDL3V7lm6WQ+XwFmlvJW
mwq4WG26ybbXRwGNat9YsVe4OPx4ycg/WAEPx6vdJG5y8Ufmxj+fJilgOMgo
dv19LmbKIxwa6hQgHG1Ys3+Mi3i+yl1VbTqSXGlNFeNckH0ux6gY0SFl9qVj
/h0Xp0oPZmyxpoO2VDUo+MSF0c7pL3XedGgnGnwvXc5DeKbfV/EcOhzTNxu/
1+QhaqxqN5mmiNxi/TtB+3k4LGZb76yqhI6hvR2XgnmgEh9Gvm5SwqxE1vPe
cB4splsnaq2VsOv3+WVb43gYrXqSVBmghHn9Bw5SKTyUOrttiqxQgmsd723Z
DR46E7WHC+0ZSJw6tTBazsNGRlhulw8Dt5VrpGnVPCQ25iWrRTEgdlTZMKGB
h4BPb9xdLjFw1/rNMftu3q89CNR/fGZA7vEhpekPPHBeXVjeeFUZ23+W6Kh9
4UGlZdNF8QZlBOsNWzrN8fBHvzMn+Zky2jJ+Xd4yPqIZ8u6mK1bjoLP4nTQy
H/Ev8+U9/Vejb/Sig64BH3i4Iy/ReA3Ofr53LPQQH5s7oxztklQQ/DTdxyOO
j896llpElgrsyvdaMo/wUSSoqYm7ogLREPFlKkl8sKwMU+xaVBA/yTvancWH
6zM+dWK5KvzG3h/WqeZDufsO8pJVYd3Y4KlYywej0e9ncJ4qVC+dtRBt4OMG
bbtm5C1VDHsYL4208rHx2Yn1jD5VsAePJ6Q+42PU9MTI4ho1bO2mx0/N8qFb
sm9srFkN9JtT7oPf+eg+fN7ca1ANX1Pum7Us8uGrOVmvMq2GcpbvQo6IE1pf
3dlxUFEda9tux+6iOmHFDxsLQZQ6SPdsY0o3O0E4RzHqtrkGpnMUXbOMnbCz
sauH4aKBrth/TY6bOWFNyKENb8I1cMok87urlRO86z2Lz5doQHTk+LamnU5o
vlzL2zCogf8Bmje6oQ==
       "]], LineBox[CompressedData["
1:eJwVj3k41PsXx2dQZMsWJUXIUqnE7XqKPqdui642/RC+s9r6jqikUhKTrYUr
kjZlWiwTJVujm+r7pWstayhxMxIjoiFl/7rz++M85znP6708Z5nXkf2+cjQa
LU42/98Go/2/Fgo8wcuFbHB1WUD2CKKn3bI84bBKiN14uw65xUIrp73CE8Sf
YzZ52WqTe9q64W2NJ7TSw/hX5bRJj0tPW1/UecLQq6SBwEYt8ug3d7m0Vk+w
4xVtnhekRQryBB6cXk+om12QHlesSU7bWSn1zMHg1vvdxk5sDVKpn7rTMg8D
pZDqmivrNUjt1AabCjUM1B00TN3UNcgVM8HsrAUYOM1a8RLC55Mepc9FvOUY
YFci/2e9W50U7XT0HdqKQaP7fNptLVWybHLR1CdHDFT0uy5W1amQtTkDifW7
MPi9sn5DfpwK+UXt8os8FwzWJCqPZSuqkNpNLdrBPhikrFbLGVJUJo96epeN
RWFAz2Gns8yVSHVBv9j+guzeuER+w1dF8lF30Oy5eAwcgr0sSnMUSUkA3141
BYPupTqvK9cpkqyIOyKjTAxqy0c9AxznkrvS32fvrMKgfICTIb2sQPb3saoT
3mKwYKvWM09Mgbxg1St514DB3pPOV/80VyBfi0ZNWW0YtLw6W1ZIypMbqjUF
x/oxWLQ4pvLkhBxpMeSUnKrKgJDr5Y9sztHJinXv8sUaDOAtqWUaudBJnxDP
huULGBDzfZ5jiDmdFNB4anlLGHCmwWaMWU8jdbVjz/+zmgFn32AREhMaqWBH
nhncx4DQyt+M4l7OEF3nbH1QCgPuUwGw0XyCMNNVTzh6kwErVsp1638ZJ/yz
JcX37jCAP+MsPndvnBhtuqWikMkAYYhrUueScULJlF5YJWLAqbIWUa7BGLG6
opa+/wMD5m7al+tm/ZMIVT4o8FnMBE7Z2ZCnj4cJIg2qUwyZsOEWrtV6aphQ
sNH/UWHChEceL11dtw4Tf2G12y1XMWFk3aLuQx1S4u4jm6FBeybgBqV2nzSk
ROUemkMIiwliH3czPZ1BQif55seL95ggHLggij4iIZRbtBuPZTDByFDJss5Q
QtD0Eiqxh0zIW9/sdrmhlxi8FVlklc+EdNWIguM2vUTF/UMJjSQT4ts/Ps+l
vhCnChy2LBIzwZ5Gd1LI/kz829j1ULiUBX1T8XZvUjuIZm3G3SRjFuR1LiSz
dDuIGtfWa6FmLDh601Bl45V2QtRWE7VrNQuU3kdq8OM+EoldhUypg8zfmvMM
i/1AbBmO0bRjsuDU5rofprebiUyNFacrU1kgvtfAWx5eQ6zu9RhXuMuCfUes
JngXq4nikoshW9JZcNdP2DSTUkVU+X098fIRC2hPL2YrFVYQ/S+Exwpfyvy3
fYK402WEFc88QNAp4wkPQ5uTi4miUlN2iAkbaBp6E1Ozccj+usunInM28A2X
FFuFXUXlAdHMkZUybo2M10hvoVa9L1igrYzXl3BHOh+gscAH7t7b2CBuTxdJ
UR7aoG+8f+9BNnBE2deiol6g0mOG28xzZHp2gkrEh2o0+fN7osETGc/t9xKH
1yDb02SHZqGMi56axJu+QcII7vHp52zIk2hFVwa+RYnxDx401rDhaLlFxuBY
HeJkmsuF9bNBicdwHvvZiGa/XWf36XDgxmXlZfKPW1Ga0V0DBxcOGIkbYxQ3
dyL1S9ZVVm4ckIbu+GHM70QRP8qCl7pzIL7kmZw80Yk4FT01sxgHEumWK+ct
FSPjgJWny7w5EGlf7V9yU4wyi0UtO4I54BNpJKqL60K5e2oT9l/hgJvakazn
WDd6GTZJwxtk+nxXNbOCXvTrSa80sokDfHkfVlxtL1rb3SS+08yBhWsT89/3
9aL0nTnkuw8cWAFKkuWGEnRJl8FHXRwgD7ttun1JgtzyXs7qjnDgquKnA2rM
PvS9m0+Va3NhixV9X7v0K1rmpDhtcoALaR1Jwj0fvyHblj+a9D24cJKvca5n
6BvazuYLNTEufEZh9wPkB5F/8ITLLIsLdfGrlHatGkRFqQOPP/rJ8qSBzs5n
B9H2b/XspJNcaIjI0YsxGEL+f914TaVwIfzgO725+76jonrLuLZmLuTtybPy
rpSiRJPSjMlWLmC7/7R+0iBFh0LcycVtXOD5NlmPtkmRseH5n8x/uZCRf2Yn
75sUXT78hf25hwuLNzo8oWkOI381ge3ALy5cD4pSH3AfRkZOOp+mF3qBojD0
flf3MIqvoNYaMb2gP8KjukEygpKDv/pzOryA3B4hLEoZRWELls439vCGH/Zp
U9vov9BDyzk1x8Xe0BNFCLIOjKG5o/RFpSwfaM7dmHz8yDjKWpsfnibxgQ8B
HYNLz0wgS6/sPTOBvnDjZ9Z9BZ9J1L7z7wLXAV+w/K3esc1sCjU9CLVCwX6g
Gplw1qNnCo0XOtsbj/iB3+aD6qdjp9Hss2NPYw8fhHSxmeST7gxyWrt9W9jP
gwCL3S/tzZxBlcw1vEPBOJh5OvrvWEYhLW+85M4JHObX3stwM6EQC7+n1hCC
AyQ5WPoup9CvY9oFtmE4RA9WKfMtKWR6YWxyJhqHd8kTmx9bU4hfQMQnXsMh
PNK6tWczhTYo7st7+jcOEcu4Nk1cCsWqXpTrK8Fh9QWTdW+8KdSoWeai/wqH
SWfVmNe+FMINbCfCy2R9m15YFPAodN164RbHNzj0dJ8IOh9EoVFM/O5jBw5/
lMe4aZ2jEHAXmal14vD4fkkpLYpCcX77T6EuHE4gYexQNIWMg/4xSO/B4aov
3bDqAoWcY4W+gUM42MyxSA5OpNDtuK5igRQHR/NLLYwrFOpL1FduGsFB6Oxy
fNtV2T+p8bnrx3Domz5epXODQm/vltN4EzjoKp8In7pJIb1Man/qFA7xdmuv
dKVSyCvn94zaGRw+vJiUr7xDody8o2OzszikJVk05wgo9B80W0QH
       "]], 
      LineBox[{{1.045293530402965, 3.257984440454694}, {3.0557610692481476`, 
       3.273178129112488}}], 
      LineBox[{{0.33199357119643075`, 3.2527423191935627`}, {
       0.9177425099948016, 3.257041469155708}}], 
      LineBox[{{0.10890748787908792`, 3.2511183198251805`}, {
       0.2044425507882675, 3.25181289071185}}]},
     Annotation[#, "Charting`Private`Tag$86638#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, 200}, {-5.332358609904164, 10.422205589087188`}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{{3.8762203703057337`*^9, 3.876220379980977*^9}, 
   3.8765050131761856`*^9, 3.8765050505154057`*^9, 3.876505166733531*^9, 
   3.876505253203033*^9, 3.8765054682425213`*^9, {3.876505736001601*^9, 
   3.876505760106862*^9}, 3.8765058280963507`*^9, 3.876505866424039*^9, 
   3.8768987397838306`*^9, 3.8769000954185047`*^9, 3.8769010854667377`*^9, 
   3.8769028434196196`*^9, 3.8769121476361966`*^9, {3.8771829651953516`*^9, 
   3.877182971722556*^9}, 3.877183021959716*^9, {3.8782890252122993`*^9, 
   3.878289063704639*^9}, 3.878292725117619*^9, 3.879568890539933*^9, 
   3.8795689360093517`*^9, {3.87956897106921*^9, 3.879568983881357*^9}, {
   3.879569032952899*^9, 3.8795690542179728`*^9}, {3.879571487079197*^9, 
   3.8795714913553343`*^9}, {3.879571532521268*^9, 3.8795715419674997`*^9}, 
   3.880084507024222*^9, 3.8800853736383677`*^9, 3.8800860259425*^9, {
   3.8800860771259003`*^9, 3.88008610210466*^9}, {3.88008616318576*^9, 
   3.880086198003504*^9}, 3.8800862293413563`*^9, {3.880086266837357*^9, 
   3.8800863385534143`*^9}, {3.880086376964778*^9, 3.880086509882822*^9}},
 CellLabel->
  "Out[363]=",ExpressionUUID->"03952c6d-0c6b-45a2-91c0-9040c75e840d"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"myU", "=", 
  RowBox[{"U", "/.", 
   RowBox[{"FindRoot", "[", 
    RowBox[{
     RowBox[{
      RowBox[{
       FractionBox[
        RowBox[{"fIn", "[", 
         RowBox[{
          FractionBox[
           SqrtBox[
            RowBox[{"2", " ", "mass", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "Esep"}], "+", "U"}], ")"}]}]], "p"], ",", 
          "range", ",", "0"}], "]"}], 
        RowBox[{"dfIn", "[", 
         RowBox[{
          FractionBox[
           SqrtBox[
            RowBox[{"2", " ", "mass", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "Esep"}], "+", "U"}], ")"}]}]], "p"], ",", "0"}], 
         "]"}]], "-", 
       FractionBox[
        RowBox[{"fOut", "[", 
         RowBox[{
          FractionBox[
           SqrtBox[
            RowBox[{"2", " ", "mass", " ", "Esep"}]], "p"], ",", "range", ",",
           "0"}], "]"}], 
        RowBox[{"dfOut", "[", 
         RowBox[{
          FractionBox[
           SqrtBox[
            RowBox[{"2", " ", "mass", " ", "Esep"}]], "p"], ",", "0"}], 
         "]"}]]}], "\[Equal]", "0"}], ",", 
     RowBox[{"{", 
      RowBox[{"U", ",", "150"}], "}"}]}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.876220514856612*^9, 3.8762205364843173`*^9}, {
   3.87641380286448*^9, 3.876413813776764*^9}, 3.8765050174731894`*^9, 
   3.876505834830249*^9, {3.8768987485820503`*^9, 3.876898765145231*^9}, {
   3.876901089974984*^9, 3.876901108003067*^9}, {3.8769028508117027`*^9, 
   3.876902864692957*^9}, {3.876912179610631*^9, 3.876912183976046*^9}, {
   3.8771829776133184`*^9, 3.8771829780132623`*^9}, {3.8782890831225157`*^9, 
   3.878289091001782*^9}, {3.878289430122672*^9, 3.878289446396935*^9}, 
   3.87957154916846*^9, {3.8800857093214073`*^9, 3.880085717549158*^9}, {
   3.8800858670665894`*^9, 3.8800858795038633`*^9}, 3.8800861705593033`*^9, 
   3.880086201872179*^9, 3.8800862943082933`*^9, {3.880086341393777*^9, 
   3.8800863417389393`*^9}, 3.8800863843913307`*^9},
 CellLabel->
  "In[364]:=",ExpressionUUID->"22f1491b-05f8-4461-9d8a-7023940b0eba"],

Cell[BoxData["124.41224361311541`"], "Output",
 CellChangeTimes->{
  3.8762205374997587`*^9, 3.8764138310515523`*^9, 3.8765050188080196`*^9, 
   3.8765050618099136`*^9, 3.8765058364421587`*^9, 3.8765058704783134`*^9, {
   3.876898758684599*^9, 3.8768987659384336`*^9}, 3.876900107897892*^9, {
   3.876901091009881*^9, 3.8769011091519423`*^9}, {3.8769028519905815`*^9, 
   3.8769028659018292`*^9}, 3.8769121525485616`*^9, 3.876912185418889*^9, {
   3.877182978736745*^9, 3.877183025647654*^9}, 3.878289092820195*^9, 
   3.878292505137838*^9, 3.8795690844906807`*^9, 3.879570642631118*^9, 
   3.879571550250393*^9, 3.880084507942923*^9, 3.88008537507251*^9, 
   3.880085724042804*^9, {3.880085881510993*^9, 3.880085919139512*^9}, 
   3.880086027259701*^9, {3.8800860790573597`*^9, 3.880086104048421*^9}, 
   3.8800861718601227`*^9, {3.880086203223221*^9, 3.88008623091987*^9}, {
   3.8800862683749933`*^9, 3.880086295540452*^9}, 3.880086343066498*^9, {
   3.880086385464319*^9, 3.8800865107303*^9}},
 CellLabel->
  "Out[364]=",ExpressionUUID->"180f3cce-0016-47de-8516-e49d3b9341e3"]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{"0.6401531831830206", " ", "2.5426173021971623", " ", 
   RowBox[{"fOut", "[", 
    RowBox[{
     FractionBox[
      SqrtBox[
       RowBox[{"2", " ", "mass", " ", "Esep"}]], "p"], ",", "r", ",", "0"}], 
    "]"}]}], "*)"}]], "Input",
 CellChangeTimes->{{3.8800859130397387`*^9, 3.88008591646206*^9}},
 CellLabel->"In[70]:=",ExpressionUUID->"0b790137-d050-4e54-820e-ce3bbd34ede5"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"myCoeff", "=", 
  RowBox[{"coeff", "/.", 
   RowBox[{"FindRoot", "[", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"fIn", "[", 
        RowBox[{
         FractionBox[
          SqrtBox[
           RowBox[{"2", " ", "mass", " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "Esep"}], "+", "myU"}], ")"}]}]], "p"], ",", 
         "range", ",", "0"}], "]"}], "-", 
       RowBox[{"coeff", " ", "*", 
        RowBox[{"fOut", "[", 
         RowBox[{
          FractionBox[
           SqrtBox[
            RowBox[{"2", " ", "mass", " ", "Esep"}]], "p"], ",", "range", ",",
           "0"}], "]"}]}]}], "\[Equal]", "0"}], ",", 
     RowBox[{"{", 
      RowBox[{"coeff", ",", "1"}], "}"}]}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.878289683612082*^9, 3.8782897289455433`*^9}, {
   3.878289801822311*^9, 3.87828982476785*^9}, 3.879569118948894*^9, {
   3.879571556146941*^9, 3.87957155754856*^9}, {3.880085925692211*^9, 
   3.8800859424236307`*^9}},
 CellLabel->
  "In[365]:=",ExpressionUUID->"444da03a-39cf-4c3d-bb1b-9d322299f630"],

Cell[BoxData["1.434098657022119`"], "Output",
 CellChangeTimes->{
  3.879570442808729*^9, 3.87957064413636*^9, 3.8795715629434137`*^9, 
   3.880084509531536*^9, 3.8800853767467327`*^9, {3.880085901406917*^9, 
   3.880085945613738*^9}, {3.88008606422783*^9, 3.880086105861314*^9}, 
   3.8800861742475863`*^9, 3.880086206862755*^9, {3.8800862701325703`*^9, 
   3.880086297423738*^9}, 3.8800863450137033`*^9, {3.8800863867709208`*^9, 
   3.880086511995571*^9}},
 CellLabel->
  "Out[365]=",ExpressionUUID->"7ab786b8-a8c7-4f5f-ab94-8d696a7c0fbe"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"myNorm", "=", 
  RowBox[{"A", "/.", 
   RowBox[{"FindRoot", "[", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{
        SuperscriptBox["A", "2"], " ", 
        RowBox[{"Integrate", "[", 
         RowBox[{
          SuperscriptBox[
           RowBox[{"fIn", "[", 
            RowBox[{
             FractionBox[
              SqrtBox[
               RowBox[{"2", " ", "mass", " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"-", "Esep"}], "+", "myU"}], ")"}]}]], "p"], ",", 
             "r", ",", "0"}], "]"}], "2"], ",", " ", 
          RowBox[{"{", 
           RowBox[{"r", ",", "0", ",", "range"}], "}"}]}], "]"}]}], "+", 
       RowBox[{
        SuperscriptBox[
         RowBox[{"(", 
          RowBox[{"myCoeff", "*", "A"}], ")"}], "2"], " ", 
        RowBox[{"Integrate", "[", 
         RowBox[{
          SuperscriptBox[
           RowBox[{"fOut", "[", 
            RowBox[{
             FractionBox[
              SqrtBox[
               RowBox[{"2", " ", "mass", " ", "Esep"}]], "p"], ",", "r", ",", 
             "0"}], "]"}], "2"], ",", " ", 
          RowBox[{"{", 
           RowBox[{"r", ",", "range", ",", "Infinity"}], "}"}]}], "]"}]}]}], 
      "\[Equal]", "1"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"A", ",", "0.5"}], "}"}]}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.878289888090169*^9, 3.87828992516182*^9}, 
   3.878290001694314*^9, {3.878290037135365*^9, 3.878290114624291*^9}, {
   3.878290154622223*^9, 3.878290162835559*^9}, {3.8782925361300364`*^9, 
   3.87829257237002*^9}, 3.879569178546549*^9, {3.879569241919948*^9, 
   3.879569245267922*^9}, {3.879569277686201*^9, 3.879569279300724*^9}, 
   3.879569335601615*^9, {3.879570630856184*^9, 3.879570632038727*^9}, {
   3.879571567743319*^9, 3.879571576536313*^9}, {3.880085952401537*^9, 
   3.8800859833274717`*^9}},
 CellLabel->
  "In[366]:=",ExpressionUUID->"7835cc6b-0406-456f-a3a3-ab6f89b509c7"],

Cell[BoxData["0.7793879594854948`"], "Output",
 CellChangeTimes->{
  3.879569148123713*^9, 3.8795692505602217`*^9, 3.879569281919847*^9, {
   3.879570637910142*^9, 3.879570646893675*^9}, 3.879571580819501*^9, 
   3.880084514207199*^9, {3.880085965594881*^9, 3.880085985349716*^9}, {
   3.880086066066136*^9, 3.880086107886919*^9}, 3.880086176698345*^9, 
   3.8800862093348007`*^9, {3.880086271816217*^9, 3.880086299605596*^9}, 
   3.8800863471623774`*^9, {3.880086388564464*^9, 3.8800865137393713`*^9}},
 CellLabel->
  "Out[366]=",ExpressionUUID->"129ca179-e6bb-47ec-97ae-7f1255d4dd2c"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"Piecewise", "[", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"myNorm", " ", 
         RowBox[{"fIn", "[", 
          RowBox[{
           FractionBox[
            SqrtBox[
             RowBox[{"2", " ", "mass", " ", 
              RowBox[{"(", 
               RowBox[{
                RowBox[{"-", "Esep"}], "+", "myU"}], ")"}]}]], "p"], ",", "r",
            ",", "0"}], "]"}]}], ",", 
        RowBox[{"r", "<", "range"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"myNorm", " ", "myCoeff", " ", 
         RowBox[{"fOut", "[", 
          RowBox[{
           FractionBox[
            SqrtBox[
             RowBox[{"2", " ", "mass", " ", "Esep"}]], "p"], ",", "r", ",", 
           "0"}], "]"}]}], ",", 
        RowBox[{"r", ">", "range"}]}], "}"}]}], "}"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"r", ",", "0", ",", "10"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.8771831750308676`*^9, 3.877183283833062*^9}, {
   3.879570725415313*^9, 3.87957082556467*^9}, {3.879571586117661*^9, 
   3.87957160585564*^9}, {3.880085993645627*^9, 3.880086011757999*^9}, 
   3.880086091143783*^9, {3.880086514826668*^9, 3.8800865176568737`*^9}},
 CellLabel->
  "In[368]:=",ExpressionUUID->"fc251146-842e-4797-bd33-46971b36a427"],

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    TagBox[
     {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[CompressedData["
1:eJwVVHk01W0XVdRbktfQpMxXZppDxf5FkZKvDK8iKkN6K6QPlYiopIgrCiWh
DJUUmhSZZbrJNZX5jj8XIa5Zvvv98ayz9jrrnOfsffY6Ss5eVm4LhYSECgTv
/zHpqvPiTa53jDzMb15jdcQaldMVLewU/8XBRSHr9Sj/xXxlUaii4mV4a6ky
x5VDkVQcqbBIMRxqIvbsMOV7OGfo6zGhkADZJRV/ViilgT5zOnZYIQuJLdvD
Pdfl4a92htqgwkdkXT3PGxYpQ0pNsne/QjWCom0lUl80Yiz+rWZ3RDVEQtLL
rkw2wtS9ltk4VQ0Thzi+4l46BoQnbAvoNZgIP3KjppMOfUPLHeFhdaAHU/eG
ijWjIWd2odrQN8SVOTc3ObUiqEZ7/6bjdEh7LM5y/96OzDfvVV/E0VGZaXxf
eVUHGuONhdfX0SFjqF8VZd8BFfcjn9cYNGExP2k+kNGBGuHrukLSzdDf4Wke
96sTqwzbpb5VtqCiykdefKIbr3LC28/p/sSeLxehzugFT4oXsf/kT2SX9HyV
nuuFmt9+I/XYnzCjC/Hoqxl4smtZCnPqJ77wH440WDAQUx1xyr6iHWxvluL7
fAb8mFEjZsc6keTY+m7uChO5psOp66M6cZBp57oylomhrEM2wqWdOCPsRZ97
wcTp85Lvi9S6sGH0bwminQn7uZiAbaNduGgn2RCqx0LystpQyb098PnZe7iN
w4KR4ZVjwU49uFCxWYwxw0KHp9a2oYs9IL2YlI8SbKxpvMOpe96DKdP3qTQD
NqjxFua3JHrROuQUQA1nI3Q9TXy+vRdPg3U0DypwoGx3levB74Wz63m9VZs4
KL6lW9whzsAl3dH3WcYczPbf9S7YzcDQ5L1/N7hx4Jv7nybfDAZE4/OWvcng
wB3fEwYuMHGbEd15RpkLodTxgvY7TJQHz5QEb+QiXkS2o+YpE8kytM4TRlx8
rT4ln9UiyAd67gw6yoW6zWyq2w4WLNOPqW+K5KL4nVK5jTUL+g4nnOISuDgi
Y8Y2OcfCnNbKttJnXNzqoqopP2YhqHS16uVCLsjTai+7FrAxf3/ptkM8LoJr
Lerr17KhMpYROTPGxRrdC78+b2FDp2LG1muei32jnzc+dGOjbsKiOUeaRFag
1Vu7ajb6xh//FWNAwrjnYospgw2PSA/rq8YkfhonTWybYeO1qLbD5gMkRJeQ
Bit0OOj2MLQccSBxJiawqCGag8iS3Oyv/iQW8lO7v2RxsF7Ku3j4GolEu69C
OaUcVItPznTfIlEjK20SOcbB9rK72dL3STgH6bsGLOeitX9e5MgjEtO9jjfO
qnKBvlkzl1QSWhmZVeZHuHC9J3aUlk2ifCmN1PfmwrY4TUw3j8Sxc6NL1W8L
6ue/+dl8IBGxychi8WcutLyegV9CQiXWxYPfxIXBoukwr0oSn8dv3WUNciHx
aof96xoSA5/oDaUKJKSson1vNJK4Lj818kafxJoW/5cyLSRkr8lLpxwm0Z45
u8/nB4l8psnW6DMkvN0tTeI6SFiY/msbFErCx143NqCbBCvzrp+ngMfozZ6t
mgwSAcvyHzi+JcEp1FV4zCKxwvPHBwsaibznpw+2cUhkN/z5sZMr4NHjUtBG
kjDdojKjKdSHS+wrJ5N5JLrizGXXyvSBYxRjpD1Awm/S03Dp5j74Xxa3DBok
Ie4Q6zS5vw+7NndFxP8ikV74MYjr0ocoOd/Ji0MkjBS7n7QE9EFjef1duWES
LSEipRVxfXg85mYZLsCebA1m/qs+LBrcu/2TAPc4/VhSUt2HwPPTwzkCbN12
S7ee1YewkOUPTwhw5WF9mx/zfQjJsNtWL+ivX8u9zF7LQwzlbAlf8P+LPQ+S
R7bxcHj83u6fgvnki0wr5g7xcLOp9YO/YP5ovXHe0nM89CuMa7YK+Am/eSax
KowHdof7oyEBfz9N2+3KqTx4ph+XKhfoQ6aJHNMt5KEyuTrGRqCfg1z+tR1t
gv6uB1RSeknQ7rtkmI7ycJlvQsvqIkFISNdbiffjSXptvGe7QOfw0t9OGv1o
WZwfMdRKIj5Ayeji8X58y7q9WrVB4OfxBpdQ/344WX8sa68lEegVHB4V148c
VJcdrBL407m7KaO2HwdW/hk5VUiiqf3udB6nHxGt6jvFBX4zszVSLF4wgO97
CPELuSS09yWdbdMbAF9GK9orncTjEgsqy2oAyzdVHF/6hITkztl3wx4DUJls
mz6eQGJcx2Hh0qcD8C/r7DC5Q6JEWibBQGIQuwz9/MmzJLZGfi3aqzWIbVt/
15xxEex38SXWYdNBfGhU31BuL/D7VMuGMwGDWLBzqPKqOYl/umMrE8lBRLtr
vVyjRoL3XGJspvgXppsXOXQxuBBVuq7WrT+MoKDztVnHBPdoF3XD1b3DuLfz
o9GgDRc3jzzWk7MaRtxxxu4tB7lYGP3BzOHsMLY4eI9VGnIx/WfAvTVpGBo7
Gn+LKXDR32Gb+W3hCNbKLYuXYnBQF6+mUVw3gtYzJq/+/MvBXYlazZQTo3B8
amakfY0N53abL8J143DZkMvyzWTCu+K6s5faNNZtPx/SM9+LN+fFrNeYzyHq
dRRdN7YHKVfO8jbaCRFJbRMZR226EJI8e3hsYAGho1Nsa062Q4jeQ3z2FyZ+
BVl7hH74AVWfU5dsNBcR+2RHn2XcbgWL/+CfZYWLieevNh9qftmM/ewzgTHe
S4iwrUE+cd/psGRFa8nLixIfuq/Mp6k2glfmHPw1exlRyDFtIj80QFeyRc7k
+HLCVtWj9HQCDaLgn+wYFScMrVTPRmrUwWzQfkLZVoL4YU16XS+sRlup6wJ7
bUnC76ntguzoKry9Ve6kyJcktqfo6+TnVeDUIt8eGk2K6PL3d1zxowwuF668
FUuUJnZH3fv2yLgUG9y1aBYnVxA5v/vOvaUX45uU9jJfg5UEw8Q1z6auCJ9e
hTzUFl5F6A24z6jaf8aVBZqH6jpXEelLJItu9H+EaIf/Ps/81cSAuEiy0ef3
YKcrxSYHrCE0rubqana+Ray6jMhHWxmiOFE2T18oH9GJp6pOKq4lyEWiz7+4
58I+1pEmxl1LRPqvz5wYf43JPc+mkovWEcExdp5l4jlgso0vlYbJEtl+ljni
CS+RarhpXf1JOeL9kT7HFNXnyL7hK+muI0+M5Z2S0RHPRJitywFhvjyR5TjX
M6eRjm0dKVlh9QpE/bG1w4VmTxGv8TC7/oAiUfc2zDf/eSqCiiPjD3xUJLLU
4/UUJVJwcJ3tsxoVJeJdNcftk1EydntWOUiGKxFHpXaGz1Y+QrKss0odX4nY
eHqPjqFVIlTKMu5bOigT655NnavQi0e0TqlmY5EyMVln7vDq8H2cGhX5Xq1J
IRQi0hMsPWLh/Oa3YngUhYh86HHekxaDcyHWs/upFELbXqzgalUMfK3zW8Xu
UYhvxMMYanEMbvF970bfpxDzyWUDRbkxeGUwNfMgiUKENpKa+g9iMFUy35L+
gkKsn5+cTTsRg+hGscjyKgpxW6ch0J5PRWKax+mb1RTiUb5MW8EvKtJ8aCb7
aimETWxtmhxJxbtV0dO1NArRIEK7w/tJRYe99Gl6M4Vo/a52I6WYCjWmjAmD
SSEam1bahkZQsTHfX/4pm0L4lQbai96kYseN9ik3LoXY8ku07F4QFRZqSa/7
eBSiMufvguwLVFw4qyg/MkIhOho3qP05QkXArmtTuaMUwvCktk68FRU3ljOa
fPgUYufkgaitFlTE56TdmZykEG17L7v6gIrUYBH3gmkKUW5R07zagIoXh92M
A2YpRNpseXLhZirylSvljP5QCMeaOyVu2lQUjapOzc9TiIyA3TskVKn4H+t7
7Vs=
       "]]},
     Annotation[#, "Charting`Private`Tag$88641#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, 10}, {0., 0.7793861446401996}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{{3.879570843622718*^9, 3.879570853774201*^9}, 
   3.879571608345409*^9, 3.880084516269706*^9, 3.880086015267064*^9, {
   3.8800860676554127`*^9, 3.8800861091655207`*^9}, 3.880086178702873*^9, 
   3.880086211051835*^9, {3.880086273266088*^9, 3.880086300605966*^9}, 
   3.88008634841398*^9, {3.880086390133037*^9, 3.8800865180731173`*^9}},
 CellLabel->
  "Out[368]=",ExpressionUUID->"8d1df92d-d2a4-4a0f-a87d-f36befdc5a77"]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{
   "R", " ", "\:0434\:043b\:044f", " ", "\:0433\:0435\:043b\:0438\:044f3"}], 
   " ", "==", " ", 
   RowBox[{"1.646", " ", "\:0444\:043c"}]}], " ", "*)"}]], "Input",
 CellChangeTimes->{{3.8800865535980377`*^9, 
  3.8800865699855547`*^9}},ExpressionUUID->"700162e1-a626-4c8d-a329-\
bbb787d8421f"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Integrate", "[", 
   RowBox[{
    RowBox[{
     SuperscriptBox["myNorm", "2"], 
     SuperscriptBox[
      RowBox[{"fIn", "[", 
       RowBox[{
        FractionBox[
         SqrtBox[
          RowBox[{"2", " ", "mass", " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "Esep"}], "+", "myU"}], ")"}]}]], "p"], ",", "r", 
        ",", "0"}], "]"}], "2"], "r"}], ",", 
    RowBox[{"{", 
     RowBox[{"r", ",", "0", ",", "range"}], "}"}]}], "]"}], "+", 
  RowBox[{"Integrate", "[", 
   RowBox[{
    RowBox[{
     SuperscriptBox["myNorm", "2"], " ", 
     SuperscriptBox["myCoeff", "2"], " ", 
     SuperscriptBox[
      RowBox[{"fOut", "[", 
       RowBox[{
        FractionBox[
         SqrtBox[
          RowBox[{"2", " ", "mass", " ", "Esep"}]], "p"], ",", "r", ",", 
        "0"}], "]"}], "2"], "r"}], ",", 
    RowBox[{"{", 
     RowBox[{"r", ",", "range", ",", "Infinity"}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.880086242737124*^9, 3.880086251815425*^9}},
 CellLabel->
  "In[370]:=",ExpressionUUID->"1a4970d6-a8d9-4293-a27c-2d4dbe2240ea"],

Cell[BoxData["1.646237974515782`"], "Output",
 CellChangeTimes->{{3.880086253659615*^9, 3.880086304157057*^9}, 
   3.880086352535798*^9, {3.8800863938671637`*^9, 3.880086521539238*^9}},
 CellLabel->
  "Out[370]=",ExpressionUUID->"4e50f98a-4679-4058-8e93-44eb471f3344"]
}, Open  ]]
},
WindowSize->{1389.75, 768.75},
WindowMargins->{{0, Automatic}, {0, Automatic}},
FrontEndVersion->"13.0 for Linux x86 (64-bit) (December 2, 2021)",
StyleDefinitions->"Default.nb",
ExpressionUUID->"ef468bfc-a077-454d-b697-d1f9ba5b95f7"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 2143, 52, 154, "Input",ExpressionUUID->"ad6df89e-d3d8-490c-ae91-dac6d80d1b5a"],
Cell[2704, 74, 2806, 75, 179, "Input",ExpressionUUID->"2fc4ef7c-d275-441b-9d3d-7d7239f74130"],
Cell[CellGroupData[{
Cell[5535, 153, 2450, 57, 75, "Input",ExpressionUUID->"c2ad0be7-22ee-4300-8ad8-45a6e113c1b0"],
Cell[7988, 212, 9178, 167, 235, "Output",ExpressionUUID->"03952c6d-0c6b-45a2-91c0-9040c75e840d"]
}, Open  ]],
Cell[CellGroupData[{
Cell[17203, 384, 2112, 53, 75, "Input",ExpressionUUID->"22f1491b-05f8-4461-9d8a-7023940b0eba"],
Cell[19318, 439, 1080, 16, 33, "Output",ExpressionUUID->"180f3cce-0016-47de-8516-e49d3b9341e3"]
}, Open  ]],
Cell[20413, 458, 420, 10, 39, "Input",ExpressionUUID->"0b790137-d050-4e54-820e-ce3bbd34ede5"],
Cell[CellGroupData[{
Cell[20858, 472, 1095, 30, 53, "Input",ExpressionUUID->"444da03a-39cf-4c3d-bb1b-9d322299f630"],
Cell[21956, 504, 541, 9, 33, "Output",ExpressionUUID->"7ab786b8-a8c7-4f5f-ab94-8d696a7c0fbe"]
}, Open  ]],
Cell[CellGroupData[{
Cell[22534, 518, 1961, 50, 53, "Input",ExpressionUUID->"7835cc6b-0406-456f-a3a3-ab6f89b509c7"],
Cell[24498, 570, 586, 9, 33, "Output",ExpressionUUID->"129ca179-e6bb-47ec-97ae-7f1255d4dd2c"]
}, Open  ]],
Cell[CellGroupData[{
Cell[25121, 584, 1347, 36, 53, "Input",ExpressionUUID->"fc251146-842e-4797-bd33-46971b36a427"],
Cell[26471, 622, 5969, 117, 237, "Output",ExpressionUUID->"8d1df92d-d2a4-4a0f-a87d-f36befdc5a77"]
}, Open  ]],
Cell[32455, 742, 348, 9, 29, "Input",ExpressionUUID->"700162e1-a626-4c8d-a329-bbb787d8421f"],
Cell[CellGroupData[{
Cell[32828, 755, 1116, 34, 53, "Input",ExpressionUUID->"1a4970d6-a8d9-4293-a27c-2d4dbe2240ea"],
Cell[33947, 791, 268, 4, 33, "Output",ExpressionUUID->"4e50f98a-4679-4058-8e93-44eb471f3344"]
}, Open  ]]
}
]
*)