Skip to content

M
mathematica
Switch branch/tag
Find file
Normal viewHistoryPermalink
bounded3He.nb 36.4 KB
Newer Older
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
1 2 3 4 5 6 7 8 9 10 11 12 
(* 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]
himyss's avatar
radii normed
9c694860
 
himyss committed
13 14 15 16 
NotebookDataLength[     36311,        828]
NotebookOptionsPosition[     34323,        787]
NotebookOutlinePosition[     34721,        803]
CellTagsIndexPosition[     34678,        800]
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 
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[{
himyss's avatar
3He
a9d10661
 
himyss committed
46 47 
  RowBox[{
   RowBox[{"Esep", "=", "5.49351"}], ";"}], " ",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
48 49 50 51 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"1", "p", " ", "separation", " ", "energy", " ", 
     SuperscriptBox[
himyss's avatar
3He
a9d10661
 
himyss committed
52 53 54 
      RowBox[{"for", " "}], "3"], "He"}], ",", " ", "MeV"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
55 56 57 58 
  RowBox[{"p", "=", "197.327"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"mass", "=", "625.411"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
himyss's avatar
radii normed
9c694860
 
himyss committed
59 
  RowBox[{"range", "=", "0.742"}], ";"}]}], "Input",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
60 61 62 63 64 65 66 67 
 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}, {
himyss's avatar
3He
a9d10661
 
himyss committed
68 69 70 
   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,
himyss's avatar
radii normed
9c694860
 
himyss committed
71 72 73 
   3.8800863244282427`*^9}, {3.880086371096949*^9, 3.8800865058217907`*^9}, {
   3.880087577450821*^9, 3.8800876498236017`*^9}},
 CellLabel->"In[77]:=",ExpressionUUID->"ad6df89e-d3d8-490c-ae91-dac6d80d1b5a"],
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 

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}},
himyss's avatar
radii normed
9c694860
 
himyss committed
149 
 CellLabel->"In[84]:=",ExpressionUUID->"2fc4ef7c-d275-441b-9d3d-7d7239f74130"],
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 

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[{"{",
himyss's avatar
3He
a9d10661
 
himyss committed
190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 
    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}},
himyss's avatar
radii normed
9c694860
 
himyss committed
209 
 CellLabel->"In[90]:=",ExpressionUUID->"c2ad0be7-22ee-4300-8ad8-45a6e113c1b0"],
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
210 211 212 213 214 

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    TagBox[
     {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
himyss's avatar
radii normed
9c694860
 
himyss committed
215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 
      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
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
269 
       "]], LineBox[CompressedData["
himyss's avatar
radii normed
9c694860
 
himyss committed
270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 
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=
himyss's avatar
3He
a9d10661
 
himyss committed
315 
       "]],
himyss's avatar
radii normed
9c694860
 
himyss committed
316 317 318 319 320 321 
      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"]& ], {}}, {}},
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 
  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]]}& )}},
himyss's avatar
radii normed
9c694860
 
himyss committed
353 
  PlotRange->{{0, 200}, {-6.959527542228844, 11.975586222486642`}},
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 
  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,
himyss's avatar
3He
a9d10661
 
himyss committed
371 372 373 374 
   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,
himyss's avatar
radii normed
9c694860
 
himyss committed
375 376 377 
   3.8800863385534143`*^9}, {3.880086376964778*^9, 3.880086509882822*^9}, {
   3.880087582924012*^9, 3.8800876543745813`*^9}},
 CellLabel->"Out[90]=",ExpressionUUID->"d721930c-928c-45cb-8146-dfcae2431ec5"]
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
378 379 380 381 382 
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
himyss's avatar
3He
a9d10661
 
himyss committed
383 384 385 
 RowBox[{"myU", "=", 
  RowBox[{"U", "/.", 
   RowBox[{"FindRoot", "[",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
386 
    RowBox[{
himyss's avatar
3He
a9d10661
 
himyss committed
387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 
     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[{"{",
himyss's avatar
radii normed
9c694860
 
himyss committed
422 
      RowBox[{"U", ",", "170"}], "}"}]}], "]"}]}]}]], "Input",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
