|
|
A059534
|
|
Beatty sequence for 1 + 1/Catalan's constant.
|
|
3
|
|
|
2, 4, 6, 8, 10, 12, 14, 16, 18, 20, 23, 25, 27, 29, 31, 33, 35, 37, 39, 41, 43, 46, 48, 50, 52, 54, 56, 58, 60, 62, 64, 66, 69, 71, 73, 75, 77, 79, 81, 83, 85, 87, 89, 92, 94, 96, 98, 100, 102, 104, 106, 108, 110, 112, 115, 117, 119, 121, 123, 125, 127, 129, 131, 133
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
LINKS
|
|
|
PROG
|
(PARI) { digits=100; default(realprecision, digits+80); s=1.0; n=5*digits; j=4*n+1; si=-1.0; for (i=3, j-2, s+=si/i^2; si=-si; i++; ); s+=0.5/j^2; ttk=4.0; d=4.0*j^3; xk=2.0; xkp=xk; for (k=2, 100000000, term=(ttk-1)*ttk*xkp; xk++; xkp*=xk; if (k>2, term*=xk; xk++; xkp*=xk; ); term*=bernreal(k)/d; sn=s+term; if (sn==s, break); s=sn; ttk*=4.0; d*=(k+1)*(k+2)*j^2; k++; ); b=1 + 1/s; for (n = 1, 2000, write("b059534.txt", n, " ", floor(n*b)); ) } \\ Harry J. Smith, Jun 27 2009
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|