(* Content-type: application/mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 7.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 145, 7] NotebookDataLength[ 32776, 831] NotebookOptionsPosition[ 30802, 766] NotebookOutlinePosition[ 31176, 783] CellTagsIndexPosition[ 31133, 780] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{ RowBox[{ "Notebook", " ", "to", " ", "demonstrate", " ", "the", " ", "use", " ", "of", " ", "the", " ", "explicit", " ", "formula", " ", "to", " ", "show", " ", "the", " ", "nonexistence", "\n", "of", " ", "elliptic", " ", "curves", " ", "with", " ", "very", " ", "small", " ", RowBox[{"conductor", ".", " ", "David"}], " ", "Farmer"}], " ", "&"}], " ", "Sally", " ", "Koutsoliotas"}], " ", "*)"}]], "Input", CellChangeTimes->{{3.507318121346743*^9, 3.50731817078178*^9}, { 3.507321657189581*^9, 3.50732167169444*^9}}], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{ "The", " ", "version", " ", "of", " ", "the", " ", "explicit", " ", "formula", " ", "used", " ", "here", " ", "is", " ", "from", " ", "Seection", " ", "3", " ", "of", " ", "Andy", " ", RowBox[{"Booker", "'"}], "s", " ", "paper", "\[IndentingNewLine]", "\"\\""}], ",", "\[IndentingNewLine]", RowBox[{"Experimental", " ", "Mathematics"}], ",", RowBox[{ RowBox[{"Vol", ".", "15"}], " ", RowBox[{"(", "2006", ")"}]}], ",", RowBox[{"No", ".4"}], ",", " ", RowBox[{ RowBox[{"p", ".", " ", "385"}], "-", "407"}]}], " ", "*)"}]], "Input", CellChangeTimes->{{3.54592200502666*^9, 3.545922015245417*^9}, { 3.545922081702731*^9, 3.545922118003262*^9}, {3.5459221565623827`*^9, 3.5459222024091578`*^9}}], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{"define", " ", "the", " ", "digamma", " ", "function"}], " ", "*)"}]], "Input", CellChangeTimes->{{3.507320659654529*^9, 3.507320696766552*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"DiGamR", "[", "xx_", "]"}], ":=", RowBox[{ RowBox[{"DiGamR", "[", "xx", "]"}], "=", RowBox[{ RowBox[{ RowBox[{"N", "[", RowBox[{ RowBox[{"PolyGamma", "[", RowBox[{"xx", "/", "2"}], "]"}], ",", "20"}], "]"}], "/", "2"}], "-", RowBox[{ RowBox[{"Log", "[", "Pi", "]"}], "/", "2", " ", RowBox[{"(*", " ", RowBox[{ RowBox[{"$", "\\", RowBox[{ RowBox[{"Gamma_R", "'"}], "/"}], "\\", "Gamma_R"}], RowBox[{"(", "xx", ")"}]}], " ", "*)"}]}]}]}]}]], "Input", CellChangeTimes->{{3.507318186134997*^9, 3.507318213989978*^9}}], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ "now", " ", "we", " ", "set", " ", "up", " ", "the", " ", "terms", " ", "on", " ", "the", " ", "right", " ", "side", " ", "of", " ", RowBox[{"(", RowBox[{"3", "-", "1"}], ")"}], " ", "in", " ", RowBox[{"Booker", "'"}], "s", " ", "paper"}], " ", "*)"}]], "Input", CellChangeTimes->{{3.507318229523122*^9, 3.50731826861794*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"digamterm", "[", RowBox[{"Aa_", ",", "mulis_", ",", "eps_", ",", "ILim_"}], "]"}], ":=", RowBox[{ RowBox[{"(", RowBox[{"1", "/", "Pi"}], ")"}], " ", RowBox[{"Re", "[", RowBox[{"Sum", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{"eps", " ", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"DiGamR", "[", RowBox[{ RowBox[{"1", "/", "2"}], "+", RowBox[{"mulis", "[", RowBox[{"[", "j", "]"}], "]"}], "+", RowBox[{"I", " ", "y"}]}], "]"}], RowBox[{"h", "[", RowBox[{"Aa", ",", "y"}], "]"}]}], ",", RowBox[{"{", RowBox[{"y", ",", RowBox[{ RowBox[{"-", "ILim"}], "+", RowBox[{"eps", "/", "2"}]}], ",", "ILim", ",", "eps"}], "}"}]}], "]"}]}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", RowBox[{"Length", "[", "mulis", "]"}]}], "}"}]}], "]"}], "]"}]}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{"the", " ", "spectral", " ", "term"}], ",", " ", RowBox[{ "where", " ", "we", " ", "evaluate", " ", "the", " ", "integral", " ", "as", " ", "a", " ", "Riemann", " ", "sum", " ", "of", " ", "stepsize", " ", "eps"}], ",", " ", RowBox[{ RowBox[{ RowBox[{"summed", " ", "from"}], " ", "-", RowBox[{"ILim", " ", "to", " ", RowBox[{"ILim", ".", " ", "Try"}], " ", "eps"}]}], "=", RowBox[{ RowBox[{ RowBox[{"1", "/", "5"}], " ", "and", " ", "ILim"}], " ", "=", RowBox[{ "500", " ", "to", " ", "get", " ", "a", " ", "few", " ", "decimal", " ", "places", " ", "of", " ", "accuracy"}]}]}]}], " ", "*)"}]}]}], "Input",\ CellChangeTimes->{{3.507235404139433*^9, 3.507235404542239*^9}, { 3.507235435717652*^9, 3.507235441257519*^9}, 3.507235581714596*^9, 3.507318222978042*^9, {3.50731827444222*^9, 3.507318330152934*^9}, { 3.507320756069751*^9, 3.507320846829998*^9}}], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{"Note", ":", " ", RowBox[{ RowBox[{"for", " ", "an", " ", "elliptic", " ", "curve", " ", "L"}], "-", "function"}]}], ",", " ", RowBox[{ "the", " ", "list", " ", "of", " ", "spectral", " ", "parameters", " ", "is", " ", RowBox[{"{", RowBox[{ RowBox[{"1", "/", "2"}], ",", " ", RowBox[{"3", "/", "2"}]}], "}"}]}]}], " ", "*)"}]], "Input", CellChangeTimes->{{3.507321013397624*^9, 3.50732107038182*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"coeffterm", "[", RowBox[{"Aa_", ",", "Sumlim_"}], "]"}], ":=", RowBox[{ RowBox[{"-", "2"}], " ", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"c", "[", "n", "]"}], RowBox[{ RowBox[{"g", "[", RowBox[{"Aa", ",", RowBox[{"Log", "[", "n", "]"}]}], "]"}], "/", RowBox[{"Sqrt", "[", "n", "]"}]}]}], ",", RowBox[{"{", RowBox[{"n", ",", "2", ",", "Sumlim"}], "}"}]}], "]"}]}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{"the", " ", "coefficient", " ", "sum"}], ",", " ", RowBox[{"where", " ", "the", " ", "upper", " ", "limit", " ", "is", " ", RowBox[{"Sumlim", "."}]}]}], " ", "*)"}]}]}], "Input", CellChangeTimes->{{3.507235516452462*^9, 3.507235596867007*^9}, { 3.507236015313023*^9, 3.507236082568724*^9}, 3.507236129344965*^9, { 3.507320893026437*^9, 3.50732093018009*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"rhsEC", "[", RowBox[{"Aa_", ",", "lev_"}], "]"}], ":=", RowBox[{ RowBox[{ RowBox[{"g", "[", RowBox[{"Aa", ",", "0"}], "]"}], RowBox[{"Log", "[", "lev", "]"}]}], "+", RowBox[{"digamterm", "[", RowBox[{"Aa", ",", RowBox[{"{", RowBox[{ RowBox[{"1", "/", "2"}], ",", RowBox[{"3", "/", "2"}]}], "}"}], ",", RowBox[{"1", "/", "15"}], ",", "700"}], "]"}], "+", "\[IndentingNewLine]", RowBox[{"coeffterm", "[", RowBox[{"Aa", ",", "1000"}], "]"}], " ", RowBox[{"(*", " ", RowBox[{ "need", " ", "to", " ", "check", " ", "how", " ", "many", " ", "significant", " ", "figures", " ", "this", " ", "gives"}], "*)"}]}]}]], "Input", CellChangeTimes->{{3.5072368400782*^9, 3.50723694459916*^9}, { 3.507237047125886*^9, 3.507237049846342*^9}, {3.507237245040582*^9, 3.507237245637454*^9}, 3.507333462445469*^9, {3.507392443621953*^9, 3.507392444177632*^9}, 3.507439484296457*^9, 3.507440181390174*^9, { 3.50750416407072*^9, 3.507504175588226*^9}, {3.508018215820088*^9, 3.508018218072677*^9}, {3.508018321914638*^9, 3.508018340267994*^9}, { 3.5462552181670237`*^9, 3.54625521936535*^9}}], Cell[BoxData[ RowBox[{" ", RowBox[{"(*", " ", RowBox[{ RowBox[{ "above", " ", "is", " ", "the", " ", "RHS", " ", "of", " ", "the", " ", "explicit", " ", "formula", " ", "for", " ", "the", " ", "specific", " ", "case", " ", "of", " ", "the", " ", "L"}], "-", RowBox[{ "function", "\[IndentingNewLine]", "of", " ", "an", " ", "elliptic", " ", RowBox[{"curve", "/", "Q"}], " ", "with", " ", RowBox[{"conductor", " ", "'"}], RowBox[{ RowBox[{"lev", "'"}], "."}]}]}], "*)"}]}]], "Input", CellChangeTimes->{{3.5462552340420847`*^9, 3.5462552573946953`*^9}, { 3.546255358940444*^9, 3.546255382189128*^9}}], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{ RowBox[{ "we", " ", "will", " ", "test", " ", "with", " ", "functions", " ", "of", " ", "the", " ", "form", " ", RowBox[{"g", "[", RowBox[{"A", ",", "x"}], "]"}]}], " ", "=", " ", RowBox[{"g", "[", RowBox[{"A", " ", "x"}], "]"}]}], ",", " ", RowBox[{ "where", " ", "we", " ", "will", " ", "decide", " ", "later", "\[IndentingNewLine]", "which", " ", "function", " ", RowBox[{"g", "[", "]"}], " ", "to", " ", RowBox[{"use", "."}]}]}], " ", "*)"}]], "Input", CellChangeTimes->{{3.5073183581236*^9, 3.507318422414816*^9}, { 3.50731852382278*^9, 3.507318539340487*^9}, {3.507320880702779*^9, 3.507320881473548*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"g", "[", RowBox[{"Aa_", ",", "t_"}], "]"}], ":=", RowBox[{"g", "[", RowBox[{"Aa", " ", "t"}], "]"}]}]], "Input", CellChangeTimes->{{3.507235664209149*^9, 3.507235682054463*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"h", "[", RowBox[{"Aa_", ",", "t_"}], "]"}], ":=", RowBox[{ RowBox[{"h", "[", RowBox[{"t", "/", "Aa"}], "]"}], "/", "Aa"}]}]], "Input", CellChangeTimes->{{3.507235708864558*^9, 3.507235718531132*^9}, { 3.507235774875444*^9, 3.507235797163828*^9}}], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{"Here", " ", "are", " ", "four", " ", RowBox[{"possibilities", ":", " ", "g0"}]}], ",", " ", "g1", ",", " ", "g2", ",", " ", RowBox[{"and", " ", RowBox[{"g3", ".", "\[IndentingNewLine]", "These"}], " ", "have", " ", "corresponding", " ", "h0"}], ",", " ", "h1", ",", " ", "h2", ",", " ", RowBox[{"and", " ", "h3"}], ",", " ", RowBox[{"in", " ", "the", " ", "notation", " ", "of", " ", RowBox[{"Booker", "'"}], "s", " ", RowBox[{"paper", "."}]}]}], " ", "*)"}]], "Input", CellChangeTimes->{{3.507318542035799*^9, 3.507318544954497*^9}, { 3.546252006839641*^9, 3.5462520715564337`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"g0", "[", "t_", "]"}], ":=", RowBox[{"With", "[", RowBox[{ RowBox[{"{", RowBox[{"Nt", "=", RowBox[{"Abs", "[", RowBox[{"N", "[", RowBox[{"t", ",", "20"}], "]"}], "]"}]}], "}"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{"Nt", ">", "1"}], ",", "0", ",", RowBox[{"1", "-", "Nt"}]}], "]"}]}], "]"}]}]], "Input", CellChangeTimes->{3.507318749018969*^9}], Cell[BoxData[ RowBox[{ RowBox[{"h0", "[", "t_", "]"}], ":=", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"Sin", "[", RowBox[{"t", "/", "2"}], "]"}], "/", RowBox[{"(", RowBox[{"t", "/", "2"}], ")"}]}], ")"}], "^", "2"}]}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"h0", "[", "0", "]"}], "=", "1"}]], "Input", CellChangeTimes->{{3.508018735965552*^9, 3.508018739145742*^9}}], Cell[BoxData["1"], "Output", CellChangeTimes->{3.508018739706375*^9, 3.546252212301498*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"g1", "[", "t_", "]"}], ":=", RowBox[{"With", "[", RowBox[{ RowBox[{"{", RowBox[{"Nt", "=", RowBox[{"Abs", "[", RowBox[{"N", "[", RowBox[{"t", ",", "20"}], "]"}], "]"}]}], "}"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{"Nt", ">", "2"}], ",", "0", ",", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"Nt", ">", "1"}], ",", RowBox[{ RowBox[{"4", "/", "3"}], "-", RowBox[{"2", "Nt"}], "+", RowBox[{"Nt", "^", "2"}], "-", RowBox[{ RowBox[{"1", "/", "6"}], " ", RowBox[{"Nt", "^", "3"}]}]}], ",", RowBox[{ RowBox[{"2", "/", "3"}], "-", RowBox[{"Nt", "^", "2"}], "+", RowBox[{ RowBox[{"1", "/", "2"}], " ", RowBox[{"Nt", "^", "3"}]}]}]}], "]"}]}], "]"}]}], "]"}]}]], "Input",\ CellChangeTimes->{{3.507235656001355*^9, 3.507235659081276*^9}, 3.507235759603612*^9, 3.507318750879644*^9}], Cell[BoxData[ RowBox[{ RowBox[{"h1", "[", "t_", "]"}], ":=", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"Sin", "[", RowBox[{"t", "/", "2"}], "]"}], "/", RowBox[{"(", RowBox[{"t", "/", "2"}], ")"}]}], ")"}], "^", "4"}]}]], "Input", CellChangeTimes->{3.507235766793907*^9}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"h1", "[", "0", "]"}], "=", "1"}]], "Input", CellChangeTimes->{3.508018745945376*^9}], Cell[BoxData["1"], "Output", CellChangeTimes->{3.508018747111911*^9, 3.5462522171305103`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Plot", "[", RowBox[{ RowBox[{"g1", "[", "t", "]"}], ",", RowBox[{"{", RowBox[{"t", ",", RowBox[{"-", "4"}], ",", "4"}], "}"}]}], "]"}]], "Input", CellChangeTimes->{{3.546262684711094*^9, 3.54626269879608*^9}}], Cell[BoxData[ GraphicsBox[{{}, {}, {Hue[0.67, 0.6, 0.6], LineBox[CompressedData[" 1:eJxN2nk0VH/4B/AhIkvJN1Ih2UOyFZJ5kLUiylayZM9SyFIRk1AkOyFFdpIl soVrT/Yt+24mkuxkCb87nd/P/OaPmXmd9/M8997PZ2bOnDNzyvT+dQtyDAZz CL0jPgoD88zu7qFKzP/e7m9wvdjYJDk3X5R3eZVkMX510x+/SJY87DPQ3Euy 3ORiY0wOydeft3wQNybZrc3H3qri4J67u2lnqu7Q75mt8v4RSVa6PX+2raU8 NEyzZ5NA0YV+7wN7Xg9zf4EA9Z6TDo4cix7fv+e4y6JOXf6Ue/a2u/PwtBrF nt33Zbo5/Sbf8x1cG3dPOtmeI8sJ0XE2mD0bPYz1lLQgmV/8mnmnMcml6cVn qXVIHgl/Wf8ASOazEV25fITkBe4f/YSDJJeMxiK4AyR7x15LOUFB8utMamb7 7V2E+FyzrF5WgMNpPHHpn3F6Q1LYssjBrL6pf9anu2mK1aBVenhw+J8nFL+6 Y8dw2ZcUO//ZmpkrGOu0dvTQ46//rH7L/C2Wwu7pQE7ZP2dwd6Zjo8ZnUgh5 RMvV9pnmY0/raTueSCMaFzUbUYYtay6/qBX3L4dp3lqshgIf9fNQojmi991r wo4VhXSV+RHd3uvD1Yl1OrP5bsn9n+ff2vViKZLMbPgdiV6v/J0yjI1iaTln ZEk0i0eA4wSWP+g8WYQB0VKnBt78wH7Zl9D8TZNoVd5Cg19YjUcHoneViB5z cDNcwI7NOZmdkyGaRifzzjLWyXxI2FbkXz5jgKxhKQaUNhN4iI42dHbfxEZd y6nrOU60Q99g7A6Wv44llI6BaK1LsqFkUHrB+7YCJdEFanG0FKCe+4vv4eYO cX2Ydiv2wxiPzvLHeaJbS/s/HgCnNxUVk3iiKzUTBuggyi9U51ob0UjH2NZh 4P+7yeFbixpXreS99h+UOprPlpYQc8ZhMiZmuDrVUrSQTbRHYow+C4zelnzG m0x0W8Kz6uPg1JmgcTuG6P6JjRtsQKFKczwsiOjP0ieoOIBffDh3+yFxfpbO KoYLStOVPcTvoZYTkPj1jhuusueq3DUj+sW14hu84Hjg2Ui3OurKZ9evUggA hddsBs0lYj6vqE4mBFErOi5yUsR5ff4xR4WhdJSfLouTeLw7DIpqouBxkexj LAvqsUS7C1ZiIBvTr+5/kNj/6YlcqDhUXg8IttzcRr06/JjmPHjnmIroLKCu pCoXviUJinQyHZd+EPNkh8VCKfha94vxVCfqBGmXxGgZeMFZm3+oATWOQZeB XRYue8Vp75SjHhMWDMzDQquUetRgBrF/17qHTB66P2Qfi3qKOpdmve6jIkRR Py/1cUPNkZEumqUEehbGBg/sUZvwP0j6pAwD7AxvNW8S+9MILH2qMB7qyEEr Sjx++dmGF+qQOKdWtcGLWtPt3XijBphf4TSdZkXtYMWFodWEaYqupDpq1AvH bX6YacGcmwQfbuwvWu8xHayhDbnddA33elBzDG7fZ9IBR1GCtWEzMZ+C6+M6 sDITmXmhGDX/n8OGL/Rgy/DPmdVg1OleFpeZDKCstK110hc1rltvk/I2eB5N v9/pTuz/vF/o720g79DPy7FCvXDAVHvDCKgVSyVsALWQUo1M3B1oTAj7rn8O daX5zCNBUwjctnFVEUQdUlJUXG4Kh4pOFHMfRe2zya//ywyYBZ7IjM5tIRiG +PfbtpbAxSCnoP0Otea7b0MH7CBIZPyxTThqaqU0m2E7WNf0zse9QF2co8P8 2R6aQ2u5PzqhxpVWzzvcB+f/1Kj2q6Bef/Z68pUjjIjPAOtF1KYsSznUTqCq /fKhmChqxfXOfF8nYIts+Wl0ArVQsKur3wOoZ77eXDS/iWAKG+1dklxARHKZ ohWPms6OK1DEFWL1ImTx/ai1FbbJK13hfnRPNkMtak9vO54fbsBy3CD0bjRq fXLTpmuPwYbNXPeEPOpW/zm6t17QLUsZLHp+E8GpHYm4dBoHWKPUryqCqBOn TG764uC/+GkpZybUYUfEYuSfQhmH/YmWnxsIZv9YfW6jNxzkdh3zDENNQW74 8JQfPFQ8eizyOephq+OLYX4wYV6s9cEDNe8Xnl+Uz6EoZbO613IDwdkmq9sv PAcTPq8UERk0vzWwttDrD58EXthMTK4juKfRWsLLr+Bkh0mJYD9qlxp5Y/cg CHSTpnZpXUf3gytJgDIYLGpnUqhK1hG5Zg+BloMhwGKkPi4UtI5UtiTdYhAN hSehjPpuUmh+Ruvz4ttwmJH8lVp5Bs23/uuynw8H3ZGa1QNcqP+amDTKR4CI oEt4HD3qFzqvT/+MgIna3taqiT9IpU4zc6JCFKiuxynRvUItS2/44lQMfH7n EqHj/QfBGTwKfvM0BjiVNCbfuaG2M7ZOnIiBrdBdnKgpWq+rMxicFgvZgqZl upKoB2ktaWXj4Igxn/j78TVEjrfyUNjreBity+M4fx51a8RIon8SeNzYUJQS WkM4hC2/FNQlAcu43N0LnGuISczoBXnyZND825aHPbiGJET6H2n2SIZKsTkF lR+rCIfk8vdCtxRIiBew0I9CXUnIVPdNA9NHSRmP11aQyp+Btj/Js2Bn/68W j9kVhIOVICp/LQtiI8SWPCfQ3GFYhCsuCzpzqi48a0VzH/um55If4dKP0aaX qSsIzlKRedw5G7hvsM7F6a4gucOXgzh8c4EgFCmGFC4jJp2Uw8Fa+SDgmTI3 /WEZqSSjzh70yYd7bZ8zGd+jPl3x92lxPqw79nBavlxGNLuctEY5CoCm6Oh/ 9CbLiEOcxEvCagGIyMcu3zqwjIjgOvst8grBXSe+YM1gCRExk8kzdSiFytQc Bw6tJaS9T0E3OqcUKNYRocvKS4hc/R0bsflSCIodS34rsoQ4PGRkj7z/BRJG T0UqUi4hmmqulmwuZVB/N9klLHsR4T+KYfIKr4DDnhnnhckWEQZWh8H+4Crw K+RvYPizgMQZPSlxz6+Czbm0m8uzC8htj7cZ6r1VMGGc6lHct4DwFzrJe5+s hk/ySdUKuQuI1In1PtFP1XCd8q2GnvECknw/dbl+rAb4jp1qwbjPIw8fMF0J N6sHofPzQinLs0gjhSDLq5UmsBT28ioZnkV0J3X9EvmaIZ73UGfL11lE86Sm 8cStZmA8etZt7c0sUvlj1my0qhn+/LlXpaI4i2DOx0tWhbZAdcmczkzELyTg pylVsVQb6MvOeZ49P4N8TuBiP6DZAc8Uf7eXPJxCrK4qqpcNdQOTdnbIzaMT iKPc0ITSziB84n1ZqrRvAkklVxmyZhkCzQ0rvOj8OOLgQth+JjYEAe9OSdF8 HUe4btP/fW85BDvTEaOlruOISolUgVvrEPx44i7M2jOGsNA/jjSJG4Z710Kj y0+PIvrz/nxDZ0chfeXybCTZKBKex6mcrjoKEzEUcvf6RxC6+HvxlndGQRfv NsXuP4LQ/7qvWBI2CvDI6DxuehiJ71Thz1sZBYYkwW6FtCFE9a357NitMfi0 VneogXsAOfz+sfMrpnEIykysxm70I2qGmmnDp8bB1sjL5XNLP6I1OnaJU3gc uOulBhNd+5FaNy35F8rjEBX5IfVJQx/C03Zlt8Z1HB5LhMqK2/UiNRXuzGPd 46DgdNv2XX43cq3/CX+v3wSw80qzM73oRg5XtVM+Dp2ArX6mjpe3uxF3vz5f +rgJKJBvO/+IshtR+s2zSpk3AbyHFci19bqQyN0wBemBCaDJ5Ys5sNWB/DcS yBd8ehJSBgmvfvq0InalBSHcXyZhwIM5/q14K7J8WjDncM0kHGRXydWaaEHy BBJPzjVOgqtJemcJtCDnp89GOA9MgjLB5mjAZhMiFcR64vP6JEzPzScIOHxD pm43H58XxQNrCMenEfZvyPpge0edFB40RbVqwloakGF6Da0AwEPxg0+ETYEG JDRMU2XoKh78150Fmgj1iNOtBqP3lngQIN/MtzWoRTpKjxTGRuHBKEmg7iRN LWJldtpeIQ4PYYoGPV3FNYjTmLBZz3s8bPqVrcsw1yCaJpkZ1Vl4aKL1kqXr qEIwnz7PZ1ThAfMxR6PSqwphlk76nVuPBwmNMWNn4SrEKaXcIqkJD3Eh8t5D LyuRO+ePml7+joc+Aeqgm68qkJUHOQNvp/AgZ383t2qmHKE10d3Y/oWH9JzG ztOq5YiUgmWL6gIeHkm8OrpJXoZILhqJv1vHw3Es4/vYhyVIpsR+ZSVqAnjj HtSQ9xQjcnQqxuu0BPhV3U2wES9GymaE7kYeIkCZymsBmblCpLHWdzaQmQCG WqwFg2YFCAt/BeY+FwFml1elXWrzEQdTdkYTXgJ4RLUhB3nykY1iVnGZ0wSI G/Rukv+Rh6hFyqnkCBNAyPO21qByHmL5go5OVZQAXzjO9zqn5SJGTjXPGsUJ UNDUWXetIxsJoYSzgVIE6Gbw6/AO+4iwsmZKNF0gwLKO9PDnG1mIKIV6yPJF AoiNxa+c6MlAorme/96VI4AWzw3MtdfpiEJDbeWkAgEcbPbTeeunIazHnr/K VSRAzood1/RAMsIY+6OQQpUArdIcZ0/EJSFPJQUIL9UIMOfZdUHDMBFpNYDc rcsEoK/1U37K/h5JlN3AaF8lgD6jc7I1Eo+k/SdYGaFOgNdxA19a5d8i/kfS e8o1CNDDK98lURuLLLteOtN2jQBH8tJmYpWjkSDH9qwGTQLckDlITvYtEqGS 4cRmahEgrM75mNWVcISKhq7xwXUCtF8bFGlpCUGkmO5I89wgAHvBsxDLhEBE JXzFuQL1mwZNRCzlOTJcGmMB2gRgGWab28nwRlKruVbSUAtOfpeOMnRHro5b HFhHzUvn8FrYyhF5pC4edFaHALhu/u4TNYZIn6yxrQZqoUPDNpRnFGFre8hf G3XPonc8w7o1/FC4MS6PeuEa39J9J1cg3y30ZUbtHq7caCuBg+QtiaBOdP6O 56M38Wd94RT55mFX1DibLLsugQBoHxTjJ0O9T3dUloo3GCiv000/QM+/hsP/ +VuGMGASzrduRa/XJ4rTQTIoAvSktasYUCvTlel30L4GdlnWAzLo+lB568jb +sdANbmUnjq6fg1/5k5TUsVB5I3uahV0fQPsXzDG+7yDuXU2MwF0/VMXQ56w Pk0ArjL5a6vofllZCVgW7iSCc1WjlSK6n6qcbreeySbDhSj/6gZ0//mHazQ0 PVJANPt73TkVAsxcN5Sc2UgDypYGzsZLBGikzxQsksqAIMtvmYvyBPjQsHbS xy0T4qxppMjR15sdNoSafTUL3jZhkBEZAqhvDP2dEc+Gi8ihqUxpAggXnF4s csqBz3KMjiaSBFiPGJ/lX8+FIMYUqhdiBCihcf9e7ZIPp/9MfW3iI0DCGaau HM4CSJfbDkjhJoCfZk57XFsB+FTJBVqcIoB21GST6+lCMAn8iEQfR/fn1NVq geFi2JCJ421D37+9ij+QowElUHDu/dERKgJUWOHKKSRLgZAVIN69jwAvPxYU j4R8gZzzlqauW3jglWbLCVOsANmP5X2RP/FAd7soy2uxAhaemxuO4/Gw7KmV afcOgbO8qRtHxvBgMlF9ncujEozateclevAgzP9Hfft7FbS1D9MaoJ9n3/KN FfP8ayFAKOj3zXA8XKwTquEbqwVjg+BGiVd4yO3ZkH93vg7m+19Hb/jhIWoj HF5O1sHln8c85R/jwRwaLljIfgUVZZMPLcZ42G0UET22+A3GkwhOi7x4ODdB xv5Uvw2YFDS/n0udBOFSvUWRjDYQY5riZns3CbxhH2vHNtrAcynZZSVyEmjf Z51LUG+H/fJ1jx/4TsJfTf2h3bl2aDPDF3w2m4Th3JzTiGgnPL+/s2TMNgnx Dka1skXdsHIs5+ey7wRwLXzZuFDdD5igezdvi46DR8sln5Xf/eBQuSFBzzMO PZlNB7OPDQCXU5z+J5ZxCLAY5DrlOACtb+OW+nfHYGFgU52KYxAGTxY2vm8Z g4o66cSuJ0OQcqKiW9hqDG69KbpsJz0KYXmqkqd5R0Etqb9FQ38C2KD0vXfG ICipM+5b/TwF6Wm/jcgDuyHozoHFjW9T4O7h2JHu2A19LpjRneEpsPl6hlpH rxts3s2VUu+fhjd+tzwRTrR+vtGJVXcaJo+3542XdEFf6LOJS2vT4MDxXM/n ZyfY9KxUhZ2fgcxGTN7Lax0QbNyLEymaBb+fVJY9T5vBb/Yi/GieBe33j2ZX NJrB81Hi9puJWbiKOZx3kbUZ7oXbP6Y6+BtiN7W71YqaQOPrPucR89+g256T YDXfCAfPilgFMs7BEc+Ei8etv0HQzgv1aft5cG7/OZboWQ9+L+do3z2bB+yX RK8OrXrwZNFuvBEzDxXPkm5K8dTDPbGTKkjtPHxa4YjPba4DvvxbHE30C7De PKSrz14HVFMLxxLNF8AkxmLlTXMNfNVgo9dkXARu8XO0v29UwRHhcWYn1kVg YdzPVCZSBSb0KRwRvItgN96EL6VHv782CUn0XViE5h5+bZWGShBSu2hgYrYI uLpvk6+hEkIuGWQ4FCwCZvFu2pp9BdySilEK1V2Cyq/nytuCSyH1qOG1fJMl aJel1BrXLoXlNY6b322WALf449rl46UQ+Dnd7pjXEiyw1QecSi0BRKw4/H0a 6sHzfcNVxcB9pnc8b30JTB7+/aV9oAh+n2Ly7IxdhsqQ2a8uPfkg5f7esSoZ 9YEVj5nUfPDuPmORm70M7TltFh/c8oH5ufLVoGrUGRvDDMfzATvrdvzyzDKI WHQVmJl+guDCgcIq6RUIcXU1XcbkgeiVhPncvhUw6XrSd2I6CzyShSbjJ1aA I56+h/tNFtRvF/cEza5ApWj6TKFGFtzK7Si3w6yC3MbPe85FH8CbaV8gP/8q VIowP1EJzoSuUQv+BLdVSJCpdp/QSwcXZ8E7wcxrwPFqwThMLRmYaebLdE+t QUJs51gmVTIUx39iYRdag8pbwUXX65Pgb6N0e5Y8Wr/7YEVdKQm8T6lCox2a 804FjqslQlCrORtlzRqMORjLrdglgIg5/6PmljXAkEVXFIsnQMfGr+7wvjXA TQpHUK7FwxGeB4Gcc2j/zsszl5njIdb96RYc+wOV0hVtnGZvIZX/Xf/j+39A rqR3fFcoFlQr7pxTeIzmIe71aY0x8PMGT+gBX9Q+vhIqd2NAyCtLNToWre95 2mCVEQ2fvpcWfa5D82rmMPpzr6HCuzdy4cQ64Er9vvQ/jAATljeLRbzrILc1 4WLFEwFk2UbqXqKo85/w93eGg2I/geKgyjpUvuMOMxMJh8azKw8EndD+feLS S6uh0DvEcN2yATWZRvO3G8Fw+aVFn3Mn6gdbbwI6gqBMutTo2RBqCiPLQs0g SIgys01YQM27f0P5xiuw1iz0GWDZAEy1jOuRuy9hs8agSOPuBuCAqivp+3Ow c8qVNXyA2nh5hufucxjhoKy1fYKaLUliedsPqp9kt78IRfsZD9taCfrBS0my meoS1LK4Ws9IH2D9kMoqSbMJGC9XoazkpxB0cytR6Qjqs/7bS1eeAoZa87Q2 +ybgNu0lLVZwgDffOOcohtrXcblGFQcf2a9e+3ALrQ960G+x8wQ4WhK+l5ij trD/dsH4CYS5rxo03EO9eqxaudoDXPveWRO8UTsz/lwOcgds2OLTkx9QZ4u3 iis8glw5ZWrhz6gXXGq1Ch4C53xs0EUEtfWjTT/+h7D/quKbm12oRa7I8x11 gzbK1wXhW6hnmjfb6FzA5NHFKaqrW4CpHOFU1XUEAfOG73O6qJnBp2XEAZY1 tGu/30Gdu9LmcdcB/Ljt3ie5oeZ/3epPcx+y2t4YQBJqE8nigQ+24FrKd5k3 G7WI+EfrERuAlHwp+hLUcv6PN47YQOfjJqbBVtTTJzi2XlnDBs9Wm+smap9q LrNsC6hmeF5hSPkXMFJu1rO75hC4xfhRkQH12Dl/Ph1z4OgQCGDkRT2dPhZA YwbKHrcUs7WI+VXnmXgTYLD6IRZ5m9gvxTtFbwL9Wk6nPKxQ40rv++8YgR1f wI7aE9Sq+uvLvbchrLOkmJBO7H8haNanD7fLldKa81HrNySLT+sBT3pHZH4F 6ugBPNeOLhQ/+en0tBs1C18540UdGOE/doYNs41ev/D00rYWpP+XcoKCFjVD Z17ARS1w3BGh+cWEGpOuloDTBIpu1aliQWK9RNuUqAYIeD18r62HmqN5lpZB DVbuUoTImBL7u6Iyx1SgXDvEk9Meteb2L8wXZdAUSDdY8CbOOxEvFqgIx5kk Lve+Qi1XgSn3vgSTu4hURTRxPv6IobcCuH7vZQrMJtaz3L4YJwdylWYUD0pQ J2Ty6NgBHPgwv3SzlpjP7D+WJwtxOKp2vgGincX/s7sANYKSATWUO+j8hTer 7yTA6Kx20uIh1Aln9w9YisO6mGPZyeOox948uCIjBkIXPsy5C6PGbL9uoRaB etkGqkwpoqeGC8iFwUSewNGnQJwnzthPIwQRqidvSOgR84VXokb8IHz1op3p HdQ4MyajdF5ouHbTN8SWmF9M8qXigW3d8KJZL2LOfb2QnBOib+W2Hw8g+nx+ 0VsOEDNq+akaQbT1ObXX7GBpQcWakk700oZp2THA3OU+1/WJOM/vfvizoxBr J69BVv7vfBfd7ZhAwsHI6uxXog2qjrj+B60P3HGGHcT+RDbd+MNg7RYd83KQ mH985dB6CMjdP38qIRC9tmEZQweS3vN4pk2iZTmy+vZDpy/dziWKXdSqA4Oq FGDnf/qo00GiL2fv/iKD/a+URRJYiI6WYsDtYhNCzNRaOYlOVW3y3cJeiMCZ /hUiuj7wMOEPtvv1W3cBSaJpKGtjVrD33pRG6MsT7ZFXE7+IpY7v/eh3heir 3JOE39jExJX6Ah2iRTQM2X9iL6YeHpswJrp4rGgMj+3JEN5gsCH672AB9zjW 4eMVRnD+V48Tqh3E0uRZC9p7Eq15/OKpHmxyga/imxdEO/ixanZgscWJht/C iM71T8I1Yvu/IK5/4ojGfUCqa7DOyFAwTxrRlXiBl2XYgzUb6TfyiMZU3GzP x6bXM1c//fIvr5fkzcQqNIoP5tQRLfeBnCsBO9SiuTLc9q9e15o3DOvaYU9P N/DP7HlNT7EM3wN4L+D/mTsXb4vN7EsD67l/ljt6+ypWaahWP2r9n3Gf1hmw C/MLUqZkGLn/+/2w3NFL6dF+kgOW6a+H0JKs5xJnlMZAMvcfAdsKJpIrNlV8 Zk+S/NKjJ2QfD8n6O+ZvjwuQvETmXagqQTLvgbLpZNX/l38VtklU/3/z/d7/ ir9Osi7F87nY2yT77mqthDiQjF+d2vWIJjlx4r+jOtMks5XZKlG+INszZpM8 w2mCfM9/tMqN6ngo9jwoYXyCx4dyzw56DIc5h/bv+VgEXp+CjXrP2QYJacsP D+xZT9ZLtaGcZs95RvacF/7S7nmd5tK+Eh36PV/vkDnPFn1wzzQytFVlb0mu Sh64cjuJZJGHj+68ySb54MmiwGN1JH+zE5tkWiJZjvp0yMGrh0jHd1w/nq1F cs7g1xR1PZLZcyy/BJqS/Fcn+QfNI5KLkk7KUqWSfI9+/mvqB5J53SquK+eR HHHZ0Nq3jGSnxdjwfV0kCxjYsCf1kTxeK52hMEKyZnRfxdOfJFOTp6udmicZ sXXrrlwh2fW7srHJJsn/+/+GPf8PnIkB9A== "]]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesOrigin->{0, 0}, PlotRange->{{-4, 4}, {0., 0.6666662876741037}}, PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Scaled[0.02]}]], "Output", CellChangeTimes->{3.546262699229405*^9}] }, Open ]], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{"Note", " ", "that", " ", RowBox[{"g1", " ", "\\", "in"}], " ", RowBox[{"C", "^", "1"}], RowBox[{"(", "R", ")"}]}], ",", " ", RowBox[{"it", " ", "is", " ", "supported", " ", RowBox[{"on", " ", "[", RowBox[{ RowBox[{"-", "2"}], ",", "2"}], "]"}]}], ",", " ", RowBox[{ RowBox[{"and", " ", "h1", " ", "is", " ", "non"}], "-", "negative"}]}], " ", "*)"}]], "Input", CellChangeTimes->{{3.546252663402267*^9, 3.546252727604809*^9}, { 3.546252852476885*^9, 3.5462528585367813`*^9}}], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{ "The", " ", "formulas", " ", "assume", " ", "that", " ", "the", " ", "functions", " ", RowBox[{"g", "[", "]"}], " ", "and", " ", RowBox[{"h", "[", "]"}], " ", "have", " ", "been", " ", "defined"}], ",", "\[IndentingNewLine]", RowBox[{"so", " ", "we", " ", "will", " ", "use", " ", RowBox[{"g1", "[", "]"}], " ", "and", " ", RowBox[{"h1", "[", "]"}]}]}], " ", "*)"}]], "Input", CellChangeTimes->{{3.5462529895493803`*^9, 3.546253001203352*^9}, { 3.5462532249433327`*^9, 3.546253252675933*^9}, {3.5462532894200373`*^9, 3.5462533278071737`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"g", "[", "t_", "]"}], ":=", RowBox[{"g1", "[", "t", "]"}]}]], "Input", CellChangeTimes->{{3.507235888686879*^9, 3.507235895251235*^9}, 3.507390732127739*^9, 3.507392719449676*^9, 3.507392858064131*^9, 3.50739295032518*^9, 3.507438969986606*^9, 3.507502562112997*^9, 3.507502810575395*^9, 3.507503053233705*^9, {3.507503093812856*^9, 3.507503109802619*^9}, {3.508018699939038*^9, 3.508018701104154*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"h", "[", "t_", "]"}], ":=", RowBox[{"h1", "[", "t", "]"}]}]], "Input", CellChangeTimes->{{3.507235905397222*^9, 3.507235910275342*^9}, 3.507390735153994*^9, 3.507392720832934*^9, 3.507392859583201*^9, 3.507392952151137*^9, 3.507438973720504*^9, 3.507502564829716*^9, 3.507502811810203*^9, 3.507503057137734*^9, 3.50801870504343*^9}], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{"Let", "'"}], "s", " ", "show", " ", "that", " ", "there", " ", "is", " ", "no", " ", "elliptic", " ", "curve", " ", "of", " ", "conductor", " ", "4.", "\[IndentingNewLine]", "Choose", " ", "the", " ", "test", " ", "function", " "}]}]], "Input", CellChangeTimes->{{3.546253712659102*^9, 3.5462537266515093`*^9}, { 3.546255190046164*^9, 3.546255192475301*^9}, {3.546255394747858*^9, 3.546255446772355*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"rhsEC", "[", RowBox[{ RowBox[{ RowBox[{"2", "/", " ", RowBox[{"Log", "[", "2", "]"}]}], "+", "0.001"}], ",", "4"}], "]"}]], "Input", CellChangeTimes->{{3.546253650118031*^9, 3.546253654996784*^9}, { 3.546256001470154*^9, 3.546256001948708*^9}}], Cell[BoxData[ RowBox[{"-", "0.11796891787963493`"}]], "Output", CellChangeTimes->{{3.546253651095337*^9, 3.546253655846133*^9}, 3.546256012390728*^9}] }, Open ]], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{ RowBox[{ "since", " ", "the", " ", "other", " ", "side", " ", "of", " ", "the", " ", "explicit", " ", "formula", " ", "is", " ", "non"}], "-", "negative"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"we", " ", "have", " ", "a", " ", RowBox[{"contradiction", ".", " ", "Therefore"}], " ", "there", " ", "is", " ", "no", " ", "L"}], "-", RowBox[{ "function", " ", "with", " ", "that", "\[IndentingNewLine]", "funcitonal", " ", "equation"}]}], ",", " ", RowBox[{"and", " ", "therefore", " ", "no", " ", "elliptic", " ", RowBox[{"curve", "/", "Q"}], " ", "with", " ", "that", " ", RowBox[{"conductor", "."}]}]}], " ", "*)"}]], "Input", CellChangeTimes->{{3.546255848289275*^9, 3.546255889900751*^9}, { 3.546255934724976*^9, 3.5462559691968727`*^9}}], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{"now", " ", RowBox[{"let", "'"}], "s", " ", "show", " ", "there", " ", "is", " ", "no", " ", "elliptic", " ", "curve", " ", "with", " ", "conductor", " ", "7"}], " ", "*)"}]], "Input", CellChangeTimes->{{3.5462559771841516`*^9, 3.546255992965218*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"rhsEC", "[", RowBox[{ RowBox[{ RowBox[{"2", "/", " ", RowBox[{"Log", "[", "2", "]"}]}], "+", "0.001"}], ",", "7"}], "]"}]], "Input", CellChangeTimes->{3.546256026732745*^9}], Cell[BoxData["0.2551082740773134`"], "Output", CellChangeTimes->{3.546256027506967*^9}] }, Open ]], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{"that", " ", RowBox[{"didn", "'"}], "t", " ", RowBox[{"work", ".", " ", "So"}], " ", RowBox[{"let", "'"}], "s", " ", "adjust", " ", "the", " ", "test", " ", "function"}], " ", "*)"}]], "Input", CellChangeTimes->{{3.546256030628879*^9, 3.546256035379528*^9}, { 3.546256136011711*^9, 3.546256147716567*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"rhs", "[", RowBox[{"1.5", ",", "7"}], "]"}], "//", "Expand"}]], "Input", CellChangeTimes->{{3.546253474166198*^9, 3.546253483259068*^9}, { 3.546253517796421*^9, 3.546253541796481*^9}, {3.546253589596895*^9, 3.546253634251959*^9}, {3.546253697964293*^9, 3.546253707836171*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"-", "0.4700487439574441`"}], "-", RowBox[{"0.20871629314320228`", " ", RowBox[{"c", "[", "2", "]"}]}], "-", RowBox[{"0.008399394572698272`", " ", RowBox[{"c", "[", "3", "]"}]}]}]], "Output", CellChangeTimes->{{3.546253499402673*^9, 3.5462535433367157`*^9}, { 3.5462535907920113`*^9, 3.546253636457519*^9}, {3.546253699479548*^9, 3.54625370867389*^9}}] }, Open ]], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{ RowBox[{"the", " ", "Hasse", " ", "bound", " ", "says"}], "\[IndentingNewLine]", "|", RowBox[{"c", "[", "2", "]"}], "|", RowBox[{ RowBox[{"\\", "le"}], " ", RowBox[{"Sqrt", "[", "2", "]"}], RowBox[{ RowBox[{"Log", "[", "2", "]"}], " ", "\\", "approx"}], " ", "0.98025"}], "\[IndentingNewLine]", "|", RowBox[{"c", "[", "3", "]"}], "|", " ", RowBox[{ RowBox[{"\\", "le"}], " ", RowBox[{"Sqrt", "[", "3", "]"}], RowBox[{ RowBox[{"Log", "[", "3", "]"}], " ", "\\", "approc"}], " ", RowBox[{"1.9028", ".", "\[IndentingNewLine]", "Check"}], " ", "that", " ", "this", " ", "implies", " ", "the", " ", "obove", " ", "quantity", " ", "is", " ", "negative"}]}], ",", "\[IndentingNewLine]", RowBox[{"so", " ", "there", " ", "is", " ", "no", " ", "elliptic", " ", RowBox[{"curve", "/", "Q"}], " ", "with", " ", "conductor", " ", "7"}]}], " ", "*)"}]], "Input", CellChangeTimes->{{3.546255068681136*^9, 3.546255082011627*^9}, { 3.546255185653804*^9, 3.546255186464798*^9}, {3.546257435380924*^9, 3.546257516019465*^9}, {3.5462575550594273`*^9, 3.5462575914498262`*^9}, 3.546257756292651*^9}] }, WindowSize->{800, 552}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, ShowSelection->True, FrontEndVersion->"8.0 for Mac OS X x86 (32-bit, 64-bit Kernel) (October 5, \ 2011)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[545, 20, 616, 12, 43, "Input"], Cell[1164, 34, 909, 19, 73, "Input"], Cell[2076, 55, 193, 4, 27, "Input"], Cell[2272, 61, 652, 20, 27, "Input"], Cell[2927, 83, 392, 8, 27, "Input"], Cell[3322, 93, 2079, 55, 88, "Input"], Cell[5404, 150, 510, 14, 43, "Input"], Cell[5917, 166, 953, 26, 58, "Input"], Cell[6873, 194, 1231, 30, 43, "Input"], Cell[8107, 226, 668, 16, 43, "Input"], Cell[8778, 244, 736, 18, 43, "Input"], Cell[9517, 264, 224, 6, 27, "Input"], Cell[9744, 272, 299, 8, 27, "Input"], Cell[10046, 282, 685, 14, 58, "Input"], Cell[10734, 298, 449, 14, 27, "Input"], Cell[11186, 314, 265, 9, 27, "Input"], Cell[CellGroupData[{ Cell[11476, 327, 144, 3, 27, "Input"], Cell[11623, 332, 92, 1, 27, "Output"] }, Open ]], Cell[11730, 336, 1041, 31, 43, "Input"], Cell[12774, 369, 307, 10, 27, "Input"], Cell[CellGroupData[{ Cell[13106, 383, 120, 3, 27, "Input"], Cell[13229, 388, 94, 1, 27, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[13360, 394, 255, 7, 27, "Input"], Cell[13618, 403, 10185, 173, 247, "Output"] }, Open ]], Cell[23818, 579, 582, 15, 27, "Input"], Cell[24403, 596, 642, 14, 43, "Input"], Cell[25048, 612, 457, 8, 27, "Input"], Cell[25508, 622, 383, 7, 27, "Input"], Cell[25894, 631, 480, 9, 43, "Input"], Cell[CellGroupData[{ Cell[26399, 644, 290, 8, 27, "Input"], Cell[26692, 654, 156, 3, 27, "Output"] }, Open ]], Cell[26863, 660, 878, 19, 58, "Input"], Cell[27744, 681, 308, 6, 27, "Input"], Cell[CellGroupData[{ Cell[28077, 691, 217, 7, 27, "Input"], Cell[28297, 700, 88, 1, 27, "Output"] }, Open ]], Cell[28400, 704, 370, 8, 27, "Input"], Cell[CellGroupData[{ Cell[28795, 716, 321, 6, 27, "Input"], Cell[29119, 724, 406, 9, 27, "Output"] }, Open ]], Cell[29540, 736, 1258, 28, 91, "Input"] } ] *) (* End of internal cache information *)