423 424 425 426 427 428 429 
 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},
himyss's avatar
3He
a9d10661
 
himyss committed
430 431 432 
   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,
himyss's avatar
radii normed
9c694860
 
himyss committed
433 434 
   3.8800863417389393`*^9}, 3.8800863843913307`*^9, 3.880087663509973*^9},
 CellLabel->"In[96]:=",ExpressionUUID->"22f1491b-05f8-4461-9d8a-7023940b0eba"],
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
435
himyss's avatar
radii normed
9c694860
 
himyss committed
436 
Cell[BoxData["178.1243246255173`"], "Output",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
437 438 439 440 441 442 443 444 
 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,
himyss's avatar
3He
a9d10661
 
himyss committed
445 446 447 448 449 
   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, {
himyss's avatar
radii normed
9c694860
 
himyss committed
450 451 452 
   3.880086385464319*^9, 3.8800865107303*^9}, {3.8800875839389973`*^9, 
   3.880087671218701*^9}},
 CellLabel->"Out[96]=",ExpressionUUID->"8446f902-e811-45cc-a8bf-ab7c70890a1f"]
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
453 454 
}, Open  ]],
himyss's avatar
3He
a9d10661
 
himyss committed
455 456 457 458 459 460 461 462 463 464 465 466 
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"],
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
467 468 469 
Cell[CellGroupData[{

Cell[BoxData[
himyss's avatar
3He
a9d10661
 
himyss committed
470 471 472 
 RowBox[{"myCoeff", "=", 
  RowBox[{"coeff", "/.", 
   RowBox[{"FindRoot", "[",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
473 
    RowBox[{
himyss's avatar
3He
a9d10661
 
himyss committed
474 
     RowBox[{
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
475 
      RowBox[{
himyss's avatar
3He
a9d10661
 
himyss committed
476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 
       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",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
494 495 
 CellChangeTimes->{{3.878289683612082*^9, 3.8782897289455433`*^9}, {
   3.878289801822311*^9, 3.87828982476785*^9}, 3.879569118948894*^9, {
himyss's avatar
3He
a9d10661
 
himyss committed
496 497 
   3.879571556146941*^9, 3.87957155754856*^9}, {3.880085925692211*^9, 
   3.8800859424236307`*^9}},
himyss's avatar
radii normed
9c694860
 
himyss committed
498 
 CellLabel->"In[97]:=",ExpressionUUID->"444da03a-39cf-4c3d-bb1b-9d322299f630"],
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
499
himyss's avatar
radii normed
9c694860
 
himyss committed
500 
Cell[BoxData["1.344522780057045`"], "Output",
himyss's avatar
3He
a9d10661
 
himyss committed
501 502 503 504 505 506 
 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,
himyss's avatar
radii normed
9c694860
 
himyss committed
507 508 
   3.880086511995571*^9}, {3.8800875856465683`*^9, 3.880087672538249*^9}},
 CellLabel->"Out[97]=",ExpressionUUID->"fbc9f7e4-1f37-4a69-b281-392c36da0aff"]
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
509 510 511 512 513 
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
himyss's avatar
3He
a9d10661
 
himyss committed
514 515 516 
 RowBox[{"myNorm", "=", 
  RowBox[{"A", "/.", 
   RowBox[{"FindRoot", "[",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
517 518 
    RowBox[{
     RowBox[{
himyss's avatar
3He
a9d10661
 
himyss committed
519 
      RowBox[{
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
520 
       RowBox[{
himyss's avatar
3He
a9d10661
 
himyss committed
521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 
        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"}], "}"}]}], "]"}]}], "+",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
536 537 
       RowBox[{
        SuperscriptBox[
himyss's avatar
3He
a9d10661
 
himyss committed
538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 
         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",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
554 555 556 557 558 559 
 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}, {
himyss's avatar
3He
a9d10661
 
himyss committed
560 561 
   3.879571567743319*^9, 3.879571576536313*^9}, {3.880085952401537*^9, 
   3.8800859833274717`*^9}},
