(* 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[ 36024, 827] NotebookOptionsPosition[ 35507, 805] NotebookOutlinePosition[ 35841, 820] CellTagsIndexPosition[ 35798, 817] WindowFrame->Normal ContainsDynamic->False*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"ctor", "[", "z_", "]"}], ":=", RowBox[{"(", RowBox[{ RowBox[{"output", "=", RowBox[{"{", "}"}]}], ";", "\[IndentingNewLine]", RowBox[{"Do", "[", RowBox[{ RowBox[{"AppendTo", "[", RowBox[{"output", ",", RowBox[{"{", RowBox[{ RowBox[{"Re", "[", RowBox[{"z", "[", RowBox[{"[", "k", "]"}], "]"}], "]"}], ",", RowBox[{"Im", "[", RowBox[{"z", "[", RowBox[{"[", "k", "]"}], "]"}], "]"}]}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"k", ",", "1", ",", RowBox[{"Length", "[", "z", "]"}]}], "}"}]}], "]"}], ";", "\[IndentingNewLine]", "output"}], ")"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]"}], "\n", RowBox[{ RowBox[{ RowBox[{ RowBox[{"plotpolylist", "[", "vertlist_", "]"}], ":=", RowBox[{"Show", "[", RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics", "[", RowBox[{ RowBox[{"{", RowBox[{"(*", RowBox[{ RowBox[{"GrayLevel", "[", "1", "]"}], ",", RowBox[{"Thickness", "[", ".01", "]"}], ",", RowBox[{"{", RowBox[{"Polygon", "[", RowBox[{"vertlist", "[", RowBox[{"[", "k", "]"}], "]"}], "]"}], "}"}], ","}], "*)"}], RowBox[{ RowBox[{"GrayLevel", "[", RowBox[{"If", "[", RowBox[{ RowBox[{"k", "\[Equal]", RowBox[{"Length", "[", "vertlist", "]"}]}], ",", "0", ",", ".1"}], "]"}], "]"}], ",", RowBox[{"Thickness", "[", RowBox[{"If", "[", RowBox[{ RowBox[{"k", "\[Equal]", RowBox[{"Length", "[", "vertlist", "]"}]}], ",", ".005", ",", ".005"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"Line", "[", RowBox[{"vertlist", "[", RowBox[{"[", "k", "]"}], "]"}], "]"}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"PlotRange", "\[Rule]", "All"}], ",", RowBox[{"Axes", "\[Rule]", "None"}], ",", RowBox[{"AspectRatio", "\[Rule]", "Automatic"}]}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"k", ",", "1", ",", RowBox[{"Length", "[", "vertlist", "]"}]}], "}"}]}], "]"}], "]"}]}], ";"}], "\n"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"plotpolylistC2", "[", "vert_", "]"}], ":=", RowBox[{"(", RowBox[{ RowBox[{"vert3", "=", RowBox[{"{", "}"}]}], ";", "\[IndentingNewLine]", RowBox[{"Do", "[", RowBox[{ RowBox[{ RowBox[{"vert1", "=", RowBox[{"ctor", "[", RowBox[{"vert", "[", RowBox[{"[", "k", "]"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"AppendTo", "[", RowBox[{"vert3", ",", "vert1"}], "]"}], ";"}], ",", RowBox[{"{", RowBox[{"k", ",", "1", ",", RowBox[{"Length", "[", "vert", "]"}]}], "}"}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"plotpolylist", "[", "vert3", "]"}]}], ")"}]}], ";"}], "\n"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"(*", RowBox[{ "compute", " ", "the", " ", "tayor", " ", "series", " ", "coefficients", " ", "of", " ", RowBox[{ RowBox[{"(", RowBox[{"1", "+", "z"}], ")"}], "^", "p"}], " ", "The", " ", "kth", " ", "item", " ", "of", " ", "the", " ", "list", " ", "the", " ", "coefficient", " ", "of", " ", RowBox[{"z", "^", RowBox[{"{", RowBox[{"k", "-", "1"}], "}"}]}]}], "*)"}], RowBox[{ RowBox[{ RowBox[{"makecoeff", "[", RowBox[{"p_", ",", "z_", ",", "n_"}], "]"}], ":=", RowBox[{"(", RowBox[{ RowBox[{"out", "=", RowBox[{"{", "1", "}"}]}], ";", "\[IndentingNewLine]", RowBox[{"Do", "[", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"AppendTo", "[", RowBox[{"out", ",", RowBox[{ RowBox[{"Last", "[", "out", "]"}], " ", RowBox[{ RowBox[{"(", RowBox[{"p", "-", "k", "+", "1"}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{"(", "k", ")"}], RowBox[{"(", RowBox[{"-", "z"}], ")"}]}], ")"}]}]}]}], "]"}], ";"}], ")"}], ",", RowBox[{"{", RowBox[{"k", ",", "1", ",", "n"}], "}"}]}], "]"}], ";", "\[IndentingNewLine]", "out"}], ")"}]}], ";"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"(*", " ", RowBox[{ "evaluate", " ", "the", " ", "poser", " ", "series", " ", "with", " ", "given", " ", "coefficients", " ", "at", " ", "given", " ", "point"}], " ", "*)"}], RowBox[{"\[AliasDelimiter]", "\[IndentingNewLine]", "\n", RowBox[{ RowBox[{ RowBox[{"sumcoeff", "[", RowBox[{"coeff_", ",", "w_"}], "]"}], ":=", RowBox[{"(", RowBox[{ RowBox[{"wout", "=", RowBox[{"First", "[", "coeff", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"wpower", "=", "1"}], ";", "\[IndentingNewLine]", RowBox[{"Do", "[", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"wpower", "=", RowBox[{"wpower", " ", "w"}]}], ";", "\[IndentingNewLine]", RowBox[{"wout", "=", RowBox[{"wout", "+", RowBox[{"wpower", " ", RowBox[{"coeff", "[", RowBox[{"[", "k", "]"}], "]"}]}]}]}], ";"}], ")"}], ",", RowBox[{"{", RowBox[{"k", ",", "2", ",", RowBox[{"Length", "[", "coeff", "]"}]}], "}"}]}], "]"}], ";", "\[IndentingNewLine]", "wout"}], ")"}]}], ";", "\[IndentingNewLine]", "\n", RowBox[{"(*", " ", RowBox[{ "compute", " ", "product", " ", "of", " ", "power", " ", "series", " ", "by", " ", "convolution", " ", "formula", " ", "only", " ", "n", " ", "terms"}], " ", "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"multcoef", "[", RowBox[{"data1_", ",", "data2_"}], "]"}], ":=", RowBox[{"(", RowBox[{ RowBox[{"mout", "=", RowBox[{"{", RowBox[{ RowBox[{"data1", "[", RowBox[{"[", "1", "]"}], "]"}], " ", RowBox[{"data2", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}]}], ";", "\[IndentingNewLine]", RowBox[{"Do", "[", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"AppendTo", "[", RowBox[{"mout", ",", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"data1", "[", RowBox[{"[", RowBox[{"1", "+", "j"}], "]"}], "]"}], " ", RowBox[{"data2", "[", RowBox[{"[", RowBox[{"1", "+", "k", "-", "j"}], "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"j", ",", "0", ",", "k"}], "}"}]}], "]"}]}], "]"}], ";"}], ")"}], ",", RowBox[{"{", RowBox[{"k", ",", "1", ",", RowBox[{ RowBox[{"Length", "[", "data1", "]"}], "-", "1"}]}], "}"}]}], "]"}], ";", "\[IndentingNewLine]", "mout"}], ")"}]}], ";", "\n", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"integrate", " ", "power", " ", "series", " ", "term"}], "-", "by", "-", "term"}], "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"intcoef", "[", "data1_", "]"}], ":=", RowBox[{"(", RowBox[{ RowBox[{"out", "=", RowBox[{"{", "0", "}"}]}], ";", "\[IndentingNewLine]", RowBox[{"Do", "[", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"AppendTo", "[", RowBox[{"out", ",", RowBox[{ RowBox[{"data1", "[", RowBox[{"[", "k", "]"}], "]"}], "/", "k"}]}], "]"}], ";"}], ")"}], ",", RowBox[{"{", RowBox[{"k", ",", "1", ",", RowBox[{"Length", "[", "data1", "]"}]}], "}"}]}], "]"}], ";", "\[IndentingNewLine]", "out"}], ")"}]}], ";", "\n", "\[IndentingNewLine]", RowBox[{ RowBox[{"makeSCcoef", "[", RowBox[{"alpha_", ",", "para_", ",", "n_"}], "]"}], ":=", RowBox[{"(", RowBox[{ RowBox[{"scout", "=", RowBox[{"makecoeff", "[", RowBox[{ RowBox[{ RowBox[{"alpha", "[", RowBox[{"[", "1", "]"}], "]"}], "-", "1"}], ",", RowBox[{"para", "[", RowBox[{"[", "1", "]"}], "]"}], ",", "n"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"Do", "[", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"next", "=", RowBox[{"makecoeff", "[", RowBox[{ RowBox[{ RowBox[{"alpha", "[", RowBox[{"[", "k", "]"}], "]"}], "-", "1"}], ",", RowBox[{"para", "[", RowBox[{"[", "k", "]"}], "]"}], ",", "n"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"scout", "=", RowBox[{"multcoef", "[", RowBox[{"scout", ",", "next"}], "]"}]}], ";"}], ")"}], ",", RowBox[{"{", RowBox[{"k", ",", "2", ",", RowBox[{"Length", "[", "alpha", "]"}]}], "}"}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"scout", "=", RowBox[{"intcoef", "[", "scout", "]"}]}], ";", "\[IndentingNewLine]", "scout"}], ")"}]}], ";", "\n", "\[IndentingNewLine]", " ", RowBox[{"(*", " ", RowBox[{ "plt", " ", "only", " ", "image", " ", "of", " ", "unit", " ", "circle"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"plotSCcoef2", "[", RowBox[{"data_", ",", "num_"}], "]"}], ":=", RowBox[{"(", RowBox[{ RowBox[{"plotpolyC", "[", RowBox[{"Table", "[", RowBox[{ RowBox[{"sumcoeff", "[", RowBox[{"data", ",", RowBox[{"Exp", "[", RowBox[{"2", " ", "Pi", " ", "I", " ", RowBox[{"k", "/", "num"}]}], "]"}]}], "]"}], ",", RowBox[{"{", RowBox[{"k", ",", "1", ",", "num"}], "}"}]}], "]"}], "]"}], ";"}], ")"}]}], ";", "\n", RowBox[{"(*", " ", RowBox[{ "plot", " ", "smaller", " ", "circles", " ", "and", " ", "radii", " ", "too"}], " ", "*)"}], " ", "\[IndentingNewLine]", RowBox[{ RowBox[{"plotSCcoef3", "[", RowBox[{"data_", ",", "para_", ",", "num_"}], "]"}], ":=", RowBox[{"(", RowBox[{ RowBox[{"circles", "=", RowBox[{"{", RowBox[{"Table", "[", RowBox[{ RowBox[{"sumcoeff", "[", RowBox[{"data", ",", RowBox[{"Exp", "[", RowBox[{"2", " ", "Pi", " ", "I", " ", RowBox[{"k", "/", "num"}]}], "]"}]}], "]"}], ",", RowBox[{"{", RowBox[{"k", ",", "0", ",", "num"}], "}"}]}], "]"}], "}"}]}], ";", "\[IndentingNewLine]", RowBox[{"num2", "=", "100"}], ";", "\[IndentingNewLine]", RowBox[{"AppendTo", "[", RowBox[{"circles", ",", RowBox[{"Table", "[", RowBox[{ RowBox[{"sumcoeff", "[", RowBox[{"data", ",", RowBox[{ RowBox[{"(", ".5", ")"}], " ", RowBox[{"Exp", "[", RowBox[{"2", " ", "Pi", " ", "I", " ", RowBox[{"k", "/", "num2"}]}], "]"}]}]}], "]"}], ",", RowBox[{"{", RowBox[{"k", ",", "0", ",", "num2"}], "}"}]}], "]"}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"AppendTo", "[", RowBox[{"circles", ",", RowBox[{"Table", "[", RowBox[{ RowBox[{"sumcoeff", "[", RowBox[{"data", ",", RowBox[{ RowBox[{"(", ".75", ")"}], " ", RowBox[{"Exp", "[", RowBox[{"2", " ", "Pi", " ", "I", " ", RowBox[{"k", "/", "num2"}]}], "]"}]}]}], "]"}], ",", RowBox[{"{", RowBox[{"k", ",", "0", ",", "num2"}], "}"}]}], "]"}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"AppendTo", "[", RowBox[{"circles", ",", RowBox[{"Table", "[", RowBox[{ RowBox[{"sumcoeff", "[", RowBox[{"data", ",", RowBox[{ RowBox[{"(", ".85", ")"}], " ", RowBox[{"Exp", "[", RowBox[{"2", " ", "Pi", " ", "I", " ", RowBox[{"k", "/", "num2"}]}], "]"}]}]}], "]"}], ",", RowBox[{"{", RowBox[{"k", ",", "0", ",", "num2"}], "}"}]}], "]"}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"num2", "=", "200"}], ";", "\[IndentingNewLine]", RowBox[{"AppendTo", "[", RowBox[{"circles", ",", RowBox[{"Table", "[", RowBox[{ RowBox[{"sumcoeff", "[", RowBox[{"data", ",", RowBox[{ RowBox[{"(", ".95", ")"}], " ", RowBox[{"Exp", "[", RowBox[{"2", " ", "Pi", " ", "I", " ", RowBox[{"k", "/", "num2"}]}], "]"}]}]}], "]"}], ",", RowBox[{"{", RowBox[{"k", ",", "0", ",", "num2"}], "}"}]}], "]"}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"radii", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{".98", "^", "k"}], ",", RowBox[{"{", RowBox[{"k", ",", "1", ",", "50"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"PrependTo", "[", RowBox[{"radii", ",", "1"}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"AppendTo", "[", RowBox[{"radii", ",", "0"}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"lines", "=", RowBox[{"{", "}"}]}], ";", "\[IndentingNewLine]", RowBox[{"Do", "[", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"AppendTo", "[", RowBox[{"lines", ",", RowBox[{"Table", "[", RowBox[{ RowBox[{"sumcoeff", "[", RowBox[{"data", ",", RowBox[{ RowBox[{"radii", "[", RowBox[{"[", "j", "]"}], "]"}], " ", RowBox[{"para", "[", RowBox[{"[", "k", "]"}], "]"}]}]}], "]"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", RowBox[{"Length", "[", "radii", "]"}]}], "}"}]}], "]"}]}], "]"}], ";"}], ")"}], ",", RowBox[{"{", RowBox[{"k", ",", "1", ",", RowBox[{"Length", "[", "para", "]"}]}], "}"}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"plotpolylistC2", "[", RowBox[{"Join", "[", RowBox[{"circles", ",", "lines"}], "]"}], "]"}]}], ")"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", "\n", RowBox[{ RowBox[{"plotSC", "[", RowBox[{"alpha_", ",", "para_", ",", "num_"}], "]"}], ":=", RowBox[{"(", RowBox[{ RowBox[{"scout", "=", RowBox[{"makeSCcoef", "[", RowBox[{"alpha", ",", "para", ",", "num"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"plotSCcoef3", "[", RowBox[{"scout", ",", "para", ",", RowBox[{"10", " ", "num"}]}], "]"}]}], ")"}]}]}]}]}]}], "Input", CellChangeTimes->{{3.46409685160068*^9, 3.464096944064831*^9}, { 3.464097330280392*^9, 3.464097331384153*^9}, {3.464097372714021*^9, 3.464097421999434*^9}, {3.464097464928156*^9, 3.464097468800885*^9}, { 3.464097501968353*^9, 3.464097532447918*^9}, 3.464097565382399*^9, { 3.464097601079987*^9, 3.464097623227595*^9}, 3.464097692005395*^9, { 3.464097782564638*^9, 3.464097788148603*^9}}], Cell[BoxData[ RowBox[{"Null", " ", "\[AliasDelimiter]"}]], "Output", CellChangeTimes->{ 3.464097034530959*^9, 3.464097334977811*^9, {3.464097374236617*^9, 3.464097426536232*^9}, 3.464097503614805*^9, 3.464097568078907*^9, { 3.464097612107178*^9, 3.464097626281115*^9}, 3.464097696955112*^9, 3.464097790245797*^9}] }, Open ]], Cell[BoxData[""], "Input", CellChangeTimes->{{3.464097616973291*^9, 3.464097619604325*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"alpha1", "=", RowBox[{"{", RowBox[{".5", ",", ".5", ",", ".5", ",", ".5"}], "}"}]}], ";"}], "\n", RowBox[{ RowBox[{"para1", "=", RowBox[{"{", RowBox[{"1", ",", "I", ",", RowBox[{"-", "1"}], ",", RowBox[{"-", "I"}]}], "}"}]}], ";"}], "\n", RowBox[{ RowBox[{"alpha2", "=", RowBox[{"{", RowBox[{ ".5", ",", ".5", ",", "1.5", ",", ".5", ",", ".5", ",", ".5", ",", "1.5", ",", ".5"}], "}"}]}], ";"}], "\n", RowBox[{ RowBox[{"para2", "=", RowBox[{"N", "[", RowBox[{"Table", "[", RowBox[{ RowBox[{"Exp", "[", RowBox[{"2", " ", "Pi", " ", "I", " ", RowBox[{"k", "/", "8"}]}], "]"}], ",", RowBox[{"{", RowBox[{"k", ",", "1", ",", "8"}], "}"}]}], "]"}], "]"}]}], ";"}], "\n"}], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"plotSC", "[", RowBox[{"alpha1", ",", "para1", ",", "10"}], "]"}], "\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]"}]], "Input", CellChangeTimes->{{3.464097014183497*^9, 3.464097027394229*^9}, { 3.464097188812593*^9, 3.46409724541642*^9}, {3.464097339081512*^9, 3.464097354595509*^9}, {3.464097657301621*^9, 3.464097678720731*^9}, { 3.464097798431574*^9, 3.464097812485444*^9}}], Cell[BoxData[ GraphicsBox[{ {GrayLevel[0.1], Thickness[0.005], LineBox[CompressedData[" 1:eJxdlV9Ik2EUxkeQdFOEBl14UQu6MI1ibkQ4OVRu9IcI7CK6EGzMoIsxyqiw AoVuuhFmkayRREmgzYl4oReuJ1ALKsJlRVRYLVdWWGpQyYJK3nO6eAbjhe/9 3vec8/zOeT5vJF7ftMzj8TT9+y+tobql34KE1zaczIW2jeRaC9cTKxck7ouk GwezUvKuczLbPC+Zy+2HG6KPpAqPByY65uTW0Wrvs4lJ8d87NOTf8E22jL38 jlUvZN057/nYplmJL+Zlr+eVBJ48/fnrzBe5Hdjf+aD4WjoSU8OVFZ8ldXFj Knd2Ssqu7Lvwp/qTjG7+sb1t8Y1kHh7w7z49I7ui2f5w4a38vhuIHhv7KFdP zbW3leQlsiK2Jzn/QWq6PUO+g+9lvKdusGymIDeK+fKR3un/qz239+yc3WP3 7tQ4FheaxxrNK6l5XtK8M1qHX+uyOq1un+pguvSoTpWqm+louraozqb71/Ln d47M1gSNywm3D9sfcOexXM/3uvthXKpcfAQ0fszlh/WaX9rlD8tfOcC4jLv6 Uar1h5w+6FN9kk4/mH5Bpy9M35tOf4yq/rqiW7k0uvdQq1yK7h6klEvaxUFY uax2eeC+ckm4PHFNuWgd6FcuygHHlYtywFblohzQp1xURwwrF9UZrcpFOUDn I0j7oPOg+0HxQfmB8gfVB6ofpA9IP5C+IP1BfED8QHxB/EH9AeofUH+B+g/U n6D+BfU3qP/RVTpdsaOr3uYDND+g+QLNH2g+QfMLmm/Q/IP8AeQfIH8B+Q/I n0D+BfI3kP+B/BHknyB/BfkvyJ9B/g3yd5D/g74P4O/HXxx2g10= "]]}, {GrayLevel[0.1], Thickness[0.005], LineBox[CompressedData[" 1:eJxdlU9IVFEUxqeCysVEMKAiQQRlQSBoA1OMeRrIoEVBf4Toj5llBTlYoWXq ysCRcrBcJJRmNZLRQm0nFHwuomyEYtKhUvNZEja1aZOLCMy657T4BoYL7757 7jnf75zvraus2V+11OfzVS3+/66lOxd/eTOyK+dYbao09vRuayQ1tnpG2oPP /Ac778nk9+KuUxlPQllD99e2Dkpj0ck1Y22exKItw9l9Q3IpGZi97ffkWndJ +a89w5Jf1H/xVXRajndt7dhU8FyiU96j7oGPksgrnA38eCm+s4W/y9JTcnV3 Mmv6xqhI/t7TKz9Nyhd5vK1n1Rvp3dI0HxufkCvBn6FAMCV1HRXxQ08+SEN6 ofLw67eyonP90f7691KTWRbyysYl90VuQd3mdzIX2Rc/kkz/X+25vWfnLI7F rdd77N7PmscOzatZ81yieT/UOqq1rgqt0+qOqw6mS5vq1KC6mY6m6x3V2XQf OdBcvGEkETYuN90+bD/szmNCz7e7+LD4MXc/Luv95S4/bNT8el3+sPyVA4zL nKsfJVp/o9MHCePi9IPpd97pC9P3q9MfOap/tluRUS7L3Xu4oFxqXRw0KpcH 7h40KZftLg98Uy4LZ/7liRblonWgT7koB5xQLsoB15WL6oRbykV1RES5qM7o US7KATofYdoHnQfFB90Pyg+UP6g+UP0gfUD6gfQF6Q/iA+IH4gviD+oPUP+A +gvUf6D+BPUvqL9B/Y9qb3Def27A5gM0P6D5As0faD5B8wuab9D8g/wB5B8g fwH5D8ifQP4F8jeQ/4H8EeSfIH8F+S/In0H+DfJ3kP+Dvg/g78cfqj4/2w== "]]}, {GrayLevel[0.1], Thickness[0.005], LineBox[CompressedData[" 1:eJxdlV9Ik2EUxocwRruIEiyLSWmgdSOZBKXBofC2i8xAyi5WbKVohghKF9VN bUEXkYKiBeUkIrwTo/2BJ8K7wIqIgY3l/LM+IlsYoRFExnsOwTMYL3zv9573 nOd3zvNVX+htjZT5fL7I5v/fGt38BXKe3K1PeTX9pXTVxKnx6CtP/Al/ui85 JS9/7S5mBjyZvicn/B0Z+RarKR4t86T/a3qukJ2V8L6K0POuz/JgcXngTPlr 6QnPTH1JFeVTdL062/lWmtY+dj/5vSKxZ4UXd7a9l+3hyUd9jStSd70uuP/W Bzl09fDPH5FlOZs73rHxLiuJkeTYruElCYTma1O+eZks1reUMotyevrp+FBl Tt7kBkvJfEEuL91sOlmVlx2h+5Ub6wtye0zavgf+r/bc3rNzFsfi+vUeu7dd 82jUvGo1z3LNO651HNG6FrROq/ux6mC6DKpOq6qb6Wi67lSdTff4TG97+tKf ZuMScPuw/aQ7Dzt/zcWHxR919+Oi3p93+eGK5qccYPkrBxiXc65+NGj9W5w+ MH1anX4w/bqcvjB9Y05/mP4VbkVcucy599CtXBIuDtqUy4S7B0HlctDlgfPK RfPEAeWidcC4KAcUlItywEPlojrhhnJRHTGrXFRn7FEuygE6H820DzoPig+6 H5QfKH9QfaD6QfqA9APpC9IfxAfED8QXxB/UH6D+AfUXqP9A/QnqX1B/g/of wbXhvaXVrceMC80PaL5A8weaT9D8guYbNP8gfwD5B8hfQP4D8ieQf4H8DeR/ IH8E+SfIX0H+C/JnkH+D/B3k/6DvA/j78Rc2XXvb "]]}, {GrayLevel[0.1], Thickness[0.005], LineBox[CompressedData[" 1:eJxdlUtIlGEUhqdpTImiRQOTVAtdaSaEZBcQDoFBi6KarQWRDGEtjESwJDKD oMQuM9DCJkMpKIgESYlIXgyC0MJLeCvF2zRT09RPI7UxuvGd0+JdDD/M93/n 8rznvH/B8dpwxO/z+SJ/f/+ejaXF2b4vGfF+LleuKFr3PJasb071ZwR1J7/t udAtpduujAYqM/KwZijeHnohBzpm77dd+iyTO5ZaQuUD0nN5TfborbSEfNV3 ok+G5XFbUWF57JMczMs/FZ54K5uWC5vCrR/l9hH/SG90XKJeTuLpjZTkpQvC p/1Tcu1e/facjqScS/TkBve+l4mqH5HOZx9koGV67cu6Gcnd3DU0OpaQml07 iwfjs+K9CzScSS9KU9mxqt3zcxKovlpbsrQgv7yRsvXe/P+n/W/v2T2LY3Ff aZ5JzXtW67iudVmdMa27XfvYqH1Zn9Z3vnIwLjPKab9ye6Acjeugcjbu0/u+ Ly70bagwXfrdOez8kbuPrXp/zsWHxQ+6/OjV/IdcfejS+uKuflj92h9Ml0bX P0yX144PxpXPCccPxu+i4wvjq/yxUvn73RO/VZev7j00qy6rXByYLmMuD96o Lq2uDpxXXW66OrFaddE+cFd1UR1wWHVRHWDzqZyQUl22OI7oVl2UM4ZVF9UB uh8VdA66D4oPyg+qD1Q/qD9Q/yA+IH4gviD+IH1A+oH0BekPmg/Q/IDmCzR/ oPkEzS9ovkHzj2RnQzCbLbH9AO0PaL9A+wfaT9D+gvYbtP8gfwD5B8hfQP4D 8ieQf4H8DeR/IH8E+SfIX0H+C/JnkH+D/B3k/6DvA/j78QeZgGre "]]}, {GrayLevel[0.1], Thickness[0.005], LineBox[CompressedData[" 1:eJxdl3tM1XUYxlmbqW2moNPWdJM/moqyvIVOkNdOeKFhMl2kY6yrZ42VIQVT S4GpKeAtRWsBXhA8eVtbuQWKPkkwzUAlhbglHEU5KsLRcFPnnLnv+/rH8wc7 2zn8vt/38nze9/mFf/TFwqUvhISE/Pn/37PPwC5v/2OtQXmr9MxXxSsjKsuy ot4oqg/K0dcqPhk246CcHB/uq/41KN+1Z5V886Rcfho9d1fG2qBMnbw8Ijyz SlJ2f11Y4QmKt3RKyfpXz0pNdntnRqBXTr25ODJ3cZ00ej4dGhnVK7/5VswZ FFov2V8+SkvM6JG1PVMHn592Se6MifP7jtyR2HNzIweObpDNl0fe3HO1Wwrq GtYtq22U8qTMvq6R3VI/xfd2jqdJduw9nRb+wW2p+P3HbS/nNUvTwzXvlv18 SxampvfPPtIiwxPf6/p40C3xXvG0vXKsVQoOtXzWsfKmPKj2D5lwuE1qvROX zOsLSOyBsGBy4b+SUzZ9/v5VAclNSom7u+mKzP8l1DvipYD0jG1uSvq2XcIS x5yJ3tMlLx6IXuKf3SGvJ/wQlx7VJXlnP1yQe6FDcvIXFc+qvSH7B6Oy532/ 9J3IbN2XckOiP787qvW+X5aHJE+KD1yXwprdvQlbr8rxw+/cu5h6Xcb7tjSc GndNLo5a9Li+s1NaHkesWlF1TYblHQ34kzqff9r39n/2nJ1j59o9MXrvfxqH xWVxWtyWRz/Ny/Ls1rytDhu0LlYnq9tfWker606ts9V9hPbB+vKP9um49s36 +Lf21fr8vfY9X3XgUV2YTtapblarjk6orhpUZ6dVd9WqQ9Ol6TRWdetTHZeo ritV5zWqe+NgXvbkIfdTJ8UYJ4fc77DfC9zzME6mu/Nh56e6+5Gs91e7+PCH xnfSxY/LGv8alx+Mk5kuf/Rq/jtcfbBF66NcoELrV+7qC6tvoqs/rP5LXX9g /dH+wTiZ6foL4yTf9R9Z2n/lAqaPAU4/MP1sdPrCRNVXidMfTH8znD5hnBQ5 /cL0G+H0DdO36h/nVf/KBYyPOvc92pQTfQ7GSZo7F8XKid4L41PjQqlyolxA 85BQlxcGKicJLm8ElROtCzYrJ+dc3SDKidYVxonOJxgnja4vME4KXN+gfTQu cEk50b5ju3KiXECUE9UNVEeiukKVcqK6g3GiukS6cqLzHdOUE+UCPuVEdY+9 yolyAd0fMfQ76HnQ+aD7QfGB4gflB8ofVB9Q/UD1BdUf1B9Q/0D9BfUfpA+Q fkD6AukPpE+QfkH6BukfxAeIHxBfIP5AfIL4BfEN4h80H0DzAzRfQPMHNJ9A 8ws030DzDzQfQfMTNF9B8xc0n0HzGzTfQfMftB8wJ2y1J/7ebNsfoP0C2j+g /QTaX6D9Btp/oP0I2p+g/Qrav6D9DNrfoP0O2v8gfwDyDyB/AfIfIH8C8i8g fwPyPyB/BPJPIH8F8l8gfwbybyB/B/J/IH8I8o8gfwnynyB/CvKvIH8L8r8g fwzyzyB/DfLfIH8O8u8gfw/y/6D3A9D7A+j9Avz+8RQzQufr "]]}, {GrayLevel[0.1], Thickness[0.005], LineBox[{{1.1416666666666668`, -1.261617073437678*^-18}, { 1.1051315697687563`, -1.010216343342756*^-18}, { 1.0710712528077977`, -8.089118971543642*^-19}, { 1.0391975572011924`, -6.477211160459941*^-19}, { 1.0092628375174142`, -5.18650604160674*^-19}, { 0.9810538474559489, -4.152997988367744*^-19}, { 0.9543865843751494, -3.325435688886887*^-19}, { 0.9291019385367228, -2.662780610993977*^-19}, { 0.9050620182798119, -2.132171915391573*^-19}, { 0.8821470432691757, -1.7072968978422738`*^-19}, { 0.8602527154596549, -1.3670861511401783`*^-19}, { 0.8392879920494216, -1.0946687403938162`*^-19}, { 0.819173196930115, -8.765355791191171*^-20}, { 0.799838417377421, -7.018695182483049*^-20}, { 0.781222141290176, -5.620089273970735*^-20}, { 0.763270097454537, -4.500181675681032*^-20}, { 0.7459342673115741, -3.6034365517875484`*^-20}, { 0.7291720417333603, -2.88538461745403*^-20}, { 0.7129455005241684, -2.3104179221666476`*^-20}, { 0.6972207958933255, -1.8500240636130363`*^-20}, { 0.6819676241061063, -1.4813722673764934`*^-20}, { 0.6671587720019858, -1.1861812166196666`*^-20}, { 0.6527697271534146, -9.49812487817902*^-21}, { 0.6387783421880712, -7.605446363294534*^-21}, { 0.6251645452676375, -6.089919339535958*^-21}, { 0.6119100899519468, -4.876389338703929*^-21}, { 0.5989983387168626, -3.9046778219623*^-21}, { 0.5864140752691981, -3.1265979466226395`*^-21}, { 0.5741433415389678, -2.5035649970506847`*^-21}, { 0.5621732958504646, -2.004682981778305*^-21}, { 0.5504920892977034, -1.6052125118244707`*^-21}, { 0.5390887577922178, -1.2853439828336804`*^-21}, { 0.5279531276250718, -1.0292152235525352`*^-21}, { 0.5170757327011898, -8.241248961674753*^-22}, { 0.5064477418718462, -6.599026412947177*^-22}, { 0.49606089501807493`, -5.284047333272769*^-22}, { 0.4859074467302848, -4.2311023585973615`*^-22}, { 0.47598011659289713`, -3.3879763067603114`*^-22}, { 0.4662720452218969, -2.7128588444204866`*^-22}, { 0.4567767553215944, -2.172271126886349*^-22}, { 0.4474881171278498, -1.739405593626487*^-22}, { 0.43840031769119536`, -1.3927965904862859`*^-22}, { 0.42950783352696037`, -1.1152558952197939`*^-22}, { 0.42080540622257656`, -8.930203594110194*^-23}, { 0.4122880206463108, -7.150693986382541*^-23}, { 0.40395088544809965`, -5.725784854514528*^-23}, { 0.39578941558308206`, -4.584815440658136*^-23}, { 0.3877992166228105, -3.671205460736713*^-23}, { 0.3799760706487839, -2.9396493074557325`*^-23}, { 0.37231592354856985`, -2.353869360689733*^-23}, { 0.3648148735569666, -1.8848169926736487`*^-23}, {0, 0}}]}, {GrayLevel[0.1], Thickness[0.005], LineBox[{{1.5139404881252136`*^-18, 1.1416666666666668`}, { 1.2576162641613903`*^-18, 1.1051315697687563`}, { 1.0440724570101524`*^-18, 1.0710712528077977`}, {8.663020113299154*^-19, 1.0391975572011924`}, {7.184165222170221*^-19, 1.0092628375174142`}, { 5.954741104209619*^-19, 0.9810538474559489}, {4.933323062772279*^-19, 0.9543865843751494}, {4.085226583631275*^-19, 0.9291019385367228}, { 3.38144126401824*^-19, 0.9050620182798119}, {2.797726498750273*^-19, 0.8821470432691757}, {2.3138459606077604`*^-19, 0.8602527154596549}, { 1.9129207921145807`*^-19, 0.8392879920494216}, { 1.5808840421731208`*^-19, 0.819173196930115}, {1.3060212294214754`*^-19, 0.799838417377421}, {1.0785840017250088`*^-19, 0.781222141290176}, { 8.904657019895402*^-20, 0.763270097454537}, {7.349292627081877*^-20, 0.7459342673115741}, {6.063792555803947*^-20, 0.7291720417333603}, { 5.0017113901212065`*^-20, 0.7129455005241684}, {4.124517957902708*^-20, 0.6972207958933255}, {3.400263552010975*^-20, 0.6819676241061063}, { 2.80247066453004*^-20, 0.6671587720019858}, {2.309206500410383*^-20, 0.6527697271534146}, {1.9023111560820464`*^-20, 0.6387783421880712}, { 1.566755122688832*^-20, 0.6251645452676375}, {1.2901048207777528`*^-20, 0.6119100899519468}, {1.0620782973409367`*^-20, 0.5989983387168626}, { 8.741761076786691*^-21, 0.5864140752691981}, {7.193748423355523*^-21, 0.5741433415389678}, {5.918728113901277*^-21, 0.5621732958504646}, { 4.868791232246378*^-21, 0.5504920892977034}, {4.004388428134773*^-21, 0.5390887577922178}, {3.2928812851614886`*^-21, 0.5279531276250718}, { 2.707342629846361*^-21, 0.5170757327011898}, {2.2255634426176285`*^-21, 0.5064477418718462}, {1.8292311390338*^-21, 0.49606089501807493`}, { 1.5032499239890992`*^-21, 0.4859074467302848}, { 1.2351788722079996`*^-21, 0.47598011659289713`}, { 1.0147675147444501`*^-21, 0.4662720452218969}, { 8.3357214799194945`*^-22, 0.4567767553215944}, { 6.8463894203699515`*^-22, 0.4474881171278498}, {5.622423040960882*^-22, 0.43840031769119536`}, {4.616689299979365*^-22, 0.42950783352696037`}, { 3.79039619015595*^-22, 0.42080540622257656`}, {3.11162290734002*^-22, 0.4122880206463108}, {2.554107738089192*^-22, 0.40395088544809965`}, { 2.096248744859162*^-22, 0.39578941558308206`}, { 1.7202801019388936`*^-22, 0.3877992166228105}, {1.411593375824936*^-22, 0.3799760706487839}, {1.158178376790965*^-22, 0.37231592354856985`}, { 9.501626224922272*^-23, 0.3648148735569666}, {0, 0}}]}, {GrayLevel[0.1], Thickness[0.005], LineBox[{{-1.1416666666666668`, 1.261617073437678*^-18}, {-1.1051315697687563`, 1.010216343342756*^-18}, {-1.0710712528077977`, 8.089118971543642*^-19}, {-1.0391975572011924`, 6.477211160459941*^-19}, {-1.0092628375174142`, 5.18650604160674*^-19}, {-0.9810538474559489, 4.152997988367744*^-19}, {-0.9543865843751494, 3.325435688886887*^-19}, {-0.9291019385367228, 2.662780610993977*^-19}, {-0.9050620182798119, 2.132171915391573*^-19}, {-0.8821470432691757, 1.7072968978422738`*^-19}, {-0.8602527154596549, 1.3670861511401783`*^-19}, {-0.8392879920494216, 1.0946687403938162`*^-19}, {-0.819173196930115, 8.765355791191171*^-20}, {-0.799838417377421, 7.018695182483049*^-20}, {-0.781222141290176, 5.620089273970735*^-20}, {-0.763270097454537, 4.500181675681032*^-20}, {-0.7459342673115741, 3.6034365517875484`*^-20}, {-0.7291720417333603, 2.88538461745403*^-20}, {-0.7129455005241684, 2.3104179221666476`*^-20}, {-0.6972207958933255, 1.8500240636130363`*^-20}, {-0.6819676241061063, 1.4813722673764934`*^-20}, {-0.6671587720019858, 1.1861812166196666`*^-20}, {-0.6527697271534146, 9.49812487817902*^-21}, {-0.6387783421880712, 7.605446363294534*^-21}, {-0.6251645452676375, 6.089919339535958*^-21}, {-0.6119100899519468, 4.876389338703929*^-21}, {-0.5989983387168626, 3.9046778219623*^-21}, {-0.5864140752691981, 3.1265979466226395`*^-21}, {-0.5741433415389678, 2.5035649970506847`*^-21}, {-0.5621732958504646, 2.004682981778305*^-21}, {-0.5504920892977034, 1.6052125118244707`*^-21}, {-0.5390887577922178, 1.2853439828336804`*^-21}, {-0.5279531276250718, 1.0292152235525352`*^-21}, {-0.5170757327011898, 8.241248961674753*^-22}, {-0.5064477418718462, 6.599026412947177*^-22}, {-0.49606089501807493`, 5.284047333272769*^-22}, {-0.4859074467302848, 4.2311023585973615`*^-22}, {-0.47598011659289713`, 3.3879763067603114`*^-22}, {-0.4662720452218969, 2.7128588444204866`*^-22}, {-0.4567767553215944, 2.172271126886349*^-22}, {-0.4474881171278498, 1.739405593626487*^-22}, {-0.43840031769119536`, 1.3927965904862859`*^-22}, {-0.42950783352696037`, 1.1152558952197939`*^-22}, {-0.42080540622257656`, 8.930203594110194*^-23}, {-0.4122880206463108, 7.150693986382541*^-23}, {-0.40395088544809965`, 5.725784854514528*^-23}, {-0.39578941558308206`, 4.584815440658136*^-23}, {-0.3877992166228105, 3.671205460736713*^-23}, {-0.3799760706487839, 2.9396493074557325`*^-23}, {-0.37231592354856985`, 2.353869360689733*^-23}, {-0.3648148735569666, 1.8848169926736487`*^-23}, {0, 0}}]}, {GrayLevel[0], Thickness[0.005], LineBox[{{4.037174635000569*^-18, -1.1416666666666668`}, { 3.2780489508469026`*^-18, -1.1051315697687563`}, { 2.661896251318881*^-18, -1.0710712528077977`}, { 2.1617442434219036`*^-18, -1.0391975572011924`}, { 1.75571773053837*^-18, -1.0092628375174142`}, { 1.4260737080945106`*^-18, -0.9810538474559489}, { 1.1584194440546054`*^-18, -0.9543865843751494}, { 9.410787805619228*^-19, -0.9291019385367228}, { 7.645785094801386*^-19, -0.9050620182798119}, { 6.212320294434821*^-19, -0.8821470432691757}, { 5.0480182628881165`*^-19, -0.8602527154596549}, { 4.1022582729022134`*^-19, -0.8392879920494216}, { 3.333955200411355*^-19, -0.819173196930115}, { 2.709760265918085*^-19, -0.799838417377421}, { 2.2026018565191555`*^-19, -0.781222141290176}, { 1.7905020371257464`*^-19, -0.763270097454537}, { 1.4556165730656975`*^-19, -0.7459342673115741}, { 1.1834561790712005`*^-19, -0.7291720417333603}, { 9.622547234454502*^-20, -0.7129455005241684}, { 7.824566085128781*^-20, -0.6972207958933255}, { 6.363008086763962*^-20, -0.6819676241061063}, { 5.1748330977693735`*^-20, -0.6671587720019858}, { 4.2088314760461865`*^-20, -0.6527697271534146}, { 3.423400428740953*^-20, -0.6387783421880712}, { 2.784738990596024*^-20, -0.6251645452676375}, { 2.2653826885185387`*^-20, -0.6119100899519468}, { 1.843013861733397*^-20, -0.5989983387168626}, { 1.499495697003197*^-20, -0.5864140752691981}, { 1.2200878417456891`*^-20, -0.5741433415389678}, { 9.928094077457888*^-21, -0.5621732958504646}, { 8.07921625589532*^-21, -0.5504920892977034}, { 6.575076393802133*^-21, -0.5390887577922178}, { 5.3513117322665586`*^-21, -0.5279531276250718}, { 4.3555924221813116`*^-21, -0.5170757327011898}, { 3.545368725207064*^-21, -0.5064477418718462}, { 2.8860406056883537`*^-21, -0.49606089501807493`}, { 2.3494703957085713`*^-21, -0.4859074467302848}, { 1.912774133560062*^-21, -0.47598011659289713`}, { 1.5573392836285472`*^-21, -0.4662720452218969}, { 1.2680263733692193`*^-21, -0.4567767553215944}, { 1.0325200607622926`*^-21, -0.4474881171278498}, { 8.408016221933454*^-22, -0.43840031769119536`}, { 6.847201090418953*^-22, -0.42950783352696037`}, { 5.576436908977989*^-22, -0.42080540622257656`}, { 4.5417617046165285`*^-22, -0.4122880206463108}, { 3.699264708992097*^-22, -0.40395088544809965`}, { 3.013211832990789*^-22, -0.39578941558308206`}, { 2.454521194086236*^-22, -0.3877992166228105}, { 1.9995232373160826`*^-22, -0.3799760706487839}, { 1.6289522489289115`*^-22, -0.37231592354856985`}, { 1.327126021026957*^-22, -0.3648148735569666}, {0, 0}}]}}, AspectRatio->Automatic, Axes->None, PlotRange->All]], "Output", CellChangeTimes->{ 3.464097701709428*^9, {3.464097795492109*^9, 3.464097814508579*^9}}] }, Open ]] }, WindowSize->{640, 750}, WindowMargins->{{150, Automatic}, {Automatic, 52}}, FrontEndVersion->"6.0 for Linux x86 (32-bit) (June 19, 2007)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[590, 23, 16310, 430, 1925, "Input"], Cell[16903, 455, 328, 6, 30, "Output"] }, Open ]], Cell[17246, 464, 92, 1, 32, "Input"], Cell[17341, 467, 829, 27, 121, "Input"], Cell[CellGroupData[{ Cell[18195, 498, 446, 8, 99, "Input"], Cell[18644, 508, 16847, 294, 374, "Output"] }, Open ]] } ] *) (* End of internal cache information *)