Skip to content

Instantly share code, notes, and snippets.

@Parcly-Taxel

Parcly-Taxel/ellint.nb

Created November 3, 2021 06:04
Show Gist options
  • Save Parcly-Taxel/7f2e3ed568aa974f63f592154e32ae0a to your computer and use it in GitHub Desktop.
Save Parcly-Taxel/7f2e3ed568aa974f63f592154e32ae0a to your computer and use it in GitHub Desktop.
Download ZIP
Trefoil knot arc length (MSE #2569180)
Raw
ellint.nb
(* 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[ 70385, 1855]
NotebookOptionsPosition[ 63663, 1740]
NotebookOutlinePosition[ 64055, 1756]
CellTagsIndexPosition[ 64012, 1753]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
RowBox[{
RowBox[{"I0", "[", "y_", "]"}], ":=",
RowBox[{"NIntegrate", "[",
RowBox[{
RowBox[{"Sqrt", "[",
RowBox[{"17", "+",
RowBox[{"8",
RowBox[{"Cos", "[",
RowBox[{"3", "t"}], "]"}]}], "+",
RowBox[{"9",
RowBox[{
RowBox[{"Cos", "[",
RowBox[{"3", "t"}], "]"}], "^", "2"}]}]}], "]"}], ",",
RowBox[{"{",
RowBox[{"t", ",", "0", ",", "y"}], "}"}], ",",
RowBox[{"WorkingPrecision", "\[Rule]", "50"}]}], "]"}]}]], "Input",
CellChangeTimes->{{3.844574363818349*^9, 3.84457439383665*^9}, {
3.844574484213416*^9, 3.8445744850146313`*^9}, {3.844574555496615*^9,
3.8445745816212177`*^9}, {3.844574627715281*^9, 3.844574671425989*^9}, {
3.844574716712289*^9, 3.844574740017115*^9}, {3.844574770185299*^9,
3.8445748304959707`*^9}, {3.844575569751794*^9, 3.844575602818334*^9}, {
3.844812280984931*^9, 3.84481228124212*^9}, {3.844812316762014*^9,
3.844812371681035*^9}, {3.8448605626432333`*^9, 3.844860624650526*^9}, {
3.84486116505408*^9, 3.8448611911881447`*^9}, 3.8448719569229927`*^9,
3.844902596064334*^9},
CellLabel->"In[1]:=",ExpressionUUID->"2c2a9424-a764-4819-aa10-799f2cd948c9"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"I0", "[", "1", "]"}]], "Input",
CellChangeTimes->{{3.8448122575046787`*^9, 3.844812265311509*^9}, {
3.844812400348383*^9, 3.8448124272508593`*^9}, {3.844814357388309*^9,
3.8448143576938753`*^9}, 3.844860633536458*^9, {3.844865686644478*^9,
3.8448657081947813`*^9}, 3.844869660241852*^9},
CellLabel->"In[2]:=",ExpressionUUID->"d08b93dc-4fd6-4d10-ad2c-5f539f664970"],
Cell[BoxData["4.\
60432427101551068986302631297185217317292643002712657214157791124573380280443`\
50."], "Output",
CellChangeTimes->{
3.844812428725093*^9, 3.844814357987006*^9, 3.844815867309887*^9,
3.844822480006954*^9, 3.844825040589384*^9, 3.844826270832842*^9,
3.8448270606828833`*^9, 3.844828676208312*^9, 3.8448361543145933`*^9,
3.844837177224538*^9, 3.844857571415897*^9, 3.844860633979177*^9, {
3.84486117394524*^9, 3.844861192875278*^9}, {3.844865687100872*^9,
3.844865708622086*^9}, 3.8448666994887*^9, {3.844869649693099*^9,
3.844869660567865*^9}, 3.844871909223892*^9, 3.8448719733656178`*^9,
3.8449026138130093`*^9, 3.844903031019682*^9},
CellLabel->"Out[2]=",ExpressionUUID->"18ca9d79-c636-45c2-8a3d-1b1b1ba5e2c3"]
}, Open ]],
Cell[BoxData[
RowBox[{
RowBox[{"I1", "[", "y_", "]"}], ":=",
RowBox[{"NIntegrate", "[",
RowBox[{
FractionBox[
RowBox[{"Sqrt", "[",
RowBox[{
RowBox[{"u", "^", "2"}], "+",
RowBox[{
RowBox[{"(",
RowBox[{"8", "/", "9"}], ")"}], "u"}], "+",
RowBox[{"(",
RowBox[{"17", "/", "9"}], ")"}]}], "]"}],
RowBox[{"Sqrt", "[",
RowBox[{"1", "-",
RowBox[{"u", "^", "2"}]}], "]"}]], ",",
RowBox[{"{",
RowBox[{"u", ",",
RowBox[{"Cos", "[",
RowBox[{"3", "y"}], "]"}], ",", "1"}], "}"}], ",",
RowBox[{"WorkingPrecision", "\[Rule]", "50"}]}], "]"}]}]], "Input",
CellChangeTimes->{{3.844814301009468*^9, 3.844814321304102*^9}, {
3.844815796644169*^9, 3.84481580142778*^9}, {3.8448159126090727`*^9,
3.844815913402286*^9}, {3.844818016625935*^9, 3.844818095733235*^9}, {
3.844818238000065*^9, 3.844818242586886*^9}, 3.844822470239521*^9, {
3.844860642324913*^9, 3.844860659379464*^9}, {3.844861178284354*^9,
3.844861203122271*^9}},
CellLabel->"In[3]:=",ExpressionUUID->"db458be3-0816-4211-97f2-064142b64d1c"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"I1", "[", "1", "]"}]], "Input",
CellChangeTimes->{{3.8448143511646233`*^9, 3.8448143607346973`*^9},
3.844818076035123*^9, {3.844860664297372*^9, 3.844860665683434*^9}, {
3.844865690795096*^9, 3.844865726410351*^9}, 3.844869663299184*^9},
CellLabel->"In[4]:=",ExpressionUUID->"2f6514a5-d960-4f90-8ca6-c54f4a5a8b42"],
Cell[BoxData["4.\
60432427101551068986302631297185217317292643002711789125783600670919398315357`\
50."], "Output",
CellChangeTimes->{{3.844814352979395*^9, 3.844814361057822*^9}, {
3.844815856654619*^9, 3.844815863003393*^9}, {3.84481589702973*^9,
3.844815902159483*^9}, {3.8448159592702627`*^9, 3.844815983083675*^9}, {
3.844818076332498*^9, 3.8448180996458693`*^9}, 3.844818245413706*^9,
3.8448224824323606`*^9, 3.844825045684621*^9, 3.844826270900231*^9,
3.844827060728867*^9, 3.8448286762629538`*^9, 3.84483615440799*^9,
3.8448371772857847`*^9, 3.844857571484582*^9, 3.844860666038741*^9, {
3.844861182055784*^9, 3.8448612046259117`*^9}, {3.844865691093285*^9,
3.8448657267796497`*^9}, 3.844866699627228*^9, {3.844869649814784*^9,
3.84486966365819*^9}, 3.844871909733181*^9, 3.844871973612968*^9,
3.844902614095231*^9, 3.844903031508366*^9},
CellLabel->"Out[4]=",ExpressionUUID->"dc204ccf-5eef-40d8-a14c-950b2f65ef69"]
}, Open ]],
Cell[BoxData[
RowBox[{
RowBox[{"tu", "=",
FractionBox[
RowBox[{
RowBox[{"A", " ", "u"}], "+", "1"}],
RowBox[{"A", "+", "u"}]]}], ";",
RowBox[{"ut", "=",
FractionBox[
RowBox[{
RowBox[{"A", " ", "t"}], "-", "1"}],
RowBox[{"A", "-", "t"}]]}], ";",
RowBox[{"dtu", "=",
FractionBox[
RowBox[{
RowBox[{"A", "^", "2"}], "-", "1"}],
RowBox[{
RowBox[{"(",
RowBox[{"A", "+", "u"}], ")"}], "^", "2"}]]}], ";"}]], "Input",
CellChangeTimes->{{3.844823239985173*^9, 3.8448232676188173`*^9}, {
3.8448234059370728`*^9, 3.844823443797614*^9}, {3.844824267281938*^9,
3.844824271582047*^9}},
CellLabel->"In[5]:=",ExpressionUUID->"1a851957-c3c7-4afc-9874-71074aa704d1"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"With", "[",
RowBox[{
RowBox[{"{",
RowBox[{"tmp", "=",
RowBox[{
RowBox[{
RowBox[{
RowBox[{
FractionBox[
RowBox[{"Sqrt", "[",
RowBox[{
RowBox[{"t", "^", "2"}], "+",
RowBox[{"X", " ", "t"}], "+", "Y"}], "]"}],
RowBox[{"Sqrt", "[",
RowBox[{"1", "-",
RowBox[{"t", "^", "2"}]}], "]"}]], "/.",
RowBox[{"t", "->",
FractionBox[
RowBox[{"1", "+",
RowBox[{"A", " ", "u"}]}],
RowBox[{"u", "+", "A"}]]}]}], "/.",
RowBox[{
RowBox[{"Power", "[",
RowBox[{"expr_", ",", "p_"}], "]"}], "->",
RowBox[{"Power", "[",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"(",
RowBox[{"A", "+", "u"}], ")"}], "^", "2"}], "expr"}], ",",
"p"}], "]"}]}]}], "//", "FullSimplify"}], "//", "ExpandAll"}]}],
"}"}], ",",
RowBox[{
RowBox[{"Sqrt", "[",
RowBox[{"Collect", "[",
RowBox[{
RowBox[{"tmp", "[",
RowBox[{"[",
RowBox[{"2", ",", "1"}], "]"}], "]"}], ",", "u"}], "]"}], "]"}], "/",
RowBox[{"Sqrt", "[",
RowBox[{"FullSimplify", "[",
RowBox[{"tmp", "[",
RowBox[{"[",
RowBox[{"1", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}]}],
"]"}]], "Input",
CellChangeTimes->{{3.8448267293702374`*^9, 3.844826806850891*^9}, {
3.844826855366386*^9, 3.844826866227957*^9}, {3.844827029055757*^9,
3.8448270307536182`*^9}, {3.844827102983036*^9, 3.844827104072507*^9}, {
3.844827158319659*^9, 3.844827166356738*^9}},
CellLabel->"In[6]:=",ExpressionUUID->"7f9843fd-d61b-4867-b7ae-2aae975a95b1"],
Cell[BoxData[
FractionBox[
SqrtBox[
RowBox[{"1", "+",
RowBox[{"A", " ", "X"}], "+",
RowBox[{
SuperscriptBox["A", "2"], " ", "Y"}], "+",
RowBox[{
SuperscriptBox["u", "2"], " ",
RowBox[{"(",
RowBox[{
SuperscriptBox["A", "2"], "+",
RowBox[{"A", " ", "X"}], "+", "Y"}], ")"}]}], "+",
RowBox[{"u", " ",
RowBox[{"(",
RowBox[{
RowBox[{"2", " ", "A"}], "+", "X", "+",
RowBox[{
SuperscriptBox["A", "2"], " ", "X"}], "+",
RowBox[{"2", " ", "A", " ", "Y"}]}], ")"}]}]}]],
SqrtBox[
RowBox[{"-",
RowBox[{"(",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"-", "1"}], "+",
SuperscriptBox["A", "2"]}], ")"}], " ",
RowBox[{"(",
RowBox[{
RowBox[{"-", "1"}], "+",
SuperscriptBox["u", "2"]}], ")"}]}], ")"}]}]]]], "Output",
CellChangeTimes->{{3.844827031649857*^9, 3.8448270609499807`*^9},
3.8448271043769627`*^9, 3.8448271672093*^9, 3.844828676474471*^9,
3.844860685601276*^9, 3.844866699908094*^9, 3.844869650091806*^9,
3.844871910058362*^9, 3.8448719739359617`*^9, 3.8449026144558153`*^9,
3.8449030324229918`*^9},
CellLabel->"Out[6]=",ExpressionUUID->"44b572ec-85eb-44b2-91e1-bb822b11f9f0"]
}, Open ]],
Cell[BoxData[
RowBox[{"(*",
RowBox[{
RowBox[{
"Together", " ", "with", " ", "the", " ", "dtu", " ", "factor", " ", "the",
" ", "above", " ", "expression", " ", "becomes", " ", "this"}], " ", "-",
" ",
RowBox[{"now", " ", "set", " ", "the", " ",
RowBox[{"numerator", "'"}], "s", " ", "u", " ", "coefficient", " ", "to",
" ", "0"}]}], "*)"}]], "Input",
CellChangeTimes->{{3.844827204397766*^9, 3.844827227791861*^9}, {
3.8448280913787327`*^9, 3.844828116350154*^9}},
CellLabel->"In[7]:=",ExpressionUUID->"487b0363-9a4b-4f76-a23a-7bedb669b210"],
Cell[BoxData[
RowBox[{
RowBox[{
FractionBox[
SqrtBox[
RowBox[{"1", "+",
RowBox[{"A", " ", "X"}], "+",
RowBox[{
SuperscriptBox["A", "2"], " ", "Y"}], "+",
RowBox[{
SuperscriptBox["u", "2"], " ",
RowBox[{"(",
RowBox[{
SuperscriptBox["A", "2"], "+",
RowBox[{"A", " ", "X"}], "+", "Y"}], ")"}]}], "+",
RowBox[{"u", " ",
RowBox[{"(",
RowBox[{
RowBox[{"2", " ", "A"}], "+", "X", "+",
RowBox[{
SuperscriptBox["A", "2"], " ", "X"}], "+",
RowBox[{"2", " ", "A", " ", "Y"}]}], ")"}]}]}]],
RowBox[{"Sqrt", "[",
RowBox[{"1", "-",
RowBox[{"u", "^", "2"}]}], "]"}]],
FractionBox[
RowBox[{"Sqrt", "[",
RowBox[{
RowBox[{"A", "^", "2"}], "-", "1"}], "]"}],
RowBox[{
RowBox[{"(",
RowBox[{"A", "+", "u"}], ")"}], "^", "2"}]]}], ";"}]], "Input",
CellChangeTimes->{{3.844827071593066*^9, 3.84482720133156*^9}, {
3.844827237467803*^9, 3.844827245715302*^9}, {3.844827314417409*^9,
3.8448273411209717`*^9}, {3.8448275115380373`*^9, 3.844827512120165*^9}},
CellLabel->"In[8]:=",ExpressionUUID->"fd06d2f8-0b02-4d43-b5a6-4892b889e532"],
Cell[BoxData[
RowBox[{
RowBox[{"tfint", "[",
RowBox[{"X_", ",", "Y_"}], "]"}], ":=",
RowBox[{"Module", "[",
RowBox[{
RowBox[{"{",
RowBox[{"A", ",", "F"}], "}"}], ",",
RowBox[{
RowBox[{"A", "=",
RowBox[{"First", "[",
RowBox[{"Select", "[",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"Root", "[",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"X",
RowBox[{"#", "^", "2"}]}], "+",
RowBox[{"2",
RowBox[{"(",
RowBox[{"1", "+", "Y"}], ")"}], "#"}], "+", "X"}], "&"}],
")"}], ",", "#"}], "]"}], "&"}], "/@",
RowBox[{"{",
RowBox[{"1", ",", "2"}], "}"}]}], ",",
RowBox[{
RowBox[{
RowBox[{"Abs", "[", "#", "]"}], ">", "1"}], "&"}]}], "]"}],
"]"}]}], ";",
RowBox[{"F", "=",
RowBox[{
RowBox[{"A", "^", "2"}], "+",
RowBox[{"A", " ", "X"}], "+", "Y"}]}], ";",
RowBox[{"{",
RowBox[{"A", ",",
RowBox[{"F",
RowBox[{"(",
RowBox[{
RowBox[{"A", "^", "2"}], "-", "1"}], ")"}]}], ",",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"A", "^", "2"}], "Y"}], "+",
RowBox[{"A", " ", "X"}], "+", "1"}], ")"}], "/", "F"}]}], "}"}]}]}],
"]"}]}]], "Input",
CellChangeTimes->{{3.844827247592165*^9, 3.844827267949626*^9}, {
3.844827499380052*^9, 3.844827580797852*^9}, {3.844827612945922*^9,
3.844827722010957*^9}, {3.844827833613723*^9, 3.8448278999417353`*^9}, {
3.844827936999802*^9, 3.844828077101706*^9}, {3.84482816619991*^9,
3.8448283223735723`*^9}, {3.844828448696837*^9, 3.844828448905661*^9}},
CellLabel->"In[9]:=",ExpressionUUID->"17f0ba9b-f3fc-4828-a621-b61d87a8440b"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{"{",
RowBox[{"A0", ",", "B0", ",", "C0"}], "}"}], "=",
RowBox[{"ToNumberField", "[",
RowBox[{"tfint", "[",
RowBox[{
RowBox[{"8", "/", "9"}], ",",
RowBox[{"17", "/", "9"}]}], "]"}], "]"}]}]], "Input",
CellChangeTimes->{{3.8448277064991302`*^9, 3.844827714270494*^9}, {
3.844827947748959*^9, 3.8448279491374187`*^9}, {3.844828336009872*^9,
3.84482835063274*^9}, {3.844828463027835*^9, 3.844828463980035*^9}},
CellLabel->"In[10]:=",ExpressionUUID->"71ff88a2-acd6-4717-b786-efb9232c7818"],
Cell[BoxData[
RowBox[{"{",
RowBox[{
InterpretationBox[
TemplateBox[{"AlgebraicNumber",
InterpretationBox[
StyleBox[
TemplateBox[{"\"-6.34\"",
DynamicBox[
FEPrivate`FrontEndResource[
"FEExpressions", "NumericalApproximationElider"]]}, "RowDefault"],
ShowStringCharacters ->
False], -6.34232921921324521719043332268483936787`15.954589770191003,
Editable -> False],
TagBox[
RowBox[{"AlgebraicNumber", "[",
RowBox[{
SqrtBox["17"], ",",
RowBox[{"{",
RowBox[{
RowBox[{"-",
FractionBox["13", "4"]}], ",",
RowBox[{"-",
FractionBox["3", "4"]}]}], "}"}]}], "]"}],
Short[#, 7]& ], -6.342329219213245},
"NumericalApproximation"],
AlgebraicNumber[17^Rational[1, 2], {
Rational[-13, 4],
Rational[-3, 4]}]], ",",
InterpretationBox[
TemplateBox[{"AlgebraicNumber",
InterpretationBox[
StyleBox[
TemplateBox[{"\"1.43\"",
DynamicBox[
FEPrivate`FrontEndResource[
"FEExpressions", "NumericalApproximationElider"]],
"\"\[ThinSpace]\[Times]\[ThinSpace]\"",
TemplateBox[{"\"10\"", "\"3\""}, "Superscript", SyntaxForm ->
SuperscriptBox]}, "RowDefault"], ShowStringCharacters -> False],
1430.7920054002679535187780857086181640625`15.954589770191003,
Editable -> False],
TagBox[
RowBox[{"AlgebraicNumber", "[",
RowBox[{
SqrtBox["17"], ",",
RowBox[{"{",
RowBox[{
FractionBox["22865", "32"], ",",
FractionBox["5559", "32"]}], "}"}]}], "]"}], Short[#, 7]& ],
1430.792005400268},
"NumericalApproximation"],
AlgebraicNumber[17^Rational[1, 2], {
Rational[22865, 32],
Rational[5559, 32]}]], ",",
InterpretationBox[
TemplateBox[{"AlgebraicNumber",
InterpretationBox[
StyleBox[
TemplateBox[{"\"1.96\"",
DynamicBox[
FEPrivate`FrontEndResource[
"FEExpressions", "NumericalApproximationElider"]]}, "RowDefault"],
ShowStringCharacters -> False],
1.95587251835637854213700848049484193325`15.954589770191003, Editable ->
False],
TagBox[
RowBox[{"AlgebraicNumber", "[",
RowBox[{
SqrtBox["17"], ",",
RowBox[{"{",
RowBox[{
FractionBox["169", "137"], ",",
FractionBox["24", "137"]}], "}"}]}], "]"}], Short[#, 7]& ],
1.9558725183563785`},
"NumericalApproximation"],
AlgebraicNumber[17^Rational[1, 2], {
Rational[169, 137],
Rational[24, 137]}]]}], "}"}]], "Output",
CellChangeTimes->{{3.844827714655978*^9, 3.84482772405648*^9}, {
3.844827887514481*^9, 3.844827903137574*^9}, {3.844827943640443*^9,
3.844827949394483*^9}, {3.844828007569727*^9, 3.844828024693095*^9}, {
3.844828330650276*^9, 3.84482835090795*^9}, 3.844828464746162*^9,
3.8448286766698017`*^9, 3.8448606857171707`*^9, 3.84486670002454*^9,
3.844869650208859*^9, 3.844871910237872*^9, 3.844871974222797*^9,
3.8449025035007477`*^9, 3.8449026148063498`*^9, 3.8449030330130577`*^9},
CellLabel->"Out[10]=",ExpressionUUID->"dbdb8512-329e-4881-b12d-927b6e332805"]
}, Open ]],
Cell[BoxData[
RowBox[{
RowBox[{"I2", "[", "y_", "]"}], ":=",
RowBox[{
RowBox[{"Sqrt", "[", "B0", "]"}],
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"-", "A0"}], " ",
RowBox[{"NIntegrate", "[",
RowBox[{
RowBox[{
FractionBox[
RowBox[{"Sqrt", "[",
RowBox[{
RowBox[{"u", "^", "2"}], "+", "C0"}], "]"}],
RowBox[{"Sqrt", "[",
RowBox[{"1", "-",
RowBox[{"u", "^", "2"}]}], "]"}]],
FractionBox[
RowBox[{"2", "u"}],
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"u", "^", "2"}], "-",
RowBox[{"A0", "^", "2"}]}], ")"}], "^", "2"}]]}], ",",
RowBox[{"{",
RowBox[{"u", ",",
FractionBox[
RowBox[{
RowBox[{"A0", " ",
RowBox[{"Cos", "[",
RowBox[{"3", "y"}], "]"}]}], "-", "1"}],
RowBox[{"A0", "-",
RowBox[{"Cos", "[",
RowBox[{"3", "y"}], "]"}]}]], ",", "1"}], "}"}], ",",
RowBox[{"WorkingPrecision", "\[Rule]", "50"}]}], "]"}]}], "+",
RowBox[{"NIntegrate", "[",
RowBox[{
RowBox[{
FractionBox[
RowBox[{"Sqrt", "[",
RowBox[{
RowBox[{"u", "^", "2"}], "+", "C0"}], "]"}],
RowBox[{"Sqrt", "[",
RowBox[{"1", "-",
RowBox[{"u", "^", "2"}]}], "]"}]],
FractionBox[
RowBox[{
RowBox[{"u", "^", "2"}], "+",
RowBox[{"A0", "^", "2"}]}],
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"u", "^", "2"}], "-",
RowBox[{"A0", "^", "2"}]}], ")"}], "^", "2"}]]}], ",",
RowBox[{"{",
RowBox[{"u", ",",
FractionBox[
RowBox[{
RowBox[{"A0", " ",
RowBox[{"Cos", "[",
RowBox[{"3", "y"}], "]"}]}], "-", "1"}],
RowBox[{"A0", "-",
RowBox[{"Cos", "[",
RowBox[{"3", "y"}], "]"}]}]], ",", "1"}], "}"}], ",",
RowBox[{"WorkingPrecision", "\[Rule]", "50"}]}], "]"}]}],
")"}]}]}]], "Input",
CellChangeTimes->{{3.844821443064513*^9, 3.844821446894187*^9}, {
3.844821480266576*^9, 3.844821513799509*^9}, {3.844821554077676*^9,
3.8448215544665337`*^9}, {3.844821693650407*^9, 3.8448217229839373`*^9}, {
3.8448217570166607`*^9, 3.844821775827737*^9}, {3.8448222430120792`*^9,
3.844822244316193*^9}, {3.84482229385002*^9, 3.8448223157861357`*^9}, {
3.844822429714842*^9, 3.844822430020014*^9}, {3.844822762505087*^9,
3.844822931529119*^9}, {3.844822973849169*^9, 3.844823015142288*^9}, {
3.844823306300263*^9, 3.8448233746470213`*^9}, {3.844823608177314*^9,
3.844823611043292*^9}, {3.844823682687092*^9, 3.844823737684701*^9}, {
3.844823821101522*^9, 3.844823891826459*^9}, {3.844823950520205*^9,
3.844823977957486*^9}, {3.844824019357082*^9, 3.8448240229601517`*^9}, {
3.844824096808169*^9, 3.844824139199308*^9}, {3.84482423548951*^9,
3.8448242414713573`*^9}, {3.844824336988799*^9, 3.84482440383946*^9}, {
3.844824433917536*^9, 3.84482445798846*^9}, {3.844824573843403*^9,
3.844824689954327*^9}, {3.844825399786707*^9, 3.844825405060638*^9}, {
3.8448261748523617`*^9, 3.8448262096122293`*^9}, {3.844826577213462*^9,
3.844826635592266*^9}, {3.8448266874845667`*^9, 3.844826695695825*^9}, {
3.844828363679172*^9, 3.844828440394864*^9}, {3.844828476014921*^9,
3.8448284788629627`*^9}, {3.8448285568605633`*^9, 3.844828569943495*^9}, {
3.8448286638673363`*^9, 3.844828666251481*^9}, {3.844828746512411*^9,
3.844828749929131*^9}, 3.8448603808250237`*^9, {3.844860689727004*^9,
3.844860691604583*^9}, {3.8448607408040657`*^9, 3.844860807279091*^9},
3.844860893551258*^9, {3.844860929731268*^9, 3.844860936770464*^9}, {
3.84486122394806*^9, 3.844861264182179*^9}, {3.844861294566547*^9,
3.844861333861546*^9}, {3.8448613657252407`*^9, 3.844861446666098*^9}},
CellLabel->"In[11]:=",ExpressionUUID->"b7f4a108-0135-4c60-aa3f-b15abe8d976d"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"I2", "[", "1", "]"}]], "Input",
CellChangeTimes->{{3.844826499413856*^9, 3.8448265003098507`*^9}, {
3.844826638720771*^9, 3.844826640583702*^9}, {3.844860754555237*^9,
3.844860755246234*^9}, {3.844865731720574*^9, 3.8448657328360643`*^9},
3.8448693664531307`*^9},
CellLabel->"In[12]:=",ExpressionUUID->"ef3517f4-ce08-4210-94b7-41051fd5e0e7"],
Cell[BoxData["4.\
60432427101551068986302631297185217317292643002711640475929259431495036323599`\
50."], "Output",
CellChangeTimes->{
3.844826640853223*^9, 3.84482669838761*^9, {3.844828414686442*^9,
3.8448284421153593`*^9}, 3.84482848069074*^9, 3.844828585364065*^9,
3.844828676741199*^9, 3.844860708644762*^9, 3.844860755503331*^9,
3.844860812021942*^9, 3.8448609069992113`*^9, 3.844861266239799*^9,
3.844861448982236*^9, 3.844865733366534*^9, 3.8448667002294407`*^9,
3.84486936697433*^9, 3.8448696504544163`*^9, 3.8448719109489594`*^9,
3.844871974794923*^9, 3.8449026153747787`*^9, 3.844903034452674*^9},
CellLabel->"Out[12]=",ExpressionUUID->"167e9b9d-c045-45dc-8db4-0fe4f8a1a011"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"alg0", "=",
RowBox[{
FractionBox["2", "3"],
RowBox[{"(",
RowBox[{
FractionBox[
SqrtBox[
RowBox[{
RowBox[{"(",
RowBox[{"1", "-",
RowBox[{"u", "^", "2"}]}], ")"}], " ",
RowBox[{"(",
RowBox[{
RowBox[{"9",
RowBox[{"u", "^", "2"}]}], "+",
RowBox[{"8", "u"}], "+", "17"}], ")"}]}]],
RowBox[{"13", "+",
RowBox[{"4", " ", "u"}]}]], "+",
RowBox[{
FractionBox["2", "3"],
RowBox[{"ArcTan", "[",
SqrtBox[
FractionBox[
RowBox[{"1", "-",
RowBox[{"u", "^", "2"}]}],
RowBox[{
RowBox[{"u", "^", "2"}], "+",
RowBox[{
RowBox[{"(",
RowBox[{"8", "/", "9"}], ")"}], "u"}], "+",
RowBox[{"(",
RowBox[{"17", "/", "9"}], ")"}]}]]], "]"}]}]}], ")"}]}]}]], "Input",
CellChangeTimes->{{3.844864152931322*^9, 3.844864205101933*^9}, {
3.844864248038286*^9, 3.844864251123369*^9}, {3.844864341837607*^9,
3.844864384716866*^9}, {3.844864538531673*^9, 3.8448645390568867`*^9}, {
3.844864818802751*^9, 3.844864844103869*^9}, 3.844869358920162*^9},
CellLabel->"In[13]:=",ExpressionUUID->"bdbd3a53-e9f5-4340-a0ab-70c4409eb89a"],
Cell[BoxData[
RowBox[{
FractionBox["2", "3"], " ",
RowBox[{"(",
RowBox[{
FractionBox[
SqrtBox[
RowBox[{
RowBox[{"(",
RowBox[{"1", "-",
SuperscriptBox["u", "2"]}], ")"}], " ",
RowBox[{"(",
RowBox[{"17", "+",
RowBox[{"8", " ", "u"}], "+",
RowBox[{"9", " ",
SuperscriptBox["u", "2"]}]}], ")"}]}]],
RowBox[{"13", "+",
RowBox[{"4", " ", "u"}]}]], "+",
RowBox[{
FractionBox["2", "3"], " ",
RowBox[{"ArcTan", "[",
SqrtBox[
FractionBox[
RowBox[{"1", "-",
SuperscriptBox["u", "2"]}],
RowBox[{
FractionBox["17", "9"], "+",
FractionBox[
RowBox[{"8", " ", "u"}], "9"], "+",
SuperscriptBox["u", "2"]}]]], "]"}]}]}], ")"}]}]], "Output",
CellChangeTimes->{
3.844864385052133*^9, 3.8448645413644743`*^9, {3.8448648232997723`*^9,
3.84486487254558*^9}, 3.8448667002432737`*^9, 3.844869362228228*^9,
3.8448696504684877`*^9, 3.84487191099831*^9, 3.844871974822027*^9,
3.844902615404888*^9, 3.844903034892516*^9},
CellLabel->"Out[13]=",ExpressionUUID->"7423d711-af75-4537-ab65-3a2fe273ea96"]
}, Open ]],
Cell[BoxData[
RowBox[{
RowBox[{"I3", "[", "y_", "]"}], ":=",
RowBox[{
RowBox[{"(",
RowBox[{"alg0", "/.",
RowBox[{"u", "->",
RowBox[{"Cos", "[",
RowBox[{"3", "y"}], "]"}]}]}], ")"}], "+",
RowBox[{
RowBox[{"Sqrt", "[", "B0", "]"}],
RowBox[{"NIntegrate", "[",
RowBox[{
RowBox[{
FractionBox["1",
RowBox[{
RowBox[{"Sqrt", "[",
RowBox[{
RowBox[{"u", "^", "2"}], "+", "C0"}], "]"}],
RowBox[{"Sqrt", "[",
RowBox[{"1", "-",
RowBox[{"u", "^", "2"}]}], "]"}]}]],
FractionBox[
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"u", "^", "2"}], "+",
RowBox[{"A0", "^", "2"}]}], ")"}],
RowBox[{"(",
RowBox[{
RowBox[{"u", "^", "2"}], "+", "C0"}], ")"}]}],
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"u", "^", "2"}], "-",
RowBox[{"A0", "^", "2"}]}], ")"}], "^", "2"}]]}], ",",
RowBox[{"{",
RowBox[{"u", ",",
FractionBox[
RowBox[{
RowBox[{"A0", " ",
RowBox[{"Cos", "[",
RowBox[{"3", "y"}], "]"}]}], "-", "1"}],
RowBox[{"A0", "-",
RowBox[{"Cos", "[",
RowBox[{"3", "y"}], "]"}]}]], ",", "1"}], "}"}], ",",
RowBox[{"WorkingPrecision", "\[Rule]", "50"}]}], "]"}]}]}]}]], "Input",
CellChangeTimes->{
3.844861470998014*^9, {3.844864618023984*^9, 3.844864633534925*^9}, {
3.844864706169724*^9, 3.8448647066461163`*^9}, {3.844864880766942*^9,
3.844864892192451*^9}, {3.844865183412695*^9, 3.8448652054471693`*^9}, {
3.8448652913312798`*^9, 3.8448652963251657`*^9}, 3.8448693606841307`*^9},
CellLabel->"In[14]:=",ExpressionUUID->"299d4625-771d-4a6a-884b-2f79fb7634c3"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"I3", "[", "1", "]"}]], "Input",
CellChangeTimes->{{3.844864634992178*^9, 3.8448646354168873`*^9}, {
3.8448647092439938`*^9, 3.844864723149314*^9}, {3.844864894792757*^9,
3.844864896419198*^9}, {3.844865680112934*^9, 3.84486568076863*^9}, {
3.8448657380523863`*^9, 3.8448657388738527`*^9}, {3.844866113376225*^9,
3.844866113945669*^9}},
CellLabel->"In[15]:=",ExpressionUUID->"e1ca574f-fa9b-4cc2-9f2f-82520996ed5b"],
Cell[BoxData["4.\
60432427101551068986302631297185217317292643002711879708452797852917578369075`\
50.00841506199381"], "Output",
CellChangeTimes->{{3.844864709792741*^9, 3.844864723652151*^9},
3.844864896773881*^9, 3.844865209207863*^9, 3.844865681082883*^9,
3.8448657392444563`*^9, 3.844866114297512*^9, 3.844866700386488*^9,
3.844869368658635*^9, 3.844869650610695*^9, 3.8448719113959513`*^9,
3.84487197511795*^9, 3.8449026157722883`*^9, 3.844903035574584*^9},
CellLabel->"Out[15]=",ExpressionUUID->"df6a278e-ca16-49d1-b9ed-479a03ca9c61"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"g", "=",
RowBox[{"Sqrt", "[",
RowBox[{"B0", "/",
RowBox[{"(",
RowBox[{"1", "+", "C0"}], ")"}]}], "]"}]}]], "Input",
CellChangeTimes->{{3.844866717973673*^9, 3.844866722322034*^9}, {
3.844866752500511*^9, 3.8448667622264023`*^9}},
CellLabel->"In[16]:=",ExpressionUUID->"bf1012a0-2be9-4b29-815d-104413e2177c"],
Cell[BoxData[
SqrtBox[
InterpretationBox[
TemplateBox[{"AlgebraicNumber",
InterpretationBox[
StyleBox[
TemplateBox[{"\"484.\"",
DynamicBox[
FEPrivate`FrontEndResource[
"FEExpressions", "NumericalApproximationElider"]]}, "RowDefault"],
ShowStringCharacters -> False],
484.05064714897247313274419866502285003662`15.954589770191003, Editable ->
False],
TagBox[
RowBox[{"AlgebraicNumber", "[",
RowBox[{
SqrtBox["17"], ",",
RowBox[{"{",
RowBox[{
FractionBox["15453", "64"], ",",
FractionBox["11297", "192"]}], "}"}]}], "]"}], Short[#, 7]& ],
484.0506471489725},
"NumericalApproximation"],
AlgebraicNumber[17^Rational[1, 2], {
Rational[15453, 64],
Rational[11297, 192]}]]]], "Output",
CellChangeTimes->{
3.8448667625572023`*^9, 3.844869650717751*^9, 3.8448719115032167`*^9,
3.8448719752876453`*^9, {3.8449024847476187`*^9, 3.844902507249858*^9},
3.8449026159132*^9, 3.844903035865369*^9},
CellLabel->"Out[16]=",ExpressionUUID->"162606d8-3ff7-4f3a-b398-87992b2cfc20"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"m", "=",
RowBox[{"1", "/",
RowBox[{"(",
RowBox[{"C0", "+", "1"}], ")"}]}]}]], "Input",
CellChangeTimes->{{3.844865516243475*^9, 3.844865524311955*^9}},
CellLabel->"In[17]:=",ExpressionUUID->"005d25f4-570c-4096-842a-ce05c17c86ae"],
Cell[BoxData[
InterpretationBox[
TemplateBox[{"AlgebraicNumber",
InterpretationBox[
StyleBox[
TemplateBox[{"\"0.338\"",
DynamicBox[
FEPrivate`FrontEndResource[
"FEExpressions", "NumericalApproximationElider"]]}, "RowDefault"],
ShowStringCharacters -> False],
0.33830958330911137244001452017982956022`15.954589770191003, Editable ->
False],
TagBox[
RowBox[{"AlgebraicNumber", "[",
RowBox[{
SqrtBox["17"], ",",
RowBox[{"{",
RowBox[{
FractionBox["1", "2"], ",",
RowBox[{"-",
FractionBox["2", "51"]}]}], "}"}]}], "]"}], Short[#, 7]& ],
0.3383095833091114},
"NumericalApproximation"],
AlgebraicNumber[17^Rational[1, 2], {
Rational[1, 2],
Rational[-2, 51]}]]], "Output",
CellChangeTimes->{
3.844866700482345*^9, 3.844869650743247*^9, 3.844871911550063*^9,
3.8448719753350153`*^9, {3.844902485692546*^9, 3.844902508666712*^9},
3.8449026160795527`*^9, 3.844903036078229*^9},
CellLabel->"Out[17]=",ExpressionUUID->"b052984b-2f52-4d1d-a7f4-8f89dd240263"]
}, Open ]],
Cell[BoxData[
RowBox[{
RowBox[{"I4", "[", "y_", "]"}], ":=",
RowBox[{
RowBox[{"(",
RowBox[{"alg0", "/.",
RowBox[{"u", "->",
RowBox[{"Cos", "[",
RowBox[{"3", "y"}], "]"}]}]}], ")"}], "+",
RowBox[{"g", " ",
RowBox[{"NIntegrate", "[",
RowBox[{
FractionBox[
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"JacobiCN", "[",
RowBox[{"u", ",", "m"}], "]"}], "^", "2"}], "+",
RowBox[{"A0", "^", "2"}]}], ")"}],
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"JacobiCN", "[",
RowBox[{"u", ",", "m"}], "]"}], "^", "2"}], "+", "C0"}], ")"}]}],
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"JacobiCN", "[",
RowBox[{"u", ",", "m"}], "]"}], "^", "2"}], "-",
RowBox[{"A0", "^", "2"}]}], ")"}], "^", "2"}]], ",",
RowBox[{"{",
RowBox[{"u", ",", "0", ",",
RowBox[{"EllipticF", "[",
RowBox[{
RowBox[{"ArcCos", "[",
FractionBox[
RowBox[{
RowBox[{"A0", " ",
RowBox[{"Cos", "[",
RowBox[{"3", "y"}], "]"}]}], "-", "1"}],
RowBox[{"A0", "-",
RowBox[{"Cos", "[",
RowBox[{"3", "y"}], "]"}]}]], "]"}], ",", "m"}], "]"}]}], "}"}],
",",
RowBox[{"WorkingPrecision", "\[Rule]", "50"}]}], "]"}]}]}]}]], "Input",
CellChangeTimes->{
3.844861470998014*^9, {3.844864618023984*^9, 3.844864633534925*^9}, {
3.844864706169724*^9, 3.8448647066461163`*^9}, {3.844864880766942*^9,
3.844864892192451*^9}, {3.844865183412695*^9, 3.8448652054471693`*^9}, {
3.8448652913312798`*^9, 3.8448652963251657`*^9}, {3.844865479676083*^9,
3.84486548915809*^9}, {3.844865543303809*^9, 3.844865591469469*^9}, {
3.844865641500799*^9, 3.844865648990294*^9}, {3.844865749267824*^9,
3.8448657542303963`*^9}, {3.844865784560482*^9, 3.844865786246904*^9},
3.844869407221527*^9},
CellLabel->"In[18]:=",ExpressionUUID->"3cef2ab3-6910-4a01-bc70-fbf276c56415"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"I4", "[", "1", "]"}]], "Input",
CellChangeTimes->{{3.8448661077003517`*^9, 3.8448661100408173`*^9}},
CellLabel->"In[19]:=",ExpressionUUID->"2cc6d49e-094d-46db-ab40-e274ab172252"],
Cell[BoxData["4.\
60432427101551068986302631297185217317292643002712657214157619838415122826182`\
50.00841506199381"], "Output",
CellChangeTimes->{3.844866110425797*^9, 3.844866142839147*^9,
3.8448663133545303`*^9, 3.844866700581058*^9, 3.844866771506757*^9,
3.8448694092198544`*^9, 3.844869650858329*^9, 3.844871911782124*^9,
3.8448719755936613`*^9, 3.8449026162775507`*^9, 3.844903036448471*^9},
CellLabel->"Out[19]=",ExpressionUUID->"266e04d5-b873-4b21-b94f-e2d6f9bac4de"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"n", "=",
InterpretationBox[
TemplateBox[{"AlgebraicNumber",
InterpretationBox[
StyleBox[
TemplateBox[{"\"-0.0255\"",
DynamicBox[
FEPrivate`FrontEndResource[
"FEExpressions", "NumericalApproximationElider"]]}, "RowDefault"],
ShowStringCharacters ->
False], -0.02549385424538811242833880044145189459`15.954589770191003,
Editable -> False],
TagBox[
RowBox[{"AlgebraicNumber", "[",
RowBox[{
SqrtBox["17"], ",",
RowBox[{"{",
RowBox[{
FractionBox["1", "2"], ",",
RowBox[{"-",
FractionBox["13", "102"]}]}], "}"}]}], "]"}],
Short[#, 7]& ], -0.025493854245388112`},
"NumericalApproximation"],
AlgebraicNumber[17^Rational[1, 2], {
Rational[1, 2],
Rational[-13, 102]}]]}]], "Input",
CellChangeTimes->{{3.8448668402426434`*^9, 3.84486684076094*^9}},
CellLabel->"In[20]:=",ExpressionUUID->"914d2159-72c9-4d2e-91f5-66d003149793"],
Cell[BoxData[
InterpretationBox[
TemplateBox[{"AlgebraicNumber",
InterpretationBox[
StyleBox[
TemplateBox[{"\"-0.0255\"",
DynamicBox[
FEPrivate`FrontEndResource[
"FEExpressions", "NumericalApproximationElider"]]}, "RowDefault"],
ShowStringCharacters ->
False], -0.02549385424538811242833880044145189459`15.954589770191003,
Editable -> False],
TagBox[
RowBox[{"AlgebraicNumber", "[",
RowBox[{
SqrtBox["17"], ",",
RowBox[{"{",
RowBox[{
FractionBox["1", "2"], ",",
RowBox[{"-",
FractionBox["13", "102"]}]}], "}"}]}], "]"}],
Short[#, 7]& ], -0.025493854245388112`},
"NumericalApproximation"],
AlgebraicNumber[17^Rational[1, 2], {
Rational[1, 2],
Rational[-13, 102]}]]], "Output",
CellChangeTimes->{3.844866842463827*^9, 3.844869650874164*^9,
3.8448719118120317`*^9, 3.844871975620245*^9, 3.8449024906178637`*^9,
3.844902616309977*^9, 3.8449030365719137`*^9},
CellLabel->"Out[20]=",ExpressionUUID->"117492f6-aca7-4ce5-8f5e-44c16f57cc41"]
}, Open ]],
Cell[BoxData[
RowBox[{
RowBox[{"alg", "=",
RowBox[{
FractionBox["1", "3"],
RowBox[{"(",
RowBox[{
FractionBox[
SqrtBox[
RowBox[{
RowBox[{"(",
RowBox[{"1", "-",
RowBox[{"u", "^", "2"}]}], ")"}], " ",
RowBox[{"(",
RowBox[{
RowBox[{"9",
RowBox[{"u", "^", "2"}]}], "+",
RowBox[{"8", "u"}], "+", "17"}], ")"}]}]],
RowBox[{"u", "-", "A0"}]], "+",
RowBox[{
FractionBox["4", "3"], " ",
RowBox[{"ArcTan", "[",
SqrtBox[
FractionBox[
RowBox[{"1", "-",
RowBox[{"u", "^", "2"}]}],
RowBox[{
RowBox[{"u", "^", "2"}], "+",
RowBox[{
RowBox[{"(",
RowBox[{"8", "/", "9"}], ")"}], "u"}], "+",
RowBox[{"(",
RowBox[{"17", "/", "9"}], ")"}]}]]], "]"}]}]}], ")"}]}]}], ";"}]],\
"Input",
CellChangeTimes->{{3.8448688079809*^9, 3.844868867746518*^9}, {
3.8448689021434383`*^9, 3.8448689030600233`*^9}, {3.844868940338352*^9,
3.844868995774519*^9}, {3.844869281912067*^9, 3.844869315593253*^9},
3.8448694892803926`*^9},
CellLabel->"In[21]:=",ExpressionUUID->"d6140dfb-24c7-4947-a556-27087f8651b5"],
Cell[BoxData[
RowBox[{
RowBox[{"I5", "[", "y_", "]"}], ":=",
RowBox[{
RowBox[{"With", "[",
RowBox[{
RowBox[{"{",
RowBox[{"\[Phi]", "=",
RowBox[{"ArcCos", "[",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"A0", " ", "u"}], "-", "1"}], ")"}], "/",
RowBox[{"(",
RowBox[{"A0", "-", "u"}], ")"}]}], "]"}]}], "}"}], ",",
RowBox[{"alg", "+",
RowBox[{"g",
RowBox[{"(",
RowBox[{
RowBox[{"n",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"EllipticF", "[",
RowBox[{"\[Phi]", ",", "m"}], "]"}], "-", " ",
RowBox[{"EllipticE", "[",
RowBox[{"\[Phi]", ",", "m"}], "]"}]}], ")"}], "/", "m"}]}], "+",
RowBox[{
RowBox[{"(",
RowBox[{"16", "/", "137"}], ")"}],
RowBox[{"(",
RowBox[{"1", "-", "m"}], ")"}], " ",
RowBox[{"EllipticPi", "[",
RowBox[{"n", ",", "\[Phi]", ",", "m"}], "]"}]}]}], ")"}]}]}]}],
"]"}], "/.",
RowBox[{"u", "->",
RowBox[{"Cos", "[",
RowBox[{"3", "y"}], "]"}]}]}]}]], "Input",
CellChangeTimes->{
3.844861470998014*^9, {3.844864618023984*^9, 3.844864633534925*^9}, {
3.844864706169724*^9, 3.8448647066461163`*^9}, {3.844864880766942*^9,
3.844864892192451*^9}, {3.844865183412695*^9, 3.8448652054471693`*^9}, {
3.8448652913312798`*^9, 3.8448652963251657`*^9}, {3.844865479676083*^9,
3.84486548915809*^9}, {3.844865543303809*^9, 3.844865591469469*^9}, {
3.844865641500799*^9, 3.844865648990294*^9}, {3.844865749267824*^9,
3.8448657542303963`*^9}, {3.844865784560482*^9, 3.844865786246904*^9}, {
3.844866136289607*^9, 3.844866303228567*^9}, {3.844866394989873*^9,
3.844866440476172*^9}, {3.84486648202693*^9, 3.844866611011279*^9}, {
3.844866679758285*^9, 3.844866680672277*^9}, {3.8448669828535767`*^9,
3.844867007310596*^9}, {3.844867218441444*^9, 3.8448672778634787`*^9}, {
3.844867308477455*^9, 3.8448673134011497`*^9}, {3.844867405051311*^9,
3.8448675039839373`*^9}, {3.844867634052408*^9, 3.844867634152132*^9}, {
3.844867699748785*^9, 3.8448677103407183`*^9}, {3.844867752498082*^9,
3.844867762367639*^9}, {3.8448678696206827`*^9, 3.8448678715394783`*^9}, {
3.844868207677558*^9, 3.84486822015469*^9}, {3.844868319969082*^9,
3.844868383129423*^9}, {3.8448684456693277`*^9, 3.844868466900239*^9},
3.844868612049697*^9, {3.844868650156987*^9, 3.844868652126088*^9}, {
3.844868873530944*^9, 3.8448689001738167`*^9}, {3.844869259947414*^9,
3.844869262106567*^9}, {3.8448695126423798`*^9, 3.844869538222821*^9}, {
3.844869569524675*^9, 3.844869585721293*^9}, 3.8448696350069857`*^9, {
3.84486998057966*^9, 3.844870003485045*^9}, {3.844870560827886*^9,
3.844870584366131*^9}},
CellLabel->"In[22]:=",ExpressionUUID->"fc7d6d4a-cf1d-4446-889d-2aef239a90ea"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"N", "[",
RowBox[{
RowBox[{"I5", "[", "1", "]"}], ",", "50"}], "]"}]], "Input",
CellChangeTimes->{{3.8448714483434467`*^9, 3.8448714530983667`*^9}, {
3.844871930435587*^9, 3.8448719403299932`*^9}},
CellLabel->"In[23]:=",ExpressionUUID->"e249b5a5-a6ee-43e9-b638-462b0a2bd512"],
Cell[BoxData["4.\
60432427101551068986302631297185217317292643002712657214148941586318041214393`\
50."], "Output",
CellChangeTimes->{
3.844871453408757*^9, {3.844871912032625*^9, 3.8448719406720533`*^9},
3.8448719757512693`*^9, {3.844902477426375*^9, 3.844902515236661*^9},
3.844902616536702*^9, 3.844903037091456*^9},
CellLabel->"Out[23]=",ExpressionUUID->"8f0cedba-c5a1-4159-8b5f-ee8a62e1e86a"]
}, Open ]],
Cell[BoxData[
RowBox[{
RowBox[{"J0", "[", "y_", "]"}], ":=",
RowBox[{"NIntegrate", "[",
RowBox[{
RowBox[{"Sqrt", "[",
RowBox[{"25", "+",
RowBox[{"16",
RowBox[{"Cos", "[",
RowBox[{"3", "t"}], "]"}]}], "+",
RowBox[{"4",
RowBox[{
RowBox[{"Cos", "[",
RowBox[{"3", "t"}], "]"}], "^", "2"}]}]}], "]"}], ",",
RowBox[{"{",
RowBox[{"t", ",", "0", ",", "y"}], "}"}], ",",
RowBox[{"WorkingPrecision", "\[Rule]", "50"}]}], "]"}]}]], "Input",
CellChangeTimes->{
3.844898161252595*^9, {3.8448981931956472`*^9, 3.844898202188374*^9}, {
3.844898236658979*^9, 3.844898239053442*^9}},
CellLabel->"In[24]:=",ExpressionUUID->"f12a59db-60f6-4605-b5e1-f6d60953ad18"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"J0", "[", "1", "]"}]], "Input",
CellChangeTimes->{{3.8448982142934513`*^9, 3.8448982279695473`*^9}},
CellLabel->"In[25]:=",ExpressionUUID->"783402ff-f7a3-4338-9491-b1687227e07f"],
Cell[BoxData["5.\
14608510651588948875063444866083409185868746013231264360420446407869091989138`\
50."], "Output",
CellChangeTimes->{{3.8448982288472*^9, 3.8448982409800167`*^9},
3.8448986857388372`*^9, 3.844902616688202*^9, 3.844903037235675*^9},
CellLabel->"Out[25]=",ExpressionUUID->"68f0a497-0b68-45de-b032-5b56f0bdfda2"]
}, Open ]],
Cell[BoxData[
RowBox[{
RowBox[{"J1", "[", "y_", "]"}], ":=",
RowBox[{
FractionBox["2", "3"],
RowBox[{"NIntegrate", "[",
RowBox[{
FractionBox[
RowBox[{"Sqrt", "[",
RowBox[{
RowBox[{"u", "^", "2"}], "+",
RowBox[{"4", "u"}], "+",
RowBox[{"(",
RowBox[{"25", "/", "4"}], ")"}]}], "]"}],
RowBox[{"Sqrt", "[",
RowBox[{"1", "-",
RowBox[{"u", "^", "2"}]}], "]"}]], ",",
RowBox[{"{",
RowBox[{"u", ",",
RowBox[{"Cos", "[",
RowBox[{"3", "y"}], "]"}], ",", "1"}], "}"}], ",",
RowBox[{"WorkingPrecision", "\[Rule]", "50"}]}], "]"}]}]}]], "Input",
CellChangeTimes->{{3.844898271161471*^9, 3.8448983012473927`*^9}, {
3.844898336972846*^9, 3.844898454735454*^9}, {3.844898698205385*^9,
3.844898706441332*^9}},
CellLabel->"In[26]:=",ExpressionUUID->"e662bf28-549b-47ce-b970-b0032324e207"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"J1", "[", "1", "]"}]], "Input",
CellChangeTimes->{{3.844898303183839*^9, 3.844898303761319*^9}},
CellLabel->"In[27]:=",ExpressionUUID->"edfea8cf-26fb-42e6-9b29-1a7f993033b7"],
Cell[BoxData["5.\
14608510651588948875063444866083409185868746013230265670335372109642357150122`\
50."], "Output",
CellChangeTimes->{
3.8448983041370773`*^9, {3.844898338525722*^9, 3.844898390499963*^9},
3.8448984577927713`*^9, {3.84489867889356*^9, 3.844898709590659*^9},
3.8449026168715153`*^9, 3.844903037759486*^9},
CellLabel->"Out[27]=",ExpressionUUID->"7a5e3764-722f-48f0-a6d5-301cc5076147"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{"{",
RowBox[{"A1", ",", "B1", ",", "C1"}], "}"}], "=",
RowBox[{"ToNumberField", "[",
RowBox[{"tfint", "[",
RowBox[{"4", ",",
RowBox[{"25", "/", "4"}]}], "]"}], "]"}]}]], "Input",
CellChangeTimes->{{3.8448984968964577`*^9, 3.844898498530471*^9}, {
3.84489854067976*^9, 3.8448985424058733`*^9}},
CellLabel->"In[28]:=",ExpressionUUID->"c59fb042-bbab-44e6-be92-674aed8698c5"],
Cell[BoxData[
RowBox[{"{",
RowBox[{
InterpretationBox[
TemplateBox[{"AlgebraicNumber",
InterpretationBox[
StyleBox[
TemplateBox[{"\"-3.32\"",
DynamicBox[
FEPrivate`FrontEndResource[
"FEExpressions", "NumericalApproximationElider"]]}, "RowDefault"],
ShowStringCharacters ->
False], -3.32417332780597796215715789003297686577`15.954589770191003,
Editable -> False],
TagBox[
RowBox[{"AlgebraicNumber", "[",
RowBox[{
SqrtBox["65"], ",",
RowBox[{"{",
RowBox[{
RowBox[{"-",
FractionBox["29", "16"]}], ",",
RowBox[{"-",
FractionBox["3", "16"]}]}], "}"}]}], "]"}],
Short[#, 7]& ], -3.324173327805978},
"NumericalApproximation"],
AlgebraicNumber[65^Rational[1, 2], {
Rational[-29, 16],
Rational[-3, 16]}]], ",",
InterpretationBox[
TemplateBox[{"AlgebraicNumber",
InterpretationBox[
StyleBox[
TemplateBox[{"\"40.2\"",
DynamicBox[
FEPrivate`FrontEndResource[
"FEExpressions", "NumericalApproximationElider"]]}, "RowDefault"],
ShowStringCharacters -> False],
40.23503546477434156258823350071907043457`15.954589770191003, Editable ->
False],
TagBox[
RowBox[{"AlgebraicNumber", "[",
RowBox[{
SqrtBox["65"], ",",
RowBox[{"{",
RowBox[{
FractionBox["145665", "8192"], ",",
FractionBox["22815", "8192"]}], "}"}]}], "]"}], Short[#, 7]& ],
40.23503546477434},
"NumericalApproximation"],
AlgebraicNumber[65^Rational[1, 2], {
Rational[145665, 8192],
Rational[22815, 8192]}]], ",",
InterpretationBox[
TemplateBox[{"AlgebraicNumber",
InterpretationBox[
StyleBox[
TemplateBox[{"\"14.2\"",
DynamicBox[
FEPrivate`FrontEndResource[
"FEExpressions", "NumericalApproximationElider"]]}, "RowDefault"],
ShowStringCharacters -> False],
14.17947552976123049006673682015389204025`15.954589770191003, Editable ->
False],
TagBox[
RowBox[{"AlgebraicNumber", "[",
RowBox[{
SqrtBox["65"], ",",
RowBox[{"{",
RowBox[{
FractionBox["57", "8"], ",",
FractionBox["7", "8"]}], "}"}]}], "]"}], Short[#, 7]& ],
14.17947552976123},
"NumericalApproximation"],
AlgebraicNumber[65^Rational[1, 2], {
Rational[57, 8],
Rational[7, 8]}]]}], "}"}]], "Output",
CellChangeTimes->{3.8448984989655867`*^9, 3.844898542830295*^9,
3.8449026170760183`*^9, 3.844903038140359*^9},
CellLabel->"Out[28]=",ExpressionUUID->"b02477d3-3227-4171-b7bc-b6c19809dbfe"]
}, Open ]],
Cell[BoxData[
RowBox[{
RowBox[{"J2", "[", "y_", "]"}], ":=",
RowBox[{
RowBox[{"(",
RowBox[{"2",
RowBox[{
RowBox[{"Sqrt", "[", "B1", "]"}], "/", "3"}]}], ")"}],
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"-", "A1"}], " ",
RowBox[{"NIntegrate", "[",
RowBox[{
RowBox[{
FractionBox[
RowBox[{"Sqrt", "[",
RowBox[{
RowBox[{"u", "^", "2"}], "+", "C1"}], "]"}],
RowBox[{"Sqrt", "[",
RowBox[{"1", "-",
RowBox[{"u", "^", "2"}]}], "]"}]],
FractionBox[
RowBox[{"2", "u"}],
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"u", "^", "2"}], "-",
RowBox[{"A1", "^", "2"}]}], ")"}], "^", "2"}]]}], ",",
RowBox[{"{",
RowBox[{"u", ",",
FractionBox[
RowBox[{
RowBox[{"A1", " ",
RowBox[{"Cos", "[",
RowBox[{"3", "y"}], "]"}]}], "-", "1"}],
RowBox[{"A1", "-",
RowBox[{"Cos", "[",
RowBox[{"3", "y"}], "]"}]}]], ",", "1"}], "}"}], ",",
RowBox[{"WorkingPrecision", "\[Rule]", "50"}]}], "]"}]}], "+",
RowBox[{"NIntegrate", "[",
RowBox[{
RowBox[{
FractionBox[
RowBox[{"Sqrt", "[",
RowBox[{
RowBox[{"u", "^", "2"}], "+", "C1"}], "]"}],
RowBox[{"Sqrt", "[",
RowBox[{"1", "-",
RowBox[{"u", "^", "2"}]}], "]"}]],
FractionBox[
RowBox[{
RowBox[{"u", "^", "2"}], "+",
RowBox[{"A1", "^", "2"}]}],
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"u", "^", "2"}], "-",
RowBox[{"A1", "^", "2"}]}], ")"}], "^", "2"}]]}], ",",
RowBox[{"{",
RowBox[{"u", ",",
FractionBox[
RowBox[{
RowBox[{"A1", " ",
RowBox[{"Cos", "[",
RowBox[{"3", "y"}], "]"}]}], "-", "1"}],
RowBox[{"A1", "-",
RowBox[{"Cos", "[",
RowBox[{"3", "y"}], "]"}]}]], ",", "1"}], "}"}], ",",
RowBox[{"WorkingPrecision", "\[Rule]", "50"}]}], "]"}]}],
")"}]}]}]], "Input",
CellChangeTimes->{{3.844898517157793*^9, 3.8448985349675703`*^9}, {
3.844898580371933*^9, 3.844898654307754*^9}, {3.844898727190341*^9,
3.8448987289338293`*^9}},
CellLabel->"In[29]:=",ExpressionUUID->"63c342f7-4715-4ee2-ae45-9b525adcea69"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"J2", "[", "1", "]"}]], "Input",
CellChangeTimes->{{3.844898634076033*^9, 3.844898635386674*^9}},
CellLabel->"In[30]:=",ExpressionUUID->"caa8606f-957b-48d0-be91-273ac08b840f"],
Cell[BoxData["5.\
14608510651588948875063444866083409185868746013229905061300573509705902037117`\
50."], "Output",
CellChangeTimes->{{3.844898635877243*^9, 3.844898659158341*^9}, {
3.844898731031642*^9, 3.844898747119618*^9}, 3.844902617703794*^9,
3.8449030394873543`*^9},
CellLabel->"Out[30]=",ExpressionUUID->"e63806f7-d0a3-4153-879c-253ab4ffabbc"]
}, Open ]],
Cell[BoxData[
RowBox[{
RowBox[{"alg1", "=",
RowBox[{
FractionBox["4", "3"],
RowBox[{"(",
RowBox[{
FractionBox[
RowBox[{"2",
RowBox[{"Sqrt", "[",
RowBox[{
RowBox[{"(",
RowBox[{"1", "-",
RowBox[{"u", "^", "2"}]}], ")"}], " ",
RowBox[{"(",
RowBox[{
RowBox[{"4",
RowBox[{"u", "^", "2"}]}], "+",
RowBox[{"16", "u"}], "+", "25"}], ")"}]}], "]"}]}],
RowBox[{"29", "+",
RowBox[{"16", " ", "u"}]}]], "+",
RowBox[{"ArcTan", "[",
SqrtBox[
FractionBox[
RowBox[{"1", "-",
RowBox[{"u", "^", "2"}]}],
RowBox[{
RowBox[{"u", "^", "2"}], "+",
RowBox[{"4", "u"}], "+",
RowBox[{"(",
RowBox[{"25", "/", "4"}], ")"}]}]]], "]"}]}], ")"}]}]}], ";"}]], \
"Input",
CellChangeTimes->{{3.8448989182513123`*^9, 3.844898949804069*^9}, {
3.8448989876275587`*^9, 3.84489905707513*^9}, {3.8448991539329433`*^9,
3.844899230368882*^9}, {3.8448992891892147`*^9, 3.84489929059109*^9}, {
3.844899325557645*^9, 3.844899374544883*^9}, {3.844899443780855*^9,
3.844899460753311*^9}, {3.844899563942567*^9, 3.844899672299459*^9},
3.844899847998584*^9},
CellLabel->"In[31]:=",ExpressionUUID->"693fe376-9cb2-456f-8dd7-7c3f87c04259"],
Cell[BoxData[
RowBox[{
RowBox[{"J3", "[", "y_", "]"}], ":=",
RowBox[{
RowBox[{"(",
RowBox[{"alg1", "/.",
RowBox[{"u", "->",
RowBox[{"Cos", "[",
RowBox[{"3", "y"}], "]"}]}]}], ")"}], "+",
RowBox[{
FractionBox["2", "3"],
RowBox[{"Sqrt", "[", "B1", "]"}],
RowBox[{"NIntegrate", "[",
RowBox[{
RowBox[{
FractionBox["1",
RowBox[{
RowBox[{"Sqrt", "[",
RowBox[{
RowBox[{"u", "^", "2"}], "+", "C1"}], "]"}],
RowBox[{"Sqrt", "[",
RowBox[{"1", "-",
RowBox[{"u", "^", "2"}]}], "]"}]}]],
FractionBox[
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"u", "^", "2"}], "+",
RowBox[{"A1", "^", "2"}]}], ")"}],
RowBox[{"(",
RowBox[{
RowBox[{"u", "^", "2"}], "+", "C1"}], ")"}]}],
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"u", "^", "2"}], "-",
RowBox[{"A1", "^", "2"}]}], ")"}], "^", "2"}]]}], ",",
RowBox[{"{",
RowBox[{"u", ",",
FractionBox[
RowBox[{
RowBox[{"A1", " ",
RowBox[{"Cos", "[",
RowBox[{"3", "y"}], "]"}]}], "-", "1"}],
RowBox[{"A1", "-",
RowBox[{"Cos", "[",
RowBox[{"3", "y"}], "]"}]}]], ",", "1"}], "}"}], ",",
RowBox[{"WorkingPrecision", "\[Rule]", "50"}]}], "]"}]}]}]}]], "Input",
CellChangeTimes->{{3.8448996640422583`*^9, 3.844899725942259*^9}},
CellLabel->"In[32]:=",ExpressionUUID->"c3f35d49-253e-4658-99b5-4f0c7140bde2"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"J3", "[", "1", "]"}]], "Input",
CellChangeTimes->{{3.844899727727159*^9, 3.8448997283643436`*^9}},
CellLabel->"In[33]:=",ExpressionUUID->"e52efa63-5fdb-4ba4-b415-7f3d34ccea18"],
Cell[BoxData["5.\
14608510651588948875063444866083409185868746013230388367653005997510912581235`\
50.01785344904433"], "Output",
CellChangeTimes->{3.84489972877203*^9, 3.844899850280012*^9,
3.844902617892692*^9, 3.844903040066387*^9},
CellLabel->"Out[33]=",ExpressionUUID->"f3453618-0a2a-4785-a590-111a0754221a"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"g1", "=",
RowBox[{"Sqrt", "[",
RowBox[{
RowBox[{"(",
RowBox[{"4", "/", "9"}], ")"}],
RowBox[{"B1", "/",
RowBox[{"(",
RowBox[{"1", "+", "C1"}], ")"}]}]}], "]"}]}]], "Input",
CellChangeTimes->{{3.844866717973673*^9, 3.844866722322034*^9}, {
3.844866752500511*^9, 3.8448667622264023`*^9}, {3.844899936551215*^9,
3.8449000060159063`*^9}},
CellLabel->"In[34]:=",ExpressionUUID->"3ecbb0d8-d5a3-4f13-a243-f7b86c1634f4"],
Cell[BoxData[
SqrtBox[
InterpretationBox[
TemplateBox[{"AlgebraicNumber",
InterpretationBox[
StyleBox[
TemplateBox[{"\"1.18\"",
DynamicBox[
FEPrivate`FrontEndResource[
"FEExpressions", "NumericalApproximationElider"]]}, "RowDefault"],
ShowStringCharacters -> False],
1.17805374429991482898572030535433441401`15.954589770191003, Editable ->
False],
TagBox[
RowBox[{"AlgebraicNumber", "[",
RowBox[{
SqrtBox["65"], ",",
RowBox[{"{",
RowBox[{
RowBox[{"-",
FractionBox["195", "512"]}], ",",
FractionBox["99", "512"]}], "}"}]}], "]"}], Short[#, 7]& ],
1.1780537442999148`},
"NumericalApproximation"],
AlgebraicNumber[65^Rational[1, 2], {
Rational[-195, 512],
Rational[99, 512]}]]]], "Output",
CellChangeTimes->{
3.8448667625572023`*^9, 3.844869650717751*^9, 3.8448719115032167`*^9,
3.8448719752876453`*^9, {3.844899972762166*^9, 3.844900006378354*^9},
3.844902617951717*^9, 3.844903040150221*^9},
CellLabel->"Out[34]=",ExpressionUUID->"d1557a2e-eef7-4691-9405-b32197a6b3af"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"m1", "=",
RowBox[{"1", "/",
RowBox[{"(",
RowBox[{"C1", "+", "1"}], ")"}]}]}]], "Input",
CellChangeTimes->{{3.844865516243475*^9, 3.844865524311955*^9}, {
3.844899991546195*^9, 3.8449000045453377`*^9}},
CellLabel->"In[35]:=",ExpressionUUID->"982bc219-5a2e-4aec-ba9d-a47e18704214"],
Cell[BoxData[
InterpretationBox[
TemplateBox[{"AlgebraicNumber",
InterpretationBox[
StyleBox[
TemplateBox[{"\"0.0659\"",
DynamicBox[
FEPrivate`FrontEndResource[
"FEExpressions", "NumericalApproximationElider"]]}, "RowDefault"],
ShowStringCharacters -> False],
0.06587842893777043351555278150044614449`15.954589770191003, Editable ->
False],
TagBox[
RowBox[{"AlgebraicNumber", "[",
RowBox[{
SqrtBox["65"], ",",
RowBox[{"{",
RowBox[{
FractionBox["1", "2"], ",",
RowBox[{"-",
FractionBox["7", "130"]}]}], "}"}]}], "]"}], Short[#, 7]& ],
0.06587842893777043},
"NumericalApproximation"],
AlgebraicNumber[65^Rational[1, 2], {
Rational[1, 2],
Rational[-7, 130]}]]], "Output",
CellChangeTimes->{
3.844866700482345*^9, 3.844869650743247*^9, 3.844871911550063*^9,
3.8448719753350153`*^9, {3.844899992221983*^9, 3.844900010089295*^9},
3.844902618013587*^9, 3.844903040194872*^9},
CellLabel->"Out[35]=",ExpressionUUID->"dbc3bbe4-05d8-4262-ae1a-91f5dc8522d4"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"n1", "=",
InterpretationBox[
TemplateBox[{"AlgebraicNumber",
InterpretationBox[
StyleBox[
TemplateBox[{"\"-0.0995\"",
DynamicBox[
FEPrivate`FrontEndResource[
"FEExpressions", "NumericalApproximationElider"]]}, "RowDefault"],
ShowStringCharacters ->
False], -0.09950121718117416191873303432657849044`15.954589770191003,
Editable -> False],
TagBox[
RowBox[{"AlgebraicNumber", "[",
RowBox[{
SqrtBox["65"], ",",
RowBox[{"{",
RowBox[{
FractionBox["1", "2"], ",",
RowBox[{"-",
FractionBox["29", "390"]}]}], "}"}]}], "]"}],
Short[#, 7]& ], -0.09950121718117416},
"NumericalApproximation"],
AlgebraicNumber[65^Rational[1, 2], {
Rational[1, 2],
Rational[-29, 390]}]]}]], "Input",
CellChangeTimes->{{3.844900240226255*^9, 3.84490024319298*^9},
3.8449022476835747`*^9},
CellLabel->"In[36]:=",ExpressionUUID->"467d6b72-a554-48d0-ace5-99b150e0b9bf"],
Cell[BoxData[
InterpretationBox[
TemplateBox[{"AlgebraicNumber",
InterpretationBox[
StyleBox[
TemplateBox[{"\"-0.0995\"",
DynamicBox[
FEPrivate`FrontEndResource[
"FEExpressions", "NumericalApproximationElider"]]}, "RowDefault"],
ShowStringCharacters ->
False], -0.09950121718117416191873303432657849044`15.954589770191003,
Editable -> False],
TagBox[
RowBox[{"AlgebraicNumber", "[",
RowBox[{
SqrtBox["65"], ",",
RowBox[{"{",
RowBox[{
FractionBox["1", "2"], ",",
RowBox[{"-",
FractionBox["29", "390"]}]}], "}"}]}], "]"}],
Short[#, 7]& ], -0.09950121718117416},
"NumericalApproximation"],
AlgebraicNumber[65^Rational[1, 2], {
Rational[1, 2],
Rational[-29, 390]}]]], "Output",
CellChangeTimes->{3.8449002436665792`*^9, 3.84490225277836*^9,
3.844902618046238*^9, 3.844903040278564*^9},
CellLabel->"Out[36]=",ExpressionUUID->"a21699ab-75b2-4150-be8a-676a99e26849"]
}, Open ]],
Cell[BoxData[
RowBox[{
RowBox[{"alg2", "=",
RowBox[{
FractionBox["1", "3"],
RowBox[{"(",
RowBox[{
FractionBox[
RowBox[{"Sqrt", "[",
RowBox[{
RowBox[{"(",
RowBox[{"1", "-",
RowBox[{"u", "^", "2"}]}], ")"}], " ",
RowBox[{"(",
RowBox[{
RowBox[{"4",
RowBox[{"u", "^", "2"}]}], "+",
RowBox[{"16", "u"}], "+", "25"}], ")"}]}], "]"}],
RowBox[{"u", "-", "A1"}]], "+",
RowBox[{"4",
RowBox[{"ArcTan", "[",
SqrtBox[
FractionBox[
RowBox[{"1", "-",
RowBox[{"u", "^", "2"}]}],
RowBox[{
RowBox[{"u", "^", "2"}], "+",
RowBox[{"4", "u"}], "+",
RowBox[{"(",
RowBox[{"25", "/", "4"}], ")"}]}]]], "]"}]}]}], ")"}]}]}],
";"}]], "Input",
CellChangeTimes->{
3.844901576257983*^9, {3.844901606779023*^9, 3.8449016432094307`*^9}, {
3.844901750794731*^9, 3.844901842543058*^9}},
CellLabel->"In[37]:=",ExpressionUUID->"1d665d02-7a2e-4866-a7f9-6abbb7cd93a3"],
Cell[BoxData[
RowBox[{
RowBox[{"J4", "[", "y_", "]"}], ":=",
RowBox[{
RowBox[{"With", "[",
RowBox[{
RowBox[{"{",
RowBox[{"\[Phi]", "=",
RowBox[{"ArcCos", "[",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"A1", " ", "u"}], "-", "1"}], ")"}], "/",
RowBox[{"(",
RowBox[{"A1", "-", "u"}], ")"}]}], "]"}]}], "}"}], ",",
RowBox[{"alg2", "+",
RowBox[{"g1",
RowBox[{"(",
RowBox[{
RowBox[{"n1",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"EllipticF", "[",
RowBox[{"\[Phi]", ",", "m1"}], "]"}], "-",
RowBox[{"EllipticE", "[",
RowBox[{"\[Phi]", ",", "m1"}], "]"}]}], ")"}], "/", "m1"}]}],
"+",
RowBox[{
RowBox[{"(",
RowBox[{"16", "/", "9"}], ")"}],
RowBox[{"(",
RowBox[{"1", "-", "m1"}], ")"}], " ",
RowBox[{"EllipticPi", "[",
RowBox[{"n1", ",", "\[Phi]", ",", "m1"}], "]"}]}]}], ")"}]}]}]}],
"]"}], "/.",
RowBox[{"u", "->",
RowBox[{"Cos", "[",
RowBox[{"3", "y"}], "]"}]}]}]}]], "Input",
CellChangeTimes->{
3.844861470998014*^9, {3.844864618023984*^9, 3.844864633534925*^9}, {
3.844864706169724*^9, 3.8448647066461163`*^9}, {3.844864880766942*^9,
3.844864892192451*^9}, {3.844865183412695*^9, 3.8448652054471693`*^9}, {
3.8448652913312798`*^9, 3.8448652963251657`*^9}, {3.844865479676083*^9,
3.84486548915809*^9}, {3.844865543303809*^9, 3.844865591469469*^9}, {
3.844865641500799*^9, 3.844865648990294*^9}, {3.844865749267824*^9,
3.8448657542303963`*^9}, {3.844865784560482*^9, 3.844865786246904*^9},
3.844869407221527*^9, {3.84490000844039*^9, 3.844900055290743*^9}, {
3.844900145628282*^9, 3.844900260163328*^9}, {3.844900291448524*^9,
3.84490051306391*^9}, {3.8449007948034163`*^9, 3.844900804079204*^9}, {
3.844900879916889*^9, 3.8449010019043903`*^9}, {3.844901200679634*^9,
3.844901210323884*^9}, {3.844901307908556*^9, 3.8449013273797493`*^9}, {
3.84490144053258*^9, 3.844901441485568*^9}, {3.844901648935479*^9,
3.844901652005869*^9}, {3.844901769476378*^9, 3.8449017722080393`*^9}, {
3.844901913588644*^9, 3.8449020083667*^9}, {3.844902066989243*^9,
3.844902166148799*^9}, {3.844902334206934*^9, 3.8449023458087187`*^9}, {
3.844902387443713*^9, 3.8449024084242477`*^9}},
CellLabel->"In[38]:=",ExpressionUUID->"3eb697fc-5b91-42d0-8898-79d14a1e5369"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"N", "[",
RowBox[{
RowBox[{"J4", "[", "1", "]"}], ",", "50"}], "]"}]], "Input",
CellChangeTimes->{{3.8449000483240356`*^9, 3.844900057429839*^9}, {
3.844901006971509*^9, 3.844901014685915*^9}},
CellLabel->"In[39]:=",ExpressionUUID->"0de3a3b9-5d7a-4c98-943b-506bf9df2fdd"],
Cell[BoxData["5.\
1460851065158894887506344486608340918586874601323126436041055590439956963499`\
50."], "Output",
CellChangeTimes->{{3.8449000513201427`*^9, 3.844900057752655*^9},
3.844900148978487*^9, {3.84490025248578*^9, 3.844900262082159*^9},
3.844900385737276*^9, {3.844900435840775*^9, 3.844900447511126*^9},
3.844900500474565*^9, {3.844901004730774*^9, 3.844901015014758*^9}, {
3.844901214131445*^9, 3.8449012314038363`*^9}, {3.844901312407398*^9,
3.844901330312584*^9}, 3.844901446717539*^9, 3.844901681562462*^9, {
3.844901774265665*^9, 3.844901793313966*^9}, 3.8449018449438334`*^9,
3.84490188061674*^9, 3.844901939866103*^9, 3.8449021411221333`*^9,
3.844902172376966*^9, 3.84490242075441*^9, 3.844902618146935*^9,
3.8449030404318867`*^9},
CellLabel->"Out[39]=",ExpressionUUID->"da9babd6-d28b-40fc-b7fd-7aa4f7630ec7"]
}, Open ]]
},
WindowSize->{808, 699},
WindowMargins->{{Automatic, 60}, {Automatic, 37}},
FrontEndVersion->"12.3 for Mac OS X x86 (64-bit) (May 11, 2021)",
StyleDefinitions->"Default.nb",
ExpressionUUID->"9e593016-283b-4f1b-a71a-067465be7445"
]
(* End of Notebook Content *)
(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 1210, 26, 30, "Input",ExpressionUUID->"2c2a9424-a764-4819-aa10-799f2cd948c9"],
Cell[CellGroupData[{
Cell[1793, 50, 401, 6, 30, "Input",ExpressionUUID->"d08b93dc-4fd6-4d10-ad2c-5f539f664970"],
Cell[2197, 58, 760, 12, 34, "Output",ExpressionUUID->"18ca9d79-c636-45c2-8a3d-1b1b1ba5e2c3"]
}, Open ]],
Cell[2972, 73, 1130, 28, 50, "Input",ExpressionUUID->"db458be3-0816-4211-97f2-064142b64d1c"],
Cell[CellGroupData[{
Cell[4127, 105, 347, 5, 30, "Input",ExpressionUUID->"2f6514a5-d960-4f90-8ca6-c54f4a5a8b42"],
Cell[4477, 112, 961, 14, 34, "Output",ExpressionUUID->"dc204ccf-5eef-40d8-a14c-950b2f65ef69"]
}, Open ]],
Cell[5453, 129, 730, 22, 49, "Input",ExpressionUUID->"1a851957-c3c7-4afc-9874-71074aa704d1"],
Cell[CellGroupData[{
Cell[6208, 155, 1749, 50, 136, "Input",ExpressionUUID->"7f9843fd-d61b-4867-b7ae-2aae975a95b1"],
Cell[7960, 207, 1273, 37, 68, "Output",ExpressionUUID->"44b572ec-85eb-44b2-91e1-bb822b11f9f0"]
}, Open ]],
Cell[9248, 247, 578, 12, 52, "Input",ExpressionUUID->"487b0363-9a4b-4f76-a23a-7bedb669b210"],
Cell[9829, 261, 1220, 35, 61, "Input",ExpressionUUID->"fd06d2f8-0b02-4d43-b5a6-4892b889e532"],
Cell[11052, 298, 1905, 54, 115, "Input",ExpressionUUID->"17f0ba9b-f3fc-4828-a621-b61d87a8440b"],
Cell[CellGroupData[{
Cell[12982, 356, 552, 12, 30, "Input",ExpressionUUID->"71ff88a2-acd6-4717-b786-efb9232c7818"],
Cell[13537, 370, 3401, 87, 43, "Output",ExpressionUUID->"dbdb8512-329e-4881-b12d-927b6e332805"]
}, Open ]],
Cell[16953, 460, 4092, 93, 177, "Input",ExpressionUUID->"b7f4a108-0135-4c60-aa3f-b15abe8d976d"],
Cell[CellGroupData[{
Cell[21070, 557, 378, 6, 30, "Input",ExpressionUUID->"ef3517f4-ce08-4210-94b7-41051fd5e0e7"],
Cell[21451, 565, 713, 11, 34, "Output",ExpressionUUID->"167e9b9d-c045-45dc-8db4-0fe4f8a1a011"]
}, Open ]],
Cell[CellGroupData[{
Cell[22201, 581, 1288, 37, 60, InheritFromParent,ExpressionUUID->"bdbd3a53-e9f5-4340-a0ab-70c4409eb89a"],
Cell[23492, 620, 1194, 35, 69, "Output",ExpressionUUID->"7423d711-af75-4537-ab65-3a2fe273ea96"]
}, Open ]],
Cell[24701, 658, 1832, 52, 114, "Input",ExpressionUUID->"299d4625-771d-4a6a-884b-2f79fb7634c3"],
Cell[CellGroupData[{
Cell[26558, 714, 449, 7, 30, "Input",ExpressionUUID->"e1ca574f-fa9b-4cc2-9f2f-82520996ed5b"],
Cell[27010, 723, 558, 8, 34, "Output",ExpressionUUID->"df6a278e-ca16-49d1-b9ed-479a03ca9c61"]
}, Open ]],
Cell[CellGroupData[{
Cell[27605, 736, 353, 8, 30, "Input",ExpressionUUID->"bf1012a0-2be9-4b29-815d-104413e2177c"],
Cell[27961, 746, 1157, 30, 45, "Output",ExpressionUUID->"162606d8-3ff7-4f3a-b398-87992b2cfc20"]
}, Open ]],
Cell[CellGroupData[{
Cell[29155, 781, 267, 6, 30, "Input",ExpressionUUID->"005d25f4-570c-4096-842a-ce05c17c86ae"],
Cell[29425, 789, 1132, 30, 41, "Output",ExpressionUUID->"b052984b-2f52-4d1d-a7f4-8f89dd240263"]
}, Open ]],
Cell[30572, 822, 2142, 57, 116, "Input",ExpressionUUID->"3cef2ab3-6910-4a01-bc70-fbf276c56415"],
Cell[CellGroupData[{
Cell[32739, 883, 204, 3, 30, "Input",ExpressionUUID->"2cc6d49e-094d-46db-ab40-e274ab172252"],
Cell[32946, 888, 486, 7, 34, "Output",ExpressionUUID->"266e04d5-b873-4b21-b94f-e2d6f9bac4de"]
}, Open ]],
Cell[CellGroupData[{
Cell[33469, 900, 1045, 28, 37, "Input",ExpressionUUID->"914d2159-72c9-4d2e-91f5-66d003149793"],
Cell[34517, 930, 1118, 29, 41, "Output",ExpressionUUID->"117492f6-aca7-4ce5-8f5e-44c16f57cc41"]
}, Open ]],
Cell[35650, 962, 1270, 38, 60, InheritFromParent,ExpressionUUID->"d6140dfb-24c7-4947-a556-27087f8651b5"],
Cell[36923, 1002, 2944, 62, 94, "Input",ExpressionUUID->"fc7d6d4a-cf1d-4446-889d-2aef239a90ea"],
Cell[CellGroupData[{
Cell[39892, 1068, 308, 6, 30, "Input",ExpressionUUID->"e249b5a5-a6ee-43e9-b638-462b0a2bd512"],
Cell[40203, 1076, 408, 7, 34, "Output",ExpressionUUID->"8f0cedba-c5a1-4159-8b5f-ee8a62e1e86a"]
}, Open ]],
Cell[40626, 1086, 744, 20, 30, "Input",ExpressionUUID->"f12a59db-60f6-4605-b5e1-f6d60953ad18"],
Cell[CellGroupData[{
Cell[41395, 1110, 204, 3, 30, "Input",ExpressionUUID->"783402ff-f7a3-4338-9491-b1687227e07f"],
Cell[41602, 1115, 331, 5, 34, "Output",ExpressionUUID->"68f0a497-0b68-45de-b032-5b56f0bdfda2"]
}, Open ]],
Cell[41948, 1123, 910, 25, 50, "Input",ExpressionUUID->"e662bf28-549b-47ce-b970-b0032324e207"],
Cell[CellGroupData[{
Cell[42883, 1152, 200, 3, 30, "Input",ExpressionUUID->"edfea8cf-26fb-42e6-9b29-1a7f993033b7"],
Cell[43086, 1157, 409, 7, 34, "Output",ExpressionUUID->"7a5e3764-722f-48f0-a6d5-301cc5076147"]
}, Open ]],
Cell[CellGroupData[{
Cell[43532, 1169, 431, 10, 30, "Input",ExpressionUUID->"c59fb042-bbab-44e6-be92-674aed8698c5"],
Cell[43966, 1181, 2862, 79, 42, "Output",ExpressionUUID->"b02477d3-3227-4171-b7bc-b6c19809dbfe"]
}, Open ]],
Cell[46843, 1263, 2488, 74, 177, "Input",ExpressionUUID->"63c342f7-4715-4ee2-ae45-9b525adcea69"],
Cell[CellGroupData[{
Cell[49356, 1341, 200, 3, 30, "Input",ExpressionUUID->"caa8606f-957b-48d0-be91-273ac08b840f"],
Cell[49559, 1346, 359, 6, 34, "Output",ExpressionUUID->"e63806f7-d0a3-4153-879c-253ab4ffabbc"]
}, Open ]],
Cell[49933, 1355, 1372, 38, 60, InheritFromParent,ExpressionUUID->"693fe376-9cb2-456f-8dd7-7c3f87c04259"],
Cell[51308, 1395, 1609, 49, 114, "Input",ExpressionUUID->"c3f35d49-253e-4658-99b5-4f0c7140bde2"],
Cell[CellGroupData[{
Cell[52942, 1448, 202, 3, 30, "Input",ExpressionUUID->"e52efa63-5fdb-4ba4-b415-7f3d34ccea18"],
Cell[53147, 1453, 317, 5, 34, "Output",ExpressionUUID->"f3453618-0a2a-4785-a590-111a0754221a"]
}, Open ]],
Cell[CellGroupData[{
Cell[53501, 1463, 478, 12, 30, "Input",ExpressionUUID->"3ecbb0d8-d5a3-4f13-a243-f7b86c1634f4"],
Cell[53982, 1477, 1182, 31, 45, "Output",ExpressionUUID->"d1557a2e-eef7-4691-9405-b32197a6b3af"]
}, Open ]],
Cell[CellGroupData[{
Cell[55201, 1513, 319, 7, 30, "Input",ExpressionUUID->"982bc219-5a2e-4aec-ba9d-a47e18704214"],
Cell[55523, 1522, 1134, 30, 41, "Output",ExpressionUUID->"dbc3bbe4-05d8-4262-ae1a-91f5dc8522d4"]
}, Open ]],
Cell[CellGroupData[{
Cell[56694, 1557, 1070, 29, 37, "Input",ExpressionUUID->"467d6b72-a554-48d0-ace5-99b150e0b9bf"],
Cell[57767, 1588, 1042, 28, 41, "Output",ExpressionUUID->"a21699ab-75b2-4150-be8a-676a99e26849"]
}, Open ]],
Cell[58824, 1619, 1095, 34, 60, "Input",ExpressionUUID->"1d665d02-7a2e-4866-a7f9-6abbb7cd93a3"],
Cell[59922, 1655, 2530, 57, 115, "Input",ExpressionUUID->"3eb697fc-5b91-42d0-8898-79d14a1e5369"],
Cell[CellGroupData[{
Cell[62477, 1716, 304, 6, 30, "Input",ExpressionUUID->"0de3a3b9-5d7a-4c98-943b-506bf9df2fdd"],
Cell[62784, 1724, 863, 13, 34, "Output",ExpressionUUID->"da9babd6-d28b-40fc-b7fd-7aa4f7630ec7"]
}, Open ]]
}
]
*)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment