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

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

(* CreatedBy='Mathematica 12.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     62484,       1433]
NotebookOptionsPosition[     59044,       1368]
NotebookOutlinePosition[     59472,       1385]
CellTagsIndexPosition[     59429,       1382]
WindowFrame->Normal*)

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

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

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{
    FractionBox[
     RowBox[{"fIn", "[", 
      RowBox[{
       FractionBox[
        SqrtBox[
         RowBox[{"2", " ", "mass", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "Esep"}], "+", "U"}], ")"}]}]], "p"], ",", "range", 
       ",", "0"}], "]"}], 
     RowBox[{"dfIn", "[", 
      RowBox[{
       FractionBox[
        SqrtBox[
         RowBox[{"2", " ", "mass", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "Esep"}], "+", "U"}], ")"}]}]], "p"], ",", "0"}], 
      "]"}]], "-", 
    FractionBox[
     RowBox[{"fOut", "[", 
      RowBox[{
       FractionBox[
        SqrtBox[
         RowBox[{"2", " ", "mass", " ", "Esep"}]], "p"], ",", "range", ",", 
       "0"}], "]"}], 
     RowBox[{"dfOut", "[", 
      RowBox[{
       FractionBox[
        SqrtBox[
         RowBox[{"2", " ", "mass", " ", "Esep"}]], "p"], ",", "0"}], "]"}]]}],
    ",", 
   RowBox[{"{", 
    RowBox[{"U", ",", "0", ",", "100"}], "}"}]}], "]"}]], "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}, {3.88053533924759*^9, 
   3.880535339536429*^9}},
 CellLabel->"In[14]:=",ExpressionUUID->"c2ad0be7-22ee-4300-8ad8-45a6e113c1b0"],

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    TagBox[
     {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[CompressedData["
1:eJwVymk4FIoCh/E5Ci2WMZulsY1thkkzGNss/iUlZSpJGoUKMyUZUqZljjok
qUhcTWWp0HYKnXjqUMpSosiNtKgI2VIuSctJTvfD+7xffuYbo/3C1QgEwoFf
/f+M9MmCPJNq0c7CSbdxng5uzGCZRsnyPUov5RJHXHRQrSguZMpKPC72f60d
ctPBfU6MVEN2x4P50K62V6iDB0eeeH+WPvYIGLvs2OOlgwb3xL+HpZ0eXxOn
QmYF6KAxoyNmQDri8Twz1lAjXgc+tzXFfdIpj752zcT2Sh00aTvPnL6eiJR3
Xj6ZIbpYkO5198JyAzixFJ61GkSEz/uuzdU2Q3VzWSVFRoQvK3yUttgSo9kv
ytd2E/FXm3JHgysT00VVKfYyPVw8ELaxewkTxdpNswq36CHPeenyfyRMhBB2
B8yJ0kOqyoA1V8lESku6n26sHjatu/7qWDUTyqcKqvZePVB6BzwDfVi4ZRmf
tuGYHnZ+WkXqD7LFp6d/Oigr9MAn2ZZMT2Bj44/CK79RSWj+TdLQlM5GgMs6
wQd9EkJHD/Vkn2HjZ8xcgw4jEpJaBvVZNWzcMNW7c9vsl0+7tM9XfS7qLXfp
57N/+dksv+zDcxEWHL6tbeEvP5352UZlD57P7EMm8SRI/zmeWxTIAd+rmzvQ
Q0KWeLh2VygHAd/vp83vJ+HuWa8hXxkHORrF53OHSDDw/uo8Ec/BovoflOBR
Eh5krW9ZdIKDgcfMpYQpElj2zGkDTzk4ONLDbjIi433obZntKi5U5oo95/zJ
MCijpf+UcJHIqLKXBJKxUFNe/mQjFy0zxBTqOjJyrjLU9sRy8fvEcHLWJjKW
fT94+mEGF/6GzH0VsWRczfRr3trCBSeJdvh8BhnR9/sdS3wdwPFZ3WvyhIy7
uBLJDnAAhrxX5LeTQayUF1wKdoBZLmGuZQcZpSXfSIXRDlBcFK7i95Ixppo9
pspwwNi5UN0zE2TsjJxX+sdTB1S3vgpYQqdgLzGeHbDeERGVzCJxJAWPUvlh
reGOsB5x0YqUU2CsTshZsc0Rm6d6Wo/soKDqW8rspQmOqLDzvdCVQAGh59SQ
x1lH+Eh/rmzMoiCxvOoi650jsqxbM67eoSBVomE1udUJQc0vb04YUiHiaE0M
xjmhuPVAXJQpFaPqpPr2vU5Y0H5oZNCSisBrxptLjzihdNDOe3weFVYzeFc3
XXFCXMSMlOWLqLhbvonX+N4JR0fHiphxVEwQq71ObOEheJl/35U2Ki733acl
xfIg+J9p88eXVKy/9ahfvpuHIkJ3jctbKuqkzw/5pPLQ3r9yWucHKjKqRpom
L/EgWWXunKRJA3urSUDYIA/Lwg22Owlo2FS/J4Irc0YgM0Dn8J80hAmsxRPR
zjBKVyUU/EVDxLUWXmW8M5ifMkPr/qZBdtpK3eugM+oYa4tNH9CwLfpxQeB5
Z+icO8USvaNht77F24ReZxBVL/OZ5vo4Lm0Mag51gd/7hz/szugj61WcZ6bM
Bd5Od6/FXtTHf1aa2gXKXcD4XD+zrlQfKn7c97cJLli97sbn5Gp95OmanPyc
54Lv9dGzj/Xo4/KNmHb6axdMq+hdW8MyQLWm4crIta5Q22cb4XvHAH0F9j4v
NrjC/5aZeKjBAFpY6Ll4iyv+GLDxO9FmgEBFNM9ijyuOfVi5gDFkgJGBe4Yd
Oa6YCt3fP0U1hHFjbLdPlyvqZE6ffsgNsevIo+22UjfYZ+wosucZIY/ZvVUV
7QZSGolotsAIdXVfwjUUbjDh1TjYLDcCcZIR2J3ihhZZx974zUa4ELVLcOqy
GySdb/yt843wVGytPuujG47fV7cd1ZkDrt6+7ME4d5hp2+fPVKOjYsYzZv1e
d/TOXGyqr0fHAoJ9ZWGSO97UyMdFZnT4j3R0hmS5IyL3fNkbER07m1xY7dfd
0V4+tiFCScetQyOVNWPuSHY7uef5FB1e+xeJz3x3h2qGM7GYaIxmRW6XUo0P
iTCj6Zy5MTqlS9XdSHzIGpNN33saQ23RBXEJl49HUa8+EVKN4a0W8va0nA/K
wUJFFt0E//1Wvn2Xgg/5SfdlKfNMIBnV0lizjw/L+Gvrzy4wQWRXhS3pGB/u
mvi8cLMJ0qpocSklfIwZy3YE3zRB2+5mjR0f+Yg787huYbApQj8J2eKtAgit
vt4RR5hBm0RoaYsTwPpLjUPhdjNUcGq2BykFaJC06BrtNwNl26JKWZoA4SV7
+mQ5ZngwIPZJLBWAnr3z5Jk2M3DehMhujgsQEBP0Y3ipOV7/MNfymBRgS/Ep
XnKQOVLn9JbcUxei8EDnbu9Ic/QGyr48oQmh2dF00+ewOVSt8uQPrkJIj4pL
A5vMQWjYX8hQCqHeaRMpkDBwtd/T+9IBIWLMd3VXbGVAoqExPC9NCC2+96Zt
CQyUeaY6CvOFeNH1ga4oYkBWdbxmTbUQ/ROc15xxBlquF3QdURfhWUX7fJNc
CyifhCeRdURYPueeIvy6BWxHbZinaCJY2Z/rGmuwQJL9legLNiLQiwIS332x
gOulsp/VS0Q4fIR9Y8UaS5zNu0f/clQE3ffnWZlWVshZx362OlsEA62vq6Jg
BZVRZkZZnghDMcqcgiArpGWHasSWiPBxvGsWO9MKu9P+GRluEeHLdKeuHg1r
7Fi28fLSFyJ8zXG3a7e0hnxWQ9jltyLU2inVeZ7WiEjOfiEdE2Fj8Zva+fut
4fc7t7ab7AHl5NOmV2o28BWolPPpHrjFWpKqsLRBVli3ldzUA5Ju52o51wb/
Ao9CgpE=
       "]], LineBox[CompressedData["
1:eJwVj3s4lHkbx+dBZEyOScIma0O9ZEu0m/xueuZkhhB6qq0sOaSyOUdWI0Kx
WAkhVLKS7U0YSTwPSuggM4nKbqOIVFLkbLzz/nFf9/W5Pt/7vq7vGu/f3Hzl
aDRahmz+v+fac7ZEDW4HlmZQU4i/LpW0W7vplowTfqp5UeyjSzG2fmubfb8d
jHTsFOb0dCltaXXv2a/bIfFivM5E+krKNMFqulweh1zWRiY/QYfakbHJ5sMP
ONC7jjxvuKBNPQvWQpvW4SCdiIzd5qVNETvHWScscBBY5QclmWhT3iuqPFVs
cPCW21d56PZyKqJgY8R6Dg4RdzQKZ99oUQVlP9YcDsRBPaLfwIbQpNpe+LQV
BeHgeMH/bPQ6TWqCnv1KHIJDufWwI39Bg+IdncVsT+AwZl2RkFGiQc3+2OLM
SMWh4wDz4K4ZdWp3vftwxd84bOKMOYtC1aiEj4lzkkoc7AtPhk+ZqFE3DepU
tYU4GHTVVY33qVJLBQabf2/E4bWctqsPV5W6jb+Lc+rEQcwW4tT6ZdSKJ5F6
o2Myn7J7a42WCuWweM3C6BsOe20Fxu1iOhVk2WfvOYODycTZ0/FZdKo1EwIa
MSboHaMX9uvQqfBdyjVpmkxov9R3fPEHZeqZJM95gxUTpg5kj3XsU6L+HG+I
OxbJBJ0O+TapuQIVJMr49ZcTTBjP67kxhSlQ/Eofe85JJixR6M7mP5enFH9T
xgyTmGCgWjq2RCBPxYy4CzqzmeCjOx2Z2yNH+fV/iLUQMuHhMwMukY1Rtp26
MR8nmEAP/DvmUY6U1L3xcW/vNBMcVOwKsZ1ScjKV3HpvngmNM+XdNapSstLR
dy5fgQWrVLjLspMWyLWtN6N5y1ngt2WPv/rJeVKjgR1VvokFikdnfhXFz5Lv
ysMjAkJYQFMscDjfN0l6Gkj/LYtgQYYvueF28SR5Pz2R/T5a9k9LMYblO0mW
hOWsDIxngaXrhuriz99Ib7u6O4ezWNBRXqZmrvyN/LdrXhokZEGeb+jqBudx
sns6Pil8lgW50ckrzAzGSPwwY6xGyoItWufdCl58Jqv/ySIm5diQe/Oax6mo
z2Rm81WzSBU2FPS+XKFfP0q6pD54eFyfDRBpFr2f+Yl8tFpFPcaODfqXBWMf
w0bIe6zM3FOn2HAwaiYz3GaQHPSrfxydyAbLCqPq2Y8DpGLSgFzoWTb8kjab
P3h5gOQ8sD7qk8kGRiz95Hu1AfIR+5U9fpkNpv5izrnRN6SYY/xhSTMbDEOX
x7c0SMg3jsJtZzAOuDwyS3LY9oKUP/w6JG4JB3prAy0HGntJ45SlZVHKMp/y
tNPYvpf0e7hHM1CDA2V1WQ6x7B5yhLf4znGNjGd6uoP3dpNf+ZwMhj0HcjdO
iHyKukhsx8s36QIOjG0xiIwyvE+2NXp8H5bAgWM5HuG5vHtkhkWXD5HMAa8s
9ZLyyBbSULVtwDCDA8UdquqXRU0kPK55V1nEAcraotbLpIEU8P4cETfK/AQR
LIm7SdI4nK86CxwQuMTW4hqFqK22ZeM8jQsCrUmpyegllGGCQiUKXIBOYXPO
0xJkqGQ9UcbgAhXecIibWo6g9fvJn/S5UGzgm974+RYS4Isze7dyQaIc+MPl
+QZEsxdixVFcUFfezHJPbEd78KGO3t+54KUx/F8lmw5UzV6ZpXGKC4ZDTR/m
3nWgAOfotfFnuWB5YEfNdfwR6tprx/PLl+V7Ihe4U0/Q5YjWrP80cMHF+/V2
BwMRwiuem9TRHEHSWmrDNuxByTpTTqIkR1jlncYAQwnKamt725ziCCWZCro/
gwQVReVFVaU7Qtjbhvs6XhIkfLXtalaOI8TcJXNiiyRosChh3vMvR0hQqtss
1e9HDqaa11+1OoLV1mkY1HyD5reYLx1cwgNc8b4e88NbdGyPT/NUPA8ijHk3
TD2GkGrRiMQ2mQd3R8YXfw4eQhVvgxfjUnngVXNaZ/MfQ2joiMCWcZ4Hp3j1
K/rvDaH9Jy8KDUt5kP36bYPepmHEL+kp57bxwNbuCW6t8h6ZjvLO5TP4ILcz
5I+O8hHUH2d1EJ3nQyZ7A9f8709o7QrVtGMX+JD3nFlQ3PIJBZYP1V66yAe8
56NI7uUnNCHKU1Eo5QPNbu6nLMVRtNQYq2oT8iFZ+E7j0/5RZNH6GHPr5UN2
e1phEuMziqb7Fx3Uc4INa6aVxszH0PJzF16eueQE3ZvUgg/6f0Gl6uuiHuQ7
Q6Ft6aD4yDhqClnNNLm+A66l2R6xcPqGCg2L9be5u4CpSlu/WdckaoiZpQU8
dQHzK8+XVayeRmt4SvPf73IFnwGpUbTXDKruNEt58cwVJAJRWH/gLEptlVoa
7nMDq2vJ9BtH59C50PeBXn1uwO/wjDvBnkcx2t+pGe3eCY8hARtRXkDXzJZ0
hEl2gs1dh6t+HQtIcQLTbdrvDl8iPZPUfpOivywrYwuH3KFafMbdSWsRmXmX
Oy8c9YBztuNnpH8tolfculseHzxAPGjAOLKcBqIr0eYo1BP+KQyjkjxoMF3l
amv01RPqWS3i3Wk0WLwdUpMYtAumfO8M3W2mAc+SxYz5tguiquS90TQNHuzb
cOhwKAFBpf66B0ww0PQJqL8YTkCdT1P7elMM9gdcWvY0koDnr9K50zKeDNG6
ZRUjyzPe1Gasw8A4eWp2IYGAYXTHlTTHQHCLTM3IJuBLvPIfWlYY/KzkcrOm
jgCuYKswFzBIZJyRG64n4CI9OvaAPQZdGs3uqxoJkE/tXLnWAYMAfauZ2GYC
hC3q925txyDnx5UOnIey+ydGU+0sDCb2SsQv+whYKdHp/8zHAH7VXbvsNQEJ
V6yfVTlhkOLndhz1E3ClsiIv0hkDo+B7+iWDBBznf7ok3YGBa2KZ79FRAm6f
Hsin78SgIKW/tmiMgO++9q9/LOPhjFV00VcCFrpF6enusj75qTespwhQ8egT
a3li8Kj4Pu3QDAHjsXrV3TLWKZW65c8RIN6WH5CzCwPv6zZXHy8QEEgGfSEI
DG7cPDa1uEgAx+Sa86rdGPwPGOM7RQ==
       "]], 
      LineBox[{{35.51863317699764, 11.596563375425081`}, {
       35.5187070404136, -6.351023695319389}}]},
     Annotation[#, "Charting`Private`Tag$4532#1"]& ]}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{
   "DefaultBoundaryStyle" -> Automatic, 
    "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "DefaultMeshStyle" -> 
    AbsolutePointSize[6], "ScalingFunctions" -> None, 
    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{{0, 100}, {-6.351023695319389, 11.596563375425081`}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{{3.8762203703057337`*^9, 3.876220379980977*^9}, 
   3.8765050131761856`*^9, 3.8765050505154057`*^9, 3.876505166733531*^9, 
   3.876505253203033*^9, 3.8765054682425213`*^9, {3.876505736001601*^9, 
   3.876505760106862*^9}, 3.8765058280963507`*^9, 3.876505866424039*^9, 
   3.8768987397838306`*^9, 3.8769000954185047`*^9, 3.8769010854667377`*^9, 
   3.8769028434196196`*^9, 3.8769121476361966`*^9, {3.8771829651953516`*^9, 
   3.877182971722556*^9}, 3.877183021959716*^9, {3.8782890252122993`*^9, 
   3.878289063704639*^9}, 3.878292725117619*^9, 3.879568890539933*^9, 
   3.8795689360093517`*^9, {3.87956897106921*^9, 3.879568983881357*^9}, {
   3.879569032952899*^9, 3.8795690542179728`*^9}, {3.879571487079197*^9, 
   3.8795714913553343`*^9}, {3.879571532521268*^9, 3.8795715419674997`*^9}, 
   3.880084507024222*^9, 3.8800853736383677`*^9, 3.8800860259425*^9, {
   3.8800860771259003`*^9, 3.88008610210466*^9}, {3.88008616318576*^9, 
   3.880086198003504*^9}, 3.8800862293413563`*^9, {3.880086266837357*^9, 
   3.8800863385534143`*^9}, {3.880086376964778*^9, 3.880086509882822*^9}, {
   3.880087582924012*^9, 3.8800876543745813`*^9}, 3.880094554518929*^9, {
   3.88009550172042*^9, 3.880095549656839*^9}, 3.880450030900723*^9, {
   3.8805353288921747`*^9, 3.880535341553361*^9}, 3.880612842330372*^9, 
   3.880950027979465*^9},
 CellLabel->"Out[14]=",ExpressionUUID->"09944d10-7523-457d-b94c-dbdbf3669e9b"]
}, Open  ]],

