(* 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[     59900,       1824]
NotebookOptionsPosition[     52345,       1585]
NotebookOutlinePosition[     52689,       1600]
CellTagsIndexPosition[     52646,       1597]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"\[CurlyPi]", "=", 
   RowBox[{"1", "/", "1168"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.579039757775041*^9, 3.5790397592201433`*^9}, {
  3.5790416585652437`*^9, 3.5790416627056007`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"BinarySearch", "[", 
    RowBox[{"f_", ",", "a_", ",", "b_"}], "]"}], "/;", 
   RowBox[{"b", "\[Equal]", 
    RowBox[{"a", "+", "1"}]}]}], ":=", "b"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"BinarySearch", "[", 
   RowBox[{"f_", ",", "a_", ",", "b_"}], "]"}], ":=", 
  RowBox[{"(", 
   RowBox[{
    RowBox[{"Print", "[", 
     RowBox[{"{", 
      RowBox[{"a", ",", "b"}], "}"}], "]"}], ";", 
    RowBox[{"With", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"m", "=", 
        RowBox[{"Floor", "[", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"a", "+", "b"}], ")"}], "/", "2"}], "]"}]}], "}"}], ",", 
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{"f", "[", "m", "]"}], ",", 
        RowBox[{"BinarySearch", "[", 
         RowBox[{"f", ",", "a", ",", "m"}], "]"}], ",", 
        RowBox[{"BinarySearch", "[", 
         RowBox[{"f", ",", "m", ",", "b"}], "]"}]}], "]"}]}], "]"}]}], 
   ")"}]}]}], "Input",
 CellChangeTimes->{{3.579055500871192*^9, 3.579055641069498*^9}, {
  3.579055723971489*^9, 3.579055788699191*^9}}],

Cell[CellGroupData[{

Cell["Before Terry' s comment #47", "Section",
 CellChangeTimes->{{3.5792855132606*^9, 3.57928552269585*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"\[Omega]", "[", 
   RowBox[{"k0_", ",", "l0_", ",", "\[CurlyPi]_"}], "]"}], ":=", 
  RowBox[{
   FractionBox["1", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"k0", "+", 
       RowBox[{"2", "l0"}]}], ")"}], "!"}]], 
   RowBox[{"Binomial", "[", 
    RowBox[{
     RowBox[{"2", "l0"}], ",", "l0"}], "]"}], 
   RowBox[{"(", 
    RowBox[{
     RowBox[{
      FractionBox[
       RowBox[{"2", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"2", "l0"}], "+", "1"}], ")"}]}], 
       RowBox[{"l0", "+", "1"}]], 
      FractionBox[
       RowBox[{"k0", 
        RowBox[{"(", 
         RowBox[{"1", "-", 
          RowBox[{"\[Kappa]2", "[", 
           RowBox[{"k0", ",", "l0", ",", "\[CurlyPi]"}], "]"}]}], ")"}]}], 
       RowBox[{"k0", "+", 
        RowBox[{"2", "l0"}], "+", "1"}]]}], "-", 
     FractionBox[
      RowBox[{"4", 
       RowBox[{"(", 
        RowBox[{"1", "+", 
         RowBox[{"\[Kappa]1", "[", 
          RowBox[{"k0", ",", "l0", ",", "\[CurlyPi]"}], "]"}]}], ")"}]}], 
      RowBox[{"1", "+", 
       RowBox[{"4", "\[CurlyPi]"}]}]]}], ")"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Kappa]1", "[", 
   RowBox[{"k0_", ",", "l0_", ",", "\[CurlyPi]_"}], "]"}], ":=", 
  RowBox[{
   RowBox[{"\[Delta]1", "[", 
    RowBox[{"k0", ",", "\[CurlyPi]"}], "]"}], 
   RowBox[{"(", 
    RowBox[{"1", "+", 
     SuperscriptBox[
      RowBox[{"\[Delta]2", "[", "k0", "]"}], "2"], "+", 
     RowBox[{
      RowBox[{"Log", "[", "293", "]"}], "k0"}]}], ")"}], 
   RowBox[{"Binomial", "[", 
    RowBox[{
     RowBox[{"k0", "+", 
      RowBox[{"2", "l0"}]}], ",", "k0"}], "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Kappa]2", "[", 
   RowBox[{"k0_", ",", "l0_", ",", "\[CurlyPi]_"}], "]"}], ":=", 
  RowBox[{
   RowBox[{"\[Delta]1", "[", 
    RowBox[{"k0", ",", "\[CurlyPi]"}], "]"}], 
   RowBox[{"(", 
    RowBox[{"1", "+", 
     RowBox[{"4", "\[CurlyPi]"}]}], ")"}], 
   RowBox[{"(", 
    RowBox[{"1", "+", 
     SuperscriptBox[
      RowBox[{"\[Delta]2", "[", "k0", "]"}], "2"], "+", 
     RowBox[{
      RowBox[{"Log", "[", "293", "]"}], "k0"}]}], ")"}], 
   RowBox[{"Binomial", "[", 
    RowBox[{
     RowBox[{"k0", "+", 
      RowBox[{"2", "l0"}], "+", "1"}], ",", 
     RowBox[{"k0", "-", "1"}]}], "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Delta]1", "[", 
   RowBox[{"k0_", ",", "\[CurlyPi]_"}], "]"}], ":=", 
  SuperscriptBox[
   RowBox[{"(", 
    RowBox[{"1", "+", 
     RowBox[{"4", "\[CurlyPi]"}]}], ")"}], 
   RowBox[{"-", "k0"}]]}]}], "Input",
 CellChangeTimes->{{3.5790417567917356`*^9, 3.579041842111093*^9}, {
   3.5790418806094007`*^9, 3.579041995269519*^9}, {3.579042036677415*^9, 
   3.579042069850539*^9}, {3.579042162504923*^9, 3.579042300278743*^9}, {
   3.579055344517562*^9, 3.579055346903365*^9}, 3.5790553835164223`*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"\[Delta]2", "[", "k0_", "]"}], ":=", 
  RowBox[{
   RowBox[{"\[Delta]2", "[", "k0", "]"}], "=", 
   RowBox[{"Sum", "[", 
    RowBox[{
     FractionBox[
      SuperscriptBox[
       RowBox[{"(", 
        RowBox[{
         RowBox[{"Log", "[", "293", "]"}], "k0"}], ")"}], "\[Nu]"], 
      RowBox[{"\[Nu]", "!"}]], ",", 
     RowBox[{"{", 
      RowBox[{"\[Nu]", ",", "0", ",", "292"}], "}"}]}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.5790417567917356`*^9, 3.579041842111093*^9}, {
  3.5790418806094007`*^9, 3.579041995269519*^9}, {3.579042036677415*^9, 
  3.579042069850539*^9}, {3.579042162504923*^9, 3.579042300278743*^9}, {
  3.579055344517562*^9, 3.579055346903365*^9}, {3.5790553835164223`*^9, 
  3.5790553852288837`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{
     RowBox[{"Terry", "'"}], "s", " ", "improved", " ", "\[Delta]2", " ", 
     "from", " ", 
     RowBox[{"http", ":"}]}], "//", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{"sbseminar", ".", "wordpress", ".", "com"}], "/", "2013"}], 
        "/", "05"}], "/", "30"}], "/", "i"}], "-", "just", "-", "cant", "-", 
     "resist", "-", "there", "-", "are", "-", "infinitely", "-", "many", "-", 
     "pairs", "-", "of", "-", "primes", "-", "at", "-", "most", "-", 
     "59470640", "-", 
     RowBox[{
      RowBox[{"apart", "/", "#"}], "comment"}], "-", "23465"}]}], " ", "*)"}],
   "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Clear", "[", "\[Delta]2", "]"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Delta]2", "[", "k0_", "]"}], ":=", 
    RowBox[{
     RowBox[{"\[Delta]2", "[", "k0", "]"}], "=", 
     RowBox[{"Product", "[", 
      RowBox[{
       RowBox[{"1", "+", 
        RowBox[{"k0", " ", 
         RowBox[{"Log", "[", 
          RowBox[{"1", "+", 
           FractionBox["1", "j"]}], "]"}]}]}], ",", 
       RowBox[{"{", 
        RowBox[{"j", ",", "1", ",", "292"}], "}"}]}], "]"}]}]}]}]}]], "Input",\

 CellChangeTimes->{{3.579055315519124*^9, 3.579055431179329*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{
     RowBox[{"Terry", "'"}], "s", " ", "next", " ", "improvement"}], ",", " ", 
    RowBox[{"based", " ", "on", " ", "estimating", " ", 
     SubscriptBox["\[CapitalSigma]", "2"], " ", "more", " ", "carefully"}]}], 
   " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"\[Delta]2", "[", 
     RowBox[{"k0_", ",", "j_"}], "]"}], ":=", 
    RowBox[{
     RowBox[{"\[Delta]2", "[", 
      RowBox[{"k0", ",", "j"}], "]"}], "=", 
     RowBox[{"Product", "[", 
      RowBox[{
       RowBox[{"1", "+", 
        RowBox[{"k0", " ", 
         RowBox[{"Log", "[", 
          RowBox[{"1", "+", 
           FractionBox["1", "i"]}], "]"}]}]}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "j"}], "}"}]}], "]"}]}]}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Kappa]1", "[", 
     RowBox[{"k0_", ",", "l0_", ",", "\[CurlyPi]_"}], "]"}], ":=", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"\[Delta]1", "[", 
        RowBox[{"k0", ",", "\[CurlyPi]"}], "]"}], "+", 
       RowBox[{"Sum", "[", 
        RowBox[{
         RowBox[{
          SuperscriptBox[
           RowBox[{"\[Delta]1", "[", 
            RowBox[{"k0", ",", "\[CurlyPi]"}], "]"}], "j"], 
          SuperscriptBox[
           RowBox[{"\[Delta]2", "[", 
            RowBox[{"k0", ",", "j"}], "]"}], "2"]}], ",", 
         RowBox[{"{", 
          RowBox[{"j", ",", "1", ",", 
           RowBox[{"1", "/", 
            RowBox[{"(", 
             RowBox[{"4", "\[CurlyPi]"}], ")"}]}]}], "}"}]}], "]"}], "+", 
       RowBox[{
        RowBox[{"\[Delta]1", "[", 
         RowBox[{"k0", ",", "\[CurlyPi]"}], "]"}], "k0", " ", 
        RowBox[{"Log", "[", 
         RowBox[{"1", "+", 
          FractionBox["1", 
           RowBox[{"4", "\[CurlyPi]"}]]}], "]"}]}]}], ")"}], 
     RowBox[{"Binomial", "[", 
      RowBox[{
       RowBox[{"k0", "+", 
        RowBox[{"2", "l0"}]}], ",", "k0"}], "]"}]}]}], "\[IndentingNewLine]", 
   
   RowBox[{
    RowBox[{"\[Kappa]2", "[", 
     RowBox[{"k0_", ",", "l0_", ",", "\[CurlyPi]_"}], "]"}], ":=", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"\[Delta]1", "[", 
         RowBox[{"k0", ",", "\[CurlyPi]"}], "]"}], 
        RowBox[{"(", 
         RowBox[{"1", "+", 
          RowBox[{"4", "\[CurlyPi]"}]}], ")"}]}], "+", 
       RowBox[{"Sum", "[", 
        RowBox[{
         RowBox[{
          SuperscriptBox[
           RowBox[{"\[Delta]1", "[", 
            RowBox[{"k0", ",", "\[CurlyPi]"}], "]"}], "j"], 
          RowBox[{"(", 
           RowBox[{"1", "+", 
            RowBox[{"4", "\[CurlyPi]"}]}], ")"}], 
          SuperscriptBox[
           RowBox[{"\[Delta]2", "[", 
            RowBox[{"k0", ",", "j"}], "]"}], "2"]}], ",", 
         RowBox[{"{", 
          RowBox[{"j", ",", "1", ",", 
           RowBox[{"1", "/", 
            RowBox[{"(", 
             RowBox[{"4", "\[CurlyPi]"}], ")"}]}]}], "}"}]}], "]"}], "+", 
       RowBox[{
        RowBox[{"\[Delta]1", "[", 
         RowBox[{"k0", ",", "\[CurlyPi]"}], "]"}], 
        RowBox[{"(", 
         RowBox[{"1", "+", 
          RowBox[{"4", "\[CurlyPi]"}]}], ")"}], "k0", " ", 
        RowBox[{"Log", "[", 
         RowBox[{"1", "+", 
          FractionBox["1", 
           RowBox[{"4", "\[CurlyPi]"}]]}], "]"}]}]}], ")"}], 
     RowBox[{"Binomial", "[", 
      RowBox[{
       RowBox[{"k0", "+", 
        RowBox[{"2", "l0"}], "+", "1"}], ",", 
       RowBox[{"k0", "-", "1"}]}], "]"}]}]}]}]}]], "Input",
 CellChangeTimes->{{3.57915658298456*^9, 3.579156880107204*^9}, {
  3.579156923755351*^9, 3.5791569277007513`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"k0", "=", 
   RowBox[{"3.5", "*", 
    SuperscriptBox["10", "6"]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.579042425475142*^9, 3.579042430342516*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Select", "[", 
  RowBox[{
   RowBox[{"Range", "[", 
    RowBox[{"140", ",", "680"}], "]"}], ",", 
   RowBox[{
    RowBox[{
     RowBox[{"\[Omega]", "[", 
      RowBox[{"k0", ",", "#", ",", "\[CurlyPi]"}], "]"}], ">", "0"}], "&"}]}],
   "]"}]], "Input",
 CellChangeTimes->{{3.579042499997485*^9, 3.579042524579956*^9}, {
  3.579055449436914*^9, 3.579055449658369*^9}, {3.57915695283364*^9, 
  3.579156997478168*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "150", ",", "151", ",", "152", ",", "153", ",", "154", ",", "155", ",", 
   "156", ",", "157", ",", "158", ",", "159", ",", "160", ",", "161", ",", 
   "162", ",", "163", ",", "164", ",", "165", ",", "166", ",", "167", ",", 
   "168", ",", "169", ",", "170", ",", "171", ",", "172", ",", "173", ",", 
   "174", ",", "175", ",", "176", ",", "177", ",", "178", ",", "179", ",", 
   "180", ",", "181", ",", "182", ",", "183", ",", "184", ",", "185", ",", 
   "186", ",", "187", ",", "188", ",", "189", ",", "190", ",", "191", ",", 
   "192", ",", "193", ",", "194", ",", "195", ",", "196", ",", "197", ",", 
   "198", ",", "199", ",", "200", ",", "201", ",", "202", ",", "203", ",", 
   "204", ",", "205", ",", "206", ",", "207", ",", "208", ",", "209", ",", 
   "210", ",", "211", ",", "212", ",", "213", ",", "214", ",", "215", ",", 
   "216", ",", "217", ",", "218", ",", "219", ",", "220", ",", "221", ",", 
   "222", ",", "223", ",", "224", ",", "225", ",", "226", ",", "227", ",", 
   "228", ",", "229", ",", "230", ",", "231", ",", "232", ",", "233", ",", 
   "234", ",", "235", ",", "236", ",", "237", ",", "238", ",", "239", ",", 
   "240", ",", "241", ",", "242", ",", "243", ",", "244", ",", "245", ",", 
   "246", ",", "247", ",", "248", ",", "249", ",", "250", ",", "251", ",", 
   "252", ",", "253", ",", "254", ",", "255", ",", "256", ",", "257", ",", 
   "258", ",", "259", ",", "260", ",", "261", ",", "262", ",", "263", ",", 
   "264", ",", "265", ",", "266", ",", "267", ",", "268", ",", "269", ",", 
   "270", ",", "271", ",", "272", ",", "273", ",", "274", ",", "275", ",", 
   "276", ",", "277", ",", "278", ",", "279", ",", "280", ",", "281", ",", 
   "282", ",", "283", ",", "284", ",", "285", ",", "286", ",", "287", ",", 
   "288", ",", "289", ",", "290", ",", "291", ",", "292", ",", "293", ",", 
   "294", ",", "295", ",", "296", ",", "297", ",", "298", ",", "299", ",", 
   "300", ",", "301", ",", "302", ",", "303", ",", "304", ",", "305", ",", 
   "306", ",", "307", ",", "308", ",", "309", ",", "310", ",", "311", ",", 
   "312", ",", "313", ",", "314", ",", "315", ",", "316", ",", "317", ",", 
   "318", ",", "319", ",", "320", ",", "321", ",", "322", ",", "323", ",", 
   "324", ",", "325", ",", "326", ",", "327", ",", "328", ",", "329", ",", 
   "330", ",", "331", ",", "332", ",", "333", ",", "334", ",", "335", ",", 
   "336", ",", "337", ",", "338", ",", "339", ",", "340", ",", "341", ",", 
   "342", ",", "343", ",", "344", ",", "345", ",", "346", ",", "347", ",", 
   "348", ",", "349", ",", "350", ",", "351", ",", "352", ",", "353", ",", 
   "354", ",", "355", ",", "356", ",", "357", ",", "358", ",", "359", ",", 
   "360", ",", "361", ",", "362", ",", "363", ",", "364", ",", "365", ",", 
   "366", ",", "367", ",", "368", ",", "369", ",", "370", ",", "371", ",", 
   "372", ",", "373", ",", "374", ",", "375", ",", "376", ",", "377", ",", 
   "378", ",", "379", ",", "380", ",", "381", ",", "382", ",", "383", ",", 
   "384", ",", "385", ",", "386", ",", "387", ",", "388", ",", "389", ",", 
   "390", ",", "391", ",", "392", ",", "393", ",", "394", ",", "395", ",", 
   "396", ",", "397", ",", "398", ",", "399", ",", "400", ",", "401", ",", 
   "402", ",", "403", ",", "404", ",", "405", ",", "406", ",", "407", ",", 
   "408", ",", "409", ",", "410", ",", "411", ",", "412", ",", "413", ",", 
   "414", ",", "415", ",", "416", ",", "417", ",", "418", ",", "419", ",", 
   "420", ",", "421", ",", "422", ",", "423", ",", "424", ",", "425", ",", 
   "426", ",", "427", ",", "428", ",", "429", ",", "430", ",", "431", ",", 
   "432", ",", "433", ",", "434", ",", "435", ",", "436", ",", "437", ",", 
   "438", ",", "439", ",", "440", ",", "441", ",", "442", ",", "443", ",", 
   "444", ",", "445", ",", "446", ",", "447", ",", "448", ",", "449", ",", 
   "450", ",", "451", ",", "452", ",", "453", ",", "454", ",", "455", ",", 
   "456", ",", "457", ",", "458", ",", "459", ",", "460", ",", "461", ",", 
   "462", ",", "463", ",", "464", ",", "465", ",", "466", ",", "467", ",", 
   "468", ",", "469", ",", "470", ",", "471", ",", "472", ",", "473", ",", 
   "474", ",", "475", ",", "476", ",", "477", ",", "478", ",", "479", ",", 
   "480", ",", "481", ",", "482", ",", "483", ",", "484", ",", "485", ",", 
   "486", ",", "487", ",", "488", ",", "489", ",", "490", ",", "491", ",", 
   "492", ",", "493", ",", "494", ",", "495", ",", "496", ",", "497", ",", 
   "498", ",", "499", ",", "500", ",", "501", ",", "502", ",", "503", ",", 
   "504", ",", "505", ",", "506", ",", "507", ",", "508", ",", "509", ",", 
   "510", ",", "511", ",", "512", ",", "513", ",", "514", ",", "515", ",", 
   "516", ",", "517", ",", "518", ",", "519", ",", "520", ",", "521", ",", 
   "522", ",", "523", ",", "524", ",", "525", ",", "526", ",", "527", ",", 
   "528", ",", "529", ",", "530", ",", "531", ",", "532", ",", "533", ",", 
   "534", ",", "535", ",", "536", ",", "537", ",", "538", ",", "539", ",", 
   "540", ",", "541", ",", "542", ",", "543", ",", "544", ",", "545", ",", 
   "546", ",", "547", ",", "548", ",", "549", ",", "550", ",", "551", ",", 
   "552", ",", "553", ",", "554", ",", "555", ",", "556", ",", "557", ",", 
   "558", ",", "559", ",", "560", ",", "561", ",", "562", ",", "563", ",", 
   "564", ",", "565", ",", "566", ",", "567", ",", "568", ",", "569", ",", 
   "570", ",", "571", ",", "572", ",", "573", ",", "574", ",", "575", ",", 
   "576", ",", "577", ",", "578", ",", "579", ",", "580", ",", "581", ",", 
   "582", ",", "583", ",", "584", ",", "585", ",", "586", ",", "587", ",", 
   "588", ",", "589", ",", "590", ",", "591", ",", "592", ",", "593", ",", 
   "594", ",", "595", ",", "596", ",", "597", ",", "598", ",", "599", ",", 
   "600", ",", "601", ",", "602", ",", "603", ",", "604", ",", "605", ",", 
   "606", ",", "607", ",", "608", ",", "609", ",", "610", ",", "611", ",", 
   "612", ",", "613", ",", "614", ",", "615", ",", "616", ",", "617", ",", 
   "618", ",", "619", ",", "620", ",", "621", ",", "622", ",", "623", ",", 
   "624", ",", "625", ",", "626", ",", "627", ",", "628", ",", "629", ",", 
   "630", ",", "631", ",", "632", ",", "633", ",", "634", ",", "635", ",", 
   "636", ",", "637", ",", "638", ",", "639", ",", "640", ",", "641", ",", 
   "642", ",", "643", ",", "644", ",", "645", ",", "646", ",", "647", ",", 
   "648", ",", "649", ",", "650", ",", "651", ",", "652", ",", "653", ",", 
   "654", ",", "655", ",", "656", ",", "657", ",", "658", ",", "659", ",", 
   "660", ",", "661", ",", "662", ",", "663", ",", "664", ",", "665", ",", 
   "666", ",", "667", ",", "668", ",", "669", ",", "670", ",", "671", ",", 
   "672"}], "}"}]], "Output",
 CellChangeTimes->{
  3.5790425254130983`*^9, {3.579055445009554*^9, 3.579055450192212*^9}, {
   3.5791569323830853`*^9, 3.579157011460746*^9}}]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"BinarySearch", "[", 
    RowBox[{"f_", ",", "a_", ",", "b_"}], "]"}], "/;", 
   RowBox[{"b", "\[Equal]", 
    RowBox[{"a", "+", "1"}]}]}], ":=", "b"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"BinarySearch", "[", 
   RowBox[{"f_", ",", "a_", ",", "b_"}], "]"}], ":=", 
  RowBox[{"(", 
   RowBox[{
    RowBox[{"Print", "[", 
     RowBox[{"{", 
      RowBox[{"a", ",", "b"}], "}"}], "]"}], ";", 
    RowBox[{"With", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"m", "=", 
        RowBox[{"Floor", "[", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"a", "+", "b"}], ")"}], "/", "2"}], "]"}]}], "}"}], ",", 
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{"f", "[", "m", "]"}], ",", 
        RowBox[{"BinarySearch", "[", 
         RowBox[{"f", ",", "a", ",", "m"}], "]"}], ",", 
        RowBox[{"BinarySearch", "[", 
         RowBox[{"f", ",", "m", ",", "b"}], "]"}]}], "]"}]}], "]"}]}], 
   ")"}]}]}], "Input",
 CellChangeTimes->{{3.579055500871192*^9, 3.579055641069498*^9}, {
  3.579055723971489*^9, 3.579055788699191*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"BinarySearch", "[", 
  RowBox[{
   RowBox[{"Function", "[", 
    RowBox[{"k", ",", 
     RowBox[{
      RowBox[{"Length", "[", 
       RowBox[{"Select", "[", 
        RowBox[{
         RowBox[{"Range", "[", 
          RowBox[{"150", ",", "672"}], "]"}], ",", 
         RowBox[{
          RowBox[{
           RowBox[{"\[Omega]", "[", 
            RowBox[{
             RowBox[{"N", "[", "k", "]"}], ",", 
             RowBox[{"N", "[", "#", "]"}], ",", "\[CurlyPi]"}], "]"}], ">", 
           "0"}], "&"}]}], "]"}], "]"}], ">", "0"}]}], "]"}], ",", "0", ",", 
   "2950000"}], "]"}]], "Input",
 CellChangeTimes->{{3.579055644274489*^9, 3.5790557466320133`*^9}, {
  3.579055815360876*^9, 3.5790558190871277`*^9}, {3.579157024734057*^9, 
  3.579157032916317*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0", ",", "2950000"}], "}"}]], "Print",
 CellChangeTimes->{{3.5790557338573112`*^9, 3.579055748193603*^9}, 
   3.579055790961153*^9, 3.5790558213953047`*^9, 3.579157035503875*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0", ",", "1475000"}], "}"}]], "Print",
 CellChangeTimes->{{3.5790557338573112`*^9, 3.579055748193603*^9}, 
   3.579055790961153*^9, 3.5790558213953047`*^9, 3.57915704880968*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"737500", ",", "1475000"}], "}"}]], "Print",
 CellChangeTimes->{{3.5790557338573112`*^9, 3.579055748193603*^9}, 
   3.579055790961153*^9, 3.5790558213953047`*^9, 3.579157061044039*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"737500", ",", "1106250"}], "}"}]], "Print",
 CellChangeTimes->{{3.5790557338573112`*^9, 3.579055748193603*^9}, 
   3.579055790961153*^9, 3.5790558213953047`*^9, 3.579157076334573*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"737500", ",", "921875"}], "}"}]], "Print",
 CellChangeTimes->{{3.5790557338573112`*^9, 3.579055748193603*^9}, 
   3.579055790961153*^9, 3.5790558213953047`*^9, 3.5791570893897943`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"829687", ",", "921875"}], "}"}]], "Print",
 CellChangeTimes->{{3.5790557338573112`*^9, 3.579055748193603*^9}, 
   3.579055790961153*^9, 3.5790558213953047`*^9, 3.5791571010273848`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"829687", ",", "875781"}], "}"}]], "Print",
 CellChangeTimes->{{3.5790557338573112`*^9, 3.579055748193603*^9}, 
   3.579055790961153*^9, 3.5790558213953047`*^9, 3.579157115121121*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"852734", ",", "875781"}], "}"}]], "Print",
 CellChangeTimes->{{3.5790557338573112`*^9, 3.579055748193603*^9}, 
   3.579055790961153*^9, 3.5790558213953047`*^9, 3.5791571289009457`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"864257", ",", "875781"}], "}"}]], "Print",
 CellChangeTimes->{{3.5790557338573112`*^9, 3.579055748193603*^9}, 
   3.579055790961153*^9, 3.5790558213953047`*^9, 3.579157143317355*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"864257", ",", "870019"}], "}"}]], "Print",
 CellChangeTimes->{{3.5790557338573112`*^9, 3.579055748193603*^9}, 
   3.579055790961153*^9, 3.5790558213953047`*^9, 3.5791571580702467`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"864257", ",", "867138"}], "}"}]], "Print",
 CellChangeTimes->{{3.5790557338573112`*^9, 3.579055748193603*^9}, 
   3.579055790961153*^9, 3.5790558213953047`*^9, 3.579157172853897*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"865697", ",", "867138"}], "}"}]], "Print",
 CellChangeTimes->{{3.5790557338573112`*^9, 3.579055748193603*^9}, 
   3.579055790961153*^9, 3.5790558213953047`*^9, 3.579157184635737*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"866417", ",", "867138"}], "}"}]], "Print",
 CellChangeTimes->{{3.5790557338573112`*^9, 3.579055748193603*^9}, 
   3.579055790961153*^9, 3.5790558213953047`*^9, 3.5791571981626263`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"866777", ",", "867138"}], "}"}]], "Print",
 CellChangeTimes->{{3.5790557338573112`*^9, 3.579055748193603*^9}, 
   3.579055790961153*^9, 3.5790558213953047`*^9, 3.5791572117951117`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"866777", ",", "866957"}], "}"}]], "Print",
 CellChangeTimes->{{3.5790557338573112`*^9, 3.579055748193603*^9}, 
   3.579055790961153*^9, 3.5790558213953047`*^9, 3.579157225703145*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"866777", ",", "866867"}], "}"}]], "Print",
 CellChangeTimes->{{3.5790557338573112`*^9, 3.579055748193603*^9}, 
   3.579055790961153*^9, 3.5790558213953047`*^9, 3.579157237419524*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"866777", ",", "866822"}], "}"}]], "Print",
 CellChangeTimes->{{3.5790557338573112`*^9, 3.579055748193603*^9}, 
   3.579055790961153*^9, 3.5790558213953047`*^9, 3.579157250671084*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"866799", ",", "866822"}], "}"}]], "Print",
 CellChangeTimes->{{3.5790557338573112`*^9, 3.579055748193603*^9}, 
   3.579055790961153*^9, 3.5790558213953047`*^9, 3.579157269500978*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"866799", ",", "866810"}], "}"}]], "Print",
 CellChangeTimes->{{3.5790557338573112`*^9, 3.579055748193603*^9}, 
   3.579055790961153*^9, 3.5790558213953047`*^9, 3.579157290410152*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"866804", ",", "866810"}], "}"}]], "Print",
 CellChangeTimes->{{3.5790557338573112`*^9, 3.579055748193603*^9}, 
   3.579055790961153*^9, 3.5790558213953047`*^9, 3.579157310083561*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"866804", ",", "866807"}], "}"}]], "Print",
 CellChangeTimes->{{3.5790557338573112`*^9, 3.579055748193603*^9}, 
   3.579055790961153*^9, 3.5790558213953047`*^9, 3.5791573305934753`*^9}]
}, Open  ]],

