(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 12.3' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 158, 7] NotebookDataLength[ 42611, 1158] NotebookOptionsPosition[ 39591, 1098] NotebookOutlinePosition[ 39990, 1114] CellTagsIndexPosition[ 39947, 1111] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Integrate", "[", RowBox[{ FractionBox[ RowBox[{ SuperscriptBox["x", RowBox[{ RowBox[{"-", "7"}], "/", "3"}]], RowBox[{"Exp", "[", RowBox[{ RowBox[{"-", "1"}], "/", SuperscriptBox["x", "2"]}], "]"}]}], RowBox[{"(", RowBox[{"x", "+", "y"}], ")"}]], ",", RowBox[{"{", RowBox[{"x", ",", "0", ",", "\[Infinity]"}], "}"}], ",", RowBox[{"Assumptions", "->", RowBox[{"{", RowBox[{ RowBox[{"x", ">", "0"}], ",", RowBox[{"y", ">", "0"}], ",", RowBox[{"yc", ">", "0"}]}], "}"}]}]}], "]"}]], "Input", CellChangeTimes->{3.8371541989870367`*^9}, CellLabel->"In[48]:=",ExpressionUUID->"fababbb8-381a-4687-8bb1-07b91e8a5e30"], Cell[BoxData[ FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox["1", SuperscriptBox["y", "2"]]}]], " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"ExpIntegralE", "[", RowBox[{ FractionBox["7", "6"], ",", RowBox[{"-", FractionBox["1", SuperscriptBox["y", "2"]]}]}], "]"}]}], " ", RowBox[{"Gamma", "[", FractionBox["1", "6"], "]"}]}], "+", RowBox[{"4", " ", "y", " ", RowBox[{"ExpIntegralE", "[", RowBox[{ FractionBox["5", "3"], ",", RowBox[{"-", FractionBox["1", SuperscriptBox["y", "2"]]}]}], "]"}], " ", RowBox[{"Gamma", "[", FractionBox["2", "3"], "]"}]}]}], ")"}]}], RowBox[{"12", " ", SuperscriptBox["y", "2"]}]]], "Output", CellChangeTimes->{3.837154204903788*^9}, CellLabel->"Out[48]=",ExpressionUUID->"18ce2878-724a-4922-9798-030bcc1e0186"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"i1", "=", RowBox[{ FractionBox[ SuperscriptBox["\[Theta]", "2"], "\[Pi]"], RowBox[{"Integrate", "[", RowBox[{ FractionBox[ RowBox[{ SuperscriptBox["x", RowBox[{ RowBox[{"-", "7"}], "/", "3"}]], RowBox[{"Exp", "[", RowBox[{ RowBox[{"-", "1"}], "/", SuperscriptBox["x", "2"]}], "]"}]}], RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"x", "+", "yc"}], ")"}], "2"], RowBox[{"(", RowBox[{"x", "+", "y"}], ")"}]}]], ",", RowBox[{"{", RowBox[{"x", ",", "0", ",", "\[Infinity]"}], "}"}], ",", RowBox[{"Assumptions", "->", RowBox[{"{", RowBox[{ RowBox[{"x", ">", "0"}], ",", RowBox[{"y", ">", "0"}], ",", RowBox[{"yc", ">", "0"}]}], "}"}]}]}], "]"}]}]}]], "Input", CellChangeTimes->{{3.837152286872525*^9, 3.8371524292104893`*^9}, { 3.837152588581934*^9, 3.8371525930135*^9}, {3.8371526561032867`*^9, 3.8371526576627903`*^9}, {3.837152833506393*^9, 3.837152859458221*^9}, { 3.8371530145983763`*^9, 3.83715307324615*^9}, {3.8371531339360533`*^9, 3.837153158839828*^9}, {3.83715327430684*^9, 3.837153280810217*^9}, { 3.837153503855143*^9, 3.837153504046052*^9}}, CellLabel->"In[37]:=",ExpressionUUID->"bf6da1af-2ff5-46b6-a7be-85f3bdd44d14"], Cell[BoxData[ RowBox[{ FractionBox["1", RowBox[{"36", " ", "\[Pi]", " ", SuperscriptBox[ RowBox[{"(", RowBox[{"y", "-", "yc"}], ")"}], "2"]}]], RowBox[{ SuperscriptBox["\[Theta]", "2"], " ", RowBox[{"(", RowBox[{ RowBox[{ FractionBox["1", SuperscriptBox["y", RowBox[{"7", "/", "3"}]]], RowBox[{"3", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox["1", SuperscriptBox["y", "2"]]}]], " ", RowBox[{"(", RowBox[{ RowBox[{"8", " ", SqrtBox["3"], " ", "\[Pi]"}], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], RowBox[{"1", "/", "6"}]], " ", RowBox[{"Gamma", "[", RowBox[{"-", FractionBox["1", "6"]}], "]"}], " ", RowBox[{"Gamma", "[", FractionBox["1", "6"], "]"}]}], "-", RowBox[{"4", " ", SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], RowBox[{"2", "/", "3"}]], " ", RowBox[{"Gamma", "[", RowBox[{"-", FractionBox["2", "3"]}], "]"}], " ", RowBox[{"Gamma", "[", FractionBox["2", "3"], "]"}]}], "+", RowBox[{"4", " ", SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], RowBox[{"2", "/", "3"}]], " ", RowBox[{"Gamma", "[", FractionBox["2", "3"], "]"}], " ", RowBox[{"Gamma", "[", RowBox[{ RowBox[{"-", FractionBox["2", "3"]}], ",", RowBox[{"-", FractionBox["1", SuperscriptBox["y", "2"]]}]}], "]"}]}], "-", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], RowBox[{"1", "/", "6"}]], " ", RowBox[{"Gamma", "[", FractionBox["1", "6"], "]"}], " ", RowBox[{"Gamma", "[", RowBox[{ RowBox[{"-", FractionBox["1", "6"]}], ",", RowBox[{"-", FractionBox["1", SuperscriptBox["y", "2"]]}]}], "]"}]}]}], ")"}]}]}], "-", RowBox[{ FractionBox["1", SuperscriptBox["yc", RowBox[{"7", "/", "3"}]]], RowBox[{"3", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox["1", SuperscriptBox["yc", "2"]]}]], " ", RowBox[{"(", RowBox[{ RowBox[{"8", " ", SqrtBox["3"], " ", "\[Pi]"}], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], RowBox[{"1", "/", "6"}]], " ", RowBox[{"Gamma", "[", RowBox[{"-", FractionBox["1", "6"]}], "]"}], " ", RowBox[{"Gamma", "[", FractionBox["1", "6"], "]"}]}], "-", RowBox[{"4", " ", SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], RowBox[{"2", "/", "3"}]], " ", RowBox[{"Gamma", "[", RowBox[{"-", FractionBox["2", "3"]}], "]"}], " ", RowBox[{"Gamma", "[", FractionBox["2", "3"], "]"}]}], "+", RowBox[{"4", " ", SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], RowBox[{"2", "/", "3"}]], " ", RowBox[{"Gamma", "[", FractionBox["2", "3"], "]"}], " ", RowBox[{"Gamma", "[", RowBox[{ RowBox[{"-", FractionBox["2", "3"]}], ",", RowBox[{"-", FractionBox["1", SuperscriptBox["yc", "2"]]}]}], "]"}]}], "-", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], RowBox[{"1", "/", "6"}]], " ", RowBox[{"Gamma", "[", FractionBox["1", "6"], "]"}], " ", RowBox[{"Gamma", "[", RowBox[{ RowBox[{"-", FractionBox["1", "6"]}], ",", RowBox[{"-", FractionBox["1", SuperscriptBox["yc", "2"]]}]}], "]"}]}]}], ")"}]}]}], "+", RowBox[{ FractionBox["1", SuperscriptBox["yc", RowBox[{"16", "/", "3"}]]], RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox["1", SuperscriptBox["yc", "2"]]}]], " ", RowBox[{"(", RowBox[{"y", "-", "yc"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "48"}], " ", SqrtBox["3"], " ", "\[Pi]"}], "+", RowBox[{"56", " ", SqrtBox["3"], " ", "\[Pi]", " ", SuperscriptBox["yc", "2"]}], "+", RowBox[{"6", " ", SuperscriptBox["\[ExponentialE]", FractionBox["1", SuperscriptBox["yc", "2"]]], " ", SuperscriptBox["yc", RowBox[{"7", "/", "3"}]], " ", RowBox[{"Gamma", "[", FractionBox["1", "6"], "]"}]}], "-", RowBox[{"6", " ", SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], RowBox[{"1", "/", "6"}]], " ", RowBox[{"Gamma", "[", RowBox[{"-", FractionBox["1", "6"]}], "]"}], " ", RowBox[{"Gamma", "[", FractionBox["1", "6"], "]"}]}], "+", RowBox[{"7", " ", SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], RowBox[{"1", "/", "6"}]], " ", SuperscriptBox["yc", "2"], " ", RowBox[{"Gamma", "[", RowBox[{"-", FractionBox["1", "6"]}], "]"}], " ", RowBox[{"Gamma", "[", FractionBox["1", "6"], "]"}]}], "-", RowBox[{"36", " ", SuperscriptBox["\[ExponentialE]", FractionBox["1", SuperscriptBox["yc", "2"]]], " ", SuperscriptBox["yc", RowBox[{"4", "/", "3"}]], " ", RowBox[{"Gamma", "[", FractionBox["2", "3"], "]"}]}], "-", RowBox[{"12", " ", SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], RowBox[{"2", "/", "3"}]], " ", SuperscriptBox["yc", "2"], " ", RowBox[{"Gamma", "[", RowBox[{"-", FractionBox["2", "3"]}], "]"}], " ", RowBox[{"Gamma", "[", FractionBox["2", "3"], "]"}]}], "-", RowBox[{"36", " ", SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], RowBox[{"2", "/", "3"}]], " ", RowBox[{"Gamma", "[", FractionBox["1", "3"], "]"}], " ", RowBox[{"Gamma", "[", FractionBox["2", "3"], "]"}]}], "+", RowBox[{"24", " ", SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], RowBox[{"2", "/", "3"}]], " ", SuperscriptBox["yc", "2"], " ", RowBox[{"Gamma", "[", FractionBox["1", "3"], "]"}], " ", RowBox[{"Gamma", "[", FractionBox["2", "3"], "]"}]}], "+", RowBox[{"12", " ", SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], RowBox[{"2", "/", "3"}]], " ", SuperscriptBox["yc", "2"], " ", RowBox[{"Gamma", "[", FractionBox["2", "3"], "]"}], " ", RowBox[{"Gamma", "[", RowBox[{ RowBox[{"-", FractionBox["2", "3"]}], ",", RowBox[{"-", FractionBox["1", SuperscriptBox["yc", "2"]]}]}], "]"}]}], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], RowBox[{"1", "/", "6"}]], " ", RowBox[{"(", RowBox[{"6", "-", RowBox[{"7", " ", SuperscriptBox["yc", "2"]}]}], ")"}], " ", RowBox[{"Gamma", "[", FractionBox["1", "6"], "]"}], " ", RowBox[{"Gamma", "[", RowBox[{ RowBox[{"-", FractionBox["1", "6"]}], ",", RowBox[{"-", FractionBox["1", SuperscriptBox["yc", "2"]]}]}], "]"}]}], "+", RowBox[{"36", " ", SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], RowBox[{"2", "/", "3"}]], " ", RowBox[{"Gamma", "[", FractionBox["2", "3"], "]"}], " ", RowBox[{"Gamma", "[", RowBox[{ FractionBox["1", "3"], ",", RowBox[{"-", FractionBox["1", SuperscriptBox["yc", "2"]]}]}], "]"}]}], "-", RowBox[{"24", " ", SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], RowBox[{"2", "/", "3"}]], " ", SuperscriptBox["yc", "2"], " ", RowBox[{"Gamma", "[", FractionBox["2", "3"], "]"}], " ", RowBox[{"Gamma", "[", RowBox[{ FractionBox["1", "3"], ",", RowBox[{"-", FractionBox["1", SuperscriptBox["yc", "2"]]}]}], "]"}]}]}], ")"}]}]}]}], ")"}]}]}]], "Output", CellChangeTimes->{ 3.83715241847911*^9, 3.837152473882375*^9, 3.837152648741281*^9, 3.837152714070775*^9, 3.8371529158052588`*^9, 3.8371530773164387`*^9, { 3.837153147250849*^9, 3.8371531813231363`*^9}, {3.8371532786171103`*^9, 3.83715328927947*^9}, 3.837153512224599*^9}, CellLabel->"Out[37]=",ExpressionUUID->"efae8b97-cac0-4dfb-b6a6-403456514c4f"] }, Open ]], Cell[BoxData[ RowBox[{"i2", "=", RowBox[{"FullSimplify", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"i1", "/.", RowBox[{"{", RowBox[{ RowBox[{"y", "->", RowBox[{"b", RowBox[{"(", RowBox[{"\[Theta]c", "-", "\[Theta]"}], ")"}]}]}], ",", RowBox[{"yc", "->", RowBox[{"b", " ", "\[Theta]c"}]}]}], "}"}]}], ")"}], "+", RowBox[{"(", RowBox[{"i1", "/.", RowBox[{"{", RowBox[{ RowBox[{"y", "->", RowBox[{"b", RowBox[{"(", RowBox[{"\[Theta]c", "+", "\[Theta]"}], ")"}]}]}], ",", RowBox[{"yc", "->", RowBox[{"b", " ", "\[Theta]c"}]}]}], "}"}]}], ")"}]}], ",", RowBox[{"Assumptions", "->", RowBox[{"{", RowBox[{ RowBox[{"b", ">", "0"}], ",", RowBox[{"\[Theta]c", ">", "0"}], ",", RowBox[{"\[Theta]", "<", "\[Theta]c"}], ",", RowBox[{"\[Theta]", ">", RowBox[{"-", "\[Theta]c"}]}]}], "}"}]}]}], "]"}]}]], "Input", CellChangeTimes->{{3.837153261498691*^9, 3.8371533873000507`*^9}, { 3.837153547543831*^9, 3.837153576103598*^9}, 3.8371549111209583`*^9}, CellLabel->"In[52]:=",ExpressionUUID->"c8a865bb-cad1-450c-8187-a809f109c3c8"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"iSing", "=", RowBox[{ RowBox[{"-", FractionBox["1", RowBox[{"12", " ", "\[Pi]"}]]}], RowBox[{"(", RowBox[{ FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox["1", RowBox[{ SuperscriptBox["b", "2"], " ", SuperscriptBox[ RowBox[{"(", RowBox[{"\[Theta]", "-", "\[Theta]c"}], ")"}], "2"]}]]}]], " ", RowBox[{"ExpIntegralE", "[", RowBox[{ FractionBox["7", "6"], ",", RowBox[{"-", FractionBox["1", RowBox[{ SuperscriptBox["b", "2"], " ", SuperscriptBox[ RowBox[{"(", RowBox[{"\[Theta]", "-", "\[Theta]c"}], ")"}], "2"]}]]}]}], "]"}], " ", RowBox[{"Gamma", "[", FractionBox["1", "6"], "]"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{"\[Theta]", "-", "\[Theta]c"}], ")"}], "2"]], "+", FractionBox[ RowBox[{"4", " ", "b", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox["1", RowBox[{ SuperscriptBox["b", "2"], " ", SuperscriptBox[ RowBox[{"(", RowBox[{"\[Theta]", "-", "\[Theta]c"}], ")"}], "2"]}]]}]], " ", RowBox[{"ExpIntegralE", "[", RowBox[{ FractionBox["5", "3"], ",", RowBox[{"-", FractionBox["1", RowBox[{ SuperscriptBox["b", "2"], " ", SuperscriptBox[ RowBox[{"(", RowBox[{"\[Theta]", "-", "\[Theta]c"}], ")"}], "2"]}]]}]}], "]"}], " ", RowBox[{"Gamma", "[", FractionBox["2", "3"], "]"}]}], RowBox[{"\[Theta]", "-", "\[Theta]c"}]]}], ")"}]}]}]], "Input", CellChangeTimes->{{3.837159003786721*^9, 3.837159023124399*^9}}, CellLabel->"In[56]:=",ExpressionUUID->"bd0a8b22-65f5-4df7-b577-0bf969c52194"], Cell[BoxData[ RowBox[{"-", FractionBox[ RowBox[{ FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox["1", RowBox[{ SuperscriptBox["b", "2"], " ", SuperscriptBox[ RowBox[{"(", RowBox[{"\[Theta]", "-", "\[Theta]c"}], ")"}], "2"]}]]}]], " ", RowBox[{"ExpIntegralE", "[", RowBox[{ FractionBox["7", "6"], ",", RowBox[{"-", FractionBox["1", RowBox[{ SuperscriptBox["b", "2"], " ", SuperscriptBox[ RowBox[{"(", RowBox[{"\[Theta]", "-", "\[Theta]c"}], ")"}], "2"]}]]}]}], "]"}], " ", RowBox[{"Gamma", "[", FractionBox["1", "6"], "]"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{"\[Theta]", "-", "\[Theta]c"}], ")"}], "2"]], "+", FractionBox[ RowBox[{"4", " ", "b", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox["1", RowBox[{ SuperscriptBox["b", "2"], " ", SuperscriptBox[ RowBox[{"(", RowBox[{"\[Theta]", "-", "\[Theta]c"}], ")"}], "2"]}]]}]], " ", RowBox[{"ExpIntegralE", "[", RowBox[{ FractionBox["5", "3"], ",", RowBox[{"-", FractionBox["1", RowBox[{ SuperscriptBox["b", "2"], " ", SuperscriptBox[ RowBox[{"(", RowBox[{"\[Theta]", "-", "\[Theta]c"}], ")"}], "2"]}]]}]}], "]"}], " ", RowBox[{"Gamma", "[", FractionBox["2", "3"], "]"}]}], RowBox[{"\[Theta]", "-", "\[Theta]c"}]]}], RowBox[{"12", " ", "\[Pi]"}]]}]], "Output", CellChangeTimes->{{3.837159019621305*^9, 3.837159023641695*^9}}, CellLabel->"Out[56]=",ExpressionUUID->"3576d1fd-5496-4abf-b6de-d73b330829ac"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"SeriesCoefficient", "[", RowBox[{"iSing", ",", RowBox[{"{", RowBox[{"\[Theta]", ",", "\[Theta]c", ",", "1"}], "}"}], ",", RowBox[{"Assumptions", "->", RowBox[{"{", RowBox[{ RowBox[{"\[Theta]c", ">", "0"}], ",", RowBox[{"\[Theta]", "<", "\[Theta]c"}], ",", RowBox[{"b", ">", "0"}]}], "}"}]}]}], "]"}], "//", "N"}]], "Input", CellChangeTimes->{{3.837159089333461*^9, 3.8371591098128653`*^9}, { 3.837159228863208*^9, 3.8371592353512363`*^9}, {3.837159280672409*^9, 3.837159285872224*^9}, {3.837159400090641*^9, 3.83715940032231*^9}, { 3.8371594403473797`*^9, 3.837159450811274*^9}}, CellLabel->"In[81]:=",ExpressionUUID->"430dc047-855a-4d97-9e15-6a1b8ed47326"], Cell[BoxData[ RowBox[{"0.14367637572608224`", " ", SuperscriptBox["b", "3"]}]], "Output", CellChangeTimes->{ 3.837159110237748*^9, {3.837159229380753*^9, 3.837159235660652*^9}, { 3.837159281262404*^9, 3.837159286279546*^9}, 3.8371594008606462`*^9, { 3.837159440963653*^9, 3.837159451210906*^9}}, CellLabel->"Out[81]=",ExpressionUUID->"4c7ac351-a724-4419-9c0c-ec0428c101c7"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"t1", "=", RowBox[{ FractionBox["1", "\[Pi]"], RowBox[{"Integrate", "[", RowBox[{ FractionBox[ RowBox[{ SuperscriptBox["x", RowBox[{ RowBox[{"-", "7"}], "/", "3"}]], RowBox[{"Exp", "[", RowBox[{ RowBox[{"-", "1"}], "/", SuperscriptBox["x", "2"]}], "]"}]}], SuperscriptBox["x", RowBox[{"n", "+", "1"}]]], ",", RowBox[{"{", RowBox[{"x", ",", "0", ",", "\[Infinity]"}], "}"}], ",", RowBox[{"Assumptions", "->", RowBox[{"{", RowBox[{ RowBox[{"x", ">", "0"}], ",", RowBox[{"y", ">", "0"}], ",", RowBox[{"yc", ">", "0"}], ",", RowBox[{"n", ">", "0"}], ",", RowBox[{"n", "\[Element]", "Integers"}]}], "}"}]}]}], "]"}]}]}]], "Input", CellChangeTimes->{{3.837159052820043*^9, 3.837159080308568*^9}, { 3.837159125069777*^9, 3.837159126237418*^9}, {3.837159159774912*^9, 3.837159167254078*^9}, {3.837159215647285*^9, 3.8371592157911263`*^9}, { 3.8371592706725616`*^9, 3.8371592715678377`*^9}}, CellLabel->"In[71]:=",ExpressionUUID->"608d5852-49d9-45d7-9763-cb10da65d830"], Cell[BoxData[ FractionBox[ RowBox[{"Gamma", "[", RowBox[{ FractionBox["7", "6"], "+", FractionBox["n", "2"]}], "]"}], RowBox[{"2", " ", "\[Pi]"}]]], "Output", CellChangeTimes->{{3.837159072418518*^9, 3.83715908203614*^9}, 3.837159127950366*^9, 3.837159168943486*^9, 3.837159217563138*^9, 3.837159273575577*^9}, CellLabel->"Out[71]=",ExpressionUUID->"6e748a34-f946-4d43-8801-0833c28048e5"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"t1", "/.", RowBox[{"n", "->", "0"}]}], "//", "N"}]], "Input", CellChangeTimes->{{3.837159127206106*^9, 3.837159130805312*^9}, { 3.837159222143116*^9, 3.837159238231101*^9}, {3.8371592776402197`*^9, 3.837159288088141*^9}, {3.837159404642499*^9, 3.8371594052183743`*^9}}, CellLabel->"In[78]:=",ExpressionUUID->"5fe4ff16-26ea-47b3-a961-9c776acf99bb"], Cell[BoxData["0.14765111774913997`"], "Output", CellChangeTimes->{ 3.837159131124867*^9, 3.8371591696139517`*^9, {3.837159218839643*^9, 3.8371592385873203`*^9}, {3.837159274439986*^9, 3.837159288334805*^9}, 3.837159405552266*^9}, CellLabel->"Out[78]=",ExpressionUUID->"cffc2fe3-b3b9-420d-ae4c-5559b2e114bb"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{ FractionBox[ RowBox[{"4", " ", "b", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox["1", RowBox[{ SuperscriptBox["b", "2"], " ", SuperscriptBox[ RowBox[{"(", RowBox[{"\[Theta]", "-", "\[Theta]c"}], ")"}], "2"]}]]}]], " ", RowBox[{"ExpIntegralE", "[", RowBox[{ FractionBox["5", "3"], ",", RowBox[{"-", FractionBox["1", RowBox[{ SuperscriptBox["b", "2"], " ", SuperscriptBox[ RowBox[{"(", RowBox[{"\[Theta]", "-", "\[Theta]c"}], ")"}], "2"]}]]}]}], "]"}], " ", RowBox[{"Gamma", "[", FractionBox["2", "3"], "]"}]}], RowBox[{"\[Theta]", "-", "\[Theta]c"}]], "+", RowBox[{ RowBox[{"(", RowBox[{"-", FractionBox[ RowBox[{"4", " ", "b", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox["1", RowBox[{ SuperscriptBox["b", "2"], " ", SuperscriptBox[ RowBox[{"(", RowBox[{"\[Theta]", "+", "\[Theta]c"}], ")"}], "2"]}]]}]], " ", RowBox[{"ExpIntegralE", "[", RowBox[{ FractionBox["5", "3"], ",", RowBox[{"-", FractionBox["1", RowBox[{ SuperscriptBox["b", "2"], " ", SuperscriptBox[ RowBox[{"(", RowBox[{"\[Theta]", "+", "\[Theta]c"}], ")"}], "2"]}]]}]}], "]"}]}], RowBox[{"\[Theta]", "+", "\[Theta]c"}]]}], ")"}], " ", RowBox[{"Gamma", "[", FractionBox["2", "3"], "]"}]}]}], "/.", RowBox[{"\[Theta]", "->", "0"}]}]], "Input", CellChangeTimes->{{3.837156673000149*^9, 3.837156678457199*^9}}, CellLabel->"In[54]:=",ExpressionUUID->"c1119876-e977-4ca0-8963-ee017dea27fd"], Cell[BoxData[ RowBox[{"-", FractionBox[ RowBox[{"8", " ", "b", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox["1", RowBox[{ SuperscriptBox["b", "2"], " ", SuperscriptBox["\[Theta]c", "2"]}]]}]], " ", RowBox[{"ExpIntegralE", "[", RowBox[{ FractionBox["5", "3"], ",", RowBox[{"-", FractionBox["1", RowBox[{ SuperscriptBox["b", "2"], " ", SuperscriptBox["\[Theta]c", "2"]}]]}]}], "]"}], " ", RowBox[{"Gamma", "[", FractionBox["2", "3"], "]"}]}], "\[Theta]c"]}]], "Output", CellChangeTimes->{{3.837156674332121*^9, 3.837156679177827*^9}}, CellLabel->"Out[54]=",ExpressionUUID->"af9d10e0-c1b2-44e5-b2ed-d286376e73b8"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"i2", "/.", RowBox[{"\[Theta]", "->", "0"}]}]], "Input", CellChangeTimes->{{3.837153471782505*^9, 3.837153475437471*^9}}, CellLabel->"In[39]:=",ExpressionUUID->"54955d9a-a1b8-49ea-8b22-6f8c2c06fec5"], Cell[BoxData["0"], "Output", CellChangeTimes->{3.837153475703751*^9, 3.837153517514927*^9}, CellLabel->"Out[39]=",ExpressionUUID->"06308438-871f-452c-aef0-07985232fef7"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"i3", "=", RowBox[{"Simplify", "[", RowBox[{"i2", "/.", RowBox[{"{", RowBox[{ RowBox[{"b", "->", "1"}], ",", RowBox[{"\[Theta]c", "->", "1"}]}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.8371541718671207`*^9, 3.837154174051682*^9}}, CellLabel->"In[46]:=",ExpressionUUID->"3b3622a7-fadc-4e42-a347-8de6c9595dd3"], Cell[BoxData[ RowBox[{"-", RowBox[{ FractionBox["1", RowBox[{"12", " ", "\[Pi]"}]], RowBox[{"(", RowBox[{ FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox["1", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Theta]"}], ")"}], "2"]]}]], " ", RowBox[{"ExpIntegralE", "[", RowBox[{ FractionBox["7", "6"], ",", RowBox[{"-", FractionBox["1", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Theta]"}], ")"}], "2"]]}]}], "]"}], " ", RowBox[{"Gamma", "[", FractionBox["1", "6"], "]"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Theta]"}], ")"}], "2"]], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", FractionBox[ RowBox[{"2", " ", RowBox[{"ExpIntegralE", "[", RowBox[{ FractionBox["7", "6"], ",", RowBox[{"-", "1"}]}], "]"}]}], "\[ExponentialE]"]}], "+", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox["1", SuperscriptBox[ RowBox[{"(", RowBox[{"1", "+", "\[Theta]"}], ")"}], "2"]]}]], " ", RowBox[{"ExpIntegralE", "[", RowBox[{ FractionBox["7", "6"], ",", RowBox[{"-", FractionBox["1", SuperscriptBox[ RowBox[{"(", RowBox[{"1", "+", "\[Theta]"}], ")"}], "2"]]}]}], "]"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{"1", "+", "\[Theta]"}], ")"}], "2"]]}], ")"}], " ", RowBox[{"Gamma", "[", FractionBox["1", "6"], "]"}]}], "+", FractionBox[ RowBox[{"4", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox["1", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Theta]"}], ")"}], "2"]]}]], " ", RowBox[{"ExpIntegralE", "[", RowBox[{ FractionBox["5", "3"], ",", RowBox[{"-", FractionBox["1", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Theta]"}], ")"}], "2"]]}]}], "]"}], " ", RowBox[{"Gamma", "[", FractionBox["2", "3"], "]"}]}], RowBox[{ RowBox[{"-", "1"}], "+", "\[Theta]"}]], "+", RowBox[{ RowBox[{"(", RowBox[{ FractionBox[ RowBox[{"8", " ", RowBox[{"ExpIntegralE", "[", RowBox[{ FractionBox["5", "3"], ",", RowBox[{"-", "1"}]}], "]"}]}], "\[ExponentialE]"], "-", FractionBox[ RowBox[{"4", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox["1", SuperscriptBox[ RowBox[{"(", RowBox[{"1", "+", "\[Theta]"}], ")"}], "2"]]}]], " ", RowBox[{"ExpIntegralE", "[", RowBox[{ FractionBox["5", "3"], ",", RowBox[{"-", FractionBox["1", SuperscriptBox[ RowBox[{"(", RowBox[{"1", "+", "\[Theta]"}], ")"}], "2"]]}]}], "]"}]}], RowBox[{"1", "+", "\[Theta]"}]]}], ")"}], " ", RowBox[{"Gamma", "[", FractionBox["2", "3"], "]"}]}]}], ")"}]}]}]], "Output", CellChangeTimes->{3.8371541745414057`*^9}, CellLabel->"Out[46]=",ExpressionUUID->"caca42d0-5916-43b6-84ea-cc91bc115d17"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Plot", "[", RowBox[{"i3", ",", RowBox[{"{", RowBox[{"\[Theta]", ",", RowBox[{"-", "1.2"}], ",", "1.2"}], "}"}], ",", RowBox[{"PlotRange", "->", "All"}], ",", RowBox[{"WorkingPrecision", "\[Rule]", "20"}]}], "]"}]], "Input", CellChangeTimes->{{3.8371534184294147`*^9, 3.837153465381342*^9}, { 3.8371538031322393`*^9, 3.8371538039795713`*^9}, {3.837154155242539*^9, 3.8371541836983547`*^9}, {3.837154249916608*^9, 3.837154260835929*^9}}, CellLabel->"In[50]:=",ExpressionUUID->"bdc138d2-dd1c-49d9-b891-f23d7c6a21a1"], Cell[BoxData[ GraphicsBox[{{{}, {}, TagBox[ {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[ 1.], LineBox[CompressedData[" 1:eJwBsQFO/iFib1JlAgAAABoAAAACAAAAaHXmkTf68r9xFzYS7cGzPwE+EGPq +fK/KsnCecvAsz8oFLKO5vbyvxbpET+AtbM/T+pTuuLz8r+ZTc+TPaqzP3bA 9eXe8PK/IZNObAOfsz/FbDk91+ryv3OiaHqoiLM/7ELbaNPn8r8z6zCZh32z PxQZfZTP5PK/zFEVDm9ysz9ixcDrx97yv+hAqs1WXLM//h1ImrjS8r/WfMBx iDCzPyX06cW0z/K/D8JvDqklsz9MyovxsMzyv81nqarRGrM/mnbPSKnG8r9b AjW4OgWzPzfPVveZuvK/vTzkCWvasj9wgGVUe6LyvzxQ5d87hrI/4eKCDj5y 8r9jRgUuXeOxP+4mFE75bvK/xfNMTZXYsT/6aqWNtGvyv3oqu3jVzbE/FPPH DCtl8r9GwSTLbbixP0gDDQsYWPK/AaR4RP2NsT+vI5cH8j3yv+h88P2OOrE/ fWSrAKYJ8r8gmrMQPZmwPxnm0/INofG/qJL99evVrj95ptylut3wv0k0rqP3 6ao/bg2pNjwe8L83VmndGJSnPw5HTHsiA/C/3y4hdDUkpz9INuaO "]], LineBox[CompressedData[" 1:eJw1l2k01XvcxUkpGZIoFGWWcEyJ6P9FMpXMbmTMVGYu4RSFrplrCsmQ4WQo 4hii8DMex+zIFFH3yCFDhjJU4ulZ63le7LVf7Jd7rb0+W+CWl7HTPjo6uq9/ 9L8+6+PC9PbHClr+3cOYc7gIOo7bX7rGs4zsoqybT47mQYfm2OIp/0WkLJM9 RarPgTU0bJY+/gX9+/eh4e9emeBeo27d89c86meTEL39KQ26s52somtpqFT0 UmSYVwoE0P3YppOdRa5W8/wjyYlgWHHMSLBqBp0spmOYPR0HLCazz/o4qWi9 dPJWwZVIOOEtmiCX/AmZNxt5+aiHAyvdJZrG9BSaFzzHt/QuBC6/URHzK5pA g768XJ5YIJy4w1PAXTyGrmbGxB+Z84Wq345sFJdhRA1E18cm3WC86HdD2c1h dIVsKjRU7wY7RmluMQbDaGHwMH4m3Q20nneSL18YRudqmrN0Td1g4rpEeC3j MOJmd1IT6neFfc9Wtp4S3qFH84dMD5HugKEG/j/nmSGkZuW2cZ/kAssR/1bv 2FDQKULFpgKTI5TMzxJHjCjIfFA93WDBAZz0VCvLNSmosQcfGt7jAFOsX8rs JCioxA8ncCHBAfpSNYraNwYRJveEo4/TAcryN57ExQ2iCs+E9Veit2Di/utd Jb8BZJ02KaxgbQfGao6/FiwGkPLF91Ve6nbQzXB0OwsGkP7NVMZ2ETt4E+O6 Rsc8gJ6pCQx4TNlCZgbfZ3JePyLtKpjFMNiCZXV4142BPsR04DuRzsIaJhcN UgMlepFMUWViubgliArXSBLYe9G2lX+/1LYFeFvxdlA2e5D14A3aZKcFMPZ9 3pBo60FTeIeYdhcLkCkP+mvqZg9Soeh+JZXcgAjvAl6N+G5Udn8lKUbxL5Db 2HzGvEpGzU+8Oai+phAsZaWsNEpGam6rXjI6pkB2aqE4NpARtals5iWfKViP xu5riiIjr2WRLVyXCUS9PuPoLUBGw312MjaCJjAdpCc2YtyJut1kOc2njCB6 N7ssp6YDyTMP+ZxavA6B0hMfI7I7kIBAwAG+59fhts1xDq9/OtAM/uALd/vr oNOUcBfMOpAi9RLLj3F9YAx5AJ++tyOjeeerHT3X4NFvO4qAQjuyd1XbE27W g5AdoY1CYitCGFzpm9ECT0k7sfjMVnTw3MW/84q1wMYqy8I/rBVJLs50tnlo waUGzqYrxq0IL2rSvb59BXbuHYikrbUgdn3OoB8cVyDoF41bXLYFbf6HIxeY XYa57yZB5rYIbf+TbK/GpQaC8m2nuM8hNCScf0IfAdj4yDW/32xCLGcLzQZM AUaX2Q9aJzShLCF9s41QDEi0vlSHxkbE26w10fFZFQrHtSu8TzYgwWCDAx69 yvDp+GsTubm3qFnZHp/krgynzES3vhHfoqqELDo3VmVIoezHAvTeIon1KdZL RkoQ1t3SExz0BuW2HOu/+Z8i2DeozMWM1aFwVmu9ICEFyPpZGnOtoA4VJ+M4 v/TLw7gSrzSbVx3yPhuLY70nD4Y1235JjHVos8ZA2ndEDtTKa/ZlnH+NJMQe 8ZQnywL/Mxw/IaUGxc84tpOtpWEyXMgMGRCRhSNNhPmDCNhVqklu8xNRipVn PtNtEaBNWzPIfq1Et650OKhtCMO6ckZlflwlatLffGB+TBgOrbKwRXRVINPq 7SohG0Fwr1uKfxJRjti5SFaeVvwQeExMlVmkHPm7RJ6+/ZQPwj3tF+63lSFp 9+jZwQ+nIFN4VMuOrgwVe+QLrTifhM5EtCeKf4G2LmEcVhncIHgn2afaoxgZ qfVvznodA+n23tOiLMUoJI2mShrhAOXTB/vTS4uQqKCmDifGAYYj9yTuzT1H bu0Sd35wHYVgdacZdXsCcmDt3urpYIWorNykqr1C5Lc+KexUxwKpW+9BJKcQ RduEDY+/YobSsutZhyYLUEuhioZBGROM8yj96T0fLajEifhTD8BnP18G9W95 yNUsVPbjr/2wMvCykpiUhyTqfk4K8u4HxggBtrT+Z8hX9jhB0WkfOH/mSqhP zEGrw2xFvG67WHBmzEJWXzbaz5GmGL64g6Ua7mk9PJyNXh96EZbt/wtraZjf 03r0FD1HIQTa823s1OO3PsMBGeiIpl0RKfgbRl779iKuLh1JHgm1LjBax/yv S9I0f6ShlfHyM4PSa1j/wWzLGvxjJC0nky13+CuGdxxJ9Xibioboz2+77ixi oi2sAyI7KSjfwytQg34Be4APuZwWnIyo9VcP92vQMMmx2uDrTUko22GrMvTB Z2xcfuU1414i+od3qfRdPxWzn0A3vtL9izoMxp8X6H3E3qTI9NmExSPi2U3j Xzc+YBz6eeoDDHHIEKpVEuLeY+0oTKKCMRq5Y64xjg7DGF/Qt5wzUZHouKOF TRHjEHZXzvFYElME4jpd/lgnZwDrXxyOpI99hF7rGgvWBvVgYoQrOz4s4eiq rPNyWm4n9tCm1psaH4qW3yxGLxxpx8ZPiM0aH3mI+AmnUyq6m7Fb+biqCIFg 1DawfoOS+xbj5Dqa5SKMR23S5gwiUIuRItcf6YgFoLBkjfibf1diAT/feZyV 8EPs68zT84GlmLhHjflhKR+kz066Na1QgFUjqQejtu6oLedyrdT5TOwEYdUw w8kZVYUvKQiZJmL4mCpBSzc7xMTLYqoUH4q9+bg0f4Fqig620u81OrpjdGaz 3AGT6khtzsvkZasmNmim+pBRURyK5UrEmzmPYfMPGCzwm/ogxhb7yErZGpPe +BVVv2IJGr/liLdyAjAFrrfrvL4OwBLwslklJwqj4LbPWLjfAZ6VUW7Rw48x T73zBunOXtCZGyVqwfIMi5Crb3Kt8IVr2NKEfeNzbBiNjUSX+UNc3jxNKqMc E9TfXCouDQTdYCtO2/4qzGeCcz+56B5UUpoHomfrsGYX+ZNzhSEwdaG4Jj2p CYv8Kc7DQHkIJzS0qzTvtGLS6nWFyQZhIGVK7MINd2Cjkdoygv3h0ORPDomP 7sKEuZy11XsiIFvDT4nyZRDrubkxNKgbBdIZcE1K4R3mm//Ixo4cDe/+Wn7S VjCCteAK/B90xIFKT+vP3PUJ7PZdOfojmgnAaM1fSKNMYeyNLXE5rf+CjiDV QsHtE4YTWaVM2ybBqvUgLrFoBgsjv04XPZ0MypRgixf2s9io2wNrz+lkUJjx W3aUnsNCiGxfflulQvBTZsHqvQWsF6To+CzTQdl/9vbn5jWs/A5Pd+7ZDHiZ HSq6Sv8NS0w5kCr4IwPaNQmda/rfMZO5KTHxjEyYVe+4Nca4hU0kJBjIj2aD dMF5dtmdHayhDs9TS8iBjqNT31ZTd7FsqvOMkn8unOC/yW96ig5knuaqnI55 Bj8Zd5dCBunh+9eAAeGz+bBY4NdDwB+AWQZPBVpiPuwc4A6g+TLCKLfjk6Lt fLgWtNxW4H0QXmsYOkiQCyD24HFS6n0mwKeJb+FcCCBxblTRdJQV6LH3/CqE YjhDVGuoJh+DdeOB8B2WErihq1vQ6s0JVJeO+Ua/EuhZvsGWeJIL2pMqiWqa peDN87vrTtBxiJqN1tKaeQHS+inWJ815gC3+oqeRwCsoDWk0j0s5Dbt5Mu84 ol+Ba4b8h375M7BSK6o0vPoK9lPHjivlnoHb3kyh5+wqwNx+xF4SLwB6mRt6 5bhKmL5UT81WEwK2r33TVf1E4Ptl0FV5RAxUmZ5a8+0R4XCRk6NMpBi4Cd/+ EIGrgud4+zHqnhiQLfdNWCRWAYfHufsz38UhrFNxZM+oGs70luiZLUvA1rPc nqsjNXDTba3hIDcORBvc9WoO1EIq7eXj7nwcmI4pd/Er1sJEkB+P/QkZILIO k9bSamHgvt9vMYIMuOMPtabfeA1DO+Ht3zpl4ZOJT93MZB2kLpSXxEgoAJsX dkGfpR4oNyJOpwUpwKUY5tpa1XrgedL9QLVLATKbCVXR2fVAteUVtnQ9D6ZS E+U42zdQ9GfI31cpQhfjZcI96lvQv3v9r+8OykCs50w+9qUJeC4e2DehBPBo c4ESy4dAdFd/L7sAwFy+5eh+YwR7Pbi3cQfV4NdLj6TvbxDYVNJHnPT7wyOG DmRXrWYgoE2lSCN1cL3c0HmVuwWYFp9XZPNdhp9nPUmsDa1QMLbPa/enFjjo WNB0R1vB+cnW0TkVbeh11mSMWG2FoZfrvpLB2pBTwKu1K9wG2tcik/vodECD r7Pja3wbSEh/uOzKogux7Gc6BmzbQcyZJ5Dr/FU4tUVpS2IgQd2QhPjdXgMg PtuvY8JHAqrSzn4KsyHo6l3o5bxAAgf9x6qpsoZwNztrON2VBOxD4yn5wYYw cNllNptCAhURqtfHE0YQlrjDWJrbCbhgktW0qTHMnRXTa1Xpgu6XHYbH10yB aHV/cP3vXqgn+18sr7KEQrbxN9uxvXDMTjjO9z9LeNwsT9gt+JMbaIfbHrkJ gcKLQYeHewGlBgZUud2ES0sWQoIKfVDPlXxEU9QKyPcu3DX63gcT9uKDLAXW MJ25frLCbwBCcWW8xD88/497dVXtRQrE1L1/qVrtALR4YYtGbQoIqErNlgw5 gM6r1N02Uwp45/+aw6392eF1P12KJwVuMevYREg7QnLA+Y+L+RR4bHMh0rrU EXLCag8LMA+BVo1lYd1zJ6jLqLOPnRwC7Q8cFepEF1hobzhid38YfiZKU0q4 3KFSnCaTFjEGkT//Fix55QsC63STgaETUPmekKi4GQDdI2WPtRunwOfeVk5f bAi4URtjKHaf4KnkdVIZKQxY55aN/RmoUEpPHqFyRoKdT3mwbcwMjPP6JU/n xEJU+dxSLtMszCvrM1/QT4Te7TuO5iU0YMfXrBmbpUAL+T5mIjUPe/jfVZOh aTBEGD9DonwBZQYhyXmbTKg7X6b9r+ciXNtZWNQsygGDqxyUnJUl4KVfXc5K yYN9ns23hJO+Qt6D9rws80L4vz8L//9n/wfRPlcP "]], LineBox[CompressedData[" 1:eJwVj30s1HEAhy+vvSwrtqgtRUkxjqQXL32iKw1/SJj3RM0ppeby0q3GGkUl L9fy0o5ydpSd391XI+uISridsoR1xo2Rl2vu903NyyL98ezZ8+djE58SfMmA w+HYr/HfZn6BQltDPUadbQTue6XoParxy+vRo36+m2/6+CUWkqP438pYnHkj qd4kk6G/4IWd5hkL4+/CuNetMjDM9Ji2kkXqW/ZC0mcZ+PPpMboaFt7kfpRe L8OQUBRsqGDhblHqLzzUgOZ8tZdrDwtRUWV4gbIBmbXY+miJxZZZ7eXlKwzC enJ7i/+y6Jj8aLkqYHBQp35QyqFYyDapWH+HwQw32kRiQvHqp9TbqZBBZFPG cos5Rd9VhxplIwOPTvnElAOFqmGlgbPKYGl8TwsviqKNu7sNpXL05n1NlcRS OMZsOMmtlqOae9fJKJ6i9nzCXnuZHIHCsar3fApGus/R5YMcYvPn93zTKHhn vbNq9XL4+liHniiiME3UnvMKUGDbD7VZlYhi80jUuqYwBWYf3u5afUoRH5DW 6RGvgGho2LNNTBF6rbE7JlOByZQK2+P1FGLV+HaLWgXyxVZ6z08UfPuR66dN CWJ5XXUVPRRJzYVN4RYEbjPpCctqCnliik6wi2DYfXCgpZ8iP2ggUnWEgNHk Fu4YWuv9KrkxjyAn+7D/LQ1FyI1OV/8gAmf1E+WxMQplns80yycwSD2VXjZB MVcSgdCbBINWv10WpygEjXXLHdkE9a2SmXAdhV201smzgCDrYoikeY4iI6V/ QllOELLRKNbq19qvLsctQEpwgCGWGX8oIoqt7ccJwUpoQt/gIkV5Cf2S844g jrS7OC9R7EzuHu1uJ/gH5cBJxg== "]]}, Annotation[#, "Charting`Private`Tag$591217#1"]& ], {}}, {}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->{True, True}, AxesLabel->{None, None}, AxesOrigin->{0, 0}, DisplayFunction->Identity, Frame->{{False, False}, {False, False}}, FrameLabel->{{None, None}, {None, None}}, FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}}, GridLines->{None, None}, GridLinesStyle->Directive[ GrayLevel[0.5, 0.4]], ImagePadding->All, Method->{ "DefaultBoundaryStyle" -> Automatic, "DefaultGraphicsInteraction" -> { "Version" -> 1.2, "TrackMousePosition" -> {True, False}, "Effects" -> { "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, "Droplines" -> { "freeformCursorMode" -> True, "placement" -> {"x" -> "All", "y" -> "None"}}}}, "DefaultMeshStyle" -> AbsolutePointSize[6], "ScalingFunctions" -> None, "CoordinatesToolOptions" -> {"DisplayFunction" -> ({ (Identity[#]& )[ Part[#, 1]], (Identity[#]& )[ Part[#, 2]]}& ), "CopiedValueFunction" -> ({ (Identity[#]& )[ Part[#, 1]], (Identity[#]& )[ Part[#, 2]]}& )}}, PlotRange->{All, All}, PlotRangeClipping->True, PlotRangePadding->{{ Scaled[0.02], Scaled[0.02]}, { Scaled[0.05], Scaled[0.05]}}, Ticks->{Automatic, Automatic}]], "Output", CellChangeTimes->{{3.8371534327807703`*^9, 3.8371534661081553`*^9}, 3.8371535186967163`*^9, {3.837153795719611*^9, 3.837153804720337*^9}, { 3.837154159258175*^9, 3.8371541842464027`*^9}, {3.837154252931675*^9, 3.8371542615178967`*^9}}, CellLabel->"Out[50]=",ExpressionUUID->"7dab72fd-cd30-44be-b73a-d23ec16b745a"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"FullSimplify", "[", "i1", "]"}]], "Input", CellChangeTimes->{{3.837153188209074*^9, 3.837153192208252*^9}}, CellLabel->"In[27]:=",ExpressionUUID->"acd277d9-1867-4c72-bd06-edff0291a14e"], Cell[BoxData["$Aborted"], "Output", CellChangeTimes->{3.837153204428977*^9}, CellLabel->"Out[27]=",ExpressionUUID->"a786f1b5-8e89-47a1-ba0e-8efd965ebaa9"] }, Open ]], Cell[BoxData[ RowBox[{"i1", "/.", RowBox[{"y", "->"}]}]], "Input", CellChangeTimes->{{3.837153084295353*^9, 3.837153086518475*^9}},ExpressionUUID->"b3b0287b-1e22-4b15-a4f5-\ 10695aa5f35c"], Cell[BoxData[" "], "Input", CellChangeTimes->{ 3.837152869906995*^9},ExpressionUUID->"072b5925-5f74-49be-9a73-\ 34916c763d7c"] }, WindowSize->{636., 350.25}, WindowMargins->{{Automatic, 2.25}, {357.75, Automatic}}, FrontEndVersion->"12.3 for Linux x86 (64-bit) (May 11, 2021)", StyleDefinitions->"Default.nb", ExpressionUUID->"72b00315-9f0b-4e16-a7c6-aa74d14cd275" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[580, 22, 747, 23, 41, "Input",ExpressionUUID->"fababbb8-381a-4687-8bb1-07b91e8a5e30"], Cell[1330, 47, 989, 31, 61, "Output",ExpressionUUID->"18ce2878-724a-4922-9798-030bcc1e0186"] }, Open ]], Cell[CellGroupData[{ Cell[2356, 83, 1367, 36, 42, "Input",ExpressionUUID->"bf6da1af-2ff5-46b6-a7be-85f3bdd44d14"], Cell[3726, 121, 9857, 285, 381, "Output",ExpressionUUID->"efae8b97-cac0-4dfb-b6a6-403456514c4f"] }, Open ]], Cell[13598, 409, 1282, 35, 41, "Input",ExpressionUUID->"c8a865bb-cad1-450c-8187-a809f109c3c8"], Cell[CellGroupData[{ Cell[14905, 448, 1994, 59, 171, InheritFromParent,ExpressionUUID->"bd0a8b22-65f5-4df7-b577-0bf969c52194"], Cell[16902, 509, 1870, 56, 70, "Output",ExpressionUUID->"3576d1fd-5496-4abf-b6de-d73b330829ac"] }, Open ]], Cell[CellGroupData[{ Cell[18809, 570, 757, 16, 24, "Input",ExpressionUUID->"430dc047-855a-4d97-9e15-6a1b8ed47326"], Cell[19569, 588, 387, 7, 25, "Output",ExpressionUUID->"4c7ac351-a724-4419-9c0c-ec0428c101c7"] }, Open ]], Cell[CellGroupData[{ Cell[19993, 600, 1177, 32, 57, "Input",ExpressionUUID->"608d5852-49d9-45d7-9763-cb10da65d830"], Cell[21173, 634, 417, 10, 47, "Output",ExpressionUUID->"6e748a34-f946-4d43-8801-0833c28048e5"] }, Open ]], Cell[CellGroupData[{ Cell[21627, 649, 393, 7, 22, "Input",ExpressionUUID->"5fe4ff16-26ea-47b3-a961-9c776acf99bb"], Cell[22023, 658, 319, 5, 25, "Output",ExpressionUUID->"cffc2fe3-b3b9-420d-ae4c-5559b2e114bb"] }, Open ]], Cell[CellGroupData[{ Cell[22379, 668, 1903, 57, 129, "Input",ExpressionUUID->"c1119876-e977-4ca0-8963-ee017dea27fd"], Cell[24285, 727, 752, 21, 55, "Output",ExpressionUUID->"af9d10e0-c1b2-44e5-b2ed-d286376e73b8"] }, Open ]], Cell[CellGroupData[{ Cell[25074, 753, 227, 4, 22, "Input",ExpressionUUID->"54955d9a-a1b8-49ea-8b22-6f8c2c06fec5"], Cell[25304, 759, 171, 2, 25, "Output",ExpressionUUID->"06308438-871f-452c-aef0-07985232fef7"] }, Open ]], Cell[CellGroupData[{ Cell[25512, 766, 373, 9, 24, "Input",ExpressionUUID->"3b3622a7-fadc-4e42-a347-8de6c9595dd3"], Cell[25888, 777, 3799, 115, 282, "Output",ExpressionUUID->"caca42d0-5916-43b6-84ea-cc91bc115d17"] }, Open ]], Cell[CellGroupData[{ Cell[29724, 897, 570, 11, 24, "Input",ExpressionUUID->"bdc138d2-dd1c-49d9-b891-f23d7c6a21a1"], Cell[30297, 910, 8540, 161, 184, "Output",ExpressionUUID->"7dab72fd-cd30-44be-b73a-d23ec16b745a"] }, Open ]], Cell[CellGroupData[{ Cell[38874, 1076, 211, 3, 24, "Input",ExpressionUUID->"acd277d9-1867-4c72-bd06-edff0291a14e"], Cell[39088, 1081, 156, 2, 25, "Output",ExpressionUUID->"a786f1b5-8e89-47a1-ba0e-8efd965ebaa9"] }, Open ]], Cell[39259, 1086, 196, 5, 22, "Input",ExpressionUUID->"b3b0287b-1e22-4b15-a4f5-10695aa5f35c"], Cell[39458, 1093, 129, 3, 22, "Input",ExpressionUUID->"072b5925-5f74-49be-9a73-34916c763d7c"] } ] *) (* End of internal cache information *)