Skip to content

M
mathematica
Commit d5a3c0bc authored by himyss's avatar himyss
Browse files
Options

added

parent d1d49668
Show whitespace changes
Inline Side-by-side
Showing with 16823 additions and 4142 deletions
bounded8He.nb
View file @ d5a3c0bc
......@@ -10,10 +10,10 @@
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 158, 7]
NotebookDataLength[ 51640, 1230]
NotebookOptionsPosition[ 48748, 1173]
NotebookOutlinePosition[ 49146, 1189]
CellTagsIndexPosition[ 49103, 1186]
NotebookDataLength[ 55741, 1330]
NotebookOptionsPosition[ 52847, 1273]
NotebookOutlinePosition[ 53245, 1289]
CellTagsIndexPosition[ 53202, 1286]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
......@@ -34,29 +34,25 @@ Cell[BoxData[{
RowBox[{"mass", " ", "of", " ", "neutron"}], ",", " ", "MeV"}],
"*)"}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"Ebind", "=", "3.925"}], " ", ";"}],
RowBox[{"Ebind1", "=", "3.925"}], " ", ";",
RowBox[{"(*",
RowBox[{
RowBox[{"binding", " ", "energy", " ",
SuperscriptBox[
RowBox[{"of", " "}], "8"], "He"}], ",", " ", "MeV"}],
"*)"}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"Esep", "=", "24.81432"}], ";"}], " ",
RowBox[{"of", " "}], "8"], "He"}], ",", " ", "MeV"}], "*)"}],
"\[IndentingNewLine]",
RowBox[{"Esep1", "=", "24.81432"}], ";", " ",
RowBox[{"(*",
RowBox[{
RowBox[{"1", "p", " ", "separation", " ", "energy", " ",
SuperscriptBox[
RowBox[{"for", " "}], "8"], "He"}], ",", " ", "MeV"}],
"*)"}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"for", " "}], "8"], "He"}], ",", " ", "MeV"}], "*)"}],
"\[IndentingNewLine]",
RowBox[{"p", "=", "197.327"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"mass", "=", "821"}], ";"}], "\[IndentingNewLine]",
RowBox[{"mass1", "=", "821"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"range", "=", "3.735"}], ";"}]}], "Input",
RowBox[{"range1", "=", "3.735"}], ";"}]}], "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,
......@@ -66,8 +62,9 @@ Cell[BoxData[{
3.879568558358831*^9, 3.879568558418344*^9}, {3.87956872817546*^9,
3.879568728576707*^9}, {3.880087700774041*^9, 3.880087701627767*^9}, {
3.8800877975749063`*^9, 3.880087987460205*^9}, 3.880694076159418*^9, {
3.880694106655357*^9, 3.88069421626336*^9}},
CellLabel->"In[2]:=",ExpressionUUID->"ad6df89e-d3d8-490c-ae91-dac6d80d1b5a"],
3.880694106655357*^9, 3.88069421626336*^9}, {3.881205548538989*^9,
3.881205552151064*^9}},
CellLabel->"In[1]:=",ExpressionUUID->"ad6df89e-d3d8-490c-ae91-dac6d80d1b5a"],
Cell[BoxData[{
RowBox[{
......@@ -106,7 +103,7 @@ Cell[BoxData[{
RowBox[{
RowBox[{"dfInR", "[",
RowBox[{"q", ",", "r", ",", "ang"}], "]"}], "/.",
RowBox[{"r", "->", "range"}]}], "]"}]}], "\[IndentingNewLine]",
RowBox[{"r", "->", "range1"}]}], "]"}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"dfOutR", "[",
RowBox[{"k_", ",", "r_", ",", "ang_"}], "]"}], ":=",
......@@ -131,7 +128,7 @@ Cell[BoxData[{
RowBox[{
RowBox[{"dfOutR", "[",
RowBox[{"q", ",", "r", ",", "ang"}], "]"}], "/.",
RowBox[{"r", "->", "range"}]}], "]"}]}]}], "Input",
RowBox[{"r", "->", "range1"}]}], "]"}]}]}], "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}, {
......@@ -142,8 +139,8 @@ Cell[BoxData[{
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[9]:=",ExpressionUUID->"2fc4ef7c-d275-441b-9d3d-7d7239f74130"],
3.8782888301901093`*^9}, {3.88120555964981*^9, 3.881205560418138*^9}},
CellLabel->"In[6]:=",ExpressionUUID->"2fc4ef7c-d275-441b-9d3d-7d7239f74130"],
Cell[CellGroupData[{
......@@ -156,33 +153,33 @@ Cell[BoxData[
RowBox[{
FractionBox[
SqrtBox[
RowBox[{"2", " ", "mass", " ",
RowBox[{"2", " ", "mass1", " ",
RowBox[{"(",
RowBox[{
RowBox[{"-", "Esep"}], "+", "U"}], ")"}]}]], "p"], ",", "range",
RowBox[{"-", "Esep1"}], "+", "U"}], ")"}]}]], "p"], ",", "range1",
",", "0"}], "]"}],
RowBox[{"dfIn", "[",
RowBox[{
FractionBox[
SqrtBox[
RowBox[{"2", " ", "mass", " ",
RowBox[{"2", " ", "mass1", " ",
RowBox[{"(",
RowBox[{
RowBox[{"-", "Esep"}], "+", "U"}], ")"}]}]], "p"], ",", "0"}],
RowBox[{"-", "Esep1"}], "+", "U"}], ")"}]}]], "p"], ",", "0"}],
"]"}]], "-",
FractionBox[
RowBox[{"fOut", "[",
RowBox[{
FractionBox[
SqrtBox[
RowBox[{"2", " ", "mass", " ", "Esep"}]], "p"], ",", "range", ",",
RowBox[{"2", " ", "mass1", " ", "Esep1"}]], "p"], ",", "range1", ",",
"0"}], "]"}],
RowBox[{"dfOut", "[",
RowBox[{
FractionBox[
SqrtBox[
RowBox[{"2", " ", "mass", " ", "Esep"}]], "p"], ",", "0"}], "]"}]]}],
",",
RowBox[{"2", " ", "mass1", " ", "Esep1"}]], "p"], ",", "0"}],
"]"}]]}], ",",
RowBox[{"{",
RowBox[{"U", ",", "0", ",", "50"}], "}"}]}], "]"}]], "Input",
CellChangeTimes->{{3.8762201382091026`*^9, 3.8762201688017263`*^9}, {
......@@ -199,8 +196,9 @@ Cell[BoxData[
3.878288877525551*^9, 3.8782889008012238`*^9}, {3.878288973320475*^9,
3.878289061721037*^9}, {3.8795689049227552`*^9, 3.87956893336374*^9}, {
3.879568964437098*^9, 3.879568981740798*^9}, {3.879569036636216*^9,
3.879569049028247*^9}, 3.880087707720909*^9, 3.880860838537876*^9},
CellLabel->"In[15]:=",ExpressionUUID->"c2ad0be7-22ee-4300-8ad8-45a6e113c1b0"],
3.879569049028247*^9}, 3.880087707720909*^9, 3.880860838537876*^9, {
3.881205565537171*^9, 3.881205582259083*^9}},
CellLabel->"In[12]:=",ExpressionUUID->"c2ad0be7-22ee-4300-8ad8-45a6e113c1b0"],
Cell[BoxData[
GraphicsBox[{{{}, {},
......@@ -299,7 +297,7 @@ B42U8787ZolIMDEOzZjtpZFX4b681hUSuvw/tTcN06ioJHBpdZWEgecX9B6O
"]],
LineBox[{{28.992520401116344`, 6.989114408516514}, {
28.992522325141426`, -3.5027620306546083`}}]},
Annotation[#, "Charting`Private`Tag$12001#1"]& ]}, {}},
Annotation[#, "Charting`Private`Tag$4194#1"]& ]}, {}},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->{True, True},
AxesLabel->{None, None},
......@@ -351,14 +349,15 @@ B42U8787ZolIMDEOzZjtpZFX4b681hUSuvw/tTcN06ioJHBpdZWEgecX9B6O
3.879569032952899*^9, 3.8795690542179728`*^9}, {3.88008674094413*^9,
3.880086745707924*^9}, 3.880087709857608*^9, {3.880087803643498*^9,
3.880087992161138*^9}, 3.88069400745238*^9, {3.88086082362481*^9,
3.880860840824017*^9}, 3.880869126504476*^9},
CellLabel->"Out[15]=",ExpressionUUID->"fc3a6f0e-96cd-440d-9384-66ff4e3f2956"]
3.880860840824017*^9}, 3.880869126504476*^9, 3.881204831143188*^9,
3.881205584395485*^9},
CellLabel->"Out[12]=",ExpressionUUID->"ab87c9e7-d22d-479e-8e60-a486d3571fb7"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"myU", "=",
RowBox[{"myU1", "=",
RowBox[{"U", "/.",
RowBox[{"FindRoot", "[",
RowBox[{
......@@ -369,32 +368,32 @@ Cell[BoxData[
RowBox[{
FractionBox[
SqrtBox[
RowBox[{"2", " ", "mass", " ",
RowBox[{"2", " ", "mass1", " ",
RowBox[{"(",
RowBox[{
RowBox[{"-", "Esep"}], "+", "U"}], ")"}]}]], "p"], ",",
"range", ",", "0"}], "]"}],
RowBox[{"-", "Esep1"}], "+", "U"}], ")"}]}]], "p"], ",",
"range1", ",", "0"}], "]"}],
RowBox[{"dfIn", "[",
RowBox[{
FractionBox[
SqrtBox[
RowBox[{"2", " ", "mass", " ",
RowBox[{"2", " ", "mass1", " ",
RowBox[{"(",
RowBox[{
RowBox[{"-", "Esep"}], "+", "U"}], ")"}]}]], "p"], ",", "0"}],
RowBox[{"-", "Esep1"}], "+", "U"}], ")"}]}]], "p"], ",", "0"}],
"]"}]], "-",
FractionBox[
RowBox[{"fOut", "[",
RowBox[{
FractionBox[
SqrtBox[
RowBox[{"2", " ", "mass", " ", "Esep"}]], "p"], ",", "range", ",",
"0"}], "]"}],
RowBox[{"2", " ", "mass1", " ", "Esep1"}]], "p"], ",", "range1",
",", "0"}], "]"}],
RowBox[{"dfOut", "[",
RowBox[{
FractionBox[
SqrtBox[
RowBox[{"2", " ", "mass", " ", "Esep"}]], "p"], ",", "0"}],
RowBox[{"2", " ", "mass1", " ", "Esep1"}]], "p"], ",", "0"}],
"]"}]]}], "\[Equal]", "0"}], ",",
RowBox[{"{",
RowBox[{"U", ",", "60"}], "}"}]}], "]"}]}]}]], "Input",
......@@ -405,8 +404,9 @@ Cell[BoxData[
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.880086673187624*^9, 3.8800877134548693`*^9, 3.880087807038721*^9},
CellLabel->"In[16]:=",ExpressionUUID->"22f1491b-05f8-4461-9d8a-7023940b0eba"],
3.880086673187624*^9, 3.8800877134548693`*^9, 3.880087807038721*^9, {
3.8812055878679543`*^9, 3.8812056017959967`*^9}},
CellLabel->"In[13]:=",ExpressionUUID->"22f1491b-05f8-4461-9d8a-7023940b0eba"],
Cell[BoxData["35.89034438287419`"], "Output",
CellChangeTimes->{
......@@ -420,14 +420,15 @@ Cell[BoxData["35.89034438287419`"], "Output",
3.880086743953726*^9, 3.8800867473357487`*^9}, 3.880087715279483*^9, {
3.880087808333634*^9, 3.880087992954145*^9}, 3.880694008901841*^9,
3.88069408135425*^9, {3.88069411322472*^9, 3.880694222245901*^9},
3.880860843173624*^9, 3.8808691277708607`*^9},
CellLabel->"Out[16]=",ExpressionUUID->"43c6b393-af8b-4903-9f47-817769bc3c06"]
3.880860843173624*^9, 3.8808691277708607`*^9, 3.8812048322740593`*^9,
3.881205604361702*^9},
CellLabel->"Out[13]=",ExpressionUUID->"d0b79513-3f3c-4a7e-9d4c-9cf52de46f51"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"myCoeff", "=",
RowBox[{"myCoeff1", "=",
RowBox[{"coeff", "/.",
RowBox[{"FindRoot", "[",
RowBox[{
......@@ -437,38 +438,40 @@ Cell[BoxData[
RowBox[{
FractionBox[
SqrtBox[
RowBox[{"2", " ", "mass", " ",
RowBox[{"2", " ", "mass1", " ",
RowBox[{"(",
RowBox[{
RowBox[{"-", "Esep"}], "+", "myU"}], ")"}]}]], "p"], ",",
"range", ",", "0"}], "]"}], "-",
RowBox[{"-", "Esep1"}], "+", "myU1"}], ")"}]}]], "p"], ",",
"range1", ",", "0"}], "]"}], "-",
RowBox[{"coeff", " ", "*",
RowBox[{"fOut", "[",
RowBox[{
FractionBox[
SqrtBox[
RowBox[{"2", " ", "mass", " ", "Esep"}]], "p"], ",", "range", ",",
"0"}], "]"}]}]}], "\[Equal]", "0"}], ",",
RowBox[{"2", " ", "mass1", " ", "Esep1"}]], "p"], ",", "range1",
",", "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.8800866773952007`*^9, 3.880086704667201*^9}},
CellLabel->"In[17]:=",ExpressionUUID->"444da03a-39cf-4c3d-bb1b-9d322299f630"],
3.8800866773952007`*^9, 3.880086704667201*^9}, {3.881205606847561*^9,
3.881205617689822*^9}},
CellLabel->"In[14]:=",ExpressionUUID->"444da03a-39cf-4c3d-bb1b-9d322299f630"],
Cell[BoxData["25.351749791847226`"], "Output",
CellChangeTimes->{
3.879570442808729*^9, 3.87957064413636*^9, 3.8800867486475277`*^9,
3.880087717563438*^9, {3.880087810097796*^9, 3.880087994342682*^9},
3.880694010579237*^9, 3.8806940829871798`*^9, {3.8806941143549423`*^9,
3.88069422353505*^9}, 3.880860847232801*^9, 3.88086912945323*^9},
CellLabel->"Out[17]=",ExpressionUUID->"c1ef58b4-3b2f-4db8-b274-f079d4bf1a4d"]
3.88069422353505*^9}, 3.880860847232801*^9, 3.88086912945323*^9,
3.881204833642619*^9, 3.8812056191396523`*^9},
CellLabel->"Out[14]=",ExpressionUUID->"f3d3808a-1e3a-4a37-b9ac-83ad7cab0a61"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"myNorm", "=",
RowBox[{"myNorm1", "=",
RowBox[{"A", "/.",
RowBox[{"FindRoot", "[",
RowBox[{
......@@ -483,17 +486,17 @@ Cell[BoxData[
RowBox[{
FractionBox[
SqrtBox[
RowBox[{"2", " ", "mass", " ",
RowBox[{"2", " ", "mass1", " ",
RowBox[{"(",
RowBox[{
RowBox[{"-", "Esep"}], "+", "myU"}], ")"}]}]], "p"], ",",
RowBox[{"-", "Esep1"}], "+", "myU1"}], ")"}]}]], "p"], ",",
"r", ",", "0"}], "]"}], "2"], ",", " ",
RowBox[{"{",
RowBox[{"r", ",", "0", ",", "range"}], "}"}]}], "]"}]}], "+",
RowBox[{"r", ",", "0", ",", "range1"}], "}"}]}], "]"}]}], "+",
RowBox[{
SuperscriptBox[
RowBox[{"(",
RowBox[{"myCoeff", "*", "A"}], ")"}], "2"], " ",
RowBox[{"myCoeff1", "*", "A"}], ")"}], "2"], " ",
RowBox[{"Integrate", "[",
RowBox[{
SuperscriptBox[
......@@ -501,10 +504,10 @@ Cell[BoxData[
RowBox[{
FractionBox[
SqrtBox[
RowBox[{"2", " ", "mass", " ", "Esep"}]], "p"], ",", "r", ",",
"0"}], "]"}], "2"], ",", " ",
RowBox[{"2", " ", "mass1", " ", "Esep1"}]], "p"], ",", "r",
",", "0"}], "]"}], "2"], ",", " ",
RowBox[{"{",
RowBox[{"r", ",", "range", ",", "Infinity"}], "}"}]}], "]"}]}]}],
RowBox[{"r", ",", "range1", ",", "Infinity"}], "}"}]}], "]"}]}]}],
"\[Equal]", "1"}], ",", " ",
RowBox[{"{",
RowBox[{"A", ",", "0.2"}], "}"}]}], "]"}]}]}]], "Input",
......@@ -514,8 +517,9 @@ Cell[BoxData[
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.8800866795052*^9, 3.880086723648281*^9}},
CellLabel->"In[18]:=",ExpressionUUID->"7835cc6b-0406-456f-a3a3-ab6f89b509c7"],
3.8800866795052*^9, 3.880086723648281*^9}, {3.881205621463208*^9,
3.881205632051371*^9}},
CellLabel->"In[15]:=",ExpressionUUID->"7835cc6b-0406-456f-a3a3-ab6f89b509c7"],
Cell[BoxData["0.6514573663189586`"], "Output",
CellChangeTimes->{
......@@ -523,8 +527,9 @@ Cell[BoxData["0.6514573663189586`"], "Output",
3.879570637910142*^9, 3.879570646893675*^9}, 3.8800867527823553`*^9,
3.880087721464933*^9, {3.880087811800939*^9, 3.8800879962360888`*^9},
3.880694013817318*^9, 3.8806940857627296`*^9, {3.8806941160551243`*^9,
3.880694225103177*^9}, 3.88086085054058*^9, 3.8808691316909847`*^9},
CellLabel->"Out[18]=",ExpressionUUID->"948665ac-1376-42e6-910e-4a23a2bb0b60"]
3.880694225103177*^9}, 3.88086085054058*^9, 3.8808691316909847`*^9,
3.881204835578549*^9, 3.881205634171431*^9},
CellLabel->"Out[15]=",ExpressionUUID->"413f2204-d5e7-4ebc-99c1-5d089819350a"]
}, Open ]],
Cell[CellGroupData[{
......@@ -537,33 +542,34 @@ Cell[BoxData[
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"myNorm", " ",
RowBox[{"myNorm1", " ",
RowBox[{"fIn", "[",
RowBox[{
FractionBox[
SqrtBox[
RowBox[{"2", " ", "mass", " ",
RowBox[{"2", " ", "mass1", " ",
RowBox[{"(",
RowBox[{
RowBox[{"-", "Esep"}], "+", "myU"}], ")"}]}]], "p"], ",", "r",
",", "0"}], "]"}]}], ",",
RowBox[{"r", "<", "range"}]}], "}"}], ",",
RowBox[{"-", "Esep1"}], "+", "myU1"}], ")"}]}]], "p"], ",",
"r", ",", "0"}], "]"}]}], ",",
RowBox[{"r", "<", "range1"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{
RowBox[{"myNorm", " ", "myCoeff", " ",
RowBox[{"myNorm1", " ", "myCoeff1", " ",
RowBox[{"fOut", "[",
RowBox[{
FractionBox[
SqrtBox[
RowBox[{"2", " ", "mass", " ", "Esep"}]], "p"], ",", "r", ",",
RowBox[{"2", " ", "mass1", " ", "Esep1"}]], "p"], ",", "r", ",",
"0"}], "]"}]}], ",",
RowBox[{"r", ">", "range"}]}], "}"}]}], "}"}], "]"}], ",",
RowBox[{"r", ">", "range1"}]}], "}"}]}], "}"}], "]"}], ",",
RowBox[{"{",
RowBox[{"r", ",", "0", ",", "10"}], "}"}]}], "]"}]], "Input",
CellChangeTimes->{{3.8771831750308676`*^9, 3.877183283833062*^9}, {
3.879570725415313*^9, 3.87957082556467*^9}, {3.8800866815822*^9,
3.88008672825279*^9}, 3.880086759707464*^9},
CellLabel->"In[19]:=",ExpressionUUID->"fc251146-842e-4797-bd33-46971b36a427"],
3.88008672825279*^9}, 3.880086759707464*^9, {3.88120563620135*^9,
3.881205647216186*^9}},
CellLabel->"In[16]:=",ExpressionUUID->"fc251146-842e-4797-bd33-46971b36a427"],
Cell[BoxData[
GraphicsBox[{{{}, {},
......@@ -722,7 +728,7 @@ oZnDCOJWTgyG/WRiRtC1zT0/DCGnoiRpfp6JLjk+SSK9hlAczeFTs8jEQ2aP
a82fGMIta+/tEX+YKJ2qPpZ1yxDuM1rlDf4xkT4aSMzkGsLTWeWF5WUm2sWZ
NDieNYT/AXAvkZQ=
"]]},
Annotation[#, "Charting`Private`Tag$14136#1"]& ]}, {}},
Annotation[#, "Charting`Private`Tag$6329#1"]& ]}, {}},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->{True, True},
AxesLabel->{None, None},
......@@ -766,8 +772,9 @@ NDieNYT/AXAvkZQ=
3.880086754549075*^9, 3.8800867612813787`*^9}, 3.880087724555195*^9, {
3.880087813143466*^9, 3.880087997646585*^9}, 3.880694015748687*^9,
3.880694087181596*^9, {3.880694118038547*^9, 3.8806942262383127`*^9},
3.880860852002902*^9, 3.880869132289547*^9},
CellLabel->"Out[19]=",ExpressionUUID->"512149a4-22e4-4c08-9068-29ce9e398f67"]
3.880860852002902*^9, 3.880869132289547*^9, 3.881204836375863*^9,
3.8812056496910753`*^9},
CellLabel->"Out[16]=",ExpressionUUID->"8082d315-b512-4039-9ab9-4f8290a4ac22"]
}, Open ]],
Cell[BoxData[
......@@ -792,43 +799,45 @@ Cell[BoxData[
RowBox[{"Integrate", "[",
RowBox[{
RowBox[{
SuperscriptBox["myNorm", "2"],
SuperscriptBox["myNorm1", "2"],
SuperscriptBox[
RowBox[{"fIn", "[",
RowBox[{
FractionBox[
SqrtBox[
RowBox[{"2", " ", "mass", " ",
RowBox[{"2", " ", "mass1", " ",
RowBox[{"(",
RowBox[{
RowBox[{"-", "Esep"}], "+", "myU"}], ")"}]}]], "p"], ",", "r",
",", "0"}], "]"}], "2"],
RowBox[{"-", "Esep1"}], "+", "myU1"}], ")"}]}]], "p"], ",",
"r", ",", "0"}], "]"}], "2"],
SuperscriptBox["r", "2"]}], ",",
RowBox[{"{",
RowBox[{"r", ",", "0", ",", "range"}], "}"}]}], "]"}], "+",
RowBox[{"r", ",", "0", ",", "range1"}], "}"}]}], "]"}], "+",
RowBox[{"Integrate", "[",
RowBox[{
RowBox[{
SuperscriptBox["myNorm", "2"], " ",
SuperscriptBox["myCoeff", "2"], " ",
SuperscriptBox["myNorm1", "2"], " ",
SuperscriptBox["myCoeff1", "2"], " ",
SuperscriptBox[
RowBox[{"fOut", "[",
RowBox[{
FractionBox[
SqrtBox[
RowBox[{"2", " ", "mass", " ", "Esep"}]], "p"], ",", "r", ",",
RowBox[{"2", " ", "mass1", " ", "Esep1"}]], "p"], ",", "r", ",",
"0"}], "]"}], "2"],
SuperscriptBox["r", "2"]}], ",",
RowBox[{"{",
RowBox[{"r", ",", "range", ",", "Infinity"}], "}"}]}], "]"}]}],
RowBox[{"r", ",", "range1", ",", "Infinity"}], "}"}]}], "]"}]}],
"]"}]}]], "Input",
CellLabel->"In[20]:=",ExpressionUUID->"6db522eb-96cb-477e-99c5-d27334b80d0b"],
CellChangeTimes->{{3.881205652720017*^9, 3.8812056628772306`*^9}},
CellLabel->"In[17]:=",ExpressionUUID->"6db522eb-96cb-477e-99c5-d27334b80d0b"],
Cell[BoxData["1.6819723510511688`"], "Output",
CellChangeTimes->{
3.880694061225527*^9, {3.880694092576046*^9, 3.8806942478489323`*^9},
3.88086085882454*^9, 3.880869134852079*^9},
CellLabel->"Out[20]=",ExpressionUUID->"96772995-7f2f-4004-8007-b9946c1988b9"]
3.88086085882454*^9, 3.880869134852079*^9, 3.881204839005999*^9,
3.881205664575473*^9},
CellLabel->"Out[17]=",ExpressionUUID->"f0141514-31a9-41e5-9f3f-c4d639d78a0e"]
}, Open ]],
Cell[CellGroupData[{
......@@ -1001,37 +1010,44 @@ Cell[BoxData[
RowBox[{"Sin", "[",
RowBox[{"a", " ", "x"}], "]"}], " ",
RowBox[{"Sin", "[",
RowBox[{"b", " ", "x"}], "]"}]}], ",",
RowBox[{"{",
RowBox[{"x", ",", "0", ",", "r"}], "}"}], ",",
RowBox[{"b", " ", "x"}], "]"}]}], ",", "x", ",",
RowBox[{"Assumptions", "\[Rule]",
RowBox[{
RowBox[{"a", ">", "0"}], " ", "&&", " ",
RowBox[{"b", ">", "0"}]}]}]}], "]"}]], "Input",
CellChangeTimes->{{3.880861521929137*^9, 3.8808615332101803`*^9}},
CellLabel->"In[22]:=",ExpressionUUID->"0773e3b1-b600-49c8-a187-fa642d2580e8"],
CellChangeTimes->{{3.880861521929137*^9, 3.8808615332101803`*^9}, {
3.88120377534515*^9, 3.881203786644293*^9}},
CellLabel->
"In[134]:=",ExpressionUUID->"0773e3b1-b600-49c8-a187-fa642d2580e8"],
Cell[BoxData[
FractionBox[
RowBox[{
RowBox[{"b", " ",
RowBox[{"Cos", "[",
RowBox[{"b", " ", "r"}], "]"}], " ",
FractionBox[
RowBox[{"Sin", "[",
RowBox[{"a", " ", "r"}], "]"}]}], "-",
RowBox[{"a", " ",
RowBox[{"Cos", "[",
RowBox[{"a", " ", "r"}], "]"}], " ",
RowBox[{
RowBox[{"(",
RowBox[{"a", "-", "b"}], ")"}], " ", "x"}], "]"}],
RowBox[{"2", " ",
RowBox[{"(",
RowBox[{"a", "-", "b"}], ")"}]}]], "-",
FractionBox[
RowBox[{"Sin", "[",
RowBox[{"b", " ", "r"}], "]"}]}]}],
RowBox[{
SuperscriptBox["a", "2"], "-",
SuperscriptBox["b", "2"]}]]], "Output",
CellChangeTimes->{3.8808615368747*^9, 3.880869149783291*^9},
CellLabel->"Out[22]=",ExpressionUUID->"8b590bfb-da17-4e96-815c-5cfc0a0f84da"]
RowBox[{"(",
RowBox[{"a", "+", "b"}], ")"}], " ", "x"}], "]"}],
RowBox[{"2", " ",
RowBox[{"(",
RowBox[{"a", "+", "b"}], ")"}]}]]}]], "Output",
CellChangeTimes->{
3.8808615368747*^9, 3.880869149783291*^9, {3.8812037807260647`*^9,
3.881203787872303*^9}},
CellLabel->
"Out[134]=",ExpressionUUID->"cc9562cb-8b47-4ffb-a0f0-a0db873c2346"]
}, Open ]],
Cell[BoxData[
RowBox[{
RowBox[{"(*",
RowBox[{
RowBox[{"PsiP", "[", "q_", "]"}], ":=", " ",
RowBox[{
......@@ -1121,13 +1137,98 @@ Cell[BoxData[
FractionBox[
SqrtBox[
RowBox[{"2", " ", "mass", " ", "Esep"}]], "p"], ")"}],
"2"]}]]}]}], ")"}]}]}]], "Input",
"2"]}]]}]}], ")"}]}]}], "*)"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"PsiP", "[", "q_", "]"}], ":=", " ",
RowBox[{
SqrtBox[
FractionBox["2", "\[Pi]"]], "myNorm1", " ",
RowBox[{"(",
RowBox[{
FractionBox[
RowBox[{"Sin", "[",
RowBox[{
RowBox[{"(",
RowBox[{"q", "-",
FractionBox[
SqrtBox[
RowBox[{"2", " ", "mass1", " ",
RowBox[{"(",
RowBox[{
RowBox[{"-", "Esep1"}], "+", "myU1"}], ")"}]}]], "p"]}],
")"}], " ", "range1"}], "]"}],
RowBox[{"2", " ",
RowBox[{"(",
RowBox[{"q", "-",
FractionBox[
SqrtBox[
RowBox[{"2", " ", "mass1", " ",
RowBox[{"(",
RowBox[{
RowBox[{"-", "Esep1"}], "+", "myU1"}], ")"}]}]], "p"]}],
")"}]}]], "-",
FractionBox[
RowBox[{"Sin", "[",
RowBox[{
RowBox[{"(",
RowBox[{"q", "+",
FractionBox[
SqrtBox[
RowBox[{"2", " ", "mass1", " ",
RowBox[{"(",
RowBox[{
RowBox[{"-", "Esep1"}], "+", "myU1"}], ")"}]}]], "p"]}],
")"}], " ", "range1"}], "]"}],
RowBox[{"2", " ",
RowBox[{"(",
RowBox[{"q", "+",
FractionBox[
SqrtBox[
RowBox[{"2", " ", "mass1", " ",
RowBox[{"(",
RowBox[{
RowBox[{"-", "Esep1"}], "+", "myU1"}], ")"}]}]], "p"]}],
")"}]}]], "+",
RowBox[{"myCoeff1",
FractionBox[
RowBox[{
SuperscriptBox["\[ExponentialE]",
RowBox[{
RowBox[{"-",
FractionBox[
SqrtBox[
RowBox[{"2", " ", "mass1", " ", "Esep1"}]], "p"]}], " ",
"range1"}]], " ",
RowBox[{"(",
RowBox[{
RowBox[{"q", " ",
RowBox[{"Cos", "[",
RowBox[{"q", " ", "range1"}], "]"}]}], "+",
RowBox[{
FractionBox[
SqrtBox[
RowBox[{"2", " ", "mass1", " ", "Esep1"}]], "p"], " ",
RowBox[{"Sin", "[",
RowBox[{"q", " ", "range1"}], "]"}]}]}], ")"}]}],
RowBox[{
SuperscriptBox["q", "2"], "+",
SuperscriptBox[
RowBox[{"(",
FractionBox[
SqrtBox[
RowBox[{"2", " ", "mass1", " ", "Esep1"}]], "p"], ")"}],
"2"]}]]}]}], ")"}]}]}]}]], "Input",
CellChangeTimes->{{3.880861581160618*^9, 3.880861704887731*^9},
3.8808617419572563`*^9, {3.880861875265856*^9, 3.8808619034776497`*^9}, {
3.880861942515173*^9, 3.8808619704782467`*^9}, {3.880875342255783*^9,
3.88087538894692*^9}, {3.881035003112468*^9,
3.881035013631791*^9}},ExpressionUUID->"c77c7be5-a519-41e6-b069-\
691e1006d400"],
3.88087538894692*^9}, {3.881035003112468*^9, 3.881035013631791*^9}, {
3.881203805263261*^9, 3.881203818283866*^9}, {3.881203924377963*^9,
3.881203938033923*^9}, {3.881204080636133*^9, 3.881204114844595*^9}, {
3.8812044164683723`*^9, 3.881204430878365*^9}, {3.8812047385493*^9,
3.881204786636528*^9}, {3.8812048642166367`*^9, 3.881204865174282*^9}, {
3.881204917757641*^9, 3.8812049207310743`*^9}, {3.8812056704257107`*^9,
3.881205704444508*^9}},
CellLabel->"In[18]:=",ExpressionUUID->"c77c7be5-a519-41e6-b069-691e1006d400"],
Cell[CellGroupData[{
......@@ -1146,8 +1247,7 @@ Cell[BoxData[
RowBox[{"(",
RowBox[{"PsiP", "[", "q", "]"}], ")"}], "2"], ",",
RowBox[{"{",
RowBox[{"q", ",", "0", ",", "Infinity"}], "}"}], ",",
RowBox[{"AccuracyGoal", "\[Rule]", "4"}]}], "]"}]}]], "Input",
RowBox[{"q", ",", "0", ",", "Infinity"}], "}"}]}], "]"}]}]], "Input",
CellChangeTimes->{{3.880862001847947*^9, 3.880862029306449*^9}, {
3.8808621852190742`*^9, 3.8808621872400637`*^9}, {3.8808622612561007`*^9,
3.8808622622720737`*^9}, {3.880869441370618*^9, 3.880869443208253*^9}, {
......@@ -1159,16 +1259,16 @@ Cell[BoxData[
3.880874814478899*^9, 3.880874818483645*^9}, {3.8808752703279543`*^9,
3.880875279874785*^9}, {3.880875321982376*^9, 3.880875327034482*^9}, {
3.880875469296236*^9, 3.880875498161496*^9}, {3.8808755743169622`*^9,
3.880875587481934*^9}},
CellLabel->
"In[101]:=",ExpressionUUID->"4e5819f2-d11b-4bce-b884-4cd47410fbd4"],
3.880875587481934*^9}, {3.881204851540709*^9, 3.88120487584007*^9}, {
3.881204929685568*^9, 3.881204938728588*^9}},
CellLabel->"In[19]:=",ExpressionUUID->"4e5819f2-d11b-4bce-b884-4cd47410fbd4"],
Cell[BoxData["0.9999996279403691`"], "Output",
Cell[BoxData["0.9999998739674564`"], "Output",
CellChangeTimes->{{3.880875317306684*^9, 3.8808753918008947`*^9}, {
3.880875476926515*^9, 3.880875499751493*^9}, {3.8808755895123568`*^9,
3.880875593014206*^9}},
CellLabel->
"Out[101]=",ExpressionUUID->"03549400-50ed-49fa-b63f-eb2e00cbf061"]
3.880875593014206*^9}, {3.881204848107025*^9, 3.881204877006847*^9}, {
3.881204927752026*^9, 3.881204939243373*^9}, 3.8812057081188097`*^9},
CellLabel->"Out[19]=",ExpressionUUID->"b9b74f50-8c3b-4f99-adc7-2201a234e737"]
}, Open ]]
},
WindowSize->{1389.75, 768.75},
......@@ -1188,49 +1288,49 @@ CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 1962, 49, 154, "Input",ExpressionUUID->"ad6df89e-d3d8-490c-ae91-dac6d80d1b5a"],
Cell[2523, 71, 2801, 74, 179, "Input",ExpressionUUID->"2fc4ef7c-d275-441b-9d3d-7d7239f74130"],
Cell[558, 20, 1962, 46, 154, "Input",ExpressionUUID->"ad6df89e-d3d8-490c-ae91-dac6d80d1b5a"],
Cell[2523, 68, 2848, 74, 179, "Input",ExpressionUUID->"2fc4ef7c-d275-441b-9d3d-7d7239f74130"],
Cell[CellGroupData[{
Cell[5349, 149, 2249, 53, 75, "Input",ExpressionUUID->"c2ad0be7-22ee-4300-8ad8-45a6e113c1b0"],
Cell[7601, 204, 8104, 150, 239, "Output",ExpressionUUID->"fc3a6f0e-96cd-440d-9384-66ff4e3f2956"]
Cell[5396, 146, 2312, 54, 75, "Input",ExpressionUUID->"c2ad0be7-22ee-4300-8ad8-45a6e113c1b0"],
Cell[7711, 202, 8151, 151, 239, "Output",ExpressionUUID->"ab87c9e7-d22d-479e-8e60-a486d3571fb7"]
}, Open ]],
Cell[CellGroupData[{
Cell[15742, 359, 1902, 49, 75, "Input",ExpressionUUID->"22f1491b-05f8-4461-9d8a-7023940b0eba"],
Cell[17647, 410, 931, 13, 33, "Output",ExpressionUUID->"43c6b393-af8b-4903-9f47-817769bc3c06"]
Cell[15899, 358, 1967, 50, 75, "Input",ExpressionUUID->"22f1491b-05f8-4461-9d8a-7023940b0eba"],
Cell[17869, 410, 981, 14, 33, "Output",ExpressionUUID->"d0b79513-3f3c-4a7e-9d4c-9cf52de46f51"]
}, Open ]],
Cell[CellGroupData[{
Cell[18615, 428, 1042, 28, 53, "Input",ExpressionUUID->"444da03a-39cf-4c3d-bb1b-9d322299f630"],
Cell[19660, 458, 431, 6, 33, "Output",ExpressionUUID->"c1ef58b4-3b2f-4db8-b274-f079d4bf1a4d"]
Cell[18887, 429, 1100, 29, 53, "Input",ExpressionUUID->"444da03a-39cf-4c3d-bb1b-9d322299f630"],
Cell[19990, 460, 481, 7, 33, "Output",ExpressionUUID->"f3d3808a-1e3a-4a37-b9ac-83ad7cab0a61"]
}, Open ]],
Cell[CellGroupData[{
Cell[20128, 469, 1903, 48, 53, "Input",ExpressionUUID->"7835cc6b-0406-456f-a3a3-ab6f89b509c7"],
Cell[22034, 519, 511, 7, 33, "Output",ExpressionUUID->"948665ac-1376-42e6-910e-4a23a2bb0b60"]
Cell[20508, 472, 1962, 49, 82, "Input",ExpressionUUID->"7835cc6b-0406-456f-a3a3-ab6f89b509c7"],
Cell[22473, 523, 559, 8, 33, "Output",ExpressionUUID->"413f2204-d5e7-4ebc-99c1-5d089819350a"]
}, Open ]],
Cell[CellGroupData[{
Cell[22582, 531, 1243, 34, 53, "Input",ExpressionUUID->"fc251146-842e-4797-bd33-46971b36a427"],
Cell[23828, 567, 11200, 202, 237, "Output",ExpressionUUID->"512149a4-22e4-4c08-9068-29ce9e398f67"]
Cell[23069, 536, 1302, 35, 53, "Input",ExpressionUUID->"fc251146-842e-4797-bd33-46971b36a427"],
Cell[24374, 573, 11249, 203, 237, "Output",ExpressionUUID->"8082d315-b512-4039-9ab9-4f8290a4ac22"]
}, Open ]],
Cell[35043, 772, 392, 10, 29, "Input",ExpressionUUID->"3894cc59-ff84-4f88-8130-0b3a87025770"],
Cell[35638, 779, 392, 10, 29, "Input",ExpressionUUID->"3894cc59-ff84-4f88-8130-0b3a87025770"],
Cell[CellGroupData[{
Cell[35460, 786, 1239, 38, 53, "Input",ExpressionUUID->"6db522eb-96cb-477e-99c5-d27334b80d0b"],
Cell[36702, 826, 265, 4, 33, "Output",ExpressionUUID->"96772995-7f2f-4004-8007-b9946c1988b9"]
Cell[36055, 793, 1317, 39, 53, "Input",ExpressionUUID->"6db522eb-96cb-477e-99c5-d27334b80d0b"],
Cell[37375, 834, 313, 5, 33, "Output",ExpressionUUID->"f0141514-31a9-41e5-9f3f-c4d639d78a0e"]
}, Open ]],
Cell[CellGroupData[{
Cell[37004, 835, 2750, 70, 117, "Input",ExpressionUUID->"16514567-4615-434d-b213-33ee847a7d1f"],
Cell[39757, 907, 1424, 41, 46, "Print",ExpressionUUID->"40fb7223-8e82-42ef-80cb-46c705e0366c"]
Cell[37725, 844, 2750, 70, 117, "Input",ExpressionUUID->"16514567-4615-434d-b213-33ee847a7d1f"],
Cell[40478, 916, 1424, 41, 46, "Print",ExpressionUUID->"40fb7223-8e82-42ef-80cb-46c705e0366c"]
}, Open ]],
Cell[CellGroupData[{
Cell[41218, 953, 645, 17, 29, "Input",ExpressionUUID->"d53a8a6a-e273-4a7d-85e5-98d445720389"],
Cell[41866, 972, 647, 19, 53, "Output",ExpressionUUID->"8b075045-5536-4e81-b7e7-6bb6c8e0532e"]
Cell[41939, 962, 645, 17, 29, "Input",ExpressionUUID->"d53a8a6a-e273-4a7d-85e5-98d445720389"],
Cell[42587, 981, 647, 19, 53, "Output",ExpressionUUID->"8b075045-5536-4e81-b7e7-6bb6c8e0532e"]
}, Open ]],
Cell[CellGroupData[{
Cell[42550, 996, 569, 15, 29, "Input",ExpressionUUID->"0773e3b1-b600-49c8-a187-fa642d2580e8"],
Cell[43122, 1013, 585, 17, 50, "Output",ExpressionUUID->"8b590bfb-da17-4e96-815c-5cfc0a0f84da"]
Cell[43271, 1005, 562, 15, 29, "Input",ExpressionUUID->"0773e3b1-b600-49c8-a187-fa642d2580e8"],
Cell[43836, 1022, 660, 22, 50, "Output",ExpressionUUID->"cc9562cb-8b47-4ffb-a0f0-a0db873c2346"]
}, Open ]],
Cell[43722, 1033, 3118, 96, 214, "Input",ExpressionUUID->"c77c7be5-a519-41e6-b069-691e1006d400"],
Cell[44511, 1047, 6271, 183, 204, "Input",ExpressionUUID->"c77c7be5-a519-41e6-b069-691e1006d400"],
Cell[CellGroupData[{
Cell[46865, 1133, 1567, 30, 52, "Input",ExpressionUUID->"4e5819f2-d11b-4bce-b884-4cd47410fbd4"],
Cell[48435, 1165, 297, 5, 33, "Output",ExpressionUUID->"03549400-50ed-49fa-b63f-eb2e00cbf061"]
Cell[50807, 1234, 1606, 29, 52, "Input",ExpressionUUID->"4e5819f2-d11b-4bce-b884-4cd47410fbd4"],
Cell[52416, 1265, 415, 5, 33, "Output",ExpressionUUID->"b9b74f50-8c3b-4f99-adc7-2201a234e737"]
}, Open ]]
}
]
......
tMatrix3He.nb
View file @ d5a3c0bc
This source diff could not be displayed because it is too large. You can view the blob instead.
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment