resonans1.nb

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

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

(* CreatedBy='Mathematica 13.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     31898,        816]
NotebookOptionsPosition[     28083,        753]
NotebookOutlinePosition[     28510,        770]
CellTagsIndexPosition[     28467,        767]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", "Input", " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"kin", " ", "=", " ", "20."}], ";"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"U0", " ", "=", " ", 
     RowBox[{"-", "40."}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"mass", "=", " ", "939."}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"range", " ", "=", " ", "3"}], ";"}]}]}]], "Input",
 CellChangeTimes->{{3.8698997367340517`*^9, 3.869899843251645*^9}, {
   3.869900162194075*^9, 3.86990022514108*^9}, {3.869900263296324*^9, 
   3.869900278820773*^9}, {3.869900716799706*^9, 3.869900721797647*^9}, 
   3.869902565634445*^9, {3.869902624866561*^9, 3.869902626057645*^9}, {
   3.869904698619644*^9, 3.8699047019298077`*^9}, 3.870142822697193*^9, {
   3.870144764448131*^9, 3.87014476799769*^9}, {3.870144888379427*^9, 
   3.870144970243932*^9}, {3.870145792817855*^9, 
   3.8701457939369917`*^9}},ExpressionUUID->"137265ab-321f-4792-b29f-\
d90a733af3c5"],

Cell[BoxData[
 RowBox[{"Clear", "[", "kIn", "]"}]], "Input",
 CellChangeTimes->{{3.870145066547662*^9, 3.870145071853197*^9}},
 CellLabel->
  "In[137]:=",ExpressionUUID->"02874061-c0e5-4c46-b890-e7a8fd8bd6e5"],

Cell[BoxData[
 RowBox[{
  RowBox[{"kIn", "[", "ee1_", "]"}], " ", ":=", " ", 
  SqrtBox[
   RowBox[{"2", "*", "mass", "*", "ee1"}]]}]], "Input",
 CellChangeTimes->{
  3.870144988515087*^9, {3.870145025281687*^9, 3.870145047838258*^9}},
 CellLabel->
  "In[160]:=",ExpressionUUID->"883ffbce-ec8a-4bea-82fd-15315e8e8450"],

Cell[BoxData[
 RowBox[{"Clear", "[", "kOut", "]"}]], "Input",
 CellChangeTimes->{{3.870145819108061*^9, 3.870145823425865*^9}},
 CellLabel->
  "In[187]:=",ExpressionUUID->"2293bf7c-5f3d-49ed-a114-653f30313e43"],

Cell[BoxData[
 RowBox[{
  RowBox[{"kOut", "[", "ee1_", "]"}], " ", ":=", " ", 
  SqrtBox[
   RowBox[{"2", " ", "mass", " ", 
    RowBox[{"(", 
     RowBox[{"ee1", "-", " ", "U0"}], ")"}]}]]}]], "Input",
 CellChangeTimes->{{3.870145798453801*^9, 3.8701458265877657`*^9}},
 CellLabel->
  "In[188]:=",ExpressionUUID->"11f78f48-8152-4fea-ae39-cefff89cb163"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"momenta", "[", "energy_", "]"}], ":=", " ", 
    FractionBox[
     SqrtBox[
      RowBox[{" ", 
       RowBox[{
        RowBox[{"energy", "*", "energy"}], " ", "+", " ", 
        RowBox[{"2", "*", "mass", "*", "energy"}]}]}]], "197.327"]}], ";"}], 
  "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{{3.8699003273146257`*^9, 3.869900338441976*^9}, {
  3.8699003751305847`*^9, 3.869900433328023*^9}, {3.869900668579546*^9, 
  3.8699007491817617`*^9}, {3.869900797210259*^9, 3.8699008787901163`*^9}},
 CellLabel->
  "In[161]:=",ExpressionUUID->"d12fe0db-0b6c-4c94-a7d8-f0defcf26c82"],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{
  "\:0441\:043e\:0437\:0434\:0430\:044e", " ", 
   "\:0444\:0443\:043d\:043a\:0446\:0438\:0438"}], " ", "*)"}]], "Input",
 CellChangeTimes->{{3.869903554214937*^9, 
  3.869903563251793*^9}},ExpressionUUID->"e3d916de-51d7-4557-9725-\
dda018af67dd"],

Cell[BoxData[
 RowBox[{
  RowBox[{"phiIn", "[", 
   RowBox[{"q_", ",", "r_", ",", "A1_", ",", "angMom_"}], "]"}], " ", ":=", 
  " ", 
  RowBox[{"A1", "*", 
   RowBox[{"SphericalBesselJ", "[", 
    RowBox[{"angMom", ",", 
     RowBox[{"q", "*", "r"}]}], "]"}], " ", "q", " ", "r"}]}]], "Input",
 CellChangeTimes->{{3.869900736428203*^9, 3.869900743569099*^9}, 
   3.8699008753099527`*^9, {3.869902402524983*^9, 3.869902436940699*^9}, {
   3.8699056338975143`*^9, 3.8699056350914783`*^9}},
 CellLabel->
  "In[189]:=",ExpressionUUID->"f52c282a-ced4-44a3-b877-cf9f404c573d"],