Cell[CellGroupData[{

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

Cell[BoxData["49.82150997759732`"], "Output",
 CellChangeTimes->{
  3.8762205374997587`*^9, 3.8764138310515523`*^9, 3.8765050188080196`*^9, 
   3.8765050618099136`*^9, 3.8765058364421587`*^9, 3.8765058704783134`*^9, {
   3.876898758684599*^9, 3.8768987659384336`*^9}, 3.876900107897892*^9, {
   3.876901091009881*^9, 3.8769011091519423`*^9}, {3.8769028519905815`*^9, 
   3.8769028659018292`*^9}, 3.8769121525485616`*^9, 3.876912185418889*^9, {
   3.877182978736745*^9, 3.877183025647654*^9}, 3.878289092820195*^9, 
   3.878292505137838*^9, 3.8795690844906807`*^9, 3.879570642631118*^9, 
   3.879571550250393*^9, 3.880084507942923*^9, 3.88008537507251*^9, 
   3.880085724042804*^9, {3.880085881510993*^9, 3.880085919139512*^9}, 
   3.880086027259701*^9, {3.8800860790573597`*^9, 3.880086104048421*^9}, 
   3.8800861718601227`*^9, {3.880086203223221*^9, 3.88008623091987*^9}, {
   3.8800862683749933`*^9, 3.880086295540452*^9}, 3.880086343066498*^9, {
   3.880086385464319*^9, 3.8800865107303*^9}, {3.8800875839389973`*^9, 
   3.880087671218701*^9}, 3.880094555917272*^9, {3.880095507297576*^9, 
   3.88009555113643*^9}, 3.880450032982963*^9, 3.880535343801962*^9, 
   3.880612844013884*^9, 3.8809500300512943`*^9},
 CellLabel->"Out[15]=",ExpressionUUID->"34307ed9-f9b6-46f0-83b9-36af98fe8641"]
}, Open  ]],

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

Cell[CellGroupData[{

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

Cell[BoxData["1.8472893558629866`"], "Output",
 CellChangeTimes->{
  3.879570442808729*^9, 3.87957064413636*^9, 3.8795715629434137`*^9, 
   3.880084509531536*^9, 3.8800853767467327`*^9, {3.880085901406917*^9, 
   3.880085945613738*^9}, {3.88008606422783*^9, 3.880086105861314*^9}, 
   3.8800861742475863`*^9, 3.880086206862755*^9, {3.8800862701325703`*^9, 
   3.880086297423738*^9}, 3.8800863450137033`*^9, {3.8800863867709208`*^9, 
   3.880086511995571*^9}, {3.8800875856465683`*^9, 3.880087672538249*^9}, 
   3.880094557765129*^9, {3.880095510469705*^9, 3.880095552895302*^9}, 
   3.88045003501324*^9, 3.8805353471630363`*^9, 3.880612846206691*^9, 
   3.880950033016953*^9},
 CellLabel->"Out[16]=",ExpressionUUID->"39e4bab8-8774-4076-8cb8-c3b0ab3ddb14"]
}, Open  ]],

Cell[CellGroupData[{

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

Cell[BoxData["0.7088388229648961`"], "Output",
 CellChangeTimes->{
  3.879569148123713*^9, 3.8795692505602217`*^9, 3.879569281919847*^9, {
   3.879570637910142*^9, 3.879570646893675*^9}, 3.879571580819501*^9, 
   3.880084514207199*^9, {3.880085965594881*^9, 3.880085985349716*^9}, {
   3.880086066066136*^9, 3.880086107886919*^9}, 3.880086176698345*^9, 
   3.8800862093348007`*^9, {3.880086271816217*^9, 3.880086299605596*^9}, 
   3.8800863471623774`*^9, {3.880086388564464*^9, 3.8800865137393713`*^9}, {
   3.880087588111286*^9, 3.88008767445716*^9}, 3.880094560272614*^9, {
   3.8800955121968946`*^9, 3.8800955547288637`*^9}, 3.880450039503292*^9, 
   3.880535349939876*^9, 3.880612851865893*^9, 3.880950035526105*^9},
 CellLabel->"Out[17]=",ExpressionUUID->"80708da8-ef7e-44aa-a0e1-3d0babfbe844"]
}, Open  ]],