himyss's avatar
radii normed
9c694860
 
himyss committed
562 
 CellLabel->"In[98]:=",ExpressionUUID->"7835cc6b-0406-456f-a3a3-ab6f89b509c7"],
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
563
himyss's avatar
radii normed
9c694860
 
himyss committed
564 
Cell[BoxData["0.8003232743677815`"], "Output",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
565 566 
 CellChangeTimes->{
  3.879569148123713*^9, 3.8795692505602217`*^9, 3.879569281919847*^9, {
himyss's avatar
3He
a9d10661
 
himyss committed
567 568 569 570 
   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},
himyss's avatar
radii normed
9c694860
 
himyss committed
571 572 573 
   3.8800863471623774`*^9, {3.880086388564464*^9, 3.8800865137393713`*^9}, {
   3.880087588111286*^9, 3.88008767445716*^9}},
 CellLabel->"Out[98]=",ExpressionUUID->"fb1c1bd4-f3ef-407c-ae27-b9cbd3bbd016"]
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
574 575 576 577 578 579 580 581 582 583 584 585 
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"Piecewise", "[", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
himyss's avatar
3He
a9d10661
 
himyss committed
586 
        RowBox[{"myNorm", " ",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
587 588 589 590 591 592 593 
         RowBox[{"fIn", "[", 
          RowBox[{
           FractionBox[
            SqrtBox[
             RowBox[{"2", " ", "mass", " ", 
              RowBox[{"(", 
               RowBox[{
himyss's avatar
3He
a9d10661
 
himyss committed
594 595 
                RowBox[{"-", "Esep"}], "+", "myU"}], ")"}]}]], "p"], ",", "r",
            ",", "0"}], "]"}]}], ",",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
596 597 598 
        RowBox[{"r", "<", "range"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{
himyss's avatar
3He
a9d10661
 
himyss committed
599 
        RowBox[{"myNorm", " ", "myCoeff", " ",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
600 601 602 603 604 605 606 607 608 609 
         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}, {
himyss's avatar
3He
a9d10661
 
himyss committed
610 611 612 
   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}},
himyss's avatar
radii normed
9c694860
 
himyss committed
613 
 CellLabel->"In[99]:=",ExpressionUUID->"fc251146-842e-4797-bd33-46971b36a427"],
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
614 615 616 617 618 619 

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    TagBox[
     {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[CompressedData["
himyss's avatar
radii normed
9c694860
 
himyss committed
620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 645 646 647 648 649 650 651 652 653 654 655 656 657 658 659 660 661 662 663 664 665 666 667 668 669 670 671 672 673 674 675 676 677 678 679 
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=
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
680 
       "]]},
himyss's avatar
radii normed
9c694860
 
himyss committed
681 
     Annotation[#, "Charting`Private`Tag$22258#1"]& ]}, {}},
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
682 683 684 685 686 687 688 689 690 691 692 693 694 695 696 697 698 699 700 701 702 703 704 705 706 707 708 709 710 711 712 
  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]]}& )}},
himyss's avatar
radii normed
9c694860
 
himyss committed
713 
  PlotRange->{{0, 10}, {0., 0.8003187560651602}},
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
714 715 716 717 718 719 720 721 
  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},
himyss's avatar
3He
a9d10661
 
himyss committed
722 723 724 
   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},
himyss's avatar
radii normed
9c694860
 
himyss committed
725 726 727 
   3.88008634841398*^9, {3.880086390133037*^9, 3.8800865180731173`*^9}, {
   3.880087591163252*^9, 3.880087675790471*^9}},
 CellLabel->"Out[99]=",ExpressionUUID->"a19709bd-12a5-47c1-be49-b078c7a3cbb6"]