Cell[BoxData[
 RowBox[{
  RowBox[{"phiOut", "[", 
   RowBox[{"k_", ",", "r_", ",", "delta_", ",", "angMom_"}], "]"}], ":=", " ", 
  RowBox[{
   RowBox[{"k", " ", "r", " ", 
    RowBox[{"SphericalBesselJ", "[", 
     RowBox[{"angMom", ",", 
      RowBox[{"k", "*", "r"}]}], "]"}], " ", 
    RowBox[{"Cos", "[", "delta", "]"}]}], " ", "+", " ", 
   RowBox[{"k", " ", "r", " ", 
    RowBox[{"SphericalBesselY", "[", 
     RowBox[{"angMom", ",", 
      RowBox[{"k", "*", "r"}]}], "]"}], " ", 
    RowBox[{"Sin", "[", "delta", "]"}]}]}]}]], "Input",
 CellChangeTimes->{{3.8699008921903763`*^9, 3.869900896872324*^9}, {
  3.8699012065049067`*^9, 3.869901327699944*^9}, {3.8699024395507*^9, 
  3.869902440239511*^9}},
 CellLabel->
  "In[190]:=",ExpressionUUID->"73ba3c33-88f7-4598-b488-9b9c4944ce11"],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{
   "\:0438", " ", "\:0438\:0445", " ", 
    "\:043f\:0440\:043e\:0438\:0437\:0432\:043e\:0434\:043d\:044b\:0435"}], 
   " ", "*)"}]}]], "Input",
 CellChangeTimes->{{3.869903571706623*^9, 
  3.869903579390106*^9}},ExpressionUUID->"e74e75f4-3e64-4825-8a63-\
3aae4485a22c"],

