(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 12.1' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 158, 7] NotebookDataLength[ 133860, 4095] NotebookOptionsPosition[ 119687, 3827] NotebookOutlinePosition[ 120200, 3845] CellTagsIndexPosition[ 120157, 3842] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ RowBox[{"Exit", "[", "]"}]], "Input", CellChangeTimes->{{3.827398416270541*^9, 3.827398417518792*^9}}, CellLabel->"In[19]:=",ExpressionUUID->"851884f2-edd5-4c2d-9f36-5bd9b7a85eec"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"<<", "Sigma.m"}]], "Input", CellChangeTimes->{{3.828588960235929*^9, 3.828588961922735*^9}, { 3.8316223191519547`*^9, 3.831622319885988*^9}, {3.831695355435175*^9, 3.8316953564021*^9}}, CellLabel->"In[1]:=",ExpressionUUID->"d1b6b1e4-eb24-487b-9e19-61e898f3794f"], Cell[BoxData[ RowBox[{ RowBox[{ StyleBox["Sigma", FontColor->RGBColor[1, 0, 0]], " ", RowBox[{ RowBox[{"-", " ", "A"}], " ", "summation", " ", "package", " ", "by", " ", "Carsten", " ", "Schneider"}]}], " ", "\[LongDash]", " ", RowBox[{"\[Copyright]", " ", "RISC"}], " ", "\[LongDash]", " ", RowBox[{"V", " ", "2.85", " ", RowBox[{"(", RowBox[{ RowBox[{"June", " ", "2"}], ",", " ", "2021"}], ")"}], " ", ButtonBox[ StyleBox["Help", "Hyperlink", FontVariations->{"Underline"->True}], Appearance->{Automatic, None, "Normal", Automatic}, ButtonData:>{"", ""}, ButtonFunction:>Sigma`Version`Private`SigmaHelp[], Evaluator->Automatic]}]}]], "Print", CellFrame->0.5, CellChangeTimes->{3.831695359479046*^9}, FontColor->RGBColor[0, 0, 0], Background->RGBColor[0.796887, 0.789075, 0.871107], ButtonBoxOptions->{ Active->True},ExpressionUUID->"8bc9d068-7ad0-4b56-b962-d121c81db79f"] }, Open ]], Cell[CellGroupData[{ Cell["Section 1: Symbolic summation \[LineSeparator](a short introduction)", \ "Title", CellChangeTimes->{{3.644643480667371*^9, 3.6446434845602093`*^9}, { 3.827385629144805*^9, 3.827385646733103*^9}, {3.924156489303849*^9, 3.924156492830974*^9}},ExpressionUUID->"704c8277-85dd-4b00-be79-\ 882101fd2c60"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"<<", "\"\\""}], ";"}]], "Input", CellChangeTimes->{{3.644652088536532*^9, 3.6446520908145857`*^9}, { 3.92294852627707*^9, 3.922948529817071*^9}, {3.9241549217239037`*^9, 3.924154921910625*^9}}, CellLabel->"In[2]:=",ExpressionUUID->"0c38568f-6427-47f1-a08e-fcd0599ea15e"], Cell[BoxData[ RowBox[{ RowBox[{ StyleBox["Sigma", FontColor->RGBColor[1, 0, 0]], " ", RowBox[{ RowBox[{"-", " ", "A"}], " ", "summation", " ", "package", " ", "by", " ", "Carsten", " ", "Schneider"}]}], " ", "\[LongDash]", " ", RowBox[{"\[Copyright]", " ", "RISC"}], " ", "\[LongDash]", " ", RowBox[{"V", " ", "2.895", " ", RowBox[{"(", RowBox[{ RowBox[{"September", " ", "20"}], ",", " ", "2022"}], ")"}], " ", "click", " ", "for", " ", ButtonBox[ StyleBox["Help", "Hyperlink", FontVariations->{"Underline"->True}], Appearance->{Automatic, None, "Normal", Automatic}, ButtonData:>{"", ""}, ButtonFunction:>Sigma`Version`Private`SigmaHelp[], Evaluator->Automatic]}]}]], "Print", CellFrame->0.5, CellChangeTimes->{3.92415492301161*^9}, FontColor->RGBColor[0, 0, 0], Background->RGBColor[0.796887, 0.789075, 0.871107], ButtonBoxOptions->{ Active->True},ExpressionUUID->"be8a0400-d84d-472f-aa39-76a649e61842"] }, Open ]], Cell[CellGroupData[{ Cell["Rational summation", "Subsection", CellChangeTimes->{{3.827386021079687*^9, 3.827386042019915*^9}},ExpressionUUID->"c0c9b550-ef2e-4668-b871-\ a8e57b157f8e"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"mySum", "=", RowBox[{"SigmaSum", "[", RowBox[{ RowBox[{"1", "/", "k"}], ",", RowBox[{"{", RowBox[{"k", ",", "1", ",", "n"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.707537463560733*^9, 3.707537495998261*^9}, { 3.827385699111377*^9, 3.827385702835271*^9}}, CellLabel->"In[3]:=",ExpressionUUID->"236e52e3-cd97-46b7-9ebd-f7fbe45f093d"], Cell[BoxData[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"k", "=", "1"}], "n"], FractionBox["1", "k"]}]], "Output", CellChangeTimes->{3.878814086856586*^9, 3.922948439896731*^9, 3.922948533724226*^9, 3.924154924930606*^9}, CellLabel->"Out[3]=",ExpressionUUID->"5d211ac4-6321-4a72-8730-84cebd5aafa0"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"SigmaReduce", "[", "mySum", "]"}]], "Input", CellChangeTimes->{{3.7075374978604116`*^9, 3.7075375024416723`*^9}}, CellLabel->"In[4]:=",ExpressionUUID->"ba8cecd1-67dd-456c-a79e-da3b767c5f7f"], Cell[BoxData[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubscriptBox["\[Iota]", "1"]]}]], "Output", CellChangeTimes->{3.878814113795476*^9, 3.922948441272003*^9, 3.9229485339312153`*^9, 3.924154926042242*^9}, CellLabel->"Out[4]=",ExpressionUUID->"81cd4355-5151-4689-9e20-f0fd57ec22a6"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"mySum", "=", RowBox[{"SigmaSum", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"x", "^", "2"}], "+", "x", "+", "1"}], ")"}], "/", RowBox[{ RowBox[{"(", RowBox[{"x", "+", "1"}], ")"}], "^", "2"}]}], "/", RowBox[{ RowBox[{"(", RowBox[{"x", "+", "2"}], ")"}], "^", "3"}]}], ",", RowBox[{"{", RowBox[{"x", ",", "1", ",", "n"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.8273857129117193`*^9, 3.827385744691577*^9}, { 3.82738610259368*^9, 3.8273861109762077`*^9}, 3.8273861549334307`*^9, { 3.878814273284793*^9, 3.87881428756211*^9}}, CellLabel->"In[5]:=",ExpressionUUID->"6353f64f-9f87-4d14-95a7-85380780b73e"], Cell[BoxData[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"x", "=", "1"}], "n"], FractionBox[ RowBox[{"1", "+", "x", "+", SuperscriptBox["x", "2"]}], RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"1", "+", "x"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"(", RowBox[{"2", "+", "x"}], ")"}], "3"]}]]}]], "Output", CellChangeTimes->{ 3.878814240020565*^9, {3.878814274325382*^9, 3.8788142881024523`*^9}, 3.922948442932727*^9, 3.922948534305649*^9, 3.924154926939382*^9}, CellLabel->"Out[5]=",ExpressionUUID->"7d693058-441d-4945-a55c-8d76ba893142"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"SigmaReduce", "[", "mySum", "]"}]], "Input", CellChangeTimes->{{3.827385747554729*^9, 3.8273857533630342`*^9}}, CellLabel->"In[6]:=",ExpressionUUID->"a519d9e9-d59b-4115-b8df-760330976bd9"], Cell[BoxData[ RowBox[{ RowBox[{"-", FractionBox[ RowBox[{"n", " ", RowBox[{"(", RowBox[{"1236", "+", RowBox[{"3538", " ", "n"}], "+", RowBox[{"4103", " ", SuperscriptBox["n", "2"]}], "+", RowBox[{"2393", " ", SuperscriptBox["n", "3"]}], "+", RowBox[{"697", " ", SuperscriptBox["n", "4"]}], "+", RowBox[{"81", " ", SuperscriptBox["n", "5"]}]}], ")"}]}], RowBox[{"8", " ", SuperscriptBox[ RowBox[{"(", RowBox[{"1", "+", "n"}], ")"}], "3"], " ", SuperscriptBox[ RowBox[{"(", RowBox[{"2", "+", "n"}], ")"}], "3"]}]]}], "+", RowBox[{"3", " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubsuperscriptBox["\[Iota]", "1", "3"]]}], ")"}]}], "+", RowBox[{"4", " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubsuperscriptBox["\[Iota]", "1", "2"]]}], ")"}]}]}]], "Output", CellChangeTimes->{ 3.878814246115511*^9, {3.878814276121427*^9, 3.878814295734398*^9}, 3.922948445207117*^9, 3.922948535651105*^9, 3.92415492821024*^9}, CellLabel->"Out[6]=",ExpressionUUID->"4f62fa47-a6e0-4fb1-91b0-671be22a7419"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Generalized summation", "Subsection", CellChangeTimes->{{3.8273860594222527`*^9, 3.827386062216543*^9}},ExpressionUUID->"5476d50b-8a19-49b0-b076-\ b62a8bb09a08"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"mySum", "=", RowBox[{ RowBox[{"SigmaSum", "[", RowBox[{ RowBox[{"SigmaHNumber", "[", "k", "]"}], ",", RowBox[{"{", RowBox[{"k", ",", "1", ",", "n"}], "}"}]}], "]"}], "//", "GetFullDefinition"}]}]], "Input", CellChangeTimes->{{3.8273857615661697`*^9, 3.827385775513524*^9}, { 3.878815624913522*^9, 3.87881564197219*^9}}, CellLabel->"In[7]:=",ExpressionUUID->"78db9115-2c33-47e1-b9f4-34989752c2ec"], Cell[BoxData[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"k", "=", "1"}], "n"], RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Omicron]", "1"], "=", "1"}], "k"], FractionBox["1", SubscriptBox["\[Omicron]", "1"]]}]}]], "Output", CellChangeTimes->{ 3.878814314109256*^9, {3.878815618457699*^9, 3.878815642344716*^9}, 3.9229484471056433`*^9, 3.922948537172584*^9, 3.9241549293857183`*^9}, CellLabel->"Out[7]=",ExpressionUUID->"a99b1899-fec8-4c7b-9617-62a02b14af7c"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"res", "=", RowBox[{"SigmaReduce", "[", "mySum", "]"}]}]], "Input", CellChangeTimes->{{3.827385777143202*^9, 3.827385781444523*^9}, { 3.878815619469287*^9, 3.878815620566596*^9}}, CellLabel->"In[8]:=",ExpressionUUID->"1bd3158d-e04a-4c6d-97a4-5190ba4b1517"], Cell[BoxData[ RowBox[{ RowBox[{"-", "n"}], "+", RowBox[{ RowBox[{"(", RowBox[{"1", "+", "n"}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubscriptBox["\[Iota]", "1"]]}], ")"}]}]}]], "Output", CellChangeTimes->{ 3.878814317588065*^9, {3.878815621489141*^9, 3.878815647708351*^9}, 3.9229484485396223`*^9, 3.922948538469617*^9, 3.9241549304474792`*^9}, CellLabel->"Out[8]=",ExpressionUUID->"59e082d6-8be1-4556-80fd-5bcdcaec6e4b"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"mySum", "=", RowBox[{"SigmaSum", "[", RowBox[{ RowBox[{ RowBox[{"SigmaPower", "[", RowBox[{ RowBox[{"-", "1"}], ",", "k"}], "]"}], RowBox[{"SigmaBinomial", "[", RowBox[{"m", ",", "k"}], "]"}], RowBox[{"SigmaHNumber", "[", "k", "]"}]}], ",", RowBox[{"{", RowBox[{"k", ",", "0", ",", "n"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.707537540419404*^9, 3.707537606982614*^9}, { 3.878814362136804*^9, 3.878814363142992*^9}, {3.878814449440351*^9, 3.878814464949389*^9}}, CellLabel->"In[9]:=",ExpressionUUID->"024523ba-9a53-46ef-9fc2-368a2cea21a4"], Cell[BoxData[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"k", "=", "0"}], "n"], RowBox[{ SubscriptBox["H", "k"], " ", SuperscriptBox[ RowBox[{"(", GridBox[{ {"m"}, {"k"} }], ")"}], "."], " ", SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], UnderscriptBox["k", "."]]}]}]], "Output", CellChangeTimes->{ 3.878814330146905*^9, 3.878814363790163*^9, {3.878814450561412*^9, 3.878814465637574*^9}, 3.922948462197253*^9, 3.922948539368691*^9, 3.9241549313064003`*^9}, CellLabel->"Out[9]=",ExpressionUUID->"c1532054-d9ea-49ca-bffb-11a7534c1566"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"g", "=", RowBox[{"SigmaReduce", "[", "mySum", "]"}]}]], "Input", CellChangeTimes->{{3.707537566294983*^9, 3.707537571356398*^9}, { 3.878814367795557*^9, 3.878814368835568*^9}}, CellLabel->"In[10]:=",ExpressionUUID->"4f783cb6-5582-445f-9b9e-922cb4ad4ecf"], Cell[BoxData[ RowBox[{ RowBox[{"-", FractionBox["1", "m"]}], "+", FractionBox[ RowBox[{ RowBox[{"(", RowBox[{"m", "-", "n"}], ")"}], " ", SuperscriptBox[ RowBox[{"(", GridBox[{ {"m"}, {"n"} }], ")"}], "."], " ", SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], UnderscriptBox["n", "."]]}], SuperscriptBox["m", "2"]], "+", FractionBox[ RowBox[{ RowBox[{"(", RowBox[{"m", "-", "n"}], ")"}], " ", SubscriptBox["H", "n"], " ", SuperscriptBox[ RowBox[{"(", GridBox[{ {"m"}, {"n"} }], ")"}], "."], " ", SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], UnderscriptBox["n", "."]]}], "m"]}]], "Output", CellChangeTimes->{{3.878814339969549*^9, 3.8788143699598703`*^9}, { 3.878814454139722*^9, 3.878814467592424*^9}, 3.9229484637779827`*^9, 3.922948541514468*^9, 3.924154933050496*^9}, CellLabel->"Out[10]=",ExpressionUUID->"1718b7eb-4ec0-4b67-8b7f-1c8b34d77660"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"g", "=", RowBox[{ RowBox[{"CollectProdSum", "[", RowBox[{ RowBox[{"NShift", "[", RowBox[{"g", ",", "n", ",", RowBox[{"-", "1"}]}], "]"}], ",", "3"}], "]"}], "/.", RowBox[{"n", "->", "k"}]}]}]], "Input", CellChangeTimes->{{3.878814373399652*^9, 3.878814420196557*^9}}, CellLabel->"In[11]:=",ExpressionUUID->"e4681773-75c4-4916-8f62-2ed964138f80"], Cell[BoxData[ RowBox[{ RowBox[{"-", FractionBox["1", "m"]}], "+", FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "k"}], "+", "m"}], ")"}], " ", SuperscriptBox[ RowBox[{"(", GridBox[{ {"m"}, {"k"} }], ")"}], "."], " ", SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], UnderscriptBox["k", "."]]}], SuperscriptBox["m", "2"]], "-", FractionBox[ RowBox[{"k", " ", SubscriptBox["H", "k"], " ", SuperscriptBox[ RowBox[{"(", GridBox[{ {"m"}, {"k"} }], ")"}], "."], " ", SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], UnderscriptBox["k", "."]]}], "m"]}]], "Output", CellChangeTimes->{{3.878814397674767*^9, 3.878814420599609*^9}, 3.878814469889485*^9, 3.9229485425494757`*^9, 3.92415493339142*^9}, CellLabel->"Out[11]=",ExpressionUUID->"a2b6656e-2628-48a2-8d67-20d972b3c1e8"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"CollectProdSum", "[", RowBox[{ RowBox[{ RowBox[{"NShift", "[", RowBox[{"g", ",", "k"}], "]"}], "-", "g"}], ",", "3"}], "]"}]], "Input", CellChangeTimes->{{3.878814423647682*^9, 3.878814439680139*^9}}, CellLabel->"In[12]:=",ExpressionUUID->"9954db1d-25e9-43b4-b08d-ffbed6373fda"], Cell[BoxData[ RowBox[{ SubscriptBox["H", "k"], " ", SuperscriptBox[ RowBox[{"(", GridBox[{ {"m"}, {"k"} }], ")"}], "."], " ", SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], UnderscriptBox["k", "."]]}]], "Output", CellChangeTimes->{{3.878814432336053*^9, 3.878814440201775*^9}, 3.8788144752971067`*^9, 3.922948543972436*^9, 3.924154935177652*^9}, CellLabel->"Out[12]=",ExpressionUUID->"f1be8f99-3b30-4819-9d82-59d5a8d0de29"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"mySum", "=", RowBox[{"SigmaSum", "[", RowBox[{ RowBox[{ RowBox[{"SigmaSum", "[", RowBox[{ RowBox[{"SigmaBinomial", "[", RowBox[{"n", ",", "k"}], "]"}], ",", RowBox[{"{", RowBox[{"k", ",", "0", ",", "r"}], "}"}]}], "]"}], "^", "2"}], ",", RowBox[{"{", RowBox[{"r", ",", "0", ",", "b"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.7075375161061907`*^9, 3.707537538647422*^9}, { 3.7075376221477757`*^9, 3.707537647072652*^9}, {3.738925681621478*^9, 3.738925682307181*^9}}, CellLabel->"In[13]:=",ExpressionUUID->"dda35b47-bcb8-491d-b219-93075c117577"], Cell[BoxData[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"r", "=", "0"}], "b"], SuperscriptBox[ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"k", "=", "0"}], "r"], SuperscriptBox[ RowBox[{"(", GridBox[{ {"n"}, {"k"} }], ")"}], "."]}], ")"}], "2"]}]], "Output", CellChangeTimes->{3.87881458155307*^9, 3.922948545630704*^9, 3.92415493620114*^9}, CellLabel->"Out[13]=",ExpressionUUID->"e79a5fe9-1243-43db-80ab-88dfde90c07a"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"res", "=", RowBox[{"SigmaReduce", "[", "mySum", "]"}]}]], "Input", CellChangeTimes->{{3.707537648524673*^9, 3.707537671751049*^9}}, CellLabel->"In[14]:=",ExpressionUUID->"9077afb3-d77f-49da-a29e-d16f411837e6"], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "b"}], "+", "n"}], ")"}], " ", SuperscriptBox[ RowBox[{"(", GridBox[{ {"n"}, {"b"} }], ")"}], "."], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "0"}], "b"], SuperscriptBox[ RowBox[{"(", GridBox[{ {"n"}, { SubscriptBox["\[Iota]", "1"]} }], ")"}], "."]}], ")"}]}], "+", RowBox[{ FractionBox["1", "2"], " ", RowBox[{"(", RowBox[{"2", "+", RowBox[{"2", " ", "b"}], "-", "n"}], ")"}], " ", SuperscriptBox[ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "0"}], "b"], SuperscriptBox[ RowBox[{"(", GridBox[{ {"n"}, { SubscriptBox["\[Iota]", "1"]} }], ")"}], "."]}], ")"}], "2"]}], "-", RowBox[{ FractionBox["1", "2"], " ", "n", " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "0"}], "b"], SuperscriptBox[ RowBox[{"(", SuperscriptBox[ RowBox[{"(", GridBox[{ {"n"}, { SubscriptBox["\[Iota]", "1"]} }], ")"}], "."], ")"}], "2"]}], ")"}]}]}]], "Output", CellChangeTimes->{3.8788145933452177`*^9, 3.922948547123674*^9, 3.924154937750104*^9}, CellLabel->"Out[14]=",ExpressionUUID->"e347dfcc-c9a2-45a8-8982-cd0eae13acce"] }, Open ]] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["Section 2: The term algebra", "Title", CellChangeTimes->{{3.827918391802244*^9, 3.8279184268944263`*^9}, { 3.828588864626494*^9, 3.82858886540082*^9}, {3.878814748468452*^9, 3.87881475025758*^9}},ExpressionUUID->"113b3435-56cb-42f5-aecf-\ 7be6a9eab464"], Cell[CellGroupData[{ Cell["Defining objects in Sigma", "Subsection", CellChangeTimes->{{3.82791899333857*^9, 3.827918996968924*^9}, { 3.828588877794149*^9, 3.828588878159463*^9}, {3.828588958025045*^9, 3.828588958589932*^9}},ExpressionUUID->"f91662f3-4ffb-4cbd-a924-\ 61fa36336114"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"P", "=", RowBox[{ RowBox[{ RowBox[{"SigmaProduct", "[", RowBox[{"k", ",", RowBox[{"{", RowBox[{"k", ",", "1", ",", "n"}], "}"}]}], "]"}], "^", "2"}], RowBox[{"SigmaPower", "[", RowBox[{ RowBox[{"-", "1"}], ",", "n"}], "]"}]}]}]], "Input", CellChangeTimes->{{3.8279186462032137`*^9, 3.827918663545827*^9}, 3.8279187966991787`*^9, {3.827918901819124*^9, 3.827918902712126*^9}}, CellLabel->"In[2]:=",ExpressionUUID->"66c7f248-105c-409f-bdb1-b6cf4f480490"], Cell[BoxData[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], UnderscriptBox["n", "."]], " ", SuperscriptBox[ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Product]", "."], RowBox[{"k", "=", "1"}], "n"], "k"}], ")"}], "2"]}]], "Output", CellChangeTimes->{3.8316953667976522`*^9}, CellLabel->"Out[2]=",ExpressionUUID->"123dcdcf-9dbd-4aba-ae36-f9dad4ae9112"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"E1", "=", RowBox[{"SigmaSum", "[", RowBox[{ RowBox[{"SigmaFactorial", "[", "k", "]"}], ",", RowBox[{"{", RowBox[{"k", ",", "1", ",", "n"}], "}"}]}], "]"}]}], "\n"}]], "Input", CellChangeTimes->{{3.8279184659994087`*^9, 3.827918476668209*^9}}, CellLabel->"In[3]:=",ExpressionUUID->"cdf2189f-2fab-401b-b60e-9bf6d818c2a8"], Cell[BoxData[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"k", "=", "1"}], "n"], RowBox[{"k", SubscriptBox["!", "."]}]}]], "Output", CellChangeTimes->{3.83169536692137*^9}, CellLabel->"Out[3]=",ExpressionUUID->"c5e36c55-d128-45d2-a655-77c3b55125c5"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"GetFullDefinition", "[", "E1", "]"}]], "Input", CellChangeTimes->{{3.827918480157784*^9, 3.827918486774253*^9}}, CellLabel->"In[4]:=",ExpressionUUID->"42a9e20e-5ee9-4507-b7e1-b963b309eaf5"], Cell[BoxData[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"k", "=", "1"}], "n"], RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Product]", "."], RowBox[{ SubscriptBox["\[Omicron]", "1"], "=", "1"}], "k"], SubscriptBox["\[Omicron]", "1"]}]}]], "Output", CellChangeTimes->{3.831695366952636*^9}, CellLabel->"Out[4]=",ExpressionUUID->"0472079d-4317-4721-9d1f-52ddae94e364"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"E1", "/.", RowBox[{"n", "\[Rule]", "5"}]}]], "Input", CellChangeTimes->{{3.8279184906551867`*^9, 3.827918494834299*^9}}, CellLabel->"In[5]:=",ExpressionUUID->"8b877440-fdfa-4e1b-ad06-8d412e7e8cc9"], Cell[BoxData["153"], "Output", CellChangeTimes->{3.831695366977838*^9}, CellLabel->"Out[5]=",ExpressionUUID->"b9aca63a-4c27-4fda-955c-30f9d50ba688"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"E2", "=", RowBox[{"SigmaSum", "[", RowBox[{ RowBox[{ RowBox[{"SigmaSum", "[", RowBox[{ RowBox[{"1", "/", "i"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "k"}], "}"}]}], "]"}], " ", RowBox[{ RowBox[{"SigmaSum", "[", RowBox[{ RowBox[{"1", "/", RowBox[{"i", "^", "3"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "k"}], "}"}]}], "]"}], "/", RowBox[{"(", RowBox[{"k", "+", "1"}], ")"}]}]}], ",", RowBox[{"{", RowBox[{"k", ",", "1", ",", "n"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.827918506114643*^9, 3.827918529145743*^9}}, CellLabel->"In[6]:=",ExpressionUUID->"81fbce5f-089e-41d5-9823-956c744a20d2"], Cell[BoxData[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"k", "=", "1"}], "n"], FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "k"], FractionBox["1", SuperscriptBox["i", "3"]]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "k"], FractionBox["1", "i"]}], ")"}]}], RowBox[{"1", "+", "k"}]]}]], "Output", CellChangeTimes->{3.831695367008788*^9}, CellLabel->"Out[6]=",ExpressionUUID->"f1dfb6bc-862c-49ba-8309-4b4f463abe1b"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"E2", "/.", RowBox[{"n", "\[Rule]", "4"}]}]], "Input", CellChangeTimes->{{3.8279185334918528`*^9, 3.8279185387826567`*^9}}, CellLabel->"In[7]:=",ExpressionUUID->"65536035-af1f-4c40-bdd4-16bb1fbe157a"], Cell[BoxData[ FractionBox["14417", "6912"]], "Output", CellChangeTimes->{3.831695367044433*^9}, CellLabel->"Out[7]=",ExpressionUUID->"4e848326-e9b8-472d-9665-0a38fe148ec3"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"E3", "=", RowBox[{"E1", "*", RowBox[{"(", RowBox[{"E1", "+", "E2"}], ")"}]}]}]], "Input", CellChangeTimes->{{3.827918567154867*^9, 3.827918597599455*^9}}, CellLabel->"In[8]:=",ExpressionUUID->"6ef87d02-6ba2-4e65-b5ff-680037702a17"], Cell[BoxData[ RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"k", "=", "1"}], "n"], RowBox[{"k", SubscriptBox["!", "."]}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"k", "=", "1"}], "n"], RowBox[{"k", SubscriptBox["!", "."]}]}], "+", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"k", "=", "1"}], "n"], FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "k"], FractionBox["1", SuperscriptBox["i", "3"]]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "k"], FractionBox["1", "i"]}], ")"}]}], RowBox[{"1", "+", "k"}]]}]}], ")"}]}]], "Output", CellChangeTimes->{3.8316953670795794`*^9}, CellLabel->"Out[8]=",ExpressionUUID->"97777b37-2ebd-4a3b-a9eb-47627b686e9a"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"E3", "/.", RowBox[{"n", "\[Rule]", "4"}]}]], "Input", CellChangeTimes->{{3.82791860025425*^9, 3.827918602903414*^9}}, CellLabel->"In[9]:=",ExpressionUUID->"75005778-1b7e-4f08-8a99-e8b09e75f474"], Cell[BoxData[ FractionBox["2667643", "2304"]], "Output", CellChangeTimes->{3.83169536711758*^9}, CellLabel->"Out[9]=",ExpressionUUID->"a8d52c27-ab66-46c8-b67b-30704d278899"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"E3", "=", RowBox[{"CollectProdSum", "[", "E3", "]"}]}]], "Input", CellChangeTimes->{{3.82791860483445*^9, 3.82791861233134*^9}}, CellLabel->"In[10]:=",ExpressionUUID->"51c5aabc-c553-471c-a146-99226db20810"], Cell[BoxData[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"k", "=", "1"}], "n"], RowBox[{"k", SubscriptBox["!", "."]}]}], ")"}], "2"], "+", RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"k", "=", "1"}], "n"], RowBox[{"k", SubscriptBox["!", "."]}]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"k", "=", "1"}], "n"], FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "k"], FractionBox["1", SuperscriptBox["i", "3"]]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "k"], FractionBox["1", "i"]}], ")"}]}], RowBox[{"1", "+", "k"}]]}], ")"}]}]}]], "Output", CellChangeTimes->{3.83169536716084*^9}, CellLabel->"Out[10]=",ExpressionUUID->"4b9b0f33-f17d-487d-bd07-e71ea0beee7f"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["The problem SigmaReduce", "Subsection", CellChangeTimes->{{3.827919006098606*^9, 3.827919009782981*^9}},ExpressionUUID->"1d2860c1-3e4a-40bc-855a-\ 7423b34aa01a"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"A1", "=", RowBox[{"SigmaSum", "[", RowBox[{ RowBox[{ RowBox[{"SigmaSum", "[", RowBox[{ RowBox[{"1", "/", "i"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "k"}], "}"}]}], "]"}], RowBox[{ RowBox[{"SigmaSum", "[", RowBox[{ RowBox[{"1", "/", RowBox[{"i", "^", "3"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "k"}], "}"}]}], "]"}], "/", RowBox[{"(", RowBox[{"k", "+", "1"}], ")"}]}]}], ",", RowBox[{"{", RowBox[{"k", ",", "1", ",", "n"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.8191765609624443`*^9, 3.819176601254879*^9}}, CellLabel->"In[11]:=",ExpressionUUID->"425a78a9-cb98-45fb-a6e6-78a30511735c"], Cell[BoxData[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"k", "=", "1"}], "n"], FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "k"], FractionBox["1", SuperscriptBox["i", "3"]]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "k"], FractionBox["1", "i"]}], ")"}]}], RowBox[{"1", "+", "k"}]]}]], "Output", CellChangeTimes->{3.831695367199843*^9}, CellLabel->"Out[11]=",ExpressionUUID->"368a0dfe-252e-4ae6-a096-0c768fdf5645"] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"A2", "=", RowBox[{ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "n"], FractionBox["1", SuperscriptBox["i", "5"]]}], "-", FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "n"], FractionBox["1", SuperscriptBox["i", "4"]]}], RowBox[{"1", "+", "n"}]], "-", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"j", "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "j"], FractionBox["1", SuperscriptBox["i", "4"]]}], "j"]}], "-", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"j", "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "j"], FractionBox["1", SuperscriptBox["i", "3"]]}], SuperscriptBox["j", "2"]]}], "+", FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"j", "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "j"], FractionBox["1", SuperscriptBox["i", "3"]]}], "j"]}], RowBox[{"1", "+", "n"}]], "-", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"j", "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "j"], FractionBox["1", "i"]}], SuperscriptBox["j", "4"]]}], "+", FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"j", "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "j"], FractionBox["1", "i"]}], SuperscriptBox["j", "3"]]}], RowBox[{"1", "+", "n"}]], "+", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"k", "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"j", "=", "1"}], "k"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "j"], FractionBox["1", SuperscriptBox["i", "3"]]}], "j"]}], "k"]}], "+", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"k", "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"j", "=", "1"}], "k"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "j"], FractionBox["1", "i"]}], SuperscriptBox["j", "3"]]}], "k"]}]}]}], ";"}]], "Input", CellChangeTimes->{{3.819176449434576*^9, 3.819176467368972*^9}}, CellLabel->"In[12]:=",ExpressionUUID->"58e7a766-98ca-4a0a-989f-688b3536f48d"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"{", RowBox[{"A1", ",", "A2"}], "}"}], "/.", RowBox[{"n", "\[Rule]", "4"}]}]], "Input", CellChangeTimes->{{3.8279182816735363`*^9, 3.82791828982588*^9}}, CellLabel->"In[13]:=",ExpressionUUID->"bfe104f0-6b23-4c4c-99d3-3ace36728362"], Cell[BoxData[ RowBox[{"{", RowBox[{ FractionBox["14417", "6912"], ",", FractionBox["14417", "6912"]}], "}"}]], "Output", CellChangeTimes->{3.8316953672577963`*^9}, CellLabel->"Out[13]=",ExpressionUUID->"9900d1b6-8b38-404b-a6d4-8e6d6fa2a77b"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"{", RowBox[{"B1", ",", "B2"}], "}"}], "=", RowBox[{"SigmaReduce", "[", RowBox[{ RowBox[{"{", RowBox[{"A1", ",", "A2"}], "}"}], ",", "n"}], "]"}]}]], "Input", CellChangeTimes->{{3.8191764272075377`*^9, 3.819176433637816*^9}, { 3.8191764743297367`*^9, 3.819176492127537*^9}}, CellLabel->"In[14]:=",ExpressionUUID->"7a5fd869-602b-493a-bccb-3fa172ee7ac3"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "3"]]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubscriptBox["\[Iota]", "2"]]}], ")"}]}], RowBox[{"1", "+", SubscriptBox["\[Iota]", "1"]}]]}], ",", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "3"]]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubscriptBox["\[Iota]", "2"]]}], ")"}]}], RowBox[{"1", "+", SubscriptBox["\[Iota]", "1"]}]]}]}], "}"}]], "Output", CellChangeTimes->{3.8316953696620197`*^9}, CellLabel->"Out[14]=",ExpressionUUID->"16d49442-504a-470a-9b88-d2ac845b79c7"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"W", "=", RowBox[{"Complement", "[", RowBox[{ RowBox[{"Variables", "[", RowBox[{"NShift", "[", RowBox[{ RowBox[{"{", RowBox[{"B1", ",", "B2"}], "}"}], ",", "n"}], "]"}], "]"}], ",", RowBox[{"{", "n", "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.8279190182305326`*^9, 3.827919056649541*^9}}, CellLabel->"In[15]:=",ExpressionUUID->"400fea04-d5f6-42ef-b976-6903a446c250"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubsuperscriptBox["\[Iota]", "1", "3"]]}], ",", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubscriptBox["\[Iota]", "1"]]}], ",", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "3"]]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubscriptBox["\[Iota]", "2"]]}], ")"}]}], RowBox[{"1", "+", SubscriptBox["\[Iota]", "1"]}]]}]}], "}"}]], "Output", CellChangeTimes->{3.83169536976283*^9}, CellLabel->"Out[15]=",ExpressionUUID->"07b473cc-8efb-491f-a597-9b22396b4176"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"{", RowBox[{"B2", ",", "B1"}], "}"}], "=", RowBox[{"SigmaReduce", "[", RowBox[{ RowBox[{"{", RowBox[{"A2", ",", "A1"}], "}"}], ",", "n"}], "]"}]}]], "Input", CellChangeTimes->{{3.8191765075745573`*^9, 3.8191765130240917`*^9}, { 3.819177279056683*^9, 3.819177287290031*^9}}, CellLabel->"In[16]:=",ExpressionUUID->"d3d91d94-9cef-4cbc-9ad1-5f13021da5c1"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubsuperscriptBox["\[Iota]", "1", "4"]]}], ")"}]}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubscriptBox["\[Iota]", "1"]]}], ")"}]}], "+", FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubsuperscriptBox["\[Iota]", "1", "3"]]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubscriptBox["\[Iota]", "1"]]}], ")"}]}], RowBox[{"1", "+", "n"}]], "+", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "4"]]}], SubscriptBox["\[Iota]", "1"]]}], "-", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "3"]]}], SubsuperscriptBox["\[Iota]", "1", "2"]]}], "+", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "3"]]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubscriptBox["\[Iota]", "2"]]}], ")"}]}], SubscriptBox["\[Iota]", "1"]]}]}], ",", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubsuperscriptBox["\[Iota]", "1", "4"]]}], ")"}]}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubscriptBox["\[Iota]", "1"]]}], ")"}]}], "+", FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubsuperscriptBox["\[Iota]", "1", "3"]]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubscriptBox["\[Iota]", "1"]]}], ")"}]}], RowBox[{"1", "+", "n"}]], "+", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "4"]]}], SubscriptBox["\[Iota]", "1"]]}], "-", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "3"]]}], SubsuperscriptBox["\[Iota]", "1", "2"]]}], "+", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "3"]]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubscriptBox["\[Iota]", "2"]]}], ")"}]}], SubscriptBox["\[Iota]", "1"]]}]}]}], "}"}]], "Output", CellChangeTimes->{3.831695371960094*^9}, CellLabel->"Out[16]=",ExpressionUUID->"76b24814-2a55-4261-bf20-3ff1e939738e"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"W", "=", RowBox[{"Complement", "[", RowBox[{ RowBox[{"Variables", "[", RowBox[{"NShift", "[", RowBox[{ RowBox[{"{", RowBox[{"B1", ",", "B2"}], "}"}], ",", "n"}], "]"}], "]"}], ",", RowBox[{"{", "n", "}"}]}], "]"}]}]], "Input", CellLabel->"In[17]:=",ExpressionUUID->"c5ac63ec-7e18-4bd9-b5de-e72811199250"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubsuperscriptBox["\[Iota]", "1", "4"]]}], ",", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubsuperscriptBox["\[Iota]", "1", "3"]]}], ",", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubscriptBox["\[Iota]", "1"]]}], ",", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "4"]]}], SubscriptBox["\[Iota]", "1"]]}], ",", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "3"]]}], SubsuperscriptBox["\[Iota]", "1", "2"]]}], ",", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "3"]]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubscriptBox["\[Iota]", "2"]]}], ")"}]}], SubscriptBox["\[Iota]", "1"]]}]}], "}"}]], "Output", CellChangeTimes->{3.831695372103735*^9}, CellLabel->"Out[17]=",ExpressionUUID->"2dc6fcd8-d9de-46e9-ab43-9ccfee7cd2c2"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"SigmaReduce", "[", RowBox[{ RowBox[{"A2", "-", "A1"}], ",", "n"}], "]"}]], "Input", CellChangeTimes->{{3.819176522007575*^9, 3.819176530720066*^9}}, CellLabel->"In[18]:=",ExpressionUUID->"652ed41e-1888-41fe-92f9-a7227b0570e3"], Cell[BoxData["0"], "Output", CellChangeTimes->{3.831695374313547*^9}, CellLabel->"Out[18]=",ExpressionUUID->"b244ca41-0c24-4f6b-b82e-759db9ff15ea"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"varL", "=", RowBox[{"Append", "[", RowBox[{ RowBox[{"Variables", "[", "A2", "]"}], ",", "A1"}], "]"}]}]], "Input", CellChangeTimes->{{3.819177999798284*^9, 3.81917805464017*^9}}, CellLabel->"In[19]:=",ExpressionUUID->"9222af92-c3f1-4086-bccb-ac14201c6226"], Cell[BoxData[ RowBox[{"{", RowBox[{"n", ",", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "n"], FractionBox["1", SuperscriptBox["i", "5"]]}], ",", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "n"], FractionBox["1", SuperscriptBox["i", "4"]]}], ",", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"j", "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "j"], FractionBox["1", SuperscriptBox["i", "4"]]}], "j"]}], ",", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"j", "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "j"], FractionBox["1", SuperscriptBox["i", "3"]]}], SuperscriptBox["j", "2"]]}], ",", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"j", "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "j"], FractionBox["1", SuperscriptBox["i", "3"]]}], "j"]}], ",", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"j", "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "j"], FractionBox["1", "i"]}], SuperscriptBox["j", "4"]]}], ",", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"j", "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "j"], FractionBox["1", "i"]}], SuperscriptBox["j", "3"]]}], ",", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"k", "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"j", "=", "1"}], "k"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "j"], FractionBox["1", SuperscriptBox["i", "3"]]}], "j"]}], "k"]}], ",", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"k", "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"j", "=", "1"}], "k"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "j"], FractionBox["1", "i"]}], SuperscriptBox["j", "3"]]}], "k"]}], ",", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"k", "=", "1"}], "n"], FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "k"], FractionBox["1", SuperscriptBox["i", "3"]]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "k"], FractionBox["1", "i"]}], ")"}]}], RowBox[{"1", "+", "k"}]]}]}], "}"}]], "Output", CellChangeTimes->{3.831695374362125*^9}, CellLabel->"Out[19]=",ExpressionUUID->"228cbb61-83a6-40b9-84a5-692821224408"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"varLRed", "=", RowBox[{"SigmaReduce", "[", RowBox[{"varL", ",", "n"}], "]"}]}]], "Input", CellChangeTimes->{{3.8191780569386044`*^9, 3.8191780652460938`*^9}, { 3.8191781369521112`*^9, 3.819178141149972*^9}, {3.8191783217720127`*^9, 3.819178322418404*^9}, {3.827918315557498*^9, 3.827918317280294*^9}}, CellLabel->"In[20]:=",ExpressionUUID->"d21be423-8a6b-42a8-9a21-fe498cc0479f"], Cell[BoxData[ RowBox[{"{", RowBox[{"n", ",", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubsuperscriptBox["\[Iota]", "1", "5"]]}], ",", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubsuperscriptBox["\[Iota]", "1", "4"]]}], ",", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "4"]]}], SubscriptBox["\[Iota]", "1"]]}], ",", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "3"]]}], SubsuperscriptBox["\[Iota]", "1", "2"]]}], ",", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "3"]]}], SubscriptBox["\[Iota]", "1"]]}], ",", RowBox[{ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubsuperscriptBox["\[Iota]", "1", "5"]]}], "+", RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubsuperscriptBox["\[Iota]", "1", "4"]]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubscriptBox["\[Iota]", "1"]]}], ")"}]}], "-", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "4"]]}], SubscriptBox["\[Iota]", "1"]]}]}], ",", RowBox[{ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubsuperscriptBox["\[Iota]", "1", "4"]]}], "+", RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubsuperscriptBox["\[Iota]", "1", "3"]]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubscriptBox["\[Iota]", "1"]]}], ")"}]}], "-", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "3"]]}], SubscriptBox["\[Iota]", "1"]]}]}], ",", RowBox[{ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "3"]]}], SubsuperscriptBox["\[Iota]", "1", "2"]]}], "+", RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubscriptBox["\[Iota]", "1"]]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "3"]]}], SubscriptBox["\[Iota]", "1"]]}], ")"}]}], "-", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "3"]]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubscriptBox["\[Iota]", "2"]]}], ")"}]}], SubscriptBox["\[Iota]", "1"]]}]}], ",", RowBox[{ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "4"]]}], SubscriptBox["\[Iota]", "1"]]}], "-", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "3"]]}], SubsuperscriptBox["\[Iota]", "1", "2"]]}], "-", RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubscriptBox["\[Iota]", "1"]]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "3"]]}], SubscriptBox["\[Iota]", "1"]]}], ")"}]}], "+", RowBox[{"2", " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "3"]]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubscriptBox["\[Iota]", "2"]]}], ")"}]}], SubscriptBox["\[Iota]", "1"]]}], ")"}]}]}], ",", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubsuperscriptBox["\[Iota]", "1", "4"]]}], ")"}]}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubscriptBox["\[Iota]", "1"]]}], ")"}]}], "+", FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubsuperscriptBox["\[Iota]", "1", "3"]]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubscriptBox["\[Iota]", "1"]]}], ")"}]}], RowBox[{"1", "+", "n"}]], "+", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "4"]]}], SubscriptBox["\[Iota]", "1"]]}], "-", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "3"]]}], SubsuperscriptBox["\[Iota]", "1", "2"]]}], "+", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "3"]]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubscriptBox["\[Iota]", "2"]]}], ")"}]}], SubscriptBox["\[Iota]", "1"]]}]}]}], "}"}]], "Output", CellChangeTimes->{3.831695376611923*^9}, CellLabel->"Out[20]=",ExpressionUUID->"1fe56d23-7a38-4012-983a-ecd4232ad48e"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"MapThread", "[", RowBox[{"Rule", ",", RowBox[{"{", RowBox[{"varL", ",", "varLRed"}], "}"}]}], "]"}], "//", "MatrixForm"}]], "Input", CellChangeTimes->{{3.819178072888908*^9, 3.8191781194516487`*^9}, { 3.819178336757028*^9, 3.819178356396052*^9}, {3.827918319762723*^9, 3.827918321048726*^9}}, CellLabel->"In[21]:=",ExpressionUUID->"4bbb293e-5b5d-4db1-bd39-ee5e7a73e766"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", TagBox[GridBox[{ { RowBox[{"n", "\[Rule]", "n"}]}, { RowBox[{ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "n"], FractionBox["1", SuperscriptBox["i", "5"]]}], "\[Rule]", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubsuperscriptBox["\[Iota]", "1", "5"]]}]}]}, { RowBox[{ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "n"], FractionBox["1", SuperscriptBox["i", "4"]]}], "\[Rule]", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubsuperscriptBox["\[Iota]", "1", "4"]]}]}]}, { RowBox[{ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"j", "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "j"], FractionBox["1", SuperscriptBox["i", "4"]]}], "j"]}], "\[Rule]", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "4"]]}], SubscriptBox["\[Iota]", "1"]]}]}]}, { RowBox[{ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"j", "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "j"], FractionBox["1", SuperscriptBox["i", "3"]]}], SuperscriptBox["j", "2"]]}], "\[Rule]", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "3"]]}], SubsuperscriptBox["\[Iota]", "1", "2"]]}]}]}, { RowBox[{ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"j", "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "j"], FractionBox["1", SuperscriptBox["i", "3"]]}], "j"]}], "\[Rule]", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "3"]]}], SubscriptBox["\[Iota]", "1"]]}]}]}, { RowBox[{ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"j", "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "j"], FractionBox["1", "i"]}], SuperscriptBox["j", "4"]]}], "\[Rule]", RowBox[{ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubsuperscriptBox["\[Iota]", "1", "5"]]}], "+", RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubsuperscriptBox["\[Iota]", "1", "4"]]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubscriptBox["\[Iota]", "1"]]}], ")"}]}], "-", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "4"]]}], SubscriptBox["\[Iota]", "1"]]}]}]}]}, { RowBox[{ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"j", "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "j"], FractionBox["1", "i"]}], SuperscriptBox["j", "3"]]}], "\[Rule]", RowBox[{ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubsuperscriptBox["\[Iota]", "1", "4"]]}], "+", RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubsuperscriptBox["\[Iota]", "1", "3"]]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubscriptBox["\[Iota]", "1"]]}], ")"}]}], "-", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "3"]]}], SubscriptBox["\[Iota]", "1"]]}]}]}]}, { RowBox[{ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"k", "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"j", "=", "1"}], "k"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "j"], FractionBox["1", SuperscriptBox["i", "3"]]}], "j"]}], "k"]}], "\[Rule]", RowBox[{ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "3"]]}], SubsuperscriptBox["\[Iota]", "1", "2"]]}], "+", RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubscriptBox["\[Iota]", "1"]]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "3"]]}], SubscriptBox["\[Iota]", "1"]]}], ")"}]}], "-", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "3"]]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubscriptBox["\[Iota]", "2"]]}], ")"}]}], SubscriptBox["\[Iota]", "1"]]}]}]}]}, { RowBox[{ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"k", "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"j", "=", "1"}], "k"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "j"], FractionBox["1", "i"]}], SuperscriptBox["j", "3"]]}], "k"]}], "\[Rule]", RowBox[{ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "4"]]}], SubscriptBox["\[Iota]", "1"]]}], "-", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "3"]]}], SubsuperscriptBox["\[Iota]", "1", "2"]]}], "-", RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubscriptBox["\[Iota]", "1"]]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "3"]]}], SubscriptBox["\[Iota]", "1"]]}], ")"}]}], "+", RowBox[{"2", " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "3"]]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubscriptBox["\[Iota]", "2"]]}], ")"}]}], SubscriptBox["\[Iota]", "1"]]}], ")"}]}]}]}]}, { RowBox[{ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"k", "=", "1"}], "n"], FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "k"], FractionBox["1", SuperscriptBox["i", "3"]]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"i", "=", "1"}], "k"], FractionBox["1", "i"]}], ")"}]}], RowBox[{"1", "+", "k"}]]}], "\[Rule]", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubsuperscriptBox["\[Iota]", "1", "4"]]}], ")"}]}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubscriptBox["\[Iota]", "1"]]}], ")"}]}], "+", FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubsuperscriptBox["\[Iota]", "1", "3"]]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubscriptBox["\[Iota]", "1"]]}], ")"}]}], RowBox[{"1", "+", "n"}]], "+", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "4"]]}], SubscriptBox["\[Iota]", "1"]]}], "-", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "3"]]}], SubsuperscriptBox["\[Iota]", "1", "2"]]}], "+", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubsuperscriptBox["\[Iota]", "2", "3"]]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubscriptBox["\[Iota]", "2"]]}], ")"}]}], SubscriptBox["\[Iota]", "1"]]}]}]}]} }, GridBoxAlignment->{"Columns" -> {{Center}}, "Rows" -> {{Baseline}}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.5599999999999999]}, Offset[0.27999999999999997`]}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}}], Column], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output", CellChangeTimes->{3.8316953766882563`*^9}, CellLabel-> "Out[21]//MatrixForm=",ExpressionUUID->"af68bb52-2b2a-4f0a-b353-\ 181325317753"] }, Open ]] }, Closed]] }, Closed]], Cell[CellGroupData[{ Cell["Section 4: Symbolic summation", "Title", CellChangeTimes->{{3.827918391802244*^9, 3.8279184268944263`*^9}, { 3.828588864626494*^9, 3.82858886540082*^9}, {3.878814748468452*^9, 3.87881475025758*^9}, {3.924156428734262*^9, 3.9241564412299547`*^9}},ExpressionUUID->"da6e2fb4-83d1-4243-8d7b-\ f595ca8fd4ba"], Cell[CellGroupData[{ Cell["Toy examples", "Section", CellChangeTimes->{{3.831622274198222*^9, 3.831622304898262*^9}},ExpressionUUID->"3ec7fa06-35d9-4269-aa2e-\ f49b20fefaa1"], Cell[CellGroupData[{ Cell["Telescoping", "Subsection", CellChangeTimes->{{3.83162427645835*^9, 3.831624277985805*^9}},ExpressionUUID->"6d9045c2-d4c1-4b53-8e37-\ 0ff603bd2631"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"mySum", "=", RowBox[{"SigmaSum", "[", RowBox[{ RowBox[{"S", "[", RowBox[{"1", ",", "k"}], "]"}], ",", RowBox[{"{", RowBox[{"k", ",", "1", ",", "n"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.831622287298766*^9, 3.8316223137350388`*^9}}, CellLabel->"In[22]:=",ExpressionUUID->"e2c8ef8e-002a-4f41-ab1b-0ba2d5a47ad7"], Cell[BoxData[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"k", "=", "1"}], "n"], RowBox[{"S", "[", RowBox[{"1", ",", "k"}], "]"}]}]], "Output", CellChangeTimes->{3.831636315810925*^9, 3.831695381778981*^9}, CellLabel->"Out[22]=",ExpressionUUID->"e79340c8-ffcf-4267-b19c-88f4c47a4253"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"SigmaReduce", "[", "mySum", "]"}]], "Input", CellChangeTimes->{{3.831624224425425*^9, 3.831624228276979*^9}}, CellLabel->"In[23]:=",ExpressionUUID->"2db5acb2-23dc-4674-9bbe-d5df57fbfd67"], Cell[BoxData[ RowBox[{ RowBox[{"-", "n"}], "+", RowBox[{ RowBox[{"(", RowBox[{"1", "+", "n"}], ")"}], " ", RowBox[{"S", "[", RowBox[{"1", ",", "n"}], "]"}]}]}]], "Output", CellChangeTimes->{3.831636323103984*^9, 3.831695382002977*^9}, CellLabel->"Out[23]=",ExpressionUUID->"62df2090-dead-4cc6-b9a4-ef6c6e4f9e78"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"mySum", "=", RowBox[{"SigmaSum", "[", RowBox[{ RowBox[{ RowBox[{"SigmaFactorial", "[", "k", "]"}], " ", "k"}], ",", RowBox[{"{", RowBox[{"k", ",", "a", ",", "b"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.831624203440839*^9, 3.831624211999302*^9}}, CellLabel->"In[24]:=",ExpressionUUID->"bc1178b4-b258-40cf-9be3-6bd7bf5f4bc2"], Cell[BoxData[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"k", "=", "a"}], "b"], RowBox[{"k", " ", RowBox[{"k", SubscriptBox["!", "."]}]}]}]], "Output", CellChangeTimes->{3.831636349065555*^9, 3.831695382034697*^9}, CellLabel->"Out[24]=",ExpressionUUID->"f7e2dcfb-f4f8-4bee-9c21-7a6dce0678d4"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"SigmaReduce", "[", "mySum", "]"}]], "Input", CellChangeTimes->{{3.831624203440839*^9, 3.8316242126586523`*^9}}, CellLabel->"In[25]:=",ExpressionUUID->"d88ff7ab-bb59-49a8-ba05-b944dce7359e"], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"-", "a"}], " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "a"}], ")"}], SubscriptBox["!", "."]}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"1", "+", "b"}], ")"}], " ", RowBox[{"b", SubscriptBox["!", "."]}]}]}]], "Output", CellChangeTimes->{3.83163635454921*^9, 3.83169538217455*^9}, CellLabel->"Out[25]=",ExpressionUUID->"10f60be6-a0ba-462b-a201-e4fe4a7ae019"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"mySum", "=", RowBox[{"SigmaSum", "[", RowBox[{ RowBox[{ RowBox[{"SigmaPower", "[", RowBox[{ RowBox[{"-", "1"}], ",", "k"}], "]"}], " ", RowBox[{"SigmaBinomial", "[", RowBox[{"n", ",", "k"}], "]"}], " ", RowBox[{"SigmaHNumber", "[", "k", "]"}]}], ",", RowBox[{"{", RowBox[{"k", ",", "a", ",", "b"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.831624203440839*^9, 3.831624213292328*^9}}, CellLabel->"In[26]:=",ExpressionUUID->"38e23792-f224-4bba-b914-3a3f637e23f4"], Cell[BoxData[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"k", "=", "a"}], "b"], RowBox[{ SubscriptBox["H", "k"], " ", SuperscriptBox[ RowBox[{"(", GridBox[{ {"n"}, {"k"} }], ")"}], "."], " ", SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], UnderscriptBox["k", "."]]}]}]], "Output", CellChangeTimes->{3.831636377164926*^9, 3.831695382207223*^9}, CellLabel->"Out[26]=",ExpressionUUID->"f7db102d-e3dc-4619-adfe-19f700ec29b4"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"SigmaReduce", "[", "mySum", "]"}]], "Input", CellChangeTimes->{{3.831624203440839*^9, 3.831624213913733*^9}}, CellLabel->"In[27]:=",ExpressionUUID->"245256e3-7de1-4f5e-8f29-50c5266baf76"], Cell[BoxData[ RowBox[{ RowBox[{"-", FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "a", "-", "n"}], ")"}], " ", RowBox[{"(", RowBox[{"a", "-", "n"}], ")"}], " ", SuperscriptBox[ RowBox[{"(", GridBox[{ {"n"}, { RowBox[{ RowBox[{"-", "1"}], "+", "a"}]} }], ")"}], "."], " ", SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], UnderscriptBox["a", "."]]}], RowBox[{"a", " ", SuperscriptBox["n", "2"]}]]}], "+", FractionBox[ RowBox[{ RowBox[{"(", RowBox[{"1", "-", "a", "+", "n"}], ")"}], " ", SubscriptBox["H", "a"], " ", SuperscriptBox[ RowBox[{"(", GridBox[{ {"n"}, { RowBox[{ RowBox[{"-", "1"}], "+", "a"}]} }], ")"}], "."], " ", SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], UnderscriptBox["a", "."]]}], "n"], "+", FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "b"}], "+", "n"}], ")"}], " ", SuperscriptBox[ RowBox[{"(", GridBox[{ {"n"}, {"b"} }], ")"}], "."], " ", SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], UnderscriptBox["b", "."]]}], SuperscriptBox["n", "2"]], "+", FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "b"}], "+", "n"}], ")"}], " ", SubscriptBox["H", "b"], " ", SuperscriptBox[ RowBox[{"(", GridBox[{ {"n"}, {"b"} }], ")"}], "."], " ", SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], UnderscriptBox["b", "."]]}], "n"]}]], "Output", CellChangeTimes->{3.831636388778604*^9, 3.831695382906118*^9}, CellLabel->"Out[27]=",ExpressionUUID->"9eca07ce-0c18-4801-b190-afcad63a17c7"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"mySum", "=", RowBox[{"SigmaSum", "[", RowBox[{ RowBox[{ RowBox[{"SigmaSum", "[", RowBox[{ RowBox[{"SigmaBinomial", "[", RowBox[{"n", ",", "k"}], "]"}], ",", RowBox[{"{", RowBox[{"k", ",", "0", ",", "r"}], "}"}]}], "]"}], "^", "2"}], ",", RowBox[{"{", RowBox[{"r", ",", "0", ",", "b"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.831624203440839*^9, 3.8316242145103703`*^9}}, CellLabel->"In[28]:=",ExpressionUUID->"af75c852-4958-4d3e-b6b3-691592942265"], Cell[BoxData[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"r", "=", "0"}], "b"], SuperscriptBox[ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"k", "=", "0"}], "r"], SuperscriptBox[ RowBox[{"(", GridBox[{ {"n"}, {"k"} }], ")"}], "."]}], ")"}], "2"]}]], "Output", CellChangeTimes->{3.83163640452861*^9, 3.831695383001484*^9}, CellLabel->"Out[28]=",ExpressionUUID->"101a6230-bc42-458a-9a69-b079b252321c"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"res", "=", RowBox[{"SigmaReduce", "[", "mySum", "]"}]}]], "Input", CellChangeTimes->{{3.831624203440839*^9, 3.8316242151100273`*^9}}, CellLabel->"In[29]:=",ExpressionUUID->"0d224b46-7d22-46ed-9dc4-6daea32b5c33"], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "b"}], "+", "n"}], ")"}], " ", SuperscriptBox[ RowBox[{"(", GridBox[{ {"n"}, {"b"} }], ")"}], "."], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "0"}], "b"], SuperscriptBox[ RowBox[{"(", GridBox[{ {"n"}, { SubscriptBox["\[Iota]", "1"]} }], ")"}], "."]}], ")"}]}], "+", RowBox[{ FractionBox["1", "2"], " ", RowBox[{"(", RowBox[{"2", "+", RowBox[{"2", " ", "b"}], "-", "n"}], ")"}], " ", SuperscriptBox[ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "0"}], "b"], SuperscriptBox[ RowBox[{"(", GridBox[{ {"n"}, { SubscriptBox["\[Iota]", "1"]} }], ")"}], "."]}], ")"}], "2"]}], "-", RowBox[{ FractionBox["1", "2"], " ", "n", " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "0"}], "b"], SuperscriptBox[ RowBox[{"(", SuperscriptBox[ RowBox[{"(", GridBox[{ {"n"}, { SubscriptBox["\[Iota]", "1"]} }], ")"}], "."], ")"}], "2"]}], ")"}]}]}]], "Output", CellChangeTimes->{3.8316364204483967`*^9, 3.831695383357109*^9}, CellLabel->"Out[29]=",ExpressionUUID->"42faf624-bc57-4fe2-8927-410a089ebe62"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"{", RowBox[{"mySum", "-", "res"}], "}"}], "/.", RowBox[{"{", RowBox[{"b", "\[Rule]", "10"}], "}"}]}], "//", "Together"}]], "Input", CellChangeTimes->{{3.831624203440839*^9, 3.8316242162011957`*^9}, { 3.831636446746657*^9, 3.8316364570174*^9}}, CellLabel->"In[30]:=",ExpressionUUID->"22cd9332-4904-4fc0-abd8-7924473c4ca0"], Cell[BoxData[ RowBox[{"{", "0", "}"}]], "Output", CellChangeTimes->{{3.8316364509481153`*^9, 3.831636457740253*^9}, 3.831695383440839*^9}, CellLabel->"Out[30]=",ExpressionUUID->"a305e6a9-8e53-4e96-bb12-56bfcfd75bc9"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Definite summation", "Subsection", CellChangeTimes->{{3.831624203440839*^9, 3.831624216849856*^9}, 3.8316242828476543`*^9},ExpressionUUID->"75de7937-3466-4af4-9301-\ 40848fbc0599"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"mySum", "=", RowBox[{"SigmaSum", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"SigmaHNumber", "[", "k", "]"}], "^", "2"}], " ", RowBox[{"SigmaBinomial", "[", RowBox[{"n", ",", "k"}], "]"}]}], ",", RowBox[{"{", RowBox[{"k", ",", "0", ",", "n"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.831624203440839*^9, 3.8316242174324083`*^9}, { 3.831695406843995*^9, 3.831695408795796*^9}}, CellLabel->"In[36]:=",ExpressionUUID->"0eac2d65-2dfc-445a-bf38-a5bc1a8bb50a"], Cell[BoxData[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"k", "=", "0"}], "n"], RowBox[{ SubsuperscriptBox["H", "k", "2"], " ", SuperscriptBox[ RowBox[{"(", GridBox[{ {"n"}, {"k"} }], ")"}], "."]}]}]], "Output", CellChangeTimes->{{3.831695383473031*^9, 3.831695409236664*^9}}, CellLabel->"Out[36]=",ExpressionUUID->"2cb7a3ec-2951-4b4d-b5e3-315735cccfc4"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"SigmaReduce", "[", "mySum", "]"}]], "Input", CellChangeTimes->{{3.831624203440839*^9, 3.831624218146708*^9}}, CellLabel->"In[37]:=",ExpressionUUID->"6396d67d-db28-48e3-a9b2-8b8bde677f01"], Cell[BoxData[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "0"}], "n"], RowBox[{ SubsuperscriptBox["H", SubscriptBox["\[Iota]", "1"], "2"], " ", SuperscriptBox[ RowBox[{"(", GridBox[{ {"n"}, { SubscriptBox["\[Iota]", "1"]} }], ")"}], "."]}]}]], "Output", CellChangeTimes->{{3.831695383656279*^9, 3.831695410404294*^9}}, CellLabel->"Out[37]=",ExpressionUUID->"aef5e8a9-cf36-4fa7-a223-01df451f5c7b"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"rec", "=", RowBox[{"GenerateRecurrence", "[", "mySum", "]"}]}]], "Input", CellChangeTimes->{{3.831624203440839*^9, 3.8316242187834053`*^9}}, CellLabel->"In[38]:=",ExpressionUUID->"82e79f20-00e6-47d1-8ba2-576f3eb0ecca"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{"-", "8"}], " ", RowBox[{"(", RowBox[{"1", "+", "n"}], ")"}], " ", RowBox[{"(", RowBox[{"3", "+", "n"}], ")"}], " ", RowBox[{"SUM", "[", "n", "]"}]}], "+", RowBox[{"4", " ", RowBox[{"(", RowBox[{"29", "+", RowBox[{"25", " ", "n"}], "+", RowBox[{"5", " ", SuperscriptBox["n", "2"]}]}], ")"}], " ", RowBox[{"SUM", "[", RowBox[{"1", "+", "n"}], "]"}]}], "-", RowBox[{"2", " ", RowBox[{"(", RowBox[{"8", "+", RowBox[{"3", " ", "n"}]}], ")"}], " ", RowBox[{"(", RowBox[{"10", "+", RowBox[{"3", " ", "n"}]}], ")"}], " ", RowBox[{"SUM", "[", RowBox[{"2", "+", "n"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"86", "+", RowBox[{"49", " ", "n"}], "+", RowBox[{"7", " ", SuperscriptBox["n", "2"]}]}], ")"}], " ", RowBox[{"SUM", "[", RowBox[{"3", "+", "n"}], "]"}]}], "-", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"4", "+", "n"}], ")"}], "2"], " ", RowBox[{"SUM", "[", RowBox[{"4", "+", "n"}], "]"}]}]}], "\[Equal]", RowBox[{"-", "1"}]}], "}"}]], "Output", CellChangeTimes->{{3.8316953839386806`*^9, 3.8316954132960663`*^9}}, CellLabel->"Out[38]=",ExpressionUUID->"0347e85a-4b7d-41d7-a55d-23c8023d20da"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"recSol", "=", RowBox[{"SolveRecurrence", "[", RowBox[{ RowBox[{"rec", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"SUM", "[", "n", "]"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.831624203440839*^9, 3.8316242194162197`*^9}}, CellLabel->"In[39]:=",ExpressionUUID->"58ed9239-ce3a-48ae-b128-4c828c74c378"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", SuperscriptBox["2", UnderscriptBox["n", "."]]}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", RowBox[{ SuperscriptBox["2", UnderscriptBox["n", "."]], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", RowBox[{ SubscriptBox["\[Iota]", "1"], " ", SuperscriptBox["2", UnderscriptBox[ SubscriptBox["\[Iota]", "1"], "."]]}]]}], ")"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", RowBox[{ SuperscriptBox["2", UnderscriptBox["n", "."]], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubscriptBox["\[Iota]", "1"]]}], ")"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", RowBox[{ SuperscriptBox["2", UnderscriptBox["n", "."]], " ", RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubsuperscriptBox["\[Iota]", "1", "2"]]}], ")"}]}], "-", SuperscriptBox[ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubscriptBox["\[Iota]", "1"]]}], ")"}], "2"], "+", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox[ SuperscriptBox["2", UnderscriptBox[ SubscriptBox["\[Iota]", "2"], "."]], SubscriptBox["\[Iota]", "2"]]}], RowBox[{ SubscriptBox["\[Iota]", "1"], " ", SuperscriptBox["2", UnderscriptBox[ SubscriptBox["\[Iota]", "1"], "."]]}]]}]}], ")"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", RowBox[{ SuperscriptBox["2", UnderscriptBox["n", "."]], " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", RowBox[{ SubsuperscriptBox["\[Iota]", "1", "2"], " ", SuperscriptBox["2", UnderscriptBox[ SubscriptBox["\[Iota]", "1"], "."]]}]]}], ")"}]}], "-", RowBox[{"2", " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubscriptBox["\[Iota]", "1"]]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", RowBox[{ SubscriptBox["\[Iota]", "1"], " ", SuperscriptBox["2", UnderscriptBox[ SubscriptBox["\[Iota]", "1"], "."]]}]]}], ")"}]}], "+", RowBox[{"3", " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubscriptBox["\[Iota]", "2"]]}], RowBox[{ SubscriptBox["\[Iota]", "1"], " ", SuperscriptBox["2", UnderscriptBox[ SubscriptBox["\[Iota]", "1"], "."]]}]]}], ")"}]}]}], ")"}]}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.831695384055179*^9, 3.8316954172457447`*^9}, CellLabel->"Out[39]=",ExpressionUUID->"a9895de0-301e-43ee-b238-d618bf8da42d"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"FindLinearCombination", "[", RowBox[{"recSol", ",", "mySum", ",", "n", ",", "4"}], "]"}]], "Input", CellChangeTimes->{{3.831624203440839*^9, 3.8316242194162197`*^9}, { 3.831695425235219*^9, 3.831695426001748*^9}}, CellLabel->"In[41]:=",ExpressionUUID->"d77910d2-4451-4f0f-9711-b71b675ffa0c"], Cell[BoxData[ RowBox[{ RowBox[{ SuperscriptBox["2", UnderscriptBox["n", "."]], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubsuperscriptBox["\[Iota]", "1", "2"]]}], ")"}]}], "+", RowBox[{ SuperscriptBox["2", UnderscriptBox["n", "."]], " ", SuperscriptBox[ RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubscriptBox["\[Iota]", "1"]]}], ")"}], "2"]}], "-", RowBox[{"2", " ", SuperscriptBox["2", UnderscriptBox["n", "."]], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", RowBox[{ SubsuperscriptBox["\[Iota]", "1", "2"], " ", SuperscriptBox["2", UnderscriptBox[ SubscriptBox["\[Iota]", "1"], "."]]}]]}], ")"}]}], "-", RowBox[{"2", " ", SuperscriptBox["2", UnderscriptBox["n", "."]], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubscriptBox["\[Iota]", "1"]]}], ")"}], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", RowBox[{ SubscriptBox["\[Iota]", "1"], " ", SuperscriptBox["2", UnderscriptBox[ SubscriptBox["\[Iota]", "1"], "."]]}]]}], ")"}]}], "+", RowBox[{"3", " ", SuperscriptBox["2", UnderscriptBox["n", "."]], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox["1", SubscriptBox["\[Iota]", "2"]]}], RowBox[{ SubscriptBox["\[Iota]", "1"], " ", SuperscriptBox["2", UnderscriptBox[ SubscriptBox["\[Iota]", "1"], "."]]}]]}], ")"}]}], "-", RowBox[{ SuperscriptBox["2", UnderscriptBox["n", "."]], " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "2"], "=", "1"}], SubscriptBox["\[Iota]", "1"]], FractionBox[ SuperscriptBox["2", UnderscriptBox[ SubscriptBox["\[Iota]", "2"], "."]], SubscriptBox["\[Iota]", "2"]]}], RowBox[{ SubscriptBox["\[Iota]", "1"], " ", SuperscriptBox["2", UnderscriptBox[ SubscriptBox["\[Iota]", "1"], "."]]}]]}], ")"}]}]}]], "Output", CellChangeTimes->{ 3.8316953841067057`*^9, {3.831695417413629*^9, 3.831695426566051*^9}}, CellLabel->"Out[41]=",ExpressionUUID->"591abeb5-30a3-413c-8fda-378079103556"] }, Open ]] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["The full machinery", "Section", CellChangeTimes->{{3.83162225869281*^9, 3.831622284289261*^9}, { 3.831624308295051*^9, 3.831624314736467*^9}},ExpressionUUID->"a71d7b28-b50d-4f6f-92ec-\ 0af2b04673e7"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"<<", "EvaluateMultiSums.m"}]], "Input", CellChangeTimes->{{3.831622591328903*^9, 3.831622598042297*^9}}, CellLabel->"In[13]:=",ExpressionUUID->"6b497308-8917-453b-a098-027cb4be0136"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{ "EvaluateMultiSums", " ", "by", " ", "Carsten", " ", "Schneider", " ", RowBox[{"(", RowBox[{ RowBox[{"still", " ", "under", " ", "development"}], "..."}], ")"}]}], " ", "\[LongDash]", " ", RowBox[{"\[Copyright]", " ", "RISC"}], " ", "\[LongDash]", " ", RowBox[{"V", " ", "0.993", " ", RowBox[{"(", RowBox[{ RowBox[{"September", " ", "18"}], ",", " ", "2020"}], ")"}]}]}]], "Print",\ CellFrame->0.5, CellChangeTimes->{3.8316371618108807`*^9}, FontColor->RGBColor[0, 0, 0], Background->RGBColor[0.796887, 0.789075, 0.871107], ButtonBoxOptions->{ Active->True},ExpressionUUID->"e324f51b-b3aa-4ddc-b35e-c0953a78dc03"], Cell[BoxData["\<\"A note for the usage of EvaluateMultiSums: If more \ complicated nested sums arise, one should also activate the package \ HarmonicSums by Jakob Ablinger (for speed ups and appropriate limit \ calculations for infinite sums).\"\>"], "Print", CellChangeTimes->{3.831637161965884*^9}, CellLabel-> "During evaluation of \ In[13]:=",ExpressionUUID->"94976cb0-0963-417a-bb79-8248d2cbb95d"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["The warm-up example", "Subsubsection", CellChangeTimes->{{3.831624321622698*^9, 3.831624327922365*^9}},ExpressionUUID->"0b5cfb3b-1194-4072-b005-\ 1c316ad5cfa4"], Cell[BoxData[ RowBox[{ RowBox[{"summand", "=", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"j", SubscriptBox["!", "."]}], " ", RowBox[{"k", SubscriptBox["!", "."]}], " ", RowBox[{"N", SubscriptBox["!", "."]}], " ", RowBox[{ RowBox[{"(", RowBox[{"j", "+", "k", "+", "N"}], ")"}], SubscriptBox["!", "."]}], " ", RowBox[{"(", RowBox[{ FractionBox[ RowBox[{"2", "+", RowBox[{"2", " ", "j"}], "+", "k", "+", "N"}], RowBox[{ RowBox[{"(", RowBox[{"1", "+", "j", "+", "k"}], ")"}], " ", RowBox[{"(", RowBox[{"1", "+", "j", "+", "N"}], ")"}]}]], "-", RowBox[{"S", "[", RowBox[{"1", ",", "j"}], "]"}], "+", RowBox[{"S", "[", RowBox[{"1", ",", RowBox[{"j", "+", "k"}]}], "]"}], "+", RowBox[{"S", "[", RowBox[{"1", ",", RowBox[{"j", "+", "N"}]}], "]"}], "-", RowBox[{"S", "[", RowBox[{"1", ",", RowBox[{"j", "+", "k", "+", "N"}]}], "]"}]}], ")"}]}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"1", "+", "j", "+", "k"}], ")"}], SubscriptBox["!", "."]}], " ", RowBox[{ RowBox[{"(", RowBox[{"1", "+", "j", "+", "N"}], ")"}], SubscriptBox["!", "."]}], " ", RowBox[{ RowBox[{"(", RowBox[{"1", "+", "k", "+", "N"}], ")"}], SubscriptBox["!", "."]}]}], ")"}]}]}], ";"}]], "Input", CellChangeTimes->{{3.695358996200685*^9, 3.695358999573455*^9}}, CellLabel->"In[14]:=",ExpressionUUID->"93f01bcf-cfc4-4b2e-8154-edbf65c117ae"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"mySum", "=", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"k", "=", "0"}], "\[Infinity]"], RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"j", "=", "0"}], "\[Infinity]"], "summand"}]}]}]], "Input", CellChangeTimes->{{3.695359054949247*^9, 3.695359058352847*^9}}, CellLabel->"In[15]:=",ExpressionUUID->"b5ec48f5-fa1f-4495-9c52-ff1ae788211e"], Cell[BoxData[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"k", "=", "0"}], "\[Infinity]"], RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{"j", "=", "0"}], "\[Infinity]"], RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"j", SubscriptBox["!", "."]}], " ", RowBox[{"k", SubscriptBox["!", "."]}], " ", RowBox[{"N", SubscriptBox["!", "."]}], " ", RowBox[{ RowBox[{"(", RowBox[{"j", "+", "k", "+", "N"}], ")"}], SubscriptBox["!", "."]}], " ", RowBox[{"(", RowBox[{ FractionBox[ RowBox[{"2", "+", RowBox[{"2", " ", "j"}], "+", "k", "+", "N"}], RowBox[{ RowBox[{"(", RowBox[{"1", "+", "j", "+", "k"}], ")"}], " ", RowBox[{"(", RowBox[{"1", "+", "j", "+", "N"}], ")"}]}]], "-", RowBox[{"S", "[", RowBox[{"1", ",", "j"}], "]"}], "+", RowBox[{"S", "[", RowBox[{"1", ",", RowBox[{"j", "+", "k"}]}], "]"}], "+", RowBox[{"S", "[", RowBox[{"1", ",", RowBox[{"j", "+", "N"}]}], "]"}], "-", RowBox[{"S", "[", RowBox[{"1", ",", RowBox[{"j", "+", "k", "+", "N"}]}], "]"}]}], ")"}]}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"1", "+", "j", "+", "k"}], ")"}], SubscriptBox["!", "."]}], " ", RowBox[{ RowBox[{"(", RowBox[{"1", "+", "j", "+", "N"}], ")"}], SubscriptBox["!", "."]}], " ", RowBox[{ RowBox[{"(", RowBox[{"1", "+", "k", "+", "N"}], ")"}], SubscriptBox["!", "."]}]}], ")"}]}]}]}]], "Output", CellChangeTimes->{3.8316371689161167`*^9}, CellLabel->"Out[15]=",ExpressionUUID->"8ed28314-4cf2-4294-9b1d-f13506fab261"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"EvaluateMultiSum", "[", RowBox[{"mySum", ",", RowBox[{"{", "}"}], ",", RowBox[{"{", "N", "}"}], ",", RowBox[{"{", "1", "}"}], ",", RowBox[{"{", "Infinity", "}"}]}], "]"}]], "Input", CellChangeTimes->{{3.695359067258122*^9, 3.695359084209755*^9}}, CellLabel->"In[16]:=",ExpressionUUID->"f1f079bf-5668-4a45-968b-62c5e0aef345"], Cell[BoxData[ RowBox[{ FractionBox[ RowBox[{"3", " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "N"], FractionBox["1", SubsuperscriptBox["\[Iota]", "1", "2"]]}], ")"}]}], RowBox[{"2", " ", "N", " ", RowBox[{"(", RowBox[{"1", "+", "N"}], ")"}]}]], "+", FractionBox[ SuperscriptBox[ RowBox[{"S", "[", RowBox[{"1", ",", "N"}], "]"}], "2"], RowBox[{"2", " ", "N", " ", RowBox[{"(", RowBox[{"1", "+", "N"}], ")"}]}]]}]], "Output", CellChangeTimes->{3.83163721235248*^9}, CellLabel->"Out[16]=",ExpressionUUID->"a3463eac-322b-4f37-9eae-9f8e2910c8f7"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["The triple sum", "Subsection", CellChangeTimes->{{3.8316229865058203`*^9, 3.831622988671309*^9}},ExpressionUUID->"3c14e639-a713-486d-b2dc-\ 63d98f3625a4"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"summand", "=", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"-", "1"}], ")"}], "^", "r"}], " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "j"}], "+", "n", "-", "2"}], ")"}], "!"}], " ", RowBox[{"Binomial", "[", RowBox[{ RowBox[{"j", "+", "1"}], ",", "r"}], "]"}], " ", RowBox[{ RowBox[{"r", "!"}], "/", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "j"}], "+", "n", "+", "r"}], ")"}], "!"}]}], " ", RowBox[{ RowBox[{"(", RowBox[{"-", "1"}], ")"}], "^", "s"}], " ", RowBox[{ RowBox[{"Binomial", "[", RowBox[{ RowBox[{ RowBox[{"-", "j"}], "+", "n", "+", "r", "-", "2"}], ",", "s"}], "]"}], "/", RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{"n", "-", "s"}], ")"}], " ", RowBox[{"(", RowBox[{"s", "+", "1"}], ")"}]}], ")"}]}]}]}]], "Input", CellChangeTimes->{{3.831622989459194*^9, 3.831622999190834*^9}}, CellLabel->"In[17]:=",ExpressionUUID->"39198eab-81af-4ded-b667-adc8ce2ef7ae"], Cell[BoxData[ FractionBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], RowBox[{"r", "+", "s"}]], " ", RowBox[{"Binomial", "[", RowBox[{ RowBox[{"1", "+", "j"}], ",", "r"}], "]"}], " ", RowBox[{"Binomial", "[", RowBox[{ RowBox[{ RowBox[{"-", "2"}], "-", "j", "+", "n", "+", "r"}], ",", "s"}], "]"}], " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "2"}], "-", "j", "+", "n"}], ")"}], "!"}], " ", RowBox[{"r", "!"}]}], RowBox[{ RowBox[{"(", RowBox[{"n", "-", "s"}], ")"}], " ", RowBox[{"(", RowBox[{"1", "+", "s"}], ")"}], " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "j"}], "+", "n", "+", "r"}], ")"}], "!"}]}]]], "Output", CellChangeTimes->{3.8316372305046787`*^9}, CellLabel->"Out[17]=",ExpressionUUID->"563d48ee-474f-465a-91ef-dd132bc74270"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"res", "=", RowBox[{ RowBox[{"EvaluateMultiSum", "[", RowBox[{"summand", ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"s", ",", "0", ",", RowBox[{ RowBox[{"-", "j"}], "+", "n", "+", "r", "-", "2"}]}], "}"}], ",", RowBox[{"{", RowBox[{"r", ",", "0", ",", RowBox[{"j", "+", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "0", ",", RowBox[{"n", "-", "2"}]}], "}"}]}], "}"}], ",", RowBox[{"{", "n", "}"}], ",", RowBox[{"{", "2", "}"}]}], "]"}], "//", "Timing"}]}]], "Input", CellChangeTimes->{{3.831622989459194*^9, 3.831623007427988*^9}}, CellLabel->"In[18]:=",ExpressionUUID->"446f258d-7c61-4e12-9dcf-f2a7e1a6f651"], Cell[BoxData[ RowBox[{"{", RowBox[{"75.601439`", ",", RowBox[{ FractionBox[ RowBox[{ RowBox[{"-", "1"}], "-", "n", "-", SuperscriptBox["n", "2"]}], RowBox[{ SuperscriptBox["n", "2"], " ", SuperscriptBox[ RowBox[{"(", RowBox[{"1", "+", "n"}], ")"}], "3"]}]], "+", FractionBox[ RowBox[{ RowBox[{"(", RowBox[{"1", "+", "n", "+", SuperscriptBox["n", "2"]}], ")"}], " ", SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], UnderscriptBox["n", "."]]}], RowBox[{ SuperscriptBox["n", "2"], " ", SuperscriptBox[ RowBox[{"(", RowBox[{"1", "+", "n"}], ")"}], "3"]}]], "-", FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubsuperscriptBox["\[Iota]", "1", "2"]]}], RowBox[{"1", "+", "n"}]], "+", FractionBox[ RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox["1", SubscriptBox["\[Iota]", "1"]]}], SuperscriptBox[ RowBox[{"(", RowBox[{"1", "+", "n"}], ")"}], "2"]], "-", FractionBox[ RowBox[{"2", " ", RowBox[{"(", RowBox[{ UnderoverscriptBox[ SuperscriptBox["\[Sum]", "."], RowBox[{ SubscriptBox["\[Iota]", "1"], "=", "1"}], "n"], FractionBox[ SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], UnderscriptBox[ SubscriptBox["\[Iota]", "1"], "."]], SubsuperscriptBox["\[Iota]", "1", "2"]]}], ")"}]}], RowBox[{"1", "+", "n"}]]}]}], "}"}]], "Output", CellChangeTimes->{3.831637308313405*^9}, CellLabel->"Out[18]=",ExpressionUUID->"d7715fd7-59b7-477b-921a-7949271bfd8d"] }, Open ]] }, Open ]] }, Closed]] }, Closed]] }, WindowSize->{1130.218487394958, 589.9159663865546}, WindowMargins->{{0, Automatic}, {0, Automatic}}, TaggingRules->{ "WelcomeScreenSettings" -> {"FEStarting" -> False}, "TryRealOnly" -> False}, FrontEndVersion->"13.0 for Linux x86 (64-bit) (December 2, 2021)", StyleDefinitions->"Default.nb", ExpressionUUID->"e559c8c7-8399-4dcc-beaf-1ef314d762f0" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[558, 20, 197, 3, 29, "Input",ExpressionUUID->"851884f2-edd5-4c2d-9f36-5bd9b7a85eec"], Cell[CellGroupData[{ Cell[780, 27, 293, 5, 29, "Input",ExpressionUUID->"d1b6b1e4-eb24-487b-9e19-61e898f3794f"], Cell[1076, 34, 954, 25, 40, "Print",ExpressionUUID->"8bc9d068-7ad0-4b56-b962-d121c81db79f"] }, Open ]], Cell[CellGroupData[{ Cell[2067, 64, 310, 5, 160, "Title",ExpressionUUID->"704c8277-85dd-4b00-be79-882101fd2c60"], Cell[CellGroupData[{ Cell[2402, 73, 322, 6, 29, "Input",ExpressionUUID->"0c38568f-6427-47f1-a08e-fcd0599ea15e"], Cell[2727, 81, 990, 26, 40, "Print",ExpressionUUID->"be8a0400-d84d-472f-aa39-76a649e61842"] }, Open ]], Cell[CellGroupData[{ Cell[3754, 112, 166, 3, 55, "Subsection",ExpressionUUID->"c0c9b550-ef2e-4668-b871-a8e57b157f8e"], Cell[CellGroupData[{ Cell[3945, 119, 391, 9, 29, "Input",ExpressionUUID->"236e52e3-cd97-46b7-9ebd-f7fbe45f093d"], Cell[4339, 130, 342, 8, 58, "Output",ExpressionUUID->"5d211ac4-6321-4a72-8730-84cebd5aafa0"] }, Open ]], Cell[CellGroupData[{ Cell[4718, 143, 216, 3, 29, "Input",ExpressionUUID->"ba8cecd1-67dd-456c-a79e-da3b767c5f7f"], Cell[4937, 148, 403, 10, 59, "Output",ExpressionUUID->"81cd4355-5151-4689-9e20-f0fd57ec22a6"] }, Open ]], Cell[CellGroupData[{ Cell[5377, 163, 753, 20, 29, "Input",ExpressionUUID->"6353f64f-9f87-4d14-95a7-85380780b73e"], Cell[6133, 185, 636, 18, 57, "Output",ExpressionUUID->"7d693058-441d-4945-a55c-8d76ba893142"] }, Open ]], Cell[CellGroupData[{ Cell[6806, 208, 214, 3, 29, "Input",ExpressionUUID->"a519d9e9-d59b-4115-b8df-760330976bd9"], Cell[7023, 213, 1454, 44, 61, "Output",ExpressionUUID->"4f62fa47-a6e0-4fb1-91b0-671be22a7419"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[8526, 263, 171, 3, 55, "Subsection",ExpressionUUID->"5476d50b-8a19-49b0-b076-b62a8bb09a08"], Cell[CellGroupData[{ Cell[8722, 270, 457, 11, 29, "Input",ExpressionUUID->"78db9115-2c33-47e1-b9f4-34989752c2ec"], Cell[9182, 283, 572, 15, 61, "Output",ExpressionUUID->"a99b1899-fec8-4c7b-9617-62a02b14af7c"] }, Open ]], Cell[CellGroupData[{ Cell[9791, 303, 286, 5, 29, "Input",ExpressionUUID->"1bd3158d-e04a-4c6d-97a4-5190ba4b1517"], Cell[10080, 310, 612, 17, 61, "Output",ExpressionUUID->"59e082d6-8be1-4556-80fd-5bcdcaec6e4b"] }, Open ]], Cell[CellGroupData[{ Cell[10729, 332, 642, 16, 29, "Input",ExpressionUUID->"024523ba-9a53-46ef-9fc2-368a2cea21a4"], Cell[11374, 350, 644, 20, 58, "Output",ExpressionUUID->"c1532054-d9ea-49ca-bffb-11a7534c1566"] }, Open ]], Cell[CellGroupData[{ Cell[12055, 375, 285, 5, 29, "Input",ExpressionUUID->"4f783cb6-5582-445f-9b9e-922cb4ad4ecf"], Cell[12343, 382, 1036, 35, 66, "Output",ExpressionUUID->"1718b7eb-4ec0-4b67-8b7f-1c8b34d77660"] }, Open ]], Cell[CellGroupData[{ Cell[13416, 422, 406, 10, 29, "Input",ExpressionUUID->"e4681773-75c4-4916-8f62-2ed964138f80"], Cell[13825, 434, 954, 33, 66, "Output",ExpressionUUID->"a2b6656e-2628-48a2-8d67-20d972b3c1e8"] }, Open ]], Cell[CellGroupData[{ Cell[14816, 472, 320, 7, 29, "Input",ExpressionUUID->"9954db1d-25e9-43b4-b08d-ffbed6373fda"], Cell[15139, 481, 479, 14, 48, "Output",ExpressionUUID->"f1be8f99-3b30-4819-9d82-59d5a8d0de29"] }, Open ]], Cell[CellGroupData[{ Cell[15655, 500, 651, 16, 29, "Input",ExpressionUUID->"dda35b47-bcb8-491d-b219-93075c117577"], Cell[16309, 518, 564, 18, 60, "Output",ExpressionUUID->"e79a5fe9-1243-43db-80ab-88dfde90c07a"] }, Open ]], Cell[CellGroupData[{ Cell[16910, 541, 238, 4, 29, "Input",ExpressionUUID->"9077afb3-d77f-49da-a29e-d16f411837e6"], Cell[17151, 547, 1656, 59, 62, "Output",ExpressionUUID->"e347dfcc-c9a2-45a8-8982-cd0eae13acce"] }, Open ]] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[18868, 613, 265, 4, 72, "Title",ExpressionUUID->"113b3435-56cb-42f5-aecf-7be6a9eab464"], Cell[CellGroupData[{ Cell[19158, 621, 267, 4, 55, "Subsection",ExpressionUUID->"f91662f3-4ffb-4cbd-a924-61fa36336114"], Cell[CellGroupData[{ Cell[19450, 629, 527, 13, 29, "Input",ExpressionUUID->"66c7f248-105c-409f-bdb1-b6cf4f480490"], Cell[19980, 644, 433, 13, 59, "Output",ExpressionUUID->"123dcdcf-9dbd-4aba-ae36-f9dad4ae9112"] }, Open ]], Cell[CellGroupData[{ Cell[20450, 662, 383, 9, 50, "Input",ExpressionUUID->"cdf2189f-2fab-401b-b60e-9bf6d818c2a8"], Cell[20836, 673, 292, 8, 58, "Output",ExpressionUUID->"c5e36c55-d128-45d2-a655-77c3b55125c5"] }, Open ]], Cell[CellGroupData[{ Cell[21165, 686, 215, 3, 29, "Input",ExpressionUUID->"42a9e20e-5ee9-4507-b7e1-b963b309eaf5"], Cell[21383, 691, 430, 12, 61, "Output",ExpressionUUID->"0472079d-4317-4721-9d1f-52ddae94e364"] }, Open ]], Cell[CellGroupData[{ Cell[21850, 708, 226, 4, 29, "Input",ExpressionUUID->"8b877440-fdfa-4e1b-ad06-8d412e7e8cc9"], Cell[22079, 714, 150, 2, 33, "Output",ExpressionUUID->"b9aca63a-4c27-4fda-955c-30f9d50ba688"] }, Open ]], Cell[CellGroupData[{ Cell[22266, 721, 775, 22, 31, "Input",ExpressionUUID->"81fbce5f-089e-41d5-9823-956c744a20d2"], Cell[23044, 745, 690, 22, 75, "Output",ExpressionUUID->"f1dfb6bc-862c-49ba-8309-4b4f463abe1b"] }, Open ]], Cell[CellGroupData[{ Cell[23771, 772, 228, 4, 29, "Input",ExpressionUUID->"65536035-af1f-4c40-bdd4-16bb1fbe157a"], Cell[24002, 778, 175, 3, 48, "Output",ExpressionUUID->"4e848326-e9b8-472d-9665-0a38fe148ec3"] }, Open ]], Cell[CellGroupData[{ Cell[24214, 786, 269, 6, 29, "Input",ExpressionUUID->"6ef87d02-6ba2-4e65-b5ff-680037702a17"], Cell[24486, 794, 1164, 38, 90, "Output",ExpressionUUID->"97777b37-2ebd-4a3b-a9eb-47627b686e9a"] }, Open ]], Cell[CellGroupData[{ Cell[25687, 837, 223, 4, 29, "Input",ExpressionUUID->"75005778-1b7e-4f08-8a99-e8b09e75f474"], Cell[25913, 843, 176, 3, 49, "Output",ExpressionUUID->"a8d52c27-ab66-46c8-b67b-30704d278899"] }, Open ]], Cell[CellGroupData[{ Cell[26126, 851, 235, 4, 29, "Input",ExpressionUUID->"51c5aabc-c553-471c-a146-99226db20810"], Cell[26364, 857, 1217, 40, 90, "Output",ExpressionUUID->"4b9b0f33-f17d-487d-bd07-e71ea0beee7f"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[27630, 903, 171, 3, 39, "Subsection",ExpressionUUID->"1d2860c1-3e4a-40bc-855a-7423b34aa01a"], Cell[CellGroupData[{ Cell[27826, 910, 773, 22, 31, "Input",ExpressionUUID->"425a78a9-cb98-45fb-a6e6-78a30511735c"], Cell[28602, 934, 691, 22, 75, "Output",ExpressionUUID->"368a0dfe-252e-4ae6-a096-0c768fdf5645"] }, Open ]], Cell[29308, 959, 3404, 111, 90, "Input",ExpressionUUID->"58e7a766-98ca-4a0a-989f-688b3536f48d"], Cell[CellGroupData[{ Cell[32737, 1074, 274, 6, 29, "Input",ExpressionUUID->"bfe104f0-6b23-4c4c-99d3-3ace36728362"], Cell[33014, 1082, 254, 6, 48, "Output",ExpressionUUID->"9900d1b6-8b38-404b-a6d4-8e6d6fa2a77b"] }, Open ]], Cell[CellGroupData[{ Cell[33305, 1093, 411, 10, 29, "Input",ExpressionUUID->"7a5fd869-602b-493a-bccb-3fa172ee7ac3"], Cell[33719, 1105, 1870, 58, 79, "Output",ExpressionUUID->"16d49442-504a-470a-9b88-d2ac845b79c7"] }, Open ]], Cell[CellGroupData[{ Cell[35626, 1168, 442, 11, 29, "Input",ExpressionUUID->"400fea04-d5f6-42ef-b976-6903a446c250"], Cell[36071, 1181, 1446, 45, 79, "Output",ExpressionUUID->"07b473cc-8efb-491f-a597-9b22396b4176"] }, Open ]], Cell[CellGroupData[{ Cell[37554, 1231, 411, 10, 29, "Input",ExpressionUUID->"d3d91d94-9cef-4cbc-9ad1-5f13021da5c1"], Cell[37968, 1243, 6156, 192, 154, "Output",ExpressionUUID->"76b24814-2a55-4261-bf20-3ff1e939738e"] }, Open ]], Cell[CellGroupData[{ Cell[44161, 1440, 374, 10, 29, "Input",ExpressionUUID->"c5ac63ec-7e18-4bd9-b5de-e72811199250"], Cell[44538, 1452, 2574, 81, 79, "Output",ExpressionUUID->"2dc6fcd8-d9de-46e9-ab43-9ccfee7cd2c2"] }, Open ]], Cell[CellGroupData[{ Cell[47149, 1538, 258, 5, 29, "Input",ExpressionUUID->"652ed41e-1888-41fe-92f9-a7227b0570e3"], Cell[47410, 1545, 149, 2, 33, "Output",ExpressionUUID->"b244ca41-0c24-4f6b-b82e-759db9ff15ea"] }, Open ]], Cell[CellGroupData[{ Cell[47596, 1552, 293, 6, 29, "Input",ExpressionUUID->"9222af92-c3f1-4086-bccb-ac14201c6226"], Cell[47892, 1560, 3688, 124, 92, "Output",ExpressionUUID->"228cbb61-83a6-40b9-84a5-692821224408"] }, Open ]], Cell[CellGroupData[{ Cell[51617, 1689, 417, 7, 29, "Input",ExpressionUUID->"d21be423-8a6b-42a8-9a21-fe498cc0479f"], Cell[52037, 1698, 12399, 388, 244, "Output",ExpressionUUID->"1fe56d23-7a38-4012-983a-ecd4232ad48e"] }, Open ]], Cell[CellGroupData[{ Cell[64473, 2091, 430, 10, 29, "Input",ExpressionUUID->"4bbb293e-5b5d-4db1-bd39-ee5e7a73e766"], Cell[64906, 2103, 19192, 545, 607, "Output",ExpressionUUID->"af68bb52-2b2a-4f0a-b353-181325317753"] }, Open ]] }, Closed]] }, Closed]], Cell[CellGroupData[{ Cell[84159, 2655, 318, 5, 72, "Title",ExpressionUUID->"da6e2fb4-83d1-4243-8d7b-f595ca8fd4ba"], Cell[CellGroupData[{ Cell[84502, 2664, 157, 3, 67, "Section",ExpressionUUID->"3ec7fa06-35d9-4269-aa2e-f49b20fefaa1"], Cell[CellGroupData[{ Cell[84684, 2671, 158, 3, 55, "Subsection",ExpressionUUID->"6d9045c2-d4c1-4b53-8e37-0ff603bd2631"], Cell[CellGroupData[{ Cell[84867, 2678, 376, 9, 29, "Input",ExpressionUUID->"e2c8ef8e-002a-4f41-ab1b-0ba2d5a47ad7"], Cell[85246, 2689, 327, 8, 58, "Output",ExpressionUUID->"e79340c8-ffcf-4267-b19c-88f4c47a4253"] }, Open ]], Cell[CellGroupData[{ Cell[85610, 2702, 213, 3, 29, "Input",ExpressionUUID->"2db5acb2-23dc-4674-9bbe-d5df57fbfd67"], Cell[85826, 2707, 337, 9, 33, "Output",ExpressionUUID->"62df2090-dead-4cc6-b9a4-ef6c6e4f9e78"] }, Open ]], Cell[CellGroupData[{ Cell[86200, 2721, 387, 9, 29, "Input",ExpressionUUID->"bc1178b4-b258-40cf-9be3-6bd7bf5f4bc2"], Cell[86590, 2732, 341, 9, 60, "Output",ExpressionUUID->"f7e2dcfb-f4f8-4bee-9c21-7a6dce0678d4"] }, Open ]], Cell[CellGroupData[{ Cell[86968, 2746, 215, 3, 29, "Input",ExpressionUUID->"d88ff7ab-bb59-49a8-ba05-b944dce7359e"], Cell[87186, 2751, 464, 15, 33, "Output",ExpressionUUID->"10f60be6-a0ba-462b-a201-e4fe4a7ae019"] }, Open ]], Cell[CellGroupData[{ Cell[87687, 2771, 555, 14, 29, "Input",ExpressionUUID->"38e23792-f224-4bba-b914-3a3f637e23f4"], Cell[88245, 2787, 520, 17, 60, "Output",ExpressionUUID->"f7db102d-e3dc-4619-adfe-19f700ec29b4"] }, Open ]], Cell[CellGroupData[{ Cell[88802, 2809, 213, 3, 29, "Input",ExpressionUUID->"245256e3-7de1-4f5e-8f29-50c5266baf76"], Cell[89018, 2814, 1869, 70, 67, "Output",ExpressionUUID->"9eca07ce-0c18-4801-b190-afcad63a17c7"] }, Open ]], Cell[CellGroupData[{ Cell[90924, 2889, 551, 14, 29, "Input",ExpressionUUID->"af75c852-4958-4d3e-b6b3-691592942265"], Cell[91478, 2905, 540, 17, 60, "Output",ExpressionUUID->"101a6230-bc42-458a-9a69-b079b252321c"] }, Open ]], Cell[CellGroupData[{ Cell[92055, 2927, 240, 4, 29, "Input",ExpressionUUID->"0d224b46-7d22-46ed-9dc4-6daea32b5c33"], Cell[92298, 2933, 1631, 58, 62, "Output",ExpressionUUID->"42faf624-bc57-4fe2-8927-410a089ebe62"] }, Open ]], Cell[CellGroupData[{ Cell[93966, 2996, 386, 9, 29, "Input",ExpressionUUID->"22cd9332-4904-4fc0-abd8-7924473c4ca0"], Cell[94355, 3007, 223, 4, 33, "Output",ExpressionUUID->"a305e6a9-8e53-4e96-bb12-56bfcfd75bc9"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[94627, 3017, 191, 3, 55, "Subsection",ExpressionUUID->"75de7937-3466-4af4-9301-40848fbc0599"], Cell[CellGroupData[{ Cell[94843, 3024, 533, 13, 29, "Input",ExpressionUUID->"0eac2d65-2dfc-445a-bf38-a5bc1a8bb50a"], Cell[95379, 3039, 425, 13, 58, "Output",ExpressionUUID->"2cb7a3ec-2951-4b4d-b5e3-315735cccfc4"] }, Open ]], Cell[CellGroupData[{ Cell[95841, 3057, 213, 3, 29, "Input",ExpressionUUID->"6396d67d-db28-48e3-a9b2-8b8bde677f01"], Cell[96057, 3062, 519, 16, 59, "Output",ExpressionUUID->"aef5e8a9-cf36-4fa7-a223-01df451f5c7b"] }, Open ]], Cell[CellGroupData[{ Cell[96613, 3083, 247, 4, 29, "Input",ExpressionUUID->"82e79f20-00e6-47d1-8ba2-576f3eb0ecca"], Cell[96863, 3089, 1424, 44, 35, "Output",ExpressionUUID->"0347e85a-4b7d-41d7-a55d-23c8023d20da"] }, Open ]], Cell[CellGroupData[{ Cell[98324, 3138, 357, 8, 29, "Input",ExpressionUUID->"58ed9239-ce3a-48ae-b128-4c828c74c378"], Cell[98684, 3148, 5021, 150, 146, "Output",ExpressionUUID->"a9895de0-301e-43ee-b238-d618bf8da42d"] }, Open ]], Cell[CellGroupData[{ Cell[103742, 3303, 322, 5, 29, "Input",ExpressionUUID->"d77910d2-4451-4f0f-9711-b71b675ffa0c"], Cell[104067, 3310, 3493, 114, 89, "Output",ExpressionUUID->"591abeb5-30a3-413c-8fda-378079103556"] }, Open ]] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[107621, 3431, 211, 4, 67, "Section",ExpressionUUID->"a71d7b28-b50d-4f6f-92ec-0af2b04673e7"], Cell[CellGroupData[{ Cell[107857, 3439, 208, 3, 29, "Input",ExpressionUUID->"6b497308-8917-453b-a098-027cb4be0136"], Cell[CellGroupData[{ Cell[108090, 3446, 695, 19, 40, "Print",ExpressionUUID->"e324f51b-b3aa-4ddc-b35e-c0953a78dc03"], Cell[108788, 3467, 406, 7, 43, "Print",ExpressionUUID->"94976cb0-0963-417a-bb79-8248d2cbb95d"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[109243, 3480, 170, 3, 45, "Subsubsection",ExpressionUUID->"0b5cfb3b-1194-4072-b005-1c316ad5cfa4"], Cell[109416, 3485, 1723, 52, 72, "Input",ExpressionUUID->"93f01bcf-cfc4-4b2e-8154-edbf65c117ae"], Cell[CellGroupData[{ Cell[111164, 3541, 444, 11, 55, "Input",ExpressionUUID->"b5ec48f5-fa1f-4495-9c52-ff1ae788211e"], Cell[111611, 3554, 1888, 58, 59, "Output",ExpressionUUID->"8ed28314-4cf2-4294-9b1d-f13506fab261"] }, Open ]], Cell[CellGroupData[{ Cell[113536, 3617, 369, 8, 29, "Input",ExpressionUUID->"f1f079bf-5668-4a45-968b-62c5e0aef345"], Cell[113908, 3627, 734, 23, 72, "Output",ExpressionUUID->"a3463eac-322b-4f37-9eae-9f8e2910c8f7"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[114691, 3656, 164, 3, 55, "Subsection",ExpressionUUID->"3c14e639-a713-486d-b2dc-63d98f3625a4"], Cell[CellGroupData[{ Cell[114880, 3663, 1065, 35, 31, "Input",ExpressionUUID->"39198eab-81af-4ded-b667-adc8ce2ef7ae"], Cell[115948, 3700, 894, 30, 52, "Output",ExpressionUUID->"563d48ee-474f-465a-91ef-dd132bc74270"] }, Open ]], Cell[CellGroupData[{ Cell[116879, 3735, 772, 20, 29, "Input",ExpressionUUID->"446f258d-7c61-4e12-9dcf-f2a7e1a6f651"], Cell[117654, 3757, 1981, 64, 73, "Output",ExpressionUUID->"d7715fd7-59b7-477b-921a-7949271bfd8d"] }, Open ]] }, Open ]] }, Closed]] }, Closed]] } ] *)