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
3He
a9d10661
 
himyss committed
13 14 15 16 
NotebookDataLength[     36218,        839]
NotebookOptionsPosition[     34231,        798]
NotebookOutlinePosition[     34629,        814]
CellTagsIndexPosition[     34586,        811]
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
3He
a9d10661
 
himyss committed
59 
  RowBox[{"range", "=", "0.912"}], ";"}]}], "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 71 72 73 
   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"],
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
3He
a9d10661
 
himyss committed
149 150 
 CellLabel->
  "In[357]:=",ExpressionUUID->"2fc4ef7c-d275-441b-9d3d-7d7239f74130"],
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
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 190 

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
191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 
    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"],
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
212 213 214 215 216 

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    TagBox[
     {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
himyss's avatar
3He
a9d10661
 
himyss committed
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 
      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==
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
267 
       "]], LineBox[CompressedData["
himyss's avatar
3He
a9d10661
 
himyss committed
268 269 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 315 316 317 318 319 320 321 322 323 324 
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"]& ], {}}, {}},
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
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 353 354 355 
  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
3He
a9d10661
 
himyss committed
356 
  PlotRange->{{0, 200}, {-5.332358609904164, 10.422205589087188`}},
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 
  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
374 375 376 377 378 379 380 
   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"]
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
381 382 383 384 385 
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
himyss's avatar
3He
a9d10661
 
himyss committed
386 387 388 
 RowBox[{"myU", "=", 
  RowBox[{"U", "/.", 
   RowBox[{"FindRoot", "[",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
389 
    RowBox[{
himyss's avatar
3He
a9d10661
 
himyss committed
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 422 423 424 425 
     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",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
426 427 428 429 430 431 432 
 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
433 434 435 436 437 438 
   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"],
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
439
himyss's avatar
3He
a9d10661
 
himyss committed
440 
Cell[BoxData["124.41224361311541`"], "Output",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
441 442 443 444 445 446 447 448 
 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
449 450 451 452 453 454 455 456 
   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"]
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
457 458 
}, Open  ]],
himyss's avatar
3He
a9d10661
 
himyss committed
459 460 461 462 463 464 465 466 467 468 469 470 
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
471 472 473 
Cell[CellGroupData[{

Cell[BoxData[
himyss's avatar
3He
a9d10661
 
himyss committed
474 475 476 
 RowBox[{"myCoeff", "=", 
  RowBox[{"coeff", "/.", 
   RowBox[{"FindRoot", "[",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
477 
    RowBox[{
himyss's avatar
3He
a9d10661
 
himyss committed
478 
     RowBox[{
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
479 
      RowBox[{
himyss's avatar
3He
a9d10661
 
himyss committed
480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 
       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
498 499 
 CellChangeTimes->{{3.878289683612082*^9, 3.8782897289455433`*^9}, {
   3.878289801822311*^9, 3.87828982476785*^9}, 3.879569118948894*^9, {
himyss's avatar
3He
a9d10661
 
himyss committed
500 501 502 503 
   3.879571556146941*^9, 3.87957155754856*^9}, {3.880085925692211*^9, 
   3.8800859424236307`*^9}},
 CellLabel->
  "In[365]:=",ExpressionUUID->"444da03a-39cf-4c3d-bb1b-9d322299f630"],
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
504
himyss's avatar
3He
a9d10661
 
himyss committed
505 506 507 508 509 510 511 512 513 514 
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"]
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
515 516 517 518 519 
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
himyss's avatar
3He
a9d10661
 
himyss committed
520 521 522 
 RowBox[{"myNorm", "=", 
  RowBox[{"A", "/.", 
   RowBox[{"FindRoot", "[",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
523 524 
    RowBox[{
     RowBox[{
himyss's avatar
3He
a9d10661
 
himyss committed
525 
      RowBox[{
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
526 
       RowBox[{
himyss's avatar
3He
a9d10661
 
himyss committed
527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 
        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
542 543 
       RowBox[{
        SuperscriptBox[
himyss's avatar
3He
a9d10661
 
himyss committed
544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 
         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
560 561 562 563 564 565 
 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
566 567 568 569 
   3.879571567743319*^9, 3.879571576536313*^9}, {3.880085952401537*^9, 
   3.8800859833274717`*^9}},
 CellLabel->
  "In[366]:=",ExpressionUUID->"7835cc6b-0406-456f-a3a3-ab6f89b509c7"],
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
570
himyss's avatar
3He
a9d10661
 
himyss committed
571 
Cell[BoxData["0.7793879594854948`"], "Output",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
572 573 
 CellChangeTimes->{
  3.879569148123713*^9, 3.8795692505602217`*^9, 3.879569281919847*^9, {
himyss's avatar
3He
a9d10661
 
himyss committed
574 575 576 577 578 579 580 
   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"]
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
581 582 583 584 585 586 587 588 589 590 591 592 
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"Piecewise", "[", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
himyss's avatar
3He
a9d10661
 
himyss committed
593 
        RowBox[{"myNorm", " ",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
594 595 596 597 598 599 600 
         RowBox[{"fIn", "[", 
          RowBox[{
           FractionBox[
            SqrtBox[
             RowBox[{"2", " ", "mass", " ", 
              RowBox[{"(", 
               RowBox[{
himyss's avatar
3He
a9d10661
 
himyss committed
601 602 
                RowBox[{"-", "Esep"}], "+", "myU"}], ")"}]}]], "p"], ",", "r",
            ",", "0"}], "]"}]}], ",",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
603 604 605 
        RowBox[{"r", "<", "range"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{
himyss's avatar
3He
a9d10661
 
himyss committed
606 
        RowBox[{"myNorm", " ", "myCoeff", " ",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
607 608 609 610 611 612 613 614 615 616 
         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
617 618 619 620 621 
   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"],
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
622 623 624 625 626 627 

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    TagBox[
     {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[CompressedData["
himyss's avatar
3He
a9d10661
 
himyss committed
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 680 681 682 683 684 685 686 687 688 689 690 691 692 
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=
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
693 
       "]]},
himyss's avatar
3He
a9d10661
 
himyss committed
694 
     Annotation[#, "Charting`Private`Tag$88641#1"]& ]}, {}},
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
695 696 697 698 699 700 701 702 703 704 705 706 707 708 709 710 711 712 713 714 715 716 717 718 719 720 721 722 723 724 725 
  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
3He
a9d10661
 
himyss committed
726 
  PlotRange->{{0, 10}, {0., 0.7793861446401996}},
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
727 728 729 730 731 732 733 734 
  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
735 736 737 738 739 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 778 779 780 781 782 783 784 785 786 787 788 789 790 791 792 793 794 795 796 
   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"]
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
797 798 799 800 801 802 803 804 805 806 807 808 809 810 811 812 813 814 815 
}, 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
3He
a9d10661
 
himyss committed
816 817 818 819 820 821 
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  ]],
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
822 
Cell[CellGroupData[{
himyss's avatar
3He
a9d10661
 
himyss committed
823 824 
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"]
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
825 
}, Open  ]],
himyss's avatar
3He
a9d10661
 
himyss committed
826 
Cell[20413, 458, 420, 10, 39, "Input",ExpressionUUID->"0b790137-d050-4e54-820e-ce3bbd34ede5"],
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
827 
Cell[CellGroupData[{
himyss's avatar
3He
a9d10661
 
himyss committed
828 829 
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"]
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
830 831 
}, Open  ]],
Cell[CellGroupData[{
himyss's avatar
3He
a9d10661
 
himyss committed
832 833 
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"]
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
834 835 
}, Open  ]],
Cell[CellGroupData[{
himyss's avatar
3He
a9d10661
 
himyss committed
836 837 
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"]
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
838 
}, Open  ]],
himyss's avatar
3He
a9d10661
 
himyss committed
839 
Cell[32455, 742, 348, 9, 29, "Input",ExpressionUUID->"700162e1-a626-4c8d-a329-bbb787d8421f"],
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
840 
Cell[CellGroupData[{
himyss's avatar
3He
a9d10661
 
himyss committed
841 842 
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"]
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
843 844 845 846 847 
}, Open  ]]
}
]
*)