Cell[CellGroupData[{

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

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    TagBox[
     {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[CompressedData["
1:eJwVlnk81Nsbx2UpF92yXJJtNNV0IxSlBZ9vVBJtSFIpYy0ibiqSFkWSksYS
IVskhKiopChkSVnKvszMt2xdzIxlcuv3/f1xXuf1/u85n+f9nHM0mT7WrqIi
IiKfqfX/PSmYOXe1S4TJ56Tfj621wkyqmmlW+2nHULWV3dFGuOL3+/IQGi0A
Ld0b9VIIfyRVRGpI0MLx7/TWEiZxFV7G/iemNO7CNBTsN4hB808P1pjGQ3BD
6wh340zM6xxgjGqUIjth+zU5gxKkfkjxHdaoRcKPOpuBP98h/3F4p5dOBzy9
nG0bt7ViSG7oxg6nDhgoXS+t8moF4/QOkxWsDtzMW9ueFd2K+0bSqeyZDrxa
3PVav6cV0bU33BzedULq6qVTjqfacJp9a9z8UDeczPhzB+59QYp0XYjs1j54
nggUK+tuhzs+3R3xY2PzGvPhAqUeiKRNlnVGsJHNKkw3MOxBvLhq14cMNu7o
5Bkl2/WgptZN/WEbGwLPhvQlMT1YYTub5rqRg6elAzajC3vx3YOR2zOHC+kF
5mdL5/ThePT58qYoEouvuDZWDfZBVJDW+/ohifM5cdvqpvqQsL9G5PFbEkpq
TPlKiX58UJU3i+ST+NAdHBiq2Q+trOxqC/tvcLbp1D1o34+RF81NbzW+41R9
+ODqyn54c/9mF+cP4oJKTveWyAH0ObZLvqkdhGiklItrwgBsvl7TaeAMoskq
7j+/rAGsr/sWwF08BJ8tWmv3vRmAWGHmQsWwIZyb1jU34A8gPkjT5MyRYUz9
si7fa8eG1GSTc0jgMJx/sDyMmWyc97kYfitmGNoZJQ5/ebPBZPa2ZNUN41qd
Pxl+lQ3t7UmeXw1HMHL/+u4rT9h4I698d8PCUbzM+ce8WJoDg8ia8q1ao9jO
o006KnLwYO5Zzt5to9i848RaAY2DGzNtuseDRnHWjr2Rs5YDu17W+4TvoxAx
/rK+zpGDoZyF/J8VP5DS7ChanMuBlOYVRu/6MbjX1urnreeixui2bvDWMcxd
arNAh+Ai1D7ZUM16DJJ3riveM+dCNOq5+UHPMaQvfx0JOy6Ev0bcvySNIWGL
pGe6HxfDXfuyP4qO4xb78un2LC5yppkF3gvGIWqWobIynwsPhZPP56uOo6Pu
93X3Yi44ltdrLNeOw6/Yvii7gouusvLv1e7jWJHLv2H2lYv6eMbfFfXjaDum
KGgTI3G92GD1kfZx+NRX141LktjetHnDL+44Ku6Zqc7MJ1E175CF0e9xXLNs
X9SgROLl6ahjz1dPwPjNr3nXV5IIjE7y3W8ygcO+MgaSuiTW5+cETO6YQLnG
sud++iSKuVXhBi4TuK8w8FzGiESe7czDgtgJyARv1be1IuF1cm7R7vQJVB4f
KTHYQ2LlDfmyH48nMD9gwz9zbEk8qFz1Qbt2Au3sI1FHDpJIXsMcyv45gWVp
g8JVx0gc2uUzYS7Jw/JX1tM+XpTXx4OEpAIPmWsclZJ9SMSmxkot0+HBxWHF
43J/EjcX1q1MPcqD/ocZ4cJLJEoN1Up2efNgR7+xZziEBMfRB7PneOCV/DTM
DyWxKU/edn8cDzYeMoWjESTcWlx7xDN5CD5jaHz8Jonon888iop4iMuw2fOR
mqvvFofOz2/kob31whybGBIKvvnzXnTyULimO8E3jgTiRaI9Bnl46fy+8Oxd
qj4y40GlOB+RIpzdRBKJN/On9Hzl+AgN+rtZJIXEiIHFC3UaH6vGVesf3idh
FjL6MWATH1d2GS8qSSfhnQMHhgUfrgQ8lDJJJHy6zWmx4yOduMtgPiDxbprt
fdmFj80aPNuYLBJjGutmdP34CNGf/29BNglV82sh3Rf4kDtcwS+m7gFz7475
EZF8MCV+uaTmkPCL0Y5fn8jHfVqTyelHVP4vg5eQ2XzU6hy4pJ9LopbdlHvn
KR9hRi1a7RTzpeiGm6v4kOGwjN3zSNDW+L/58YmPJzvnPe2m2PJAteW9Xj5S
9QpYRvkkTl9UbrMY5cO6Wbs9hOK0LM+jU0I+vknrhzyhuKHx1VCGpAAhxO3o
OoqnBQv8rRUFuOyuM6+eYroa8/dvugDh86T6iynetaU4PG+1ADcMjWihFAd6
zlU4CAGeJAuqQHFmtH2y5E4BGKcufOyn6mkqzVnx1EGA98GbCC+KZ/tmi5w9
BKhSyFDqpc7DkNxtLHtaABUVK7uNFFvrplaXhwggJdE7e57KI9iOt9frtgDm
zTSZHCqvh+e3dimnCLBL/seVcirPlow4t+pcAY71bPJ4QeUtUj84dqpMgKdf
lZ7dp/qhxdt0bkmNAAv0sv28qX7ZLb4p0dQqQF+C1r2lVD/zPNYs1hoXwJRx
yWRHGomvt65kfP0lgJatZVkp5YP4szadUJlJKGdPP5SlfHGQOGfWz5iExD/C
4cBEEqHaDQ03107C+D8p0xuUb4U2GvZGZpPQ+HRaPpTyUTKt0ivOcRJDFVw9
/Tsk9Gv/mtriNYmWUo9LJOWz45j7pYmASRQWbrO7TPleYiITuzNmEvfX0+L+
CSfB7LR9LVY/CfqqmQqjYBJSif0nfNsnsadCpdv0HIkiB2/VXnISSYfueK05
S0KsIyygbM4U7CPY4w2+JDK+lun7Gk5BbSJzfNKZhFW8+UDPlino3TkpZnGU
8ml/S5SV9RQWCJPTrxwisfXL6CjjxBQOvXVIfk7dD2QrLbsnbQpKrWkj8luo
+Y7Js7MqmIJ2uGlfEUis27dRouzVFMycBJc3bSIR1mLDjPk6Ba9lsd/+WENi
RXOoqtWf03CSe1KQqE7ieNNIVGnANAy+89PieFzIRQWAETaNdReWfm4e5aJs
99wfLNY0NuVeFRV+o97VjxqWJx9P46yuV7B4Fxe5DdYSDO40+sNDb8hVcvHj
w/MA1p4ZrG76FLcokgvfd1eYPgwhBl9v2btekQt1zra5k/pCfHxQs2THAi7q
RP/ICSKEsBN33+QoycUyInI8/IAQnz7yTSKFHLSXxVzMiBCiaIl34NxeDjYX
PEhu/1eIBmVfVddsDmQTazrMnv9E42ZTssGIg8KTMjaLLP7DHFnl5cu82Eg9
5zmkt1+EOPpYo0S6qh+XU2b38kfmEJ10fes2+z6INPcRLwPFiGdxdXFszR4s
P+V21nalBOEwcOSVPasTHEGcnfSrucS0mCDvtks7dnCPn4/2lSTylUv/aDT8
gl2cKC11dSnCTDjSn2pG/ZcqmRdr8qSJJ5b/1GjFNUNHtk3N7Mh84raEqeeF
kU+QgsCpi/cnoRsRDNkLTTAfdZhasm8hEfPJ1bqC0Yivb13mOGjLEmNcp7Nu
iXUouVblSBPIEoG/1DYke9XCTcK/r7FRjjhKHlCIJarh7HeuRCZBnrhDnx1y
c3kHXXetRisnBcIlu1Hl9tNKfJTTlvbf8Bexu0DOwUbqLV7kX07UFlMkLK2y
xMMTK3Buzso99d2KRCD7/qPMw+WQ6grc7l2sRFwOCWL4Nr4A94EmKyVoEcHN
W1uul1QK1gpl8dJ9yoSSocRiX+9niEpwq3aiLSZCJX4yv98qgQPrcKPMt8VE
RPHAld7eJ5jekjmTUq5CNKfUrzNeXQQ21/Ts2zBV4nD59hGnigKkGa9WaXBS
I5amvXjESctH3lV/WfdV6sTTpFztEzq5CNvnbCkmUCfaSvK+rJ56iLVdqQ/D
GjSI28KgVxvfZiH+78S8BksacYadqnWgJxMXKiLjLUtphInog0r7/9KxU2Vf
5oelmsTh5y02djvTsNm7+qBsuCYhczRP36LzPlJUmUvrBZrErLOG1c/WZCyt
zIrddXAJIafTSTKO3EPUqrcrP5cvIS6KsaQLlyXAjSf+qXYlnZAI0ukw+hEH
ZuEELfwWndAI9unRl4yF12Wb2R236cSZUKVYeZFY+NsUf5G5Qyc+npHYOTkV
g2sC/5tRsXQiLsS8uOpbDPI3zPyMS6ITHguvW1ytjsHMm99tDx7Rie5/I+6y
wmIQ9VkmsqqaTlhkF6zn/BGDhPQTHqG1dKJLPOGlpmgM0k81mm2voxNBn/XU
mEIWnipGCesa6cR1A52Tw0MsdDnIezS30okSsQnQ6llgsJXNBth0Yp3jmnuy
t1jQKw5Uz+DSiWsFGieCwljYeLVzxvUbnUhWlZAausCCFSOpYHCIOq9DfHrD
SRb8PGnq4+N0Inq/4oan1iwEGV2aKeLRiRVlU1aGlixcnT/QckpAJ0jTxKUv
zFiIf5weMT1NJ+4KFVTqDFhIuyjuXiakE3IGA4b7V7HwaK+radAsnTiefUSB
XMZC8ZL3aia/6MS7+vDHZ9RZKOctn6EePGLR6cOSMkos/A9M949z
       "]]},
     Annotation[#, "Charting`Private`Tag$6686#1"]& ]}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{
   "DefaultBoundaryStyle" -> Automatic, 
    "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "DefaultMeshStyle" -> 
    AbsolutePointSize[6], "ScalingFunctions" -> None, 
    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{{0, 10}, {0., 0.7088375225314254}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{{3.879570843622718*^9, 3.879570853774201*^9}, 
   3.879571608345409*^9, 3.880084516269706*^9, 3.880086015267064*^9, {
   3.8800860676554127`*^9, 3.8800861091655207`*^9}, 3.880086178702873*^9, 
   3.880086211051835*^9, {3.880086273266088*^9, 3.880086300605966*^9}, 
   3.88008634841398*^9, {3.880086390133037*^9, 3.8800865180731173`*^9}, {
   3.880087591163252*^9, 3.880087675790471*^9}, 3.88009456141636*^9, {
   3.880095513424275*^9, 3.880095556137561*^9}, 3.88045004068221*^9, 
   3.880535351094419*^9, 3.8806128530849047`*^9, 3.880950036573413*^9},
 CellLabel->"Out[18]=",ExpressionUUID->"7811fcb1-4c6e-4b69-84c5-2a4539505a01"]
}, Open  ]],

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

Cell[CellGroupData[{

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

Cell[BoxData["1.990237974515781`"], "Output",
 CellChangeTimes->{{3.880086253659615*^9, 3.880086304157057*^9}, 
   3.880086352535798*^9, {3.8800863938671637`*^9, 3.880086521539238*^9}, {
   3.880087593758895*^9, 3.880087646136713*^9}, 3.8800876779743423`*^9, 
   3.8800945643277493`*^9, {3.8800955160531816`*^9, 3.8800955589398813`*^9}, 
   3.880450043503827*^9},
 CellLabel->"Out[19]=",ExpressionUUID->"347850e7-84fd-46b7-86c9-4d8581f3799d"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  FractionBox["2", "3"], 
  RowBox[{"Sqrt", "[", 
   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"], 
       SuperscriptBox["r", "2"]}], ",", 
      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"], 
       SuperscriptBox["r", "2"]}], ",", 
      RowBox[{"{", 
       RowBox[{"r", ",", "range", ",", "Infinity"}], "}"}]}], "]"}]}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.880095066018161*^9, 3.8800950760469503`*^9}, {
  3.880095479252779*^9, 3.880095482250162*^9}},
 CellLabel->"In[19]:=",ExpressionUUID->"c32958ba-7ddc-4867-b2dd-1fa3dd43d8ea"],

Cell[BoxData["1.5627722778696047`"], "Output",
 CellChangeTimes->{
  3.880095484411269*^9, {3.8800955177242403`*^9, 3.8800955604897947`*^9}, 
   3.880450045304723*^9, 3.880535359224338*^9, 3.880612865473521*^9, 
   3.8809500412555923`*^9},
 CellLabel->"Out[19]=",ExpressionUUID->"0d862d62-21ad-40be-bf86-361b15cd5a3e"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"PsiR", "[", "r_", "]"}], ":=", 
  RowBox[{"Piecewise", "[", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"myNorm", " ", 
        RowBox[{"fIn", "[", 
         RowBox[{
          FractionBox[
           SqrtBox[
            RowBox[{"2", " ", "mass", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "Esep"}], "+", "myU"}], ")"}]}]], "p"], ",", "r", 
          ",", "0"}], "]"}]}], ",", 
       RowBox[{"r", "<", "range"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"myNorm", " ", "myCoeff", " ", 
        RowBox[{"fOut", "[", 
         RowBox[{
          FractionBox[
           SqrtBox[
            RowBox[{"2", " ", "mass", " ", "Esep"}]], "p"], ",", "r", ",", 
          "0"}], "]"}]}], ",", 
       RowBox[{"r", ">", "range"}]}], "}"}]}], "}"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.880450259849861*^9, 3.88045026245728*^9}, {
  3.880450342004263*^9, 3.880450428007826*^9}, {3.8804505083503447`*^9, 
  3.88045054654538*^9}, {3.8804509533799887`*^9, 3.880450954981093*^9}},
 CellLabel->"In[20]:=",ExpressionUUID->"75241ed7-817c-4778-b63c-fc8ddf258a79"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"PsiR", "[", "r", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"r", ",", "0", ",", "10"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.880450434981985*^9, 3.880450437069388*^9}, {
   3.880450467622654*^9, 3.880450475414855*^9}, 3.880450957309256*^9},
 CellLabel->"In[21]:=",ExpressionUUID->"0c53d04a-d321-40e7-bd01-78ab6fb12897"],

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    TagBox[
     {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[CompressedData["
1:eJwVlnk81Nsbx2UpF92yXJJtNNV0IxSlBZ9vVBJtSFIpYy0ibiqSFkWSksYS
IVskhKiopChkSVnKvszMt2xdzIxlcuv3/f1xXuf1/u85n+f9nHM0mT7WrqIi
IiKfqfX/PSmYOXe1S4TJ56Tfj621wkyqmmlW+2nHULWV3dFGuOL3+/IQGi0A
Ld0b9VIIfyRVRGpI0MLx7/TWEiZxFV7G/iemNO7CNBTsN4hB808P1pjGQ3BD
6wh340zM6xxgjGqUIjth+zU5gxKkfkjxHdaoRcKPOpuBP98h/3F4p5dOBzy9
nG0bt7ViSG7oxg6nDhgoXS+t8moF4/QOkxWsDtzMW9ueFd2K+0bSqeyZDrxa
3PVav6cV0bU33BzedULq6qVTjqfacJp9a9z8UDeczPhzB+59QYp0XYjs1j54
nggUK+tuhzs+3R3xY2PzGvPhAqUeiKRNlnVGsJHNKkw3MOxBvLhq14cMNu7o
5Bkl2/WgptZN/WEbGwLPhvQlMT1YYTub5rqRg6elAzajC3vx3YOR2zOHC+kF
5mdL5/ThePT58qYoEouvuDZWDfZBVJDW+/ohifM5cdvqpvqQsL9G5PFbEkpq
TPlKiX58UJU3i+ST+NAdHBiq2Q+trOxqC/tvcLbp1D1o34+RF81NbzW+41R9
+ODqyn54c/9mF+cP4oJKTveWyAH0ObZLvqkdhGiklItrwgBsvl7TaeAMoskq
7j+/rAGsr/sWwF08BJ8tWmv3vRmAWGHmQsWwIZyb1jU34A8gPkjT5MyRYUz9
si7fa8eG1GSTc0jgMJx/sDyMmWyc97kYfitmGNoZJQ5/ebPBZPa2ZNUN41qd
Pxl+lQ3t7UmeXw1HMHL/+u4rT9h4I698d8PCUbzM+ce8WJoDg8ia8q1ao9jO
o006KnLwYO5Zzt5to9i848RaAY2DGzNtuseDRnHWjr2Rs5YDu17W+4TvoxAx
/rK+zpGDoZyF/J8VP5DS7ChanMuBlOYVRu/6MbjX1urnreeixui2bvDWMcxd
arNAh+Ai1D7ZUM16DJJ3riveM+dCNOq5+UHPMaQvfx0JOy6Ev0bcvySNIWGL
pGe6HxfDXfuyP4qO4xb78un2LC5yppkF3gvGIWqWobIynwsPhZPP56uOo6Pu
93X3Yi44ltdrLNeOw6/Yvii7gouusvLv1e7jWJHLv2H2lYv6eMbfFfXjaDum
KGgTI3G92GD1kfZx+NRX141LktjetHnDL+44Ku6Zqc7MJ1E175CF0e9xXLNs
X9SgROLl6ahjz1dPwPjNr3nXV5IIjE7y3W8ygcO+MgaSuiTW5+cETO6YQLnG
sud++iSKuVXhBi4TuK8w8FzGiESe7czDgtgJyARv1be1IuF1cm7R7vQJVB4f
KTHYQ2LlDfmyH48nMD9gwz9zbEk8qFz1Qbt2Au3sI1FHDpJIXsMcyv45gWVp
g8JVx0gc2uUzYS7Jw/JX1tM+XpTXx4OEpAIPmWsclZJ9SMSmxkot0+HBxWHF
43J/EjcX1q1MPcqD/ocZ4cJLJEoN1Up2efNgR7+xZziEBMfRB7PneOCV/DTM
DyWxKU/edn8cDzYeMoWjESTcWlx7xDN5CD5jaHz8Jonon888iop4iMuw2fOR
mqvvFofOz2/kob31whybGBIKvvnzXnTyULimO8E3jgTiRaI9Bnl46fy+8Oxd
qj4y40GlOB+RIpzdRBKJN/On9Hzl+AgN+rtZJIXEiIHFC3UaH6vGVesf3idh
FjL6MWATH1d2GS8qSSfhnQMHhgUfrgQ8lDJJJHy6zWmx4yOduMtgPiDxbprt
fdmFj80aPNuYLBJjGutmdP34CNGf/29BNglV82sh3Rf4kDtcwS+m7gFz7475
EZF8MCV+uaTmkPCL0Y5fn8jHfVqTyelHVP4vg5eQ2XzU6hy4pJ9LopbdlHvn
KR9hRi1a7RTzpeiGm6v4kOGwjN3zSNDW+L/58YmPJzvnPe2m2PJAteW9Xj5S
9QpYRvkkTl9UbrMY5cO6Wbs9hOK0LM+jU0I+vknrhzyhuKHx1VCGpAAhxO3o
OoqnBQv8rRUFuOyuM6+eYroa8/dvugDh86T6iynetaU4PG+1ADcMjWihFAd6
zlU4CAGeJAuqQHFmtH2y5E4BGKcufOyn6mkqzVnx1EGA98GbCC+KZ/tmi5w9
BKhSyFDqpc7DkNxtLHtaABUVK7uNFFvrplaXhwggJdE7e57KI9iOt9frtgDm
zTSZHCqvh+e3dimnCLBL/seVcirPlow4t+pcAY71bPJ4QeUtUj84dqpMgKdf
lZ7dp/qhxdt0bkmNAAv0sv28qX7ZLb4p0dQqQF+C1r2lVD/zPNYs1hoXwJRx
yWRHGomvt65kfP0lgJatZVkp5YP4szadUJlJKGdPP5SlfHGQOGfWz5iExD/C
4cBEEqHaDQ03107C+D8p0xuUb4U2GvZGZpPQ+HRaPpTyUTKt0ivOcRJDFVw9
/Tsk9Gv/mtriNYmWUo9LJOWz45j7pYmASRQWbrO7TPleYiITuzNmEvfX0+L+
CSfB7LR9LVY/CfqqmQqjYBJSif0nfNsnsadCpdv0HIkiB2/VXnISSYfueK05
S0KsIyygbM4U7CPY4w2+JDK+lun7Gk5BbSJzfNKZhFW8+UDPlino3TkpZnGU
8ml/S5SV9RQWCJPTrxwisfXL6CjjxBQOvXVIfk7dD2QrLbsnbQpKrWkj8luo
+Y7Js7MqmIJ2uGlfEUis27dRouzVFMycBJc3bSIR1mLDjPk6Ba9lsd/+WENi
RXOoqtWf03CSe1KQqE7ieNNIVGnANAy+89PieFzIRQWAETaNdReWfm4e5aJs
99wfLNY0NuVeFRV+o97VjxqWJx9P46yuV7B4Fxe5DdYSDO40+sNDb8hVcvHj
w/MA1p4ZrG76FLcokgvfd1eYPgwhBl9v2btekQt1zra5k/pCfHxQs2THAi7q
RP/ICSKEsBN33+QoycUyInI8/IAQnz7yTSKFHLSXxVzMiBCiaIl34NxeDjYX
PEhu/1eIBmVfVddsDmQTazrMnv9E42ZTssGIg8KTMjaLLP7DHFnl5cu82Eg9
5zmkt1+EOPpYo0S6qh+XU2b38kfmEJ10fes2+z6INPcRLwPFiGdxdXFszR4s
P+V21nalBOEwcOSVPasTHEGcnfSrucS0mCDvtks7dnCPn4/2lSTylUv/aDT8
gl2cKC11dSnCTDjSn2pG/ZcqmRdr8qSJJ5b/1GjFNUNHtk3N7Mh84raEqeeF
kU+QgsCpi/cnoRsRDNkLTTAfdZhasm8hEfPJ1bqC0Yivb13mOGjLEmNcp7Nu
iXUouVblSBPIEoG/1DYke9XCTcK/r7FRjjhKHlCIJarh7HeuRCZBnrhDnx1y
c3kHXXetRisnBcIlu1Hl9tNKfJTTlvbf8Bexu0DOwUbqLV7kX07UFlMkLK2y
xMMTK3Buzso99d2KRCD7/qPMw+WQ6grc7l2sRFwOCWL4Nr4A94EmKyVoEcHN
W1uul1QK1gpl8dJ9yoSSocRiX+9niEpwq3aiLSZCJX4yv98qgQPrcKPMt8VE
RPHAld7eJ5jekjmTUq5CNKfUrzNeXQQ21/Ts2zBV4nD59hGnigKkGa9WaXBS
I5amvXjESctH3lV/WfdV6sTTpFztEzq5CNvnbCkmUCfaSvK+rJ56iLVdqQ/D
GjSI28KgVxvfZiH+78S8BksacYadqnWgJxMXKiLjLUtphInog0r7/9KxU2Vf
5oelmsTh5y02djvTsNm7+qBsuCYhczRP36LzPlJUmUvrBZrErLOG1c/WZCyt
zIrddXAJIafTSTKO3EPUqrcrP5cvIS6KsaQLlyXAjSf+qXYlnZAI0ukw+hEH
ZuEELfwWndAI9unRl4yF12Wb2R236cSZUKVYeZFY+NsUf5G5Qyc+npHYOTkV
g2sC/5tRsXQiLsS8uOpbDPI3zPyMS6ITHguvW1ytjsHMm99tDx7Rie5/I+6y
wmIQ9VkmsqqaTlhkF6zn/BGDhPQTHqG1dKJLPOGlpmgM0k81mm2voxNBn/XU
mEIWnipGCesa6cR1A52Tw0MsdDnIezS30okSsQnQ6llgsJXNBth0Yp3jmnuy
t1jQKw5Uz+DSiWsFGieCwljYeLVzxvUbnUhWlZAausCCFSOpYHCIOq9DfHrD
SRb8PGnq4+N0Inq/4oan1iwEGV2aKeLRiRVlU1aGlixcnT/QckpAJ0jTxKUv
zFiIf5weMT1NJ+4KFVTqDFhIuyjuXiakE3IGA4b7V7HwaK+radAsnTiefUSB
XMZC8ZL3aia/6MS7+vDHZ9RZKOctn6EePGLR6cOSMkos/A9M949z
       "]]},
     Annotation[#, "Charting`Private`Tag$12839#1"]& ]}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{
   "DefaultBoundaryStyle" -> Automatic, 
    "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "DefaultMeshStyle" -> 
    AbsolutePointSize[6], "ScalingFunctions" -> None, 
    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{{0, 10}, {0., 0.7088375225314254}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{
  3.880450477448555*^9, {3.8804505437974977`*^9, 3.880450549561881*^9}, 
   3.880450904266193*^9, 3.880450960572151*^9, 3.8805353762901163`*^9, 
   3.880613017272277*^9},
 CellLabel->"Out[21]=",ExpressionUUID->"bfcaf053-8156-415b-a766-14cb460d5213"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"myPsiR", "[", "r_", "]"}], ":=", 
  RowBox[{"Simplify", "[", 
   RowBox[{"PsiR", "[", "r", "]"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.880451016501692*^9, 3.880451021730836*^9}, {
  3.880451090050208*^9, 3.8804511103964*^9}, {3.880451390037879*^9, 
  3.8804513904750967`*^9}, {3.880613601521275*^9, 3.8806136353054867`*^9}},
 CellLabel->"In[21]:=",ExpressionUUID->"726dafeb-13ea-42c8-bb65-5a6d68e2f75b"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"NIntegrate", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"myPsiR", "[", "r", "]"}], 
    RowBox[{"Sin", "[", 
     RowBox[{"10", " ", "r"}], "]"}]}], ",", 
   RowBox[{"{", 
    RowBox[{"r", ",", "0", ",", "Infinity"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.8806136371930943`*^9, 3.880613751137649*^9}, {
  3.880613827324993*^9, 3.8806138290488358`*^9}},
 CellLabel->"In[22]:=",ExpressionUUID->"717762d6-1e69-485a-ab75-de90dc3644c9"],

Cell[BoxData["0.0010508523874437625`"], "Output",
 CellChangeTimes->{3.880613685926548*^9, 3.880950054868312*^9},
 CellLabel->"Out[22]=",ExpressionUUID->"26f52ff1-a852-4dff-9697-265208673679"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"PsiP", "[", "q_", "]"}], ":=", 
  RowBox[{"Simplify", "[", 
   RowBox[{
    SqrtBox[
     FractionBox["2", "\[Pi]"]], 
    RowBox[{"Integrate", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"myPsiR", "[", "r", "]"}], 
       RowBox[{"Sin", "[", 
        RowBox[{"q", " ", "r"}], "]"}]}], ",", 
      RowBox[{"{", 
       RowBox[{"r", ",", "0", ",", "Infinity"}], "}"}]}], "]"}]}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"Print", "[", 
  RowBox[{"PsiP", "[", "q", "]"}], "]"}]}], "Input",
 CellChangeTimes->{{3.8804512707170753`*^9, 3.880451307528953*^9}, {
  3.880451370661892*^9, 3.880451405277438*^9}, {3.8804517177532797`*^9, 
  3.88045175856171*^9}, {3.880451907905079*^9, 3.8804519184783583`*^9}, {
  3.880451961747505*^9, 3.8804519620769672`*^9}, {3.880451997308036*^9, 
  3.880452058978894*^9}, {3.880452125277264*^9, 3.880452125385812*^9}, {
  3.880536060261842*^9, 3.880536067939273*^9}, {3.880536216471127*^9, 
  3.880536222232499*^9}, {3.880613145434863*^9, 3.880613180632429*^9}},
 CellLabel->"In[25]:=",ExpressionUUID->"c1ceda09-1e18-4730-a8a4-798bf3a10de8"],