Cell[BoxData["866805"], "Output",
 CellChangeTimes->{{3.579055733875759*^9, 3.579055752002152*^9}, {
   3.5790557959167957`*^9, 3.579055822049107*^9}, 3.579157349674293*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Select", "[", 
  RowBox[{
   RowBox[{"Range", "[", 
    RowBox[{"150", ",", "208"}], "]"}], ",", 
   RowBox[{
    RowBox[{
     RowBox[{"\[Omega]", "[", 
      RowBox[{"866805.", ",", 
       RowBox[{"N", "[", "#", "]"}], ",", "\[CurlyPi]"}], "]"}], ">", "0"}], 
    "&"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.5790558392696342`*^9, 3.579055839788232*^9}, {
   3.5791573551823606`*^9, 3.579157356964855*^9}, 3.579157618077054*^9}],

Cell[BoxData[
 RowBox[{"{", "164", "}"}]], "Output",
 CellChangeTimes->{
  3.579055840321251*^9, {3.579157620373219*^9, 3.579157627176448*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"N", "[", 
  RowBox[{"\[Omega]", "[", 
   RowBox[{"866805", ",", "164", ",", "\[CurlyPi]"}], "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.579043024039431*^9, 3.579043033447566*^9}, {
   3.579043088864394*^9, 3.579043103990138*^9}, 3.579055862676907*^9, 
   3.5791574188191223`*^9, {3.579157647825068*^9, 3.579157648295678*^9}}],

Cell[BoxData["$Aborted"], "Output",
 CellChangeTimes->{3.579157876931909*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"N", "[", 
  RowBox[{"\[Omega]", "[", 
   RowBox[{
    RowBox[{"866805", "-", "1"}], ",", "164", ",", "\[CurlyPi]"}], "]"}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.5790431206296873`*^9, 3.579043120740706*^9}, {
  3.5790558694376698`*^9, 3.579055870191731*^9}, {3.5791576734403257`*^9, 
  3.5791576753335667`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"2", 
  RowBox[{"Prime", "[", 
   RowBox[{"36716", "+", 
    RowBox[{"Ceiling", "[", 
     RowBox[{"2618607", "/", "2"}], "]"}], "+", "1"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.5790558812410383`*^9, 3.5790559023640947`*^9}}],

Cell[BoxData["42543038"], "Output",
 CellChangeTimes->{{3.57905588294319*^9, 3.579055904823296*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Prime", "[", 
   RowBox[{"11651", "+", 
    RowBox[{"Ceiling", "[", 
     RowBox[{"866805", "/", "2"}], "]"}], "+", "1"}], "]"}], "+", 
  RowBox[{"Prime", "[", 
   RowBox[{"11651", "+", 
    RowBox[{"Floor", "[", 
     RowBox[{"866805", "/", "2"}], "]"}], "+", "1"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.579056746470296*^9, 3.57905674737389*^9}, {
  3.579157880673614*^9, 3.5791578938723783`*^9}, {3.579161875996298*^9, 
  3.579161878172415*^9}}],

Cell[BoxData["13008612"], "Output",
 CellChangeTimes->{3.579056748151313*^9, 3.579157896683292*^9, 
  3.579161884455652*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"PrimePi", "[", "866805", "]"}]], "Input",
 CellChangeTimes->{{3.579160218036357*^9, 3.579160222059648*^9}}],

Cell[BoxData["68850"], "Output",
 CellChangeTimes->{3.579160228323195*^9}]
}, Open  ]]
}, Closed]],

Cell[CellGroupData[{

Cell["After Terry' s comment #47", "Section",
 CellChangeTimes->{{3.5792855132606*^9, 3.579285537153631*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"\[Kappa]", "[", 
   RowBox[{"k0_", ",", "l0_", ",", "\[CurlyPi]_"}], "]"}], ":=", 
  RowBox[{"Sum", "[", 
   RowBox[{
    RowBox[{
     SuperscriptBox[
      RowBox[{"(", 
       RowBox[{"1", "-", 
        FractionBox[
         RowBox[{"2", "n", " ", "\[CurlyPi]"}], 
         RowBox[{"1", "+", 
          RowBox[{"4", "\[CurlyPi]"}]}]]}], ")"}], 
      RowBox[{
       RowBox[{"k0", "/", "2"}], "+", "l0"}]], 
     RowBox[{"Product", "[", 
      RowBox[{
       RowBox[{"1", "+", 
        RowBox[{"3", "k0", " ", 
         RowBox[{"Log", "[", 
          RowBox[{"1", "+", 
           FractionBox["1", "j"]}], "]"}]}]}], ",", 
       RowBox[{"{", 
        RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ",", 
    RowBox[{"{", 
     RowBox[{"n", ",", "1", ",", 
      RowBox[{"2", "+", 
       FractionBox["1", 
        RowBox[{"2", "\[CurlyPi]"}]], "-", "1"}]}], "}"}]}], "]"}]}]], "Input",\

 CellChangeTimes->{{3.579285544121467*^9, 3.5792856043331957`*^9}, {
  3.579285634869767*^9, 3.579285680945569*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"f", "[", 
   RowBox[{"k0_", ",", "l0_", ",", "\[CurlyPi]_"}], "]"}], ":=", 
  RowBox[{
   FractionBox[
    RowBox[{"1", "+", 
     RowBox[{"4", "\[CurlyPi]"}]}], 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"1", "+", 
       FractionBox["1", 
        RowBox[{
         RowBox[{"2", "l0"}], "+", "1"}]]}], ")"}], 
     RowBox[{"(", 
      RowBox[{"1", "+", 
       FractionBox[
        RowBox[{
         RowBox[{"2", "l0"}], "+", "1"}], "k0"]}], ")"}], 
     RowBox[{"(", 
      RowBox[{"1", "+", 
       RowBox[{"\[Kappa]", "[", 
        RowBox[{"k0", ",", "l0", ",", "\[CurlyPi]"}], "]"}]}], ")"}]}]], "-", 
   "1"}]}]], "Input",
 CellChangeTimes->{{3.579285793971397*^9, 3.579285844312901*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"g", "[", 
   RowBox[{"k0_", ",", "\[CurlyPi]_"}], "]"}], ":=", 
  RowBox[{"With", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"P", "=", "50"}], "}"}], ",", 
    RowBox[{"FindMaximum", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"f", "[", 
         RowBox[{
          RowBox[{"N", "[", 
           RowBox[{"k0", ",", "P"}], "]"}], ",", "l0", ",", "\[CurlyPi]"}], 
         "]"}], ",", 
        RowBox[{"l0", "\[GreaterEqual]", "0"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"l0", ",", "800"}], "}"}], ",", 
      RowBox[{"WorkingPrecision", "\[Rule]", "P"}]}], "]"}]}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.579285857362918*^9, 3.5792859001438427`*^9}, {
  3.5792859486695633`*^9, 3.57928595347364*^9}, {3.579286126022748*^9, 
  3.579286127780036*^9}, {3.579286169039331*^9, 3.57928617327332*^9}, {
  3.579286263796997*^9, 3.579286277228628*^9}, {3.579286354060665*^9, 
  3.5792863615472307`*^9}}],

Cell[CellGroupData[{

Cell["\[CurlyPi] = 1/1168;", "Subsection",
 CellChangeTimes->{{3.579039757775041*^9, 3.5790397592201433`*^9}, {
  3.5790416585652437`*^9, 3.5790416627056007`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"\[CurlyPi]", "=", 
   RowBox[{"1", "/", "1168"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.579039757775041*^9, 3.5790397592201433`*^9}, {
  3.5790416585652437`*^9, 3.5790416627056007`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"g", "[", 
  RowBox[{"3500000", ",", "\[CurlyPi]"}], "]"}]], "Input",
 CellChangeTimes->{{3.579285911907754*^9, 3.579285923379943*^9}, {
  3.57928614758582*^9, 3.5792861666068907`*^9}, {3.579286246639627*^9, 
  3.579286261922916*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "0.0023528109194062497496980395395230147214904504882306460795652316693646174\
9146`50.", ",", 
   RowBox[{"{", 
    RowBox[{
    "l0", "\[Rule]", 
     "934.91434669348616383545840511492385047653436461751291831633606702794774\
195927487`50."}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.579285921698407*^9, 3.579285963561235*^9}, {
   3.579286131884048*^9, 3.579286157050905*^9}, {3.5792862442790213`*^9, 
   3.579286249856029*^9}, 3.5792862820486927`*^9, 3.579379601247984*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"g", "[", 
  RowBox[{"400000", ",", "\[CurlyPi]"}], "]"}]], "Input",
 CellChangeTimes->{{3.5792862847573423`*^9, 3.579286289089642*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "0.0002590600067561710110331174237546869952395431028567097559074316616059028\
1761`50.", ",", 
   RowBox[{"{", 
    RowBox[{
    "l0", "\[Rule]", 
     "315.72776601683843259721715356107971665366223849802997048328627171926736\
314024163`50."}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.579286300856514*^9, 3.5792863694498997`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"g", "[", 
  RowBox[{"200000", ",", "\[CurlyPi]"}], "]"}]], "Input",
 CellChangeTimes->{{3.5792863791675777`*^9, 3.57928637923061*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{
   "-", "0.0010477873344049464965698275491370296591201639705588788650549472269\
1341322072`50."}], ",", 
   RowBox[{"{", 
    RowBox[{
    "l0", "\[Rule]", 
     "223.10679774997897862188563798282240984175992054555952795869793982298190\
607689417`50."}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.579286384049242*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"BinarySearch", "[", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"g", "[", 
       RowBox[{"#", ",", "\[CurlyPi]"}], "]"}], "\[LeftDoubleBracket]", "1", 
      "\[RightDoubleBracket]"}], ">", "0"}], "&"}], ",", "200000", ",", 
   "400000"}], "]"}]], "Input",
 CellChangeTimes->{{3.579286389308769*^9, 3.579286418064632*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"200000", ",", "400000"}], "}"}]], "Print",
 CellChangeTimes->{3.579286418612815*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"300000", ",", "400000"}], "}"}]], "Print",
 CellChangeTimes->{3.579286424913472*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"300000", ",", "350000"}], "}"}]], "Print",
 CellChangeTimes->{3.5792864313425817`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"325000", ",", "350000"}], "}"}]], "Print",
 CellChangeTimes->{3.579286438208674*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"337500", ",", "350000"}], "}"}]], "Print",
 CellChangeTimes->{3.579286450478159*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"337500", ",", "343750"}], "}"}]], "Print",
 CellChangeTimes->{3.579286459894164*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"340625", ",", "343750"}], "}"}]], "Print",
 CellChangeTimes->{3.5792864680590477`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"340625", ",", "342187"}], "}"}]], "Print",
 CellChangeTimes->{3.579286474916254*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"341406", ",", "342187"}], "}"}]], "Print",
 CellChangeTimes->{3.579286482766721*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"341406", ",", "341796"}], "}"}]], "Print",
 CellChangeTimes->{3.579286490391204*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"341601", ",", "341796"}], "}"}]], "Print",
 CellChangeTimes->{3.57928649838844*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"341601", ",", "341698"}], "}"}]], "Print",
 CellChangeTimes->{3.5792865065166893`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"341601", ",", "341649"}], "}"}]], "Print",
 CellChangeTimes->{3.579286514211459*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"341625", ",", "341649"}], "}"}]], "Print",
 CellChangeTimes->{3.579286522024201*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"341637", ",", "341649"}], "}"}]], "Print",
 CellChangeTimes->{3.579286530220573*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"341637", ",", "341643"}], "}"}]], "Print",
 CellChangeTimes->{3.579286538510601*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"341637", ",", "341640"}], "}"}]], "Print",
 CellChangeTimes->{3.579286545363533*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"341638", ",", "341640"}], "}"}]], "Print",
 CellChangeTimes->{3.579286552103311*^9}]
}, Open  ]],

Cell[BoxData["341640"], "Output",
 CellChangeTimes->{3.579286558819166*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"g", "[", 
  RowBox[{"%", ",", "\[CurlyPi]"}], "]"}]], "Input",
 CellChangeTimes->{{3.579286432507218*^9, 3.579286434664695*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "1.2476767753529155011545618741486182064925569372710316088116358956130361141\
76963161120955`50.*^-9", ",", 
   RowBox[{"{", 
    RowBox[{
    "l0", "\[Rule]", 
     "291.74989307098171801846152913931073901951938706798888767547146797840699\
864179031`50."}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.579286570247089*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"f", "[", 
  RowBox[{"341640.", ",", "291", ",", "\[CurlyPi]"}], "]"}]], "Input",
 CellChangeTimes->{{3.579286619373083*^9, 3.579286627953705*^9}}],

