[math-fun] Simple formula for "Roth's constant"
(A014578 <http://oeis.org/A014578>) Out[665]= k (-1) Sum[--------, {k, 0, ∞}] k 3 2 - 1 In[666]:= RealDigits[% /. ∞] -> 4, 8] Out[666]= {{{6, 6, 7, 6, 6, 7, 6, 6, 6, 6, 6, 7, 6, 6, 7, 6, 6, 6, 6, 6, 7, 6, 6, 7, 6, 6, 7}}, 0} --Bill (As computed by the Life pattern x = 50, y = 50, rule = B3/S23 26b3o11b3o$26bo2bo10bo2bo$26bo6b3o4bo$26bo5bo2bo4bo$27bo4bo3bo4bo$33b 4o$bo8b2o23bo$o8b4o$o3bo3b2ob2o28bo$4o5b2o29b3o$40bob2o$41b3o$41b2o$4b 2o$2b2ob2o$2bo2bo$2bo2bo$3b2o$24b3o$24bo2bo$bo22bo$o23bo$o3bo14bo5bo$ 4o14bo$18bo3bo5b4o$18b4o5bo3bo$31bo14b4o$24bo5bo14bo3bo$25bo23bo$25bo 22bo$22bo2bo$23b3o$45b2o$44bo2bo$44bo2bo$43b2ob2o$44b2o$7b2o$6b3o$6b2o bo$7b3o29b2o5b4o$8bo28b2ob2o3bo3bo$37b4o8bo$14bo23b2o8bo$13b4o$8bo4bo 3bo4bo$9bo4bo2bo5bo$9bo4b3o6bo$6bo2bo10bo2bo$7b3o11b3o! misnamed thue.)
What the heck is this ⬇? I can't delete it.
I'll edit the sequence. If you do not see the edit within 10 days, remind me. Note the L-system I gave for the binary expansion, giving a divisionless method for computation. Best, jj * Bill Gosper <billgosper@gmail.com> [Apr 16. 2014 20:05]:
(A014578 <http://oeis.org/A014578>) Out[665]= k (-1) Sum[--------, {k, 0, ∞}] k 3 2 - 1
In[666]:= RealDigits[% /. ∞] -> 4, 8]
Out[666]= {{{6, 6, 7, 6, 6, 7, 6, 6, 6, 6, 6, 7, 6, 6, 7, 6, 6, 6, 6, 6, 7, 6, 6, 7, 6, 6, 7}}, 0} --Bill (As computed by the Life pattern x = 50, y = 50, rule = B3/S23 26b3o11b3o$26bo2bo10bo2bo$26bo6b3o4bo$26bo5bo2bo4bo$27bo4bo3bo4bo$33b 4o$bo8b2o23bo$o8b4o$o3bo [wow, great intro!] 3b2ob2o28bo$4o5b2o29b3o$40bob2o$41b3o$41b2o$4b 2o$2b2ob2o$2bo2bo$2bo2bo$3b2o$24b3o$24bo2bo$bo22bo$o23bo$o3bo14bo5bo$ 4o14bo$18bo3bo5 [...amazing development... ] b4o$18b4o5bo3bo$31bo14b4o$24bo5bo14bo3bo$25bo23bo$25bo 22bo$22bo2bo$23b3o$45b2o$44bo2bo$44bo2bo$43b2ob2o$44b2o$7b2o$6b3o$6b2o bo$7b3o29b2o5b4o$8b [but this part could stand polishing; it's not clear what b02bo is doing at all] o28b2ob2o3bo3bo$37b4o8bo$14bo23b2o8bo$13b4o$8bo4bo 3bo4bo$9bo4bo2bo5bo$9bo4b3o6bo$6bo2bo10bo2bo$7b3o11b3o! [but that grand finale saves everything!] misnamed thue.)
What the heck is this ⬇? I can't delete it.
Char: ⬇ (11015, #o25407, #x2b07, file ...) point=1577 of 1803 (87%) column=28 Doesn't render as character in emacs (version 24.3.4), invalid UTF-8?
_______________________________________________ math-fun mailing list math-fun@mailman.xmission.com http://mailman.xmission.com/cgi-bin/mailman/listinfo/math-fun
On Tue, Apr 15, 2014 at 3:13 PM, Bill Gosper <billgosper@gmail.com> wrote:
(A014578 <http://oeis.org/A014578>) Out[665]= k (-1) Sum[--------, {k, 0, ∞}] k 3 2 - 1
In[666]:= RealDigits[% /. ∞-> 4, 8]
Out[666]= {{{6, 6, 7, 6, 6, 7, 6, 6, 6, 6, 6, 7, 6, 6, 7, 6, 6, 6, 6, 6, 7, 6, 6, 7, 6, 6, 7}}, 0}
Note that this and the other partial sums are actual exemplars of Roth's criterion. In[685]:= FromDigits[%666, 8]
Out[685]= 2077175852086666748456887/2417851639229258349412351 In[686]:= $MaxExtraPrecision = 9999; N[Round[1/(% - %665)], 239] N::meprec: Internal precision limit $MaxExtraPrecision = 9999.` reached while evaluating Round[1/((2077175852086666748456887/2417851639229258349412351)-\!\(\*UnderoverscriptBox[\(\[Sum]\), \(k = 0\), \(\[Infinity]\)]\*FractionBox[SuperscriptBox[\((\(-1\))\), \(k\)], \(\(-1\) + Power[<<1>>]\)]\))]. >> (I don't understand this warning. 239 and 9999 are ridiculous overkill. Varying them (above a safe minimum) has no effect on the answer or the warning.) Out[686]= 14134776518227074636666380005943348126619871175004951664972849610340958207 In[697]:= Denominator[%685]^3 Out[697]= \ 14134776518227074636666362467923700155784852730794212624494301514198679551 where all that is needed is Denominator^(2+ϵ) (infinitely often). (As computed by the Life pattern
x = 50, y = 50, rule = B3/S23
[...]) Tighter: #CXRLE #C Thue (really Roth) generator #C #C Four puffers produce gliders, whose collisions produce four #C lines of blinkers, each representing the transcendental number #C .110110111110110111110110110..., where the n'th bit is the #C the mod 2 of 1 + the largest power of 3 dividing n. #C Built by Bill Gosper, in or before 1987 x = 50, y = 50, rule = B3/S23 26b3o11b3o$26bo2bo10bo2bo$26bo6b3o4bo$26bo5bo2bo4bo$27bo7bo5bo2$bo8b2o $o8b4o$o3bo3b2ob2o28bo$4o5b2o29b3o$40bob2o$41b3o$41b2o2$2b3o$2bo$2bo$ 3bo$24b3o$24bo2bo$bo22bo$o23bo$o3bo14bo5bo$4o14bo$18bo3bo5b4o$18b4o5bo 3bo$31bo14b4o$24bo5bo14bo3bo$25bo23bo$25bo22bo$22bo2bo$23b3o$46bo$47bo $47bo$45b3o2$7b2o$6b3o$6b2obo$7b3o29b2o5b4o$8bo28b2ob2o3bo3bo$37b4o8bo $38b2o8bo2$8bo5bo7bo$9bo4bo2bo5bo$9bo4b3o6bo$6bo2bo10bo2bo$7b3o11b3o! --rwg
participants (2)
-
Bill Gosper -
Joerg Arndt