Cell[BoxData[
 RowBox[{
  SqrtBox[
   FractionBox["2", "\[Pi]"]], " ", 
  RowBox[{"(", 
   TagBox[GridBox[{
      {"\[Piecewise]", GridBox[{
         {
          RowBox[{
           FractionBox[
            RowBox[{"2.901469663670466`*^7", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"4.3395022`*^7", " ", "q", " ", 
                RowBox[{"Cos", "[", 
                 RowBox[{"1.6`", " ", "q"}], "]"}]}], "+", 
               RowBox[{"1.8229557`*^7", " ", 
                RowBox[{"Sin", "[", 
                 RowBox[{"1.6`", " ", "q"}], "]"}]}]}], ")"}]}], 
            RowBox[{"3.32316748416249`*^14", "+", 
             RowBox[{"1.883127934380484`*^15", " ", 
              SuperscriptBox["q", "2"]}]}]], "+", 
           RowBox[{"1.3270979362006529`*^8", " ", 
            RowBox[{"(", 
             RowBox[{
              FractionBox[
               RowBox[{"Sin", "[", 
                RowBox[{"1.6`", " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"-", "1.1933018414019887`"}], "+", "q"}], ")"}]}], 
                "]"}], 
               RowBox[{
                RowBox[{"-", "4.46823272`*^8"}], "+", 
                RowBox[{"3.74442791`*^8", " ", "q"}]}]], "-", 
              FractionBox[
               RowBox[{"1.`", " ", 
                RowBox[{"Sin", "[", 
                 RowBox[{"1.6`", " ", 
                  RowBox[{"(", 
                   RowBox[{
                   "1.1933018414019887`", "\[VeryThinSpace]", "+", "q"}], 
                   ")"}]}], "]"}]}], 
               RowBox[{"4.46823272`*^8", "+", 
                RowBox[{"3.74442791`*^8", " ", "q"}]}]]}], ")"}]}]}], 
          RowBox[{
           RowBox[{"Im", "[", "q", "]"}], "<", "0.42008405941123844`"}]},
         {
          RowBox[{"Integrate", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              TagBox[GridBox[{
                 {"\[Piecewise]", GridBox[{
                    {
                    RowBox[{"0.8458586727012287`", " ", "r", " ", 
                    RowBox[{"SphericalBesselJ", "[", 
                    RowBox[{"0.`", ",", 
                    RowBox[{"1.1933018414019887`", " ", "r"}]}], "]"}]}], 
                    RowBox[{"r", "<", "1.6`"}]},
                    {
                    RowBox[{"1.3094304126855005`", " ", 
                    SuperscriptBox["\[ExponentialE]", 
                    RowBox[{
                    RowBox[{"-", "0.42008405941123844`"}], " ", "r"}]]}], 
                    RowBox[{"r", ">", "1.6`"}]},
                    {"0.`", 
                    TagBox["True",
                    "PiecewiseDefault",
                    AutoDelete->True]}
                   },
                   AllowedDimensions->{2, Automatic},
                   Editable->True,
                   
                   GridBoxAlignment->{
                    "Columns" -> {{Left}}, "Rows" -> {{Baseline}}},
                   
                   GridBoxItemSize->{
                    "Columns" -> {{Automatic}}, "Rows" -> {{1.}}},
                   GridBoxSpacings->{"Columns" -> {
                    Offset[0.27999999999999997`], {
                    Offset[0.84]}, 
                    Offset[0.27999999999999997`]}, "Rows" -> {
                    Offset[0.2], {
                    Offset[0.4]}, 
                    Offset[0.2]}},
                   Selectable->True]}
                },
                
                GridBoxAlignment->{
                 "Columns" -> {{Left}}, "Rows" -> {{Baseline}}},
                
                GridBoxItemSize->{
                 "Columns" -> {{Automatic}}, "Rows" -> {{1.}}},
                GridBoxSpacings->{"Columns" -> {
                    Offset[0.27999999999999997`], {
                    Offset[0.35]}, 
                    Offset[0.27999999999999997`]}, "Rows" -> {
                    Offset[0.2], {
                    Offset[0.4]}, 
                    Offset[0.2]}}],
               "Piecewise",
               DeleteWithContents->True,
               Editable->False,
               SelectWithContents->True,
               Selectable->False,
               StripWrapperBoxes->True], ")"}], " ", 
             RowBox[{"Sin", "[", 
              RowBox[{"q", " ", "r"}], "]"}]}], ",", 
            RowBox[{"{", 
             RowBox[{"r", ",", "0.`", ",", "\[Infinity]"}], "}"}], ",", 
            RowBox[{"Assumptions", "\[Rule]", 
             RowBox[{
              RowBox[{"Im", "[", "q", "]"}], "\[GreaterEqual]", 
              "0.42008405941123844`"}]}]}], "]"}], 
          TagBox["True",
           "PiecewiseDefault",
           AutoDelete->True]}
        },
        AllowedDimensions->{2, Automatic},
        Editable->True,
        GridBoxAlignment->{"Columns" -> {{Left}}, "Rows" -> {{Baseline}}},
        GridBoxItemSize->{"Columns" -> {{Automatic}}, "Rows" -> {{1.}}},
        GridBoxSpacings->{"Columns" -> {
            Offset[0.27999999999999997`], {
             Offset[0.84]}, 
            Offset[0.27999999999999997`]}, "Rows" -> {
            Offset[0.2], {
             Offset[0.4]}, 
            Offset[0.2]}},
        Selectable->True]}
     },
     GridBoxAlignment->{"Columns" -> {{Left}}, "Rows" -> {{Baseline}}},
     GridBoxItemSize->{"Columns" -> {{Automatic}}, "Rows" -> {{1.}}},
     GridBoxSpacings->{"Columns" -> {
         Offset[0.27999999999999997`], {
          Offset[0.35]}, 
         Offset[0.27999999999999997`]}, "Rows" -> {
         Offset[0.2], {
          Offset[0.4]}, 
         Offset[0.2]}}],
    "Piecewise",
    DeleteWithContents->True,
    Editable->False,
    SelectWithContents->True,
    Selectable->False,
    StripWrapperBoxes->True], ")"}]}]], "Print",
 CellChangeTimes->{3.880536153618195*^9, 3.880536286816221*^9, 
  3.8806131692510443`*^9, 3.88061329435061*^9},
 CellLabel->
  "During evaluation of \
In[25]:=",ExpressionUUID->"f909ae20-7144-47db-9cdf-70f17a037f9b"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Integrate", "[", 
  RowBox[{
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{"PsiP", "[", "q", "]"}], ")"}], "2"], ",", 
   RowBox[{"{", 
    RowBox[{"q", ",", "0", ",", "Infinity"}], "}"}], ",", 
   RowBox[{"Assumptions", "\[Rule]", " ", 
    RowBox[{"q", ">", "0"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.880536344175078*^9, 3.8805363711858873`*^9}},
 CellLabel->"In[30]:=",ExpressionUUID->"5b2f680a-fc53-4293-bc32-a5a09105f2c2"],

