

A036792


Decimal expansion of Integral(sin(x)/x, x=0..Pi).


7



1, 8, 5, 1, 9, 3, 7, 0, 5, 1, 9, 8, 2, 4, 6, 6, 1, 7, 0, 3, 6, 1, 0, 5, 3, 3, 7, 0, 1, 5, 7, 9, 9, 1, 3, 6, 3, 3, 4, 5, 8, 0, 9, 7, 2, 8, 9, 8, 1, 1, 5, 4, 9, 0, 9, 8, 0, 4, 7, 8, 3, 7, 8, 1, 8, 7, 6, 9, 8, 1, 8, 9, 0, 1, 6, 6, 3, 4, 8, 3, 5, 8, 5, 3, 2, 7, 1, 0, 3, 3, 6, 5, 0, 2, 9, 5, 4, 7, 5, 7, 7, 0, 1, 6, 8
(list;
constant;
graph;
refs;
listen;
history;
text;
internal format)



OFFSET

1,2


COMMENTS

Also known as Gibbs Constant.
Integral(sin(x)/x dx) = x  x^3/(3*3!) + x^5/(5*5!)  x^7/(7*7!) + ... [From Harry J. Smith, May 01 2009]


LINKS

Harry J. Smith, Table of n, a(n) for n=1,...,20000
Simon Plouffe, Gibbs
Eric Weisstein's World of Mathematics, WilbrahamGibbs Constant


MATHEMATICA

RealDigits[ N[ SinIntegral[Pi], 110]] [[1]]


PROG

(PARI) { default(realprecision, 20080); y=0; x=Pi; m=x; x2=x*x; n=1; nf=1; s=1; while (x!=y, y=x; n++; nf*=n; n++; nf*=n; m*=x2; s=s; x+=s*m/(n*nf)); for (n=1, 20000, d=floor(x); x=(xd)*10; write("b036792.txt", n, " ", d)); } [From Harry J. Smith, May 01 2009]


CROSSREFS

Cf. A036790.
Sequence in context: A086730 A139721 A021058 * A198581 A154735 A010147
Adjacent sequences: A036789 A036790 A036791 * A036793 A036794 A036795


KEYWORD

nonn,cons


AUTHOR

N. J. A. Sloane.


EXTENSIONS

Fixed my PARI program, had n Harry J. Smith, May 19 2009


STATUS

approved