Cell[BoxData[
 RowBox[{"-", "1.0007093376174225`*^-8"}]], "Output",
 CellChangeTimes->{{3.5792866262749166`*^9, 3.579286629444607*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"f", "[", 
  RowBox[{"341640.", ",", "291.7", ",", "\[CurlyPi]"}], "]"}]], "Input",
 CellChangeTimes->{{3.579286670826133*^9, 3.57928667548426*^9}}],

Cell[BoxData["1.1979741643131092`*^-9"], "Output",
 CellChangeTimes->{{3.579286673106853*^9, 3.579286676847947*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"N", "[", 
  RowBox[{"f", "[", 
   RowBox[{"341640", ",", 
    FractionBox["2917", "10"], ",", "\[CurlyPi]"}], "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.579286712902433*^9, 3.579286727188678*^9}}],

Cell[BoxData["$Aborted"], "Output",
 CellChangeTimes->{3.5792868221967373`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"f", "[", 
  RowBox[{
   RowBox[{"N", "[", 
    RowBox[{"341640", ",", "50"}], "]"}], ",", "292", ",", "\[CurlyPi]"}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.579286633232902*^9, 3.5792866439504347`*^9}}],

Cell[BoxData["0``52.76715586608218"], "Output",
 CellChangeTimes->{{3.57928663487868*^9, 3.579286646461955*^9}}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{
   RowBox[{
   "This", " ", "is", " ", "with", " ", "an", " ", "older", " ", "value", " ",
     "of", " ", "m"}], ",", " ", 
   RowBox[{
    RowBox[{"that", " ", "sufficed", " ", "for", " ", "k_", "0"}], " ", "=", 
    " ", "866"}], ",", 
   RowBox[{"805", ":"}]}], " ", "*)"}]], "Input",
 CellChangeTimes->{{3.579286912272828*^9, 3.579286931571913*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"2", 
  RowBox[{"Prime", "[", 
   RowBox[{"11651", "+", 
    RowBox[{"341640", "/", "2"}], "-", "1"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.5792867958510427`*^9, 3.579286821051568*^9}}],

Cell[BoxData["4982086"], "Output",
 CellChangeTimes->{3.579286822849308*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
\[CurlyPi] = 1/1040 (* Conjectured feasible by v08ltu at \
http://mathoverflow.net/questions/132632/tightening-zhangs-bound-closed *)\
\>", "Subsection",
 CellChangeTimes->{{3.579039757775041*^9, 3.5790397592201433`*^9}, {
  3.5790416585652437`*^9, 3.5790416627056007`*^9}, {3.5793795333589163`*^9, 
  3.579379533479011*^9}, {3.579379667581636*^9, 3.579379690960455*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"\[CurlyPi]", "=", 
   RowBox[{"1", "/", "1040"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.579039757775041*^9, 3.5790397592201433`*^9}, {
  3.5790416585652437`*^9, 3.5790416627056007`*^9}, {3.579379538064554*^9, 
  3.5793795383210287`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"g", "[", 
  RowBox[{"400000", ",", "\[CurlyPi]"}], "]"}]], "Input",
 CellChangeTimes->{{3.5792862847573423`*^9, 3.579286289089642*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "0.0006792265848634983603682849186032094065231409647265030286322478277076514\
5577`50.", ",", 
   RowBox[{"{", 
    RowBox[{
    "l0", "\[Rule]", 
     "315.72776601683843235775236336539921676863061035549075867932892504114998\
114205919`50."}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.579286300856514*^9, 3.5792863694498997`*^9, 
  3.579379622703681*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"g", "[", 
  RowBox[{"200000", ",", "\[CurlyPi]"}], "]"}]], "Input",
 CellChangeTimes->{{3.5792863791675777`*^9, 3.57928637923061*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{
   "-", "0.0006281697076617193803699093879605054391141027068103423644693790574\
7669005373`50."}], ",", 
   RowBox[{"{", 
    RowBox[{
    "l0", "\[Rule]", 
     "223.10679774997897862140677917182108927179020606973332801160785139472292\
35149187`50."}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.579286384049242*^9, 3.579379625086128*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"BinarySearch", "[", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"g", "[", 
       RowBox[{"#", ",", "\[CurlyPi]"}], "]"}], "\[LeftDoubleBracket]", "1", 
      "\[RightDoubleBracket]"}], ">", "0"}], "&"}], ",", "200000", ",", 
   "400000"}], "]"}]], "Input",
 CellChangeTimes->{{3.579286389308769*^9, 3.579286418064632*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"200000", ",", "400000"}], "}"}]], "Print",
 CellChangeTimes->{3.5793796340315313`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"200000", ",", "300000"}], "}"}]], "Print",
 CellChangeTimes->{3.5793796362181683`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"250000", ",", "300000"}], "}"}]], "Print",
 CellChangeTimes->{3.579379638592743*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"250000", ",", "275000"}], "}"}]], "Print",
 CellChangeTimes->{3.579379641514044*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"262500", ",", "275000"}], "}"}]], "Print",
 CellChangeTimes->{3.579379644252029*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"268750", ",", "275000"}], "}"}]], "Print",
 CellChangeTimes->{3.579379646361641*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"268750", ",", "271875"}], "}"}]], "Print",
 CellChangeTimes->{3.5793796487588673`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"270312", ",", "271875"}], "}"}]], "Print",
 CellChangeTimes->{3.57937965111607*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"270312", ",", "271093"}], "}"}]], "Print",
 CellChangeTimes->{3.579379653753745*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"270702", ",", "271093"}], "}"}]], "Print",
 CellChangeTimes->{3.579379655981289*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"270897", ",", "271093"}], "}"}]], "Print",
 CellChangeTimes->{3.5793796583318777`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"270897", ",", "270995"}], "}"}]], "Print",
 CellChangeTimes->{3.579379660739868*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"270897", ",", "270946"}], "}"}]], "Print",
 CellChangeTimes->{3.579379663019311*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"270897", ",", "270921"}], "}"}]], "Print",
 CellChangeTimes->{3.5793796654241133`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"270909", ",", "270921"}], "}"}]], "Print",
 CellChangeTimes->{3.579379667790627*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"270915", ",", "270921"}], "}"}]], "Print",
 CellChangeTimes->{3.579379670229705*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"270918", ",", "270921"}], "}"}]], "Print",
 CellChangeTimes->{3.579379672586233*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"270919", ",", "270921"}], "}"}]], "Print",
 CellChangeTimes->{3.579379675412138*^9}]
}, Open  ]],

Cell[BoxData["270920"], "Output",
 CellChangeTimes->{3.579379677766898*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
\[CurlyPi] = 1/224; (* Conjectured feasible by v08ltu at \
http://polymathprojects.org/2013/06/04/polymath-proposal-bounded-gaps-between-\
primes/#comment-19961 *)\
\>", "Subsection",
 CellChangeTimes->{{3.579039757775041*^9, 3.5790397592201433`*^9}, {
  3.5790416585652437`*^9, 3.5790416627056007`*^9}, {3.5793795333589163`*^9, 
  3.579379545296225*^9}, {3.579379697473724*^9, 3.579379707714003*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"\[CurlyPi]", "=", 
   RowBox[{"1", "/", "224"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.579039757775041*^9, 3.5790397592201433`*^9}, {
  3.5790416585652437`*^9, 3.5790416627056007`*^9}, {3.579379538064554*^9, 
  3.579379547849513*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"g", "[", 
  RowBox[{"100000", ",", "\[CurlyPi]"}], "]"}]], "Input",
 CellChangeTimes->{{3.5792862847573423`*^9, 3.579286289089642*^9}, {
  3.57937965063918*^9, 3.579379651621005*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "0.0114500565203569149889806262553175809899701749231498782607027008341371295\
1683`50.", ",", 
   RowBox[{"{", 
    RowBox[{
    "l0", "\[Rule]", 
     "157.61388300841896683054138842096061931477984715242570304086047510091712\
190400327`50."}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.579286300856514*^9, 3.5792863694498997`*^9, 
  3.579379622703681*^9, 3.579379754063883*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"g", "[", 
  RowBox[{"10000", ",", "\[CurlyPi]"}], "]"}]], "Input",
 CellChangeTimes->{{3.5792863791675777`*^9, 3.57928637923061*^9}, 
   3.579379654477789*^9, {3.5793797740591307`*^9, 3.579379793557247*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{
   "-", "0.0021986639965279409797628072642093066359255388217321437614590057184\
0389514976`50."}], ",", 
   RowBox[{"{", 
    RowBox[{
    "l0", "\[Rule]", 
     "49.500000000001115251308308776823641057547330144040666900059461327732845\
27322259`50."}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.579286384049242*^9, 3.579379625086128*^9, {3.579379756558496*^9, 
   3.579379794626955*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"BinarySearch", "[", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"g", "[", 
       RowBox[{"#", ",", "\[CurlyPi]"}], "]"}], "\[LeftDoubleBracket]", "1", 
      "\[RightDoubleBracket]"}], ">", "0"}], "&"}], ",", "10000", ",", 
   "100000"}], "]"}]], "Input",
 CellChangeTimes->{{3.579379717701388*^9, 3.579379720612649*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"10000", ",", "100000"}], "}"}]], "Print",
 CellChangeTimes->{3.579379721098937*^9, 3.57937980210065*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"10000", ",", "55000"}], "}"}]], "Print",
 CellChangeTimes->{3.579379721098937*^9, 3.579379802447838*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"10000", ",", "32500"}], "}"}]], "Print",
 CellChangeTimes->{3.579379721098937*^9, 3.579379802833311*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"10000", ",", "21250"}], "}"}]], "Print",
 CellChangeTimes->{3.579379721098937*^9, 3.5793798032476807`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"10000", ",", "15625"}], "}"}]], "Print",
 CellChangeTimes->{3.579379721098937*^9, 3.5793798036424093`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"10000", ",", "12812"}], "}"}]], "Print",
 CellChangeTimes->{3.579379721098937*^9, 3.5793798040339947`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"11406", ",", "12812"}], "}"}]], "Print",
 CellChangeTimes->{3.579379721098937*^9, 3.5793798044130363`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"12109", ",", "12812"}], "}"}]], "Print",
 CellChangeTimes->{3.579379721098937*^9, 3.579379804800035*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"12460", ",", "12812"}], "}"}]], "Print",
 CellChangeTimes->{3.579379721098937*^9, 3.579379805166738*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"12636", ",", "12812"}], "}"}]], "Print",
 CellChangeTimes->{3.579379721098937*^9, 3.579379805547946*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"12636", ",", "12724"}], "}"}]], "Print",
 CellChangeTimes->{3.579379721098937*^9, 3.579379805919415*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"12636", ",", "12680"}], "}"}]], "Print",
 CellChangeTimes->{3.579379721098937*^9, 3.579379806288227*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"12636", ",", "12658"}], "}"}]], "Print",
 CellChangeTimes->{3.579379721098937*^9, 3.57937980666541*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"12647", ",", "12658"}], "}"}]], "Print",
 CellChangeTimes->{3.579379721098937*^9, 3.579379807047282*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"12652", ",", "12658"}], "}"}]], "Print",
 CellChangeTimes->{3.579379721098937*^9, 3.5793798074227333`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"12655", ",", "12658"}], "}"}]], "Print",
 CellChangeTimes->{3.579379721098937*^9, 3.579379807793285*^9}]
}, Open  ]],

Cell[BoxData["12656"], "Output",
 CellChangeTimes->{3.579379750325243*^9, 3.5793798081945477`*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Ignoring \[Kappa]", "Subsection",
 CellChangeTimes->{{3.5793806750266237`*^9, 3.57938067953253*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"\[Kappa]", "[", 
   RowBox[{"k0_", ",", "l0_", ",", "\[CurlyPi]_"}], "]"}], ":=", 
  "0"}]], "Input",
 CellChangeTimes->{{3.579379915158923*^9, 3.579379922927245*^9}, {
  3.57937995969657*^9, 3.579379961918908*^9}}],

Cell[BoxData[
 RowBox[{"k0", "=."}]], "Input",
 CellChangeTimes->{{3.57937998322309*^9, 3.579379986169091*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Solve", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"D", "[", 
     RowBox[{
      RowBox[{"f", "[", 
       RowBox[{"k0", ",", "l0", ",", "\[CurlyPi]0"}], "]"}], ",", "l0"}], 
     "]"}], "\[Equal]", "0"}], ",", "l0"}], "]"}]], "Input",
 CellChangeTimes->{{3.5793799637541237`*^9, 3.579380013262562*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"l0", "\[Rule]", 
     RowBox[{
      FractionBox["1", "2"], " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "1"}], "-", 
        SqrtBox["k0"]}], ")"}]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"l0", "\[Rule]", 
     RowBox[{
      FractionBox["1", "2"], " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "1"}], "+", 
        SqrtBox["k0"]}], ")"}]}]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.579379973658626*^9, 3.579380013673251*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Solve", "[", 
  RowBox[{
   RowBox[{"f", "[", 
    RowBox[{"k0", ",", 
     RowBox[{
      FractionBox["1", "2"], " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "1"}], "+", 
        SqrtBox["k0"]}], ")"}]}], ",", "\[CurlyPi]0"}], "]"}], "\[Equal]", 
   "0"}], "]"}]], "Input",
 CellChangeTimes->{{3.579380025423558*^9, 3.579380039352467*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Solve", "::", "\<\"svars\"\>"}], 
  RowBox[{
  ":", " "}], "\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"Equations may not give \
