(* Content-type: application/mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 6.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 145, 7] NotebookDataLength[ 138739, 3779] NotebookOptionsPosition[ 125442, 3336] NotebookOutlinePosition[ 125955, 3356] CellTagsIndexPosition[ 125912, 3353] WindowFrame->Normal ContainsDynamic->False*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ RowBox[{" ", RowBox[{"(*", "\[IndentingNewLine]", " ", "\[IndentingNewLine]", " ", RowBox[{ "A", " ", "program", " ", "to", " ", "solve", " ", "the", " ", "truncated", " ", "Hamburger", " ", "moment", " ", "problem", "\[IndentingNewLine]", "\[IndentingNewLine]", " ", "Lawrence", " ", "Fialkow"}], "\[IndentingNewLine]", " ", "\[IndentingNewLine]", " ", "*)"}], "\[IndentingNewLine]", " ", "\[IndentingNewLine]", " ", "\[IndentingNewLine]", " ", RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{ "Compute", " ", "the", " ", "densities", " ", "and", " ", "nodes", " ", "\[IndentingNewLine]", " ", "for", " ", "a", " ", "one"}], "-", RowBox[{ "dimensional", " ", "real", " ", "truncated", " ", "power", " ", "moment", " ", RowBox[{"problem", ".", "\[IndentingNewLine]", "Moments"}], " ", "of", " ", "orders", " ", "0"}]}], ",", RowBox[{"...", "."}], ",", " ", RowBox[{ RowBox[{"m", " ", "must", " ", "be", " ", "pre"}], "-", RowBox[{"loaded", " ", "into", "\[IndentingNewLine]", " ", RowBox[{"bb", "[", "0", "]"}]}]}], ",", "...", ",", RowBox[{"bb", "[", "m", "]"}], ",", " ", RowBox[{ RowBox[{ RowBox[{"with", " ", RowBox[{"bb", "[", "0", "]"}]}], " ", ">", " ", RowBox[{"0", " ", RowBox[{ RowBox[{"(", RowBox[{"the", " ", "mass"}], ")"}], ".", " ", "The"}], " ", "program", " ", "will", "\[IndentingNewLine]", "determine", " ", "whether", " ", "or", " ", "not", " ", "there", " ", "exists", " ", "a", " ", "positive", " ", "Borel", " ", "measure", " ", "u", " ", "on", " ", "the", " ", "real", " ", "line", " ", "such", " ", "that", "\[IndentingNewLine]", "\[IndentingNewLine]", " ", RowBox[{"bb", "[", "i", "]"}]}]}], " ", "=", " ", RowBox[{ RowBox[{"Integral", "[", RowBox[{"x", "^", "i"}], "]"}], "du", " ", RowBox[{ RowBox[{"(", RowBox[{"0", "\[LessEqual]", "i", "\[LessEqual]", "m"}], ")"}], ".", "\[IndentingNewLine]", "\[IndentingNewLine]", "If"}], " ", "such", " ", "a", " ", "representing", " ", "measure", " ", "exists"}]}], ",", "\[IndentingNewLine]", " ", RowBox[{ RowBox[{ "then", " ", "the", " ", "program", " ", "will", " ", "compute", " ", "a", " ", "finitely", " ", "atomic", " ", "representing", " ", "\[IndentingNewLine]", " ", "measure", " ", "with", " ", "the", " ", "fewest", " ", "possible", " ", RowBox[{ "atoms", ".", "\[IndentingNewLine]", " ", "\[IndentingNewLine]", " ", "\[IndentingNewLine]", "If"}], " ", "m"}], " ", "=", " ", RowBox[{"2", "n"}]}], ",", " ", RowBox[{ RowBox[{ RowBox[{ "the", " ", "moment", " ", "data", " ", "will", " ", "fill", " ", "out", " ", "a", " ", "Hankel", " ", "matrix", " ", "\[IndentingNewLine]", " ", "H"}], " ", "\[Equal]", " ", RowBox[{"H", RowBox[{"(", "n", ")"}]}]}], ":=", " ", RowBox[{ RowBox[{"(", RowBox[{"bb", "[", RowBox[{"i", "+", "j"}], "]"}], ")"}], " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"0", "\[LessEqual]", "i"}], ",", RowBox[{"j", "\[LessEqual]", "n"}]}], ")"}], ".", "\[IndentingNewLine]", "If"}], " ", "H", " ", "is", " ", "positive", " ", "definite"}]}], ",", RowBox[{ RowBox[{ "there", " ", "are", " ", "infinitely", " ", "many", " ", "representing", " ", "measures"}], ";", " ", "\[IndentingNewLine]", " ", RowBox[{"to", " ", "fix", " ", "one"}]}], ",", " ", RowBox[{ RowBox[{ "an", " ", "extra", " ", "\"\\"", " ", "should", " ", "be", " ", "pre"}], "-", RowBox[{"loaded", " ", "into", " ", "\[IndentingNewLine]", " ", RowBox[{ RowBox[{"bb", "[", RowBox[{ RowBox[{"2", "n"}], "+", "1"}], "]"}], ".", " ", "An"}], " ", "auxiliary", " ", "column", " ", "vector", " ", "v", " ", "will", " ", "then", " ", "contain", " ", "\[IndentingNewLine]", " ", RowBox[{"bb", "[", RowBox[{"n", "+", "1"}], "]"}]}]}], ",", "...", ",", RowBox[{ RowBox[{ RowBox[{ RowBox[{"bb", "[", RowBox[{ RowBox[{"2", "n"}], "+", "1"}], "]"}], ".", "\[IndentingNewLine]", "\[IndentingNewLine]", "If"}], " ", "m"}], " ", "=", " ", RowBox[{ RowBox[{"2", "n"}], "+", "1"}]}], ",", " ", RowBox[{"the", " ", "moment", " ", "data", " ", RowBox[{"bb", "[", "0", "]"}]}], ",", RowBox[{ RowBox[{"...", RowBox[{"bb", "[", RowBox[{"2", "n"}], "]"}], " ", "will", " ", "fill", " ", "out"}], "\[IndentingNewLine]", " ", RowBox[{"the", " ", "Hankel", " ", "matrix", " ", "H", RowBox[{"(", "n", ")"}], " ", RowBox[{"(", RowBox[{"as", " ", "above"}], ")"}], "and"}], " ", "\[IndentingNewLine]", " ", RowBox[{"bb", "[", RowBox[{"n", "+", "1"}], "]"}]}], ",", "...", ",", RowBox[{ RowBox[{ RowBox[{"bb", "[", RowBox[{ RowBox[{"2", "n"}], "+", "1"}], "]"}], " ", "\[IndentingNewLine]", " ", "will", " ", "reside", " ", "in", " ", "the", " ", "auxiliary", " ", "column", " ", RowBox[{"v", ".", " ", "\[IndentingNewLine]", "If"}], " ", "m"}], " ", "=", " ", RowBox[{ RowBox[{"2", "n"}], "+", RowBox[{ "1", " ", "and", " ", "the", " ", "moments", " ", "correspond", " ", "to", " ", "Lebesgue", " ", "measure", " ", RowBox[{"on", " ", "[", RowBox[{"a", ",", "b"}], "]"}]}]}]}], ",", " ", "\[IndentingNewLine]", " ", RowBox[{ "the", " ", "resulting", " ", "nodes", " ", "and", " ", "weights", " ", "correspond", " ", "to", " ", "classical", " ", "Gaussian", " ", "\[IndentingNewLine]", " ", RowBox[{ "quadrature", ".", "\[IndentingNewLine]", " ", "\[IndentingNewLine]", " ", "The"}], " ", "method", " ", "we", " ", "use", " ", "is", " ", "based", " ", RowBox[{"on", " ", "[", "CF", "]"}], " ", RowBox[{"R", ".", "Curto"}], " ", "and", " ", RowBox[{"L", ".", "Fialkow"}]}], ",", "\[IndentingNewLine]", " ", RowBox[{"Houston", " ", RowBox[{"J", ".", " ", "Math", ".", " ", "17"}], RowBox[{"(", "1991", ")"}]}], ",", " ", RowBox[{"603", "-", "635"}], ",", " ", RowBox[{ RowBox[{"MR", " ", "93", RowBox[{"a", ":", "47016."}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"For", " ", "m"}]}], " ", "=", " ", RowBox[{"2", "n"}]}], ",", " ", RowBox[{ RowBox[{ "a", " ", "representing", " ", "measure", " ", "exists", " ", "if", " ", "and", " ", "only", " ", "if", "\[IndentingNewLine]", " ", "the", " ", "Hankel", " ", "H", RowBox[{"(", "n", ")"}], " ", "matrix", " ", "is", " ", "positive", " ", "semidefinite", " ", "and", " ", "\[IndentingNewLine]", " ", "\"\\""}], ";", " ", "equivalently"}], ",", " ", RowBox[{ RowBox[{"H", RowBox[{"(", "n", ")"}], " ", "admits", " ", "a", " ", "\[IndentingNewLine]", " ", "positive", " ", "semidefinite", " ", "Hankel", " ", "extension", " ", "H", RowBox[{ RowBox[{"(", RowBox[{"n", "+", "1"}], ")"}], ".", "\[IndentingNewLine]", "\[IndentingNewLine]", " ", "For"}], " ", "m"}], "=", RowBox[{ RowBox[{"2", "n"}], "+", "1"}]}], ",", " ", RowBox[{ RowBox[{ "a", " ", "representing", " ", "measure", " ", "exists", " ", "if", " ", "and", " ", "only", " ", "if", "\[IndentingNewLine]", " ", "H", RowBox[{"(", "n", ")"}], " ", "is", " ", "positive", " ", "semidefinite", " ", "and", " ", "v", " ", "belongs", " ", "to", " ", "Ran", " ", "H", RowBox[{"(", "n", ")"}]}], ";", "\[IndentingNewLine]", " ", "equivalently"}], ",", " ", RowBox[{ "the", " ", "partially", " ", "defined", " ", "matrix", "\[IndentingNewLine]", "\[IndentingNewLine]", " ", RowBox[{"(", " ", RowBox[{"H", RowBox[{"(", "n", ")"}], " ", "v"}], " ", ")"}], " ", "\[IndentingNewLine]", " ", RowBox[{"(", " ", RowBox[{ RowBox[{"v", "^", "T"}], " ", "*"}], " ", ")"}], "\[IndentingNewLine]", " ", "\[IndentingNewLine]", " ", "admits", " ", "a", " ", "positive", " ", "semidefinite", " ", "Hankel", " ", "completion", " ", "H", RowBox[{ RowBox[{"(", RowBox[{"n", "+", "1"}], ")"}], ".", "\[IndentingNewLine]", "\[IndentingNewLine]", " ", "The"}], " ", "columns", " ", "of", " ", "the", " ", "above", " ", "Hankel", " ", "matrices", " ", "are", " ", "denoted", " ", "\[IndentingNewLine]", " ", "1"}], ",", "t", ",", RowBox[{"t", "^", "2"}], ",", "...", ",", RowBox[{ RowBox[{ RowBox[{"t", "^", RowBox[{"(", RowBox[{"n", "+", "1"}], ")"}]}], ".", "\[IndentingNewLine]", " ", "\[IndentingNewLine]", " ", "When"}], " ", "a", " ", "measure", " ", "exists"}], ",", " ", RowBox[{ "we", " ", "find", " ", "a", " ", "minimal", " ", "i", " ", "for", " ", "which", "\[IndentingNewLine]", " ", "column", " ", RowBox[{"t", "^", "i"}], " ", "is", " ", "linearly", " ", "dependent", " ", "on", " ", "columns", " ", "to", " ", "its", " ", "left"}], ",", "\[IndentingNewLine]", "\[IndentingNewLine]", " ", RowBox[{ RowBox[{"t", "^", "i"}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"a0", " ", "1"}], " ", "+", " ", RowBox[{"a1", " ", "t"}], " ", "+"}], " ", "..."}], ".", " ", RowBox[{"+", " ", "a_"}]}], RowBox[{"{", "i", "}"}], " ", RowBox[{ RowBox[{"t", "^", RowBox[{"{", RowBox[{"i", "-", "1"}], "}"}]}], ".", "\[IndentingNewLine]", "\[IndentingNewLine]", "The"}], " ", "nodes", " ", "of", " ", "a", " ", "representing", " ", "measure", " ", RowBox[{"(", RowBox[{"with", " ", "the", " ", "fewest", " ", "nodes"}], ")"}], " ", "\[IndentingNewLine]", " ", RowBox[{ "consist", " ", "of", " ", "the", " ", "i", " ", "distinct", " ", "real", " ", "roots", " ", "of", " ", "the", " ", "polynomial"}], "\[IndentingNewLine]", "\[IndentingNewLine]", " ", RowBox[{"p", RowBox[{"(", "t", ")"}]}]}], " ", "=", " ", RowBox[{ RowBox[{"t", "^", "i"}], " ", "-", " ", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"a0", " ", "1"}], " ", "+", " ", RowBox[{"a1", " ", "t"}], " ", "+"}], " ", "..."}], ".", " ", RowBox[{"+", " ", "a_"}]}], RowBox[{"{", "i", "}"}], " ", RowBox[{"t", "^", RowBox[{"{", RowBox[{"i", "-", "1"}], "}"}]}]}], ")"}], ".", "\[IndentingNewLine]", " ", "\[IndentingNewLine]", "The"}], " ", "densities", " ", "of", " ", "the", " ", "measure", " ", "are", " ", "then", " ", "computed", " ", "using", " ", "a", " ", "Vandermonde", "\[IndentingNewLine]", " ", RowBox[{"matrix", " ", "based", " ", "on", " ", "the", " ", RowBox[{"nodes", ".", " ", RowBox[{"See", " ", "[", "CF", "]"}]}], " ", "for", " ", "complete", " ", RowBox[{"details", "."}]}]}]}]}]}]}], "\[IndentingNewLine]", " ", "\[IndentingNewLine]", "\[IndentingNewLine]", "*)"}], " "}]], "Input"], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"g", "[", "x_", "]"}], ":=", RowBox[{"Sin", "[", "x", "]"}]}]}]], "Input", CellChangeTimes->{{3.42115225078125*^9, 3.4211522514375*^9}, { 3.42115296884375*^9, 3.42115296921875*^9}, {3.42115335840625*^9, 3.421153369296875*^9}, {3.421153481625*^9, 3.421153487515625*^9}, { 3.422261704873625*^9, 3.422261722779875*^9}, {3.42226181320175*^9, 3.422261825654875*^9}, {3.4222619910455*^9, 3.4222619939205*^9}, { 3.422262104654875*^9, 3.422262120967375*^9}, {3.422262179108*^9, 3.422262180061125*^9}, {3.42226223001425*^9, 3.42226227538925*^9}, { 3.422262528858*^9, 3.42226258245175*^9}, {3.42226261720175*^9, 3.42226264151425*^9}, {3.422262696686125*^9, 3.4222627611705*^9}, { 3.422264873483*^9, 3.4222650056705*^9}, {3.422265066529875*^9, 3.42226510988925*^9}, {3.42226516913925*^9, 3.422265183029875*^9}, { 3.42226537838925*^9, 3.42226538438925*^9}, {3.422265434529875*^9, 3.42226544182675*^9}, 3.422265894983*^9, {3.422266508279875*^9, 3.422266560217375*^9}, 3.4222667097955*^9, 3.422267366092375*^9, { 3.422267496061125*^9, 3.422267528748625*^9}, 3.422354018248625*^9, { 3.42235409901425*^9, 3.422354231498625*^9}, {3.4223580858346252`*^9, 3.4223580932565002`*^9}, 3.4223587077721252`*^9, {3.42279220596875*^9, 3.422792207796875*^9}, {3.422792446734375*^9, 3.42279244946875*^9}, { 3.422793091859375*^9, 3.422793132453125*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{"choose", " ", "f"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"f", "[", "xx_", "]"}], " ", "=", "1"}], ";"}], "\[IndentingNewLine]"}]}]], "Input", CellChangeTimes->{{3.421148905671875*^9, 3.42114891371875*^9}, { 3.42115085215625*^9, 3.42115088084375*^9}, {3.42115094278125*^9, 3.4211509455625*^9}, {3.421150998484375*^9, 3.4211510215625*^9}, { 3.421153815578125*^9, 3.421153824515625*^9}, {3.421406399*^9, 3.4214064035625*^9}, {3.422266363404875*^9, 3.422266363779875*^9}, { 3.4223581016002502`*^9, 3.4223581021002502`*^9}, {3.4227917664375*^9, 3.422791831203125*^9}, 3.422794296829151*^9}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"h", "[", "rr_", "]"}], " ", "=", "1"}], ";"}]], "Input", CellChangeTimes->{{3.421148905671875*^9, 3.42114891371875*^9}, { 3.42115085215625*^9, 3.42115088084375*^9}, {3.42115094278125*^9, 3.4211509455625*^9}, {3.421150998484375*^9, 3.4211510215625*^9}, { 3.421153815578125*^9, 3.421153824515625*^9}, {3.421406399*^9, 3.4214064035625*^9}, {3.422266363404875*^9, 3.422266363779875*^9}, { 3.4223581016002502`*^9, 3.4223581021002502`*^9}, {3.4227917664375*^9, 3.422791846671875*^9}, {3.4227920576875*^9, 3.42279206134375*^9}, { 3.422792099515625*^9, 3.422792106375*^9}, 3.422794300063526*^9}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"h", "[", "2", "]"}]], "Input", CellChangeTimes->{{3.421148905671875*^9, 3.42114891371875*^9}, { 3.42115085215625*^9, 3.42115088084375*^9}, {3.42115094278125*^9, 3.4211509455625*^9}, {3.421150998484375*^9, 3.4211510215625*^9}, { 3.421153815578125*^9, 3.421153824515625*^9}, {3.421406399*^9, 3.4214064035625*^9}, {3.422266363404875*^9, 3.422266363779875*^9}, { 3.4223581016002502`*^9, 3.4223581021002502`*^9}, {3.4227917664375*^9, 3.42279183896875*^9}}], Cell[BoxData["1"], "Output", CellChangeTimes->{ 3.4227918158125*^9, 3.422791852671875*^9, {3.42279204396875*^9, 3.42279206575*^9}, 3.42279211765625*^9, 3.422792167921875*^9, { 3.422792457984375*^9, 3.42279249528125*^9}, 3.42279256290625*^9, 3.42279297975*^9, 3.422793101640625*^9, 3.4227931349375*^9, 3.422793218765625*^9, 3.422793380859375*^9, 3.422793505359375*^9, 3.422794029984375*^9, 3.422794194454151*^9, 3.422794314501026*^9, 3.422794584751026*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"f", "[", "2", "]"}]], "Input", CellChangeTimes->{{3.421151026453125*^9, 3.42115102921875*^9}}], Cell[BoxData["1"], "Output", CellChangeTimes->{{3.422262708467375*^9, 3.422262741092375*^9}, { 3.4222631977955*^9, 3.422263220529875*^9}, {3.422263288623625*^9, 3.4222633087955*^9}, 3.422263768858*^9, 3.422264068983*^9, 3.422264807248625*^9, {3.422264879029875*^9, 3.422264899998625*^9}, 3.422265454983*^9, 3.4222665134205*^9, 3.422266728186125*^9, 3.422354024498625*^9, 3.4223578172408752`*^9, 3.4223581087721252`*^9, 3.4223583115533752`*^9, 3.4223587161002502`*^9, 3.42279182121875*^9, 3.422791856234375*^9, 3.42279204621875*^9, 3.422792119875*^9, 3.422792169203125*^9, {3.4227924594375*^9, 3.4227924959375*^9}, 3.422792563328125*^9, 3.42279298034375*^9, 3.422793102140625*^9, 3.4227931355*^9, 3.422793219296875*^9, 3.422793381453125*^9, 3.422793507203125*^9, 3.422794030421875*^9, 3.422794194813526*^9, 3.422794315204151*^9, 3.422794585141651*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Plot", "[", RowBox[{ RowBox[{"g", "[", "x", "]"}], ",", RowBox[{"{", RowBox[{"x", ",", "0", ",", "10"}], "}"}]}], "]"}]], "Input", CellChangeTimes->{{3.42115298003125*^9, 3.42115298059375*^9}, { 3.422261735561125*^9, 3.422261745592375*^9}, {3.422262060498625*^9, 3.422262063998625*^9}, {3.422262319623625*^9, 3.422262320654875*^9}, 3.42226259788925*^9, {3.42226521288925*^9, 3.422265223983*^9}, { 3.422354049404875*^9, 3.422354078998625*^9}, {3.4223582972565002`*^9, 3.4223582978815002`*^9}, {3.422792222546875*^9, 3.422792222875*^9}, { 3.4227922606875*^9, 3.422792261328125*^9}}], Cell[BoxData[ GraphicsBox[{{}, {}, {Hue[0.67, 0.6, 0.6], LineBox[CompressedData[" 1:eJwt1nk0lN//AHBLZAszYynLzKCQilZFPt6XdimlSJIKSVkTRWlRIiRLZKmI CGVpUSJ6rCHJEiKSfZ+5dsb6e77n/P6Ze17n3Gfm3vf2jIKVs/F5Hi4uLk3y 43/r81tW/JtsAnUv/P9a8otpeIJ5EQyr1bKPk+vSt6/3mExPENNct1GKXJ8X BDH4mP4g2dJQ28LwB4f/3B2nGdHQgKJ2GJPrrzm78BFGKqQOnKg7JZ8Ky1s6 VViMHIjo2rTyqWQOxH+PuzzEqIBTzyJvvp8qh4xM/xYH9T/gVVlRLuTeDBeg NnrYtQueXIO13Hs7gSthKrclsAvqmi8kXDHvhKhlcq3fE7vAZGli6o9TJ5RX 2NJTG7tARC+b60FkJ6gen084r90NfTfT/zPv74R+O5W0Nu4eEDyW8Lrerwsu hd38WhPSC5YCehpV+d3g1LO2KytjAEYsKHt28/ZBu2WzQGHFAOgJZB5bLtUH x5oeqFd1D4B+HaU7R7UPdlT2efbIDMKQcslx1qE+4H2XJC7lNwghEYU8DyP7 IMpLQffamSEYyYxiPFfph0LaqmgtcRY4nR878lF7ALYGlX/ds44FFwqbXhgd HIBX/B7dR/eyIET9TdvvUwPwkNOoccmLBT/rOifeeQ2A6b/wbzH9LMjjkdkd mD8Ag6/FJ+YK2OBS4mFRvnMQhBR8VP7tGAGBL56JdLUhKNcJ1bi1ZwQslr2u ddw+BL5msdvljUfA+9d5sbTdQ8AT8nnfKfsR8DwS/nTScghmF4cv/H4+Aut1 31DzwoZgqNUkpZpnFIomYm+KTQ/BjyiVtQU/RsEt0j4g5v0wBGRt3XSmeRSM ikTi8/KHYX+NntZizyi0nz7V8KN8GEqWWxzQWRqF+3l3b2S3DUPe1ZCLnzeN wbrZ6aUsQRakH+ekvn0yBuaa2+kHLVjwSLxSLf7sOPC2f0nOHmdBznb5j4ed xkFa5L+yQwss6LZ0hvkb46C7d8N8FT8bdqbTjp+IHIcKuaSAh6vY0H/A4uaK n+Owc9brBrcuG3bdY1V77pyAHBE+dsldNji9BnOVAxOwre72Xk4AG2JqQ7vr TSfA+HyHofRjNowwNDkarhOwqTrjoPhLNsTm3VLsTZmAk/xnnCcL2TAzKeZu LDUJjxL2uPnNsUFJ3mppSWkSErpUlx/nwXB4d5Z/+qZJUAnJGREUxJAUZhYr cGgStMdNsjZJYTDWiC/7em8SMg+Dh4sGhlum40cdQifh6N02bb9tGFJv7mld FTcJZ2wL1P12YuD6MTDiljsJfwR/Renvw5But1lm3egkZIgELL06jaEp2Cex aXESxjyZDfzWGJZlN6r7ikyBgn9OrpEdBnO+G7s6VKZgICD01TNXDAIJxQ6R llNwhdMYpeiLYUuF5PRuhylwFhX2/R6AwXLkgveY5xTYXt1seyoYw0ddkSeH IqagH2hj26MwWLUcJ3h/TEH5moa2lhQMQk87HC83T8Hx+bz2z2kY3ps7yf3r nYJVXDNht99i4P3j55nLPQ1Hw34tVGRjSGzK3XJ5+zQcbr7/yrgUg2HUvs62 3dOgK7PY5VWOYeJEfYih8TR4fdX+GlKJYc9vFkvFcRqqVc4j91oMvQ3MlLaE aRBYTNff1IrhUUS6qeHbabC9LbQhpQ2Dpok2X27+NFRUBiUKdGDwqz9mFdE0 DS+22m717cGwMbxdnLd3Gt6Ehx542UfG75gj4TI+DasPhPSkDmBQ/eUrZyg6 A0pHrNKusDDUhNEqc2RnINU/oF4TY/AwfuGpsnYGbvdp3OoawcCkblCN0JyB oWurUzzHMFTU5jTy7J4BpGhwlDOOwTV0732XozOQnBzmYTWJQebory1tljPw 9+WS3McpDJdqhkNyPGcg5/OzZikOBmqIJ6j4zUD9McbImlkMuUb87PDwGfht b3+HPkfGX+zxM56EGeAfqvXnmscgXM046JI5A8XRH4R/kM56lMb5mzcDvAxv Pu8FDBaHtVIOfp+BDK4vVxmLZH2IfjPN+T0Dcy9HLZNJp1UZ86n0zABxqDR/ 5RKG40H/PoSPzcA7xd6wK6TnDR2seLg4EKci0ZZNOlGEI+6yggOxeoy4HtKH ftwn/spwYL1mX8M86alAqtNBVQ7seXPq7gLp2INxcjnbOOCs4pXYR3qv8PpK 5V0c+EIx0/lCmv39s2f4EdLvRA08SD8J2KPKY8kBO43qKjppXYO6Rmd7DvA/ birLIM/bK3jm/l8PDkyWeO5QJR1cMbTloC8HrgVRmIHkfXf4e3R+fswBt3DR 27/JeLTv5wtVjufAmnPYZAXpBwJhEJ7Bga8n9ZPWk/HcWE5nc+dxQLA4yXEb Ge8mvzfPnCs4MEuUZaqQ+bizb8fBv40c2HDU1oF3hqyP5aUcg24OXDhulFhO 5q/m29GUz6Mc4PE5a3aNzK+Hb5up8hIHjK/6+YtNkPWx154vXGQW5p8V7ggl 6+NyqY+Vs8osHAnLkjxK1hO9ey//1JZZoIRe2pZB1lslj+BrLzQLe2KfFE8P YViDgkb9T85CKhGscaQfQ52lUYS47SzI7IwRO9OL4fZNilak6yzIhUVYmXdj aM6NuJMYOAvb3E/bSrVj8G0+sXpd1Cwk7S5b9fsv2e8zq8rfJc7CFzeLPT4t GIK2xYoR+bOwqF4mlNiIQe/tq9hmPAtLd509OWT/sX/a6Z+dn4WNC+9W0isw PGWp9fYKzEHqtl86Gt8wTKplbphQmIPycLaMdAE5v5Ky80SPz8FP2k/rhfcY TEs8z0acnYMnQ9OWvpnkPOjauUzOcQ4CZjLZc28wnGYWHFzrS35fiU99WiIG ytPyP7s+z4FJxGXOmggMX3MCbn0vmQPJdZeaaKEY7JsMFY/WzsEYPcaS/RBD qWTtRcvBOcjXabx00QeDZ0jzjIfcPKSer/1xjZx3HfcHpTLuzIP+rY3qMwbk vEhMy90aNA++r65UvtuDYWexk+WX6HkYeKuga4YwRHCNJ5W/nwfLjL639uQ8 PnBjbmtX9zxkC1WppdMxvHMRObbywAJQhMcDn7PY8COG5pluugBSjyXionvZ 0FsiE6dvswBXBi463v7HBplVa4ccbi9AqobWEYFaNvgU7rlX+HEBuq6+KRT7 wAZT6u0PlxQXIU7kqOKsK/ke1/FtXtJYBG2n5svhl9gQYBu0FP7fItQWTjas tGIDkfv0IGG2COZsx4l/R9igav25ixa8CPtq7bTsNrBh9sMoLX9uEUYOugXb drJAom1Gy1hwCeYntqomNbFAXYDrbJ/UEigmZen+/MkCawvRNMrmJQiLP83d nsuCH8vW7ba1WwL9O5nJPGEsiDWxcRNrWILo/PgFAx0WxN+wH9x4ggvdQY4y DO9h2FggOON3mgux/z18Gew+DAXLUvj+WXOhrdathuyLw9Ae1M0McuFCDcqH 5a8aDwP9xekTA/5cSLmkxOet4jDElBiVxOdxocaxUpHfxBCErdgaS1XkRnLp UVqKrEG4Gzd/dGKYG0XXv6Y/khyAIxs/X143zo1SS26USvEPAL3wSqgVhxt9 +hd1NGSqH3I7B2tq+HgQNUwHTH/3w9iaZqN0Og8q/l3uoR7VD1ZpHw/ZHuFB KzTcy1ZK94NejtOBpg88SPSlzudVQn3A9asd5V3nRVPjKevEy7pB2c3W47ga H8pty0za6dYG3ZORpsL5/EjT4sm3uLY6MOi5dDPssgDq+37Rmcn/CQ53h6yj 04UQ2zRFdsGZIAaLre6Upwujj988Wn/x1xPqlEb5XWdWIN5anV1l69uJ3qDe uxHnVyCjNT79j/XbiedC03199ivQqI6dralZOyGybOX7hx4rkEbR+49599qJ wamTextDV6A0PQeVFX/aiaTWv04XS1aglUdSIt75dBByKT1E8FpRtA/fwm9q OwkhmDzXOi6KZEqzb0Xt6yE89FQFPs2Koms/3gQtt+ghenedygjmFkMj2S7j ti49RNH+wll9MTH0efSz31RUD3HdOOhxqpoY4s1v3XJioIcYOr+m9Oo5MZQ6 +KqU4tdL/HhookqpFkOl6n3Pln/qI7SDH/wcbBBDoYbo+vnyPiIl9ItbSasY er09uO/dnz7C54lC4bVBMaT8yUVHaqmPoNdMrcc84sjZoshceF8/0dhcVGpC F0dP7mmPPP3VT+xjmU8rmoij5+O8FxfaBohYiXu/bpuJI9nvZy9XsAeIiZ1v MltPiaO1676vubc4QLwImLvwxEocGfaY2xLyg8SsyrMmQRdxxLp8jWp+apBI s/qbMxIgjna3fIkX/TVIiDed8fpaII6ycg9nrssaImyX/MxkS8TRzxc33kYV DhF5ym+3epSJo9RT32ljP4eIi+5crE0/xdHibnrLxYEhoogafzqpRRzx9XZc cJAfJtwOdeo+nBJHn/Y1en/0Hiaaimy4zddTUIKKQZvMThYxovdH5qYGBfE4 Gejs38MilhcabX2xmYJWOqFoayMWoUloX+jdQUELf7S321iziPAv4lWueyio ZoJm9jCARRhl5UUFWlJQR7/etGsDi7Ddsvl9xjkKGpqZdb3ZxiJuvU+urLWh oP7NjsHufSwi7W3YorQ9BZ3cci1cncMihNLtbBKvURCfrOv2Ylk28S1JYmN+ CAV5mTbudz3NJtpWBxxof0xBL+6NGx8/zyYmXy5Z8UaS+5ObXys5sonVCYMR B55TkL9Ww8YALzZxN7ZgriGVgjrf/aVcjWETupEOZewiCmIzizaW1rEJU6nO dso3Cio00xnxbGYTjhEnZrdWUJDFW9FPsu1s4tlj/Q03qilI1fDqbmUWm5gN Xvl4eSsFmRdz5k35MfHxQYklc5KC9D8HPlDahokUg8qVfjNkfPrCGuR2YiJG pK6ONUdBu7WrgpYQJu6E/Nv7hYeKBNg96t6GmDCMmlU3Faci3Z/3C9ysMAHm 3AN5NCp6KR02/vkCJjbLCbxUkqaifwYhl3sdMCH9QlJ6VJ6Kth31hIWrmOhM 3rgUuI7cn6T5dVUAJhoubs8ZU6eiK2ujupIfYaJ8ne6Vk5upiF+idDvjMSbS Mw/2KWtR0ePhV2cqnmLCI9u2unAvFUlv33hCOA0T9p6O/qoGVKQgxnjWmYkJ y51uu4IPUZF4ZN10wgdM7CK8sy2OU9HW0q/9bbmYEC17Fjd9joqS71fZ7SrH BLf/y5OW58nfawuWta7ExITBa1qpHRWdO7a11v4nJpp/ZvuFOVORYqOjzo56 TPwI+arHcaWioMha1lwjJr4al86duUpFfB1Z0SnNmEhsrHPecJOKYu8dbs5r w0RkVPPa8DtUdDXX7OLqDkwEmLd3zd4j77s9YNitCxPObawT5QFUdEuoqLSm DxNWLyYoGo+oyD0sfmXrACaOW81VRoRS0atv5mY1Q5jQ7hVA1lFUtFGqJOoK xsSGFLHZiqdUJNla90RxFBPMS1JZG+OoSDPJwyt3DBPU9fJOkQlUZO7zaP/O CUzwsZVUF5OoSBQtchInMTGTqdZpk0pFtseKgqenMDF4edOzyjQqymP3CGyZ wcTfLTtMN7+lovMzPjYnOJiomdQVj/5ARUYrfsTbzGKiOHvP96VPVPTCe4Yw n8PEJ09DH9tcKtI7Z/ZVa56sr53HdKvyqeiUi8Ez7gVMPF04ObOlkIrq1ume fE/6EXH2fUwJFaVWfJg4uEjWm/cFB+5yKrLU0nH6SfrKLidlu0oq+m+1Z6nW EiZs+dzbf/6kokSW3GwQaRlGxe+5Wir604o+lZH+uUO+WrWBigq81uxjkb5r fPmbSRMVNTE+pC2S1nQozb/bQkUS1W29c6QH76/6mNlGRa0HXk92k46Nc0xr 7aAiQf7VzTmkjXMKXwr2UFGPxoXQ66T5f0k+1ewn6+v1U7oK6dzhi2HWQ1Tk vGLqDkGe14n/q38Im4p4U7M/7yKtyKR6549SUbm3+rcP5H0btWw9BieoyMv4 Y6oY6YBjuc7SM1TkFPjB+gQZL11H0Qu758h+MX4+EUDGc8zXyvLyIhXZe49b vCbjnfTik0ksNw1ZpOx6/pHMh1mu0KHKZTT0Wkj1UxqZL5F6y90zy2kIx+1L CZ7GhNvy5VuMRWlIupgmu4rMt6rCKbXbFBpSLDj0nBjHRKt2pkKaBA15Vf2d O0rWy26nE+L8sjQU++zVUW2yvmb83izfQqchJfM6gxAWJtLilxbPKNDQZrr2 6jqyHiUaklk5KjTkoqPhtaofE907Z747bqWhaKuQzk//MBFlYlgUs52G3CRX L9n/JeeH84ucMm0aKjjL5F/RgomshP0pCno0VPj+ZDOV7Ccfwej7DYdoyP5M SeEGsh93KLG8eI7S0L+DWpb6ZL8O6+i5aRynoWe+jmP6pWR/uAxY+ZvTkKdJ pqAAgYnVv7XQfxfI8/XHrEt4h4mmkaDtly7RUOClrHypDEw8FOpUj3SkoaFZ YbNrr8l+/y9AfvQKDbUuvMyYT8BESeKf2cQ7NKSef0t5JTl/rF1vfBSJoaHO dIO74c5k/Wxu+7j2OQ1FwmOBrkuYEBpDn/a+oKHmxzG+dFtMvL3Mn33nFQ0t u33L3cICE/MuoZ8n3tPQDTkj14l9mHjinPyltZKGQgQvvImUxcRFDaE8zk8a QmInFLIlMaGDHfKk6sjnD8cWFYthostpc/6RJhoK3ZwakciLiY1O+V9LummI 8dixx2eQTVQ4/CpMW6Ch6gPtj4qzyPm/XrOogksC7ZPIvpedziZchqOKenkl kPKuns6oV2xCysGymClE+sHjUaUoNmFlP1ASLiWBHkgQlaeus4n5i0tlXhoS 6EjkETfTnWxC48K6n4bnJNDvB1290Zksou0iwzndRgJxSXqVeCSxiIcONHFR OwlEaVe7rf+URQxcnjta7SSBcuM4KU98WcRLr8qGo14SaIXp4HN9CxYhHWb/ 1zRSAi3+S8zN5GURS3lvhs9WSaCgiP2a0ruGiQwiLqiwRgLdvO40dlZzmLAo eqyuWC+BAt9YhT5ZO0zklN1w6fojgfgUdZ+Vig0TV+oOTpzvl0Cbu2P6t7QM EX19Q3OXeCVR1AqT+JNOQ0Q1db2wu5YkcszYFCfqO0j2l+PAUx1JxLPffMPQ tUEioyyjrAgkUfx5b8Psi4PEI69NPuJ7JdE3muw5+qFBwrBXc+GNsSTSut/d qkobJL7n6OFOe0nE2KMuvP3ZAFF69sSvI7GSiDdG+q1nQj/xJePu0/W8Uqil 9/F/xd69BL+hiUwkvxTSMWaMHHLqJY4OqERzC0kho+bPe8vNe4k+paonjeJS KCRCctXDLb0ELVo67A5dCq1vEug+2t1D2N9Le1CvJYXU3DpLO/V7CFmzRncv Fyl0WtTqjfdEF3GDW+3Ij79SiH9y7wXpDR2EUOv1/U5Z0uiDy8mes6cbiZ5X CuFxXisRdc09yzvhJUS46qplOSarkE2pXwh/ayqExNiWnWPKoPXc8nRjz2ow Dz/9U6RPBtVICvt32LbAzO4kTtxXWWS3MqtQ+HIXeH2Ju1dSKIteHd61uODT BVybY1YMlMgifsZm9e6oLljOCFbYXCmL5Hxp2x8WdIEUx+NAyW9Z5C/y48kB sW7Ykm4Y3T8ii1a5b5pyet0NTpITOzYpySGq9tsVlN89MBbILjJRlkOVVrtv XRroAXeeAcPra+XQA8p6u09zPeCF/54t1pBDT4drrq9h9kJgRdkDEx05ZJKy 69wZu15Iufm0ydNEDuk66t3cMdULXT36HkV+cmiU57nw1uX9oB/3M0EtUA5p nKx3WLeqH16YmVeFPZJDkq+yjSnr+sGy8rKidYQcsrGJ80g53A9/3r74wftS Dql5B5TwRPRD3Y1F5p6vcogroniXBX0Aiii5FWUTcuiw9rtvhkqDoFC5Z0Jj hjyP3r+Q+c2DcMenlh41R+6PCTsZrT8IutP9V+x45JHXlQM50ecGIbdVmi4o Lo+utSCPrNhBeJ/s7mqwTh5RqReO75ccgoT/NslWnZNHcVyPHeLGh6DosNC9 ORt59P3zGq/ZpSHoPNM1uNZOHgW1FjbsEhkGpbsRub5O8qjPW3h//OphSPzG OYluyCPhWyuMfY8PQ5JRUdSHcHl0e1Dt0tjbYfh29ilXZ6Q8qj9t/0g9bxh6 L7vZiT+VR9c7U++blQ2DcrjyDsd4eRT9klXj/ncYXjUF/FbOkEcK1unSDQIs SDl3TCqmTB7941twY1qwoNx1/a3y7/JILUr4saYtC/rv8fVOVckjAdfAldou LFB9lf3xWL08yvwajgR9WJA6KGuyokMenfUyNy17zYLXV7of35mVRytWTxzM n2BB+n13yoUNdPT9ab+93XU2DOucvnt/Ix011L4K4nizYf3E7vGXW+gotoap 5OHPhtdWEg3tWnTUmLgq8r8oNqTAhyjzvXQUVHhUeSCLDYmcEcbhM3RUKv2S SQyyofttU4iDFR3JN8u6rRllg5JdAXfgeToKvJkz7THNhoTG4K4yezrysRa5 08eLIe6DerKeBx0FV21c3ymLIcbRQX1bKB1ZmSvVGuzH8Gf18bhj4XQ0ZOXp evcQhlWtO8VdI+lIs+CHcLIxhigD4bGM53R0aPvkqmwLDBEqrz+qvqajPvEY ozFnDPVtocr70uko41vjqVA3DBJPPCPPv6WjmC9H98t7YghbdsAz4RMdvfZc +jTnjSGko09HrpiORPPevagLw1ATVZ2u/Y2OnPz4bZqeYBA7kk0/WUFH38KK xr7HYHj01ZfrSTUd6dzvtLuagOHh0zUlYq10dIvufK33LYYfxiu2qf+jo7ob J02MszAIC00mGXbS0e5TxPDrbAwBHiV+/v10pLitFxS/YvAzsT7IO0lH8RuD 9NZ+x5Bg1GpxfoaOZHDZsGYVhrwDJk5lc3TU+cy6aUMNhrH/9oUG8jCQaZCx 3e8GDCI7ChJYfAy0tS9KIKAJg8pmrSwjQQZ6d065Y00LBgvldb9p4gzkI/d2 y6p2DNeYif3uNAZyqVF47dJJxkNGfva3FANlPCo9+q4bQ5momPwzeQb6d1tP nd2PoVPggfoik4GM31DODg5imOfhRmdXM9B/9bz51cMYNk+PWa1ex0CimaNj JiMYDo3au/mqM5AIt13F1CgGu6Hu+/2bGOhe2GTpnXEMd3tORxpsY6C9wXf6 xicwPP/XmJK2g4FCvJbWH5nC8LnZKFdUh4GytDwjIqYx1P0qr3QBBmrl4mOW zGBgVen9rdNnoFcXvle2cDAIlOeyt+5loOIz3FGtsxiUirZwRR5gIL/tY/e+ zWHQzUujcAwZ6GLM5KPoeQxmn9YonTrCQDdD3T6bLGC48jZ2a/4xBuq81780 S/rRa+m9jBMMJLe9xcZ/EUNKYsgJb3MGcq7u7eVawlAcK3ix6zQDran5eP8c 6baou9f3nGOgSllevTekOWFzgck2DGQ9Eb2yjbREkNtzQTsGCpwzEZklre7H yrC3Z6CwzUpyS6QPeNsWVDkxUPDswv/+b4PNjX+1Gq4MZC4xHlxC+pa7WVeo O/n8C8aED+lo59qJcQ8GsmS8c1MnnXXRgN/Ui4GiXXooheR5q62LpT/fZiCL X0I/dEgPntZZK3OPgQZ6br+MJ++3zOyjtpcveb6kD1GYjAfDWN2wzZ+BPB3X pCmT1jZMPo2CGIg2Y9O2l4zf8b1M54QQBuqxoG48QsbXGUXfWRbOQPPibfH6 ZPwDtKlhtpEMtK/3wVY6mZ+krYEvy2MYiFe4u7+TzGeB+rKParEMpCFYS4RM YphUnPrNTmSg5Ds7mpPHMIjLOw8cSWGg3LoyeQpZL2rS/bPv3zDQoQyqvzXG cFb4j/y19wzkffdzU/EQhut8xzSaPzLQB4k3RO0A2f9LlWhnDgOlnHtQV96H 4ft4vvUSwUC3+Xu9XbswbGuNT/WrYiAn9/EtWc3k76nbzWTWkP3SY0Fo/cYw fFt9X9MvBvJK/nHlTT2GRKUv3ap/GKhWV83yWDWZP/t6xvdecv/B23YfizGM 5MU4jQ0wUPMRW/tnBWS/i57Ll2Ex0J/xS0nO+Rh83rPM7ccZ6Aez+NufTxhm vXgP8nMxUe9BWU5xCob+JFvXGG4mKpfjwyJJGBp/VkSr8zKRaumODfrxGN4z Q/pN+ZnI6txC4OVoDJdK5HyTRZjIxnftyu4HGJqFNQv3r2KiNK13eTI2GL5t je5vlWEigd4wx7VnyHo4PS92WY6JLG8ZhK02J+dbZpFlNIOJuqYqdrUbYdh/ zGh+YA0TGTdHxTtqk/0WY7c9cDMT2QYX7u4SxvCquNKSsZWJejb/KcnlwxA+ rO77YRsT3d1zd+TWEhtcdCfrW3Yw0WFVKfOqMTaodnq7rgfy+TPLeu82sSFq 7dP0qoNM5E6l+B5/wYaJ8b1J9w4x0aYU3U3m5PvE6OvYMy0jJvrw1dzzQAgb +I8ZPEwyZqKLE4Iz3bfZ4O7FuXTzJBNFPX9gJm7JBuPqE6obLjDR9PJna3Sl 2JARzcvssmOid0EP/dpF2CBokykdfYmJrgtcH3bgZUPBDP9yPicmIqr4WnaN sEBD8VNPqxsTiYx4iA+Vs0DEXeLlw7tMtGaeR3fNVRbYQUGMvg8TKRSmbbS3 Z0GxoEPYzH0maqm69uzpWRZcjy32tvFnouyaNUtfDFgwUOZ6VieEiex8/6ub k2fBN5la+eHnTCQT68u6VTAMtwuCog7mMBEzM16qZnIIBv8WhHNymUhNqEct q38ITOfGQ5LzmOjcaskQn5YhWK9p7s9bQN4nwMt1sGAImt4oX//yjYn4zLfF VQYMweZI4tS6eiaa/VvPOywzBL2Oo3RhzET97hP0ofWDcEjWJOn7agXk+VB2 A8j0Q1S3esEpZQUUWbDeOFewH7rTBVqGVRQQJ8jvkzKnD7xQnrjYOgUUn1p/ /tfvPkg7r+R1bJMCElFqv3kzvA+E3o4at/6ngBwYL1LlhfugbPcjLrapAmKs NbVvYvWAnlPZKYq/AhLWfLI+ObAL4uSsVv+YVED+U9eL9R60wOri5CeHTymi tHe/dOe0qiFkQ5Fa3VdFpJHwKUA3+TXYji+rrVBTQjZzy+/e/K+EsHo3xvQP VkKKrAxW9L5Govy7wL28eSUUeGyfUbZOB+GwqyOkgvQnqty+C/91EOJ5ObGN pOstqDrSuh2EWfql3BHSP+Ws5a9BB9EX/GN09YISapMOy9fU7yD4TULOBpGW tBRx+7ivg9j1T1rXclEJxVkmfX5r3EH0nhg5aE/6Q4Xr8bPHOoiAmvKTHqTf SVuyxI53EHWFnu5hpC37gqScTToIq8SWtG+kbZSM9m4w6yDuXIyT1VhSQjW3 dGVTT3cQqzuvrdUhLWvb+sDMsoMoNz+y/QBptsWT8eVnyPMbch+zJn1xmWWB 7dkOIk7dOiCS9DlBre2rrTsI/eSdUYmkn+7MCvpFupch8eodafzYoOOuTQeh TiktrCTNr19+u+M8eT7/59VNpC+fzvweYttBXOW5+reHtKbaOyq60EHI3Dg8 NEZaLbDqBCb9dVyZs0RawXpZTKxdB/F/76Phhg== "]]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesOrigin->{0, 0}, PlotRange->{{0, 10}, {-0.9999993680039507, 0.9999995755730018}}, PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Scaled[0.02]}]], "Output", CellChangeTimes->{{3.422261727811125*^9, 3.4222617471705*^9}, { 3.422261815858*^9, 3.422261828233*^9}, {3.422261981654875*^9, 3.4222619967955*^9}, 3.422262027108*^9, 3.422262064623625*^9, { 3.422262108404875*^9, 3.42226212738925*^9}, 3.422262183248625*^9, { 3.422262235404875*^9, 3.422262322108*^9}, {3.42226253207675*^9, 3.422262644358*^9}, {3.42226271045175*^9, 3.422262764529875*^9}, { 3.42226319963925*^9, 3.42226322232675*^9}, {3.422263290279875*^9, 3.422263310279875*^9}, 3.422263483358*^9, 3.422263770498625*^9, 3.422264071748625*^9, 3.422264809108*^9, {3.42226488032675*^9, 3.42226500845175*^9}, {3.422265079654875*^9, 3.422265112108*^9}, { 3.42226517351425*^9, 3.422265186748625*^9}, {3.422265219217375*^9, 3.422265224686125*^9}, 3.42226538832675*^9, 3.422265456561125*^9, 3.422265902029875*^9, {3.422266515436125*^9, 3.422266564998625*^9}, { 3.42226671520175*^9, 3.422266728233*^9}, 3.422267372904875*^9, 3.42226741201425*^9, {3.422267501592375*^9, 3.422267533358*^9}, { 3.4223540255455*^9, 3.4223541552955*^9}, {3.4223541935455*^9, 3.422354257092375*^9}, 3.4223578190690002`*^9, 3.4223581096002502`*^9, { 3.4223582999752502`*^9, 3.4223583121158752`*^9}, 3.4223587169283752`*^9, 3.422791871484375*^9, 3.42279203528125*^9, 3.4227921236875*^9, 3.422792173359375*^9, {3.422792214125*^9, 3.42279222590625*^9}, 3.422792261859375*^9, 3.422792497*^9, 3.422792563828125*^9, 3.42279298090625*^9, {3.422793102640625*^9, 3.422793136125*^9}, 3.42279321996875*^9, 3.422793262734375*^9, 3.42279338209375*^9, 3.42279350809375*^9, 3.42279403096875*^9, 3.422794195344776*^9, 3.422794315954151*^9, 3.422794585626026*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{ "initialization", " ", "for", " ", "Gaussian", " ", "quadrature", " ", RowBox[{"over", " ", "[", RowBox[{"a", ",", "b"}], "]"}], " ", "of", " ", "degree", " ", "m"}], "=", RowBox[{ RowBox[{"2", "n"}], "+", "1"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"m", "=", "7"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"a", "=", " ", "1"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"b", "=", "2"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"c", "=", "7"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"d", "=", "8"}], ";"}], "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "0"}], ",", RowBox[{"i", "\[LessEqual]", "m"}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"bb", "[", "i", "]"}], " ", "=", " ", RowBox[{ RowBox[{"NIntegrate", "[", RowBox[{ RowBox[{ RowBox[{"x", "^", "i"}], " ", RowBox[{"f", "[", "x", "]"}]}], " ", ",", RowBox[{"{", RowBox[{"x", ",", "a", ",", "b"}], "}"}]}], " ", "]"}], " ", "+", " ", RowBox[{"NIntegrate", "[", RowBox[{ RowBox[{ RowBox[{"r", "^", "i"}], " ", RowBox[{"h", "[", "r", "]"}]}], " ", ",", RowBox[{"{", RowBox[{"r", ",", "c", ",", "d"}], "}"}]}], " ", "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"Print", "[", RowBox[{"bb", "[", "i", "]"}], "]"}]}]}], "]"}], "\[IndentingNewLine]"}]}]], "Input", CellChangeTimes->{{3.421148905671875*^9, 3.42114891371875*^9}, { 3.42115085215625*^9, 3.42115088084375*^9}, {3.42115094278125*^9, 3.4211509455625*^9}, 3.421150998484375*^9, {3.421151117796875*^9, 3.421151132140625*^9}, {3.421151188140625*^9, 3.421151188515625*^9}, { 3.421153126046875*^9, 3.4211531284375*^9}, {3.42115364790625*^9, 3.421153745171875*^9}, {3.42115378421875*^9, 3.42115379796875*^9}, { 3.42140674778125*^9, 3.421406759328125*^9}, {3.421408312875*^9, 3.421408315171875*^9}, {3.4222617846705*^9, 3.42226178526425*^9}, 3.42226318851425*^9, {3.42226329688925*^9, 3.422263300186125*^9}, 3.42226347563925*^9, 3.42226375888925*^9, {3.422264056998625*^9, 3.42226405738925*^9}, {3.42226416388925*^9, 3.422264167436125*^9}, { 3.422264232873625*^9, 3.4222642330455*^9}, 3.4222652336705*^9, 3.422265520983*^9, {3.422265578483*^9, 3.422265578623625*^9}, { 3.422265681123625*^9, 3.422265681248625*^9}, 3.422265918873625*^9, { 3.422266098108*^9, 3.42226609876425*^9}, 3.422266148279875*^9, 3.422266859654875*^9, 3.422266889998625*^9, {3.422266929217375*^9, 3.422266929842375*^9}, 3.422267630873625*^9, 3.422267702967375*^9, { 3.422267786358*^9, 3.42226778707675*^9}, {3.422267822108*^9, 3.422267822186125*^9}, {3.422267863904875*^9, 3.422267864108*^9}, { 3.422354281342375*^9, 3.422354282529875*^9}, {3.4223578516471252`*^9, 3.4223578524752502`*^9}, {3.4223579122721252`*^9, 3.4223579127721252`*^9}, {3.4223581660533752`*^9, 3.4223581664127502`*^9}, {3.4223582099440002`*^9, 3.4223582106471252`*^9}, {3.4223582475377502`*^9, 3.4223582799127502`*^9}, {3.4223584024752502`*^9, 3.4223584028971252`*^9}, {3.4223584400690002`*^9, 3.4223584404752502`*^9}, {3.4223584795690002`*^9, 3.4223584805846252`*^9}, {3.4223589404596252`*^9, 3.4223589417096252`*^9}, {3.4223589810846252`*^9, 3.4223589906627502`*^9}, {3.42279189478125*^9, 3.422791968578125*^9}, { 3.422792074390625*^9, 3.4227921095*^9}, {3.42279228565625*^9, 3.422792362265625*^9}, {3.42279320015625*^9, 3.422793210359375*^9}, { 3.422794016265625*^9, 3.42279401834375*^9}, {3.422794183079151*^9, 3.422794183954151*^9}, {3.422794574704151*^9, 3.422794576094776*^9}}], Cell[CellGroupData[{ Cell[BoxData["2.000000000000002`"], "Print", CellChangeTimes->{{3.422792353796875*^9, 3.422792362765625*^9}, 3.422792567015625*^9, 3.42279281675*^9, 3.422792997125*^9, 3.42279315946875*^9, 3.42279322509375*^9, 3.422793269484375*^9, 3.422793383109375*^9, 3.422793509203125*^9, 3.422794031859375*^9, 3.422794195907276*^9, 3.422794316860401*^9, 3.422794586282276*^9}], Cell[BoxData["9.00000000000001`"], "Print", CellChangeTimes->{{3.422792353796875*^9, 3.422792362765625*^9}, 3.422792567015625*^9, 3.42279281675*^9, 3.422792997125*^9, 3.42279315946875*^9, 3.42279322509375*^9, 3.422793269484375*^9, 3.422793383109375*^9, 3.422793509203125*^9, 3.422794031859375*^9, 3.422794195907276*^9, 3.422794316860401*^9, 3.422794586313526*^9}], Cell[BoxData["58.66666666666673`"], "Print", CellChangeTimes->{{3.422792353796875*^9, 3.422792362765625*^9}, 3.422792567015625*^9, 3.42279281675*^9, 3.422792997125*^9, 3.42279315946875*^9, 3.42279322509375*^9, 3.422793269484375*^9, 3.422793383109375*^9, 3.422793509203125*^9, 3.422794031859375*^9, 3.422794195907276*^9, 3.422794316860401*^9, 3.422794586360401*^9}], Cell[BoxData["427.5000000000005`"], "Print", CellChangeTimes->{{3.422792353796875*^9, 3.422792362765625*^9}, 3.422792567015625*^9, 3.42279281675*^9, 3.422792997125*^9, 3.42279315946875*^9, 3.42279322509375*^9, 3.422793269484375*^9, 3.422793383109375*^9, 3.422793509203125*^9, 3.422794031859375*^9, 3.422794195907276*^9, 3.422794316860401*^9, 3.422794586407276*^9}], Cell[BoxData["3198.4000000000037`"], "Print", CellChangeTimes->{{3.422792353796875*^9, 3.422792362765625*^9}, 3.422792567015625*^9, 3.42279281675*^9, 3.422792997125*^9, 3.42279315946875*^9, 3.42279322509375*^9, 3.422793269484375*^9, 3.422793383109375*^9, 3.422793509203125*^9, 3.422794031859375*^9, 3.422794195907276*^9, 3.422794316860401*^9, 3.422794586454151*^9}], Cell[BoxData["24093.00000000003`"], "Print", CellChangeTimes->{{3.422792353796875*^9, 3.422792362765625*^9}, 3.422792567015625*^9, 3.42279281675*^9, 3.422792997125*^9, 3.42279315946875*^9, 3.42279322509375*^9, 3.422793269484375*^9, 3.422793383109375*^9, 3.422793509203125*^9, 3.422794031859375*^9, 3.422794195907276*^9, 3.422794316860401*^9, 3.422794586501026*^9}], Cell[BoxData["181962.28571428594`"], "Print", CellChangeTimes->{{3.422792353796875*^9, 3.422792362765625*^9}, 3.422792567015625*^9, 3.42279281675*^9, 3.422792997125*^9, 3.42279315946875*^9, 3.42279322509375*^9, 3.422793269484375*^9, 3.422793383109375*^9, 3.422793509203125*^9, 3.422794031859375*^9, 3.422794195907276*^9, 3.422794316860401*^9, 3.422794586547901*^9}], Cell[BoxData["1.3765837500000016`*^6"], "Print", CellChangeTimes->{{3.422792353796875*^9, 3.422792362765625*^9}, 3.422792567015625*^9, 3.42279281675*^9, 3.422792997125*^9, 3.42279315946875*^9, 3.42279322509375*^9, 3.422793269484375*^9, 3.422793383109375*^9, 3.422793509203125*^9, 3.422794031859375*^9, 3.422794195907276*^9, 3.422794316860401*^9, 3.422794586610401*^9}] }, Open ]] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{"THMP", "[", "m_", "]"}], ":=", RowBox[{"(", " ", RowBox[{"(*", RowBox[{ RowBox[{ "Solve", " ", "the", " ", "Hamburger", " ", "moment", " ", "problem", " ", "for", " ", "moment", " ", "data", " ", "\[IndentingNewLine]", " ", RowBox[{"bb", "[", "0", "]"}]}], ",", "...", ",", RowBox[{"bb", "[", "m", "]"}]}], "*)"}], " "}]}], "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ "build", " ", "the", " ", "moment", " ", "table", " ", "from", " ", "the", " ", "moment", " ", RowBox[{"data", ":"}]}], " ", "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"n", "=", " ", RowBox[{"Floor", "[", RowBox[{"m", "/", "2"}], "]"}]}], ";"}], " ", "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "0"}], ",", RowBox[{"i", "\[LessEqual]", RowBox[{"n", "+", "1"}]}], ",", RowBox[{"i", "=", RowBox[{"i", "+", "1"}]}], ",", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", "0"}], ",", RowBox[{"j", "\[LessEqual]", RowBox[{"n", "+", "1"}]}], ",", RowBox[{"j", "=", RowBox[{"j", "+", "1"}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"Mtab", "[", RowBox[{"i", ",", "j"}], "]"}], "=", RowBox[{"bb", "[", RowBox[{"i", "+", "j"}], "]"}]}]}], "]"}]}], "]"}], ";"}], " "}]}]], "Input"], Cell[BoxData[""], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{"build", " ", "the", " ", "moment", " ", RowBox[{"(", "Hankel", ")"}], " ", "matrix", " ", "with", " ", "an", " ", "auxiliary", " ", "column", " ", "on", " ", "the", " ", RowBox[{"right", ".", " ", "In"}], " ", "the", " ", "case", " ", "m"}], " ", "=", " ", RowBox[{"2", "n"}]}], ",", " ", "\[IndentingNewLine]", " ", RowBox[{ RowBox[{ "this", " ", "extra", " ", "column", " ", "is", " ", "used", " ", "to", " ", "compute", " ", "a", " ", "representing", " ", "\[IndentingNewLine]", "measure", " ", "with", " ", "a", " ", "prescribed", " ", "moment", " ", "of", " ", "order", " ", "m"}], "+", RowBox[{ "1", " ", "\[IndentingNewLine]", "if", " ", "the", " ", "Hankel", " ", "matrix", " ", "is", " ", "positive", " ", RowBox[{"definite", ".", " ", "In"}], " ", "that", " ", "case"}]}], ",", " ", "\[IndentingNewLine]", " ", RowBox[{"the", " ", "free", " ", "choice", " ", RowBox[{"bb", "[", RowBox[{"m", "+", "1"}], "]"}], " ", "is", " ", "the", " ", "final", " ", "entry", " ", "in", " ", "the", " ", "auxiliary", " ", RowBox[{"column", "."}]}]}], " ", "*)"}], "\[IndentingNewLine]", " ", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"H", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Mtab", "[", RowBox[{"r", ",", "s"}], "]"}], ",", RowBox[{"{", RowBox[{"r", ",", "0", ",", "n"}], "}"}], ",", RowBox[{"{", RowBox[{"s", ",", "0", ",", RowBox[{"n", "+", "1"}]}], "}"}]}], "]"}]}], ";"}], " ", RowBox[{"(*", " ", RowBox[{ RowBox[{"Hankel", " ", "matrix"}], " ", "+", " ", RowBox[{"extra", " ", "column"}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", "\"\\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"m", "\[Equal]", RowBox[{"2", "n"}]}], ",", "\[IndentingNewLine]", RowBox[{"Print", "[", " ", RowBox[{"MatrixForm", "[", RowBox[{"H", "[", RowBox[{"[", RowBox[{ RowBox[{"Range", "[", RowBox[{"1", ",", RowBox[{"n", "+", "1"}]}], "]"}], ",", RowBox[{"Range", "[", RowBox[{"1", ",", RowBox[{"n", "+", "1"}]}], "]"}]}], "]"}], "]"}], "]"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Print", "[", RowBox[{"MatrixForm", "[", RowBox[{"H", "[", RowBox[{"[", RowBox[{ RowBox[{"Range", "[", RowBox[{"1", ",", RowBox[{"n", "+", "1"}]}], "]"}], ",", RowBox[{"Range", "[", RowBox[{"1", ",", RowBox[{"n", "+", "2"}]}], "]"}]}], "]"}], "]"}], "]"}], "]"}]}], "]"}], ";"}], "\[IndentingNewLine]"}]}]], "Input"], Cell[CellGroupData[{ Cell[BoxData["\<\"Moment matrix :\"\>"], "Print", CellChangeTimes->{3.422263344436125*^9, 3.422263494811125*^9, 3.422263788936125*^9, 3.422264079858*^9, 3.422264174436125*^9, 3.422264237561125*^9, 3.422265526686125*^9, 3.42226558413925*^9, 3.422265686358*^9, 3.422265925592375*^9, 3.422266103858*^9, 3.422266151154875*^9, 3.422266806998625*^9, 3.42226686238925*^9, 3.4222668929205*^9, 3.42226693338925*^9, 3.422267633404875*^9, 3.422267705123625*^9, 3.42226778876425*^9, 3.4222678242955*^9, 3.422267865717375*^9, 3.42235428726425*^9, 3.4223578584127502`*^9, 3.4223579204596252`*^9, 3.4223581796940002`*^9, 3.4223582131627502`*^9, 3.4223583168502502`*^9, 3.4223584066940002`*^9, 3.4223584450221252`*^9, 3.4223584874596252`*^9, 3.4223587290846252`*^9, 3.4223589451940002`*^9, 3.4223589949908752`*^9, 3.422792571265625*^9, 3.422793001625*^9, 3.422793161984375*^9, 3.422793230609375*^9, 3.422793386671875*^9, 3.42279351359375*^9, 3.422794035625*^9, 3.422794197766651*^9, 3.422794319204151*^9, 3.422794588172901*^9}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"2.000000000000002`", "9.00000000000001`", "58.66666666666673`", "427.5000000000005`", "3198.4000000000037`"}, {"9.00000000000001`", "58.66666666666673`", "427.5000000000005`", "3198.4000000000037`", "24093.00000000003`"}, {"58.66666666666673`", "427.5000000000005`", "3198.4000000000037`", "24093.00000000003`", "181962.28571428594`"}, {"427.5000000000005`", "3198.4000000000037`", "24093.00000000003`", "181962.28571428594`", "1.3765837500000016`*^6"} }, GridBoxAlignment->{ "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Print", CellChangeTimes->{3.422263344436125*^9, 3.422263494811125*^9, 3.422263788936125*^9, 3.422264079858*^9, 3.422264174436125*^9, 3.422264237561125*^9, 3.422265526686125*^9, 3.42226558413925*^9, 3.422265686358*^9, 3.422265925592375*^9, 3.422266103858*^9, 3.422266151154875*^9, 3.422266806998625*^9, 3.42226686238925*^9, 3.4222668929205*^9, 3.42226693338925*^9, 3.422267633404875*^9, 3.422267705123625*^9, 3.42226778876425*^9, 3.4222678242955*^9, 3.422267865717375*^9, 3.42235428726425*^9, 3.4223578584127502`*^9, 3.4223579204596252`*^9, 3.4223581796940002`*^9, 3.4223582131627502`*^9, 3.4223583168502502`*^9, 3.4223584066940002`*^9, 3.4223584450221252`*^9, 3.4223584874596252`*^9, 3.4223587290846252`*^9, 3.4223589451940002`*^9, 3.4223589949908752`*^9, 3.422792571265625*^9, 3.422793001625*^9, 3.422793161984375*^9, 3.422793230609375*^9, 3.422793386671875*^9, 3.42279351359375*^9, 3.422794035625*^9, 3.422794197766651*^9, 3.422794319204151*^9, 3.422794588188526*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{ "Check", " ", "that", " ", "the", " ", "Hankel", " ", "matrix", " ", "is", " ", "positive", " ", "semidefinite", " ", "and", " ", "recursively", " ", RowBox[{"generated", ".", " ", "First"}]}], ",", " ", RowBox[{ "use", " ", "nested", " ", "determinants", " ", "to", " ", "find", " ", "the", " ", "largest", " ", "positive", " ", "definite", " ", "corner"}], ",", " ", RowBox[{"with", " ", "size", " ", RowBox[{"i", "."}]}]}], " ", "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"det", " ", "=", " ", "1"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"i", " ", "=", " ", "1"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"While", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"det", " ", ">", " ", "0"}], ")"}], " ", "&&", " ", RowBox[{"(", RowBox[{"i", "\[LessEqual]", " ", RowBox[{"n", "+", "1"}]}], ")"}]}], ",", "\[IndentingNewLine]", " ", RowBox[{ RowBox[{"det", " ", "=", " ", RowBox[{"Det", "[", RowBox[{"H", "[", RowBox[{"[", RowBox[{ RowBox[{"Range", "[", RowBox[{"1", ",", "i"}], "]"}], ",", RowBox[{"Range", "[", RowBox[{"1", ",", "i"}], "]"}]}], "]"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"Print", "[", RowBox[{"N", "[", "det", "]"}], "]"}], ";", "\[IndentingNewLine]", " ", RowBox[{"i", " ", "=", " ", RowBox[{"i", "+", " ", "1"}]}], ";"}]}], "]"}], ";"}]}]}]], "Input"], Cell[CellGroupData[{ Cell[BoxData["2.000000000000002`"], "Print", CellChangeTimes->{3.422263346561125*^9, 3.422263497217375*^9, 3.422263791936125*^9, 3.42226408176425*^9, 3.42226417613925*^9, 3.422264241373625*^9, 3.42226552951425*^9, 3.4222655857955*^9, 3.422265688029875*^9, 3.422265926936125*^9, 3.42226610438925*^9, 3.422266151467375*^9, 3.422266807998625*^9, 3.422266862686125*^9, 3.42226689326425*^9, 3.422266934467375*^9, 3.422267633748625*^9, 3.422267705436125*^9, 3.422267789092375*^9, 3.422267824358*^9, 3.42226786576425*^9, 3.422354288358*^9, 3.4223578602252502`*^9, 3.4223579215846252`*^9, 3.4223581805533752`*^9, 3.4223582136783752`*^9, 3.4223583174283752`*^9, 3.4223584072721252`*^9, 3.4223584456783752`*^9, 3.4223584889752502`*^9, 3.4223587298033752`*^9, 3.4223589459440002`*^9, 3.4223589956158752`*^9, 3.422792572765625*^9, 3.42279300246875*^9, 3.422793162890625*^9, 3.422793231515625*^9, 3.422793387359375*^9, 3.42279351696875*^9, 3.422794036703125*^9, 3.422794198297901*^9, 3.422794319844776*^9, 3.422794588813526*^9}], Cell[BoxData["36.33333333333339`"], "Print", CellChangeTimes->{3.422263346561125*^9, 3.422263497217375*^9, 3.422263791936125*^9, 3.42226408176425*^9, 3.42226417613925*^9, 3.422264241373625*^9, 3.42226552951425*^9, 3.4222655857955*^9, 3.422265688029875*^9, 3.422265926936125*^9, 3.42226610438925*^9, 3.422266151467375*^9, 3.422266807998625*^9, 3.422266862686125*^9, 3.42226689326425*^9, 3.422266934467375*^9, 3.422267633748625*^9, 3.422267705436125*^9, 3.422267789092375*^9, 3.422267824358*^9, 3.42226786576425*^9, 3.422354288358*^9, 3.4223578602252502`*^9, 3.4223579215846252`*^9, 3.4223581805533752`*^9, 3.4223582136783752`*^9, 3.4223583174283752`*^9, 3.4223584072721252`*^9, 3.4223584456783752`*^9, 3.4223584889752502`*^9, 3.4223587298033752`*^9, 3.4223589459440002`*^9, 3.4223589956158752`*^9, 3.422792572765625*^9, 3.42279300246875*^9, 3.422793162890625*^9, 3.422793231515625*^9, 3.422793387359375*^9, 3.42279351696875*^9, 3.422794036703125*^9, 3.422794198297901*^9, 3.422794319844776*^9, 3.422794588844776*^9}], Cell[BoxData["218.40370370367148`"], "Print", CellChangeTimes->{3.422263346561125*^9, 3.422263497217375*^9, 3.422263791936125*^9, 3.42226408176425*^9, 3.42226417613925*^9, 3.422264241373625*^9, 3.42226552951425*^9, 3.4222655857955*^9, 3.422265688029875*^9, 3.422265926936125*^9, 3.42226610438925*^9, 3.422266151467375*^9, 3.422266807998625*^9, 3.422266862686125*^9, 3.42226689326425*^9, 3.422266934467375*^9, 3.422267633748625*^9, 3.422267705436125*^9, 3.422267789092375*^9, 3.422267824358*^9, 3.42226786576425*^9, 3.422354288358*^9, 3.4223578602252502`*^9, 3.4223579215846252`*^9, 3.4223581805533752`*^9, 3.4223582136783752`*^9, 3.4223583174283752`*^9, 3.4223584072721252`*^9, 3.4223584456783752`*^9, 3.4223584889752502`*^9, 3.4223587298033752`*^9, 3.4223589459440002`*^9, 3.4223589956158752`*^9, 3.422792572765625*^9, 3.42279300246875*^9, 3.422793162890625*^9, 3.422793231515625*^9, 3.422793387359375*^9, 3.42279351696875*^9, 3.422794036703125*^9, 3.422794198297901*^9, 3.422794319844776*^9, 3.422794588891651*^9}], Cell[BoxData["11752.840002632342`"], "Print", CellChangeTimes->{3.422263346561125*^9, 3.422263497217375*^9, 3.422263791936125*^9, 3.42226408176425*^9, 3.42226417613925*^9, 3.422264241373625*^9, 3.42226552951425*^9, 3.4222655857955*^9, 3.422265688029875*^9, 3.422265926936125*^9, 3.42226610438925*^9, 3.422266151467375*^9, 3.422266807998625*^9, 3.422266862686125*^9, 3.42226689326425*^9, 3.422266934467375*^9, 3.422267633748625*^9, 3.422267705436125*^9, 3.422267789092375*^9, 3.422267824358*^9, 3.42226786576425*^9, 3.422354288358*^9, 3.4223578602252502`*^9, 3.4223579215846252`*^9, 3.4223581805533752`*^9, 3.4223582136783752`*^9, 3.4223583174283752`*^9, 3.4223584072721252`*^9, 3.4223584456783752`*^9, 3.4223584889752502`*^9, 3.4223587298033752`*^9, 3.4223589459440002`*^9, 3.4223589956158752`*^9, 3.422792572765625*^9, 3.42279300246875*^9, 3.422793162890625*^9, 3.422793231515625*^9, 3.422793387359375*^9, 3.42279351696875*^9, 3.422794036703125*^9, 3.422794198297901*^9, 3.422794319844776*^9, 3.422794588922901*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{ "If", " ", "the", " ", "largest", " ", "positive", " ", "definite", " ", "corner", " ", "is", " ", "followed", " ", "by", " ", "a", " ", "corner", " ", "with", " ", "negative", " ", "determinant"}], ",", " ", RowBox[{ "then", " ", "H", " ", "is", " ", "not", " ", "positive", " ", "semidefinite"}], ",", " ", RowBox[{"so", " ", "there", " ", "is", " ", "no", " ", RowBox[{"measure", ":"}]}]}], " ", "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"det", "<", "0"}], ",", " ", RowBox[{ "Print", "[", " ", "\"\\"", "]"}], " ", ",", " ", "\[IndentingNewLine]", " ", "\[IndentingNewLine]", " ", RowBox[{"(*", RowBox[{ "The", " ", "rest", " ", "of", " ", "the", " ", "code", " ", "for", " ", "this", " ", "function", " ", "is", " ", "the", " ", "\"\\"", " ", "part", " ", "of", " ", "this", " ", "\"\\""}], "*)"}], " ", "\[IndentingNewLine]", "\[IndentingNewLine]", " ", RowBox[{"(*", RowBox[{ "Test", " ", "to", " ", "see", " ", "if", " ", "the", " ", "moment", " ", "data", " ", "is", " ", "recursively", " ", "generated"}], "*)"}], " ", "\[IndentingNewLine]", "\[IndentingNewLine]", " ", RowBox[{"(*", RowBox[{ RowBox[{"At", " ", "this", " ", "point"}], ",", " ", RowBox[{ "we", " ", "know", " ", "that", " ", "either", " ", "the", " ", "Hankel", " ", "matrix", " ", "is", " ", "positive", " ", "definite"}], ",", " ", RowBox[{ "or", " ", "the", " ", "largest", " ", "positive", " ", "definite", " ", "corner", " ", "is", " ", "followed", " ", "by", " ", "a", " ", "singular", " ", RowBox[{"corner", ".", " ", "If"}], " ", "the", " ", "largest", " ", "positive", " ", "definite", " ", "corner", " ", "is", " ", "followed", " ", "by", " ", "a", " ", "singular", "\[IndentingNewLine]", "corner"}], ",", " ", RowBox[{ "K", " ", "denotes", " ", "the", " ", "largest", " ", "positive", " ", "definite", " ", RowBox[{"corner", ".", " ", "Otherwise"}]}], ",", "\[IndentingNewLine]", " ", RowBox[{"K", " ", "denotes", " ", "the", " ", RowBox[{"(", RowBox[{"positive", " ", "definite"}], ")"}], " ", "Hankel", " ", "matrix", " ", "H", RowBox[{ RowBox[{"(", "n", ")"}], "."}]}]}], "*)"}], " ", "\[IndentingNewLine]", "\[IndentingNewLine]", " ", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"det", " ", "==", " ", "0"}], ",", " ", "\[IndentingNewLine]", RowBox[{ RowBox[{"K", " ", "=", " ", RowBox[{"H", "[", RowBox[{"[", RowBox[{ RowBox[{"Range", "[", RowBox[{"1", ",", RowBox[{"i", "-", "2"}]}], "]"}], ",", RowBox[{"Range", "[", RowBox[{"1", ",", RowBox[{"i", "-", "2"}]}], "]"}]}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"i", "=", RowBox[{"i", "-", "1"}]}]}], " ", RowBox[{"(*", RowBox[{ "size", " ", "of", " ", "the", " ", "largest", " ", "positive", " ", "definite", " ", "block", " ", "when", " ", "H", " ", "is", " ", "singular"}], "*)"}], ",", "\[IndentingNewLine]", RowBox[{"K", " ", "=", " ", RowBox[{"H", "[", RowBox[{"[", RowBox[{ RowBox[{"Range", "[", RowBox[{"1", ",", RowBox[{"i", "-", "1"}]}], "]"}], ",", RowBox[{"Range", "[", RowBox[{"1", ",", RowBox[{"i", "-", "1"}]}], "]"}]}], "]"}], "]"}]}]}], "]"}], " ", RowBox[{"(*", RowBox[{"K", "=", RowBox[{ "H", " ", "when", " ", "H", " ", "is", " ", "positive", " ", "definite"}]}], "*)"}], ";"}]}], "]"}], "\[IndentingNewLine]", " ", RowBox[{"MatrixForm", "[", "K", "]"}]}]}]], "Input"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"2.000000000000002`", "9.00000000000001`", "58.66666666666673`", "427.5000000000005`"}, {"9.00000000000001`", "58.66666666666673`", "427.5000000000005`", "3198.4000000000037`"}, {"58.66666666666673`", "427.5000000000005`", "3198.4000000000037`", "24093.00000000003`"}, {"427.5000000000005`", "3198.4000000000037`", "24093.00000000003`", "181962.28571428594`"} }, GridBoxAlignment->{ "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output", CellChangeTimes->{3.422263350154875*^9, 3.422263499217375*^9, 3.4222637945455*^9, 3.422264083779875*^9, 3.422264178342375*^9, 3.422264247858*^9, 3.42226553151425*^9, 3.422265587733*^9, 3.4222656901705*^9, 3.422265927811125*^9, 3.422266104967375*^9, 3.422266151748625*^9, 3.4222668090455*^9, 3.422266862998625*^9, 3.42226689357675*^9, 3.422266935154875*^9, 3.42226763407675*^9, 3.422267706592375*^9, 3.422267789467375*^9, 3.422267824498625*^9, 3.422267865842375*^9, 3.422354289404875*^9, 3.4223578613502502`*^9, 3.4223579225065002`*^9, 3.4223581813346252`*^9, 3.4223582141627502`*^9, 3.4223583179752502`*^9, 3.4223584077721252`*^9, 3.4223584466002502`*^9, 3.4223584904283752`*^9, 3.4223587308502502`*^9, 3.4223589468033752`*^9, 3.4223589964908752`*^9, 3.422792573921875*^9, 3.422793003921875*^9, 3.42279316384375*^9, 3.422793232359375*^9, 3.4227933881875*^9, 3.42279351853125*^9, 3.42279403775*^9, 3.422794198891651*^9, 3.422794320516651*^9, 3.422794589516651*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ "Htest", " ", "is", " ", "the", " ", "moment", " ", "data", " ", "to", " ", "be", " ", "tested", " ", "for", " ", RowBox[{"recursiveness", ":"}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"2", "n"}], "\[Equal]", "m"}], ",", " ", RowBox[{"Htest", " ", "=", " ", RowBox[{"H", "[", RowBox[{"[", RowBox[{ RowBox[{"Range", "[", RowBox[{"1", ",", RowBox[{"n", "+", "1"}]}], "]"}], ",", RowBox[{"Range", "[", RowBox[{"1", ",", RowBox[{"n", "+", "1"}]}], "]"}]}], "]"}], "]"}]}], ",", "\[IndentingNewLine]", RowBox[{"Htest", " ", "=", " ", RowBox[{"H", "[", RowBox[{"[", RowBox[{ RowBox[{"Range", "[", RowBox[{"1", ",", RowBox[{"n", "+", "1"}]}], "]"}], ",", RowBox[{"Range", "[", RowBox[{"1", ",", RowBox[{"n", "+", "2"}]}], "]"}]}], "]"}], "]"}]}]}], "]"}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"MatrixForm", "[", "Htest", "]"}]}]}]], "Input"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"2.000000000000002`", "9.00000000000001`", "58.66666666666673`", "427.5000000000005`", "3198.4000000000037`"}, {"9.00000000000001`", "58.66666666666673`", "427.5000000000005`", "3198.4000000000037`", "24093.00000000003`"}, {"58.66666666666673`", "427.5000000000005`", "3198.4000000000037`", "24093.00000000003`", "181962.28571428594`"}, {"427.5000000000005`", "3198.4000000000037`", "24093.00000000003`", "181962.28571428594`", "1.3765837500000016`*^6"} }, GridBoxAlignment->{ "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output", CellChangeTimes->{3.42226335332675*^9, 3.4222635011705*^9, 3.42226379682675*^9, 3.422264085998625*^9, 3.422264180654875*^9, 3.422264250436125*^9, 3.422265533561125*^9, 3.422265589873625*^9, 3.42226569201425*^9, 3.422265928623625*^9, 3.422266105358*^9, 3.422266152061125*^9, 3.422266809873625*^9, 3.422266863279875*^9, 3.422266893904875*^9, 3.422266935873625*^9, 3.42226763438925*^9, 3.42226770801425*^9, 3.422267789842375*^9, 3.422267824717375*^9, 3.422267866311125*^9, 3.422354290483*^9, 3.4223578626471252`*^9, 3.4223579240533752`*^9, 3.4223581820065002`*^9, 3.4223582146940002`*^9, 3.4223583185221252`*^9, 3.4223584082565002`*^9, 3.4223584473346252`*^9, 3.4223584915221252`*^9, 3.4223587318971252`*^9, 3.4223589475533752`*^9, 3.4223589971471252`*^9, 3.4227925750625*^9, 3.422793006125*^9, 3.42279316475*^9, 3.422793233203125*^9, 3.42279338890625*^9, 3.422793520375*^9, 3.422794038765625*^9, 3.422794199376026*^9, 3.422794321172901*^9, 3.422794590079151*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{"Recursiveness", " ", RowBox[{"test", ":"}]}], " ", "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", " ", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"MatrixRank", "[", "Htest", "]"}], " ", ">", " ", RowBox[{"MatrixRank", "[", "K", "]"}]}], ",", "\[IndentingNewLine]", " ", RowBox[{ "Print", "[", " ", "\"\\"", "]"}], " ", ",", "\[IndentingNewLine]", " ", RowBox[{"Print", "[", "\"\\"", "]"}]}], "]"}], " "}]], "Input"], Cell[BoxData["\<\"there is a measure\"\>"], "Print", CellChangeTimes->{3.422263355998625*^9, 3.4222635031705*^9, 3.422263799029875*^9, 3.422264088467375*^9, 3.42226418251425*^9, 3.42226425288925*^9, 3.422265535529875*^9, 3.42226559170175*^9, 3.422265693842375*^9, 3.422265929311125*^9, 3.422266105733*^9, 3.4222661526705*^9, 3.4222668107955*^9, 3.422266863561125*^9, 3.422266894217375*^9, 3.422266936561125*^9, 3.42226763470175*^9, 3.42226770857675*^9, 3.422267790217375*^9, 3.422267824811125*^9, 3.422267866811125*^9, 3.42235429188925*^9, 3.4223578638971252`*^9, 3.4223579249752502`*^9, 3.4223581826158752`*^9, 3.4223582152565002`*^9, 3.4223583190846252`*^9, 3.4223584086940002`*^9, 3.4223584480533752`*^9, 3.4223584927721252`*^9, 3.4223587326940002`*^9, 3.4223589481940002`*^9, 3.4223589978033752`*^9, 3.42279257625*^9, 3.42279300953125*^9, 3.42279316559375*^9, 3.42279323396875*^9, 3.422793389546875*^9, 3.422793521640625*^9, 3.422794039875*^9, 3.422794199844776*^9, 3.422794321860401*^9, 3.422794590688526*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{"data", " ", "is", " ", "recursively", " ", "generated"}], ";", " ", RowBox[{"a", " ", "measure", " ", "exists"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Print", "[", "\"\\"", "]"}], ";"}], "\[IndentingNewLine]", " ", RowBox[{ RowBox[{"i", " ", "=", " ", RowBox[{"i", "-", "1"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"express", " ", "column", " ", RowBox[{"t", "^", "i"}], " ", RowBox[{"(", RowBox[{"the", " ", "first", " ", "column", " ", "beyond", " ", "K"}], ")"}], "as", " ", "a", " ", "linear", " ", "combination", " ", "of", " ", "columns", " ", "1"}], ",", "1", ",", "...", ",", RowBox[{"t", "^", RowBox[{"(", RowBox[{"i", "-", "1"}], ")"}]}]}], "*)"}], "\[IndentingNewLine]", " ", RowBox[{ RowBox[{"J", " ", "=", " ", RowBox[{"N", "[", RowBox[{ RowBox[{"Inverse", "[", "K", "]"}], ",", "100"}], "]"}]}], ";"}], "\[IndentingNewLine]", " ", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", "0"}], ",", " ", RowBox[{"j", "\[LessEqual]", RowBox[{"i", "-", "1"}]}], ",", RowBox[{"j", "=", RowBox[{"j", "+", "1"}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"v", "[", "j", "]"}], " ", "=", " ", RowBox[{"bb", "[", RowBox[{"i", "+", "j"}], "]"}]}]}], "]"}], ";"}], " ", "\[IndentingNewLine]", RowBox[{ RowBox[{"w", " ", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"v", "[", "j", "]"}], ",", RowBox[{"{", RowBox[{"j", ",", "0", ",", RowBox[{"i", "-", "1"}]}], "}"}]}], "]"}]}], ";"}], " ", RowBox[{"(*", RowBox[{"colunmn", " ", RowBox[{"t", "^", "i"}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"z", " ", "=", " ", RowBox[{"J", ".", "w"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"N", "[", "z", "]"}], "]"}], ";"}], "\[IndentingNewLine]"}]}]], "Input"], Cell[CellGroupData[{ Cell[BoxData["\<\"A measure exists\"\>"], "Print", CellChangeTimes->{3.42226335845175*^9, 3.422263505404875*^9, 3.422263800686125*^9, 3.422264090373625*^9, 3.422264185108*^9, 3.42226425738925*^9, 3.4222655375455*^9, 3.422265593733*^9, 3.422265695748625*^9, 3.422265929873625*^9, 3.422266106217375*^9, 3.42226615301425*^9, 3.42226681138925*^9, 3.422266863858*^9, 3.422266894529875*^9, 3.4222669371705*^9, 3.422267635029875*^9, 3.422267709154875*^9, 3.422267790561125*^9, 3.42226782532675*^9, 3.422267866858*^9, 3.422354292967375*^9, 3.4223578650221252`*^9, 3.4223579257408752`*^9, 3.4223581832252502`*^9, 3.4223582157877502`*^9, 3.4223583196158752`*^9, 3.4223584091002502`*^9, 3.4223584486158752`*^9, 3.4223584940377502`*^9, 3.4223587336940002`*^9, 3.4223589488502502`*^9, 3.4223589984440002`*^9, 3.422792577453125*^9, 3.4227930111875*^9, 3.422793166390625*^9, 3.42279323478125*^9, 3.422793390171875*^9, 3.422793522734375*^9, 3.422794040875*^9, 3.422794200360401*^9, 3.422794322438526*^9, 3.422794591297901*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"-", "119.4587008185456`"}], ",", "199.2025349883603`", ",", RowBox[{"-", "103.13361499869751`"}], ",", "18.000000000004952`"}], "}"}]], "Print", CellChangeTimes->{3.42226335845175*^9, 3.422263505404875*^9, 3.422263800686125*^9, 3.422264090373625*^9, 3.422264185108*^9, 3.42226425738925*^9, 3.4222655375455*^9, 3.422265593733*^9, 3.422265695748625*^9, 3.422265929873625*^9, 3.422266106217375*^9, 3.42226615301425*^9, 3.42226681138925*^9, 3.422266863858*^9, 3.422266894529875*^9, 3.4222669371705*^9, 3.422267635029875*^9, 3.422267709154875*^9, 3.422267790561125*^9, 3.42226782532675*^9, 3.422267866858*^9, 3.422354292967375*^9, 3.4223578650221252`*^9, 3.4223579257408752`*^9, 3.4223581832252502`*^9, 3.4223582157877502`*^9, 3.4223583196158752`*^9, 3.4223584091002502`*^9, 3.4223584486158752`*^9, 3.4223584940377502`*^9, 3.4223587336940002`*^9, 3.4223589488502502`*^9, 3.4223589984440002`*^9, 3.422792577453125*^9, 3.4227930111875*^9, 3.422793166390625*^9, 3.42279323478125*^9, 3.422793390171875*^9, 3.422793522734375*^9, 3.422794040875*^9, 3.422794200360401*^9, 3.422794322438526*^9, 3.422794591329151*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{"build", " ", "the", " ", "characteristic", " ", "polynomial"}], ",", " ", RowBox[{ "whose", " ", "roots", " ", "form", " ", "the", " ", "support", " ", "of", " ", "the", " ", "representing", " ", "measure"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"p", " ", "=", " ", "0"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", "0"}], ",", RowBox[{"j", "\[LessEqual]", RowBox[{"i", "-", "1"}]}], ",", RowBox[{"j", "=", RowBox[{"j", "+", "1"}]}], ",", "\[IndentingNewLine]", RowBox[{"p", " ", "=", " ", RowBox[{"p", " ", "+", " ", RowBox[{ RowBox[{"z", "[", RowBox[{"[", RowBox[{"j", "+", "1"}], "]"}], "]"}], " ", RowBox[{"t", "^", "j"}]}]}]}]}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"q", " ", "=", " ", RowBox[{ RowBox[{"t", "^", RowBox[{"(", "i", ")"}]}], "-", "p"}]}], ";"}], " ", RowBox[{"(*", RowBox[{"characteristic", " ", "polynomial"}], "*)"}], "\[IndentingNewLine]", RowBox[{"Print", "[", RowBox[{ RowBox[{"N", "[", "q", "]"}], ",", " ", "\"\< char. poly.\>\""}], "]"}], "\[IndentingNewLine]"}]}]], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{"119.4587008185456`", "\[InvisibleSpace]", "-", RowBox[{"199.2025349883603`", " ", "t"}], "+", RowBox[{"103.13361499869751`", " ", SuperscriptBox["t", "2"]}], "-", RowBox[{"18.000000000004952`", " ", SuperscriptBox["t", "3"]}], "+", RowBox[{ SuperscriptBox["t", "4"], "\[InvisibleSpace]", "\<\" char. poly.\"\>"}]}], SequenceForm[ 119.4587008185456 - 199.2025349883603 $CellContext`t + 103.13361499869751` $CellContext`t^2 - 18.000000000004952` $CellContext`t^3 + $CellContext`t^4, " char. poly."], Editable->False]], "Print", CellChangeTimes->{ 3.422263361483*^9, 3.422263507373625*^9, {3.422263802858*^9, 3.422263808748625*^9}, 3.422264092811125*^9, 3.42226418738925*^9, 3.42226425970175*^9, 3.422265539686125*^9, 3.42226559576425*^9, 3.42226569763925*^9, 3.422265930717375*^9, 3.422266106983*^9, 3.42226615507675*^9, 3.422266811904875*^9, 3.422266864186125*^9, 3.4222668949205*^9, 3.422266937873625*^9, 3.422267635358*^9, 3.422267709936125*^9, 3.422267790936125*^9, 3.42226782545175*^9, 3.422267867108*^9, 3.42235429420175*^9, 3.4223578662565002`*^9, 3.4223579265377502`*^9, 3.4223581838658752`*^9, 3.4223582163190002`*^9, 3.4223583201783752`*^9, 3.4223584095065002`*^9, 3.4223584491158752`*^9, 3.4223584961471252`*^9, 3.4223587350221252`*^9, 3.4223589495221252`*^9, 3.4223589993346252`*^9, 3.422792578640625*^9, 3.42279301246875*^9, 3.42279316725*^9, 3.422793235546875*^9, 3.422793390796875*^9, 3.422793523921875*^9, 3.4227940424375*^9, 3.422794200860401*^9, 3.422794323063526*^9, 3.422794591969776*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"roots", " ", "=", " ", RowBox[{"Solve", "[", RowBox[{ RowBox[{"q", "\[Equal]", "0"}], ",", "t"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Print", "[", RowBox[{ RowBox[{"N", "[", "roots", "]"}], ",", " ", "\"\< roots of char. poly.\>\""}], "]"}]}], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"t", "\[Rule]", "1.1955000529356674`"}], "}"}], ",", RowBox[{"{", RowBox[{"t", "\[Rule]", "1.771142161794588`"}], "}"}], ",", RowBox[{"{", RowBox[{"t", "\[Rule]", "7.228857838084334`"}], "}"}], ",", RowBox[{"{", RowBox[{"t", "\[Rule]", "7.804499947190363`"}], "}"}]}], "}"}], "\[InvisibleSpace]", "\<\" roots of char. poly.\"\>"}], SequenceForm[{{$CellContext`t -> 1.1955000529356674`}, {$CellContext`t -> 1.771142161794588}, {$CellContext`t -> 7.228857838084334}, {$CellContext`t -> 7.804499947190363}}, " roots of char. poly."], Editable->False]], "Print", CellChangeTimes->{3.422263364248625*^9, 3.42226350951425*^9, 3.42226381095175*^9, 3.422264094904875*^9, 3.42226418920175*^9, 3.422264262436125*^9, 3.422265541748625*^9, 3.42226559895175*^9, 3.422265699842375*^9, 3.4222659315455*^9, 3.422266107592375*^9, 3.42226615538925*^9, 3.422266812529875*^9, 3.42226686457675*^9, 3.42226689526425*^9, 3.422266939061125*^9, 3.422267635654875*^9, 3.422267710467375*^9, 3.422267791358*^9, 3.422267825623625*^9, 3.422267867217375*^9, 3.42235429570175*^9, 3.4223578673971252`*^9, 3.4223579273658752`*^9, 3.4223581844283752`*^9, 3.4223582168502502`*^9, 3.4223583207408752`*^9, 3.4223584099127502`*^9, 3.4223584496002502`*^9, 3.4223584976158752`*^9, 3.4223587368815002`*^9, 3.4223589505221252`*^9, 3.4223590001940002`*^9, 3.422792579671875*^9, 3.422793014546875*^9, 3.422793168140625*^9, 3.422793236265625*^9, 3.422793391421875*^9, 3.4227935251875*^9, 3.4227940435*^9, 3.422794201360401*^9, 3.422794323735401*^9, 3.422794592626026*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", "1"}], ",", RowBox[{"j", "\[LessEqual]", "i"}], ",", RowBox[{"j", "=", RowBox[{"j", "+", "1"}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"x", "[", RowBox[{"j", "-", "1"}], "]"}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"roots", "[", RowBox[{"[", "j", "]"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], "[", RowBox[{"[", "2", "]"}], "]"}]}]}], "]"}], ";"}], " ", RowBox[{"(*", RowBox[{ "nodes", " ", "of", " ", "the", " ", "representing", " ", "measure"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", " ", RowBox[{ RowBox[{"ii", "=", "1"}], ",", RowBox[{"ii", "\[LessEqual]", RowBox[{"i", "-", "1"}]}], ",", RowBox[{"ii", "=", RowBox[{"ii", "+", "1"}]}], ",", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", "0"}], ",", RowBox[{"j", "\[LessEqual]", RowBox[{"i", "-", "1"}]}], ",", RowBox[{"j", "=", RowBox[{"j", "+", "1"}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"TTab", "[", RowBox[{"ii", ",", "j"}], "]"}], "=", " ", RowBox[{ RowBox[{"x", "[", "j", "]"}], "^", "ii"}]}]}], "]"}]}], "]"}], ";"}], "\n", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", "0"}], ",", RowBox[{"j", "\[LessEqual]", RowBox[{"i", "-", "1"}]}], ",", RowBox[{"j", "=", RowBox[{"j", "+", "1"}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"TTab", "[", RowBox[{"0", ",", "j"}], "]"}], "=", "1"}]}], "]"}], ";"}], "\n", RowBox[{ RowBox[{"V", "=", " ", RowBox[{"Table", "[", RowBox[{ RowBox[{"TTab", "[", RowBox[{"rr", ",", "ss"}], "]"}], ",", RowBox[{"{", RowBox[{"rr", ",", "0", ",", RowBox[{"i", "-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"ss", ",", "0", ",", RowBox[{"i", "-", "1"}]}], "}"}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"MatrixForm", "[", RowBox[{"N", "[", "V", "]"}], "]"}], RowBox[{"(*", "Vandermonde", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", " ", "\"\\"", "]"}], ";"}]}], "Input"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"1.`", "1.`", "1.`", "1.`"}, {"1.1955000529356674`", "1.771142161794588`", "7.228857838084334`", "7.804499947190363`"}, {"1.4292203765691835`", "3.1369445572864065`", "52.25638564323331`", "60.910219425694386`"}, {"1.7086330358451935`", "5.555974764622012`", "377.7539829470448`", "475.37380429118525`"} }, GridBoxAlignment->{ "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output", CellChangeTimes->{3.422263366279875*^9, 3.42226351157675*^9, 3.422263812608*^9, 3.422264096404875*^9, 3.422264191061125*^9, 3.422264264217375*^9, 3.422265543279875*^9, 3.42226560045175*^9, 3.4222657015455*^9, 3.422265932233*^9, 3.42226610813925*^9, 3.422266155858*^9, 3.422266812904875*^9, 3.4222668649205*^9, 3.42226689563925*^9, 3.42226693963925*^9, 3.422267635998625*^9, 3.422267711217375*^9, 3.422267791748625*^9, 3.422267825686125*^9, 3.422267867717375*^9, 3.422354296592375*^9, 3.4223578684440002`*^9, 3.4223579280846252`*^9, 3.4223581851158752`*^9, 3.4223582173658752`*^9, 3.4223583213815002`*^9, 3.4223584103033752`*^9, 3.4223584500846252`*^9, 3.4223584998658752`*^9, 3.4223587378971252`*^9, 3.4223589511627502`*^9, 3.422792580984375*^9, 3.42279301596875*^9, 3.42279316909375*^9, 3.422793237015625*^9, 3.422793392046875*^9, 3.42279352625*^9, 3.4227940451875*^9, 3.422794201829151*^9, 3.422794324391651*^9, 3.422794593485401*^9}], Cell[BoxData["\<\"Vandermonde\"\>"], "Print", CellChangeTimes->{3.4222633662955*^9, 3.42226351157675*^9, 3.422263812608*^9, 3.422264096404875*^9, 3.422264191061125*^9, 3.422264264217375*^9, 3.422265543279875*^9, 3.42226560045175*^9, 3.4222657015455*^9, 3.422265932248625*^9, 3.422266108154875*^9, 3.422266155873625*^9, 3.4222668129205*^9, 3.4222668649205*^9, 3.42226689563925*^9, 3.422266939654875*^9, 3.422267635998625*^9, 3.422267711217375*^9, 3.422267791748625*^9, 3.422267825686125*^9, 3.422267867733*^9, 3.422354296608*^9, 3.4223578685065002`*^9, 3.4223579281471252`*^9, 3.4223581851627502`*^9, 3.4223582174283752`*^9, 3.4223583214440002`*^9, 3.4223584103502502`*^9, 3.4223584501627502`*^9, 3.4223584999440002`*^9, 3.4223587379596252`*^9, 3.4223589512252502`*^9, 3.422792581015625*^9, 3.422793016*^9, 3.422793169125*^9, 3.422793237046875*^9, 3.42279339209375*^9, 3.42279352628125*^9, 3.42279404521875*^9, 3.422794201876026*^9, 3.422794324422901*^9, 3.422794593516651*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"ii", "=", "0"}], ",", " ", RowBox[{"ii", "\[LessEqual]", RowBox[{"i", "-", "1"}]}], ",", RowBox[{"ii", "=", RowBox[{"ii", "+", "1"}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"vvv", "[", "ii", "]"}], " ", "=", " ", RowBox[{"bb", "[", "ii", "]"}]}]}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"w", " ", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"vvv", "[", "ii", "]"}], ",", RowBox[{"{", RowBox[{"ii", ",", "0", ",", RowBox[{"i", "-", "1"}]}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"weights", "=", RowBox[{"N", "[", RowBox[{ RowBox[{ RowBox[{"Inverse", "[", "V", "]"}], ".", "w"}], ",", "100"}], "]"}]}], ";"}], RowBox[{"(*", RowBox[{"densities", " ", "of", " ", "the", " ", "measure"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"PP", "[", "1", "]"}], ";"}], RowBox[{"(*", RowBox[{"print", " ", "the", " ", "nodes", " ", "and", " ", "densities"}], "*)"}]}]}], "Input"], Cell[CellGroupData[{ Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"nodes\"\>", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"weights\"\>"}], SequenceForm["nodes", " ", "weights"], Editable->False]], "Print", CellChangeTimes->{3.422792592625*^9, 3.422793017296875*^9, 3.42279316996875*^9, 3.422793237734375*^9, 3.4227933926875*^9, 3.4227935279375*^9, 3.42279404684375*^9, 3.422794202344776*^9, 3.422794325110401*^9, 3.422794594313526*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "1.1955000529356674`", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "0.47124745417920205`"}], SequenceForm[1.1955000529356674`, " ", 0.47124745417920205`], Editable->False]], "Print", CellChangeTimes->{3.422792592625*^9, 3.422793017296875*^9, 3.42279316996875*^9, 3.422793237734375*^9, 3.4227933926875*^9, 3.4227935279375*^9, 3.42279404684375*^9, 3.422794202344776*^9, 3.422794325110401*^9, 3.422794594344776*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "1.771142161794588`", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "0.5287525458148301`"}], SequenceForm[1.771142161794588, " ", 0.5287525458148301], Editable->False]], "Print", CellChangeTimes->{3.422792592625*^9, 3.422793017296875*^9, 3.42279316996875*^9, 3.422793237734375*^9, 3.4227933926875*^9, 3.4227935279375*^9, 3.42279404684375*^9, 3.422794202344776*^9, 3.422794325110401*^9, 3.422794594391651*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "7.228857838084334`", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "0.5287525458745922`"}], SequenceForm[7.228857838084334, " ", 0.5287525458745922], Editable->False]], "Print", CellChangeTimes->{3.422792592625*^9, 3.422793017296875*^9, 3.42279316996875*^9, 3.422793237734375*^9, 3.4227933926875*^9, 3.4227935279375*^9, 3.42279404684375*^9, 3.422794202344776*^9, 3.422794325110401*^9, 3.422794594454151*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "7.804499947190363`", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "0.4712474541313019`"}], SequenceForm[7.804499947190363, " ", 0.4712474541313019], Editable->False]], "Print", CellChangeTimes->{3.422792592625*^9, 3.422793017296875*^9, 3.42279316996875*^9, 3.422793237734375*^9, 3.4227933926875*^9, 3.4227935279375*^9, 3.42279404684375*^9, 3.422794202344776*^9, 3.422794325110401*^9, 3.422794594501026*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"N", "[", RowBox[{ RowBox[{ RowBox[{"weights", "[", RowBox[{"[", "1", "]"}], "]"}], " ", RowBox[{ RowBox[{"x", "[", "0", "]"}], "^", "5"}]}], "+", " ", RowBox[{ RowBox[{"weights", "[", RowBox[{"[", "2", "]"}], "]"}], " ", RowBox[{ RowBox[{"x", "[", "1", "]"}], "^", "5"}]}], "+", " ", RowBox[{ RowBox[{"weights", "[", RowBox[{"[", "3", "]"}], "]"}], " ", RowBox[{ RowBox[{"x", "[", "2", "]"}], "^", "5"}]}]}], "]"}]], "Input", CellChangeTimes->{{3.421152432515625*^9, 3.42115247015625*^9}, { 3.42115253740625*^9, 3.42115255953125*^9}, {3.421152596828125*^9, 3.421152608203125*^9}, {3.42115265034375*^9, 3.421152706921875*^9}, { 3.421152759375*^9, 3.4211528155625*^9}, {3.4211528465625*^9, 3.421152858453125*^9}}], Cell[BoxData["10447.972130169212`"], "Output", CellChangeTimes->{ 3.422263376779875*^9, 3.422263515983*^9, 3.42226381713925*^9, 3.422264102733*^9, 3.42226419557675*^9, 3.42226426895175*^9, 3.42226554776425*^9, 3.422265612233*^9, 3.4222657064205*^9, 3.422265933998625*^9, 3.42226610963925*^9, 3.42226615695175*^9, 3.42226681451425*^9, 3.4222668662955*^9, 3.42226689632675*^9, 3.422267636623625*^9, 3.422267712529875*^9, 3.42226779238925*^9, 3.422267826233*^9, {3.42235429907675*^9, 3.422354304717375*^9}, 3.4223578725065002`*^9, 3.4223579308971252`*^9, 3.4223581870533752`*^9, 3.4223582188346252`*^9, 3.4223583226158752`*^9, 3.4223584115690002`*^9, 3.4223584517565002`*^9, 3.4223585031158752`*^9, 3.4223587404752502`*^9, 3.4223589530221252`*^9, 3.422792598703125*^9, 3.422793030484375*^9, 3.42279317203125*^9, 3.422793238703125*^9, 3.422793393484375*^9, 3.422793608359375*^9, 3.4227940654375*^9, 3.422794205735401*^9, 3.422794715876026*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"bb", "[", "5", "]"}]], "Input", CellChangeTimes->{{3.42115261309375*^9, 3.421152620578125*^9}, { 3.4211527156875*^9, 3.421152716078125*^9}, {3.421152831359375*^9, 3.421152831609375*^9}, {3.42115287409375*^9, 3.4211528743125*^9}}], Cell[BoxData["24093.00000000003`"], "Output", CellChangeTimes->{ 3.422263379311125*^9, 3.422263518842375*^9, 3.422263822436125*^9, 3.422264104483*^9, 3.4222641977955*^9, 3.422264270248625*^9, 3.422265549373625*^9, 3.422265614123625*^9, 3.42226570857675*^9, 3.422265935279875*^9, 3.422266116733*^9, 3.42226615745175*^9, 3.42226681513925*^9, {3.422266866654875*^9, 3.422266896654875*^9}, 3.422267636936125*^9, 3.4222677135455*^9, 3.422267793092375*^9, 3.4222678264205*^9, {3.42235430045175*^9, 3.42235430626425*^9}, 3.4223578744908752`*^9, 3.4223579318346252`*^9, 3.4223581872721252`*^9, 3.4223582196627502`*^9, 3.4223583232096252`*^9, 3.4223584120377502`*^9, 3.4223584525533752`*^9, 3.4223585042408752`*^9, 3.4223587417252502`*^9, 3.4223589536940002`*^9, 3.422792601359375*^9, 3.422793031734375*^9, 3.42279317284375*^9, 3.42279323978125*^9, 3.4227933941875*^9, 3.422793609328125*^9, 3.42279406596875*^9, 3.422794206563526*^9, 3.422794717079151*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"N", "[", "%", "]"}]], "Input", CellChangeTimes->{{3.42115262378125*^9, 3.421152625015625*^9}}], Cell[BoxData["24093.00000000003`"], "Output", CellChangeTimes->{ 3.422263381342375*^9, 3.422263520342375*^9, 3.422263824248625*^9, 3.422264106123625*^9, 3.42226427120175*^9, 3.42226555157675*^9, 3.422265615733*^9, 3.42226571038925*^9, 3.42226593582675*^9, 3.422266117436125*^9, 3.4222661577955*^9, 3.422266815733*^9, { 3.422266867311125*^9, 3.422266896998625*^9}, 3.42226763770175*^9, 3.422267714233*^9, 3.4222677934205*^9, 3.42226782651425*^9, 3.422354307029875*^9, 3.4223578756627502`*^9, 3.4223579326315002`*^9, 3.4223581886158752`*^9, 3.4223582203346252`*^9, 3.4223583238033752`*^9, 3.4223584126783752`*^9, 3.4223584531315002`*^9, 3.4223585051940002`*^9, 3.4223587431940002`*^9, 3.4223589543502502`*^9, 3.4227926026875*^9, 3.422793033359375*^9, 3.422793173828125*^9, 3.4227932413125*^9, 3.422793394890625*^9, 3.42279361075*^9, 3.422794066578125*^9, 3.422794207094776*^9, 3.422794717797901*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"MatrixForm", "[", "H", "]"}]], "Input", CellChangeTimes->{{3.421151299953125*^9, 3.421151307609375*^9}}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"2.000000000000002`", "9.00000000000001`", "58.66666666666673`", "427.5000000000005`", "3198.4000000000037`"}, {"9.00000000000001`", "58.66666666666673`", "427.5000000000005`", "3198.4000000000037`", "24093.00000000003`"}, {"58.66666666666673`", "427.5000000000005`", "3198.4000000000037`", "24093.00000000003`", "181962.28571428594`"}, {"427.5000000000005`", "3198.4000000000037`", "24093.00000000003`", "181962.28571428594`", "1.3765837500000016`*^6"} }, GridBoxAlignment->{ "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output", CellChangeTimes->{ 3.422263384279875*^9, 3.422263523108*^9, 3.422263833936125*^9, 3.42226410913925*^9, 3.422264201936125*^9, 3.422264272311125*^9, 3.422265553373625*^9, 3.422265619404875*^9, 3.42226571188925*^9, 3.4222659364205*^9, 3.422266118092375*^9, 3.422266158358*^9, 3.422266816404875*^9, {3.422266867654875*^9, 3.422266897358*^9}, 3.42226763826425*^9, 3.422267714936125*^9, 3.422267793983*^9, 3.422267826717375*^9, 3.422354307983*^9, 3.4223578771002502`*^9, 3.4223579336002502`*^9, 3.4223581892877502`*^9, 3.4223582209440002`*^9, 3.4223583243502502`*^9, 3.4223584132408752`*^9, 3.4223584537565002`*^9, 3.4223585060533752`*^9, 3.4223587442721252`*^9, 3.4223589549908752`*^9, 3.4227926044375*^9, 3.42279303434375*^9, 3.422793174515625*^9, 3.42279324225*^9, 3.42279339559375*^9, 3.422793611375*^9, 3.42279406721875*^9, 3.422794207547901*^9, 3.422794718407276*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"___", "___", "___", "___", "___", "___", "___", "___", "___", "___", "___", "___", "___", "_"}]], "Input", CellChangeTimes->{{3.422263690498625*^9, 3.42226369257675*^9}}], Cell[BoxData[ RowBox[{"_", " ", SuperscriptBox["___", "13"]}]], "Output", CellChangeTimes->{ 3.422265937467375*^9, 3.4222661194205*^9, 3.422266159108*^9, 3.42226681913925*^9, {3.42226686826425*^9, 3.422266897779875*^9}, 3.422267638936125*^9, 3.422267715623625*^9, 3.422267794717375*^9, 3.422267826733*^9, 3.422354309123625*^9, 3.4223578780690002`*^9, 3.4223579346002502`*^9, 3.4223581901315002`*^9, 3.4223582214596252`*^9, 3.4223583248658752`*^9, 3.4223584137877502`*^9, 3.4223584542721252`*^9, 3.4223585070846252`*^9, 3.4223587471940002`*^9, 3.4223589556940002`*^9, 3.42279260565625*^9, 3.4227930353125*^9, 3.42279317540625*^9, 3.422793243125*^9, 3.42279339628125*^9, 3.422793612578125*^9, 3.42279406825*^9, 3.422794208188526*^9, 3.422794719266651*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"G", "[", "1", "]"}]], "Input"], Cell[BoxData["1.8560557022658646`"], "Output", CellChangeTimes->{ 3.422263651592375*^9, 3.42226386782675*^9, 3.42226411520175*^9, 3.422264209858*^9, 3.422264278779875*^9, 3.4222655606705*^9, 3.4222656224205*^9, 3.4222657144205*^9, 3.422265938342375*^9, 3.42226612013925*^9, 3.422266159936125*^9, 3.422266820186125*^9, { 3.422266868904875*^9, 3.422266898483*^9}, 3.422267639498625*^9, 3.422267716342375*^9, 3.42226779526425*^9, 3.422267826858*^9, 3.422354310279875*^9, 3.4223578789596252`*^9, 3.4223579356627502`*^9, 3.4223581908190002`*^9, 3.4223582220690002`*^9, 3.4223583253971252`*^9, 3.4223584143346252`*^9, 3.4223584547877502`*^9, 3.4223585078190002`*^9, 3.4223587480690002`*^9, 3.4223589563033752`*^9, 3.422792608609375*^9, 3.422793041078125*^9, 3.42279307865625*^9, 3.422793176109375*^9, 3.42279324396875*^9, 3.422793396984375*^9, 3.4227936134375*^9, 3.422794068765625*^9, 3.422794208688526*^9, 3.422794719844776*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{ "Gaussian", " ", "quadrature", " ", "estimate", " ", "of", " ", "degree", " ", "m"}], " ", "=", " ", RowBox[{ RowBox[{"2", "n"}], "+", RowBox[{ "1", " ", "for", " ", "integral", " ", "of", " ", "\[IndentingNewLine]", "predefined", " ", "function", " ", RowBox[{"f", "[", "x", "]"}], RowBox[{"over", " ", "[", RowBox[{"a", ",", "b"}], "]"}]}]}]}], ";", " ", RowBox[{ RowBox[{ "uses", " ", "results", " ", "of", " ", "\[IndentingNewLine]", "previous", " ", "call", " ", "to", " ", RowBox[{"THMP", "[", "m", "]"}], " ", "with", " ", "moment", " ", "data", " ", "up", " ", "to", " ", "degree", "\[IndentingNewLine]", "2", "n"}], "+", RowBox[{ "1", " ", "corresponding", " ", "to", " ", "Lebesgue", " ", "measure", " ", RowBox[{"on", " ", "[", RowBox[{"a", ",", "b"}], "]"}]}]}]}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"G", "[", "j_", "]"}], ":=", RowBox[{"(", " ", RowBox[{"(*", RowBox[{"j", " ", "is", " ", "ignored"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"ss", "=", "0"}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "0"}], ",", " ", RowBox[{"i", "\[LessEqual]", "n"}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{"ss", " ", "=", " ", RowBox[{"ss", " ", "+", " ", RowBox[{ RowBox[{"weights", "[", RowBox[{"[", RowBox[{"i", "+", "1"}], "]"}], "]"}], " ", RowBox[{"g", "[", RowBox[{"x", "[", "i", "]"}], "]"}]}]}]}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"ss", "=", " ", RowBox[{"N", "[", "ss", "]"}]}]}], ")"}]}], ";"}]}]], "Input", CellChangeTimes->{{3.4211530040625*^9, 3.42115300446875*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ "Print", " ", "the", " ", "nodes", " ", "and", " ", "weights", " ", "of", " ", "the", " ", "measure", " ", "determined", " ", "by", " ", "THMP"}], "*)"}], "\n", RowBox[{ RowBox[{"PP", "[", "mm_", "]"}], ":=", RowBox[{"(", " ", RowBox[{"(*", " ", RowBox[{"mm", " ", "is", " ", "ignored"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{ "\"\\"", ",", " ", "\"\< \>\"", ",", " ", "\"\\""}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", "0"}], ",", RowBox[{"j", "<", "i"}], ",", RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", RowBox[{"Print", "[", RowBox[{ RowBox[{"N", "[", RowBox[{"x", "[", "j", "]"}], "]"}], ",", " ", "\"\< \>\"", ",", " ", RowBox[{"N", "[", RowBox[{"weights", "[", RowBox[{"[", RowBox[{"j", "+", "1"}], "]"}], "]"}], "]"}]}], "]"}]}], "]"}], ";"}], ")"}]}]}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"G", "[", "1", "]"}]], "Input"], Cell[BoxData["1.8560557022658646`"], "Output", GeneratedCell->False, CellAutoOverwrite->False, CellChangeTimes->{{3.422793400203125*^9, 3.4227934126875*^9}, 3.422793620125*^9, 3.42279407296875*^9, 3.422794210907276*^9, { 3.422794722501026*^9, 3.422794744110401*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"N", "[", RowBox[{"Integrate", "[", RowBox[{ RowBox[{ RowBox[{"g", "[", "x", "]"}], " ", RowBox[{"f", "[", "x", "]"}]}], " ", ",", RowBox[{"{", RowBox[{"x", ",", "a", ",", "b"}], "}"}]}], "]"}], "]"}], " ", "+", " ", RowBox[{"NIntegrate", "[", RowBox[{ RowBox[{ RowBox[{"g", "[", "x", "]"}], " ", RowBox[{"h", "[", "x", "]"}]}], " ", ",", RowBox[{"{", RowBox[{"x", ",", "c", ",", "d"}], "}"}]}], "]"}], " "}]], "Input", CellChangeTimes->{{3.422263438342375*^9, 3.422263440358*^9}, 3.4222636715455*^9, 3.422263874108*^9, 3.42226412538925*^9, 3.422264310623625*^9, 3.4222656315455*^9, 3.422265726436125*^9, 3.422265944436125*^9, 3.422266123873625*^9, 3.42226616232675*^9, 3.4222668256705*^9, {3.422266872373625*^9, 3.422266900592375*^9}, 3.422267798311125*^9, 3.42226783401425*^9, 3.422354314217375*^9, 3.4223578821471252`*^9, 3.4223579395065002`*^9, 3.4223581933658752`*^9, 3.4223582244440002`*^9, 3.4223583286627502`*^9, 3.4223584581315002`*^9, 3.4223585118815002`*^9, 3.4223587513033752`*^9, 3.4223589582565002`*^9, 3.422792613625*^9, 3.422793082828125*^9, 3.422793249140625*^9, { 3.42279331228125*^9, 3.42279335275*^9}, 3.4227936966875*^9}], Cell[BoxData["1.8558514305672014`"], "Output", CellChangeTimes->{3.422794074703125*^9, 3.422794212985401*^9, 3.422794724751026*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"N", "[", "%", "]"}]], "Input", CellChangeTimes->{{3.422263438342375*^9, 3.422263440358*^9}, 3.4222636715455*^9, 3.422263874108*^9, 3.42226412538925*^9, 3.422264310623625*^9, 3.4222656315455*^9, 3.422265726436125*^9, 3.422265944436125*^9, 3.422266123873625*^9, 3.42226616232675*^9, 3.4222668256705*^9, {3.422266872373625*^9, 3.422266900592375*^9}, 3.422267798311125*^9, 3.42226783401425*^9, 3.422354314217375*^9, 3.4223578821471252`*^9, 3.4223579395065002`*^9, 3.4223581933658752`*^9, 3.4223582244440002`*^9, 3.4223583286627502`*^9, 3.4223584581315002`*^9, 3.4223585118815002`*^9, 3.4223587513033752`*^9, 3.4223589582565002`*^9, 3.422792613625*^9, 3.422793082828125*^9, 3.422793249140625*^9, { 3.42279331228125*^9, 3.42279335275*^9}, 3.4227936966875*^9}], Cell[BoxData["65.15302334728507`"], "Output", CellChangeTimes->{3.422793404609375*^9, 3.42279362484375*^9, 3.422793678578125*^9}], Cell[BoxData["65.15302334728507`"], "Output", CellChangeTimes->{3.422793404609375*^9, 3.42279362484375*^9, 3.422793678578125*^9}] }, Open ]], Cell[BoxData[""], "Input", CellChangeTimes->{{3.422263438342375*^9, 3.422263440358*^9}, 3.4222636715455*^9, 3.422263874108*^9, 3.42226412538925*^9, 3.422264310623625*^9, 3.4222656315455*^9, 3.422265726436125*^9, 3.422265944436125*^9, 3.422266123873625*^9, 3.42226616232675*^9, 3.4222668256705*^9, {3.422266872373625*^9, 3.422266900592375*^9}, 3.422267798311125*^9, 3.42226783401425*^9, 3.422354314217375*^9, 3.4223578821471252`*^9, 3.4223579395065002`*^9, 3.4223581933658752`*^9, 3.4223582244440002`*^9, 3.4223583286627502`*^9, 3.4223584581315002`*^9, 3.4223585118815002`*^9, 3.4223587513033752`*^9, 3.4223589582565002`*^9, 3.422792613625*^9, 3.422793082828125*^9, 3.422793249140625*^9, { 3.42279331228125*^9, 3.42279334584375*^9}}], Cell[BoxData[ RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{ RowBox[{ "##", "##", "##", "##", "##", "##", "##", "##", "##", "##", "cqwecqewcqewcweceqwcqc", "##", "##"}], "//", RowBox[{"/@", RowBox[{ RowBox[{"@", "sdcsdcwec"}], "@@@", "@@@"}]}]}], "//"}], "//"}], "*)"}]], "Input", CellChangeTimes->{{3.42226369838925*^9, 3.422263728608*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"THMP", "[", "m", "]"}]], "Input"], Cell[BoxData[ RowBox[{"THMP", "[", "9", "]"}]], "Output", CellChangeTimes->{3.422263835873625*^9, 3.4223589597408752`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{".55", "+", ".55", "+", ".89"}]], "Input"], Cell[BoxData["1.99`"], "Output", CellChangeTimes->{3.4223589603815002`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"x", "[", "0", "]"}]], "Input"], Cell[BoxData["1.217651892505655`"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"weights", "[", "1", "]"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"{", RowBox[{ "0.8957908187568364`", ",", "4.141768136045698`", ",", "8.371987626394912`", ",", "6.590453418802383`"}], "}"}], "[", "1", "]"}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"N", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"x", "[", "0", "]"}], "^", "2"}], " ", RowBox[{"Sin", "[", RowBox[{"x", "[", "0", "]"}], "]"}], " ", RowBox[{"weights", "[", RowBox[{"[", "1", "]"}], "]"}]}], " ", "+", " ", RowBox[{ RowBox[{ RowBox[{"x", "[", "1", "]"}], "^", "2"}], " ", RowBox[{"Sin", "[", RowBox[{"x", "[", "1", "]"}], "]"}], " ", RowBox[{"weights", "[", RowBox[{"[", "2", "]"}], "]"}]}], "+", RowBox[{ RowBox[{ RowBox[{"x", "[", "2", "]"}], "^", "2"}], " ", RowBox[{"Sin", "[", RowBox[{"x", "[", "2", "]"}], "]"}], " ", RowBox[{"weights", "[", RowBox[{"[", "3", "]"}], "]"}]}]}], "]"}]], "Input"], Cell[BoxData["46.5268752089569`"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"N", "[", RowBox[{"Integrate", "[", RowBox[{ RowBox[{ RowBox[{"x", "^", "2"}], " ", RowBox[{"Sin", "[", "x", "]"}]}], ",", RowBox[{"{", RowBox[{"x", ",", "1", ",", "3"}], "}"}]}], "]"}], "]"}], "\[IndentingNewLine]"}]], "Input"], Cell[BoxData["5.553423249078388`"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"-", RowBox[{"Cos", "[", "3", "]"}]}], "+", RowBox[{"Cos", "[", "1", "]"}]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"Cos", "[", "1", "]"}], "-", RowBox[{"Cos", "[", "3", "]"}]}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"N", "[", "%", "]"}]], "Input"], Cell[BoxData["1.5302948024685852`"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"N", "[", RowBox[{"364", "/", "3"}], "]"}]], "Input"], Cell[BoxData["121.33333333333333`"], "Output"] }, Open ]], Cell[BoxData[ RowBox[{"N", "["}]], "Input"], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{"(*", RowBox[{ "test", " ", "the", " ", "rule", " ", "with", " ", "some", " ", "examples"}], "*)"}]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"f", "[", "x_", "]"}], ":=", " ", RowBox[{"Sin", "[", "x", "]"}]}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"G", "[", "1", "]"}]], "Input"], Cell[BoxData["1.530350795286579`"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"NIntegrate", "[", RowBox[{ RowBox[{"f", "[", "x", "]"}], ",", RowBox[{"{", RowBox[{"x", ",", "1", ",", "3"}], "}"}]}], "]"}]], "Input"], Cell[BoxData["1.5302948024685852`"], "Output"] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"f", "[", "x_", "]"}], ":=", " ", RowBox[{"x", "^", "6"}], " ", RowBox[{"(*", RowBox[{ RowBox[{"polynomial", " ", "of", " ", "degree"}], " ", ">", " ", "5"}], "*)"}]}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"G", "[", "1", "]"}]], "Input"], Cell[BoxData[ RowBox[{"G", "[", "1", "]"}]], "Output", CellChangeTimes->{3.4222634226705*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"NIntegrate", "[", RowBox[{ RowBox[{"f", "[", "x", "]"}], ",", RowBox[{"{", RowBox[{"x", ",", "1", ",", "3"}], "}"}]}], "]"}]], "Input"], Cell[BoxData["312.2857142857143`"], "Output"] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"f", "[", "x_", "]"}], ":=", " ", RowBox[{"Exp", "[", RowBox[{"x", "^", "2"}], "]"}], " ", RowBox[{"(*", RowBox[{"rapidly", " ", "growing", " ", "function"}], "*)"}]}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"G", "[", "1", "]"}]], "Input"], Cell[BoxData["1275.9078061041184`"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"NIntegrate", "[", RowBox[{ RowBox[{"f", "[", "x", "]"}], ",", RowBox[{"{", RowBox[{"x", ",", "1", ",", "3"}], "}"}]}], "]"}]], "Input"], Cell[BoxData["1443.0824711468072`"], "Output"] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{"increase", " ", "the", " ", "degree"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"m", "=", "9"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"a", "=", "1"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"b", "=", "3"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "0"}], ",", RowBox[{"i", "\[LessEqual]", "m"}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"bb", "[", "i", "]"}], " ", "=", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"b", "^", RowBox[{"(", RowBox[{"i", "+", "1"}], ")"}]}], "-", RowBox[{"a", "^", RowBox[{"(", RowBox[{"i", "+", "1"}], ")"}]}]}], ")"}], "/", RowBox[{"(", RowBox[{"i", "+", "1"}], ")"}]}]}]}], "]"}], ";"}], "\[IndentingNewLine]"}]}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"THMP", "[", "m", "]"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData["\<\"Moment matrix :\"\>"], "Print"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"2", "4", FractionBox["26", "3"], "20", FractionBox["242", "5"], FractionBox["364", "3"]}, {"4", FractionBox["26", "3"], "20", FractionBox["242", "5"], FractionBox["364", "3"], FractionBox["2186", "7"]}, { FractionBox["26", "3"], "20", FractionBox["242", "5"], FractionBox["364", "3"], FractionBox["2186", "7"], "820"}, {"20", FractionBox["242", "5"], FractionBox["364", "3"], FractionBox["2186", "7"], "820", FractionBox["19682", "9"]}, { FractionBox["242", "5"], FractionBox["364", "3"], FractionBox["2186", "7"], "820", FractionBox["19682", "9"], FractionBox["29524", "5"]} }], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Print"], Cell[BoxData["\<\"A measure exists\"\>"], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"nodes\"\>", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"weights\"\>"}], SequenceForm["nodes", " ", "weights"], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "1.093820154061336`", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "0.23692688505618908`"}], SequenceForm[1.093820154061336, " ", 0.23692688505618908`], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "1.461530689894317`", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "0.47862867049936647`"}], SequenceForm[1.461530689894317, " ", 0.47862867049936647`], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "2.`", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "0.5688888888888889`"}], SequenceForm[2., " ", 0.5688888888888889], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "2.538469310105683`", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "0.47862867049936647`"}], SequenceForm[2.538469310105683, " ", 0.47862867049936647`], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "2.906179845938664`", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "0.23692688505618908`"}], SequenceForm[2.906179845938664, " ", 0.23692688505618908`], Editable->False]], "Print"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"NIntegrate", "[", RowBox[{ RowBox[{"f", "[", "x", "]"}], ",", RowBox[{"{", RowBox[{"x", ",", "1", ",", "3"}], "}"}]}], "]"}], " "}]], "Input"], Cell[BoxData["1443.0824711468072`"], "Output"] }, Open ]], Cell[BoxData[ RowBox[{"(*", RowBox[{"much", " ", "better", " ", "estimate"}], "*)"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ "Some", " ", "examples", " ", "with", " ", "moment", " ", "problems"}], "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"m", "=", "4"}]}]], "Input"], Cell[BoxData["4"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"bb", "[", "0", "]"}], "=", RowBox[{ RowBox[{"bb", "[", "1", "]"}], "=", RowBox[{ RowBox[{"bb", "[", "2", "]"}], "=", RowBox[{ RowBox[{"bb", "[", "3", "]"}], "=", "1"}]}]}]}]], "Input"], Cell[BoxData["1"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"bb", "[", "4", "]"}], "=", "2"}]], "Input"], Cell[BoxData["2"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"THMP", "[", "m", "]"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData["\<\"Moment matrix :\"\>"], "Print"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"1", "1", "1"}, {"1", "1", "1"}, {"1", "1", "2"} }], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Print"], Cell[BoxData["\<\"no measure: moment data not recursively generated\"\>"], \ "Print"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"bb", "[", "4", "]"}], "=", "1"}]], "Input"], Cell[BoxData["1"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"THMP", "[", "m", "]"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData["\<\"Moment matrix :\"\>"], "Print"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"1", "1", "1"}, {"1", "1", "1"}, {"1", "1", "1"} }], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Print"], Cell[BoxData["\<\"A measure exists\"\>"], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"nodes\"\>", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"weights\"\>"}], SequenceForm["nodes", " ", "weights"], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "1.`", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.`"}], SequenceForm[1., " ", 1.], Editable->False]], "Print"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"bb", "[", "2", "]"}], "=", "2"}]}]], "Input"], Cell[BoxData["2"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"bb", "[", "3", "]"}], "=", "3"}]], "Input"], Cell[BoxData["3"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"bb", "[", "4", "]"}], "=", "5"}]], "Input"], Cell[BoxData["5"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"THMP", "[", "m", "]"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData["\<\"Moment matrix :\"\>"], "Print"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"1", "1", "2"}, {"1", "2", "3"}, {"2", "3", "5"} }], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Print"], Cell[BoxData["\<\"A measure exists\"\>"], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"nodes\"\>", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"weights\"\>"}], SequenceForm["nodes", " ", "weights"], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ RowBox[{"-", "0.6180339887498949`"}], "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "0.276393202250021`"}], SequenceForm[-0.6180339887498949, " ", 0.276393202250021], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "1.618033988749895`", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "0.7236067977499789`"}], SequenceForm[1.618033988749895, " ", 0.7236067977499789], Editable->False]], "Print"] }, Open ]] }, Open ]], Cell[BoxData[ RowBox[{"(*", RowBox[{ RowBox[{"H", " ", "positive"}], ",", " ", RowBox[{ RowBox[{"rank", " ", "2"}], ";", " ", RowBox[{"2", "-", RowBox[{"atomic", " ", "measure"}]}]}]}], "*)"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"bb", "[", "4", "]"}], "=", "4"}]], "Input"], Cell[BoxData["4"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"THMP", "[", "m", "]"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData["\<\"Moment matrix :\"\>"], "Print"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"1", "1", "2"}, {"1", "2", "3"}, {"2", "3", "4"} }], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Print"], Cell[BoxData["\<\"no measure: moment matrix not positive semidefinite\"\>"], \ "Print"] }, Open ]] }, Open ]], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"bb", "[", "4", "]"}], "=", "6"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"bb", "[", "5", "]"}], "=", " ", "0"}], ";"}], " ", RowBox[{"(*", RowBox[{ "supply", " ", "an", " ", "extra", " ", "moment", " ", "to", " ", "determine", " ", "a", " ", "particular", " ", "representing", " ", "measure"}], "*)"}]}]}], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"THMP", "[", "m", "]"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData["\<\"Moment matrix :\"\>"], "Print"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"1", "1", "2"}, {"1", "2", "3"}, {"2", "3", "6"} }], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Print"], Cell[BoxData["\<\"A measure exists\"\>"], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"nodes\"\>", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"weights\"\>"}], SequenceForm["nodes", " ", "weights"], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ RowBox[{"-", "10.099901985167726`"}], "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "0.00008019774726564291`"}], SequenceForm[-10.099901985167726`, " ", 0.00008019774726564291], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ RowBox[{"-", "0.5425433326913511`"}], "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "0.293595853642608`"}], SequenceForm[-0.5425433326913511, " ", 0.293595853642608], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "1.6424453178590765`", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "0.7063239486101264`"}], SequenceForm[1.6424453178590765`, " ", 0.7063239486101264], Editable->False]], "Print"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"bb", "[", "5", "]"}], " ", "=", " ", RowBox[{"-", "100"}]}]], "Input"], Cell[BoxData[ RowBox[{"-", "100"}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"THMP", "[", "m", "]"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData["\<\"Moment matrix :\"\>"], "Print"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"1", "1", "2"}, {"1", "2", "3"}, {"2", "3", "6"} }], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Print"], Cell[BoxData["\<\"A measure exists\"\>"], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"nodes\"\>", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"weights\"\>"}], SequenceForm["nodes", " ", "weights"], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ RowBox[{"-", "110.0090909022619`"}], "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "6.705960509706324`*^-9"}], SequenceForm[-110.0090909022619, " ", 6.705960509706324*^-9], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ RowBox[{"-", "0.611426464591715`"}], "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "0.27801623388884633`"}], SequenceForm[-0.611426464591715, " ", 0.27801623388884633`], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "1.6205173668536184`", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "0.7219837594051931`"}], SequenceForm[1.6205173668536184`, " ", 0.7219837594051931], Editable->False]], "Print"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"m", "=", "5"}]], "Input"], Cell[BoxData["5"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"bb", "[", "4", "]"}], "=", "5"}]], "Input"], Cell[BoxData["5"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"bb", "[", "5", "]"}], "=", "8"}]], "Input"], Cell[BoxData["8"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"THMP", "[", "m", "]"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData["\<\"Moment matrix :\"\>"], "Print"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"1", "1", "2", "3"}, {"1", "2", "3", "5"}, {"2", "3", "5", "8"} }], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Print"], Cell[BoxData["\<\"A measure exists\"\>"], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"nodes\"\>", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"weights\"\>"}], SequenceForm["nodes", " ", "weights"], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ RowBox[{"-", "0.6180339887498949`"}], "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "0.276393202250021`"}], SequenceForm[-0.6180339887498949, " ", 0.276393202250021], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "1.618033988749895`", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "0.7236067977499789`"}], SequenceForm[1.618033988749895, " ", 0.7236067977499789], Editable->False]], "Print"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"bb", "[", "5", "]"}], "=", "9"}]], "Input"], Cell[BoxData["9"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"THMP", "[", "m", "]"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData["\<\"Moment matrix :\"\>"], "Print"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"1", "1", "2", "3"}, {"1", "2", "3", "5"}, {"2", "3", "5", "9"} }], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Print"], Cell[BoxData["\<\"no measure: moment data not recursively generated\"\>"], \ "Print"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"bb", "[", "4", "]"}], "=", "4"}]], "Input"], Cell[BoxData["4"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"THMP", "[", "m", "]"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData["\<\"Moment matrix :\"\>"], "Print"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"1", "1", "2", "3"}, {"1", "2", "3", "4"}, {"2", "3", "4", "9"} }], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Print"], Cell[BoxData["\<\"no measure: moment matrix not positive semidefinite\"\>"], \ "Print"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"bb", "[", "4", "]"}], "=", "6"}]], "Input"], Cell[BoxData["6"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData["m"], "Input"], Cell[BoxData["5"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"THMP", "[", "m", "]"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData["\<\"Moment matrix :\"\>"], "Print"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"1", "1", "2", "3"}, {"1", "2", "3", "6"}, {"2", "3", "6", "9"} }], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Print"], Cell[BoxData["\<\"A measure exists\"\>"], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"nodes\"\>", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"weights\"\>"}], SequenceForm["nodes", " ", "weights"], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ RowBox[{"-", "1.7320508075688772`"}], "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "0.04465819873852045`"}], SequenceForm[-1.7320508075688772`, " ", 0.04465819873852045], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "0.`", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "0.3333333333333333`"}], SequenceForm[0., " ", 0.3333333333333333], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "1.7320508075688772`", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "0.6220084679281462`"}], SequenceForm[1.7320508075688772`, " ", 0.6220084679281462], Editable->False]], "Print"] }, Open ]] }, Open ]], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"bb", "[", "0", "]"}], "=", "1"}], ";", " ", RowBox[{ RowBox[{"bb", "[", "1", "]"}], "=", "1"}], ";", RowBox[{ RowBox[{"bb", "[", "2", "]"}], "=", "2"}], ";", RowBox[{ RowBox[{"bb", "[", "3", "]"}], "=", "3"}], ";", RowBox[{ RowBox[{"bb", "[", "4", "]"}], "=", "6"}], ";", " ", RowBox[{ RowBox[{"bb", "[", "5", "]"}], "=", "1"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"bb", "[", "6", "]"}], "=", "82"}], ";", RowBox[{ RowBox[{"bb", "[", "7", "]"}], "=", " ", RowBox[{"-", "597"}]}], ";", " ", RowBox[{ RowBox[{"bb", "[", "8", "]"}], "=", "5686"}], ";", RowBox[{ RowBox[{"bb", "[", "9", "]"}], "=", RowBox[{"-", "51399"}]}], ";", RowBox[{ RowBox[{"bb", "[", "10", "]"}], "=", " ", "468962"}], ";"}], " "}]}], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"m", "=", "10"}]], "Input"], Cell[BoxData["10"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"THMP", "[", "m", "]"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData["\<\"Moment matrix :\"\>"], "Print"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"1", "1", "2", "3", "6", "1"}, {"1", "2", "3", "6", "1", "82"}, {"2", "3", "6", "1", "82", RowBox[{"-", "597"}]}, {"3", "6", "1", "82", RowBox[{"-", "597"}], "5686"}, {"6", "1", "82", RowBox[{"-", "597"}], "5686", RowBox[{"-", "51399"}]}, {"1", "82", RowBox[{"-", "597"}], "5686", RowBox[{"-", "51399"}], "468962"} }], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Print"], Cell[BoxData["\<\"A measure exists\"\>"], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"nodes\"\>", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"weights\"\>"}], SequenceForm["nodes", " ", "weights"], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ RowBox[{"-", "9.110962423069871`"}], "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "0.00011895988516730375`"}], SequenceForm[-9.110962423069871, " ", 0.00011895988516730375`], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ RowBox[{"-", "0.5338408190584355`"}], "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "0.2953780972966052`"}], SequenceForm[-0.5338408190584355, " ", 0.2953780972966052], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{ "1.6448032421283072`", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "0.7045029428182275`"}], SequenceForm[1.6448032421283072`, " ", 0.7045029428182275], Editable->False]], "Print"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"bb", "[", "10", "]"}], "=", " ", "468963"}]], "Input"], Cell[BoxData["468963"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"THMP", "[", "m", "]"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData["\<\"Moment matrix :\"\>"], "Print"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"1", "1", "2", "3", "6", "1"}, {"1", "2", "3", "6", "1", "82"}, {"2", "3", "6", "1", "82", RowBox[{"-", "597"}]}, {"3", "6", "1", "82", RowBox[{"-", "597"}], "5686"}, {"6", "1", "82", RowBox[{"-", "597"}], "5686", RowBox[{"-", "51399"}]}, {"1", "82", RowBox[{"-", "597"}], "5686", RowBox[{"-", "51399"}], "468963"} }], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Print"], Cell[BoxData["\<\"no measure: moment data not recursively generated\"\>"], \ "Print"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"bb", "[", "10", "]"}], "=", "468962"}]], "Input"], Cell[BoxData["468962"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"bb", "[", "6", "]"}], "=", "81"}]], "Input"], Cell[BoxData["81"], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"THMP", "[", "m", "]"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData["\<\"Moment matrix :\"\>"], "Print"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"1", "1", "2", "3", "6", "1"}, {"1", "2", "3", "6", "1", "81"}, {"2", "3", "6", "1", "81", RowBox[{"-", "597"}]}, {"3", "6", "1", "81", RowBox[{"-", "597"}], "5686"}, {"6", "1", "81", RowBox[{"-", "597"}], "5686", RowBox[{"-", "51399"}]}, {"1", "81", RowBox[{"-", "597"}], "5686", RowBox[{"-", "51399"}], "468962"} }], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Print"], Cell[BoxData["\<\"no measure: moment matrix not positive semidefinite\"\>"], \ "Print"] }, Open ]] }, Open ]] }, WindowSize->{1272, 905}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, DockedCells->(FrontEndExecute[{ FrontEnd`NotebookApply[ FrontEnd`InputNotebook[], #, Placeholder]}]& ), PrintingCopies->1, PrintingPageRange->{Automatic, Automatic}, FrontEndVersion->"6.0 for Microsoft Windows (32-bit) (June 19, 2007)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[568, 21, 12002, 275, 1452, "Input"], Cell[12573, 298, 1555, 24, 132, "Input"], Cell[14131, 324, 709, 15, 72, "Input"], Cell[14843, 341, 659, 11, 31, "Input"], Cell[CellGroupData[{ Cell[15527, 356, 491, 8, 31, "Input"], Cell[16021, 366, 481, 8, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[16539, 379, 119, 2, 31, "Input"], Cell[16661, 383, 898, 13, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[17596, 401, 641, 12, 31, "Input"], Cell[18240, 415, 13757, 228, 231, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[32034, 648, 3944, 83, 212, "Input"], Cell[CellGroupData[{ Cell[36003, 735, 381, 5, 23, "Print"], Cell[36387, 742, 380, 5, 23, "Print"], Cell[36770, 749, 381, 5, 23, "Print"], Cell[37154, 756, 381, 5, 23, "Print"], Cell[37538, 763, 382, 5, 23, "Print"], Cell[37923, 770, 381, 5, 23, "Print"], Cell[38307, 777, 382, 5, 23, "Print"], Cell[38692, 784, 385, 5, 23, "Print"] }, Open ]] }, Open ]], Cell[39104, 793, 1636, 46, 192, "Input"], Cell[40743, 841, 26, 0, 31, "Input"], Cell[CellGroupData[{ Cell[40794, 845, 2996, 75, 252, "Input"], Cell[CellGroupData[{ Cell[43815, 924, 1054, 15, 23, "Print"], Cell[44872, 941, 2057, 38, 69, "Print"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[46978, 985, 1713, 46, 192, "Input"], Cell[CellGroupData[{ Cell[48716, 1035, 1058, 15, 23, "Print"], Cell[49777, 1052, 1058, 15, 23, "Print"], Cell[50838, 1069, 1059, 15, 23, "Print"], Cell[51900, 1086, 1059, 15, 23, "Print"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[53008, 1107, 4369, 103, 392, "Input"], Cell[57380, 1212, 1966, 38, 86, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[59383, 1255, 1246, 35, 112, "Input"], Cell[60632, 1292, 2057, 38, 90, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[62726, 1335, 619, 17, 112, "Input"], Cell[63348, 1354, 1057, 15, 23, "Print"] }, Open ]], Cell[CellGroupData[{ Cell[64442, 1374, 2176, 63, 232, "Input"], Cell[CellGroupData[{ Cell[66643, 1441, 1053, 15, 23, "Print"], Cell[67699, 1458, 1206, 20, 23, "Print"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[68954, 1484, 1385, 41, 152, "Input"], Cell[70342, 1527, 1693, 32, 23, "Print"] }, Open ]], Cell[CellGroupData[{ Cell[72072, 1564, 343, 10, 52, "Input"], Cell[72418, 1576, 1767, 33, 23, "Print"] }, Open ]], Cell[CellGroupData[{ Cell[74222, 1614, 2333, 71, 212, "Input"], Cell[76558, 1687, 1876, 37, 86, "Output"], Cell[78437, 1726, 1016, 14, 23, "Print"] }, Open ]], Cell[CellGroupData[{ Cell[79490, 1745, 1162, 38, 132, "Input"], Cell[CellGroupData[{ Cell[80677, 1787, 453, 9, 23, "Print"], Cell[81133, 1798, 494, 10, 23, "Print"], Cell[81630, 1810, 488, 10, 23, "Print"], Cell[82121, 1822, 488, 10, 23, "Print"], Cell[82612, 1834, 488, 10, 23, "Print"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[83149, 1850, 820, 22, 31, "Input"], Cell[83972, 1874, 994, 15, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[85003, 1894, 259, 4, 31, "Input"], Cell[85265, 1900, 1002, 15, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[86304, 1920, 119, 2, 31, "Input"], Cell[86426, 1924, 947, 14, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[87410, 1943, 129, 2, 31, "Input"], Cell[87542, 1947, 1972, 37, 90, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[89551, 1989, 198, 3, 31, "Input"], Cell[89752, 1994, 796, 13, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[90585, 2012, 54, 1, 31, "Input"], Cell[90642, 2015, 974, 14, 30, "Output"] }, Open ]], Cell[91631, 2032, 2013, 53, 192, "Input"], Cell[93647, 2087, 1134, 32, 112, "Input"], Cell[CellGroupData[{ Cell[94806, 2123, 54, 1, 31, "Input"], Cell[94863, 2126, 276, 5, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[95176, 2136, 1295, 28, 31, "Input"], Cell[96474, 2166, 135, 2, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[96646, 2173, 817, 12, 31, "Input"], Cell[97466, 2187, 133, 2, 30, "Output"], Cell[97602, 2191, 133, 2, 30, "Output"] }, Open ]], Cell[97750, 2196, 772, 11, 31, "Input"], Cell[98525, 2209, 387, 12, 31, "Input"], Cell[CellGroupData[{ Cell[98937, 2225, 57, 1, 31, "Input"], Cell[98997, 2228, 124, 2, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[99158, 2235, 65, 1, 31, "Input"], Cell[99226, 2238, 76, 1, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[99339, 2244, 54, 1, 31, "Input"], Cell[99396, 2247, 45, 0, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[99478, 2252, 60, 1, 31, "Input"], Cell[99541, 2255, 203, 6, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[99781, 2266, 743, 23, 31, "Input"], Cell[100527, 2291, 44, 0, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[100608, 2296, 305, 10, 52, "Input"], Cell[100916, 2308, 45, 0, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[100998, 2313, 128, 4, 31, "Input"], Cell[101129, 2319, 110, 3, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[101276, 2327, 54, 1, 31, "Input"], Cell[101333, 2330, 46, 0, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[101416, 2335, 79, 2, 31, "Input"], Cell[101498, 2339, 46, 0, 30, "Output"] }, Open ]], Cell[101559, 2342, 44, 1, 31, "Input"], Cell[101606, 2345, 181, 5, 52, "Input"], Cell[101790, 2352, 114, 3, 31, "Input"], Cell[CellGroupData[{ Cell[101929, 2359, 54, 1, 31, "Input"], Cell[101986, 2362, 45, 0, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[102068, 2367, 174, 5, 31, "Input"], Cell[102245, 2374, 46, 0, 30, "Output"] }, Open ]], Cell[102306, 2377, 233, 7, 31, "Input"], Cell[CellGroupData[{ Cell[102564, 2388, 54, 1, 31, "Input"], Cell[102621, 2391, 95, 2, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[102753, 2398, 174, 5, 31, "Input"], Cell[102930, 2405, 45, 0, 30, "Output"] }, Open ]], Cell[102990, 2408, 230, 6, 31, "Input"], Cell[CellGroupData[{ Cell[103245, 2418, 54, 1, 31, "Input"], Cell[103302, 2421, 46, 0, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[103385, 2426, 174, 5, 31, "Input"], Cell[103562, 2433, 46, 0, 30, "Output"] }, Open ]], Cell[103623, 2436, 1016, 31, 152, "Input"], Cell[CellGroupData[{ Cell[104664, 2471, 57, 1, 31, "Input"], Cell[CellGroupData[{ Cell[104746, 2476, 49, 0, 23, "Print"], Cell[104798, 2478, 894, 30, 125, "Print"], Cell[105695, 2510, 50, 0, 23, "Print"], Cell[105748, 2512, 213, 5, 23, "Print"], Cell[105964, 2519, 251, 6, 23, "Print"], Cell[106218, 2527, 251, 6, 23, "Print"], Cell[106472, 2535, 218, 6, 23, "Print"], Cell[106693, 2543, 251, 6, 23, "Print"], Cell[106947, 2551, 251, 6, 23, "Print"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[107247, 2563, 200, 6, 31, "Input"], Cell[107450, 2571, 46, 0, 30, "Output"] }, Open ]], Cell[107511, 2574, 99, 2, 31, "Input"], Cell[CellGroupData[{ Cell[107635, 2580, 225, 6, 72, "Input"], Cell[107863, 2588, 28, 0, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[107928, 2593, 243, 8, 31, "Input"], Cell[108174, 2603, 28, 0, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[108239, 2608, 78, 2, 31, "Input"], Cell[108320, 2612, 28, 0, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[108385, 2617, 57, 1, 31, "Input"], Cell[CellGroupData[{ Cell[108467, 2622, 49, 0, 23, "Print"], Cell[108519, 2624, 218, 8, 51, "Print"], Cell[108740, 2634, 85, 1, 23, "Print"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[108874, 2641, 78, 2, 31, "Input"], Cell[108955, 2645, 28, 0, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[109020, 2650, 57, 1, 31, "Input"], Cell[CellGroupData[{ Cell[109102, 2655, 49, 0, 23, "Print"], Cell[109154, 2657, 218, 8, 51, "Print"], Cell[109375, 2667, 50, 0, 23, "Print"], Cell[109428, 2669, 213, 5, 23, "Print"], Cell[109644, 2676, 182, 5, 23, "Print"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[109875, 2687, 138, 3, 72, "Input"], Cell[110016, 2692, 28, 0, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[110081, 2697, 78, 2, 31, "Input"], Cell[110162, 2701, 28, 0, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[110227, 2706, 78, 2, 31, "Input"], Cell[110308, 2710, 28, 0, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[110373, 2715, 57, 1, 31, "Input"], Cell[CellGroupData[{ Cell[110455, 2720, 49, 0, 23, "Print"], Cell[110507, 2722, 218, 8, 51, "Print"], Cell[110728, 2732, 50, 0, 23, "Print"], Cell[110781, 2734, 213, 5, 23, "Print"], Cell[110997, 2741, 265, 6, 23, "Print"], Cell[111265, 2749, 248, 6, 23, "Print"] }, Open ]] }, Open ]], Cell[111540, 2759, 229, 7, 31, "Input"], Cell[CellGroupData[{ Cell[111794, 2770, 78, 2, 31, "Input"], Cell[111875, 2774, 28, 0, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[111940, 2779, 57, 1, 31, "Input"], Cell[CellGroupData[{ Cell[112022, 2784, 49, 0, 23, "Print"], Cell[112074, 2786, 218, 8, 51, "Print"], Cell[112295, 2796, 87, 1, 23, "Print"] }, Open ]] }, Open ]], Cell[112409, 2801, 420, 12, 52, "Input"], Cell[CellGroupData[{ Cell[112854, 2817, 57, 1, 31, "Input"], Cell[CellGroupData[{ Cell[112936, 2822, 49, 0, 23, "Print"], Cell[112988, 2824, 218, 8, 51, "Print"], Cell[113209, 2834, 50, 0, 23, "Print"], Cell[113262, 2836, 213, 5, 23, "Print"], Cell[113478, 2843, 276, 6, 23, "Print"], Cell[113757, 2851, 265, 6, 23, "Print"], Cell[114025, 2859, 251, 6, 23, "Print"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[114325, 2871, 108, 3, 31, "Input"], Cell[114436, 2876, 47, 1, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[114520, 2882, 57, 1, 31, "Input"], Cell[CellGroupData[{ Cell[114602, 2887, 49, 0, 23, "Print"], Cell[114654, 2889, 218, 8, 51, "Print"], Cell[114875, 2899, 50, 0, 23, "Print"], Cell[114928, 2901, 213, 5, 23, "Print"], Cell[115144, 2908, 271, 6, 23, "Print"], Cell[115418, 2916, 268, 6, 23, "Print"], Cell[115689, 2924, 251, 6, 23, "Print"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[115989, 2936, 49, 1, 31, "Input"], Cell[116041, 2939, 28, 0, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[116106, 2944, 78, 2, 31, "Input"], Cell[116187, 2948, 28, 0, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[116252, 2953, 78, 2, 31, "Input"], Cell[116333, 2957, 28, 0, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[116398, 2962, 57, 1, 31, "Input"], Cell[CellGroupData[{ Cell[116480, 2967, 49, 0, 23, "Print"], Cell[116532, 2969, 233, 8, 51, "Print"], Cell[116768, 2979, 50, 0, 23, "Print"], Cell[116821, 2981, 213, 5, 23, "Print"], Cell[117037, 2988, 265, 6, 23, "Print"], Cell[117305, 2996, 248, 6, 23, "Print"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[117602, 3008, 78, 2, 31, "Input"], Cell[117683, 3012, 28, 0, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[117748, 3017, 57, 1, 31, "Input"], Cell[CellGroupData[{ Cell[117830, 3022, 49, 0, 23, "Print"], Cell[117882, 3024, 233, 8, 51, "Print"], Cell[118118, 3034, 85, 1, 23, "Print"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[118252, 3041, 78, 2, 31, "Input"], Cell[118333, 3045, 28, 0, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[118398, 3050, 57, 1, 31, "Input"], Cell[CellGroupData[{ Cell[118480, 3055, 49, 0, 23, "Print"], Cell[118532, 3057, 233, 8, 51, "Print"], Cell[118768, 3067, 87, 1, 23, "Print"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[118904, 3074, 78, 2, 31, "Input"], Cell[118985, 3078, 28, 0, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[119050, 3083, 27, 0, 31, "Input"], Cell[119080, 3085, 28, 0, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[119145, 3090, 57, 1, 31, "Input"], Cell[CellGroupData[{ Cell[119227, 3095, 49, 0, 23, "Print"], Cell[119279, 3097, 233, 8, 51, "Print"], Cell[119515, 3107, 50, 0, 23, "Print"], Cell[119568, 3109, 213, 5, 23, "Print"], Cell[119784, 3116, 270, 6, 23, "Print"], Cell[120057, 3124, 218, 6, 23, "Print"], Cell[120278, 3132, 251, 6, 23, "Print"] }, Open ]] }, Open ]], Cell[120556, 3142, 876, 28, 52, "Input"], Cell[CellGroupData[{ Cell[121457, 3174, 50, 1, 31, "Input"], Cell[121510, 3177, 29, 0, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[121576, 3182, 57, 1, 31, "Input"], Cell[CellGroupData[{ Cell[121658, 3187, 49, 0, 23, "Print"], Cell[121710, 3189, 541, 17, 95, "Print"], Cell[122254, 3208, 50, 0, 23, "Print"], Cell[122307, 3210, 213, 5, 23, "Print"], Cell[122523, 3217, 274, 6, 23, "Print"], Cell[122800, 3225, 267, 6, 23, "Print"], Cell[123070, 3233, 251, 6, 23, "Print"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[123370, 3245, 89, 2, 31, "Input"], Cell[123462, 3249, 33, 0, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[123532, 3254, 57, 1, 31, "Input"], Cell[CellGroupData[{ Cell[123614, 3259, 49, 0, 23, "Print"], Cell[123666, 3261, 541, 17, 95, "Print"], Cell[124210, 3280, 85, 1, 23, "Print"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[124344, 3287, 84, 2, 31, "Input"], Cell[124431, 3291, 33, 0, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[124501, 3296, 79, 2, 31, "Input"], Cell[124583, 3300, 29, 0, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[124649, 3305, 57, 1, 31, "Input"], Cell[CellGroupData[{ Cell[124731, 3310, 49, 0, 23, "Print"], Cell[124783, 3312, 541, 17, 95, "Print"], Cell[125327, 3331, 87, 1, 23, "Print"] }, Open ]] }, Open ]] } ] *) (* End of internal cache information *)