Cell[BoxData[
 TemplateBox[{
  "Integrate", "idiv", 
   "\"Integral of \\!\\(\\*SuperscriptBox[RowBox[{\\\"(\\\", \
RowBox[{FractionBox[RowBox[{\\\"0.03000292302147212`\\\", \\\" \\\", \
\\\"q\\\"}], RowBox[{\\\"1\\\", \\\"+\\\", RowBox[{\\\"0.0001455306836108173`\
\\\", \\\" \\\", SuperscriptBox[\\\"q\\\", \\\"2\\\"]}]}]], \\\"+\\\", \
FractionBox[RowBox[{SuperscriptBox[\\\"\[ExponentialE]\\\", \
RowBox[{RowBox[{\\\"(\\\", RowBox[{RowBox[{\\\"0.`\\\", \
\\\"\[VeryThinSpace]\\\"}], \\\"-\\\", \
RowBox[{RowBox[{\\\"\[LeftSkeleton]\\\", \\\"21\\\", \
\\\"\[RightSkeleton]\\\"}], \\\" \\\", \\\"\[ImaginaryI]\\\"}]}], \
\\\")\\\"}], \\\" \\\", \\\"q\\\"}]], \\\" \\\", RowBox[{\\\"(\\\", \
RowBox[{RowBox[{\\\"(\\\", RowBox[{RowBox[{\\\"0.`\\\", \
\\\"\[VeryThinSpace]\\\"}], \\\"-\\\", \
RowBox[{RowBox[{\\\"\[LeftSkeleton]\\\", \\\"18\\\", \
\\\"\[RightSkeleton]\\\"}], \\\" \\\", \\\"\[ImaginaryI]\\\"}]}], \
\\\")\\\"}], \\\"+\\\", RowBox[{\\\"\[LeftSkeleton]\\\", \\\"1\\\", \\\"\
\[RightSkeleton]\\\"}], \\\"-\\\", RowBox[{RowBox[{\\\"\[LeftSkeleton]\\\", \
\\\"18\\\", \\\"\[RightSkeleton]\\\"}], \\\" \\\", \\\"q\\\"}]}], \
\\\")\\\"}]}], RowBox[{RowBox[{\\\"-\\\", \\\"55446.43758797824`\\\"}], \\\"+\
\\\", RowBox[{\\\"1.`\\\", \\\" \\\", SuperscriptBox[\\\"q\\\", \
\\\"2\\\"]}]}]], \\\"+\\\", FractionBox[RowBox[{RowBox[{RowBox[{\\\"-\\\", \\\
\"206.16218021559547`\\\"}], \\\" \\\", \\\"q\\\"}], \\\"+\\\", \
RowBox[{\\\"105.27001771946456`\\\", \\\" \\\", \\\"q\\\", \\\" \\\", \
RowBox[{\\\"Cos\\\", \\\"[\\\", RowBox[{\\\"0.008108368342902897`\\\", \\\" \
\\\", \\\"q\\\"}], \\\"]\\\"}]}], \\\"+\\\", RowBox[{\\\"8726.245184279049`\\\
\", \\\" \\\", RowBox[{\\\"Sin\\\", \\\"[\\\", \
RowBox[{\\\"0.008108368342902897`\\\", \\\" \\\", \\\"q\\\"}], \
\\\"]\\\"}]}]}], RowBox[{RowBox[{\\\"6871.403165219998`\\\", \\\"\
\[VeryThinSpace]\\\"}], \\\"+\\\", RowBox[{\\\"1.`\\\", \\\" \\\", \
SuperscriptBox[\\\"q\\\", \\\"2\\\"]}]}]]}], \\\")\\\"}], \\\"2\\\"]\\) does \
not converge on \\!\\(\\*RowBox[{\\\"{\\\", RowBox[{\\\"0\\\", \\\",\\\", \
\\\"\[Infinity]\\\"}], \\\"}\\\"}]\\).\"", 2, 30, 1, 26470880187150754667, 
   "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.880536439018804*^9},
 CellLabel->
  "During evaluation of \
In[30]:=",ExpressionUUID->"aa1683e9-4807-4c48-a219-fedc79848a4f"],

Cell[BoxData[
 RowBox[{"Integrate", "[", 
  RowBox[{
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      FractionBox[
       RowBox[{"0.03000292302147212`", " ", "q"}], 
       RowBox[{"1", "+", 
        RowBox[{"0.0001455306836108173`", " ", 
         SuperscriptBox["q", "2"]}]}]], "+", 
      FractionBox[
       RowBox[{
        SuperscriptBox["\[ExponentialE]", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"0.`", "\[VeryThinSpace]", "-", 
            RowBox[{"0.008108368342902897`", " ", "\[ImaginaryI]"}]}], ")"}], 
          " ", "q"}]], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"0.`", "\[VeryThinSpace]", "-", 
            RowBox[{"4363.122592139528`", " ", "\[ImaginaryI]"}]}], ")"}], 
          "+", 
          RowBox[{
           SuperscriptBox["\[ExponentialE]", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"0.`", "\[VeryThinSpace]", "+", 
               RowBox[{"0.016216736685805793`", " ", "\[ImaginaryI]"}]}], 
              ")"}], " ", "q"}]], " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"0.`", "\[VeryThinSpace]", "+", 
               RowBox[{"4363.122592139528`", " ", "\[ImaginaryI]"}]}], ")"}], 
             "-", 
             RowBox[{"52.63500885973228`", " ", "q"}]}], ")"}]}], "-", 
          RowBox[{"52.63500885973228`", " ", "q"}]}], ")"}]}], 
       RowBox[{
        RowBox[{"-", "55446.43758797824`"}], "+", 
        RowBox[{"1.`", " ", 
         SuperscriptBox["q", "2"]}]}]], "+", 
      FractionBox[
       RowBox[{
        RowBox[{
         RowBox[{"-", "206.16218021559547`"}], " ", "q"}], "+", 
        RowBox[{"105.27001771946456`", " ", "q", " ", 
         RowBox[{"Cos", "[", 
          RowBox[{"0.008108368342902897`", " ", "q"}], "]"}]}], "+", 
        RowBox[{"8726.245184279049`", " ", 
         RowBox[{"Sin", "[", 
          RowBox[{"0.008108368342902897`", " ", "q"}], "]"}]}]}], 
       RowBox[{"6871.403165219998`", "\[VeryThinSpace]", "+", 
        RowBox[{"1.`", " ", 
         SuperscriptBox["q", "2"]}]}]]}], ")"}], "2"], ",", 
   RowBox[{"{", 
    RowBox[{"q", ",", "0", ",", "\[Infinity]"}], "}"}], ",", 
   RowBox[{"Assumptions", "\[Rule]", 
    RowBox[{"q", ">", "0"}]}]}], "]"}]], "Output",
 CellChangeTimes->{3.88053643909369*^9},
 CellLabel->"Out[30]=",ExpressionUUID->"a7707aa6-c791-47d1-aafc-aca2444dec48"]
}, Open  ]]
},
WindowSize->{1389.75, 768.75},
WindowMargins->{{0, Automatic}, {0, Automatic}},
Magnification:>0.9 Inherited,
FrontEndVersion->"13.0 for Linux x86 (64-bit) (December 2, 2021)",
StyleDefinitions->"Default.nb",
ExpressionUUID->"ef468bfc-a077-454d-b697-d1f9ba5b95f7"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 2286, 53, 140, "Input",ExpressionUUID->"ad6df89e-d3d8-490c-ae91-dac6d80d1b5a"],
Cell[2847, 75, 2801, 74, 161, "Input",ExpressionUUID->"2fc4ef7c-d275-441b-9d3d-7d7239f74130"],
Cell[CellGroupData[{
Cell[5673, 153, 2495, 57, 68, "Input",ExpressionUUID->"c2ad0be7-22ee-4300-8ad8-45a6e113c1b0"],
Cell[8171, 212, 9109, 166, 211, "Output",ExpressionUUID->"09944d10-7523-457d-b94c-dbdbf3669e9b"]
}, Open  ]],
Cell[CellGroupData[{
Cell[17317, 383, 2179, 53, 68, "Input",ExpressionUUID->"22f1491b-05f8-4461-9d8a-7023940b0eba"],
Cell[19499, 438, 1292, 18, 30, "Output",ExpressionUUID->"34307ed9-f9b6-46f0-83b9-36af98fe8641"]
}, Open  ]],
Cell[20806, 459, 420, 10, 36, "Input",ExpressionUUID->"0b790137-d050-4e54-820e-ce3bbd34ede5"],
Cell[CellGroupData[{
Cell[21251, 473, 1091, 29, 49, "Input",ExpressionUUID->"444da03a-39cf-4c3d-bb1b-9d322299f630"],
Cell[22345, 504, 755, 11, 30, "Output",ExpressionUUID->"39e4bab8-8774-4076-8cb8-c3b0ab3ddb14"]
}, Open  ]],
Cell[CellGroupData[{
Cell[23137, 520, 1957, 49, 49, "Input",ExpressionUUID->"7835cc6b-0406-456f-a3a3-ab6f89b509c7"],
Cell[25097, 571, 799, 11, 30, "Output",ExpressionUUID->"80708da8-ef7e-44aa-a0e1-3d0babfbe844"]
}, Open  ]],
Cell[CellGroupData[{
Cell[25933, 587, 1343, 35, 49, "Input",ExpressionUUID->"fc251146-842e-4797-bd33-46971b36a427"],
Cell[27279, 624, 6653, 126, 213, "Output",ExpressionUUID->"7811fcb1-4c6e-4b69-84c5-2a4539505a01"]
}, Open  ]],
Cell[33947, 753, 397, 10, 26, "Input",ExpressionUUID->"700162e1-a626-4c8d-a329-bbb787d8421f"],
Cell[CellGroupData[{
Cell[34369, 767, 1112, 33, 49, "Input",ExpressionUUID->"1a4970d6-a8d9-4293-a27c-2d4dbe2240ea"],
Cell[35484, 802, 442, 6, 30, "Output",ExpressionUUID->"347850e7-84fd-46b7-86c9-4d8581f3799d"]
}, Open  ]],
Cell[CellGroupData[{
Cell[35963, 813, 1356, 40, 49, "Input",ExpressionUUID->"c32958ba-7ddc-4867-b2dd-1fa3dd43d8ea"],
Cell[37322, 855, 318, 5, 30, "Output",ExpressionUUID->"0d862d62-21ad-40be-bf86-361b15cd5a3e"]
}, Open  ]],
Cell[37655, 863, 1185, 32, 49, "Input",ExpressionUUID->"75241ed7-817c-4778-b63c-fc8ddf258a79"],
Cell[CellGroupData[{
Cell[38865, 899, 389, 8, 28, "Input",ExpressionUUID->"0c53d04a-d321-40e7-bd01-78ab6fb12897"],
Cell[39257, 909, 6274, 122, 213, "Output",ExpressionUUID->"bfcaf053-8156-415b-a766-14cb460d5213"]
}, Open  ]],
Cell[45546, 1034, 442, 8, 28, "Input",ExpressionUUID->"726dafeb-13ea-42c8-bb65-5a6d68e2f75b"],
Cell[CellGroupData[{
Cell[46013, 1046, 463, 11, 28, "Input",ExpressionUUID->"717762d6-1e69-485a-ab75-de90dc3644c9"],
Cell[46479, 1059, 192, 2, 30, "Output",ExpressionUUID->"26f52ff1-a852-4dff-9697-265208673679"]
}, Open  ]],
Cell[CellGroupData[{
Cell[46708, 1066, 1118, 25, 73, "Input",ExpressionUUID->"c1ceda09-1e18-4730-a8a4-798bf3a10de8"],
Cell[47829, 1093, 5936, 151, 82, "Print",ExpressionUUID->"f909ae20-7144-47db-9cdf-70f17a037f9b"]
}, Open  ]],
Cell[CellGroupData[{
Cell[53802, 1249, 462, 11, 28, "Input",ExpressionUUID->"5b2f680a-fc53-4293-bc32-a5a09105f2c2"],
Cell[54267, 1262, 2335, 39, 42, "Message",ExpressionUUID->"aa1683e9-4807-4c48-a219-fedc79848a4f"],
Cell[56605, 1303, 2423, 62, 79, "Output",ExpressionUUID->"a7707aa6-c791-47d1-aafc-aca2444dec48"]
}, Open  ]]
}
]
*)