(* 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[ 34923, 927] NotebookOptionsPosition[ 30269, 851] NotebookOutlinePosition[ 30667, 867] CellTagsIndexPosition[ 30624, 864] 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}}, 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[{"aa_", ",", "A1_", ",", "angMom_"}], "]"}], " ", ":=", " ", RowBox[{"phiIn", "[", RowBox[{ RowBox[{"kIn", "[", "aa", "]"}], ",", "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}}, CellLabel-> "In[197]:=",ExpressionUUID->"b19389b2-33a5-43a0-ab12-563d4c8b3b62"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Print", "[", RowBox[{"fIn", "[", RowBox[{"aa", ",", "A1", ",", "angMom"}], "]"}], "]"}]], "Input", CellChangeTimes->{{3.869904794400037*^9, 3.869904804165381*^9}, { 3.870142779338613*^9, 3.87014277958138*^9}, {3.870145767563078*^9, 3.8701457742342577`*^9}}, CellLabel-> "In[199]:=",ExpressionUUID->"e09b75c5-e7ee-41f2-92ef-6d99a9250df2"], Cell[BoxData[ RowBox[{"130.0076920801227`", " ", "A1", " ", SqrtBox["aa"], " ", RowBox[{"SphericalBesselJ", "[", RowBox[{"angMom", ",", RowBox[{"130.0076920801227`", " ", SqrtBox["aa"]}]}], "]"}]}]], "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}, CellLabel-> "During evaluation of \ In[199]:=",ExpressionUUID->"fd042a95-017b-41dc-970a-caa921531275"] }, Open ]], Cell[BoxData[ RowBox[{"Clear", "[", "fOut", "]"}]], "Input", CellChangeTimes->{{3.870145889497493*^9, 3.8701458916378307`*^9}}, CellLabel-> "In[201]:=",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[{"aa_", ",", "delta_", ",", "angMom_"}], "]"}], " ", ":=", " ", RowBox[{"phiOut", "[", RowBox[{ RowBox[{"kOut", "[", "aa", "]"}], ",", "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}}, CellLabel-> "In[202]:=",ExpressionUUID->"a1d36f07-1393-4716-b2a5-7dce8c9b5c36"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Print", "[", RowBox[{"fOut", "[", RowBox[{"aa", ",", "delta", ",", "angMom"}], "]"}], "]"}]], "Input", CellChangeTimes->{{3.869904937474552*^9, 3.869904941820552*^9}, { 3.870145876466538*^9, 3.870145877221031*^9}}, CellLabel-> "In[203]:=",ExpressionUUID->"05cf6d8c-1da7-4252-a566-4d8205b8e2d9"], Cell[BoxData[ RowBox[{ RowBox[{"130.0076920801227`", " ", SqrtBox[ RowBox[{"40.`", "\[VeryThinSpace]", "+", "aa"}]], " ", RowBox[{"Cos", "[", "delta", "]"}], " ", RowBox[{"SphericalBesselJ", "[", RowBox[{"angMom", ",", RowBox[{"130.0076920801227`", " ", SqrtBox[ RowBox[{"40.`", "\[VeryThinSpace]", "+", "aa"}]]}]}], "]"}]}], "+", RowBox[{"130.0076920801227`", " ", SqrtBox[ RowBox[{"40.`", "\[VeryThinSpace]", "+", "aa"}]], " ", RowBox[{"Sin", "[", "delta", "]"}], " ", RowBox[{"SphericalBesselY", "[", RowBox[{"angMom", ",", RowBox[{"130.0076920801227`", " ", SqrtBox[ RowBox[{"40.`", "\[VeryThinSpace]", "+", "aa"}]]}]}], "]"}]}]}]], "Print", CellChangeTimes->{3.869904945883834*^9, 3.869905674060755*^9, 3.870142433312495*^9, 3.8701427877730618`*^9, 3.8701458954244537`*^9}, CellLabel-> "During evaluation of \ In[203]:=",ExpressionUUID->"3e65c34a-a20d-49b9-a06a-b22c57e55b93"] }, 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_", ",", "q_", ",", "r_"}], "]"}], " ", ":=", " ", 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"}]}], "]"}]}]}], ")"}]}]}]}]}]], "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}}, CellLabel-> "In[206]:=",ExpressionUUID->"3c13eff2-f5aa-4a65-b4fe-877985d5db36"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Print", "[", RowBox[{"fIn1", "[", RowBox[{"A1", ",", "angMom", ",", FractionBox["kIn", "197.327"], ",", "3"}], "]"}], "]"}]], "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}}, CellLabel->"In[61]:=",ExpressionUUID->"3ce8683e-f553-40fc-b9dc-ca771f703664"], Cell[BoxData[ RowBox[{"0.49107325587986`", " ", "A1", " ", RowBox[{"(", RowBox[{ RowBox[{"2.94643953527916`", " ", RowBox[{"SphericalBesselJ", "[", RowBox[{ RowBox[{ RowBox[{"-", "1"}], "+", "angMom"}], ",", "2.94643953527916`"}], "]"}]}], "+", RowBox[{"SphericalBesselJ", "[", RowBox[{"angMom", ",", "2.94643953527916`"}], "]"}], "-", RowBox[{"2.94643953527916`", " ", RowBox[{"SphericalBesselJ", "[", RowBox[{ RowBox[{"1", "+", "angMom"}], ",", "2.94643953527916`"}], "]"}]}]}], ")"}]}]], "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}, CellLabel-> "During evaluation of \ In[61]:=",ExpressionUUID->"71adf521-737e-4433-9e68-6d2a46404b8b"] }, 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_", ",", "k_", ",", "r_"}], "]"}], " ", ":=", " ", 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"}]}], "]"}]}]}], ")"}]}]}]}]}]], "Input", CellChangeTimes->{{3.869905973732863*^9, 3.869906046015753*^9}, { 3.8701431002750063`*^9, 3.870143137123361*^9}}, CellLabel->"In[62]:=",ExpressionUUID->"4e60cda4-270f-41c1-a4f4-78d03b351537"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Print", "[", RowBox[{"fOut1", "[", RowBox[{"delta", ",", "angMom", ",", FractionBox["kOut", "197.327"], ",", "3"}], "]"}], "]"}]], "Input", CellChangeTimes->{{3.869906066697554*^9, 3.8699060728682833`*^9}, { 3.870143144626992*^9, 3.870143158426641*^9}}, CellLabel->"In[64]:=",ExpressionUUID->"107cc588-8aa8-49dd-a831-be09f57f8eda"], Cell[BoxData[ RowBox[{"0.8505638294221894`", " ", RowBox[{"(", RowBox[{ RowBox[{"5.103382976533137`", " ", RowBox[{"Cos", "[", "delta", "]"}], " ", RowBox[{"SphericalBesselJ", "[", RowBox[{ RowBox[{ RowBox[{"-", "1"}], "+", "angMom"}], ",", "5.103382976533137`"}], "]"}]}], "+", RowBox[{ RowBox[{"Cos", "[", "delta", "]"}], " ", RowBox[{"SphericalBesselJ", "[", RowBox[{"angMom", ",", "5.103382976533137`"}], "]"}]}], "-", RowBox[{"5.103382976533137`", " ", RowBox[{"Cos", "[", "delta", "]"}], " ", RowBox[{"SphericalBesselJ", "[", RowBox[{ RowBox[{"1", "+", "angMom"}], ",", "5.103382976533137`"}], "]"}]}], "+", RowBox[{"5.103382976533137`", " ", RowBox[{"Sin", "[", "delta", "]"}], " ", RowBox[{"SphericalBesselY", "[", RowBox[{ RowBox[{ RowBox[{"-", "1"}], "+", "angMom"}], ",", "5.103382976533137`"}], "]"}]}], "+", RowBox[{ RowBox[{"Sin", "[", "delta", "]"}], " ", RowBox[{"SphericalBesselY", "[", RowBox[{"angMom", ",", "5.103382976533137`"}], "]"}]}], "-", RowBox[{"5.103382976533137`", " ", RowBox[{"Sin", "[", "delta", "]"}], " ", RowBox[{"SphericalBesselY", "[", RowBox[{ RowBox[{"1", "+", "angMom"}], ",", "5.103382976533137`"}], "]"}]}]}], ")"}]}]], "Print", CellChangeTimes->{3.869906075725178*^9, 3.870142456725*^9, 3.8701431647662573`*^9}, CellLabel-> "During evaluation of \ In[64]:=",ExpressionUUID->"0a0cafd0-cdc3-47ee-924c-3fa7aeb43452"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"NSolve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"fIn", "[", RowBox[{"aa", ",", "1"}], "]"}], "==", RowBox[{"fOut", "[", RowBox[{"dd", ",", "1"}], "]"}]}], ",", RowBox[{ RowBox[{"fIn1", "[", RowBox[{"aa", ",", "1", ",", FractionBox["kIn", "197.327"], ",", "range"}], "]"}], "==", RowBox[{"fOut1", "[", RowBox[{"dd", ",", "1", ",", FractionBox["kOut", "197.327"], ",", "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}}, CellLabel->"In[79]:=",ExpressionUUID->"2259a25f-7509-401a-afac-84d646e2336f"], Cell[BoxData[ TemplateBox[{ "NSolve", "ifun", "\"Inverse functions are being used by \ \\!\\(\\*RowBox[{\\\"NSolve\\\"}]\\), so some solutions may not be found; use \ Reduce for complete solution information.\"", 2, 79, 15, 26402770535264617284, "Local"}, "MessageTemplate"]], "Message", "MSG", CellChangeTimes->{{3.870143784640518*^9, 3.870143787476328*^9}, { 3.870143828863248*^9, 3.870143835748345*^9}, {3.870144453264967*^9, 3.870144481932513*^9}}, CellLabel-> "During evaluation of \ In[79]:=",ExpressionUUID->"9244a4dd-143f-45ab-a98d-28422173269e"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"aa", "\[Rule]", RowBox[{"-", "0.9699407033214296`"}]}], ",", RowBox[{"dd", "\[Rule]", RowBox[{"-", "1.0709952516350376`"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"aa", "\[Rule]", "0.9699407033214296`"}], ",", RowBox[{"dd", "\[Rule]", "2.0705974019547555`"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{ 3.870143787488616*^9, {3.870143828888331*^9, 3.870143835761471*^9}, { 3.870144453275545*^9, 3.870144481957519*^9}}, CellLabel->"Out[79]=",ExpressionUUID->"63a47552-f5f5-426b-be4c-082d226082c7"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"NSolve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"fIn", "[", RowBox[{"aa", ",", "2"}], "]"}], "==", RowBox[{"fOut", "[", RowBox[{"dd", ",", "2"}], "]"}]}], ",", RowBox[{ RowBox[{"fIn1", "[", RowBox[{"aa", ",", "2", ",", FractionBox["kIn", "197.327"], ",", "range"}], "]"}], "==", RowBox[{"fOut1", "[", RowBox[{"dd", ",", "2", ",", FractionBox["kOut", "197.327"], ",", "range"}], "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"aa", ",", "dd"}], "}"}]}], "]"}]], "Input", CellChangeTimes->{{3.870144470121334*^9, 3.870144492624196*^9}}, CellLabel->"In[80]:=",ExpressionUUID->"f01799d3-d019-42fd-bcfe-c3836f30df57"], Cell[BoxData[ TemplateBox[{ "NSolve", "ifun", "\"Inverse functions are being used by \ \\!\\(\\*RowBox[{\\\"NSolve\\\"}]\\), so some solutions may not be found; use \ Reduce for complete solution information.\"", 2, 80, 16, 26402770535264617284, "Local"}, "MessageTemplate"]], "Message", "MSG", CellChangeTimes->{3.870144496668459*^9}, CellLabel-> "During evaluation of \ In[80]:=",ExpressionUUID->"4f280a88-130e-41c1-aa29-c30aeeab0b71"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"aa", "\[Rule]", RowBox[{"-", "1.1434403243410791`"}]}], ",", RowBox[{"dd", "\[Rule]", RowBox[{"-", "1.811811372234778`"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"aa", "\[Rule]", "1.1434403243410791`"}], ",", RowBox[{"dd", "\[Rule]", "1.3297812813550152`"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.8701444966805964`*^9}, CellLabel->"Out[80]=",ExpressionUUID->"c02235fe-7202-4083-a1a7-5c1cef001eb7"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"NSolve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"fIn", "[", RowBox[{"aa", ",", "4"}], "]"}], "==", RowBox[{"fOut", "[", RowBox[{"dd", ",", "4"}], "]"}]}], ",", RowBox[{ RowBox[{"fIn1", "[", RowBox[{"aa", ",", "4", ",", FractionBox["kIn", "197.327"], ",", "range"}], "]"}], "==", RowBox[{"fOut1", "[", RowBox[{"dd", ",", "4", ",", FractionBox["kOut", "197.327"], ",", "range"}], "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"aa", ",", "dd"}], "}"}]}], "]"}]], "Input", CellChangeTimes->{{3.870144527416617*^9, 3.870144532663507*^9}}, CellLabel->"In[81]:=",ExpressionUUID->"20d99c70-d5e0-40e2-ae06-43ab0b254118"], Cell[BoxData[ TemplateBox[{ "NSolve", "ifun", "\"Inverse functions are being used by \ \\!\\(\\*RowBox[{\\\"NSolve\\\"}]\\), so some solutions may not be found; use \ Reduce for complete solution information.\"", 2, 81, 17, 26402770535264617284, "Local"}, "MessageTemplate"]], "Message", "MSG", CellChangeTimes->{3.870144533940234*^9}, CellLabel-> "During evaluation of \ In[81]:=",ExpressionUUID->"cf720439-438e-4f44-bb06-85b9c13e93d6"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"aa", "\[Rule]", RowBox[{"-", "4.403018165298792`"}]}], ",", RowBox[{"dd", "\[Rule]", RowBox[{"-", "2.851922695522528`"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"aa", "\[Rule]", "4.403018165298792`"}], ",", RowBox[{"dd", "\[Rule]", "0.28966995806726537`"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.87014453395214*^9}, CellLabel->"Out[81]=",ExpressionUUID->"2d577fe4-6ba2-485b-ac27-870aa1bdca53"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"NSolve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"fIn", "[", RowBox[{"aa", ",", "6"}], "]"}], "==", RowBox[{"fOut", "[", RowBox[{"dd", ",", "6"}], "]"}]}], ",", RowBox[{ RowBox[{"fIn1", "[", RowBox[{"aa", ",", "6", ",", FractionBox["kIn", "197.327"], ",", "range"}], "]"}], "==", RowBox[{"fOut1", "[", RowBox[{"dd", ",", "6", ",", FractionBox["kOut", "197.327"], ",", "range"}], "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"aa", ",", "dd"}], "}"}]}], "]"}]], "Input", CellChangeTimes->{{3.8701445405223303`*^9, 3.8701445529369593`*^9}}, CellLabel->"In[83]:=",ExpressionUUID->"7471309e-5ddf-4c6c-952f-ca3808195e99"], Cell[BoxData[ TemplateBox[{ "NSolve", "ifun", "\"Inverse functions are being used by \ \\!\\(\\*RowBox[{\\\"NSolve\\\"}]\\), so some solutions may not be found; use \ Reduce for complete solution information.\"", 2, 83, 19, 26402770535264617284, "Local"}, "MessageTemplate"]], "Message", "MSG", CellChangeTimes->{{3.8701445455037537`*^9, 3.870144553990568*^9}}, CellLabel-> "During evaluation of \ In[83]:=",ExpressionUUID->"952a5a99-dfb9-49d0-8501-dadec2d5ae08"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"aa", "\[Rule]", RowBox[{"-", "21.353155624646753`"}]}], ",", RowBox[{"dd", "\[Rule]", RowBox[{"-", "3.1255541542481433`"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"aa", "\[Rule]", "21.353155624646753`"}], ",", RowBox[{"dd", "\[Rule]", "0.016038499341649987`"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{{3.87014454551573*^9, 3.8701445540064983`*^9}}, CellLabel->"Out[83]=",ExpressionUUID->"e403747e-08d3-4041-a50a-f86d6a57de1e"] }, Open ]] }, WindowSize->{1389.75, 768.75}, WindowMargins->{{0, Automatic}, {0, Automatic}}, FrontEndVersion->"13.0 for Linux x86 (64-bit) (December 2, 2021)", StyleDefinitions->"Default.nb", ExpressionUUID->"ece92cc8-ee08-4698-89c1-d978b1813711" ] (* 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, 534, 14, 29, "Input",ExpressionUUID->"6f155ea1-2306-4704-a233-560f70b5d462"], Cell[CellGroupData[{ Cell[5937, 170, 339, 7, 29, "Input",ExpressionUUID->"785e23c6-adf0-4e04-99e7-8de8b70d6714"], Cell[6279, 179, 999, 26, 38, "Print",ExpressionUUID->"a2f6bfed-b9aa-420b-b27b-846a66fae3eb"] }, Open ]], Cell[7293, 208, 489, 13, 29, "Input",ExpressionUUID->"30b0d700-2070-478c-bc3c-d787b0cb0c36"], Cell[CellGroupData[{ Cell[7807, 225, 299, 7, 29, "Input",ExpressionUUID->"5f2128b6-cab8-40ec-bfa2-267ec247a695"], Cell[8109, 234, 1776, 47, 59, "Print",ExpressionUUID->"ec888d6c-c887-43e1-b6e3-22e6c4edb4ec"] }, Open ]], Cell[9900, 284, 311, 7, 29, "Input",ExpressionUUID->"b7fb5240-5d57-4dc2-a5cf-2ef2edd84a37"], Cell[10214, 293, 1141, 26, 77, "Input",ExpressionUUID->"b19389b2-33a5-43a0-ab12-563d4c8b3b62"], Cell[CellGroupData[{ Cell[11380, 323, 377, 8, 29, "Input",ExpressionUUID->"e09b75c5-e7ee-41f2-92ef-6d99a9250df2"], Cell[11760, 333, 623, 14, 25, "Print",ExpressionUUID->"fd042a95-017b-41dc-970a-caa921531275"] }, Open ]], Cell[12398, 350, 212, 4, 29, "Input",ExpressionUUID->"5b2110d3-922f-468b-93f7-55f20889fabc"], Cell[12613, 356, 412, 11, 35, "Input",ExpressionUUID->"51bcd62e-b916-47fb-b730-ed6a9ea5a103"], Cell[13028, 369, 662, 14, 29, "Input",ExpressionUUID->"a1d36f07-1393-4716-b2a5-7dce8c9b5c36"], Cell[CellGroupData[{ Cell[13715, 387, 331, 7, 29, "Input",ExpressionUUID->"05cf6d8c-1da7-4252-a566-4d8205b8e2d9"], Cell[14049, 396, 981, 25, 25, "Print",ExpressionUUID->"3e65c34a-a20d-49b9-a06a-b22c57e55b93"] }, Open ]], Cell[15045, 424, 1736, 42, 91, "Input",ExpressionUUID->"3c13eff2-f5aa-4a65-b4fe-877985d5db36"], Cell[CellGroupData[{ Cell[16806, 470, 466, 8, 46, "Input",ExpressionUUID->"3ce8683e-f553-40fc-b9dc-ca771f703664"], Cell[17275, 480, 906, 23, 23, "Print",ExpressionUUID->"71adf521-737e-4433-9e68-6d2a46404b8b"] }, Open ]], Cell[18196, 506, 2244, 58, 134, "Input",ExpressionUUID->"4e60cda4-270f-41c1-a4f4-78d03b351537"], Cell[CellGroupData[{ Cell[20465, 568, 371, 7, 46, "Input",ExpressionUUID->"107cc588-8aa8-49dd-a831-be09f57f8eda"], Cell[20839, 577, 1580, 42, 44, "Print",ExpressionUUID->"0a0cafd0-cdc3-47ee-924c-3fa7aeb43452"] }, Open ]], Cell[CellGroupData[{ Cell[22456, 624, 1130, 25, 46, "Input",ExpressionUUID->"2259a25f-7509-401a-afac-84d646e2336f"], Cell[23589, 651, 576, 13, 28, "Message",ExpressionUUID->"9244a4dd-143f-45ab-a98d-28422173269e"], Cell[24168, 666, 642, 17, 33, "Output",ExpressionUUID->"63a47552-f5f5-426b-be4c-082d226082c7"] }, Open ]], Cell[CellGroupData[{ Cell[24847, 688, 756, 20, 46, "Input",ExpressionUUID->"f01799d3-d019-42fd-bcfe-c3836f30df57"], Cell[25606, 710, 454, 11, 28, "Message",ExpressionUUID->"4f280a88-130e-41c1-aa29-c30aeeab0b71"], Cell[26063, 723, 544, 15, 33, "Output",ExpressionUUID->"c02235fe-7202-4083-a1a7-5c1cef001eb7"] }, Open ]], Cell[CellGroupData[{ Cell[26644, 743, 756, 20, 46, "Input",ExpressionUUID->"20d99c70-d5e0-40e2-ae06-43ab0b254118"], Cell[27403, 765, 454, 11, 28, "Message",ExpressionUUID->"cf720439-438e-4f44-bb06-85b9c13e93d6"], Cell[27860, 778, 540, 15, 33, "Output",ExpressionUUID->"2d577fe4-6ba2-485b-ac27-870aa1bdca53"] }, Open ]], Cell[CellGroupData[{ Cell[28437, 798, 760, 20, 46, "Input",ExpressionUUID->"7471309e-5ddf-4c6c-952f-ca3808195e99"], Cell[29200, 820, 480, 11, 28, "Message",ExpressionUUID->"952a5a99-dfb9-49d0-8501-dadec2d5ae08"], Cell[29683, 833, 570, 15, 33, "Output",ExpressionUUID->"e403747e-08d3-4041-a50a-f86d6a57de1e"] }, Open ]] } ] *)