(* 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[     36311,        828]
NotebookOptionsPosition[     34323,        787]
NotebookOutlinePosition[     34721,        803]
CellTagsIndexPosition[     34678,        800]
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.742"}], ";"}]}], "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}, {
   3.880087577450821*^9, 3.8800876498236017`*^9}},
 CellLabel->"In[77]:=",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[84]:=",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[90]:=",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.0993440325837347`}, {
       0.00499839233763319, 3.09936427718092}}], 
      LineBox[{{3.8501838081184587`, 3.1153700731572442`}, {4.497804249705671,
        3.118150852990969}}], LineBox[CompressedData["
1:eJwVU3s01AkDnUIiZN5ovIaZScnm0R4VuiRpPZJU855fKSE2SURbSQ9l16NZ
j80jFJLaymOFUrSlSWvLs3dIptUWwhprw+f745577rn3j3vuOdd8x17/XXNJ
JFLqLP7PxbQmc0k8Baf04xSGo/Phlrqu4dLGxSiTmNTGrNTCrm8mdW11l+NS
cgN54XYt+Fjt+sJY74BgHe7U/RQtVHQcPvDI0RHXAgYOe97WwmrKkuvqR52g
fCZo0fmihd3/yfOK+cCbphO6tYu1sbfpg/11Hzc82HvqXTyhjSThPM5UmDvi
mpVV1wu1EfjwUJBtsAcKxp+ID/Vpo1HTcNMegSeai/ta7KwWwJYcnzkQ9R2i
t8z3izy0AMSos7VvmDdOCbfJnZ4sQOH5ByxVsi+qqEu299rpoDJl1as2fz+Y
iXOez83XQVSds1+e9yYcrJD/l6Kni7vNxSNnRP6Q12RtXJ2hi5D1K+tD+Jtx
cMfCWhpbD7mltr/tCQ1ALuOaafSfejg7Wn8sImYLPnWsealMXQhl2YHo4Mit
OFQxoLnnG33c95D/kpCwDZlrd1okletjzsaX71Lj+WD3BpxdDTIUd7ZYRJ3g
ozOcychxIyPNpjWQf5oP5obcZV/dyTDTU7w3S+ODO+fbRb9vIAMtvynL8/lw
6zYx3rmZjHivsx/b7/Ah1zrCHQsig+TpOcKc4oPldd2InzqrXavnFMQKUCiw
7s7oJ0Po/qH5+WEBdp1RrSv+i4yq9Qbp5AQBoi6cpN78m4xg3zju8SQBfIM8
Bnu/kNEqcvEKyhFgbNuLkYBpMi5EN6Vb1wsQKtn1uIJJgfvVLl4tSYhwm9JN
gd4UnGaqfNoShaiLqlPfeoeCdIWi796PQhS3s9c8b6QgPzY7tjJViNyaeZsl
DyiofuVcnJ4lhOCgPT3iDwr680983XpJiCTLuUE1LylwW0y58qpJiOn/7gc+
U1Hw1XHZ/H4NEYLD6oJCHKiY/3E6r1NLhEo/krqvIxXUnKf2TboiGI4Zjq9w
omLJ1H7ZJboIgs7xgYXuVAga66pDOCJkdq9O+rSZiuoNnrsG3UUY11iS4R1F
RYQw8J7quAgmTx3qFTep0Mv/2ON0WgQNl4o1R25TcbVv38yxn0QQX56QOTZS
8SEs3kknQwS1WKuMumYqpEfzqs1KRJjer5+uekOFd9Gzsg0KER4qo82FGjQs
HvT6OUdHjDLHiMjxbTQ02bWX9+iLoVy2ctxHQsPOGOFTDl2M6gJC78oOGvJJ
Ibo3jMV4u3PIJDqcBgb1VOJ9GzFuX2bd8jtOg7pjw6HPfmIoWryGTG/Q0HvM
YeeaDDEmNw9PhOrRwWXopUScEyP61ltlKI2O0LIPNwvzxPjebSxknxEdY23Z
C9RLxNjzE1sk59Ix33JOpaJaDIOrtaacNXTYNLXM8X8uht+HpBW/7KMjTnt3
/s5FErzP/dNe+zUdd8/jUYapBEcM7R6PvqND3d5otMlCgq6TRGL/AB3JohYP
K2sJJvZNXepV0VFw1X7ws5MEuQV9XBcaAw99Sc4xUgmifTKiM3wZoP187uWZ
QgkYWSrL5c0MaHdSWyOLZ/Nyf05FKwMkZspD0WUJbMs5us4vGPicnVC1rFwC
irWyLfwvBpou7ElpbZBgj93CKSdN5uyvnd0Me2b75KXEF3gw8aa193KpiRQp
zR6Vk38w0UEVF5xlS7GipPHHW51MNG/pyozjSnHt15nsxLdMVL9oPu5tI0W/
e2SV6zATab2VkmFnKaqG3VxAM4Dbl5NkR4kUkpnnTXFSA5ToL4l9mCNFS3Hy
iZ+/GsBGKZhQL5BCeXfe60eahrh560yMW5EU57veWZKphlAEDRyovyqFzUEx
p93KEB9vl0ZW1kuRuSBxhRXfEMtCeGH53VL4Vb5uca0xRFWjpSzGQobaDNe+
HYlGcMoKeFvFk+GIKHRmQYYRHoSdkIwslWEbb+TD4wtG6GK+F4U7yDCyYlrt
7B0jqMIv8gPXyWC0fOzVtyojrDJi+2/cLQOJpunHC1+ExkjTdbwrsz4/5tgP
4SxM/jOUxrouQ1yADVV8mAWH2IbX5EoZ/E5bu2xLZqH06Paor3UyvN8/fqzw
VxbSfrp4sbVZhgmNLmXyEAtECW/uDx9l6PnVpV071hjZVhM++wZlEBbfNVMk
GaPjquJc0IgMcxQrdUpyjeFZGWK7aVKGX3oLeQ8ajGHQs2U5VYNALfNEUba2
CWY+Zcn+ohHgWvtMThSZINcmXqOPQcCiM/vk5xoTOEYEX3ljQKA7a2B6fosJ
IkcdVW0sAl+Wtp+q/scEyn9fpNRbEtiutOye8DTFn/NYd+QOBApP9SWfVpni
vFkByzmAQIMhx1h42wx6SbaKZVsJbN5C7FZ/ZIajo/f2m/AJDH5/zaO50wxE
U3/zjIjA7h768rIhM7DDlsbeCyTA9lWp7luao+Rmdef6/QRsVpNc76Sbg2Hu
eczxAAG5/MnQxovmSEx6YW0VQ0CnNL9MrcIcwZKvx7UPEfi925Tb/sQcS9Tc
7FoSCPgM8HtMddm45tuS4i8nMNcw/m5/ChsmNdJVa9MJaH5am9eZz0aq+XC/
fSaBaUcvm+EbbESMUZzp2QTeFMjHL7WzYZfN//tZIYF//SNvOCyyQHV/33pJ
OYFe7xB6YKUFuBsPjPpUErCdHBjY9NACWTXz8l1+I+CXM+UV88oCB3+0Gjet
nd27LsYtRN0Sq2z3Fr1rIFCqq9mhJrBE/Q+TpOCnBEpIav5juhyMX1cOJ7QR
CD55jvSYx8HyvraevA4CfMm6a/+4clC04UpD+/PZvYsiEnxiOEhiiOPX9BIo
6Cro+FbJwdYb9TOMEQLuQ2kbMru4SOu7PGQ7RiBMvFX79DgXjxiZ3d7jBM6H
CnZ8ZvDgdPj7uwmTBPR/73nEEvDA/s706ODc7ahYX/4Nv4+Hpem3HnSobYer
eppW/CgP/wPj1umt
       "]], LineBox[CompressedData["
1:eJwVj3k4lXkbxw8RIVt2TQnz0ttiHTU0/e5MSumlbENU9DzPOQ+SpUhx0GYY
hVfa1FjKlqRDJIXnV5hk6aqTcZDerCHbsW/FnPeP+7qvz3Uv3+93w/EAR0qc
xWJxRfX/XjUi1axgRUD1mAuR0KKN/2AyC3pFvLq/7+RMuTZWSjZ9IviFgK4y
1dKWM9pYZ5szrtwl2o+c67ae0cLo/M32uH0EvNaSP3DguyaOUF2noOtOwOPa
QOfrxhq4qWjVYeRBwFCoZ2DGtDpeZz+d5XmEABPdK5WnX6hjJqbR8qY3AZf1
FWR89qlj8blzbDlfAhJ+OjipdlIN/94uqJw5R0Bu6Df/OIEKbg2tllbhEpAf
LDDh5qpgwzWPnUyiCFDLe71OMkwF19vFDPpdFM2LVnz+S0sFy1Waq3bFE+Bb
G+fWwF6Dk9P/e6LhLgEWN5eAr6aMe6y4ZQNpBJz6JHZtcFQJm7fS4iszRf53
7tULeq2E/1aEW5BNgIRZc8Gbs0pY48JodekjAvIWCl0CehRxGmmnnVFFQFKS
saFNsgJ+sFGy/nQnAeaRhkc/58li/tpdi/7dBHgptgRpnZLF3xS4mzm9BIS5
UhlNv8hi+5mpBPcBAhzDNZVvNMvgqVe9TjuForl++PrLUjIYPGs+SbFIGCiQ
8U8/L419HVgKYuIkaMzuji50lMbXrHfAwgoSwl6yr6/Xl8ZfDEvuDUuR8Hxw
zwq9Oil8Zfo+570CCYE67oUJKlK4NeHieOp6EqZuL1usrJLE4hcY3ZQNJFSV
l+/gpUjiTSGLTlf1SKiINd/z0E8SR3qcehplQIKaYUVxkpYk1jckwkkjEoqt
wzuWwyVwwEtrya2IhE6fY0+LDqzAK6fENF8eJeFZ/GpJZR0xPMe7ldDrRQJO
2dGf/I2Fh/yNJKQJErScO/td21j4Xb/nmD2HhCxr8mD2NRZO/VhW2xFAgteX
c/l7ZVl4a7V/8Fw0CeYqd0tU5JYY1+S2BqP7JDRPZBRE7V9gbB0Cdzlli/Ix
1Yt1SguMlZxUWWiuKL9EoFFw2zyjE2N+r/KhSD8vidvpM898PZsYZldCQmhv
7JRy4hwT6W3zI6eWhBRbSTnZrzNMrnFRZFo/CSd/3un7pGOSOak6Pl82SMJc
I/dH08eTjPmCccj7IRKi5cVNp89PMria5ychJOGSWEW7jeEk0+rKc/OdI4Gc
X1YUC5tgpLmPTS1WUaAjcafOW2+cod88+tL0bwqi9+70Yx8YYbYWjnj3b6bg
0lKKGckaYaaSt3xiGVHAitxoV18yzER7PvpgZkaBYYelUuO6YebWWAG+bUUB
XZWpnjj7lalTLUilDoj+HQrSJl4MMBuP59t/9xfpLV5Am3/vZer+tORNBFBw
ujUrWt6ul2G31SsNBFFg7C+Xs02hl8k6ONTMDxHtJ0rle97qYXTQZo88LgX6
Z2o39z7qZtTWFnKcrlKg8XIT+3xfJyPeUhT9oIACxQkT4zapj0ymknV3WiEF
t8o5T1PT2xn4D//XFB4FYQvZpW+3tTPcmomVUSUUbKfuGMX7tDFzxWZXnSso
eOb286ddfAEzmlh6e7lB5H+TTvyqZ83Mx33lxS5DFGR8XyriZzQx1g5ZbiEj
Ij2DX90X2huZfOfEpZQxCjrZHi5Lqo3MuWPU/uZJEQfUJV1PqGe0QpS6Hb+J
7m/45rTFvWbcMnyUDq5mQ2DT/faY3JdMy4xm4H4jNoDJtptLbQ8Y/v1zW9Ap
EXcFZa4sf4GUWU/LLUNE7FF8KOrvCuTkOW5jcYYN0RF9flPCStSsQh/dEs4G
r0aFkhPKGAkuuyRqX2QDj6z3PAyvUAfHeHw2WcRDbbYjl/5CfZu+lPCK2ZD0
WdrWs/QtmntyaIfuBBswpeKVM9iC5LZbvF0zxQbFe1e1GzcIkE6FlpfkjIif
e9TWuAuQbU3PxYF5NsR2FVYYvRGg280hDYViHHBb5aEykd2KLKdvH7ZS4sBd
oz5919/aUYRFd5izMQfIoEe/zKd+QsvPgktjTnIg9Q+d4A+CLtSxsjj+XiAH
tl8Pc58TdqHnzuPeVcEceL5l8MgqmW50WhggPxPKgd7trmuFVt1owMCfQ0Vx
YKuq+iI/rRvxb3A0dydxoGPwjIG4dw/KPu3JFSvmgHWXPZUm6EV2xntsIqY5
kP+DbLphcj96fcTIx+8UDXzdcqcA+WGkTNAv/gyhgc2kB3frDKOjdObqd2do
+F+VgG9rNoxmgtcUm0fQwAtRvDL22zDSj51d+H6JBnHXCIecjGEUXcxcSbpB
w8RXnn2N0QiylDrIKy2nQUJCz0nDZhTFyMWJD7ygocUh/qus6yh6r/TKWauK
BkM/l31T7FFErzWfj3xFg45aTtz92FF000TD2raBhhM9evu5DaNoyqPzQ3sH
DaN2PGsDhzEE3pr/Wv2Zhi+Ts34njo2heLZjGOqiIfRjpHRuwBjSDapZm9VH
g/KyZ823xDF0KCaP8h+lwfjoDLHp3Ri6G99Vli6kQd7RdPcPnWNoIElLhj9B
w7rYiAEJ4RiKvnOl0GKWBvcnrvx2KSFqzKhl+czTYHpPjWYUhEg9Z8nxziIN
1xZDJ9PUhej4w23ZTd9p+MnX/uzZ9UJUyAucXV6mwdrHf8neQIj+AW9jz+g=

       "]], 
      LineBox[{{0.7314343869017764, 3.102323629352612}, {1.9628681609693666`, 
       3.107407951595617}, {3.722632787710295, 3.1148253544083313`}}], 
      LineBox[{{0.13254941274579646`, 3.099881775754108}, {
       0.24536209154974226`, 3.100340227512465}, {0.4907201014668314, 
       3.101339758144879}, {0.6038833664936131, 3.101801887304794}}]},
     Annotation[#, "Charting`Private`Tag$18876#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}, {-6.959527542228844, 11.975586222486642`}},
  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}, {
   3.880087582924012*^9, 3.8800876543745813`*^9}},
 CellLabel->"Out[90]=",ExpressionUUID->"d721930c-928c-45cb-8146-dfcae2431ec5"]
}, 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", ",", "170"}], "}"}]}], "]"}]}]}]], "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, 3.880087663509973*^9},
 CellLabel->"In[96]:=",ExpressionUUID->"22f1491b-05f8-4461-9d8a-7023940b0eba"],

Cell[BoxData["178.1243246255173`"], "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}, {3.8800875839389973`*^9, 
   3.880087671218701*^9}},
 CellLabel->"Out[96]=",ExpressionUUID->"8446f902-e811-45cc-a8bf-ab7c70890a1f"]
}, 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[97]:=",ExpressionUUID->"444da03a-39cf-4c3d-bb1b-9d322299f630"],

Cell[BoxData["1.344522780057045`"], "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}, {3.8800875856465683`*^9, 3.880087672538249*^9}},
 CellLabel->"Out[97]=",ExpressionUUID->"fbc9f7e4-1f37-4a69-b281-392c36da0aff"]
}, 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[98]:=",ExpressionUUID->"7835cc6b-0406-456f-a3a3-ab6f89b509c7"],

Cell[BoxData["0.8003232743677815`"], "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}, {
   3.880087588111286*^9, 3.88008767445716*^9}},
 CellLabel->"Out[98]=",ExpressionUUID->"fb1c1bd4-f3ef-407c-ae27-b9cbd3bbd016"]
}, 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[99]:=",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:eJwVU3k01QsXVZoQj8gzT9cUEc/zUi/2TwohlaFBPaREkVCGhEqITJeISEIl
ohIyJtFgnl/JzB1/0jNXV8p3vz/OOmuvtdc6+5y9j6LLORvXlTw8PEXc+n/P
DHVZo3syxiicsdHnY8dtozc9ClaHFE5jNsV9xqkoAMvvaq8pKFyECv/fn3iL
riOzLk5+tUI0embcNcwLU+Fp6Hf2m/xtFATvopXk56Hnh3vytHw+7u9tvZWY
VY61A+NqX+Qrwedq8tkg4D366l76TbIrUc04euTclkYUPExr+FxchdqXnyRt
WY2w9rV2mjCpQYRwdp3EwWbc4qu6xXJ7hVQVOa3DW9qw4U7jypmP9XCmjn01
cetCdnOWz2f5Jojv46eFfPyA+bQyjZHYJtjo7vvvT72PMHVroXVzmuCxm1Mf
H/8Rk7zf7Kt6mnH6xydBlV19MDC03h59vRVbDzeeKC34hM6nSyvVpjrQ4b1X
SOXkIC43b7bQderBlqkcGe9Xo3hUXK76OKUHsYrTKcrMUXSn7eRVae2BUuDd
TSnrx6DsdrhGYlsvNoZ7bH10eAzNvOHaPKL/QtTtjlXElzGIGw5s6Hj3AZHv
annT+Wh48jR6wFO7H5pSjvrOaxmY2DARa3G8H+/uxemmKjGg5m9hpJ7cjy/+
L5TDDRm4t0Mgm8bph57Y84BQXwaSmmJPObwdQPLDtlcTfQz40xJmzI4N4bXM
QI1dKhPPTadzVBKGcIW9X/91ERNT+fvteOuH8Ohm2I7hBibcvUXKa9WG0fns
UL36FBMOP5OC9eeGwXeuK+CbMQtZAi3XRHaPAiXdsYZ9LBgZXjp2xXEU92ba
BRJJFga9NPWnAkYxom78nrrIgkR3DLO1YBQCqqFrkqXZSEyz2hMlPIZSmxyP
3CNs6DQvKX7bNIb77ikxNu5stP8oXHQ1GUPNVY5moT8b/E6ChSb+Yzjeoa3q
mMTGNZV2oeWBMQS09zSXv2VD6VAo6+zCGK7cffSsrIuNuijtukGhcSgONjw4
NsTG0ud4nyrjcWxRqpJ+MMdGuiwsNh0dx8GqSGvTX2wY7JtSSrswjtnpkPmE
dST8nu/r9csbRwnB4ZeRISFGXy6k143jpvmqrY4qJJ5vfBZh2z8Oixrzveba
JP4LFN6qK0hDc21zjShIxBbU/XZPlYvT9AbmTUloDHqzhQga/jTR8gqxJuGG
rtuTvjQ0vQ2KDj9Ggifna9VADA0RQRXyNBcSaatkBpvv03D3jKDhpDsJHTfj
n5UvaRhxUWvL8iLR2HRKLv8DDfe299b8PE/CeXMs0qZoCBFbubj6Ionv8cXO
19fRIZvx6nxFCAl1u6Uc1+10KNo1DktHkqh7ofjGzpaOBN9/S1qjSRyWNGOY
eNLx48LXWLU4EtOXPNfoRdAhUDXiuJlKImo4UU3pLh2Vu93kB5JIKBiXm4uU
0xFDC3+rl0KiInfwNE8nHbO7bPb8lUpi/5qVMVNsOuakBe6Pp5Fgu6sVDq9g
oF9XsNMgncSVFqu2NikGzIL6Gg0ySEho+/5Xo8fAfqXfYsa5+Bk19bdCKwbe
ZviJ/XWHhPlcjU6GKwPv7s+e0OPiUfvxAzdCGbhkmXCxn8sPrFh7/mIqAyNf
/A9qcrGwtFay+zMGWPrfllS48/JDbMoONTHwXOSKdzNXz87RgA+m4wz8ysst
k+Dq7d+Z+U3/BwMN5uLdwtx9fB/US6iIMbFw5Enrc+6+/OvY28S0mLhBnCz6
lUAi54zgUV5TJvJ224TMx5LY3vZH8KwjE9MBpabJ3Ht2bzmcORbAhEg0W2w4
gsSZpJDaTioTEXln57uukli5kDPyKp+JMDv67DmuP+mHGnme1jPB85eQfEUg
iWYZUZO4eSZmKlSNzbh+u1w2OBksyIJb7A6HeG4eFsf+ifBQZUE6zJZ1mZsX
zbxH7/ccZoGxTf+Ziz2JN3ztbAMfFtos1CTsuHk75jnHp36DheVId6tpbh5j
dY2s1tSwwGoIDt20lcRkdU9nvTwboRZVVX+LkwiX48wUG7DhUX+CXCNEQuaq
nGj2ATZWXTR537uahJXpafvL19jYV6UUF7jARlHnr09/s9joNR1UH+1lw4ux
iVb6hMuz4RvSTmZj1PHTutdNJGQHVhv7x7Bh2xel3UYnkZyu87IujPt/LayL
DKkJJLQUchx92OAtfiAsfn0C1YrpJ+us2UgLVjQKcPoMVwcnIXF+Nl6LSt7e
JvwFfHznVVdfY4FfMVxtxGAaltYdVuUhTMQLt2hkO89Buv38ubkSOlwG7F7x
tn5Fwt24nYOT4/B5G+5yTm0RNH/aKym7MRR7r7eV2PMT1RdWVk59HkH2JY8J
nUM8hPt0vCL5xxDCspYOzE+uIK5Zm2z5OtQPnp5RoiaIl5i3C9sj1NAH1Qun
Au00VhP8H8/gV/4H0BdSDwq8XEMIcQr7fgz1woJxJiTJZx1xUHOt83mxHljT
qZpycvyE2+LxfZ4XuzDR4HKlsUiAiOy2dLaV7IS2yAdZEydBYjJ6r66rXxv4
sXB8cE6IqHaL92yRaYHZF4dvSvbCRJ2wbszI+0b01Z9c4bBZhDCycVRpfvoO
ZVFvHBUWRAjb81Tl2/Q3OLXab7S9fQNxMy6irkC6ASd8L5WtTxclBuzuaDsk
v8YWN812q+NiRK+N+DYFvTp0bNgs4LdtIyHVIaDt8fAlqp+EZWzmFSfaaR90
UiyrcWmFxv7WIXEiZ+xq56b5CvAPBpl7lf5OnLLVXqz+9AKMh4rJWcESxD8J
DUHX1pchWV1yVaW9JLEj0yX2iHEJqOmn3h9XkCJaGm/kSlUUwyH5n/b1LCmi
zdqXI2L/DN93PeBk1UoTJpX9/zkxi0Bj7Aysvy5DWGgcTlaOeowcQ13ptuOy
xM66gmMiWvkoivATcdOSI7yHIxPvqOfhuv0JS94FOWIsY4An6OAD6A9m519v
kyfKzmoZ/xmZi7RNGUVtlgqEqFmIH/V7Ni7XxaVZVioQNycjfCzO3sNeafsH
zcqKhMaJFIV4sbsw9np/VCRakZiRFZM/VpqBLBkX5dYFRSK4d3iF4KHbUG7I
u2V9VImI6LKczdufCqpWvUZ3rRLh1MkpfxqcglNzq7qaNCjErZJdZx0LbsKl
eFYhOoFCPDlaHboglQTPMNsli0QKUVZotmgmlgQ/29KP629SiO8LRybvCCYh
asEvnnqLQrTEeulY8yThyTbOj9RMChFr8VOunZEIzuvlDw8fU4jK6vmWncWJ
oHavj3vznkK4HHypztmTiPTcs+6RTRSiePtYRZpJInIvtJuYt1AILwnZ3O2G
iXghTl1saacQbk0nwyJ0EjHoIOre8y+Xr/rLduvviVCjSZqM0yjElFRZGD+D
Cp3SILn7DAohf3rWqmOYiu0RAxxXFoWwLI49cauPCiu1zGfkBIVoDrzgq9lK
ha+HgtzMDIWwc3Zf9i6hInjHVc7zOQoR2TguZl5ERYTgeO+FBQqR/3HxnGIe
FWlPc2O+f6cQ+9QK5vrTqci5ssqtapFC9FmmSFUnU/H4gOvO4CUKITny9lpm
PBWlSu9kjX5RiK1pp7TDoqionVPlLC9TCNGQHMnTYVT8D5IsW4M=
       "]]},
     Annotation[#, "Charting`Private`Tag$22258#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.8003187560651602}},
  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}, {
   3.880087591163252*^9, 3.880087675790471*^9}},
 CellLabel->"Out[99]=",ExpressionUUID->"a19709bd-12a5-47c1-be49-b078c7a3cbb6"]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{
   "R", " ", "\:0434\:043b\:044f", " ", "\:0433\:0435\:043b\:0438\:044f3"}], 
   " ", "==", " ", 
   RowBox[{"1.561", " ", "\:0444\:043c"}]}], " ", "*)"}]], "Input",
 CellChangeTimes->{{3.8800865535980377`*^9, 3.8800865699855547`*^9}, {
  3.880087570285039*^9, 
  3.880087570627646*^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[100]:=",ExpressionUUID->"1a4970d6-a8d9-4293-a27c-2d4dbe2240ea"],

Cell[BoxData["1.5612379745157818`"], "Output",
 CellChangeTimes->{{3.880086253659615*^9, 3.880086304157057*^9}, 
   3.880086352535798*^9, {3.8800863938671637`*^9, 3.880086521539238*^9}, {
   3.880087593758895*^9, 3.880087646136713*^9}, 3.8800876779743423`*^9},
 CellLabel->
  "Out[100]=",ExpressionUUID->"f2064ae6-052c-4c29-85b6-4122f8ffe9c5"]
}, 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, 2191, 52, 154, "Input",ExpressionUUID->"ad6df89e-d3d8-490c-ae91-dac6d80d1b5a"],
Cell[2752, 74, 2802, 74, 179, "Input",ExpressionUUID->"2fc4ef7c-d275-441b-9d3d-7d7239f74130"],
Cell[CellGroupData[{
Cell[5579, 152, 2446, 56, 75, "Input",ExpressionUUID->"c2ad0be7-22ee-4300-8ad8-45a6e113c1b0"],
Cell[8028, 210, 9175, 166, 235, "Output",ExpressionUUID->"d721930c-928c-45cb-8146-dfcae2431ec5"]
}, Open  ]],
Cell[CellGroupData[{
Cell[17240, 381, 2130, 52, 75, "Input",ExpressionUUID->"22f1491b-05f8-4461-9d8a-7023940b0eba"],
Cell[19373, 435, 1127, 16, 33, "Output",ExpressionUUID->"8446f902-e811-45cc-a8bf-ab7c70890a1f"]
}, Open  ]],
Cell[20515, 454, 420, 10, 39, "Input",ExpressionUUID->"0b790137-d050-4e54-820e-ce3bbd34ede5"],
Cell[CellGroupData[{
Cell[20960, 468, 1091, 29, 53, "Input",ExpressionUUID->"444da03a-39cf-4c3d-bb1b-9d322299f630"],
Cell[22054, 499, 585, 8, 33, "Output",ExpressionUUID->"fbc9f7e4-1f37-4a69-b281-392c36da0aff"]
}, Open  ]],
Cell[CellGroupData[{
Cell[22676, 512, 1957, 49, 53, "Input",ExpressionUUID->"7835cc6b-0406-456f-a3a3-ab6f89b509c7"],
Cell[24636, 563, 631, 9, 33, "Output",ExpressionUUID->"fb1c1bd4-f3ef-407c-ae27-b9cbd3bbd016"]
}, Open  ]],
Cell[CellGroupData[{
Cell[25304, 577, 1343, 35, 53, "Input",ExpressionUUID->"fc251146-842e-4797-bd33-46971b36a427"],
Cell[26650, 614, 5758, 112, 237, "Output",ExpressionUUID->"a19709bd-12a5-47c1-be49-b078c7a3cbb6"]
}, Open  ]],
Cell[32423, 729, 397, 10, 29, "Input",ExpressionUUID->"700162e1-a626-4c8d-a329-bbb787d8421f"],
Cell[CellGroupData[{
Cell[32845, 743, 1116, 34, 53, "Input",ExpressionUUID->"1a4970d6-a8d9-4293-a27c-2d4dbe2240ea"],
Cell[33964, 779, 343, 5, 33, "Output",ExpressionUUID->"f2064ae6-052c-4c29-85b6-4122f8ffe9c5"]
}, Open  ]]
}
]
*)