Cell[BoxData[
 RowBox[{
  RowBox[{"phiIn1", "[", 
   RowBox[{"q_", ",", "r_", ",", "A1_", ",", "angMom_"}], "]"}], " ", ":=", 
  " ", 
  RowBox[{"Simplify", "[", 
   RowBox[{"D", "[", 
    RowBox[{
     RowBox[{"phiIn", "[", 
      RowBox[{"q", ",", "r", ",", "A1", ",", "angMom"}], "]"}], ",", "r"}], 
    "]"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.869903338277128*^9, 3.869903408044917*^9}, {
  3.8699056459714823`*^9, 3.8699056488705063`*^9}, {3.870148040875915*^9, 
  3.8701480410409946`*^9}},
 CellLabel->
  "In[191]:=",ExpressionUUID->"6f155ea1-2306-4704-a233-560f70b5d462"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Print", "[", 
  RowBox[{"phiIn1", "[", 
   RowBox[{"q", ",", "r", ",", "A1", ",", "angMom"}], "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.869903416610919*^9, 3.869903447278545*^9}, {
  3.870142381422407*^9, 3.870142381722019*^9}},
 CellLabel->
  "In[192]:=",ExpressionUUID->"785e23c6-adf0-4e04-99e7-8de8b70d6714"],

Cell[BoxData[
 RowBox[{
  FractionBox["1", "2"], " ", "A1", " ", "q", " ", 
  RowBox[{"(", 
   RowBox[{
    RowBox[{"q", " ", "r", " ", 
     RowBox[{"SphericalBesselJ", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"-", "1"}], "+", "angMom"}], ",", 
       RowBox[{"q", " ", "r"}]}], "]"}]}], "+", 
    RowBox[{"SphericalBesselJ", "[", 
     RowBox[{"angMom", ",", 
      RowBox[{"q", " ", "r"}]}], "]"}], "-", 
    RowBox[{"q", " ", "r", " ", 
     RowBox[{"SphericalBesselJ", "[", 
      RowBox[{
       RowBox[{"1", "+", "angMom"}], ",", 
       RowBox[{"q", " ", "r"}]}], "]"}]}]}], ")"}]}]], "Print",
 CellChangeTimes->{
  3.869903451021017*^9, 3.869904723106035*^9, 3.8699056602620077`*^9, {
   3.870142358907596*^9, 3.870142383684188*^9}, 3.870142834609494*^9, 
   3.870144787570512*^9, 3.870144856120544*^9, 3.870145267960792*^9, 
   3.87014574889021*^9, 3.870145833367424*^9},
 CellLabel->
  "During evaluation of \
In[192]:=",ExpressionUUID->"a2f6bfed-b9aa-420b-b27b-846a66fae3eb"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"phiOut1", "[", 
   RowBox[{"k_", ",", "r_", ",", "delta_", ",", "angMom_"}], "]"}], " ", ":=",
   " ", 
  RowBox[{"Simplify", "[", 
   RowBox[{"D", "[", 
    RowBox[{
     RowBox[{"phiOut", "[", 
      RowBox[{"k", ",", "r", ",", "delta", ",", "angMom"}], "]"}], ",", "r"}],
     "]"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.869903466278475*^9, 3.869903509164291*^9}},
 CellLabel->
  "In[193]:=",ExpressionUUID->"30b0d700-2070-478c-bc3c-d787b0cb0c36"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Print", "[", 
  RowBox[{"phiOut1", "[", 
   RowBox[{"k", ",", "r", ",", "delta", ",", "angMom"}], "]"}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.869903518460868*^9, 3.8699035336187468`*^9}},
 CellLabel->
  "In[194]:=",ExpressionUUID->"5f2128b6-cab8-40ec-bfa2-267ec247a695"],

Cell[BoxData[
 RowBox[{
  FractionBox["1", "2"], " ", "k", " ", 
  RowBox[{"(", 
   RowBox[{
    RowBox[{"k", " ", "r", " ", 
     RowBox[{"Cos", "[", "delta", "]"}], " ", 
     RowBox[{"SphericalBesselJ", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"-", "1"}], "+", "angMom"}], ",", 
       RowBox[{"k", " ", "r"}]}], "]"}]}], "+", 
    RowBox[{
     RowBox[{"Cos", "[", "delta", "]"}], " ", 
     RowBox[{"SphericalBesselJ", "[", 
      RowBox[{"angMom", ",", 
       RowBox[{"k", " ", "r"}]}], "]"}]}], "-", 
    RowBox[{"k", " ", "r", " ", 
     RowBox[{"Cos", "[", "delta", "]"}], " ", 
     RowBox[{"SphericalBesselJ", "[", 
      RowBox[{
       RowBox[{"1", "+", "angMom"}], ",", 
       RowBox[{"k", " ", "r"}]}], "]"}]}], "+", 
    RowBox[{"k", " ", "r", " ", 
     RowBox[{"Sin", "[", "delta", "]"}], " ", 
     RowBox[{"SphericalBesselY", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"-", "1"}], "+", "angMom"}], ",", 
       RowBox[{"k", " ", "r"}]}], "]"}]}], "+", 
    RowBox[{
     RowBox[{"Sin", "[", "delta", "]"}], " ", 
     RowBox[{"SphericalBesselY", "[", 
      RowBox[{"angMom", ",", 
       RowBox[{"k", " ", "r"}]}], "]"}]}], "-", 
    RowBox[{"k", " ", "r", " ", 
     RowBox[{"Sin", "[", "delta", "]"}], " ", 
     RowBox[{"SphericalBesselY", "[", 
      RowBox[{
       RowBox[{"1", "+", "angMom"}], ",", 
       RowBox[{"k", " ", "r"}]}], "]"}]}]}], ")"}]}]], "Print",
 CellChangeTimes->{3.8699035356894627`*^9, 3.8699047276536207`*^9, 
  3.869905665596455*^9, 3.870142391094543*^9, 3.870142837374173*^9, 
  3.870144790961524*^9, 3.870144864272018*^9, 3.8701452718887568`*^9, 
  3.870145753593706*^9, 3.870145836469759*^9},
 CellLabel->
  "During evaluation of \
In[194]:=",ExpressionUUID->"ec888d6c-c887-43e1-b6e3-22e6c4edb4ec"]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{
  "\:041f\:043e\:0434\:0441\:0442\:0430\:0432\:043b\:044f\:044e", " ", 
   "\:0434\:0430\:043d\:043d\:044b\:0435"}], " ", "*)"}]], "Input",
 CellChangeTimes->{{3.869903587185812*^9, 
  3.869903599028606*^9}},ExpressionUUID->"b7fb5240-5d57-4dc2-a5cf-\
2ef2edd84a37"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"fIn", "[", 
     RowBox[{"A1_", ",", "angMom_"}], "]"}], " ", ":=", " ", 
    RowBox[{"phiIn", "[", 
     RowBox[{
      FractionBox["kIn", "197.327"], ",", "range", ",", "A1", ",", "angMom"}],
      "]"}]}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Clear", "[", "fIn", "]"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"fIn", "[", 
     RowBox[{"ee_", ",", "A1_", ",", "angMom_"}], "]"}], " ", ":=", " ", 
    RowBox[{"phiIn", "[", 
     RowBox[{
      RowBox[{"kIn", "[", "ee", "]"}], ",", "range", ",", "A1", ",", 
      "angMom"}], "]"}]}]}]}]], "Input",
 CellChangeTimes->{{3.869902345120473*^9, 3.869902393018056*^9}, {
  3.8699024490439787`*^9, 3.869902499539846*^9}, {3.8699033471531887`*^9, 
  3.869903353290122*^9}, {3.869903605734902*^9, 3.869903674108346*^9}, {
  3.870142411001401*^9, 3.870142411233304*^9}, {3.8701427696545973`*^9, 
  3.870142769924192*^9}, {3.870144806361532*^9, 3.8701448384469233`*^9}, {
  3.8701451322440577`*^9, 3.8701452339645367`*^9}, {3.870148029511754*^9, 
  3.870148054266265*^9}},
 CellLabel->
  "In[235]:=",ExpressionUUID->"b19389b2-33a5-43a0-ab12-563d4c8b3b62"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Print", "[", 
  RowBox[{"fIn", "[", 
   RowBox[{"ee", ",", "A1", ",", "angMom"}], "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.869904794400037*^9, 3.869904804165381*^9}, {
  3.870142779338613*^9, 3.87014277958138*^9}, {3.870145767563078*^9, 
  3.8701457742342577`*^9}, {3.870148021397421*^9, 3.8701480215655212`*^9}},
 CellLabel->
  "In[237]:=",ExpressionUUID->"e09b75c5-e7ee-41f2-92ef-6d99a9250df2"],

