1093 lines
40 KiB
Mathematica
1093 lines
40 KiB
Mathematica
(* Content-type: application/vnd.wolfram.mathematica *)
|
|
|
|
(*** Wolfram Notebook File ***)
|
|
(* http://www.wolfram.com/nb *)
|
|
|
|
(* CreatedBy='Mathematica 12.2' *)
|
|
|
|
(*CacheID: 234*)
|
|
(* Internal cache information:
|
|
NotebookFileLineBreakTest
|
|
NotebookFileLineBreakTest
|
|
NotebookDataPosition[ 158, 7]
|
|
NotebookDataLength[ 40718, 1085]
|
|
NotebookOptionsPosition[ 37729, 1034]
|
|
NotebookOutlinePosition[ 38219, 1052]
|
|
CellTagsIndexPosition[ 38176, 1049]
|
|
WindowFrame->Normal*)
|
|
|
|
(* Beginning of Notebook Content *)
|
|
Notebook[{
|
|
Cell[BoxData[
|
|
RowBox[{"(*",
|
|
RowBox[{"<<", "NumericalCalculus`"}], "*)"}]], "Input",
|
|
CellChangeTimes->{{3.833368642280746*^9, 3.833368654634403*^9}, {
|
|
3.83336891101714*^9, 3.8333689126310177`*^9}},
|
|
CellLabel->
|
|
"In[136]:=",ExpressionUUID->"bf2e892a-0008-4529-bacd-4d3619ac745e"],
|
|
|
|
Cell["Sierpinski", "Text",
|
|
CellChangeTimes->{{3.833367411622897*^9, 3.833367415650799*^9}, {
|
|
3.833368638331826*^9,
|
|
3.833368641207601*^9}},ExpressionUUID->"a2a33caa-af3c-4a3a-8054-\
|
|
70d79e05c1fa"],
|
|
|
|
Cell[BoxData[
|
|
RowBox[{"(*",
|
|
RowBox[{
|
|
RowBox[{
|
|
SubscriptBox["A", "1"], "=",
|
|
RowBox[{"(", GridBox[{
|
|
{
|
|
FractionBox["1",
|
|
SqrtBox["2"]],
|
|
FractionBox["1",
|
|
SqrtBox["2"]]},
|
|
{"0",
|
|
SqrtBox["2"]}
|
|
}], ")"}]}], ";", "\[IndentingNewLine]",
|
|
RowBox[{
|
|
SubscriptBox["A", "2"], "=",
|
|
RowBox[{"(", GridBox[{
|
|
{
|
|
FractionBox["1",
|
|
SqrtBox["2"]],
|
|
RowBox[{
|
|
RowBox[{"-",
|
|
FractionBox["1",
|
|
RowBox[{"2",
|
|
SqrtBox["2"]}]]}], "+",
|
|
RowBox[{
|
|
FractionBox[
|
|
SqrtBox["3"],
|
|
RowBox[{"2",
|
|
SqrtBox["2"]}]], "\[ImaginaryI]"}]}]},
|
|
{"0",
|
|
SqrtBox["2"]}
|
|
}], ")"}]}], ";", "\[IndentingNewLine]",
|
|
RowBox[{
|
|
SubscriptBox["A", "3"], "=",
|
|
RowBox[{"(", GridBox[{
|
|
{
|
|
FractionBox["1",
|
|
SqrtBox["2"]],
|
|
RowBox[{
|
|
RowBox[{"-",
|
|
FractionBox["1",
|
|
RowBox[{"2",
|
|
SqrtBox["2"]}]]}], "-",
|
|
RowBox[{
|
|
FractionBox[
|
|
SqrtBox["3"],
|
|
RowBox[{"2",
|
|
SqrtBox["2"]}]], "\[ImaginaryI]"}]}]},
|
|
{"0",
|
|
SqrtBox["2"]}
|
|
}], ")"}]}], ";"}], "*)"}]], "Input",
|
|
CellChangeTimes->{{3.8332357196432323`*^9, 3.833235742060638*^9},
|
|
3.833235802856276*^9, {3.833236070029366*^9, 3.83323607752209*^9}, {
|
|
3.83336400179309*^9, 3.833364026698101*^9}, {3.833367367535881*^9,
|
|
3.8333673947696733`*^9}},
|
|
CellLabel->
|
|
"In[137]:=",ExpressionUUID->"be9cc703-c8ea-4b90-862f-061a2b1498c5"],
|
|
|
|
Cell["Apollonian (McMullen\[CloseCurlyQuote]s to unit)", "Text",
|
|
CellChangeTimes->{{3.83336742412843*^9, 3.833367465118862*^9}, {
|
|
3.833373828770335*^9,
|
|
3.833373830363236*^9}},ExpressionUUID->"48ecbeab-626b-4ed1-bffb-\
|
|
20691c3748f2"],
|
|
|
|
Cell[BoxData[
|
|
RowBox[{"(*",
|
|
RowBox[{
|
|
RowBox[{"M", "=",
|
|
RowBox[{"(", GridBox[{
|
|
{
|
|
RowBox[{
|
|
RowBox[{"-", "1"}], "-",
|
|
RowBox[{"2", "\[ImaginaryI]"}]}],
|
|
RowBox[{"1", "+", "\[ImaginaryI]"}]},
|
|
{
|
|
RowBox[{"4", "+",
|
|
RowBox[{"2", "\[ImaginaryI]"}]}],
|
|
RowBox[{"2", "+",
|
|
RowBox[{"2", "\[ImaginaryI]"}]}]}
|
|
}], ")"}]}], ";", "\[IndentingNewLine]",
|
|
RowBox[{
|
|
SubscriptBox["A", "1"], "=",
|
|
RowBox[{"M", ".",
|
|
RowBox[{"(", GridBox[{
|
|
{"1", "1"},
|
|
{"0", "1"}
|
|
}], ")"}], ".",
|
|
RowBox[{"Inverse", "[", "M", "]"}]}]}], ";", "\[IndentingNewLine]",
|
|
RowBox[{
|
|
SubscriptBox["A", "2"], "=",
|
|
RowBox[{"M", ".",
|
|
RowBox[{"(", GridBox[{
|
|
{"1", "0"},
|
|
{"1", "1"}
|
|
}], ")"}], ".",
|
|
RowBox[{"Inverse", "[", "M", "]"}]}]}], ";", "\[IndentingNewLine]",
|
|
RowBox[{
|
|
SubscriptBox["A", "3"], "=",
|
|
RowBox[{"M", ".",
|
|
RowBox[{"(", GridBox[{
|
|
{
|
|
RowBox[{"1", "+", "\[ImaginaryI]"}], "1"},
|
|
{"1",
|
|
RowBox[{"1", "-", "\[ImaginaryI]"}]}
|
|
}], ")"}], ".",
|
|
RowBox[{"Inverse", "[", "M", "]"}]}]}], ";"}], "*)"}]], "Input",
|
|
CellChangeTimes->{{3.8333674000285597`*^9, 3.8333674021150637`*^9}, {
|
|
3.8333675407218227`*^9, 3.833367599438445*^9}, {3.833373841265559*^9,
|
|
3.833373851104115*^9}, {3.833374173441412*^9, 3.833374175210209*^9}, {
|
|
3.8333767815304327`*^9, 3.8333767821355343`*^9}},
|
|
CellLabel->
|
|
"In[138]:=",ExpressionUUID->"86a5d30e-95b8-493a-a4bf-7613b620acfa"],
|
|
|
|
Cell[BoxData[
|
|
RowBox[{"(*",
|
|
RowBox[{
|
|
RowBox[{
|
|
SubscriptBox["f", "A_"], "[", "z_", "]"}], ":=",
|
|
RowBox[{
|
|
RowBox[{
|
|
FractionBox[
|
|
RowBox[{
|
|
RowBox[{
|
|
RowBox[{"A", "\[LeftDoubleBracket]",
|
|
RowBox[{"1", ",", "1"}], "\[RightDoubleBracket]"}], "z"}], "+",
|
|
RowBox[{"A", "\[LeftDoubleBracket]",
|
|
RowBox[{"1", ",", "2"}], "\[RightDoubleBracket]"}]}],
|
|
RowBox[{
|
|
RowBox[{
|
|
RowBox[{"A", "\[LeftDoubleBracket]",
|
|
RowBox[{"2", ",", "1"}], "\[RightDoubleBracket]"}], "z"}], "+",
|
|
RowBox[{"A", "\[LeftDoubleBracket]",
|
|
RowBox[{"2", ",", "2"}], "\[RightDoubleBracket]"}]}]],
|
|
"\[IndentingNewLine]",
|
|
RowBox[{
|
|
SubscriptBox[
|
|
SuperscriptBox["f",
|
|
RowBox[{"-", "1"}]], "A_"], "[", "z_", "]"}]}], ":=",
|
|
RowBox[{
|
|
RowBox[{
|
|
FractionBox[
|
|
RowBox[{
|
|
RowBox[{
|
|
RowBox[{
|
|
RowBox[{"Inverse", "[", "A", "]"}], "\[LeftDoubleBracket]",
|
|
RowBox[{"1", ",", "1"}], "\[RightDoubleBracket]"}], "z"}], "+",
|
|
RowBox[{
|
|
RowBox[{"Inverse", "[", "A", "]"}], "\[LeftDoubleBracket]",
|
|
RowBox[{"1", ",", "2"}], "\[RightDoubleBracket]"}]}],
|
|
RowBox[{
|
|
RowBox[{
|
|
RowBox[{
|
|
RowBox[{"Inverse", "[", "A", "]"}], "\[LeftDoubleBracket]",
|
|
RowBox[{"2", ",", "1"}], "\[RightDoubleBracket]"}], "z"}], "+",
|
|
RowBox[{
|
|
RowBox[{"Inverse", "[", "A", "]"}], "\[LeftDoubleBracket]",
|
|
RowBox[{"2", ",", "2"}], "\[RightDoubleBracket]"}]}]],
|
|
"\[IndentingNewLine]",
|
|
RowBox[{
|
|
SubscriptBox["M",
|
|
RowBox[{"A_", ",", "q_"}]], "[",
|
|
RowBox[{"n_", ",", "s_"}], "]"}]}], ":=",
|
|
RowBox[{
|
|
RowBox[{
|
|
RowBox[{
|
|
RowBox[{"D", "[",
|
|
RowBox[{
|
|
FractionBox[
|
|
SuperscriptBox[
|
|
RowBox[{"(",
|
|
RowBox[{
|
|
RowBox[{
|
|
RowBox[{"A", "\[LeftDoubleBracket]",
|
|
RowBox[{"1", ",", "1"}], "\[RightDoubleBracket]"}], "z"}],
|
|
"+",
|
|
RowBox[{"A", "\[LeftDoubleBracket]",
|
|
RowBox[{"1", ",", "2"}], "\[RightDoubleBracket]"}]}], ")"}],
|
|
"n"],
|
|
SuperscriptBox[
|
|
RowBox[{"(",
|
|
RowBox[{
|
|
RowBox[{
|
|
RowBox[{"A", "\[LeftDoubleBracket]",
|
|
RowBox[{"2", ",", "1"}], "\[RightDoubleBracket]"}], "z"}],
|
|
"+",
|
|
RowBox[{"A", "\[LeftDoubleBracket]",
|
|
RowBox[{"2", ",", "2"}], "\[RightDoubleBracket]"}]}], ")"}],
|
|
RowBox[{"n", "+", "q"}]]], ",",
|
|
RowBox[{"{",
|
|
RowBox[{"z", ",", "s"}], "}"}]}], "]"}],
|
|
RowBox[{"(",
|
|
FractionBox["1",
|
|
RowBox[{"s", "!"}]], ")"}]}], "/.",
|
|
RowBox[{
|
|
RowBox[{"{",
|
|
RowBox[{"z", "\[Rule]", "0"}], "}"}], "\[IndentingNewLine]",
|
|
RowBox[{
|
|
SubscriptBox["L",
|
|
RowBox[{"A_", ",", "q_"}]], "[",
|
|
RowBox[{"m_", ",", "n_", ",", "r_", ",", "s_"}], "]"}]}]}], ":=",
|
|
RowBox[{
|
|
RowBox[{
|
|
RowBox[{
|
|
SubscriptBox["M",
|
|
RowBox[{"A", ",", "q"}]], "[",
|
|
RowBox[{"m", ",", "r"}], "]"}], "\[Conjugate]"}],
|
|
RowBox[{
|
|
SubscriptBox["M",
|
|
RowBox[{"A", ",", "q"}]], "[",
|
|
RowBox[{"n", ",", "s"}], "]"}]}]}]}]}]}], "*)"}]], "Input",
|
|
CellChangeTimes->{{3.833376783564871*^9, 3.833376785303887*^9},
|
|
3.83344967241963*^9},
|
|
CellLabel->
|
|
"In[139]:=",ExpressionUUID->"cae51aa1-774a-45bc-8ddb-2d9d2bac6186"],
|
|
|
|
Cell["Apollonian (Bai-Finch)", "Text",
|
|
CellChangeTimes->{{3.833373910951479*^9, 3.833373938990329*^9}, {
|
|
3.833447554710703*^9,
|
|
3.833447570194236*^9}},ExpressionUUID->"aa6c1252-995d-4016-89e4-\
|
|
9cfcf4fb8ba5"],
|
|
|
|
Cell[BoxData[{
|
|
RowBox[{
|
|
RowBox[{"M", "=",
|
|
RowBox[{"(", GridBox[{
|
|
{
|
|
RowBox[{
|
|
RowBox[{"-", "1"}], "-",
|
|
RowBox[{"2", "\[ImaginaryI]"}]}],
|
|
RowBox[{"1", "+", "\[ImaginaryI]"}]},
|
|
{
|
|
RowBox[{"4", "+",
|
|
RowBox[{"2", "\[ImaginaryI]"}]}],
|
|
RowBox[{"2", "+",
|
|
RowBox[{"2", "\[ImaginaryI]"}]}]}
|
|
}], ")"}]}], ";"}], "\[IndentingNewLine]",
|
|
RowBox[{
|
|
RowBox[{"R", "=",
|
|
RowBox[{
|
|
FractionBox["1", "12"],
|
|
RowBox[{"(", GridBox[{
|
|
{
|
|
RowBox[{
|
|
RowBox[{"-", "6"}], "+",
|
|
RowBox[{"8", "\[ImaginaryI]"}]}],
|
|
RowBox[{"11", "\[ImaginaryI]"}]},
|
|
{
|
|
RowBox[{"4", "\[ImaginaryI]"}],
|
|
RowBox[{
|
|
RowBox[{"-", "6"}], "-",
|
|
RowBox[{"8", "\[ImaginaryI]"}]}]}
|
|
}], ")"}]}]}], ";"}], "\[IndentingNewLine]",
|
|
RowBox[{
|
|
RowBox[{"G", "=",
|
|
RowBox[{
|
|
FractionBox["1", "12"],
|
|
RowBox[{"(", GridBox[{
|
|
{
|
|
RowBox[{
|
|
RowBox[{"-", "2"}], "-",
|
|
RowBox[{"2", "\[ImaginaryI]"}]}],
|
|
RowBox[{
|
|
RowBox[{"-", "1"}], "-",
|
|
RowBox[{"5", "\[ImaginaryI]"}]}]},
|
|
{
|
|
RowBox[{"4", "-",
|
|
RowBox[{"4", "\[ImaginaryI]"}]}],
|
|
RowBox[{
|
|
RowBox[{"-", "2"}], "-",
|
|
RowBox[{"10", "\[ImaginaryI]"}]}]}
|
|
}], ")"}]}]}], ";"}], "\[IndentingNewLine]",
|
|
RowBox[{
|
|
RowBox[{
|
|
SubscriptBox["A", "k_"], ":=",
|
|
RowBox[{"R", "+",
|
|
RowBox[{"k", " ", "G"}]}]}], ";"}], "\[IndentingNewLine]",
|
|
RowBox[{
|
|
RowBox[{
|
|
SubscriptBox["k", "0"], "=", "100"}], ";"}], "\[IndentingNewLine]",
|
|
RowBox[{
|
|
RowBox[{
|
|
SubscriptBox["N", "c"], "=", "10"}], ";"}], "\[IndentingNewLine]",
|
|
RowBox[{
|
|
RowBox[{"Lc", "=", "6"}], ";"}]}], "Input",
|
|
CellChangeTimes->{{3.833373852642791*^9, 3.833373856486021*^9},
|
|
3.833373912862501*^9, {3.833373953244506*^9, 3.833374145561005*^9},
|
|
3.8333742186253242`*^9, {3.833374267333*^9, 3.833374267492844*^9},
|
|
3.833374327175881*^9, {3.833376466438765*^9, 3.833376502268091*^9}, {
|
|
3.8333772003399982`*^9, 3.833377201970336*^9}, {3.833378691497251*^9,
|
|
3.833378691642603*^9}, 3.833378785570231*^9, {3.8333788481941338`*^9,
|
|
3.833378848267768*^9}, {3.8334474384956913`*^9, 3.8334474386980953`*^9}, {
|
|
3.833448076160449*^9, 3.833448076672225*^9}, {3.8334481145341053`*^9,
|
|
3.833448125168035*^9}, {3.833448248753934*^9, 3.833448250632118*^9},
|
|
3.833450009474032*^9, {3.8334500717143087`*^9, 3.833450084607336*^9}, {
|
|
3.8334501284393167`*^9, 3.8334501323161707`*^9}, {3.833450163301908*^9,
|
|
3.83345016444783*^9}, {3.833451466703332*^9, 3.833451505172805*^9}, {
|
|
3.83345187734093*^9, 3.833451885320603*^9}, {3.833457215301762*^9,
|
|
3.833457215768548*^9}, {3.8334574826300993`*^9, 3.833457485156622*^9}, {
|
|
3.8334593722986507`*^9, 3.833459376867346*^9}, {3.833459729519393*^9,
|
|
3.833459729578371*^9}},
|
|
CellLabel->
|
|
"In[140]:=",ExpressionUUID->"746ed5e5-49cf-4db7-8474-08917e9d577c"],
|
|
|
|
Cell[BoxData[
|
|
RowBox[{
|
|
RowBox[{"(*",
|
|
RowBox[{
|
|
RowBox[{
|
|
SubscriptBox["F", "q_"], "[",
|
|
RowBox[{"n_", ",", "s_", ",", "l_"}], "]"}], ":=",
|
|
RowBox[{
|
|
RowBox[{"(",
|
|
RowBox[{
|
|
RowBox[{"(",
|
|
FractionBox["1",
|
|
RowBox[{
|
|
RowBox[{"l", "!"}],
|
|
RowBox[{"s", "!"}]}]], ")"}],
|
|
RowBox[{"D", "[",
|
|
RowBox[{
|
|
RowBox[{"(",
|
|
RowBox[{
|
|
RowBox[{"D", "[",
|
|
RowBox[{
|
|
FractionBox[
|
|
SuperscriptBox[
|
|
RowBox[{"(",
|
|
RowBox[{
|
|
RowBox[{
|
|
RowBox[{"(", " ",
|
|
RowBox[{
|
|
RowBox[{"G", "\[LeftDoubleBracket]",
|
|
RowBox[{"1", ",", "1"}], "\[RightDoubleBracket]"}], "+",
|
|
RowBox[{"t", " ",
|
|
RowBox[{"R", "\[LeftDoubleBracket]",
|
|
RowBox[{"1", ",", "1"}], "\[RightDoubleBracket]"}]}]}],
|
|
")"}], "z"}], "+", " ",
|
|
RowBox[{"G", "\[LeftDoubleBracket]",
|
|
RowBox[{"1", ",", "2"}], "\[RightDoubleBracket]"}], "+",
|
|
RowBox[{"t", " ",
|
|
RowBox[{"R", "\[LeftDoubleBracket]",
|
|
RowBox[{"1", ",", "2"}], "\[RightDoubleBracket]"}]}]}],
|
|
")"}], "n"],
|
|
SuperscriptBox[
|
|
RowBox[{"(",
|
|
RowBox[{
|
|
RowBox[{
|
|
RowBox[{"(",
|
|
RowBox[{
|
|
RowBox[{"G", "\[LeftDoubleBracket]",
|
|
RowBox[{"2", ",", "1"}], "\[RightDoubleBracket]"}], "+",
|
|
RowBox[{"t", " ",
|
|
RowBox[{"R", "\[LeftDoubleBracket]",
|
|
RowBox[{"2", ",", "1"}], "\[RightDoubleBracket]"}]}]}],
|
|
")"}], "z"}], "+", " ",
|
|
RowBox[{"G", "\[LeftDoubleBracket]",
|
|
RowBox[{"2", ",", "2"}], "\[RightDoubleBracket]"}], "+",
|
|
RowBox[{"t", " ",
|
|
RowBox[{"R", "\[LeftDoubleBracket]",
|
|
RowBox[{"2", ",", "2"}], "\[RightDoubleBracket]"}]}]}],
|
|
")"}],
|
|
RowBox[{"n", "+", "q"}]]], ",",
|
|
RowBox[{"{",
|
|
RowBox[{"z", ",", "s"}], "}"}]}], "]"}], "/.",
|
|
RowBox[{"{",
|
|
RowBox[{"z", "\[Rule]", "0"}], "}"}]}], ")"}], ",",
|
|
RowBox[{"{",
|
|
RowBox[{"t", ",", "l"}], "}"}]}], "]"}]}], ")"}], "/.",
|
|
RowBox[{"{",
|
|
RowBox[{"t", "\[Rule]", "0"}], "}"}]}]}], "*)"}], "\[IndentingNewLine]",
|
|
RowBox[{
|
|
RowBox[{
|
|
RowBox[{
|
|
SubscriptBox["F", "q_"], "[",
|
|
RowBox[{"n_", ",", "s_", ",", "l_"}], "]"}], "=",
|
|
RowBox[{
|
|
RowBox[{"SeriesCoefficient", "[",
|
|
RowBox[{
|
|
RowBox[{"Series", "[",
|
|
RowBox[{
|
|
FractionBox[
|
|
SuperscriptBox[
|
|
RowBox[{"(",
|
|
RowBox[{
|
|
RowBox[{
|
|
RowBox[{"(", " ",
|
|
RowBox[{
|
|
RowBox[{"G", "\[LeftDoubleBracket]",
|
|
RowBox[{"1", ",", "1"}], "\[RightDoubleBracket]"}], "+",
|
|
RowBox[{"t", " ",
|
|
RowBox[{"R", "\[LeftDoubleBracket]",
|
|
RowBox[{"1", ",", "1"}], "\[RightDoubleBracket]"}]}]}],
|
|
")"}], "z"}], "+", " ",
|
|
RowBox[{"G", "\[LeftDoubleBracket]",
|
|
RowBox[{"1", ",", "2"}], "\[RightDoubleBracket]"}], "+",
|
|
RowBox[{"t", " ",
|
|
RowBox[{"R", "\[LeftDoubleBracket]",
|
|
RowBox[{"1", ",", "2"}], "\[RightDoubleBracket]"}]}]}], ")"}],
|
|
"n"],
|
|
SuperscriptBox[
|
|
RowBox[{"(",
|
|
RowBox[{
|
|
RowBox[{
|
|
RowBox[{"(",
|
|
RowBox[{
|
|
RowBox[{"G", "\[LeftDoubleBracket]",
|
|
RowBox[{"2", ",", "1"}], "\[RightDoubleBracket]"}], "+",
|
|
RowBox[{"t", " ",
|
|
RowBox[{"R", "\[LeftDoubleBracket]",
|
|
RowBox[{"2", ",", "1"}], "\[RightDoubleBracket]"}]}]}],
|
|
")"}], "z"}], "+", " ",
|
|
RowBox[{"G", "\[LeftDoubleBracket]",
|
|
RowBox[{"2", ",", "2"}], "\[RightDoubleBracket]"}], "+",
|
|
RowBox[{"t", " ",
|
|
RowBox[{"R", "\[LeftDoubleBracket]",
|
|
RowBox[{"2", ",", "2"}], "\[RightDoubleBracket]"}]}]}], ")"}],
|
|
|
|
RowBox[{"n", "+", "q"}]]], ",",
|
|
RowBox[{"{",
|
|
RowBox[{"z", ",", "0", ",",
|
|
SubscriptBox["N", "c"]}], "}"}], ",",
|
|
RowBox[{"{",
|
|
RowBox[{"t", ",", "0", ",", "Lc"}], "}"}]}], "]"}], ",",
|
|
RowBox[{"{",
|
|
RowBox[{"s", ",", "l"}], "}"}]}], "]"}], ":=",
|
|
RowBox[{"SeriesCoefficient", "[",
|
|
RowBox[{
|
|
RowBox[{"Series", "[",
|
|
RowBox[{
|
|
FractionBox[
|
|
SuperscriptBox[
|
|
RowBox[{"(",
|
|
RowBox[{
|
|
RowBox[{
|
|
RowBox[{"(", " ",
|
|
RowBox[{
|
|
RowBox[{"G", "\[LeftDoubleBracket]",
|
|
RowBox[{"1", ",", "1"}], "\[RightDoubleBracket]"}], "+",
|
|
RowBox[{"t", " ",
|
|
RowBox[{"R", "\[LeftDoubleBracket]",
|
|
RowBox[{"1", ",", "1"}], "\[RightDoubleBracket]"}]}]}],
|
|
")"}], "z"}], "+", " ",
|
|
RowBox[{"G", "\[LeftDoubleBracket]",
|
|
RowBox[{"1", ",", "2"}], "\[RightDoubleBracket]"}], "+",
|
|
RowBox[{"t", " ",
|
|
RowBox[{"R", "\[LeftDoubleBracket]",
|
|
RowBox[{"1", ",", "2"}], "\[RightDoubleBracket]"}]}]}], ")"}],
|
|
"n"],
|
|
SuperscriptBox[
|
|
RowBox[{"(",
|
|
RowBox[{
|
|
RowBox[{
|
|
RowBox[{"(",
|
|
RowBox[{
|
|
RowBox[{"G", "\[LeftDoubleBracket]",
|
|
RowBox[{"2", ",", "1"}], "\[RightDoubleBracket]"}], "+",
|
|
RowBox[{"t", " ",
|
|
RowBox[{"R", "\[LeftDoubleBracket]",
|
|
RowBox[{"2", ",", "1"}], "\[RightDoubleBracket]"}]}]}],
|
|
")"}], "z"}], "+", " ",
|
|
RowBox[{"G", "\[LeftDoubleBracket]",
|
|
RowBox[{"2", ",", "2"}], "\[RightDoubleBracket]"}], "+",
|
|
RowBox[{"t", " ",
|
|
RowBox[{"R", "\[LeftDoubleBracket]",
|
|
RowBox[{"2", ",", "2"}], "\[RightDoubleBracket]"}]}]}], ")"}],
|
|
|
|
RowBox[{"n", "+", "q"}]]], ",",
|
|
RowBox[{"{",
|
|
RowBox[{"z", ",", "0", ",",
|
|
SubscriptBox["N", "c"]}], "}"}], ",",
|
|
RowBox[{"{",
|
|
RowBox[{"t", ",", "0", ",", "Lc"}], "}"}]}], "]"}], ",",
|
|
RowBox[{"{",
|
|
RowBox[{"s", ",", "l"}], "}"}]}], "]"}]}]}], "\[IndentingNewLine]",
|
|
RowBox[{
|
|
RowBox[{
|
|
SubscriptBox["L",
|
|
RowBox[{"k_", ",", "q_"}]], "[",
|
|
RowBox[{"m_", ",", "n_", ",", "r_", ",", "s_"}], "]"}], ":=",
|
|
RowBox[{
|
|
UnderoverscriptBox["\[Sum]",
|
|
RowBox[{"l", "=", "0"}], "Lc"],
|
|
RowBox[{"(",
|
|
RowBox[{
|
|
SuperscriptBox["k",
|
|
RowBox[{"-",
|
|
RowBox[{"(",
|
|
RowBox[{"l", "+",
|
|
RowBox[{"2", "q"}]}], ")"}]}]],
|
|
RowBox[{
|
|
UnderoverscriptBox["\[Sum]",
|
|
RowBox[{"lp", "=", "0"}], "l"], " ",
|
|
RowBox[{
|
|
RowBox[{
|
|
RowBox[{
|
|
SubscriptBox["F", "q"], "[",
|
|
RowBox[{"m", ",", "r", ",",
|
|
RowBox[{"l", "-", "lp"}]}], "]"}], "\[Conjugate]"}],
|
|
RowBox[{
|
|
SubscriptBox["F", "q"], "[",
|
|
RowBox[{"n", ",", "s", ",", "lp"}], "]"}]}]}]}],
|
|
")"}]}]}]}]}]], "Input",
|
|
CellChangeTimes->{{3.833375448945825*^9, 3.8333755889655933`*^9}, {
|
|
3.833375808794491*^9, 3.8333758831917677`*^9}, {3.8333759170572987`*^9,
|
|
3.8333761729999723`*^9}, {3.833376203100834*^9, 3.833376208114214*^9}, {
|
|
3.83337624142187*^9, 3.833376371027378*^9}, {3.833376418142404*^9,
|
|
3.8333764184979477`*^9}, {3.833377668413629*^9, 3.8333776897111683`*^9}, {
|
|
3.8333778053908663`*^9, 3.833377807319936*^9}, {3.833377848815303*^9,
|
|
3.833377849244471*^9}, {3.833378311144019*^9, 3.833378399389579*^9}, {
|
|
3.833378431474545*^9, 3.833378445400144*^9}, {3.833378507575803*^9,
|
|
3.8333785681349707`*^9}, {3.833378601004755*^9, 3.833378617100129*^9}, {
|
|
3.833378910585396*^9, 3.833378922369632*^9}, {3.8333790087934523`*^9,
|
|
3.833379010710129*^9}, {3.833379889197473*^9, 3.833379941043254*^9}, {
|
|
3.8333806621701603`*^9, 3.833380699161913*^9}, {3.833408479553835*^9,
|
|
3.8334084818164463`*^9}, {3.833408554618367*^9, 3.833408577788725*^9}, {
|
|
3.833408777006459*^9, 3.833408799828104*^9}, 3.833408835310032*^9, {
|
|
3.833408907930086*^9, 3.833408936850506*^9}, 3.833408977217951*^9, {
|
|
3.833409009042762*^9, 3.833409011452094*^9}, {3.8334090414606667`*^9,
|
|
3.8334090689933653`*^9}, {3.833409116263338*^9, 3.833409146904228*^9}, {
|
|
3.833410275342947*^9, 3.833410338150243*^9}, {3.833410502432825*^9,
|
|
3.833410536539865*^9}, {3.833410677414356*^9, 3.833410680111587*^9},
|
|
3.833410734400564*^9, {3.833410802999261*^9, 3.833410808571859*^9}, {
|
|
3.833410978148294*^9, 3.833411010767108*^9}, {3.8334469266966467`*^9,
|
|
3.833446944879003*^9}, {3.833446992683023*^9, 3.833446993858767*^9}, {
|
|
3.833447065678322*^9, 3.8334470659001217`*^9}, {3.833447678159583*^9,
|
|
3.833447721262374*^9}, {3.8334477615011787`*^9, 3.83344779379137*^9}, {
|
|
3.833447841087942*^9, 3.833447942430419*^9}, {3.8334480115433617`*^9,
|
|
3.833448015497703*^9}, {3.833449621867939*^9, 3.83344962404976*^9}, {
|
|
3.833449804774026*^9, 3.833449846882372*^9}, {3.8334511543683767`*^9,
|
|
3.833451158127261*^9}, {3.833451325016913*^9, 3.833451326898633*^9}, {
|
|
3.833461548653658*^9,
|
|
3.8334615493688726`*^9}},ExpressionUUID->"784ec4e3-0e47-4303-898e-\
|
|
2c424985901f"],
|
|
|
|
Cell[BoxData[{
|
|
RowBox[{
|
|
RowBox[{
|
|
SubscriptBox["M",
|
|
RowBox[{"A_", ",", "q_"}]], "[",
|
|
RowBox[{"n_", ",", "s_"}], "]"}], ":=",
|
|
RowBox[{"SeriesCoefficient", "[",
|
|
RowBox[{
|
|
RowBox[{"Series", "[",
|
|
RowBox[{
|
|
FractionBox[
|
|
SuperscriptBox[
|
|
RowBox[{"(",
|
|
RowBox[{
|
|
RowBox[{
|
|
RowBox[{"A", "\[LeftDoubleBracket]",
|
|
RowBox[{"1", ",", "1"}], "\[RightDoubleBracket]"}], "z"}], "+",
|
|
RowBox[{"A", "\[LeftDoubleBracket]",
|
|
RowBox[{"1", ",", "2"}], "\[RightDoubleBracket]"}]}], ")"}], "n"],
|
|
|
|
SuperscriptBox[
|
|
RowBox[{"(",
|
|
RowBox[{
|
|
RowBox[{
|
|
RowBox[{"A", "\[LeftDoubleBracket]",
|
|
RowBox[{"2", ",", "1"}], "\[RightDoubleBracket]"}], "z"}], "+",
|
|
RowBox[{"A", "\[LeftDoubleBracket]",
|
|
RowBox[{"2", ",", "2"}], "\[RightDoubleBracket]"}]}], ")"}],
|
|
RowBox[{"n", "+", "q"}]]], ",",
|
|
RowBox[{"{",
|
|
RowBox[{"z", ",", "0", ",", "s"}], "}"}]}], "]"}], ",", "s"}],
|
|
"]"}]}], "\[IndentingNewLine]",
|
|
RowBox[{
|
|
RowBox[{
|
|
SubscriptBox["L",
|
|
RowBox[{"A_", ",", "q_"}]], "[",
|
|
RowBox[{"m_", ",", "n_", ",", "r_", ",", "s_"}], "]"}], ":=",
|
|
RowBox[{
|
|
RowBox[{
|
|
RowBox[{
|
|
SubscriptBox["M",
|
|
RowBox[{"A", ",", "q"}]], "[",
|
|
RowBox[{"m", ",", "r"}], "]"}], "\[Conjugate]"}],
|
|
RowBox[{
|
|
SubscriptBox["M",
|
|
RowBox[{"A", ",", "q"}]], "[",
|
|
RowBox[{"n", ",", "s"}], "]"}]}]}]}], "Input",
|
|
CellChangeTimes->{{3.833449793340254*^9, 3.833449813260859*^9}, {
|
|
3.8334500897364483`*^9, 3.833450115321815*^9}},
|
|
CellLabel->
|
|
"In[149]:=",ExpressionUUID->"f8e74c1e-8b8a-4bd8-9de6-2e4d3a09a42c"],
|
|
|
|
Cell[BoxData[
|
|
RowBox[{
|
|
RowBox[{
|
|
SubscriptBox["F",
|
|
RowBox[{"l_", ",", "q_"}]], "[",
|
|
RowBox[{"m_", ",", "n_", ",", "r_", ",", "s_"}], "]"}], ":=",
|
|
RowBox[{
|
|
UnderoverscriptBox["\[Sum]",
|
|
RowBox[{"lp", "=", "0"}], "l"], " ",
|
|
RowBox[{
|
|
RowBox[{
|
|
RowBox[{
|
|
SubscriptBox["F", "q"], "[",
|
|
RowBox[{"m", ",", "r", ",",
|
|
RowBox[{"l", "-", "lp"}]}], "]"}], "\[Conjugate]"}],
|
|
RowBox[{
|
|
SubscriptBox["F", "q"], "[",
|
|
RowBox[{"n", ",", "s", ",", "lp"}], "]"}]}]}]}]], "Input",
|
|
CellChangeTimes->{{3.83345124461913*^9, 3.8334512487293158`*^9}, {
|
|
3.833451297537081*^9, 3.8334513125181313`*^9}, 3.83345140471535*^9},
|
|
CellLabel->
|
|
"In[151]:=",ExpressionUUID->"90715246-bc56-45e7-8672-f47995804b03"],
|
|
|
|
Cell[BoxData[{
|
|
RowBox[{
|
|
SubscriptBox["L", "q_"], ":=",
|
|
RowBox[{"2",
|
|
RowBox[{"Re", "[",
|
|
RowBox[{
|
|
RowBox[{"(",
|
|
RowBox[{
|
|
UnderoverscriptBox["\[Sum]",
|
|
RowBox[{"k", "=", "1"}],
|
|
RowBox[{
|
|
SubscriptBox["k", "0"], "-", "1"}]],
|
|
RowBox[{"Table", "[",
|
|
RowBox[{
|
|
RowBox[{
|
|
SubscriptBox["L",
|
|
RowBox[{
|
|
SubscriptBox["A", "k"], ",", "q"}]], "[",
|
|
RowBox[{"m", ",", "n", ",", "r", ",", "s"}], "]"}], ",",
|
|
RowBox[{"{",
|
|
RowBox[{"m", ",", "0", ",",
|
|
SubscriptBox["N", "c"]}], "}"}], ",",
|
|
RowBox[{"{",
|
|
RowBox[{"n", ",", "0", ",",
|
|
SubscriptBox["N", "c"]}], "}"}], ",",
|
|
RowBox[{"{",
|
|
RowBox[{"r", ",", "0", ",",
|
|
SubscriptBox["N", "c"]}], "}"}], ",",
|
|
RowBox[{"{",
|
|
RowBox[{"s", ",", "0", ",",
|
|
SubscriptBox["N", "c"]}], "}"}]}], "]"}]}], ")"}], "+",
|
|
RowBox[{"(",
|
|
RowBox[{
|
|
UnderoverscriptBox["\[Sum]",
|
|
RowBox[{"l", "=", "0"}], "Lc"],
|
|
RowBox[{
|
|
RowBox[{"Zeta", "[",
|
|
RowBox[{
|
|
RowBox[{"l", "+",
|
|
RowBox[{"2", "q"}]}], ",",
|
|
SubscriptBox["k", "0"]}], "]"}], " ",
|
|
RowBox[{"Table", "[",
|
|
RowBox[{
|
|
RowBox[{
|
|
SubscriptBox["F",
|
|
RowBox[{"l", ",", "q"}]], "[",
|
|
RowBox[{"m", ",", "n", ",", "r", ",", "s"}], "]"}], ",",
|
|
RowBox[{"{",
|
|
RowBox[{"m", ",", "0", ",",
|
|
SubscriptBox["N", "c"]}], "}"}], ",",
|
|
RowBox[{"{",
|
|
RowBox[{"n", ",", "0", ",",
|
|
SubscriptBox["N", "c"]}], "}"}], ",",
|
|
RowBox[{"{",
|
|
RowBox[{"r", ",", "0", ",",
|
|
SubscriptBox["N", "c"]}], "}"}], ",",
|
|
RowBox[{"{",
|
|
RowBox[{"s", ",", "0", ",",
|
|
SubscriptBox["N", "c"]}], "}"}]}], "]"}]}]}], ")"}]}],
|
|
"]"}]}]}], "\[IndentingNewLine]",
|
|
RowBox[{
|
|
RowBox[{
|
|
SubscriptBox["\[Phi]L",
|
|
RowBox[{"A_", ",", "q_"}]], "[",
|
|
RowBox[{"m_", ",", "n_"}], "]"}], ":=",
|
|
RowBox[{"Table", "[",
|
|
RowBox[{
|
|
RowBox[{
|
|
SubscriptBox["L",
|
|
RowBox[{"A", ",", "q"}]], "[",
|
|
RowBox[{"m", ",", "n", ",", "r", ",", "s"}], "]"}], ",",
|
|
RowBox[{"{",
|
|
RowBox[{"r", ",", "0", ",",
|
|
SubscriptBox["N", "c"]}], "}"}], ",",
|
|
RowBox[{"{",
|
|
RowBox[{"s", ",", "0", ",",
|
|
SubscriptBox["N", "c"]}], "}"}]}], "]"}]}], "\[IndentingNewLine]",
|
|
RowBox[{
|
|
RowBox[{
|
|
SubscriptBox["\[Phi]L",
|
|
RowBox[{"A_", ",", "q_"}]], "[", "\[CapitalPhi]_", "]"}], ":=",
|
|
RowBox[{"Total", "[",
|
|
RowBox[{"Flatten", "[",
|
|
RowBox[{
|
|
RowBox[{"Table", "[",
|
|
RowBox[{
|
|
RowBox[{
|
|
RowBox[{"\[CapitalPhi]", "[",
|
|
RowBox[{"[",
|
|
RowBox[{
|
|
RowBox[{"m", "+", "1"}], ",",
|
|
RowBox[{"n", "+", "1"}]}], "]"}], "]"}],
|
|
RowBox[{
|
|
SubscriptBox["\[Phi]L",
|
|
RowBox[{"A", ",", "q"}]], "[",
|
|
RowBox[{"m", ",", "n"}], "]"}]}], ",",
|
|
RowBox[{"{",
|
|
RowBox[{"m", ",", "0", ",",
|
|
SubscriptBox["N", "c"]}], "}"}], ",",
|
|
RowBox[{"{",
|
|
RowBox[{"n", ",", "0", ",",
|
|
SubscriptBox["N", "c"]}], "}"}]}], "]"}], ",", "1"}], "]"}],
|
|
"]"}]}], "\[IndentingNewLine]",
|
|
RowBox[{
|
|
RowBox[{"\[Phi]L", "[",
|
|
RowBox[{"\[CapitalPhi]_", ",", "\[CapitalLambda]_"}], "]"}], ":=",
|
|
RowBox[{"N", "[",
|
|
RowBox[{
|
|
UnderoverscriptBox["\[Sum]",
|
|
RowBox[{"m", "=", "0"}],
|
|
SubscriptBox["N", "c"]],
|
|
RowBox[{"(",
|
|
RowBox[{
|
|
UnderoverscriptBox["\[Sum]",
|
|
RowBox[{"n", "=", "0"}],
|
|
SubscriptBox["N", "c"]],
|
|
RowBox[{"(",
|
|
RowBox[{
|
|
RowBox[{"\[CapitalPhi]", "\[LeftDoubleBracket]",
|
|
RowBox[{
|
|
RowBox[{"m", "+", "1"}], ",",
|
|
RowBox[{"n", "+", "1"}]}], "\[RightDoubleBracket]"}],
|
|
RowBox[{"Table", "[",
|
|
RowBox[{
|
|
RowBox[{"\[CapitalLambda]", "\[LeftDoubleBracket]",
|
|
RowBox[{
|
|
RowBox[{"m", "+", "1"}], ",",
|
|
RowBox[{"n", "+", "1"}], ",",
|
|
RowBox[{"r", "+", "1"}], ",",
|
|
RowBox[{"s", "+", "1"}]}], "\[RightDoubleBracket]"}], ",",
|
|
RowBox[{"{",
|
|
RowBox[{"r", ",", "0", ",",
|
|
SubscriptBox["N", "c"]}], "}"}], ",",
|
|
RowBox[{"{",
|
|
RowBox[{"s", ",", "0", ",",
|
|
SubscriptBox["N", "c"]}], "}"}]}], "]"}]}], ")"}]}], ")"}]}],
|
|
"]"}]}]}], "Input",
|
|
CellChangeTimes->{{3.8333759184285297`*^9, 3.8333759188814993`*^9}, {
|
|
3.8333764326110353`*^9, 3.833376461246565*^9}, {3.833376496239792*^9,
|
|
3.833376512901948*^9}, {3.8333765504422493`*^9, 3.833376626171432*^9}, {
|
|
3.8333766645571127`*^9, 3.833376722534096*^9}, {3.833377023407977*^9,
|
|
3.833377049285466*^9}, {3.833449482195262*^9, 3.833449499000392*^9}, {
|
|
3.833449808794199*^9, 3.8334498110908213`*^9}, {3.8334499720213223`*^9,
|
|
3.833449988045641*^9}, {3.8334511689858828`*^9, 3.833451200828492*^9},
|
|
3.833451275308282*^9, {3.833451422211905*^9, 3.833451451429522*^9}, {
|
|
3.833457211761485*^9, 3.833457212652108*^9}},
|
|
CellLabel->
|
|
"In[152]:=",ExpressionUUID->"af8b7ff6-d757-4322-aa94-8cb2c99b4992"],
|
|
|
|
Cell[BoxData[
|
|
RowBox[{
|
|
RowBox[{"\[Lambda]func", "[", "q_", "]"}], ":=",
|
|
RowBox[{"Block", "[",
|
|
RowBox[{
|
|
RowBox[{"{",
|
|
RowBox[{
|
|
"\[CapitalPhi]0", ",", "\[CapitalPhi]1", ",", "\[CapitalLambda]"}],
|
|
"}"}], ",", "\[IndentingNewLine]",
|
|
RowBox[{
|
|
RowBox[{"\[CapitalPhi]0", "=",
|
|
RowBox[{"Table", "[",
|
|
RowBox[{
|
|
RowBox[{"If", "[",
|
|
RowBox[{
|
|
RowBox[{
|
|
RowBox[{"n", "\[Equal]", "0"}], "&&",
|
|
RowBox[{"m", "\[Equal]", "0"}]}], ",", "1", ",", "0"}], "]"}], ",",
|
|
RowBox[{"{",
|
|
RowBox[{"m", ",", "0", ",",
|
|
SubscriptBox["N", "c"]}], "}"}], ",",
|
|
RowBox[{"{",
|
|
RowBox[{"n", ",", "0", ",",
|
|
SubscriptBox["N", "c"]}], "}"}]}], "]"}]}], ";",
|
|
"\[IndentingNewLine]",
|
|
RowBox[{"\[CapitalLambda]", "=",
|
|
RowBox[{"N", "[",
|
|
SubscriptBox["L", "q"], "]"}]}], ";", "\[IndentingNewLine]",
|
|
RowBox[{"For", "[",
|
|
RowBox[{
|
|
RowBox[{"i", "=", "0"}], ",",
|
|
RowBox[{"i", "<", "52"}], ",",
|
|
RowBox[{"i", "++"}], ",",
|
|
RowBox[{"(",
|
|
RowBox[{"\[CapitalPhi]0", "=",
|
|
RowBox[{"\[Phi]L", "[",
|
|
RowBox[{"\[CapitalPhi]0", ",", "\[CapitalLambda]"}], "]"}]}],
|
|
")"}]}], "]"}], ";", "\[IndentingNewLine]",
|
|
RowBox[{"\[CapitalPhi]1", "=",
|
|
RowBox[{"\[Phi]L", "[",
|
|
RowBox[{"\[CapitalPhi]0", ",", "\[CapitalLambda]"}], "]"}]}], ";",
|
|
"\[IndentingNewLine]",
|
|
RowBox[{
|
|
FractionBox["\[CapitalPhi]1", "\[CapitalPhi]0"], "\[LeftDoubleBracket]",
|
|
RowBox[{"1", ",", "1"}], "\[RightDoubleBracket]"}]}]}], "]"}]}]], \
|
|
"Input",
|
|
CellChangeTimes->{{3.833457257689906*^9, 3.833457332822703*^9}, {
|
|
3.8334574139515743`*^9, 3.8334574309823647`*^9}, 3.833459172551968*^9},
|
|
CellLabel->
|
|
"In[156]:=",ExpressionUUID->"f0aa2bab-3547-4b25-8ec6-841323d25f13"],
|
|
|
|
Cell[BoxData[
|
|
RowBox[{
|
|
RowBox[{"secant", "[",
|
|
RowBox[{"x0_", ",", "y0_", ",", "x1_", ",", "y1_", ",", "z_"}], "]"}], ":=",
|
|
RowBox[{"x0", "-",
|
|
RowBox[{
|
|
RowBox[{"(",
|
|
RowBox[{"y0", "-", "z"}], ")"}],
|
|
RowBox[{"(",
|
|
FractionBox[
|
|
RowBox[{"x1", "-", "x0"}],
|
|
RowBox[{"y1", "-", "y0"}]], ")"}]}]}]}]], "Input",
|
|
CellChangeTimes->{{3.833451590804988*^9, 3.833451639877801*^9}, {
|
|
3.833457585325695*^9, 3.833457589181333*^9}},
|
|
CellLabel->
|
|
"In[157]:=",ExpressionUUID->"e43a65db-cc5c-4bc2-9396-e7ad73829ad3"],
|
|
|
|
Cell[BoxData[
|
|
RowBox[{
|
|
RowBox[{"secantmethod", "[",
|
|
RowBox[{
|
|
"f_", ",", "z_", ",", "x1_", ",", "x2_", ",", "\[Epsilon]_", ",",
|
|
"iterations_"}], "]"}], ":=",
|
|
RowBox[{"Block", "[",
|
|
RowBox[{
|
|
RowBox[{"{",
|
|
RowBox[{"k1", ",", "k2", ",", "k3", ",", "y1", ",", "y2", ",", "count"}],
|
|
"}"}], ",", "\[IndentingNewLine]",
|
|
RowBox[{
|
|
RowBox[{"k1", "=", "x1"}], ";", "\[IndentingNewLine]",
|
|
RowBox[{"k2", "=", "x2"}], ";", "\[IndentingNewLine]",
|
|
RowBox[{"y1", "=",
|
|
RowBox[{"f", "[", "k1", "]"}]}], ";", "\[IndentingNewLine]",
|
|
RowBox[{"y2", "=",
|
|
RowBox[{"f", "[", "k2", "]"}]}], ";", "\[IndentingNewLine]",
|
|
RowBox[{"count", "=", "0"}], ";", "\[IndentingNewLine]",
|
|
"\[IndentingNewLine]",
|
|
RowBox[{"While", "[",
|
|
RowBox[{
|
|
RowBox[{
|
|
RowBox[{
|
|
RowBox[{"Abs", "[",
|
|
RowBox[{"y1", "-", "z"}], "]"}], ">", "\[Epsilon]"}], "&&",
|
|
RowBox[{"count", "<", "iterations"}]}], ",", "\[IndentingNewLine]",
|
|
RowBox[{
|
|
RowBox[{"k3", "=",
|
|
RowBox[{"secant", "[",
|
|
RowBox[{"k1", ",", "y1", ",", "k2", ",", "y2", ",", "z"}], "]"}]}],
|
|
";", "\[IndentingNewLine]",
|
|
RowBox[{"k1", "=", "k2"}], ";", "\[IndentingNewLine]",
|
|
RowBox[{"y1", "=", "y2"}], ";", "\[IndentingNewLine]",
|
|
"\[IndentingNewLine]",
|
|
RowBox[{"k2", "=", "k3"}], ";", "\[IndentingNewLine]",
|
|
RowBox[{"y2", "=",
|
|
RowBox[{"f", "[", "k2", "]"}]}], ";", "\[IndentingNewLine]",
|
|
"\[IndentingNewLine]",
|
|
RowBox[{"count", "=",
|
|
RowBox[{"count", "+", "1"}]}], ";", "\[IndentingNewLine]",
|
|
RowBox[{"Print", "[",
|
|
RowBox[{"{",
|
|
RowBox[{"count", ",", "k1", ",", "y1"}], "}"}], "]"}], ";"}]}],
|
|
"\[IndentingNewLine]", "\[IndentingNewLine]", "]"}], ";"}]}],
|
|
"\[IndentingNewLine]", "]"}]}]], "Input",
|
|
CellChangeTimes->{{3.833457602153057*^9, 3.8334581013528223`*^9}, {
|
|
3.8334581580506687`*^9, 3.833458162457984*^9}, 3.833458274570757*^9, {
|
|
3.833458430593354*^9, 3.833458431967393*^9}, {3.8334584935635223`*^9,
|
|
3.833458521581132*^9}, {3.833458580962482*^9, 3.833458584101678*^9}, {
|
|
3.8334586878557177`*^9, 3.8334587505923033`*^9}, {3.83345904108678*^9,
|
|
3.833459050637452*^9}, 3.833459197693523*^9},
|
|
CellLabel->
|
|
"In[158]:=",ExpressionUUID->"84443b1d-b013-4655-81e3-152b9f45907d"],
|
|
|
|
Cell[CellGroupData[{
|
|
|
|
Cell[BoxData[
|
|
RowBox[{"secantmethod", "[",
|
|
RowBox[{
|
|
RowBox[{
|
|
RowBox[{"\[Lambda]func", "[", "#", "]"}], "&"}], ",", "1", ",", "1.3",
|
|
",", "1.31", ",",
|
|
SuperscriptBox["10",
|
|
RowBox[{"-", "20"}]], ",", "4"}], "]"}]], "Input",
|
|
CellChangeTimes->{{3.8334581450355663`*^9, 3.833458257040001*^9}, {
|
|
3.833458352159052*^9, 3.833458353133132*^9}, {3.833458441817086*^9,
|
|
3.833458447261086*^9}, {3.833458593108589*^9, 3.833458596697485*^9}, {
|
|
3.833458628329327*^9, 3.8334586753425503`*^9}, {3.833458752601663*^9,
|
|
3.833458829913876*^9}, {3.833458942692945*^9, 3.833458942754899*^9}},
|
|
CellLabel->
|
|
"In[159]:=",ExpressionUUID->"f95f417a-f1b8-48b0-aff2-4e1d3ba30c1c"],
|
|
|
|
Cell[CellGroupData[{
|
|
|
|
Cell[BoxData[
|
|
RowBox[{"{",
|
|
RowBox[{"1", ",", "1.31`", ",", "0.9938231242404644`"}], "}"}]], "Print",
|
|
CellChangeTimes->{3.83345908660977*^9, 3.833459254750489*^9,
|
|
3.833459555272402*^9, 3.8334606452685013`*^9},
|
|
CellLabel->
|
|
"During evaluation of \
|
|
In[159]:=",ExpressionUUID->"8cfa9596-8635-4deb-8559-014274eacf1d"],
|
|
|
|
Cell[BoxData[
|
|
RowBox[{"{",
|
|
RowBox[{"2", ",", "1.3065724836475656`", ",", "0.9999489377526637`"}],
|
|
"}"}]], "Print",
|
|
CellChangeTimes->{3.83345908660977*^9, 3.833459254750489*^9,
|
|
3.833459555272402*^9, 3.833460960959877*^9},
|
|
CellLabel->
|
|
"During evaluation of \
|
|
In[159]:=",ExpressionUUID->"62506904-a90c-46bc-bcc6-dd8c02c6f70c"],
|
|
|
|
Cell[BoxData[
|
|
RowBox[{"{",
|
|
RowBox[{"3", ",", "1.3065439132891326`", ",", "1.0000002220772635`"}],
|
|
"}"}]], "Print",
|
|
CellChangeTimes->{3.83345908660977*^9, 3.833459254750489*^9,
|
|
3.833459555272402*^9, 3.8334612775745773`*^9},
|
|
CellLabel->
|
|
"During evaluation of \
|
|
In[159]:=",ExpressionUUID->"0db1bbf7-ad36-4861-b0d2-678a09fa57d7"]
|
|
}, Open ]],
|
|
|
|
Cell[BoxData["$Aborted"], "Output",
|
|
CellChangeTimes->{3.83346143227352*^9},
|
|
CellLabel->
|
|
"Out[159]=",ExpressionUUID->"1358b658-b2ee-41e9-b4c2-bac808ea179f"]
|
|
}, Open ]],
|
|
|
|
Cell[CellGroupData[{
|
|
|
|
Cell[BoxData["\[IndentingNewLine]"], "Input",
|
|
CellChangeTimes->{{3.833459676746973*^9, 3.833459693261793*^9}, {
|
|
3.833461436640419*^9, 3.833461454673786*^9},
|
|
3.83346439676086*^9},ExpressionUUID->"c982c960-1df8-42e6-862e-\
|
|
6525c6293482"],
|
|
|
|
Cell[BoxData[
|
|
RowBox[{"SeriesCoefficient", "[",
|
|
RowBox[{
|
|
RowBox[{"Sin", "[", "x", "]"}], ",", "4"}], "]"}]], "Output",
|
|
CellChangeTimes->{3.833461455025072*^9},
|
|
CellLabel->
|
|
"Out[160]=",ExpressionUUID->"241c4e7a-5c6a-40df-b6ca-0ed2407b3a89"]
|
|
}, Open ]],
|
|
|
|
Cell[CellGroupData[{
|
|
|
|
Cell[BoxData[""], "Input",
|
|
CellChangeTimes->{{3.833463735672291*^9, 3.833463779711793*^9}, {
|
|
3.833463828292696*^9, 3.833463924942574*^9}, 3.833464270935576*^9, {
|
|
3.833464314261936*^9, 3.8334643230200853`*^9},
|
|
3.833464397852789*^9},ExpressionUUID->"559d4b36-36eb-476d-8986-\
|
|
55ef85a5aaa6"],
|
|
|
|
Cell[BoxData[
|
|
FractionBox[
|
|
RowBox[{
|
|
RowBox[{"-", "3"}], "+",
|
|
RowBox[{"2", " ", "x"}]}],
|
|
RowBox[{
|
|
RowBox[{
|
|
RowBox[{"-", "5"}], " ", "\[ImaginaryI]"}], "+",
|
|
RowBox[{"4", " ", "x"}]}]]], "Output",
|
|
CellChangeTimes->{{3.833463751888822*^9, 3.8334637802243357`*^9}, {
|
|
3.8334638294823227`*^9, 3.8334638317917*^9}, 3.8334639296836243`*^9,
|
|
3.833464272681696*^9, 3.833464334264291*^9},
|
|
CellLabel->
|
|
"Out[180]=",ExpressionUUID->"8b1ef0a9-765e-4869-8831-565b5674e9fb"]
|
|
}, Open ]],
|
|
|
|
Cell[CellGroupData[{
|
|
|
|
Cell[BoxData[""], "Input",
|
|
CellChangeTimes->{{3.833463785188533*^9, 3.83346378639987*^9}, {
|
|
3.833463822594265*^9, 3.833463838879908*^9}, {3.833463932523161*^9,
|
|
3.833463973059655*^9},
|
|
3.8334643990105267`*^9},ExpressionUUID->"d8004888-be25-4ebb-9a1a-\
|
|
7e6a7799816a"],
|
|
|
|
Cell[BoxData[
|
|
RowBox[{
|
|
RowBox[{"-", "0.48`"}], "+",
|
|
RowBox[{"0.4`", " ", "\[ImaginaryI]"}]}]], "Output",
|
|
CellChangeTimes->{{3.833463824356659*^9, 3.8334638417855797`*^9}, {
|
|
3.83346393714958*^9, 3.833463973391265*^9}, 3.833464277936246*^9,
|
|
3.833464336081966*^9},
|
|
CellLabel->
|
|
"Out[181]=",ExpressionUUID->"f002827f-7f21-45c6-837b-7ccf943398b8"]
|
|
}, Open ]]
|
|
},
|
|
WindowSize->{1440, 851},
|
|
WindowMargins->{{0, Automatic}, {Automatic, 0}},
|
|
TaggingRules->{
|
|
"WelcomeScreenSettings" -> {"FEStarting" -> False}, "TryRealOnly" -> False},
|
|
FrontEndVersion->"12.2 for Mac OS X x86 (64-bit) (December 12, 2020)",
|
|
StyleDefinitions->"Default.nb",
|
|
ExpressionUUID->"af9a09c9-df62-422c-a988-a490d98ae1b6"
|
|
]
|
|
(* End of Notebook Content *)
|
|
|
|
(* Internal cache information *)
|
|
(*CellTagsOutline
|
|
CellTagsIndex->{}
|
|
*)
|
|
(*CellTagsIndex
|
|
CellTagsIndex->{}
|
|
*)
|
|
(*NotebookFileOutline
|
|
Notebook[{
|
|
Cell[558, 20, 286, 6, 30, "Input",ExpressionUUID->"bf2e892a-0008-4529-bacd-4d3619ac745e"],
|
|
Cell[847, 28, 201, 4, 35, "Text",ExpressionUUID->"a2a33caa-af3c-4a3a-8054-70d79e05c1fa"],
|
|
Cell[1051, 34, 1607, 57, 181, "Input",ExpressionUUID->"be9cc703-c8ea-4b90-862f-061a2b1498c5"],
|
|
Cell[2661, 93, 238, 4, 35, "Text",ExpressionUUID->"48ecbeab-626b-4ed1-bffb-20691c3748f2"],
|
|
Cell[2902, 99, 1572, 47, 156, "Input",ExpressionUUID->"86a5d30e-95b8-493a-a4bf-7613b620acfa"],
|
|
Cell[4477, 148, 3608, 100, 123, "Input",ExpressionUUID->"cae51aa1-774a-45bc-8ddb-2d9d2bac6186"],
|
|
Cell[8088, 250, 213, 4, 35, "Text",ExpressionUUID->"aa6c1252-995d-4016-89e4-9cfcf4fb8ba5"],
|
|
Cell[8304, 256, 2975, 80, 211, "Input",ExpressionUUID->"746ed5e5-49cf-4db7-8474-08917e9d577c"],
|
|
Cell[11282, 338, 9874, 229, 189, "Input",ExpressionUUID->"784ec4e3-0e47-4303-898e-2c424985901f"],
|
|
Cell[21159, 569, 1713, 50, 77, "Input",ExpressionUUID->"f8e74c1e-8b8a-4bd8-9de6-2e4d3a09a42c"],
|
|
Cell[22875, 621, 752, 21, 60, "Input",ExpressionUUID->"90715246-bc56-45e7-8672-f47995804b03"],
|
|
Cell[23630, 644, 5287, 147, 187, "Input",ExpressionUUID->"af8b7ff6-d757-4322-aa94-8cb2c99b4992"],
|
|
Cell[28920, 793, 1887, 49, 162, InheritFromParent,ExpressionUUID->"f0aa2bab-3547-4b25-8ec6-841323d25f13"],
|
|
Cell[30810, 844, 546, 15, 50, "Input",ExpressionUUID->"e43a65db-cc5c-4bc2-9396-e7ad73829ad3"],
|
|
Cell[31359, 861, 2404, 53, 430, "Input",ExpressionUUID->"84443b1d-b013-4655-81e3-152b9f45907d"],
|
|
Cell[CellGroupData[{
|
|
Cell[33788, 918, 686, 14, 33, "Input",ExpressionUUID->"f95f417a-f1b8-48b0-aff2-4e1d3ba30c1c"],
|
|
Cell[CellGroupData[{
|
|
Cell[34499, 936, 322, 7, 24, "Print",ExpressionUUID->"8cfa9596-8635-4deb-8559-014274eacf1d"],
|
|
Cell[34824, 945, 337, 8, 24, "Print",ExpressionUUID->"62506904-a90c-46bc-bcc6-dd8c02c6f70c"],
|
|
Cell[35164, 955, 339, 8, 24, "Print",ExpressionUUID->"0db1bbf7-ad36-4861-b0d2-678a09fa57d7"]
|
|
}, Open ]],
|
|
Cell[35518, 966, 159, 3, 34, "Output",ExpressionUUID->"1358b658-b2ee-41e9-b4c2-bac808ea179f"]
|
|
}, Open ]],
|
|
Cell[CellGroupData[{
|
|
Cell[35714, 974, 243, 4, 52, "Input",ExpressionUUID->"c982c960-1df8-42e6-862e-6525c6293482"],
|
|
Cell[35960, 980, 250, 6, 34, "Output",ExpressionUUID->"241c4e7a-5c6a-40df-b6ca-0ed2407b3a89"]
|
|
}, Open ]],
|
|
Cell[CellGroupData[{
|
|
Cell[36247, 991, 299, 5, 30, "Input",ExpressionUUID->"559d4b36-36eb-476d-8986-55ef85a5aaa6"],
|
|
Cell[36549, 998, 490, 13, 51, "Output",ExpressionUUID->"8b1ef0a9-765e-4869-8831-565b5674e9fb"]
|
|
}, Open ]],
|
|
Cell[CellGroupData[{
|
|
Cell[37076, 1016, 276, 5, 30, "Input",ExpressionUUID->"d8004888-be25-4ebb-9a1a-7e6a7799816a"],
|
|
Cell[37355, 1023, 358, 8, 34, "Output",ExpressionUUID->"f002827f-7f21-45c6-837b-7ccf943398b8"]
|
|
}, Open ]]
|
|
}
|
|
]
|
|
*)
|
|
|