Skip to content

M
mathematica
Switch branch/tag
Find file
Normal viewHistoryPermalink
bounded8He.nb 46.9 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[     46895,       1013]
NotebookOptionsPosition[     45096,        974]
NotebookOutlinePosition[     45494,        990]
CellTagsIndexPosition[     45451,        987]
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 46 47 48 49 50 51 52 53 54 55 56 57 58 
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", "=", "3.925"}], " ", ";"}], 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"binding", " ", "energy", " ", 
     SuperscriptBox[
      RowBox[{"of", " "}], "8"], "He"}], ",", " ", "MeV"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Esep", "=", "24.81432"}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"1", "p", " ", "separation", " ", "energy", " ", 
     SuperscriptBox[
      RowBox[{"for", " "}], "8"], "He"}], ",", " ", "MeV"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"p", "=", "197.327"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"mass", "=", "821"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
himyss's avatar
radii normed
9c694860
 
himyss committed
59 
  RowBox[{"range", "=", "2.399"}], ";"}]}], "Input",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
60 61 62 63 64 65 66 
 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,
himyss's avatar
radii normed
9c694860
 
himyss committed
67 68 69 70 
   3.879568728576707*^9}, {3.880087700774041*^9, 3.880087701627767*^9}, {
   3.8800877975749063`*^9, 3.880087987460205*^9}},
 CellLabel->
  "In[218]:=",ExpressionUUID->"ad6df89e-d3d8-490c-ae91-dac6d80d1b5a"],
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
71 72 73 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 

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
146 147 
 CellLabel->
  "In[225]:=",ExpressionUUID->"2fc4ef7c-d275-441b-9d3d-7d7239f74130"],
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
148 149 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 

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
radii normed
9c694860
 
himyss committed
188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 
    RowBox[{"U", ",", "0", ",", "150"}], "}"}]}], "]"}]], "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.880087707720909*^9},
 CellLabel->
  "In[231]:=",ExpressionUUID->"c2ad0be7-22ee-4300-8ad8-45a6e113c1b0"],
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
206 207 208 209 210 

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    TagBox[
     {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
himyss's avatar
radii normed
9c694860
 
himyss committed
211 212 213 214 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 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 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 
      1.], LineBox[{{3.061224489795918*^-6, 1.940807922738254}, {
       3.061224489795918*^-6, 1.940807922738254}}], 
      LineBox[{{2.282297965209684, 1.9845387764422011`}, {2.9941734674575162`,
        1.9993610965838908`}, {2.99539175357331, 1.9993869892252918`}}], 
      LineBox[{{3.8625306041212233`, 2.018295309200041}, {6.136253562818797, 
       2.0729010424770893`}, {6.182822765517297, 2.0741037793717307`}, {
       6.229391968215797, 2.075310175845261}, {6.322530373612798, 
       2.0777340178746426`}, {6.508807184406798, 2.0826263971026}, {
       6.881360805994801, 2.0925940258106626`}, {7.626468049170805, 
       2.113295145936034}, {8.843893520444967, 2.14950944837228}}], 
      LineBox[{{10.791098607617133`, 2.2145823816793007`}, {
       11.308156126004256`, 2.2335674387661397`}, {12.038647322831405`, 
       2.261777314434328}, {12.088172683001606`, 2.2637520106504767`}, {
       12.137698043171804`, 2.265734881901272}, {12.236748763512203`, 
       2.2697253612238386`}, {12.434850204193001`, 2.2778067751671274`}, {
       12.831053085554595`, 2.2943838922288315`}, {13.623458848277785`, 
       2.3293019531410413`}, {13.672984208447986`, 2.3315669181473053`}, {
       13.722509568618184`, 2.3338420034147838`}, {13.821560288958583`, 
       2.338422816820904}, {13.979240340272746`, 2.3458005797934325`}}], 
      LineBox[{{17.861058649914543`, 2.568496428309608}, {18.182932399399863`,
        2.5914197321642227`}}], 
      LineBox[{{20.612359074679755`, 2.7965701464575736`}, {
       21.384177652356154`, 2.8766877426854762`}}], 
      LineBox[{{23.74486483122288, 3.188517520757271}, {24.5194156385197, 
       3.320930813286054}, {24.545732075444665`, 3.3257699818548376`}}], 
      LineBox[CompressedData["
1:eJwVy31QkwUAx/Gx22xsbD4Pe2luS4FDwESyJZO9PM9+jgWdoCODAaHkuNuJ
gAdYm8hUVIidIA0QbQiIZrwYJBfd5RRFC08pIenAViBnIy44okA5a3qW2R/f
+/z1Dc0p3GZlMhgM1Yv+Nz1pJDf+Rhy6pTEfbpzjgn7zKus+1Khks6sOXeaB
YTDQp5PVyDs7FeAY5uEbenBfRoYamzlPGzt8PCTGjc2NFaoh8ny8eIIbhJS1
8yPDZ9QwpzfmDGcHwUIqWq4/U6PC9FVgMclHuOCTHw9yNKisj8/8LoqP37hR
hE6kgafV1rwTfOSxNpR71mqgUQ7eSSziY+/fSbu/eFeDo8w19YpRPiomHbHn
PRqEFx/OHmkTIGE8oNByUwP7zv4j9q8F4HidnatGNPCfCizInxSgZuSEvHlG
g72Dm5cSJMtxaqAr4KREi4Z2nmy+ejk6LkwMV36ghc/35EmYkUCI5eg5wxEt
kmJKDa+ZCbila2z/1mixyrzscUkugSqnXWHv0GKlzez+toZAgZXMt45rkRzZ
NFQxQWBa4aFDZ7RgTJnCzi0Q2DGWHTy5pEXIRe4Ym0Via3z35VS+DqwdzsqL
60isD03kGKGD/N6ejK4KEp0//THxPFmHp5leW1wTidDahp6+TB26+4Vt8l4S
wc995jfe14GtmhPP/0Li8f0DbWHtOvz68q0bC4ZgXHH3GgKCKLBcIUSdTAhl
SqbkmpSCkT1vdm0UouslxlzJagqnfQ4nL02IZvuWukWawqafP7teVSdEWers
gwfFL36H37uMFMHPc33ZdIjCYsF2plkpQtFArDO9mgL5T3o5UkXIUZbH3P2U
gvWeacbfKIKRWHmw30shyhQRW7pOjKu3b6aUTlNQNQrIiFQxYsvyw1UPKWyr
7mltcIgR8afnzudcGotT0Vk/DIkROPSOrIWi0SVRDjr2S/C989glx3kaU/sY
fa8wpch6K/hhdg8N9y5lS+cGKWY5Ta8a+mj4rSjZkisFs6q7hTNKw8i/MM0d
lUJ1/G75SaYedNyzWdulFRhIzrhWItDDcng2Y2FpBVL4vr+yZHqobrvaxl+X
Ie+jR7vDlHq4hN0JVK8MrbXilB6LHu8dOCv13pIj+u0zx+r36NGexvy9jFDg
Chk5YNuvhxv8aPl2BdLa6xiRtXqs37SryPRIgf8APFhbCQ==
       "]], LineBox[CompressedData["
1:eJwVjHc41f0fxu2shKzsRFl5tJD1eeucYxyjiLPzIFJCRSGrEk0rRROJjBPy
EBXqe5KQpCOEyMpWOlkd+9fvj/u6r9f1uq97s9cJFx8+Hh6eyb/5fw90BQ9G
SuGgMndeZrvORlZIRoFz6UYcmI+ZF/cIbmSltJ5tElDEQWCwXlxftjTrg4Vi
TcFWHITgwg0zx6RY5tL0ohmEA09fL8adG5Isio2etrU1DkhmV5/WWEmyzkQu
PbrjgAPXWsvwXWMbWMUj9+9YUHEw1DAcPPFKgqVW1XPxchAOvF9KV/qkibPM
OYWrXWE4uCN5sUhWVpxF0Yo6q38OByM3+WWn08RYyUkqJ1qu46CTE6RUlS3K
4vN2pyk9xsH5J0NGj3uFWap3DdoDnuBAffHQcm2wMMu0eXU/6z8cCJU5p9uI
CbOCTTIJPq9x8NBXSE7Lah1rSLzf8GkHDnKb1bMs6gRZ16yz3FO/4SDsvbKk
e5Aga8d5r/iI7ziYjbhnuUldkBU9MzRq8wsHR/ZWuu84L8BS+Dqe0S+Eh2jd
778rnflZrzc+aaoTx0PaH2fb8yL8LB9H/8VCaTzcwR9rTajhY5WyptzOquJh
a1L2irsJH8sxb0Zc2ggP/tRjltwdvKzZ/memXDM8HPDjuaI0w8O6pxhytNcK
D9LlS4Z25TyssQTuW6YjHiQ6KgILzHlYMadXwnE+eGgReFiRXbSKVVitGz99
Cw9vaWYXln2XMEZEgxz9Hh6W3JPGdwktYfzlV/FWD/Hghqfv7ny8iB3QFn+4
vhAPItXx5TbjC9iEhBQ59y0efHa9fsc5x8VUexTfdUzjIdvFpCJuZA6rle2Z
fsXFw/3xMAu4M4cd35+unrOKh/3rn0atEeewFzVqkSdFCeC7P6k9tGwWcynQ
3CWiQYDUW6eaG5JmsMshBlmmzgQYHn3VKRP0G7uynudBJpkAP20O6WTq/8au
5bSkCbgToL7NKF5sjIMlfj4d3+xHALWVDVZnPTlYqkFVqFcsATgq20b8qqew
7BEbp2vPCfBiGFfW/mcCy4naZPfrFQG0DZv9r8ZNYLkykzjXWgI4mBHNE2Qm
MOa+xL2qLQSARLP3AUbj2H8ZbZqlEwQIVP5UPH9lFMNInktdKtbQb2JgW+Q9
hLF+7py31LQGQ7Jn06jAEFYTK/A7W9camoJRp0bud6yuNG8kwNgabO0PG1+b
HMQ+Sky18Dpbg6b9SmbJhQGsuy48XyfWGjx0bsccLvmG/TG65XZ2whrgRKnU
YcM27DTB05PFsYZaX4gNyW3FOAcNAtb9sYZ8bgq7WqUVGzvZEJvKbwMeMyrX
WzZ8xr4yl8v+U7YBirPIvUJeNlat6iM94WQDB3wfx7buaMAuCBl9opXZgPoz
D7LXYhm2KsPfnfXSBkqKRIdFnEux8C3skTHsL/fqHvnOLMGCwG81pOnv3swn
X4zNxLzDM7bfHLYBnjw7RwnlB5j1lFD8B3lbUK8TfG2J3UJiXzpszSP/8q93
bk15VUh6PHWdfIwtKGyd2Gg1VI0Ulg/W/b5sC5zZUzdrVV+jrRpsXN5NWzCc
3ZMvaMNCENhgKf3EFqqTvfdLPKlBp4Ve7h7rsoU7b0IFNGzrUY/RPfVbxnag
+fzf32+GPqFBO0pfoIUdrFKGFMR3s9E4Qy7DDmcH8Ymv7Ohn2Gg+JkVpzckO
DE5StMbn2Ujy0zU5vyN2wGpN2N0114IIvhHikGoHpNTSrIrvraj4NoM7OW0H
Q9NScuTYDiQprBdLWLADduyNoEesDhQUtrAhc80OHtEes+8vd6A91LRtLuJE
YK5C1P3gTlStyCa90CKCdveGm4OHulBjOq48lkIETSWtbm3NbjSSrRuk8poI
t17o7F+J6EW2MgsrIbVEsDlKtnB62ouYsfVX2Y1EaDFxm/cY7EUBPt5ZFzuI
8KUruP2jTR+a3ZrBHucQofb4Yeb4un7Ex5T6p2KLPTCFJV9GNvYj1RLupNNV
e0g1N2YfMxtEolovdM2T7KEiWntq24FBNH839JhOqj18DiyQbvEeRM0X50f4
s+wh4VmfyHDCIIomzw68eGEPDYeiOuT7BlHv6q8OjTF7GJksF3sf9R1lOI6+
nbd2gGV+tTBi0RBSmWy/nyHgCJapgd/8ukeQ96vHeaJijtDYfeKj/c8RVJh0
pixEyhH8Q58MK6+NINPdsh+cVB0h5Esjb7LGKCJFuS6uGjuCk7YBz6ujoyh5
Qyvl3+OOcCBmTZI4PYoEdrNl1T87gqjxpU3Di2NoIrIxKSvTCYhmEcOaoxOo
KMTJ89amAzCXkVoVbTGFtkMnpbn9AFhJhNm+ieKgjeUB0wNZzmBq2sB16/qN
rKq2TBz3cIGcXyM6K2Iz6E2JLa+J0kFIlv/ng/TOWZTnfLby2/BB2H487ILD
6TnUfrGxZeGRKzyc3TOgXTyP/Dy3Xf91wg3GXEynv6lwkbPBzOFte0lgWFfT
975sATX2NSvLmZNA2e0eSbNyAe27wWwXQCRY/My4EMNaQLtnPW0G8STg2c8v
Svy4gOQq2brpB0iQvTk6QWh0AfUQin9LHyXBTv4zBaZKi+jooWPnVtNIED3W
n1N1aRGdi++71z5LAj52Oif/+BLK9cdV8nNJ8Mhss3x98BJqcsjr2rFEAm6o
8vPJiCWksP6EQiIvGYJN1t+wjV9CpQmrqdYSZAjz6+c4Fy2hoUSVG8+3kaHx
nX9zEWcJEW/QL9+lkeFnzNbHg5HL6OQpLLf+EBluRZ7oHY1bRmnOW+rmPMgg
57/JYC5xGX2XmhRw8SWDdEm6vXbWMopMCY8RPU0GDdOuLdPvltHTm3ejIhLI
4LAt8+tmqRUkk9YRfIhFBsuZO44+xSvoPLfM3+0tGfh6ch36nq+gSVqyj2Md
GUJAycDrzQqqUbMjWzSRIUasuSOmbQWdLKg0VekkA0Fz1/UdSyvoQ3U677df
ZHD71iR5kLiKYr4fTmKoUYArd+b+Fs4qmiLAFVcNCginPY3pXlxF1HzlCw5a
FCBYDNzMEFxDhgHtQeZ6FEhxfnHFWmkN9c5bk5SNKZAkaC6vbrOG9orqqvQ4
UUDOguS4PmsNcQw5THo0BfBCJ9fM/uGBTfuJhfgLFJDQ0hXZt4cH9gXkFG2P
pcBgi/64kxkP3CyglKxdpUCYga1XlA0PGGm8Kc9OpYA6WbHawIMHIjemsCYL
KaCbvTvzVwoPrJvf9SWimwJtYckvSas8YCiT2OHdSwG9SY10V0FeoOwc63Qc
+Pvn9P4gXZwXmIEPutVGKbAevYq7qMgLDqMCA2+nKXClGJO3MuaFlM72STFR
Kng128bFBPGCalUo7wMTKlw3lTv+nMMLs+tox2rNqHBqF6e7cIEXGl3NW35Y
UqE42pi3hI8PQqd4siwIVBgzip7pk+ED9uZr0OdMhafvSNSJvXxw8cqDcxrH
qRBSat/8/hIfjLm9WclPp0JT1XunNn1+eP0o26flIRVcz90cbzLmh9RfcR8X
sqnQw/BK6tzHD/uuEtOJTCq0vY1v0aHyw/3qVvMfFVSIb4GdlZf4wWnLSMQ/
bCoYFR41eDfMD2Uc0cVyPhpIkH8Q7MsFoBnzMcsXpEGxydylf2oEYCKRFXlP
mAaq/HdL9T4JgPr2kLVoCRoMLZh99x8XgMRjAwK2ijT44ZJ1V1FNEPwGn2/4
uoMGpLPuOUeSBOHSf9LOTbtp4FgoN/EjXRCyzgekvDamQYlwsXByoSB0qm6R
zbaggVLjo70GjYJgTU9U9LejgfYeiYw2ISHY3OattepBA9cBEZz6JSEwz8aO
cA7TILnSJS8rTQgoQYr5g0doMNFlirfOE4JESbZuvT8NxHr+rC42CMGyg5lh
chgNFCTHPqpJrIPOd1JmGsk0GCN+ZDdkroOk568PWGM0oO2u//fpkDBQyrO4
oTU0OMoXoxK3LAwaZbEPC97R4I754PH7MiLwrJjIEW+iQX969kICQQQ6H39J
au2kQc7hryMFTBFQvzXV7DFNA6EKPtyHGFEYv9ESkjJHg6/9D34YZYpCadIz
1VouDby4l43+VIkC/np44LY1GkhTPr2NmBeFozFCElNidBDiiFx5c0IMSk6q
OkZo0YFHLDi24bQ4nA3knSvUpoNo1EIU8bY47PMfetCrRwePR6c8LSvFoc2X
+cNqJx2MbgVyY3nXA9fdKF4Y0eHSlzN0dup6sHJ0+pBKoUOORHnIfx0SIGq/
I7ieTofpivzuOL4N0GYro7TgToe0tup87vYNcAT/1Y/hQwcNFce0+ksb4JrZ
EdEtQXSI6smv2asuCZ91ou1KrtPh/ZUnSp2TktBf0NX4JJEOVXXja26CUjCl
s8ch7wYdFKpKXxmrSYGo7g+njNt0ML/vo2foKgU4XZprfDYdbjqzO0zfSoEL
s7z9ci4dHq132aTSLwUeulLkiwV0sJt8KRW2IgWRug3U8Kd0OMnMJuTslYZn
ukb/Hq2iw05mQ4DhC2moYd7oP/yaDk95lUJ3dkkDW/en579v6OCe1ufcsCgN
P3VzvEn1dDiv+MYLQxthq560H76VDrXFS1ddWjbCbb2p0+qTdGieLvlAlZeF
3Cd280pTdFD8qLOnGScLz/Qeh8r//uvF99bonpIFth4jXOIPHcAsNeBlsywI
6zeeW+ZjwHVq4bhomhzIFWrxcgUZcCHV5fSDRjnQ1L8QMyPMgBa/h0DlkQcr
fZO4CQkGyKRHPIw6IQ/h+rnXOhUZoBmzJFR4SAGuFPKKt6kwoD9lC/XTAwVI
0z+U8EmdAcZfpDKivilAqb5Mct1WBgzovUiqOLwJJvVjUp/tZMARLFrNN04R
Fgu/yZXsYUBgtlgR87MiCG/fe+eJCQPitV39hTSUQHM7594jSwaoMq19POqV
YOR4cc41WwYYWvme3D+tDP8DgMXNgg==
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
340 
       "]], LineBox[CompressedData["
himyss's avatar
radii normed
9c694860
 
himyss committed
341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 
1:eJwVkHtYzPkXxytNpTTd1Q+lhCQVbaVEn0MtRbGrtG30vUkSEhKtVLpoRYmS
SqjtvopQKKPPd+lidHEZoptmrDJDRjS6SPxm/zjPeV7Peb+f8z7HjNmzcZuS
goLCSXn917dYqSdNDAZCXcmIvrWlHiueTNHxlAaCZM7/Audy9NhtL1yL7WWB
sGlT2D+vt+uyREpJq4YCAQ8aZEU/3HXYWQfjshSnEFDdnpDTqKDDdm0NoEY5
BHyLndRwrNdm/VZwZa81CPiUEbO6c4U2u2HowMw6QwL2zbq0wuMCl+X2bui/
OoOAQ8vDbwYv4rKtfMurxcYERJiGnFKt12Q9CntXnTaX7zsp8IkWT2PB7+fQ
EFsC1n+1lE3312C/r5ptT9gRYNRyzWz+uDrLsx2f9HEgIN9nqcG8PHXWaWrl
aeRCgNBrovDm26nsEp5+7fQ1BHDHpUMlGWrsUJk0ftpaApSnW7tfWaXGXjn7
wEvJm4A/kl0XLx9WZRfuiRZ+2EiA1Gz+BtZflZ0z541qI0GAQnD7q/X2KqyQ
W/+kjiZA23aX0W4ph704ce58VRABWVG74xvKOezM5+ts80IJGLMO2v3dnMPq
/3nDb38kAa993aMezlNmOdL4YrM0AmziNTMl/kps+NT+536nCXiouWZXj4kS
2zN3jcrJTPk9TwsaGvsV2erNGttHcgkQt86p2hqpyAbxMxa0lBLwstDa3/ov
BfZeUdHliH8IMMu+HzcTfcc2rEpPeQMB+93CrbvfTOLc7pBpfc0EJP3yy9bB
E5N4r671bs92Ajx1LqgP9H7DpnE1NibdBAi8xSgifQLHbG661iwjYPNhlZ3J
C8fxu8gFr7+NErDdfEa4z+sx7HcmRddugoAWh3W52ufHsDV//f48RRJcou6n
r9Mawz2OL+z3ckm4Kxbor1Eaxct0xbdmLCDBNvg3v2buF1xivVa8wYoEi5VZ
9hUCGdb1rDBKsiFhdnbeerccGZbEhkd9tCdhr+HkpcPzZTj7w9iyhpUk7KAf
X2Y9hvHIg6n1uzeT8CQ++/BPf3/Cs/MlIxkECc78fsuWqE/Y4yDfto4mwfpz
eli/5yecO+94gUoICVt/9YvIeT+EXePVki5FkHB84Udxs/0QPuai6vUklQRq
b13e6toP+Kru28TRdBKy4+4PKG3/gF9Kmu4aZ5KglGVie8LgA16Yfcw2NJeE
ptzyR8MRg7hdxtGbUkrCxbpa3efL32PDK8pd9iwJJglWavFSMYakN7qb78v9
WndtakrFOHRLw7qjTSRcb3hwPYUWY5564t22VhIyfR6r6L94i6ntUwqCO0lI
T5t7MIg/gMtMlUKyP5Nwe05gK+K/wXn8H1+P6FLQ0LxVrdRShH1GrJRH9CjQ
ndZpbjEqxOrm/pphBhTUJ5IGnEYhPhRdNZswkvObUd8OWoh9bWl3VxMKWj2j
rINQH9Y8ey/1uyUFvEByyWe/HhxDJZnGrqTgpbTtc/LXDuyQes1ydBUFiaJ1
Zw+Xd+DB2l67Pe5yTvCe+cW/AwfoOf5MrqEg7mhwckTtc7y0aWAH8qbAXy+1
LS7uGR6y8rjxw5+CBROKzrx5TzE9qr46LpyCMqUjvN7jLTjM5bCS0T4Kwv9y
uJ4w/hBHx76vv7qfAlZ8RKoT+hBnq7Yu7YukID/VZcUrbz5+bJhq6XqEgpB1
gr6x2c0YOWlpfjtOQdUm5xhm8B42jtJ7drBQrp9689DOJTfxwrsJ6VrFFJia
3b/HdavBTooyr9ISCoaOJc3S8KvGvsefNXaUy/P/XtCvFHcdn8w5e8uhigLh
CDHHqbkST9Qanh/mUaCwc1fX+/OXcOfXmcyeDgpA5d9gs0eZSD3OQ2Pwhdxv
ylMdvnMOLeMcqA7plOsnF8+3Kc5FuVrtKkyPPE9rjE2MYz7aMjfusu9rOb/j
1ARIipHQ+43MWUrB4j832uvGVqG3+RXJyio0vJy8rHDqWx0ymt+5+KgqDVV3
FFsML95BHpc5XZNqch4fT3znykNlNcSiUQ25fluTYF/MXRT6UPuJRIcG4R81
KUu+YCQdjpjxyJgGNm/KbWf6HvqyekVljgMNXg59Y9JPTUj577Iu7lIaBvNn
6GgtbUb6mvpqiU40lJUGXBAfbkZ2gndMmAsNrYs6u/jKD1AYmW3ktpKGQ9JO
u1/1+Gjg4HDCoBcNFZVZWU/NW9DLsvIAFERDw1SDgl0Gj1Cd+nS1/jM0ZDpd
aAqzECDTpEZxVCYNTm63b0y4C9CxHxF8bhYNPZ8KvHcxArRRJkhxyqHB5pT5
0vE8AZL0ntY8eYkG+2jD7CbtZ8jomqbeTxU0RLYKbHyHn6HI3zgmRxv/my/r
Xl/cgeyKZD8Zj8n/8ehi4YGGTlRHUornxuVzb+Puo12dyG1ma7v2BA1FvbYj
0UOdyPdM4Q7l7zR8Naq2gVldKDJu46X3UxgIdsxUmdzXhWq3VKnXaTGwSOLe
xBh3IzDYJfKzYEDbRtF8Z3APWn/s37R0PwaMVqUxYW2v0Jq1aV6O/gxsNj6w
6IPoFQKus3rP7wwkMdHuW0ZeoSVZaUkWgQzwJGZVGrP7kH6Jc3T9VgbcempK
ru3pQ92Np3YMhjMwfdz4rIGqEIVwlrt7nGDA4taDIKt8IaL5A4rSkww0zas2
iy8XooDU0zgjjYHYiimLnl4XIi+Dt8v6TjPwY/yeJ9EoRLbzziyJzGZgBi85
1eKdEI24i02Kihl42TZgmWQnQh/VMno8SxngVt6JiXARIXHrityPZQyENHy8
GOguQp2+GQYuFQw4ZruTxn4ixAtynfb0OgNjy+9U+kSJUM0CCf9gNQNPv45a
6caL0JXBjGTjm3I/T7nqYYoI5UdIlHbUMtBu9rHfMk+EcpwyWe4dBkJzRcWP
i0TozDfXmGoeA3vftyXvrRShFFbiElDPQKJn4ZlpN0UoITFz/AdmYK716raC
ehH6P9mIOYE=
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
391 
       "]],
himyss's avatar
radii normed
9c694860
 
himyss committed
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 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 
      LineBox[{{0.580471590483438, 1.9514200875838648`}, {1.4721511207270248`,
        1.9683836516713138`}, {2.186634699903562, 1.9825925710320407`}}], 
      LineBox[{{8.939556785751089, 2.152490616868101}, {9.207993935126206, 
       2.160967792009999}, {9.299305334729599, 2.163889625808596}, {
       9.481928133936387, 2.16979276003634}, {9.847173732349962, 
       2.181842875836804}, {10.577664929177109`, 2.2069690219426485`}, {
       10.59036934281167, 2.2074186517273064`}}], 
      LineBox[{{15.835754719747317`, 2.4414299437024716`}, {
       16.687319194227936`, 2.4915397677971853`}, {16.73197212591572, 
       2.4942919725966415`}}], 
      LineBox[{{22.763523401115037`, 3.0445179991141096`}, {
       22.947972382630205`, 3.069820622958715}}], 
      LineBox[{{3.091055018879432, 2.0014259047659824`}, {3.343293478053214, 
       2.006857035698921}, {3.7422877758768687`, 2.015614922736117}, {
       3.7668673388151013`, 2.016161269426587}}], 
      LineBox[{{18.917112571488154`, 2.6470421147113985`}, {
       19.747774596335077`, 2.7162552749872737`}}], 
      LineBox[{{21.479840917662276`, 2.887241944107288}, {21.568829651773182`,
        2.8971921120413873`}, {21.76553538422295, 2.9196521759320224`}, {
       22.158946849122486`, 2.96659757335866}, {22.200608173424442`, 
       2.9717348968934307`}}], 
      LineBox[{{14.074903605578868`, 2.35032866174145}, {15.180775289873466`, 
       2.405721642181061}}], 
      LineBox[{{17.05545921690165, 2.5146262200643115`}, {17.42684360447982, 
       2.5388652813832806`}, {17.76539538460842, 2.561843832651462}}], 
      LineBox[{{0.09344920708285748, 1.9424942553287703`}, {
       0.18402156866230668`, 1.9441367148350859`}, {0.36804007610012357`, 
       1.947498215336329}, {0.48480832517731554`, 1.949648438817852}}], 
      LineBox[{{15.52807047584685, 2.4243668761537593`}, {15.740091454441195`,
        2.4360684247050886`}}], 
      LineBox[{{18.492341839908118`, 2.6142740417093675`}, {
       18.821449306182032`, 2.6395162373133374`}}], 
      LineBox[{{19.94153124592807, 2.733470483245029}, {20.169965455101522`, 
       2.754330941876069}, {20.516695809373633`, 2.7872225167896643`}}], 
      LineBox[{{23.226047178497186`, 3.1094414545487137`}, {23.64920156591676,
        3.17340116189556}}], 
      LineBox[{{22.30769504889355, 2.985091982002679}, {22.552358314022023`, 
       3.0164564905150737`}, {22.667860135808915`, 3.031688294707287}}], 
      LineBox[{{15.343333012503038`, 2.4143697715089045`}, {
       15.432407210540728`, 2.419167306163315}}], 
      LineBox[{{18.278595664705986`, 2.5983979935088914`}, {18.36672775876869,
        2.6048959453198033`}, {18.396678574601996`, 2.607119481454969}}], 
      LineBox[{{16.827635391221843`, 2.500232297983178}, {16.959795951595527`,
        2.5085393163924277`}}], 
      LineBox[{{24.641395340750787`, 3.343565604471415}, {24.655365320941417`,
        3.3461915693103883`}}], 
      LineBox[{{10.686032608117792`, 2.210818688735391}, {10.695435342311011`,
        2.2111542530000747`}}], 
      LineBox[{{23.043635647936327`, 3.083246216359192}, {23.05083928655206, 
       3.084265758591256}}], 
      LineBox[{{15.25449064936491, 2.409626312171736}, {15.25449064936491, 
       2.409626312171736}}], 
      LineBox[{{19.8193359030368, 2.7225634976773403`}, {19.845867980621946`, 
       2.7249171149694478`}}]},
     Annotation[#, "Charting`Private`Tag$48818#1"]& ], {}}, {}},
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 
  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
478 
  PlotRange->{{0, 150}, {-4.790211036257171, 7.236631161024249}},
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 
  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}, {
himyss's avatar
radii normed
9c694860
 
himyss committed
495 496 497 498 499 
   3.879569032952899*^9, 3.8795690542179728`*^9}, {3.88008674094413*^9, 
   3.880086745707924*^9}, 3.880087709857608*^9, {3.880087803643498*^9, 
   3.880087992161138*^9}},
 CellLabel->
  "Out[231]=",ExpressionUUID->"3ccd7abc-262b-44b5-b38b-4667bce6baf7"]
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
500 501 502 503 504 
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
himyss's avatar
radii normed
9c694860
 
himyss committed
505 506 507 
 RowBox[{"myU", "=", 
  RowBox[{"U", "/.", 
   RowBox[{"FindRoot", "[",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
508 
    RowBox[{
himyss's avatar
radii normed
9c694860
 
himyss committed
509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 
     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", ",", "60"}], "}"}]}], "]"}]}]}]], "Input",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
545 546 547 548 549 550 
 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,
himyss's avatar
radii normed
9c694860
 
himyss committed
551 552 553 554 
   3.878289091001782*^9}, {3.878289430122672*^9, 3.878289446396935*^9}, 
   3.880086673187624*^9, 3.8800877134548693`*^9, 3.880087807038721*^9},
 CellLabel->
  "In[232]:=",ExpressionUUID->"22f1491b-05f8-4461-9d8a-7023940b0eba"],
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
555
himyss's avatar
radii normed
9c694860
 
himyss committed
556 
Cell[BoxData["48.01652525001887`"], "Output",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
557 558 559 560 561 562 563 
 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,
himyss's avatar
radii normed
9c694860
 
himyss committed
564 565 566 567 568 
   3.878292505137838*^9, 3.8795690844906807`*^9, 3.879570642631118*^9, {
   3.880086743953726*^9, 3.8800867473357487`*^9}, 3.880087715279483*^9, {
   3.880087808333634*^9, 3.880087992954145*^9}},
 CellLabel->
  "Out[232]=",ExpressionUUID->"b4205776-3774-42a9-8a33-cf7915e0b844"]
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
569 570 571 572 573 
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
himyss's avatar
radii normed
9c694860
 
himyss committed
574 575 576 
 RowBox[{"myCoeff", "=", 
  RowBox[{"coeff", "/.", 
   RowBox[{"FindRoot", "[",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
577 
    RowBox[{
himyss's avatar
radii normed
9c694860
 
himyss committed
578 
     RowBox[{
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
579 
      RowBox[{
himyss's avatar
radii normed
9c694860
 
himyss committed
580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 
       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
598 
 CellChangeTimes->{{3.878289683612082*^9, 3.8782897289455433`*^9}, {
himyss's avatar
radii normed
9c694860
 
himyss committed
599 600 601 602 
   3.878289801822311*^9, 3.87828982476785*^9}, 3.879569118948894*^9, {
   3.8800866773952007`*^9, 3.880086704667201*^9}},
 CellLabel->
  "In[233]:=",ExpressionUUID->"444da03a-39cf-4c3d-bb1b-9d322299f630"],
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
603
himyss's avatar
radii normed
9c694860
 
himyss committed
604 605 606 607 608 609 
Cell[BoxData["8.088071130275875`"], "Output",
 CellChangeTimes->{
  3.879570442808729*^9, 3.87957064413636*^9, 3.8800867486475277`*^9, 
   3.880087717563438*^9, {3.880087810097796*^9, 3.880087994342682*^9}},
 CellLabel->
  "Out[233]=",ExpressionUUID->"e978557d-baa8-4635-867c-ce6bf4ad2127"]
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
610 611 612 613 614 
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
himyss's avatar
radii normed
9c694860
 
himyss committed
615 616 617 
 RowBox[{"myNorm", "=", 
  RowBox[{"A", "/.", 
   RowBox[{"FindRoot", "[",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
618 619 
    RowBox[{
     RowBox[{
himyss's avatar
radii normed
9c694860
 
himyss committed
620 
      RowBox[{
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
621 
       RowBox[{
himyss's avatar
radii normed
9c694860
 
himyss committed
622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 
        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
637 638 
       RowBox[{
        SuperscriptBox[
himyss's avatar
radii normed
9c694860
 
himyss committed
639 640 641 642 643 644 645 646 647 648 649 650 651 652 653 654 
         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.2"}], "}"}]}], "]"}]}]}]], "Input",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
655 656 657 658 659 
 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},
himyss's avatar
radii normed
9c694860
 
himyss committed
660 661 662 663 
   3.879569335601615*^9, {3.879570630856184*^9, 3.879570632038727*^9}, {
   3.8800866795052*^9, 3.880086723648281*^9}},
 CellLabel->
  "In[234]:=",ExpressionUUID->"7835cc6b-0406-456f-a3a3-ab6f89b509c7"],
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
664
himyss's avatar
radii normed
9c694860
 
himyss committed
665 
Cell[BoxData["0.7696211451653487`"], "Output",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
666 667 
 CellChangeTimes->{
  3.879569148123713*^9, 3.8795692505602217`*^9, 3.879569281919847*^9, {
himyss's avatar
radii normed
9c694860
 
himyss committed
668 669 670 671 
   3.879570637910142*^9, 3.879570646893675*^9}, 3.8800867527823553`*^9, 
   3.880087721464933*^9, {3.880087811800939*^9, 3.8800879962360888`*^9}},
 CellLabel->
  "Out[234]=",ExpressionUUID->"46e6a3ad-c2c6-4495-b04b-21996d66b414"]
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
672 673 674 675 676 677 678 679 680 681 682 683 
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"Piecewise", "[", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
himyss's avatar
radii normed
9c694860
 
himyss committed
684 
        RowBox[{"myNorm", " ",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
685 686 687 688 689 690 691 
         RowBox[{"fIn", "[", 
          RowBox[{
           FractionBox[
            SqrtBox[
             RowBox[{"2", " ", "mass", " ", 
              RowBox[{"(", 
               RowBox[{
himyss's avatar
radii normed
9c694860
 
himyss committed
692 693 
                RowBox[{"-", "Esep"}], "+", "myU"}], ")"}]}]], "p"], ",", "r",
            ",", "0"}], "]"}]}], ",",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
694 695 696 
        RowBox[{"r", "<", "range"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{
himyss's avatar
radii normed
9c694860
 
himyss committed
697 
        RowBox[{"myNorm", " ", "myCoeff", " ",
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
698 699 700 701 702 703 704 705 706 707 
         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
radii normed
9c694860
 
himyss committed
708 709 710 711 
   3.879570725415313*^9, 3.87957082556467*^9}, {3.8800866815822*^9, 
   3.88008672825279*^9}, 3.880086759707464*^9},
 CellLabel->
  "In[235]:=",ExpressionUUID->"fc251146-842e-4797-bd33-46971b36a427"],
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
712 713 714 715 716 717 

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
718 719 720 721 722 723 724 725 726 727 728 729 730 731 732 733 734 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 797 798 799 800 801 802 803 804 805 806 807 808 809 810 811 812 813 814 815 816 817 818 819 820 821 822 823 824 825 826 827 828 829 830 831 832 833 834 835 836 837 838 839 840 841 842 843 844 845 846 847 848 849 850 851 852 853 854 855 856 857 858 859 860 861 862 863 864 865 866 867 868 869 870 871 872 
1:eJwVlnc4Fd4fx+3M7L0uV8gOyag+x6qkbCozI3tHIdIXWaWQrWyVsgpllJXs
mRGi7HuJIuRexO/+/jrP6znPc57zeb8+5/McEQcfkxsUZGRkreRkZP9fn911
oDnh9ODstOZbkZOTUWfbhjGXrmDcoIIpRakd4wiH7Y2RGEwwjEa71KdhbsKz
5gRhakwchHI597ljIsHzTKDXjnAmyC0f46LCPIHhPdeUdeES2NPQaE0QKoQj
3+Yk1oTrwH6QyrKDpwryu3P9fgp3waFlzoUW8k9QXhH3zVNuEiY1pqQuOA2D
CwxlrvrPg9P27WqlqW9AVvC3/tuDeWheLT0SsvcNMqgEprqL5mEv705DBd8U
dHY5C5WMzcNNk45d/NUpkDTbL7ihvgAdol7DiSNTgHeVKP1OvggpFxvxxZ3T
4J4c1jiYuATFZSX5LoU/wHvx+Hx1+TJkZWkGND6ehRnbCdqWrmUYrJabkciZ
BdPxWLm+hWXw4wzZiyqdBdUeXPAi3wpkkTHw8HXNAuWbYhaumBWY3+f8Uk4+
BxmhImdv2/0EcYn48mDfOWhh581UY1kDvb5aChOdeVBO6GzUlV4DG8+IJCHj
eXhOE7RgfG4Nxuxyc6dt5uEhcUzePXQN4kTUuJRuz4PFj5T2LPwa9Dl2BqqX
zMPKK5atveZf4LOe6VRFvwD0IlESP1TXQXb/uwL7pwXoPJ0kf1d3HTzy1acc
+hYg+mrOKUGTdcC9tJku/roAFIm156081uEzC1Un9eoC7B6sunx9tg713D1J
MhyL8HPK/OUAxQYo16a+ZHFYhN4MiePNvRugoCrDQL+1CPHVyifsJjbAydJy
tXZ/ES4MaqodLG6AxJkCd0vqJWg7Yq13+nAD7pdEHA/kWoIPtxLdak/8gbTX
LsJ0qktQZkYsqUz7A2o+rWHPg5bgEUuPVP71TWAP7T+f+nsJ6k4J1hh4b5I8
AiN2ZwkWbH1g/84mDAgvOeUdLIFGGbvZlfRN0PEP/OjGhAO8nnUYU/8mRF9L
CFKRwoF25NpAsMYWDDTuFVVdx4H3K7CU0NuCE/FEbVUXHGQNJS2MWGzBxzxa
iVIvHKwLqxDl/bdgkPbPc/sQHOR8uCu69HIL4t9XM1M+wQFhmznQhGsb/Inx
9fwtOMAKOhweYrfhs+iZB8ntODDQqY4rO7ENDx2+WBF6cFCcfDWH9vI2EDBa
a4ljODCRz+9ojNwGjIuY6/YKDu5abBp7Jm1DIGNs//xvHJSE6U7x5m4D+Wtt
7pYtHJD1Lq8H1G/DecdFzQsHOChzVeST3tiGs/w3n8Wy4GH8cVTR+ME25OvN
vl/jwAPV+zG5aMa/sBkvnaXJiwdL6jvasxJ/wW3hsKRGBA+0BZ88023/wtNa
hsL+E3hw+GbWRNn7F7Rar0a3GOCBPnvWy2/iLyi87gukMMHDW0tvgR9LfyE4
8EqBojkeKCdjguvJd4Ds9M+S61Z4KBqvV/I7tQNrhzdUN53xcCnj/Nx3nR1Q
eRw5WuaGh60rI4mXTHZgcvVG7hVPPOh+XVuT8NoBrJjs2yA/PCyNYl5+L9iB
wVmNwvUQPDxKLbO4VLkDtULFgWJheFAxV6eu/7gDPLJuLufD8RAzYuqQOr4D
ASnHCy0j8aCQMsNCubQDVusWKxfvk/Iw9Wry3dyBxGOa+lIxeJAcjha4dJQA
cmpy5u/i8TCYzN5Tx0+AK+eb9hwe4iHIJC9Y4jgBmBKPvz9IwEPXUN0YhQ4B
pr3EXSmT8OCfdO6+rzEBHoUUWXkm44HPeFjpuy0BjN5+sWt7ggf3wdXEumAC
TP5OS9NMwwNbYjBIxBBgpSewwzkdD/WGNL9SUgjQeS+N+k4GKW/mJ08pCghw
tFTd6G4mHhgGhPV9KwjAllXx3DcLD9WPSonTHwjges6A1iQbD9YGai/1uwlw
2fpOgOhTkt+j7RZ1XwnAeerFyiyJS/tMqCUWCZA4d9k98RkezBJ+VKX8IYA1
/84f2Rw87F/ydKAgI4K+oUpUPYmLGIksvkxEyHj2SlAlFw+Xe+83TfMR4dTX
g8Y8Ev99wOatL0kExvBll10S5+jnCtSdJAJvtwyPdh4ezjHI9IhrE6EmN20g
hMS/umuDU4yIUJxJ9iifxGnxupIUtkRQpjIyfU/isxe/jPl4EMHwe6DwBxIv
0dndnw4iQsKE958KEj/u+qmkH00EAUvt3mQSq8YFzdU+IcL8rSNlTiSeuUCd
JJ5PhKvEyZRjJI6lTYaUciJ42CxHjpLup9Ap9Iv8AxEChtxCAkg8HvP6qU8X
EWhVcoMoSHzvvKr+9BgRWgXH74aT6pc88pl4cYEIIn3eCT9JeQ22G7+s3SCC
onpX0TkSB0V/txA/JIIV47XPiaR8Mec8qFMYd+EwO/t3F8mH3+coBx+JXQju
dHWkI/kSWjhH81dpF06GxB1gSD57KOhehaJdkLodNnGc5PsYStiIu7YLjPuq
xsyk/vhia5jK4rwLYQHav3+m4CE8jFUt3X8XbvK8GKgj9dNEfeq9oge7QENp
b3mc1H/RE1fEpDN2IXAD0fQ8xoMSgbfzTdEuqFgm7Ng8wkPCyRzmpo+7kMk5
U3blAR40K5/nTPzeBdyf+MisKJKfflet6/u7EJtVHtAegYfsNamlJdo9EJZi
qpq/h4dtqQrZLZE9sOqiE/gTioeS4vcfjprtgVb3aLNtAB4s2oKvp17fA8sG
lzOc/qT3P69BJeC1Bxycp1g/+uDBBtOsfzx6D2ZdVdpm3PHAmt05qV27B3s6
Z/8Luo6Hxrr4u91te/Ath2q2zgYPHuOXRI2H9mAl4EzRqiUePnMOudmu7IFg
S4eVOGm+BCdOEIIE9iFMx6955QIexCuynx5K7gN5QrPne108DPfZoOiT+/AS
i48P1MKDLMNczBODfRiUzBRr18DD7P0VrvJ7+6Bjdo2mTpY0T4pK65UT9iGv
jCqEVQoPGp+8bRsy9+HvBwFfK3E8pJJtFne+3Yek8IjVT0J40Luzpzy/sA/3
JNOdcUfx8MaX0ZRH7x/QSS+Gh6zioDeLPbjM4h98YMuVjcPhYKmNL1fL6R8k
s7uTxczhgI/3+E/P8H9w+oEd1mgcB1EtupEtNf/g4FdqG88nHFiwhVe5ix7A
xBvfruhUHPiejp44lD8AXu/VuLjHOIh3TjhMOXMAWuX3y+7E4aCpPlu/6eoB
/KiWl5cLw4GkY+08++MD8LbqSut1xMFu1Qb7x70DiFmLW/8jiwOO7wQ1E7pD
wCsV7q+J40COluw6jusQkm8/rRkXxoGj9dFSVsVDyK1Q1gpjJdVDJa3j7HoI
RfvPzSi3liDH3CmAefQQhm1uEo+8W4L8Ox4rClfIkEaqoRef3BIoNNMRYmzI
kKfpraiKY0vQTPWS+ocjGQrXDjtUEVyCmYQFTIIvGRq+SaMiyLgEQnk2V5bj
SCwyonp+eRGy2gzb8j+Qoc5/NIZc+YuQzKScwyZKjm5891HvpFmEiNx9461V
chT1YMp86PU8GCnU+klvkiOvCpHLrnnzINRyM8mBSI4sKXjLNlPmoX5uZXCQ
mgItDii34sLm4c+xCcMyIQpk+e002TnDeXAorbnsbESBHONyJmd+z4Fmnbfe
eBUFumtD/N4iOQdkwzPoQwgl4jxm7hF2bwailhvzQ/6jRLX/3Dsu+M8AHdkz
CrVYSvTpAo8txmkG2GSvtdWkUiIqb3HWrfMzcCxm6HxFJSUK0I1xk2GZgYsa
rQYFS5SoPaaxtkL3BzwpLLSONaFClp3bDc/zp0E8wDnITIoa+aTf2IrLmgRZ
9R+muQrUiJc8eEjq/iQok12VX1GhJuX9dmLKZxK0HuothWtTowb1IOUA3Umw
LZIxe21NjYJ/f33ivTEBaSN/5CkeUyO694Ki0ZcmgPpkOK5ikxp9JPfdmD86
Dgvb6RYMH2nQVVvdb+kto1CgqXRz+hMNsjAof/K8fBTsEvofV3TTIPV/Kt2f
s0dhUoy6y3ScBml1agedvDUKQ6b+Gk83aVCQ9EW9celRaHxzCSMrdQQFCwX4
CGeNQLoX+YpB+hFkliLRZRoxDBcX3cOS/WhRFA3bk5C7Q3AuDFuqE0SLxujC
QqJ9hkCLc2ry711atCDOJpV7fQjUdS+rWj2kRRriLO/+ag2BVLH8JvYFLcpv
qjc+QzsEDDe2XGumaNE/lgNV4dRB6F0IMxs/T4eyL8+ZTfgMgMFCorSQED3i
d9M5SQm90Cu9MGItRo/ozr5TuivZC/o3T93NlqJH9ZZ24gxsvXCB4vsgzyl6
1H9RpMJysQe0MNK32I3oEdMp5+b+hz1w0vpzM20EPYo12Us7+N4NAiNE881F
evR1NotMMLELVj453OssY0BK3+wZbot2gHDSuPujagbU1CATeutoB5jZGpib
NTAgHy4l1ZjddmgkqEnNdDKgLhoFi6kv7ZAsyzqyM8+A6H55xptGtYNGepOE
BB8j0qYfI3db+QwP3QUG7kczoh8PB6R6W9tAjnVMUNuOCZm12BaE5LfCUsJS
ROoNJuQeLID8klvhGf0ODufBhLgqCskCIluBkYrn7cMgJiT8Bt/67EYrrPy9
dm4siQkJyP3T+E+qFYqnpr3d2pjQC3L39PB3LSDwcrHp8fGjiI2s+wzP12ag
h237qc2jqK68PjLGrxGCNCVp3+0eRcYX88rwFo2wpG1V/picGb0ISxKyP90I
rRdadrWYmZFSypmo9CONEGKS8KREihk9lJY7zpP7EX7eOPb5lj0zUqRLdL7y
5QP0PjSXZB1gRhe4JHN09Rrg/Jrljqg5CzL6JRsxFlALORyRw+FXWRCj2VGX
r3a1sKXxumLKigWNfnlqsnKxFvLi91zSHFjQEPPKoa5ILexKPB2n82VBfp+U
qYf730Opw3TdejwLIsSRB6rIvgeWcbvQxmYWdJviq23D3xoYb3Uit5RhRV8S
mzxL6qpgXXOSL0yeFdUItiZPlVTBkRZD5TxFVlSd1nxdIKsKVJrUXZZUWVGU
qUxvS0gVpDSw9PnrsqJXIbPtBaerwLD6Q8YDW1Yk2SyaXt36FtqLORQ+JrKi
46dqsoa/voGa2DZbzDYrilyl/riDrYSXF3t4YgisKDdv98MmZyVkMX75srbH
iqrCUvB7RyrhXuKPcw0UbEhoc1vBfryC9J/elbNgYUNekiPuZu4VMPdC4fCB
NBtSL/F/+Sq5HI52PM3dsWdD0oe7m7FbpeBMHTjT38+G5rMGRyYIJcAn3PV1
b4gN0cKrsu75EuhXFRyQHGVDjDdmMnv7S0DF8/PHiG9s6PYDhkDK4hKgGebM
VsGzIeXzfjUaJiVQnPfOPIecHYlZBtzprngJCxqEbi9lduRnE3KW6vYLcPS/
U8OYxY4Mei2IWKViUFH8XnP8GTtqw9S90xMtBvo/6N25PHY0r/dmIpS1GCr9
aN7fe86OwvS+HGNYL4J936Tarbfs6B6urry9rAjSfF40TPWwo2KRAptK6SLo
8hxuKf3Hjuz5dN/SnCgEeRfp/kv2HEjr4VXZZ3b58N1N2KfMiQM506tyFRjl
w0NPdpajrhwobYZypUozH5b99owHvDkQXnnq2AE2HwpDe0aNQzlQwDJtnfRy
HnAne0xbpHOgp0EKzbOBeXD44fXq9T4ONCYmmjutkAsDbDIMgWqciPztXfHS
w2w46+W1nH2aEwUtOX/FLmdDeUd5RytwosdVAa/KvmTDo9ATUSznOFGSX1no
TnE2XFpS+ffahBNZPJzQcLicDd11mr/nPDhRn4PZs5jcLPh8/cqwUQ4nEra8
SllplAkN5RHZMpRc6ORdwvqD+TSguWTOl07DhaYdXQZSB9LAeFkik5yeC1UG
OXO9aUgDHLYvbYyFCx3hds3nT0kD9kzu5HtCXGjWOY+6STcNPCJLY0fUuNAb
2U0KxdJU4L86Fhjqy4XkdyKNHO+nwB1yKaPeaS70vMjZ9feNZFjSYRgXn+VC
d6H/vZxFMhjFrtr9t8CFxizoou6cSwYx5goflZ9cKP7eF6czEsnQK6D8OI/A
hYK122jsl5NASPVsfwA7N9oJaFoQ9E2CVm+TS4J63Oj9oe4nTEwi0E+FXPCu
5ka7b5qm/r1JAJXbJUyP33OjV7s4hXNZCWDPNv6lop4bZQ4G38qPSIDaCydt
Npq50TVP+8xE0wRwrvntF9DHjYa3UzMN/z6ElkdO2SFL3EjWwvG+NjwkzSuD
X9E8PEjTVCX82fd4WHwukpIbyoMWDkSS2IxiYTVqlkh5jweVVhkZU2nEwh+H
fDvXSB4kTp5/mVY8FsiERaRPxJP2meMazu3HAF86prU1nQd9xxfVZJfEgFGM
8PriWx6UzFWXkXEkBj64CF6SWeZBxySIoQGD9yFFkpeqzpwXreQO2iclR4K1
rfiL+au86Dne+XRjRCSIpShdPGrNi1TNPOd2/COh+vByoqMDL2r+6eGSZBoJ
I2MRAkd9eNGoe9FJa85IYL+/quwYy4vuaWqRX8iOgOTZphtMDbwoicF2XbHy
P0jMcu6wx/AhsbiOY3GEcLh9K+JmKJYPlT69hjvxIxxsTXKE08X5EDFssvHX
53CQoR+73SvDh7jGaXyLnoRDV7Cu5Ck1PnTmuQ3mQD4cKK6JxTGa8KEpt3uM
Np53IZB77uL7SD4U/YRGwXYzFCxTbPoZcXyoZtRf4sypEBCoouf3WCGd73C9
6Q0mBH4MvXfpWuNDp/vC6tXoQ8CJmY08eosP2f0tTo+cDgbv+HalQ3J+pEdu
LfQvKhgi7slnbgjwI2YhXdbgr0FQ4knuNGbCj/KTJCmPxd8Ggk4xMbeRH9Vy
6dvtiwZCaENuZFsLP8IHJVrQsAUCmWIW03IbP4p04V4RJA+EI8KPRRR7+NHp
Da36wB8BwEUM0mv7yo86j7sNtD8NAKWyS5n4dX7kq8/tucgbAN6cW6onsALo
+0LEsR7BmzC/qBXUGiOAHjS0cNBd9AOt3P4CqQcCiMlP3uFQ3Q/yrlr2JT8S
QDFR1G0UMn5g2+Mn6pgqgB5x30xXOOoHk5V5vZSFAkjRS0VTYtgXvtw5wOg2
CqB/uzd/3LLzhVbW+q6OLQF09edvt+cUPlBw5gR/n70gcrdsib0z4AGtBvSR
e06CqKEsI2mh1gPm7OZXjrsKohHZrk27Ag/ARqTWR3sLIqvKwtiEQA8oaide
Q3cEEcfNBfNWAQ8oNmzNqEoRRBjH7bf3vd3hpb0pV1aHILq8Y/Yfs7AblN0P
ZHWRFULmUV41ym+cYfW0TcR9BSHk2JfrLJrrDDJbOpuFSkJojYo3VCjBGV45
cIzOqAkhe37W55puzvASqjIszwkhtUCP94qizlBEXBc2sBNCRj+2dmrSbkCW
l6fcySQh1EhfQb4V5wQx5o76lNtCiPGN0yuKdAcoMJyyvkEQQjw6tme3oh3g
g565d8eeEDqkTa7cvOUAf86cT3pAIYxo9jsP+a84gLW49Fd2FmFUWJu9z8Lj
AIo7fxzEpIVRJjZaKuKpPXzPiAjRtRdG7Ee91E5UXIeTU/klMX3CKGSrS6/l
jw2wyLkSKgaFkZ7frytPJ21gNVzu/PiwMMpQk5yOb7WBImzDguSkMKqaEj1a
kGwDHB4jwt1LwqjzAG/+XskGdkMp9WnIMKiUYEbpG2wNEwwqLRd4MUiwxTGh
ks0KMo5nl/XpY1CT+PeWi8FXYWvzXHHkZQwa7k7sWXa8CoaNf56qGWKQcAlP
RLbBVaAxvfiw2ASDFJ0Dq1XErkJgKNE97BoGkXnyu10cvAImA1ckZV0wKPjv
NueE9BVgDOQofBiBQf9lzniLr5tDeHNChn4dBt3Mv/JT+5UprEw3pxDrMchm
bdrYOdMULPY2E198wKBB7mvy6bGmIKNiGUfZjEF0npmCYi6mMP5aPKShHYPC
n0Z6L4uZgmJ6k5X0CAad7Lx1lrvQBJa8NoQYfmPQ8jHutNISY7jMb17cLSaC
RkWiNh7WG0LGglyzlbgIcttAU8fzDWGhjPbbqoQIEonm5v4WYwih6AMLs7QI
CqMQOB1ibgilN7ChpidEkLxFTnjJhgHQV26YTJ0RQUvy0TsfZQ2gQ+cR2S8L
ERSz/DHduPoSaHp3WLHGiaBEXG/T3I4e7I4vn2uLF0GPTat+W83qQbU2o+Lt
hyJocFXgFb5bD8R5jWmnH4ugONP6kNM5esDYNln9Mk0EcZpyV/fo6sE47y8m
VCSC5kPfHQRnXADfzxzN3k0i6OVQbXqR/nnIFXAQ690m1TN388Xwog6oK0nL
OeyIoA89yfFLIzoword1ikAQQRe6jbBH2nSA9na0vti+CCq/e0HPr0AHfIde
+YdSiKJffP6KeDsdQNGbzTLMoih9SmLDYlob5n5H2SRIiqJSzpGitVktEPv0
Is3AShS9kxr5r4NZE/plCEVT1qKoOVYXNCg0ISjtQpW7rSja+lDdWLeFoMdt
eSDaXhQdFXedHZtA4M8iTdfkIopWeXM6E4sQNNmU35EPEEWvork79zQQWBKq
r7MkiCKXVeUbKwtnIVG2VepLoyjqXqm0CSdowEtVK7qIZlG0Js9h6TKjAU3a
W7gTraLII+Uxt1OnBvy6Jl6c9FkU9ftbPcrO0AD96Hhh415RdIs6NIRMXQNo
fphwDE6Iov0mC0LyPXW4kzj/r29LFM0fTP1xEFAD502qoS4pLDJd+vgX80AF
bF/p9FvLYNGev1gJebAKWNhH9fyWxSIujoPVdWcVODdA2c55AouqJzR2d7VU
QPw1RYP9KSwqWT6/2LB3EhYdyIqJ2liUo2Xye8nnJDh92Qs6boNFdg844YmD
Mji8+YOJe4xFBB359VVvRfCMMN2/mIRFC1Ix7I62ihBoWv2V8QkWvXoZw4S/
rAix24GPEtNI52ddNhKSVYRyNeJe+jMsklLiylVfPQHElsOx56+x6AvjPaK4
1wlI/MKY0NaBRbdPpn4MDFaArEIv1+guLBr10u8tc1eAwoB+7Qs9WMRw8UBw
00oB3nEl7vb0YxHNo17L4rMKMGXJ7jo8ikWh9v/9fk6lABLzvNpz81g08Jvi
haS4PChUhwgVLWIRu8EjKRV6eVC//414A4dFhurSLda/5OCSxLPK5RUsmlyv
7Nt6Jwf+HhihjQ0s+liw9eeynhyEnv6P+HYTi57vRD8TkZOD+0xzIwHbWKQu
6FDCxC4HGRWFDwgELGqMjUmQm5aFgntULvW7WPRJ7LK6Y6ssvDa+oRW6j0VX
+Ew2Kl7IQrVou+DZAyw61kb1jCNBFho3xYmHh1hE8XVJIslfFv4HFzk/YQ==
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
873 
       "]]},
himyss's avatar
radii normed
9c694860
 
himyss committed
874 
     Annotation[#, "Charting`Private`Tag$51114#1"]& ]}, {}},
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
875 876 877 878 879 880 881 882 883 884 885 886 887 888 889 890 891 892 893 894 895 896 897 898 899 900 901 902 903 904 905 
  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
906 
  PlotRange->{{0, 10}, {0., 0.7696209539853472}},
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
907 908 909 910 911 912 913 
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
himyss's avatar
radii normed
9c694860
 
himyss committed
914 915 916 917 918 919 920 921 922 923 924 925 926 927 928 929 930 931 932 933 934 935 936 937 938 939 940 941 942 943 944 945 946 947 948 949 950 951 952 953 954 955 956 957 958 959 960 961 962 963 964 965 966 967 968 969 970 971 972 
 CellChangeTimes->{{3.879570843622718*^9, 3.879570853774201*^9}, {
   3.880086754549075*^9, 3.8800867612813787`*^9}, 3.880087724555195*^9, {
   3.880087813143466*^9, 3.880087997646585*^9}},
 CellLabel->
  "Out[235]=",ExpressionUUID->"7832bb28-3d46-47d5-a351-fd05f230e3ed"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
    "R", " ", "\:0434\:043b\:044f", " ", "\:0433\:0435\:043b\:0438\:044f8"}], 
    " ", "==", " ", 
    RowBox[{"1.688", " ", "\:0444\:043c"}]}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  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.880087741923903*^9, 3.88008775501427*^9}},
 CellLabel->
  "In[236]:=",ExpressionUUID->"3894cc59-ff84-4f88-8130-0b3a87025770"],

Cell[BoxData["1.6882859930125582`"], "Output",
 CellChangeTimes->{
  3.88008775773094*^9, {3.880087791981393*^9, 3.8800879996867847`*^9}},
 CellLabel->
  "Out[236]=",ExpressionUUID->"02d33532-4691-458b-888b-cfeaf1cd9444"]
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
973 974 975 976 977 978 979 980 981 982 983 984 985 986 987 988 989 990 991 
}, 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->"64e28415-be8d-47eb-8c2a-84b3776a3911"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
himyss's avatar
radii normed
9c694860
 
himyss committed
992 993 994 995 996 997 
Cell[558, 20, 1896, 49, 154, "Input",ExpressionUUID->"ad6df89e-d3d8-490c-ae91-dac6d80d1b5a"],
Cell[2457, 71, 2806, 75, 179, "Input",ExpressionUUID->"2fc4ef7c-d275-441b-9d3d-7d7239f74130"],
Cell[CellGroupData[{
Cell[5288, 150, 2232, 54, 75, "Input",ExpressionUUID->"c2ad0be7-22ee-4300-8ad8-45a6e113c1b0"],
Cell[7523, 206, 16900, 292, 235, "Output",ExpressionUUID->"3ccd7abc-262b-44b5-b38b-4667bce6baf7"]
}, Open  ]],
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
998 
Cell[CellGroupData[{
himyss's avatar
radii normed
9c694860
 
himyss committed
999 1000 
Cell[24460, 503, 1906, 50, 75, "Input",ExpressionUUID->"22f1491b-05f8-4461-9d8a-7023940b0eba"],
Cell[26369, 555, 793, 12, 33, "Output",ExpressionUUID->"b4205776-3774-42a9-8a33-cf7915e0b844"]
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
1001 1002 
}, Open  ]],
Cell[CellGroupData[{
himyss's avatar
radii normed
9c694860
 
himyss committed
1003 1004 
Cell[27199, 572, 1046, 29, 53, "Input",ExpressionUUID->"444da03a-39cf-4c3d-bb1b-9d322299f630"],
Cell[28248, 603, 290, 5, 33, "Output",ExpressionUUID->"e978557d-baa8-4635-867c-ce6bf4ad2127"]
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
1005 1006 
}, Open  ]],
Cell[CellGroupData[{
himyss's avatar
radii normed
9c694860
 
himyss committed
1007 1008 
Cell[28575, 613, 1907, 49, 53, "Input",ExpressionUUID->"7835cc6b-0406-456f-a3a3-ab6f89b509c7"],
Cell[30485, 664, 368, 6, 33, "Output",ExpressionUUID->"46e6a3ad-c2c6-4495-b04b-21996d66b414"]
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
1009 1010 
}, Open  ]],
Cell[CellGroupData[{
himyss's avatar
radii normed
9c694860
 
himyss committed
1011 1012 
Cell[30890, 675, 1247, 35, 53, "Input",ExpressionUUID->"fc251146-842e-4797-bd33-46971b36a427"],
Cell[32140, 712, 11288, 205, 239, "Output",ExpressionUUID->"7832bb28-3d46-47d5-a351-fd05f230e3ed"]
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
1013 1014 
}, Open  ]],
Cell[CellGroupData[{
himyss's avatar
radii normed
9c694860
 
himyss committed
1015 1016 
Cell[43465, 922, 1391, 43, 75, "Input",ExpressionUUID->"3894cc59-ff84-4f88-8130-0b3a87025770"],
Cell[44859, 967, 221, 4, 33, "Output",ExpressionUUID->"02d33532-4691-458b-888b-cfeaf1cd9444"]
himyss's avatar
wf in systems
1d727bf3
 
himyss committed
1017 1018 1019 1020 1021 
}, Open  ]]
}
]
*)