Cell[BoxData[
 RowBox[{"130.0076920801227`", " ", "A1", " ", 
  SqrtBox["ee"], " ", 
  RowBox[{"SphericalBesselJ", "[", 
   RowBox[{"angMom", ",", 
    RowBox[{"130.0076920801227`", " ", 
     SqrtBox["ee"]}]}], "]"}]}]], "Print",
 CellChangeTimes->{
  3.8699048057582283`*^9, 3.869905670676991*^9, 3.870142417093445*^9, 
   3.870142781188665*^9, 3.8701448450752087`*^9, {3.870145147694181*^9, 
   3.870145205575617*^9}, 3.8701452780754757`*^9, {3.870145760135851*^9, 
   3.870145776341847*^9}, 3.870145841956135*^9, 3.870148063131756*^9},
 CellLabel->
  "During evaluation of \
In[237]:=",ExpressionUUID->"c89acb3d-0a4c-49d7-80b0-b609b986aa24"]
}, Open  ]],

Cell[BoxData[
 RowBox[{"Clear", "[", "fOut", "]"}]], "Input",
 CellChangeTimes->{{3.870145889497493*^9, 3.8701458916378307`*^9}},
 CellLabel->
  "In[281]:=",ExpressionUUID->"5b2110d3-922f-468b-93f7-55f20889fabc"],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{"fOut", "[", 
    RowBox[{"delta_", ",", "angMom_"}], "]"}], " ", ":=", " ", 
   RowBox[{"phiOut", "[", 
    RowBox[{
     FractionBox["kOut", "197.327"], ",", "range", ",", "delta", ",", 
     "angMom"}], "]"}]}], "*)"}]], "Input",
 CellChangeTimes->{{3.8701459979441843`*^9, 
  3.870146017022655*^9}},ExpressionUUID->"51bcd62e-b916-47fb-b730-\
ed6a9ea5a103"],

Cell[BoxData[
 RowBox[{
  RowBox[{"fOut", "[", 
   RowBox[{"ee_", ",", "delta_", ",", "angMom_"}], "]"}], " ", ":=", " ", 
  RowBox[{"phiOut", "[", 
   RowBox[{
    RowBox[{"kOut", "[", "ee", "]"}], ",", "range", ",", "delta", ",", 
    "angMom"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.869902502776773*^9, 3.869902522389965*^9}, {
  3.869902654400511*^9, 3.869902659743734*^9}, {3.869902852863531*^9, 
  3.869902859644326*^9}, {3.869903356231689*^9, 3.8699033585596457`*^9}, {
  3.869903634073707*^9, 3.869903650690649*^9}, {3.870145856912568*^9, 
  3.8701458723720713`*^9}, {3.87014920356913*^9, 3.870149206587562*^9}},
 CellLabel->
  "In[282]:=",ExpressionUUID->"a1d36f07-1393-4716-b2a5-7dce8c9b5c36"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Print", "[", 
  RowBox[{"fOut", "[", 
   RowBox[{"ee", ",", "delta", ",", "angMom"}], "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.869904937474552*^9, 3.869904941820552*^9}, {
  3.870145876466538*^9, 3.870145877221031*^9}, {3.870149208544265*^9, 
  3.870149213499585*^9}},
 CellLabel->
  "In[283]:=",ExpressionUUID->"05cf6d8c-1da7-4252-a566-4d8205b8e2d9"],

Cell[BoxData[
 RowBox[{
  RowBox[{"130.0076920801227`", " ", 
   SqrtBox[
    RowBox[{"40.`", "\[VeryThinSpace]", "+", "ee"}]], " ", 
   RowBox[{"Cos", "[", "delta", "]"}], " ", 
   RowBox[{"SphericalBesselJ", "[", 
    RowBox[{"angMom", ",", 
     RowBox[{"130.0076920801227`", " ", 
      SqrtBox[
       RowBox[{"40.`", "\[VeryThinSpace]", "+", "ee"}]]}]}], "]"}]}], "+", 
  RowBox[{"130.0076920801227`", " ", 
   SqrtBox[
    RowBox[{"40.`", "\[VeryThinSpace]", "+", "ee"}]], " ", 
   RowBox[{"Sin", "[", "delta", "]"}], " ", 
   RowBox[{"SphericalBesselY", "[", 
    RowBox[{"angMom", ",", 
     RowBox[{"130.0076920801227`", " ", 
      SqrtBox[
       RowBox[{"40.`", "\[VeryThinSpace]", "+", "ee"}]]}]}], 
    "]"}]}]}]], "Print",
 CellChangeTimes->{3.869904945883834*^9, 3.869905674060755*^9, 
  3.870142433312495*^9, 3.8701427877730618`*^9, 3.8701458954244537`*^9, 
  3.870148230384247*^9, 3.8701492137612953`*^9},
 CellLabel->
  "During evaluation of \
In[283]:=",ExpressionUUID->"83bdcd50-0ad4-45e5-bad5-d8cfb17b57d5"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"fIn1", "[", 
     RowBox[{"A1_", ",", "angMom_"}], "]"}], " ", ":=", 
    RowBox[{"phiIn1", "[", 
     RowBox[{
      FractionBox["kIn", "197.327"], ",", "range", ",", "A1", ",", "angMom"}],
      "]"}]}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Clear", "[", "fIn1", "]"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"fIn1", "[", 
     RowBox[{"A1_", ",", "angMom_", ",", "ee_", ",", "r_"}], "]"}], " ", ":=",
     " ", 
    RowBox[{"Simplify", "[", 
     RowBox[{
      FractionBox["1", "2"], " ", "A1", "  ", 
      RowBox[{"kIn", "[", "ee", "]"}], " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{"kIn", "[", "ee", "]"}], "*", "r", "*", " ", 
         RowBox[{"SphericalBesselJ", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"-", "1"}], "+", "angMom"}], ",", 
           RowBox[{
            RowBox[{"kIn", "[", "ee", "]"}], " ", "r"}]}], "]"}]}], "+", 
        RowBox[{"SphericalBesselJ", "[", 
         RowBox[{"angMom", ",", 
          RowBox[{
           RowBox[{"kIn", "[", "ee", "]"}], " ", "r"}]}], "]"}], "-", 
        RowBox[{
         RowBox[{"kIn", "[", "ee", "]"}], " ", "r", "*", " ", 
         RowBox[{"SphericalBesselJ", "[", 
          RowBox[{
           RowBox[{"1", "+", "angMom"}], ",", 
           RowBox[{
            RowBox[{"kIn", "[", "ee", "]"}], " ", "r"}]}], "]"}]}]}], ")"}]}],
      "]"}]}]}]}]], "Input",
 CellChangeTimes->{{3.869903692263486*^9, 3.8699037607967243`*^9}, {
   3.869903805919382*^9, 3.8699038065397987`*^9}, {3.869904993102446*^9, 
   3.869905015751034*^9}, 3.869905077887982*^9, {3.869905420758458*^9, 
   3.869905459460041*^9}, {3.869905505117428*^9, 3.869905529205653*^9}, 
   3.8699055762480593`*^9, 3.869905727288272*^9, {3.8699057582146273`*^9, 
   3.8699057833774157`*^9}, {3.869905907239922*^9, 3.869905935578825*^9}, {
   3.870142907640519*^9, 3.870142982023192*^9}, {3.87014693443699*^9, 
   3.870146952468712*^9}, {3.870147963882876*^9, 3.870147968392056*^9}, {
   3.870148237519827*^9, 3.87014824701538*^9}, {3.870148362718192*^9, 
   3.870148413818166*^9}, {3.870148513473283*^9, 3.870148525308412*^9}},
 CellLabel->
  "In[272]:=",ExpressionUUID->"3c13eff2-f5aa-4a65-b4fe-877985d5db36"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Print", "[", 
  RowBox[{"fIn1", "[", 
   RowBox[{"A1", ",", "angMom", ",", "ee", ",", "range"}], "]"}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.8699050365884542`*^9, 3.8699050389173403`*^9}, 
   3.8699057283160553`*^9, 3.86990646066066*^9, {3.870142469522188*^9, 
   3.8701424697043743`*^9}, {3.870143006086361*^9, 3.870143019666407*^9}, {
   3.870146975514615*^9, 3.8701469756687326`*^9}, {3.870148250255692*^9, 
   3.870148250944021*^9}, {3.8701484005239563`*^9, 3.870148401023958*^9}, {
   3.870148528349903*^9, 3.8701485291199427`*^9}},
 CellLabel->
  "In[274]:=",ExpressionUUID->"3ce8683e-f553-40fc-b9dc-ca771f703664"],

Cell[BoxData[
 RowBox[{
  RowBox[{"2817.`", " ", "A1", " ", "ee", " ", 
   RowBox[{"SphericalBesselJ", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"-", "1"}], "+", "angMom"}], ",", 
     RowBox[{"130.0076920801227`", " ", 
      SqrtBox["ee"]}]}], "]"}]}], "+", 
  RowBox[{"21.66794868002045`", " ", "A1", " ", 
   SqrtBox["ee"], " ", 
   RowBox[{"SphericalBesselJ", "[", 
    RowBox[{"angMom", ",", 
     RowBox[{"130.0076920801227`", " ", 
      SqrtBox["ee"]}]}], "]"}]}], "-", 
  RowBox[{"2817.`", " ", "A1", " ", "ee", " ", 
   RowBox[{"SphericalBesselJ", "[", 
    RowBox[{
     RowBox[{"1", "+", "angMom"}], ",", 
     RowBox[{"130.0076920801227`", " ", 
      SqrtBox["ee"]}]}], "]"}]}]}]], "Print",
 CellChangeTimes->{
  3.869905960563758*^9, {3.869906456305697*^9, 3.869906462492515*^9}, {
   3.870142445346938*^9, 3.8701424733919773`*^9}, 3.8701427938593597`*^9, 
   3.870142858281001*^9, 3.870143021952096*^9, 3.8701469065602903`*^9, 
   3.8701469813749847`*^9, 3.8701479785357533`*^9, 3.870148257772119*^9, {
   3.870148406395288*^9, 3.870148420319944*^9}, 3.8701485366968107`*^9},
 CellLabel->
  "During evaluation of \
In[274]:=",ExpressionUUID->"5fcd0e1b-e82c-4a52-b9fa-9d34f7748276"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"fOut1", "[", 
     RowBox[{"delta_", ",", "angMom_"}], "]"}], " ", ":=", 
    RowBox[{"phiOut1", "[", 
     RowBox[{
      FractionBox["kOut", "197.327"], ",", "range", ",", "delta", ",", 
      "angMom"}], "]"}]}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Clear", "[", "fOut1", "]"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"fOut1", "[", 
     RowBox[{"delta_", ",", "angMom_", ",", "ee_", ",", "r_"}], "]"}], " ", ":=",
     " ", 
    RowBox[{"Simplify", "[", 
     RowBox[{
      FractionBox["1", "2"], " ", 
      RowBox[{"kOut", "[", "ee", "]"}], " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{"kOut", "[", "ee", "]"}], " ", "r", " ", 
         RowBox[{"Cos", "[", "delta", "]"}], " ", 
         RowBox[{"SphericalBesselJ", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"-", "1"}], "+", "angMom"}], ",", 
           RowBox[{
            RowBox[{"kOut", "[", "ee", "]"}], " ", "r"}]}], "]"}]}], "+", 
        RowBox[{
         RowBox[{"Cos", "[", "delta", "]"}], " ", 
         RowBox[{"SphericalBesselJ", "[", 
          RowBox[{"angMom", ",", 
           RowBox[{
            RowBox[{"kOut", "[", "ee", "]"}], " ", "r"}]}], "]"}]}], "-", 
        RowBox[{
         RowBox[{"kOut", "[", "ee", "]"}], " ", "r", " ", 
         RowBox[{"Cos", "[", "delta", "]"}], " ", 
         RowBox[{"SphericalBesselJ", "[", 
          RowBox[{
           RowBox[{"1", "+", "angMom"}], ",", 
           RowBox[{
            RowBox[{"kOut", "[", "ee", "]"}], " ", "r"}]}], "]"}]}], "+", 
        RowBox[{
         RowBox[{"kOut", "[", "ee", "]"}], " ", "r", " ", 
         RowBox[{"Sin", "[", "delta", "]"}], " ", 
         RowBox[{"SphericalBesselY", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"-", "1"}], "+", "angMom"}], ",", 
           RowBox[{
            RowBox[{"kOut", "[", "ee", "]"}], " ", "r"}]}], "]"}]}], "+", 
        RowBox[{
         RowBox[{"Sin", "[", "delta", "]"}], " ", 
         RowBox[{"SphericalBesselY", "[", 
          RowBox[{"angMom", ",", 
           RowBox[{
            RowBox[{"kOut", "[", "ee", "]"}], " ", "r"}]}], "]"}]}], "-", 
        RowBox[{
         RowBox[{"kOut", "[", "ee", "]"}], " ", "r", " ", 
         RowBox[{"Sin", "[", "delta", "]"}], " ", 
         RowBox[{"SphericalBesselY", "[", 
          RowBox[{
           RowBox[{"1", "+", "angMom"}], ",", 
           RowBox[{
            RowBox[{"kOut", "[", "ee", "]"}], " ", "r"}]}], "]"}]}]}], 
       ")"}]}], "]"}]}]}]}]], "Input",
 CellChangeTimes->{{3.869905973732863*^9, 3.869906046015753*^9}, {
  3.8701431002750063`*^9, 3.870143137123361*^9}, {3.870147001538087*^9, 
  3.870147030213578*^9}, {3.8701479542811413`*^9, 3.870147959423932*^9}, {
  3.870148268064863*^9, 3.870148281708652*^9}},
 CellLabel->
  "In[277]:=",ExpressionUUID->"4e60cda4-270f-41c1-a4f4-78d03b351537"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Print", "[", 
  RowBox[{"fOut1", "[", 
   RowBox[{"delta", ",", "angMom", ",", "ee", ",", "range"}], "]"}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.869906066697554*^9, 3.8699060728682833`*^9}, {
  3.870143144626992*^9, 3.870143158426641*^9}, {3.8701470403211727`*^9, 
  3.870147040458294*^9}, {3.8701482909067717`*^9, 3.8701482947231627`*^9}, {
  3.87014845357169*^9, 3.870148454689481*^9}},
 CellLabel->
  "In[279]:=",ExpressionUUID->"107cc588-8aa8-49dd-a831-be09f57f8eda"],