solutions for all \\\\\\\\\\\\\\\"solve\\\\\\\\\\\\\\\" \
variables.\\\\\\\"\\\", \\\"MT\\\"]\\) \
\\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", \
ButtonFrame->None, ButtonData:>\\\"paclet:ref/message/Solve/svars\\\", \
ButtonNote -> \\\"Solve::svars\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{3.579380040019038*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"k0", "\[Rule]", 
     FractionBox[
      RowBox[{"1", "+", 
       RowBox[{"2", " ", "\[CurlyPi]0"}], "-", 
       SqrtBox[
        RowBox[{"1", "+", 
         RowBox[{"4", " ", "\[CurlyPi]0"}]}]]}], 
      RowBox[{"8", " ", 
       SuperscriptBox["\[CurlyPi]0", "2"]}]]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"k0", "\[Rule]", 
     FractionBox[
      RowBox[{"1", "+", 
       RowBox[{"2", " ", "\[CurlyPi]0"}], "+", 
       SqrtBox[
        RowBox[{"1", "+", 
         RowBox[{"4", " ", "\[CurlyPi]0"}]}]]}], 
      RowBox[{"8", " ", 
       SuperscriptBox["\[CurlyPi]0", "2"]}]]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.579380026360524*^9, 3.579380040020384*^9}}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"k", "[", "\[CurlyPi]_", "]"}], ":=", 
  RowBox[{"Ceiling", "[", 
   FractionBox[
    RowBox[{"1", "+", 
     RowBox[{"2", " ", "\[CurlyPi]"}], "+", 
     SqrtBox[
      RowBox[{"1", "+", 
       RowBox[{"4", " ", "\[CurlyPi]"}]}]]}], 
    RowBox[{"8", " ", 
     SuperscriptBox["\[CurlyPi]", "2"]}]], "]"}]}]], "Input",
 CellChangeTimes->{{3.579380699242011*^9, 3.579380712920763*^9}, {
  3.5793807434935093`*^9, 3.579380745868534*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"k", "[", "\[CurlyPi]", "]"}]], "Input",
 CellChangeTimes->{{3.579380725044923*^9, 3.57938072623335*^9}}],

Cell[BoxData["12656"], "Output",
 CellChangeTimes->{{3.5793807266443367`*^9, 3.5793807475703583`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"k", "[", 
  RowBox[{"1", "/", "42"}], "]"}]], "Input",
 CellChangeTimes->{{3.5793807525443*^9, 3.5793807538924227`*^9}}],