himyss's avatar
3He
a9d10661
 
himyss committed
728 729 730 731 732 733 734 735 
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{
   "R", " ", "\:0434\:043b\:044f", " ", "\:0433\:0435\:043b\:0438\:044f3"}], 
   " ", "==", " ",
himyss's avatar
radii normed
9c694860
 
himyss committed
736 737 738 739 
   RowBox[{"1.561", " ", "\:0444\:043c"}]}], " ", "*)"}]], "Input",
 CellChangeTimes->{{3.8800865535980377`*^9, 3.8800865699855547`*^9}, {
  3.880087570285039*^9, 
  3.880087570627646*^9}},ExpressionUUID->"700162e1-a626-4c8d-a329-\
himyss's avatar
3He
a9d10661
 
himyss committed
740 741 742 743 744 745 746 747 748 749 750 751 752 753 754 755 756 757 758 759 760 761 762 763 764 765 766 767 768 769 770 771 772 773 774 775 776 777 
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->
himyss's avatar
radii normed
9c694860
 
himyss committed
778 
  "In[100]:=",ExpressionUUID->"1a4970d6-a8d9-4293-a27c-2d4dbe2240ea"],
himyss's avatar
3He
a9d10661
 
himyss committed
779
himyss's avatar
radii normed
9c694860
 
himyss committed
780 
Cell[BoxData["1.5612379745157818`"], "Output",
himyss's avatar
3He
a9d10661
 
himyss committed
781 
 CellChangeTimes->{{3.880086253659615*^9, 3.880086304157057*^9},
himyss's avatar
radii normed
9c694860
 
himyss committed
782 783 
   3.880086352535798*^9, {3.8800863938671637`*^9, 3.880086521539238*^9}, {
   3.880087593758895*^9, 3.880087646136713*^9}, 3.8800876779743423`*^9},
himyss's avatar
3He
a9d10661
 
himyss committed
784 
 CellLabel->
himyss's avatar
radii normed
9c694860
 
himyss committed
785 
  "Out[100]=",ExpressionUUID->"f2064ae6-052c-4c29-85b6-4122f8ffe9c5"]
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
786 787 788 789 790 791 792 793 794 795 796 797 798 799 800 801 802 803 804 
}, 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[{
himyss's avatar
radii normed
9c694860
 
himyss committed
805 806 
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"],
himyss's avatar
3He
a9d10661
 
himyss committed
807 
Cell[CellGroupData[{
himyss's avatar
radii normed
9c694860
 
himyss committed
808 809 
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"]
himyss's avatar
3He
a9d10661
 
himyss committed
810 
}, Open  ]],
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
811 
Cell[CellGroupData[{
himyss's avatar
radii normed
9c694860
 
himyss committed
812 813 
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"]
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
814 
}, Open  ]],
himyss's avatar
radii normed
9c694860
 
himyss committed
815 
Cell[20515, 454, 420, 10, 39, "Input",ExpressionUUID->"0b790137-d050-4e54-820e-ce3bbd34ede5"],
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
816 
Cell[CellGroupData[{
himyss's avatar
radii normed
9c694860
 
himyss committed
817 818 
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"]
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
819 820 
}, Open  ]],
Cell[CellGroupData[{
himyss's avatar
radii normed
9c694860
 
himyss committed
821 822 
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"]
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
823 824 
}, Open  ]],
Cell[CellGroupData[{
himyss's avatar
radii normed
9c694860
 
himyss committed
825 826 
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"]
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
827 
}, Open  ]],
himyss's avatar
radii normed
9c694860
 
himyss committed
828 
Cell[32423, 729, 397, 10, 29, "Input",ExpressionUUID->"700162e1-a626-4c8d-a329-bbb787d8421f"],
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
829 
Cell[CellGroupData[{
himyss's avatar
radii normed
9c694860
 
himyss committed
830 831 
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"]
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
832 833 834 835 836 
}, Open  ]]
}
]
*)