|
| |
|
|
A014571
|
|
Consider the Morse-Thue sequence (A010060) as defining a binary constant and convert it to decimal.
|
|
10
| |
|
|
4, 1, 2, 4, 5, 4, 0, 3, 3, 6, 4, 0, 1, 0, 7, 5, 9, 7, 7, 8, 3, 3, 6, 1, 3, 6, 8, 2, 5, 8, 4, 5, 5, 2, 8, 3, 0, 8, 9, 4, 7, 8, 3, 7, 4, 4, 5, 5, 7, 6, 9, 5, 5, 7, 5, 7, 3, 3, 7, 9, 4, 1, 5, 3, 4, 8, 7, 9, 3, 5, 9, 2, 3, 6, 5, 7, 8, 2, 5, 8, 8, 9, 6, 3, 8, 0, 4, 5, 4, 0, 4, 8, 6, 2, 1, 2, 1, 3, 3, 3, 9, 6, 2, 5, 6
(list; constant; graph; refs; listen; history; internal format)
|
|
|
|
OFFSET
| 0,1
|
|
|
COMMENTS
| The constant is Sum A010060(n)*2^(-n).
|
|
|
LINKS
| Harry J. Smith, Table of n, a(n) for n = 0..20000
Eric Weisstein's World of Mathematics, Thue-Morse Constant
|
|
|
EXAMPLE
| 0.412454033640107597783361368258455283089...
|
|
|
MAPLE
| A010060 := proc(n) add(i, i=convert(n, base, 2)) mod 2 ; end: A014571 := proc() local nlim, aold, a ; nlim := ilog2(10^Digits) ; aold := add( A010060(n)/2^n, n=0..nlim) ; a := 0.0 ; while abs(a-aold) > abs(a)/10^(Digits-3) do aold := a; nlim := nlim+200 ; a := add( A010060(n)/2^n, n=0..nlim) ; od: evalf(%/2) ; end: A014571() ; - R. J. Mathar (mathar(AT)strw.leidenuniv.nl), Mar 03 2008
|
|
|
PROG
| (PARI) { default(realprecision, 20080); x=0.0; m=67000; for (n=1, m-1, x=x+x; x=x+sum(k=0, length(binary(n))-1, bittest(n, k))%2); x=10*x/2^m; for (n=0, 20000, d=floor(x); x=(x-d)*10; write("b014571.txt", n, " ", d)); } [From Harry J. Smith (hjsmithh(AT)sbcglobal.net), Apr 25 2009]
|
|
|
CROSSREFS
| Cf. A010060, A058631.
Sequence in context: A063987 A010126 A021712 * A152523 A082903 A154589
Adjacent sequences: A014568 A014569 A014570 * A014572 A014573 A014574
|
|
|
KEYWORD
| nonn,cons
|
|
|
AUTHOR
| Eric Weisstein (eric(AT)weisstein.com)
|
|
|
EXTENSIONS
| Corrected and extended by R. J. Mathar (mathar(AT)strw.leidenuniv.nl), Mar 03 2008
Fixed my PARI program, had -n Harry J. Smith (hjsmithh(AT)sbcglobal.net), May 19 2009
|
| |
|
|
