|
Search: id:A014571
|
|
|
| A014571 |
|
Consider the Morse-Thue sequence (A010060) as defining a binary constant and convert it to decimal. |
|
+0
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; cons; graph; listen)
|
|
|
OFFSET
|
0,1
|
|
|
COMMENT
|
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, Link to a section of The World of Mathematics.
|
|
EXAMPLE
|
0.4124540336401...
|
|
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
|
|
PROGRAM
|
(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 this_sequence A152523 A082903 A154589
Adjacent sequences: A014568 A014569 A014570 this_sequence 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
|
|
|
Search completed in 0.002 seconds
|
