Homec4science
Diffusion COSOlver (new_coulomb) Browse

cosolver/GKCoulomb/CurlyNpjmn/src/curlyNpjmn.nb
73330a622fd4new_coulomb

Last Change
Raw File

curlyNpjmn.nb
View Options

(* Content-type: application/vnd.wolfram.mathematica *)
(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)
(* CreatedBy='Mathematica 11.1' *)
(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 158, 7]
NotebookDataLength[ 30933, 896]
NotebookOptionsPosition[ 28901, 840]
NotebookOutlinePosition[ 29239, 855]
CellTagsIndexPosition[ 29196, 852]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell["\<\
Verification of curlNpjmn
\
\>", "Section",
CellChangeTimes->{{3.809949360514091*^9, 3.809949386179058*^9}, {
3.809949421271123*^9,
3.809949442516075*^9}},ExpressionUUID->"ce2d4ef1-cf31-4bde-b7b2-\
8cf83050e13c"],
Cell[CellGroupData[{
Cell["Basis Transformation=", "Chapter",
CellChangeTimes->{{3.8099494626270027`*^9,
3.809949482469351*^9}},ExpressionUUID->"10e49ebd-8512-48fe-8a88-\
b02847721655"],
Cell[BoxData[{
RowBox[{
RowBox[{"\[Tau]1", "[",
RowBox[{"l_", ",", "k_", ",", "m_", ",", "p_", ",", "j_"}], "]"}], ":=",
RowBox[{"Sum", "[",
RowBox[{
RowBox[{
RowBox[{"Binomial", "[",
RowBox[{"l", ",",
RowBox[{
RowBox[{"2", "i"}], "+", "m"}]}], "]"}],
RowBox[{"Binomial", "[",
RowBox[{
RowBox[{"i", "+",
RowBox[{
RowBox[{"(",
RowBox[{"l", "+", "m", "-", "1"}], ")"}], "/", "2"}]}], ",", "l"}],
"]"}],
RowBox[{"Binomial", "[",
RowBox[{
RowBox[{"r", "+",
RowBox[{
RowBox[{"(",
RowBox[{"l", "-", "m"}], ")"}], "/", "2"}], "-", "i"}], ",", "s"}],
"]"}],
RowBox[{"Binomial", "[",
RowBox[{
RowBox[{"r", "+",
RowBox[{
RowBox[{"(",
RowBox[{"l", "-", "m"}], ")"}], "/", "2"}], "-", "i", "-", "s"}],
",", "j"}], "]"}],
RowBox[{"Factorial", "[",
RowBox[{
RowBox[{"2", "i"}], "+", "m"}], "]"}],
RowBox[{"Factorial", "[",
RowBox[{"2",
RowBox[{"(",
RowBox[{"i", "+", "s"}], ")"}]}], "]"}],
RowBox[{
RowBox[{"Factorial", "[",
RowBox[{"r", "+",
RowBox[{
RowBox[{"(",
RowBox[{"l", "-", "m"}], ")"}], "/", "2"}], "-", "i", "-", "s"}],
"]"}], "/",
RowBox[{"(",
RowBox[{
RowBox[{"Factorial", "[",
RowBox[{"2", "i"}], "]"}],
RowBox[{"2", "^",
RowBox[{"(",
RowBox[{"2",
RowBox[{"(",
RowBox[{"i", "+", "s"}], ")"}]}], ")"}]}],
RowBox[{"Factorial", "[",
RowBox[{"i", "+", "s", "-",
RowBox[{"p", "/", "2"}]}], "]"}]}], ")"}]}],
RowBox[{"Boole", "[",
RowBox[{
RowBox[{"2",
RowBox[{"(",
RowBox[{"i", "+", "s"}], ")"}]}], ">=", "p"}], "]"}],
RowBox[{"Boole", "[",
RowBox[{
RowBox[{"r", "+",
RowBox[{
RowBox[{"(",
RowBox[{"l", "-", "m"}], ")"}], "/", "2"}], "-", "i", "-", "s"}], ">=",
"j"}], "]"}]}], ",",
RowBox[{"{",
RowBox[{"i", ",", "0", ",",
RowBox[{
RowBox[{"(",
RowBox[{"l", "-", "m"}], ")"}], "/", "2"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"s", ",", "0", ",",
RowBox[{"r", "+",
RowBox[{
RowBox[{"(",
RowBox[{"l", "-", "m"}], ")"}], "/", "2"}], "-", "i"}]}], "}"}]}],
"]"}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"\[Tau]2", "[",
RowBox[{"l_", ",", "k_", ",", "m_", ",", "p_", ",", "j_"}], "]"}], ":=",
RowBox[{
RowBox[{
RowBox[{"Sqrt", "[", "Pi", "]"}], "/", "2"}], "*",
RowBox[{"Sum", "[",
RowBox[{
RowBox[{
RowBox[{"Binomial", "[",
RowBox[{"l", ",",
RowBox[{
RowBox[{"2", "i"}], "+", "m", "+", "1"}]}], "]"}],
RowBox[{"Binomial", "[",
RowBox[{
RowBox[{"i", "+",
RowBox[{
RowBox[{"(",
RowBox[{"l", "+", "m"}], ")"}], "/", "2"}]}], ",", "l"}], "]"}],
RowBox[{"Binomial", "[",
RowBox[{
RowBox[{"r", "+",
RowBox[{
RowBox[{"(",
RowBox[{"l", "-", "m"}], ")"}], "/", "2"}], "-", "i"}], ",",
"s"}], "]"}],
RowBox[{"Binomial", "[",
RowBox[{
RowBox[{"r", "+",
RowBox[{
RowBox[{"(",
RowBox[{"l", "-", "m"}], ")"}], "/", "2"}], "-", "i", "+",
"\[Eta]", "-", "s"}], ",", "j"}], "]"}],
RowBox[{"Factorial", "[",
RowBox[{
RowBox[{"2", "i"}], "+", "m", "+", "1"}], "]"}],
RowBox[{"Factorial", "[",
RowBox[{"2",
RowBox[{"(",
RowBox[{"i", "+", "s", "+", "1", "-", "\[Eta]"}], ")"}]}], "]"}],
RowBox[{
RowBox[{"Factorial", "[",
RowBox[{"r", "+",
RowBox[{
RowBox[{"(",
RowBox[{"l", "-", "m"}], ")"}], "/", "2"}], "-", "i", "+",
"\[Eta]", "-", "s"}], "]"}], "/",
RowBox[{"(",
RowBox[{
RowBox[{"Factorial", "[",
RowBox[{
RowBox[{"2", "i"}], "+", "1"}], "]"}],
RowBox[{"Factorial", "[",
RowBox[{
RowBox[{"1", "/", "2"}], "-", "\[Eta]"}], "]"}],
RowBox[{"Factorial", "[",
RowBox[{"i", "+", "s", "+", "1", "-", "\[Eta]", "-",
RowBox[{"p", "/", "2"}]}], "]"}],
RowBox[{"2", "^",
RowBox[{"(",
RowBox[{"2",
RowBox[{"(",
RowBox[{"i", "+", "s", "+", "1", "-", "\[Eta]"}], ")"}]}],
")"}]}]}], ")"}]}],
RowBox[{"Boole", "[",
RowBox[{
RowBox[{"r", "+",
RowBox[{
RowBox[{"(",
RowBox[{"l", "-", "m"}], ")"}], "/", "2"}], "-", "i", "+",
"\[Eta]", "-", "s"}], ">=", "j"}], "]"}]}], ",",
RowBox[{"{",
RowBox[{"i", ",", "0", ",",
RowBox[{
RowBox[{
RowBox[{"(",
RowBox[{"l", "-", "m"}], ")"}], "/", "2"}], "-", "1"}]}], "}"}],
",",
RowBox[{"{",
RowBox[{"s", ",", "0", ",",
RowBox[{"r", "+",
RowBox[{
RowBox[{"(",
RowBox[{"l", "-", "m"}], ")"}], "/", "2"}], "-", "i"}]}], "}"}],
",",
RowBox[{"{",
RowBox[{"\[Eta]", ",", "0", ",",
RowBox[{"i", "+", "s", "+", "1", "-",
RowBox[{"p", "/", "2"}]}]}], "}"}]}], "]"}]}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"TP", "[",
RowBox[{"l_", ",", "k_", ",", "m_", ",", "p_", ",", "j_"}], "]"}], ":=",
RowBox[{
RowBox[{"Boole", "[",
RowBox[{"EvenQ", "[", "p", "]"}], "]"}],
RowBox[{
RowBox[{"(",
RowBox[{"-", "1"}], ")"}], "^",
RowBox[{"(",
RowBox[{"m", "+", "j"}], ")"}]}],
RowBox[{
RowBox[{"2", "^", "l"}], "/",
RowBox[{"Factorial", "[", "p", "]"}]}],
RowBox[{"Sum", "[",
RowBox[{
RowBox[{
RowBox[{
RowBox[{
RowBox[{"(",
RowBox[{"-", "1"}], ")"}], "^", "r"}], "/",
RowBox[{"Factorial", "[", "r", "]"}]}], "*",
RowBox[{"Binomial", "[",
RowBox[{
RowBox[{"k", "+", "l", "+",
RowBox[{"1", "/", "2"}]}], ",",
RowBox[{"k", "-", "r"}]}], "]"}],
RowBox[{"(",
RowBox[{
RowBox[{"\[Tau]1", "[",
RowBox[{"l", ",", "k", ",", "m", ",", "p", ",", "j"}], "]"}], "+",
RowBox[{"\[Tau]2", "[",
RowBox[{"l", ",", "k", ",", "m", ",", "p", ",", "j"}], "]"}]}],
")"}]}], ",",
RowBox[{"{",
RowBox[{"r", ",", "0", ",", "k"}], "}"}]}],
"]"}]}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"\[Tau]3", "[",
RowBox[{"l_", ",", "k_", ",", "m_", ",", "p_", ",", "j_"}], "]"}], ":=",
RowBox[{
RowBox[{
RowBox[{"Sqrt", "[", "Pi", "]"}], "/", "2"}], "*",
RowBox[{"Sum", "[",
RowBox[{
RowBox[{
RowBox[{"Binomial", "[",
RowBox[{"l", ",",
RowBox[{
RowBox[{"2", "i"}], "+", "m"}]}], "]"}],
RowBox[{"Binomial", "[",
RowBox[{
RowBox[{"i", "+",
RowBox[{
RowBox[{"(",
RowBox[{"l", "+", "m", "-", "1"}], ")"}], "/", "2"}]}], ",", "l"}],
"]"}],
RowBox[{"Binomial", "[",
RowBox[{
RowBox[{"r", "+",
RowBox[{
RowBox[{"(",
RowBox[{"l", "-", "m", "-", "1"}], ")"}], "/", "2"}], "-", "i"}],
",", "s"}], "]"}],
RowBox[{"Binomial", "[",
RowBox[{
RowBox[{"r", "+",
RowBox[{
RowBox[{"(",
RowBox[{"l", "-", "m", "-", "1"}], ")"}], "/", "2"}], "+",
"\[Eta]", "-", "i", "-", "s"}], ",", "j"}], "]"}],
RowBox[{"Factorial", "[",
RowBox[{
RowBox[{"2", "i"}], "+", "m"}], "]"}],
RowBox[{"Factorial", "[",
RowBox[{"1", "+",
RowBox[{"2",
RowBox[{"(",
RowBox[{"i", "+", "s", "-", "\[Eta]"}], ")"}]}]}], "]"}],
RowBox[{
RowBox[{"Factorial", "[",
RowBox[{"r", "+",
RowBox[{
RowBox[{"(",
RowBox[{"l", "-", "m", "-", "1"}], ")"}], "/", "2"}], "+",
"\[Eta]", "-", "i", "-", "s"}], "]"}], "/",
RowBox[{"(",
RowBox[{
RowBox[{"Factorial", "[",
RowBox[{"2", "i"}], "]"}],
RowBox[{"Factorial", "[",
RowBox[{
RowBox[{"1", "/", "2"}], "-", "\[Eta]"}], "]"}],
RowBox[{"2", "^",
RowBox[{"(",
RowBox[{"1", "+",
RowBox[{"2",
RowBox[{"(",
RowBox[{"i", "+", "s", "-", "\[Eta]"}], ")"}]}]}], ")"}]}],
RowBox[{"Factorial", "[",
RowBox[{"i", "+", "s", "-", "\[Eta]", "-",
RowBox[{
RowBox[{"(",
RowBox[{"p", "-", "1"}], ")"}], "/", "2"}]}], "]"}]}], ")"}]}],
RowBox[{"Boole", "[",
RowBox[{
RowBox[{"r", "+",
RowBox[{
RowBox[{"(",
RowBox[{"l", "-", "m", "-", "1"}], ")"}], "/", "2"}], "+",
"\[Eta]", "-", "i", "-", "s"}], ">=", "j"}], "]"}]}], ",",
RowBox[{"{",
RowBox[{"i", ",", "0", ",",
RowBox[{
RowBox[{"(",
RowBox[{"l", "-", "m", "-", "1"}], ")"}], "/", "2"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"s", ",", "0", ",",
RowBox[{"r", "+",
RowBox[{
RowBox[{"(",
RowBox[{"l", "-", "m", "-", "1"}], ")"}], "/", "2"}], "-", "i"}]}],
"}"}], ",",
RowBox[{"{",
RowBox[{"\[Eta]", ",", "0", ",",
RowBox[{"Floor", "[",
RowBox[{"i", "+", "s", "+",
RowBox[{
RowBox[{"(",
RowBox[{"1", "-", "p"}], ")"}], "/", "2"}]}], "]"}]}], "}"}]}],
"]"}]}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"\[Tau]4", "[",
RowBox[{"l_", ",", "k_", ",", "m_", ",", "p_", ",", "j_"}], "]"}], ":=",
RowBox[{"Sum", "[",
RowBox[{
RowBox[{
RowBox[{"Binomial", "[",
RowBox[{"l", ",",
RowBox[{
RowBox[{"2", "i"}], "+", "m", "+", "1"}]}], "]"}],
RowBox[{"Binomial", "[",
RowBox[{
RowBox[{"i", "+",
RowBox[{
RowBox[{"(",
RowBox[{"l", "+", "m"}], ")"}], "/", "2"}]}], ",", "l"}], "]"}],
RowBox[{"Binomial", "[",
RowBox[{
RowBox[{"r", "+",
RowBox[{
RowBox[{"(",
RowBox[{"l", "-", "m", "-", "1"}], ")"}], "/", "2"}], "-", "i"}],
",", "s"}], "]"}],
RowBox[{"Binomial", "[",
RowBox[{
RowBox[{"r", "+",
RowBox[{
RowBox[{"(",
RowBox[{"l", "-", "m", "-", "1"}], ")"}], "/", "2"}], "-", "i", "-",
"s"}], ",", "j"}], "]"}],
RowBox[{"Factorial", "[",
RowBox[{
RowBox[{"2", "i"}], "+", "m", "+", "1"}], "]"}],
RowBox[{"Factorial", "[",
RowBox[{"1", "+",
RowBox[{"2",
RowBox[{"(",
RowBox[{"i", "+", "s"}], ")"}]}]}], "]"}],
RowBox[{
RowBox[{"Factorial", "[",
RowBox[{"r", "+",
RowBox[{
RowBox[{"(",
RowBox[{"l", "-", "m", "-", "1"}], ")"}], "/", "2"}], "-", "i", "-",
"s"}], "]"}], "/",
RowBox[{"(",
RowBox[{
RowBox[{"Factorial", "[",
RowBox[{
RowBox[{"2", "i"}], "+", "1"}], "]"}],
RowBox[{"2", "^",
RowBox[{"(",
RowBox[{"1", "+",
RowBox[{"2",
RowBox[{"(",
RowBox[{"i", "+", "s"}], ")"}]}]}], ")"}]}],
RowBox[{"Factorial", "[",
RowBox[{"i", "+", "s", "-",
RowBox[{
RowBox[{"(",
RowBox[{"p", "-", "1"}], ")"}], "/", "2"}]}], "]"}]}], ")"}]}],
RowBox[{"Boole", "[",
RowBox[{
RowBox[{"1", "+",
RowBox[{"2",
RowBox[{"(",
RowBox[{"i", "+", "s"}], ")"}]}]}], ">=", "p"}], "]"}],
RowBox[{"Boole", "[",
RowBox[{
RowBox[{"r", "+",
RowBox[{
RowBox[{"(",
RowBox[{"l", "-", "m", "-", "1"}], ")"}], "/", "2"}], "-", "i", "-",
"s"}], ">=", "j"}], "]"}]}], ",",
RowBox[{"{",
RowBox[{"i", ",", "0", ",",
RowBox[{
RowBox[{"(",
RowBox[{"l", "-", "m", "-", "1"}], ")"}], "/", "2"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"s", ",", "0", ",",
RowBox[{"r", "+",
RowBox[{
RowBox[{"(",
RowBox[{"l", "-", "m", "-", "1"}], ")"}], "/", "2"}], "-", "i"}]}],
"}"}]}], "]"}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"TI", "[",
RowBox[{"l_", ",", "k_", ",", "m_", ",", "p_", ",", "j_"}], "]"}], ":=",
RowBox[{
RowBox[{"Boole", "[",
RowBox[{"OddQ", "[", "p", "]"}], "]"}],
RowBox[{
RowBox[{"(",
RowBox[{"-", "1"}], ")"}], "^",
RowBox[{"(",
RowBox[{"m", "+", "j"}], ")"}]}],
RowBox[{
RowBox[{"2", "^", "l"}], "/",
RowBox[{"Factorial", "[", "p", "]"}]}],
RowBox[{"Sum", "[",
RowBox[{
RowBox[{
RowBox[{
RowBox[{
RowBox[{"(",
RowBox[{"-", "1"}], ")"}], "^", "r"}], "/",
RowBox[{"Factorial", "[", "r", "]"}]}], "*",
RowBox[{"Binomial", "[",
RowBox[{
RowBox[{"l", "+", "k", "+",
RowBox[{"1", "/", "2"}]}], ",",
RowBox[{"k", "-", "r"}]}], "]"}],
RowBox[{"(",
RowBox[{
RowBox[{"\[Tau]3", "[",
RowBox[{"l", ",", "k", ",", "m", ",", "p", ",", "j"}], "]"}], "+",
RowBox[{"\[Tau]4", "[",
RowBox[{"l", ",", "k", ",", "m", ",", "p", ",", "j"}], "]"}]}],
")"}]}], ",",
RowBox[{"{",
RowBox[{"r", ",", "0", ",", "k"}], "}"}]}],
"]"}]}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"T", "[",
RowBox[{"l_", ",", "k_", ",", "m_", ",", "p_", ",", "j_"}], "]"}], ":=",
RowBox[{
RowBox[{
RowBox[{"Boole", "[",
RowBox[{"EvenQ", "[",
RowBox[{"l", "-", "m"}], "]"}], "]"}], "*",
RowBox[{"TP", "[",
RowBox[{"l", ",", "k", ",", "m", ",", "p", ",", "j"}], "]"}]}], "+",
RowBox[{
RowBox[{"Boole", "[",
RowBox[{"OddQ", "[",
RowBox[{"l", "-", "m"}], "]"}], "]"}], "*",
RowBox[{"TI", "[",
RowBox[{"l", ",", "k", ",", "m", ",", "p", ",", "j"}],
"]"}]}]}]}]}], "Input",
CellChangeTimes->{{3.809949457224197*^9,
3.809949484917872*^9}},ExpressionUUID->"79f0a12d-39c6-4f4d-aa8f-\
ae303b0591b4"],
Cell[BoxData[""], "Input",ExpressionUUID->"55d663e2-929b-44de-ba52-35dd6cf816b0"]
}, Open ]],
Cell[CellGroupData[{
Cell["Laguerre-Laguerre-Associated Laguerre-x^m 2 Laguerre", "Chapter",
CellChangeTimes->{{3.809949491503059*^9, 3.809949501905436*^9}, {
3.809949541643442*^9,
3.809949566633795*^9}},ExpressionUUID->"2eaf49d1-81ca-441f-9d45-\
dc66b7f48283"],
Cell[BoxData[{
RowBox[{
SubscriptBox["L",
RowBox[{"p_", ",", "j_", ",", "l_", " "}]], ":=",
RowBox[{
SuperscriptBox[
RowBox[{"(",
RowBox[{"-", "1"}], ")"}], "l"], "*",
FractionBox[
RowBox[{"Factorial", "[",
RowBox[{"p", " ", "+", " ", "j", " ", "+",
RowBox[{"1", "/", "2"}]}], "]"}],
RowBox[{
RowBox[{"Factorial", "[",
RowBox[{"j", "-", "l"}], "]"}], "*",
RowBox[{"Factorial", "[",
RowBox[{"p", "+", "l", " ", "+",
RowBox[{"1", "/", "2"}]}], "]"}], "*",
RowBox[{"Factorial", "[", "l", "]"}]}]]}]}], "\[IndentingNewLine]",
RowBox[{
SubscriptBox["d",
RowBox[{"m_", ",", "n_", ",", "k_", ",", "s_"}]], ":=",
RowBox[{
SubsuperscriptBox["\[Sum]",
RowBox[{"np", "=", "0"}], "n"],
RowBox[{
SubsuperscriptBox["\[Sum]",
RowBox[{"kp", "=", "0"}], "k"],
RowBox[{
SubsuperscriptBox["\[Sum]",
RowBox[{"sp", "=", "0"}], "s"],
RowBox[{
SubscriptBox["L",
RowBox[{
RowBox[{
RowBox[{"-", "1"}], "/", "2"}], ",", "k", ",", "kp"}]], "*",
SubscriptBox["L",
RowBox[{
RowBox[{"m", "-",
RowBox[{"1", "/", "2"}]}], ",", "n", ",", "np"}]], "*",
SubscriptBox["L",
RowBox[{
RowBox[{
RowBox[{"-", "1"}], "/", "2"}], ",", "s", ",", "sp"}]], " ", "*",
RowBox[{"Factorial", "[",
RowBox[{
"np", " ", "+", " ", "kp", " ", "+", " ", "sp", " ", "+", " ", "m"}],
"]"}]}]}]}]}]}]}], "Input",
CellChangeTimes->{{3.809950373710847*^9,
3.809950401007883*^9}},ExpressionUUID->"e04009ca-0afa-4338-b4f6-\
5056f8bd7989"]
}, Open ]],
Cell[CellGroupData[{
Cell["Curly Npjmn", "Chapter",
CellChangeTimes->{{3.809950236940247*^9, 3.809950280487936*^9}, {
3.809950333947772*^9, 3.8099503476636753`*^9}, {3.809950410316083*^9,
3.8099504138476562`*^9}, {3.8099507399873743`*^9,
3.80995075803596*^9}},ExpressionUUID->"dfddeaf4-8bb9-43d7-99a8-\
724e6a38d875"],
Cell[BoxData[
RowBox[{"\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"Pmax", " ", "=", "15"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"Jmax", " ", "=", " ", "6"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"curlyN", "[",
RowBox[{"p_", ",", "j_", ",", "m_", ",", "n_"}], "]"}], " ", ":=",
RowBox[{"Sum", "[",
RowBox[{
RowBox[{"Sum", "[",
RowBox[{
RowBox[{"Sum", "[",
RowBox[{
RowBox[{
SubscriptBox["d",
RowBox[{"m", ",", "n", ",", "h", ",", "s1"}]],
RowBox[{"T", "[",
RowBox[{"p", ",", "j", ",", "m", ",", "g", ",", "h"}], "]"}],
SqrtBox[
RowBox[{
SuperscriptBox["2", "g"],
RowBox[{"g", "!"}], " "}]],
SubscriptBox["N",
RowBox[{"g", ",", "s1"}]]}], ",",
RowBox[{"{",
RowBox[{"s1", ",", "0", ",",
RowBox[{"Min", "[",
RowBox[{"Jmax", ",",
RowBox[{"n", " ", "+", "m", "+", "h"}]}], "]"}]}], "}"}]}],
"]"}], ",",
RowBox[{"{",
RowBox[{"h", ",", "0", ",",
RowBox[{"j", " ", "+", " ",
RowBox[{"Floor", "[",
RowBox[{
RowBox[{"(",
RowBox[{"m", "+", "p"}], ")"}], "/", "2"}], "]"}]}]}], "}"}]}],
"]"}], ",",
RowBox[{"{",
RowBox[{"g", ",", "0", ",",
RowBox[{"Min", "[",
RowBox[{"Pmax", ",",
RowBox[{"p", " ", "+",
RowBox[{"2", "j"}], "+", "m"}]}], "]"}]}], "}"}]}],
"]"}]}]}]}]], "Input",
CellChangeTimes->{{3.809950759463616*^9, 3.809951059819352*^9}, {
3.8099511283234663`*^9, 3.8099511441926737`*^9}, {3.8099513008409557`*^9,
3.809951394384892*^9}, {3.809952885451585*^9, 3.8099528868892517`*^9}, {
3.809952983594968*^9, 3.809952984672966*^9}, {3.809953511103022*^9,
3.809953570764791*^9}, {3.809953661590807*^9, 3.809953687740341*^9}, {
3.8099537321265078`*^9, 3.8099537375399113`*^9}, {3.809954027964505*^9,
3.809954035505652*^9}, {3.8100155750949697`*^9, 3.810015592707633*^9}, {
3.810031920315482*^9,
3.810031922226462*^9}},ExpressionUUID->"61247607-09e1-4380-ac7e-\
a5907aab6555"],
Cell[BoxData[""], "Input",
CellChangeTimes->{{3.809956295780773*^9,
3.809956296657847*^9}},ExpressionUUID->"eadb081c-3ade-4bec-8ace-\
6405d9823a8e"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"N", "[",
RowBox[{
RowBox[{"curlyN", "[",
RowBox[{"11", ",", "2", ",", "2", ",", "3"}], "]"}], ",", "15"}],
"]"}]], "Input",
CellChangeTimes->CompressedData["
1:eJwdz1toEgAYBWBbusW6kjVzjELXWgVRxNCKipwPLRIZs9pcpWUWC9lksswS
x1i1mcvWXI4YlkWj1E26rMa6IiQRXQhqoxCxbSxSM6dtD24a5Pkffj7Oy+H8
PJW26mQOg8Eoyx78xL438ksVE+3+pxiFco03AiXHcqMwHdDHYVBoT8KdsXgK
usSGWRi+pJgfzprpSefDq9XlHPiW5+RCfr9bClmDCTkc2JY6CpcMzJEhg0EN
G3eESdmd0Ub4+6BXB+P8KybY4DY2w4Qo2ErZFiL7Jg0O+EVmJN94mC7YG7ve
D2Ndhx7QLnsNaVBOPYe219U++CzD8cPxNJes9A99gM2hkY/QeqMkCO2TtaRK
NzQLC+s+kw0xDSuCfxP15LSXw4bLOWJyap+DA1dxneScOVkEGZkZkqfv5EFJ
wkbKF4Q3wr+CKOnIDW+Fj7XdInjq/KO9sHhFfgXc0FMhgbr6/WSZx6aHzDW3
STd7sxG+bGu+AGu6LWaYvHnAAoO3Wq/B40pzL3SZTO9gQ2UL6R9kfYdbPHlk
7bKVEzDQwicl/c7t0ayCnylS6P4hgYsKJsgOZUQOx45wa2F75Ss13MP0keKc
ojPQdPa9AUp10cvQLBBYYV/BvC7q0zyxQUvON27eiZhoplxcCOvOja+DTXcD
pdDRPnx4YdbiiFcBH+6SqKH/j5Ts7ChrgvflQpL3ws1dnNWqHSCHTydXw4vs
aXLM93Qt1LfJSqBlk7J0adaqr8Xr4X9sxcJK
"],ExpressionUUID->"dc7bd711-4f20-4581-8768-aa7b440a2510"],
Cell[BoxData[
RowBox[{
RowBox[{"1.00332618619426304695403508`15.*^8", " ",
SubscriptBox["N",
RowBox[{"1", ",", "3"}]]}], "-",
RowBox[{"8.026609489554104375632280641`15.*^8", " ",
SubscriptBox["N",
RowBox[{"1", ",", "4"}]]}], "+",
RowBox[{"3.1855606411667851740790613795`15.*^9", " ",
SubscriptBox["N",
RowBox[{"1", ",", "5"}]]}], "-",
RowBox[{"7.8761105616249649185891753793`15.*^9", " ",
SubscriptBox["N",
RowBox[{"1", ",", "6"}]]}], "+",
RowBox[{"1.075216275765017945623194526`15.*^8", " ",
SubscriptBox["N",
RowBox[{"3", ",", "2"}]]}], "-",
RowBox[{"8.601730206120143564985556205`15.*^8", " ",
SubscriptBox["N",
RowBox[{"3", ",", "3"}]]}], "+",
RowBox[{"3.5482137100245592205565419345`15.*^9", " ",
SubscriptBox["N",
RowBox[{"3", ",", "4"}]]}], "-",
RowBox[{"9.2468599715791543323594729203`15.*^9", " ",
SubscriptBox["N",
RowBox[{"3", ",", "5"}]]}], "+",
RowBox[{"1.61282441364752691843479178843`15.*^10", " ",
SubscriptBox["N",
RowBox[{"3", ",", "6"}]]}], "+",
RowBox[{"1.66825973931104395748210635`15.*^7", " ",
SubscriptBox["N",
RowBox[{"5", ",", "1"}]]}], "-",
RowBox[{"1.334607791448835165985685076`15.*^8", " ",
SubscriptBox["N",
RowBox[{"5", ",", "2"}]]}], "+",
RowBox[{"5.922322074554206049061477525`15.*^8", " ",
SubscriptBox["N",
RowBox[{"5", ",", "3"}]]}], "-",
RowBox[{"1.6849423367041543970569274085`15.*^9", " ",
SubscriptBox["N",
RowBox[{"5", ",", "4"}]]}], "+",
RowBox[{"3.2113999981737596181530547141`15.*^9", " ",
SubscriptBox["N",
RowBox[{"5", ",", "5"}]]}], "-",
RowBox[{"4.1039189587051681354059816086`15.*^9", " ",
SubscriptBox["N",
RowBox[{"5", ",", "6"}]]}], "-",
RowBox[{"5.7540508956106170116372551`15.*^6", " ",
SubscriptBox["N",
RowBox[{"7", ",", "0"}]]}], "+",
RowBox[{"4.60324071648849360930980412`15.*^7", " ",
SubscriptBox["N",
RowBox[{"7", ",", "1"}]]}], "-",
RowBox[{"2.474241885112565315004019716`15.*^8", " ",
SubscriptBox["N",
RowBox[{"7", ",", "2"}]]}], "+",
RowBox[{"8.400914307591500836990392524`15.*^8", " ",
SubscriptBox["N",
RowBox[{"7", ",", "3"}]]}], "-",
RowBox[{"1.8700665410734505287821079249`15.*^9", " ",
SubscriptBox["N",
RowBox[{"7", ",", "4"}]]}], "+",
RowBox[{"2.7389282263106536975393334531`15.*^9", " ",
SubscriptBox["N",
RowBox[{"7", ",", "5"}]]}], "-",
RowBox[{"2.5490445467555033361553040331`15.*^9", " ",
SubscriptBox["N",
RowBox[{"7", ",", "6"}]]}], "-",
RowBox[{"1.07071800199441957220402131`15.*^7", " ",
SubscriptBox["N",
RowBox[{"9", ",", "1"}]]}], "+",
RowBox[{"6.42430801196651743322412784`15.*^7", " ",
SubscriptBox["N",
RowBox[{"9", ",", "2"}]]}], "-",
RowBox[{"2.034364203789397187187640482`15.*^8", " ",
SubscriptBox["N",
RowBox[{"9", ",", "3"}]]}], "+",
RowBox[{"3.854584807179910459934476703`15.*^8", " ",
SubscriptBox["N",
RowBox[{"9", ",", "4"}]]}], "-",
RowBox[{"4.389943808177120246036487356`15.*^8", " ",
SubscriptBox["N",
RowBox[{"9", ",", "5"}]]}], "+",
RowBox[{"2.783866805185490887730455396`15.*^8", " ",
SubscriptBox["N",
RowBox[{"9", ",", "6"}]]}], "+",
RowBox[{"3.59353124603780229882528354`15.*^7", " ",
SubscriptBox["N",
RowBox[{"11", ",", "2"}]]}], "-",
RowBox[{"1.437412498415120919530113416`15.*^8", " ",
SubscriptBox["N",
RowBox[{"11", ",", "3"}]]}], "+",
RowBox[{"2.515471872226461609177698478`15.*^8", " ",
SubscriptBox["N",
RowBox[{"11", ",", "4"}]]}], "-",
RowBox[{"2.156118747622681379295170124`15.*^8", " ",
SubscriptBox["N",
RowBox[{"11", ",", "5"}]]}], "+",
RowBox[{"7.18706249207560459765056708`15.*^7", " ",
SubscriptBox["N",
RowBox[{"11", ",", "6"}]]}], "+",
RowBox[{"1.87013295322401476579665868`15.*^7", " ",
SubscriptBox["N",
RowBox[{"13", ",", "3"}]]}], "-",
RowBox[{"3.74026590644802953159331736`15.*^7", " ",
SubscriptBox["N",
RowBox[{"13", ",", "4"}]]}], "+",
RowBox[{"1.87013295322401476579665868`15.*^7", " ",
SubscriptBox["N",
RowBox[{"13", ",", "5"}]]}]}]], "Output",
CellChangeTimes->{{3.810024127419029*^9, 3.810024189941394*^9},
3.810024234995727*^9, 3.810024362444035*^9, 3.810027734021021*^9,
3.8100277764440823`*^9, {3.8100278227676907`*^9, 3.810027831906974*^9}, {
3.8100319403672733`*^9,
3.810031955672827*^9}},ExpressionUUID->"3ad22726-c0b6-426b-9615-\
6ae3f7b083cf"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{"Table", "[",
RowBox[{
RowBox[{"N", "[",
RowBox[{"T", "[",
RowBox[{"14", ",", "4", ",", "4", ",", "2", ",", "i"}], "]"}], "]"}],
",",
RowBox[{"{",
RowBox[{"i", ",", "0", ",", "9"}], "}"}]}],
"]"}], "\[IndentingNewLine]",
RowBox[{"N", "[",
RowBox[{"T", "[",
RowBox[{"14", ",", "4", ",", "4", ",", "2", ",", "3"}], "]"}],
"]"}]}], "Input",
CellChangeTimes->CompressedData["
1:eJxTTMoPSmViYGAQA2IQXXZoiyd78htHLadV3iB6Xd+PaBCtwSkZA6JDTsim
gOinZvJg+k/f2mwQ7f/CIw9Eq+WfbAbR91Z8aQPRQTnTekC03t4ZYJonPHIC
iLbujgbTuyYoTQPRPY0qYNpPTHIWiG4+8GsRiHbQ5VgBor+W8K0D0awXkzeB
6MvfZ+0D0UmVcqdBtN2NQ2D6Vv7y8yD6S8dKMO0So3cVRCvcMQfTDG9/tHEB
admIwG4Q/XXvh34Q/an60SQQreLGdBhEf9nEAqYfbVx2HESzlK0E057zJK+C
6AtOimB6QuTsmyBa8P8mMH3mjqkCN5AuXjBXBUSb/eSbDaJX7Zo/D0QDAIyJ
o5Q=
"],ExpressionUUID->"b9257fc9-2adb-4647-bb19-7b538a7c70a9"],
Cell[BoxData[
RowBox[{"{",
RowBox[{
"5.348210445703125`*^10", ",", "3.889607596875`*^10", ",",
"2.7227253178125`*^10", ",", "1.815150211875`*^10", ",",
"1.134468882421875`*^10", ",", "6.482679328125`*^9", ",",
"3.2413396640625`*^9", ",", "1.296535865625`*^9", ",",
"3.2413396640625`*^8", ",", "0.`"}], "}"}]], "Output",
CellChangeTimes->{{3.8100166053061523`*^9, 3.810016624976615*^9}, {
3.810016662704342*^9, 3.810016757824242*^9}, {3.810016858644104*^9,
3.810016862753151*^9}, 3.810016894423894*^9, 3.810016938436727*^9, {
3.810022456741439*^9, 3.810022491290533*^9}, {3.810022528253839*^9,
3.8100225517952967`*^9}, 3.81002293697961*^9, 3.8100229740478773`*^9,
3.810023081702433*^9, 3.810023118354491*^9, 3.810023644598068*^9, {
3.8100236824324923`*^9, 3.810023718900343*^9}, {3.8100246476214867`*^9,
3.810024693695239*^9}},ExpressionUUID->"e96cded0-dc4c-4635-a97a-\
6adb2ef68b93"],
Cell[BoxData["1.815150211875`*^10"], "Output",
CellChangeTimes->{{3.8100166053061523`*^9, 3.810016624976615*^9}, {
3.810016662704342*^9, 3.810016757824242*^9}, {3.810016858644104*^9,
3.810016862753151*^9}, 3.810016894423894*^9, 3.810016938436727*^9, {
3.810022456741439*^9, 3.810022491290533*^9}, {3.810022528253839*^9,
3.8100225517952967`*^9}, 3.81002293697961*^9, 3.8100229740478773`*^9,
3.810023081702433*^9, 3.810023118354491*^9, 3.810023644598068*^9, {
3.8100236824324923`*^9, 3.810023718900343*^9}, {3.8100246476214867`*^9,
3.810024693795895*^9}},ExpressionUUID->"12b75c2f-b224-40ed-837a-\
393fc1e7e0f3"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"NumberForm", "[",
RowBox[{"1.815150211875`*^10", ",", "16"}], "]"}]], "Input",
NumberMarks->False,ExpressionUUID->"3987dbd5-d974-41d0-a8b8-1b067dd2eba9"],
Cell[BoxData[
TagBox[
InterpretationBox[
RowBox[{"\<\"1.815150211875\"\>", "\[Times]",
SuperscriptBox["10", "\<\"10\"\>"]}],
1.815150211875*^10,
AutoDelete->True],
NumberForm[#, 16]& ]], "Output",
CellChangeTimes->{
3.81002372572234*^9},ExpressionUUID->"6f671564-40bc-4003-bc36-e23292f4b5b3"]
}, Open ]]
}, Open ]]
},
WindowSize->{1272, 687},
WindowMargins->{{251, Automatic}, {Automatic, 114}},
FrontEndVersion->"11.1 for Linux x86 (64-bit) (April 18, 2017)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)
(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 228, 7, 101, "Section", "ExpressionUUID" -> \
"ce2d4ef1-cf31-4bde-b7b2-8cf83050e13c"],
Cell[CellGroupData[{
Cell[811, 31, 168, 3, 55, "Chapter", "ExpressionUUID" -> \
"10e49ebd-8512-48fe-8a88-b02847721655"],
Cell[982, 36, 14501, 446, 580, "Input", "ExpressionUUID" -> \
"79f0a12d-39c6-4f4d-aa8f-ae303b0591b4"],
Cell[15486, 484, 81, 0, 33, "Input", "ExpressionUUID" -> \
"55d663e2-929b-44de-ba52-35dd6cf816b0"]
}, Open ]],
Cell[CellGroupData[{
Cell[15604, 489, 246, 4, 66, "Chapter", "ExpressionUUID" -> \
"2eaf49d1-81ca-441f-9d45-dc66b7f48283"],
Cell[15853, 495, 1649, 50, 88, "Input", "ExpressionUUID" -> \
"e04009ca-0afa-4338-b4f6-5056f8bd7989"]
}, Open ]],
Cell[CellGroupData[{
Cell[17539, 550, 305, 5, 66, "Chapter", "ExpressionUUID" -> \
"dfddeaf4-8bb9-43d7-99a8-724e6a38d875"],
Cell[17847, 557, 2252, 57, 174, "Input", "ExpressionUUID" -> \
"61247607-09e1-4380-ac7e-a5907aab6555"],
Cell[20102, 616, 152, 3, 33, "Input", "ExpressionUUID" -> \
"eadb081c-3ade-4bec-8ace-6405d9823a8e"],
Cell[CellGroupData[{
Cell[20279, 623, 1031, 21, 35, "Input", "ExpressionUUID" -> \
"dc7bd711-4f20-4581-8768-aa7b440a2510"],
Cell[21313, 646, 4540, 115, 285, "Output", "ExpressionUUID" -> \
"3ad22726-c0b6-426b-9615-6ae3f7b083cf"]
}, Open ]],
Cell[CellGroupData[{
Cell[25890, 766, 875, 22, 59, "Input", "ExpressionUUID" -> \
"b9257fc9-2adb-4647-bb19-7b538a7c70a9"],
Cell[26768, 790, 931, 16, 37, "Output", "ExpressionUUID" -> \
"e96cded0-dc4c-4635-a97a-6adb2ef68b93"],
Cell[27702, 808, 636, 9, 33, "Output", "ExpressionUUID" -> \
"12b75c2f-b224-40ed-837a-393fc1e7e0f3"]
}, Open ]],
Cell[CellGroupData[{
Cell[28375, 822, 181, 3, 37, "Input", "ExpressionUUID" -> \
"3987dbd5-d974-41d0-a8b8-1b067dd2eba9"],
Cell[28559, 827, 314, 9, 47, "Output", "ExpressionUUID" -> \
"6f671564-40bc-4003-bc36-e23292f4b5b3"]
}, Open ]]
}, Open ]]
}
]
*)
c4science · Help