Cell[BoxData[
 RowBox[{"21.66794868002045`", " ", 
  SqrtBox[
   RowBox[{"40.`", "\[VeryThinSpace]", "+", "ee"}]], " ", 
  RowBox[{"(", 
   RowBox[{
    RowBox[{"130.0076920801227`", " ", 
     SqrtBox[
      RowBox[{"40.`", "\[VeryThinSpace]", "+", "ee"}]], " ", 
     RowBox[{"Cos", "[", "delta", "]"}], " ", 
     RowBox[{"SphericalBesselJ", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"-", "1"}], "+", "angMom"}], ",", 
       RowBox[{"130.0076920801227`", " ", 
        SqrtBox[
         RowBox[{"40.`", "\[VeryThinSpace]", "+", "ee"}]]}]}], "]"}]}], "+", 
    RowBox[{
     RowBox[{"Cos", "[", "delta", "]"}], " ", 
     RowBox[{"SphericalBesselJ", "[", 
      RowBox[{"angMom", ",", 
       RowBox[{"130.0076920801227`", " ", 
        SqrtBox[
         RowBox[{"40.`", "\[VeryThinSpace]", "+", "ee"}]]}]}], "]"}]}], "-", 
    RowBox[{"130.0076920801227`", " ", 
     SqrtBox[
      RowBox[{"40.`", "\[VeryThinSpace]", "+", "ee"}]], " ", 
     RowBox[{"Cos", "[", "delta", "]"}], " ", 
     RowBox[{"SphericalBesselJ", "[", 
      RowBox[{
       RowBox[{"1", "+", "angMom"}], ",", 
       RowBox[{"130.0076920801227`", " ", 
        SqrtBox[
         RowBox[{"40.`", "\[VeryThinSpace]", "+", "ee"}]]}]}], "]"}]}], "+", 
    RowBox[{"130.0076920801227`", " ", 
     SqrtBox[
      RowBox[{"40.`", "\[VeryThinSpace]", "+", "ee"}]], " ", 
     RowBox[{"Sin", "[", "delta", "]"}], " ", 
     RowBox[{"SphericalBesselY", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"-", "1"}], "+", "angMom"}], ",", 
       RowBox[{"130.0076920801227`", " ", 
        SqrtBox[
         RowBox[{"40.`", "\[VeryThinSpace]", "+", "ee"}]]}]}], "]"}]}], "+", 
    RowBox[{
     RowBox[{"Sin", "[", "delta", "]"}], " ", 
     RowBox[{"SphericalBesselY", "[", 
      RowBox[{"angMom", ",", 
       RowBox[{"130.0076920801227`", " ", 
        SqrtBox[
         RowBox[{"40.`", "\[VeryThinSpace]", "+", "ee"}]]}]}], "]"}]}], "-", 
    RowBox[{"130.0076920801227`", " ", 
     SqrtBox[
      RowBox[{"40.`", "\[VeryThinSpace]", "+", "ee"}]], " ", 
     RowBox[{"Sin", "[", "delta", "]"}], " ", 
     RowBox[{"SphericalBesselY", "[", 
      RowBox[{
       RowBox[{"1", "+", "angMom"}], ",", 
       RowBox[{"130.0076920801227`", " ", 
        SqrtBox[
         RowBox[{"40.`", "\[VeryThinSpace]", "+", "ee"}]]}]}], "]"}]}]}], 
   ")"}]}]], "Print",
 CellChangeTimes->{3.869906075725178*^9, 3.870142456725*^9, 
  3.8701431647662573`*^9, 3.870147042050164*^9, 3.870147987933353*^9, 
  3.870148297603366*^9, 3.870148466801008*^9, 3.870148543467256*^9},
 CellLabel->
  "During evaluation of \
In[279]:=",ExpressionUUID->"683f7d8f-5bdf-445b-b5fa-cda0d3db7032"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"NSolve", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{"fIn", "[", 
       RowBox[{"ee", ",", "aa", ",", "1"}], "]"}], "==", 
      RowBox[{"fOut", "[", 
       RowBox[{"ee", ",", "dd", ",", "1"}], "]"}]}], ",", 
     RowBox[{
      RowBox[{"fIn1", "[", 
       RowBox[{"aa", ",", "1", ",", "ee", ",", "range"}], "]"}], "==", 
      RowBox[{"fOut1", "[", 
       RowBox[{"dd", ",", "1", ",", "ee", ",", "range"}], "]"}]}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{"aa", ",", "dd"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.8699038602251577`*^9, 3.869904101667548*^9}, {
   3.8699041723855553`*^9, 3.869904210530024*^9}, 3.869904252311246*^9, {
   3.869904396536236*^9, 3.869904399027741*^9}, {3.8701431958772173`*^9, 
   3.870143225476993*^9}, {3.870143652479889*^9, 3.870143678067491*^9}, {
   3.870143757110928*^9, 3.870143786639409*^9}, {3.870143823852461*^9, 
   3.8701438346853113`*^9}, {3.8701444368885393`*^9, 3.870144480820072*^9}, {
   3.8701471863516808`*^9, 3.870147281566517*^9}, {3.870148494940176*^9, 
   3.8701485080902767`*^9}, {3.870149182748822*^9, 3.8701491897996817`*^9}, {
   3.870149242903923*^9, 3.870149243039599*^9}, {3.870149473732283*^9, 
   3.8701494836969976`*^9}, {3.870149514821795*^9, 
   3.8701495238721933`*^9}},ExpressionUUID->"2259a25f-7509-401a-afac-\
84d646e2336f"],

Cell[BoxData["$Aborted"], "Output",
 CellChangeTimes->{
  3.870143787488616*^9, {3.870143828888331*^9, 3.870143835761471*^9}, {
   3.870144453275545*^9, 3.870144481957519*^9}, 3.870147408223731*^9, 
   3.870147928822309*^9, 3.8701480182199383`*^9, 3.870149179379559*^9, 
   3.870149245263503*^9, 3.870149492543207*^9},
 CellLabel->
  "Out[285]=",ExpressionUUID->"7679836d-076b-40df-b7cb-b2861ef0c59e"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"phase", "[", "ee_", "]"}], ":=", " "}]], "Input",
 CellChangeTimes->{{3.870148760144835*^9, 
  3.870148768918832*^9}},ExpressionUUID->"ea8a29ab-c052-471d-ae15-\
ee682cfe1269"],