Cell[BoxData["462"], "Output",
 CellChangeTimes->{3.5793807544051123`*^9}]
}, Open  ]]
}, Open  ]]
}, Open  ]]
},
WindowSize->{1170, 920},
WindowMargins->{{Automatic, -1229}, {Automatic, 71}},
FrontEndVersion->"7.0 for Mac OS X x86 (32-bit) (February 18, 2009)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[545, 20, 224, 5, 27, "Input"],
Cell[772, 27, 1110, 32, 43, "Input"],
Cell[CellGroupData[{
Cell[1907, 63, 109, 1, 67, "Section"],
Cell[2019, 66, 2831, 85, 109, "Input"],
Cell[4853, 153, 764, 19, 49, "Input"],
Cell[5620, 174, 1309, 37, 96, "Input"],
Cell[6932, 213, 3641, 105, 153, "Input"],
Cell[10576, 320, 187, 5, 30, "Input"],
Cell[CellGroupData[{
Cell[10788, 329, 439, 12, 27, "Input"],
Cell[11230, 343, 6825, 93, 283, "Output"]
}, Open  ]],
Cell[18070, 439, 1110, 32, 43, "Input"],
Cell[CellGroupData[{
Cell[19205, 475, 783, 21, 27, "Input"],
Cell[CellGroupData[{
Cell[20013, 500, 218, 4, 20, "Print"],
Cell[20234, 506, 217, 4, 20, "Print"],
Cell[20454, 512, 223, 4, 20, "Print"],
Cell[20680, 518, 223, 4, 20, "Print"],
Cell[20906, 524, 224, 4, 20, "Print"],
Cell[21133, 530, 224, 4, 20, "Print"],
Cell[21360, 536, 222, 4, 20, "Print"],
Cell[21585, 542, 224, 4, 20, "Print"],
Cell[21812, 548, 222, 4, 20, "Print"],
Cell[22037, 554, 224, 4, 20, "Print"],
Cell[22264, 560, 222, 4, 20, "Print"],
Cell[22489, 566, 222, 4, 20, "Print"],
Cell[22714, 572, 224, 4, 20, "Print"],
Cell[22941, 578, 224, 4, 20, "Print"],
Cell[23168, 584, 222, 4, 20, "Print"],
Cell[23393, 590, 222, 4, 20, "Print"],
Cell[23618, 596, 222, 4, 20, "Print"],
Cell[23843, 602, 222, 4, 20, "Print"],
Cell[24068, 608, 222, 4, 20, "Print"],
Cell[24293, 614, 222, 4, 20, "Print"],
Cell[24518, 620, 224, 4, 20, "Print"]
}, Open  ]],
Cell[24757, 627, 173, 2, 27, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[24967, 634, 458, 12, 27, "Input"],
Cell[25428, 648, 143, 3, 27, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[25608, 656, 351, 6, 27, "Input"],
Cell[25962, 664, 77, 1, 27, "Output"]
}, Open  ]],
Cell[26054, 668, 339, 8, 27, "Input"],
Cell[CellGroupData[{
Cell[26418, 680, 252, 6, 27, "Input"],
Cell[26673, 688, 100, 1, 27, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[26810, 694, 486, 12, 27, "Input"],
Cell[27299, 708, 124, 2, 27, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[27460, 715, 131, 2, 27, "Input"],
Cell[27594, 719, 74, 1, 27, "Output"]
}, Open  ]]
}, Closed]],
Cell[CellGroupData[{
Cell[27717, 726, 109, 1, 37, "Section"],
Cell[27829, 729, 1050, 32, 50, "Input"],
Cell[28882, 763, 733, 24, 58, "Input"],
Cell[29618, 789, 981, 26, 27, "Input"],
Cell[CellGroupData[{
Cell[30624, 819, 163, 2, 34, "Subsection"],
Cell[30790, 823, 224, 5, 27, "Input"],
Cell[CellGroupData[{
Cell[31039, 832, 257, 5, 27, "Input"],
Cell[31299, 839, 527, 12, 27, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[31863, 856, 159, 3, 27, "Input"],
Cell[32025, 861, 379, 10, 27, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[32441, 876, 158, 3, 27, "Input"],
Cell[32602, 881, 375, 11, 27, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[33014, 897, 368, 10, 27, "Input"],
Cell[CellGroupData[{
Cell[33407, 911, 124, 3, 20, "Print"],
Cell[33534, 916, 124, 3, 20, "Print"],
Cell[33661, 921, 126, 3, 20, "Print"],
Cell[33790, 926, 124, 3, 20, "Print"],
Cell[33917, 931, 124, 3, 20, "Print"],
Cell[34044, 936, 124, 3, 20, "Print"],
Cell[34171, 941, 126, 3, 20, "Print"],
Cell[34300, 946, 124, 3, 20, "Print"],
Cell[34427, 951, 124, 3, 20, "Print"],
Cell[34554, 956, 124, 3, 20, "Print"],
Cell[34681, 961, 123, 3, 20, "Print"],
Cell[34807, 966, 126, 3, 20, "Print"],
Cell[34936, 971, 124, 3, 20, "Print"],
Cell[35063, 976, 124, 3, 20, "Print"],
Cell[35190, 981, 124, 3, 20, "Print"],
Cell[35317, 986, 124, 3, 20, "Print"],
Cell[35444, 991, 124, 3, 20, "Print"],
Cell[35571, 996, 124, 3, 20, "Print"]
}, Open  ]],
Cell[35710, 1002, 75, 1, 27, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[35822, 1008, 152, 3, 27, "Input"],
Cell[35977, 1013, 369, 10, 33, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[36383, 1028, 170, 3, 27, "Input"],
Cell[36556, 1033, 135, 2, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[36728, 1040, 171, 3, 27, "Input"],
Cell[36902, 1045, 116, 1, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[37055, 1051, 223, 5, 47, "Input"],
Cell[37281, 1058, 79, 1, 27, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[37397, 1064, 229, 6, 27, "Input"],
Cell[37629, 1072, 112, 1, 30, "Output"]
}, Open  ]],
Cell[37756, 1076, 400, 10, 27, "Input"],
Cell[CellGroupData[{
Cell[38181, 1090, 212, 5, 27, "Input"],
Cell[38396, 1097, 76, 1, 27, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[38521, 1104, 381, 6, 35, "Subsection"],
Cell[38905, 1112, 275, 6, 27, "Input"],
Cell[CellGroupData[{
Cell[39205, 1122, 159, 3, 27, "Input"],
Cell[39367, 1127, 404, 11, 27, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[39808, 1143, 158, 3, 27, "Input"],
Cell[39969, 1148, 396, 11, 27, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[40402, 1164, 368, 10, 27, "Input"],
Cell[CellGroupData[{
Cell[40795, 1178, 126, 3, 20, "Print"],
Cell[40924, 1183, 126, 3, 20, "Print"],
Cell[41053, 1188, 124, 3, 20, "Print"],
Cell[41180, 1193, 124, 3, 20, "Print"],
Cell[41307, 1198, 124, 3, 20, "Print"],
Cell[41434, 1203, 124, 3, 20, "Print"],
Cell[41561, 1208, 126, 3, 20, "Print"],
Cell[41690, 1213, 123, 3, 20, "Print"],
Cell[41816, 1218, 124, 3, 20, "Print"],
Cell[41943, 1223, 124, 3, 20, "Print"],
Cell[42070, 1228, 126, 3, 20, "Print"],
Cell[42199, 1233, 124, 3, 20, "Print"],
Cell[42326, 1238, 124, 3, 20, "Print"],
Cell[42453, 1243, 126, 3, 20, "Print"],
Cell[42582, 1248, 124, 3, 20, "Print"],
Cell[42709, 1253, 124, 3, 20, "Print"],
Cell[42836, 1258, 124, 3, 20, "Print"],
Cell[42963, 1263, 124, 3, 20, "Print"]
}, Open  ]],
Cell[43102, 1269, 75, 1, 27, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[43226, 1276, 411, 7, 35, "Subsection"],
Cell[43640, 1285, 272, 6, 27, "Input"],
Cell[CellGroupData[{
Cell[43937, 1295, 207, 4, 27, "Input"],
Cell[44147, 1301, 426, 11, 27, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[44610, 1317, 231, 4, 27, "Input"],
Cell[44844, 1323, 449, 13, 27, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[45330, 1341, 367, 10, 27, "Input"],
Cell[CellGroupData[{
Cell[45722, 1355, 144, 3, 20, "Print"],
Cell[45869, 1360, 144, 3, 20, "Print"],
Cell[46016, 1365, 144, 3, 20, "Print"],
Cell[46163, 1370, 146, 3, 20, "Print"],
Cell[46312, 1375, 146, 3, 20, "Print"],
Cell[46461, 1380, 146, 3, 20, "Print"],
Cell[46610, 1385, 146, 3, 20, "Print"],
Cell[46759, 1390, 144, 3, 20, "Print"],
Cell[46906, 1395, 144, 3, 20, "Print"],
Cell[47053, 1400, 144, 3, 20, "Print"],
Cell[47200, 1405, 144, 3, 20, "Print"],
Cell[47347, 1410, 144, 3, 20, "Print"],
Cell[47494, 1415, 143, 3, 20, "Print"],
Cell[47640, 1420, 144, 3, 20, "Print"],
Cell[47787, 1425, 146, 3, 20, "Print"],
Cell[47936, 1430, 144, 3, 20, "Print"]
}, Open  ]],
Cell[48095, 1436, 98, 1, 27, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[48242, 1443, 106, 1, 35, "Subsection"],
Cell[48351, 1446, 250, 6, 27, "Input"],
Cell[48604, 1454, 111, 2, 27, "Input"],
Cell[CellGroupData[{
Cell[48740, 1460, 326, 9, 27, "Input"],
Cell[49069, 1471, 549, 19, 46, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[49655, 1495, 382, 12, 46, "Input"],
Cell[50040, 1509, 507, 10, 23, "Message"],
Cell[50550, 1521, 753, 23, 53, "Output"]
}, Open  ]],
Cell[51318, 1547, 471, 13, 53, "Input"],
Cell[CellGroupData[{
Cell[51814, 1564, 128, 2, 27, "Input"],
Cell[51945, 1568, 102, 1, 27, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[52084, 1574, 144, 3, 27, "Input"],
Cell[52231, 1579, 74, 1, 27, "Output"]
}, Open  ]]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)