Cell[BoxData[
 RowBox[{
  RowBox[{"sMatrix", "[", "ee_", "]"}], " ", ":=", " ", 
  RowBox[{"Exp", "[", 
   RowBox[{"2", "*", "I", "*", 
    RowBox[{"phase", "[", "ee", "]"}]}], "]"}], "  "}]], "Input",
 CellChangeTimes->{{3.853486741956423*^9, 3.853486748574054*^9}, 
   3.8534869884309907`*^9, {3.8701487152125597`*^9, 
   3.870148750228079*^9}},ExpressionUUID->"720974a2-687b-47f2-a4c9-\
1c81aa8e5e65"]
},
WindowSize->{1389.75, 768.75},
WindowMargins->{{0, Automatic}, {0, Automatic}},
Magnification:>1. Inherited,
FrontEndVersion->"13.0 for Linux x86 (64-bit) (December 2, 2021)",
StyleDefinitions->"Default.nb",
ExpressionUUID->"938b4164-1517-40b9-833c-732c6639372d"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 1047, 23, 113, "Input",ExpressionUUID->"137265ab-321f-4792-b29f-d90a733af3c5"],
Cell[1608, 45, 209, 4, 29, "Input",ExpressionUUID->"02874061-c0e5-4c46-b890-e7a8fd8bd6e5"],
Cell[1820, 51, 318, 8, 30, "Input",ExpressionUUID->"883ffbce-ec8a-4bea-82fd-15315e8e8450"],
Cell[2141, 61, 210, 4, 29, "Input",ExpressionUUID->"2293bf7c-5f3d-49ed-a114-653f30313e43"],
Cell[2354, 67, 353, 9, 31, "Input",ExpressionUUID->"11f78f48-8152-4fea-ae39-cefff89cb163"],
Cell[2710, 78, 646, 16, 74, "Input",ExpressionUUID->"d12fe0db-0b6c-4c94-a7d8-f0defcf26c82"],
Cell[3359, 96, 293, 7, 29, "Input",ExpressionUUID->"e3d916de-51d7-4557-9725-dda018af67dd"],
Cell[3655, 105, 570, 13, 29, "Input",ExpressionUUID->"f52c282a-ced4-44a3-b877-cf9f404c573d"],
Cell[4228, 120, 793, 19, 29, "Input",ExpressionUUID->"73ba3c33-88f7-4598-b488-9b9c4944ce11"],
Cell[5024, 141, 351, 9, 51, "Input",ExpressionUUID->"e74e75f4-3e64-4825-8a63-3aae4485a22c"],
Cell[5378, 152, 585, 15, 29, "Input",ExpressionUUID->"6f155ea1-2306-4704-a233-560f70b5d462"],
Cell[CellGroupData[{
Cell[5988, 171, 339, 7, 29, "Input",ExpressionUUID->"785e23c6-adf0-4e04-99e7-8de8b70d6714"],
Cell[6330, 180, 999, 26, 38, "Print",ExpressionUUID->"a2f6bfed-b9aa-420b-b27b-846a66fae3eb"]
}, Open  ]],
Cell[7344, 209, 489, 13, 29, "Input",ExpressionUUID->"30b0d700-2070-478c-bc3c-d787b0cb0c36"],
Cell[CellGroupData[{
Cell[7858, 226, 299, 7, 29, "Input",ExpressionUUID->"5f2128b6-cab8-40ec-bfa2-267ec247a695"],
Cell[8160, 235, 1776, 47, 59, "Print",ExpressionUUID->"ec888d6c-c887-43e1-b6e3-22e6c4edb4ec"]
}, Open  ]],
Cell[9951, 285, 311, 7, 29, "Input",ExpressionUUID->"b7fb5240-5d57-4dc2-a5cf-2ef2edd84a37"],
Cell[10265, 294, 1190, 27, 77, "Input",ExpressionUUID->"b19389b2-33a5-43a0-ab12-563d4c8b3b62"],
Cell[CellGroupData[{
Cell[11480, 325, 425, 8, 29, "Input",ExpressionUUID->"e09b75c5-e7ee-41f2-92ef-6d99a9250df2"],
Cell[11908, 335, 645, 14, 25, "Print",ExpressionUUID->"c89acb3d-0a4c-49d7-80b0-b609b986aa24"]
}, Open  ]],
Cell[12568, 352, 212, 4, 29, "Input",ExpressionUUID->"5b2110d3-922f-468b-93f7-55f20889fabc"],
Cell[12783, 358, 412, 11, 35, "Input",ExpressionUUID->"51bcd62e-b916-47fb-b730-ed6a9ea5a103"],
Cell[13198, 371, 707, 14, 29, "Input",ExpressionUUID->"a1d36f07-1393-4716-b2a5-7dce8c9b5c36"],
Cell[CellGroupData[{
Cell[13930, 389, 380, 8, 29, "Input",ExpressionUUID->"05cf6d8c-1da7-4252-a566-4d8205b8e2d9"],
Cell[14313, 399, 1030, 26, 25, "Print",ExpressionUUID->"83bdcd50-0ad4-45e5-bad5-d8cfb17b57d5"]
}, Open  ]],
Cell[15358, 428, 2282, 53, 111, "Input",ExpressionUUID->"3c13eff2-f5aa-4a65-b4fe-877985d5db36"],
Cell[CellGroupData[{
Cell[17665, 485, 649, 12, 29, "Input",ExpressionUUID->"3ce8683e-f553-40fc-b9dc-ca771f703664"],
Cell[18317, 499, 1203, 29, 25, "Print",ExpressionUUID->"5fcd0e1b-e82c-4a52-b9fa-9d34f7748276"]
}, Open  ]],
Cell[19535, 531, 2924, 74, 187, "Input",ExpressionUUID->"4e60cda4-270f-41c1-a4f4-78d03b351537"],
Cell[CellGroupData[{
Cell[22484, 609, 500, 10, 29, "Input",ExpressionUUID->"107cc588-8aa8-49dd-a831-be09f57f8eda"],
Cell[22987, 621, 2648, 68, 67, "Print",ExpressionUUID->"683f7d8f-5bdf-445b-b5fa-cda0d3db7032"]
}, Open  ]],
Cell[CellGroupData[{
Cell[25672, 694, 1371, 29, 29, "Input",ExpressionUUID->"2259a25f-7509-401a-afac-84d646e2336f"],
Cell[27046, 725, 401, 7, 33, "Output",ExpressionUUID->"7679836d-076b-40df-b7cb-b2861ef0c59e"]
}, Open  ]],
Cell[27462, 735, 210, 5, 29, "Input",ExpressionUUID->"ea8a29ab-c052-471d-ae15-ee682cfe1269"],
Cell[27675, 742, 404, 9, 29, "Input",ExpressionUUID->"720974a2-687b-47f2-a4c9-1c81aa8e5e65"]
}
]
*)