[go: up one dir, main page]

login
Decimal expansion of Sum_{n >= 2} zeta(n)/n!.
6

%I #24 Mar 15 2017 04:02:24

%S 1,0,7,8,1,8,8,7,2,9,5,7,5,8,1,8,4,8,2,7,5,8,2,6,5,4,3,6,7,6,9,8,3,2,

%T 3,8,1,7,0,7,2,1,9,6,0,9,6,7,2,3,4,7,1,6,0,0,3,7,1,7,0,2,0,7,8,0,0,7,

%U 1,5,2,3,0,0,3,2,7,8,4,3,4,8,6,5,6,7,6,7,6,8,0,8,8,5,8,2,9,0,1

%N Decimal expansion of Sum_{n >= 2} zeta(n)/n!.

%H G. C. Greubel, <a href="/A093720/b093720.txt">Table of n, a(n) for n = 1..10000</a>

%H J. Sondow and E. W. Weisstein, <a href="http://mathworld.wolfram.com/RiemannZetaFunction.html">MathWorld: Riemann Zeta Function</a>

%F Equals Sum_{k>=1} (exp(1/k) - 1 - 1/k). - _Vaclav Kotesovec_, Mar 04 2016

%e 1.078188729575818482758265436769832381707219...

%p evalf(Sum(exp(1/n)-1-1/n, n=1..infinity), 120); # _Vaclav Kotesovec_, Mar 04 2016

%t digits = 99; ClearAll[z, rd]; z[k_] := z[k] = z[k-1] + N[Sum[Zeta[n]/n!, {n, 2^(k-1) + 1, 2^k}], digits]; z[0] = 0; rd[k_] := rd[k] = RealDigits[z[k]][[1]]; rd[0]; rd[k = 1]; While[ rd[k] != rd[k-1], k++]; rd[k] (* _Jean-François Alcover_, Nov 09 2012 *)

%o (PARI) suminf(n=2, zeta(n)/n!) \\ _Michel Marcus_, Mar 15 2017

%Y Cf. A076813, A093721, A269574, A269611, A269720, A269768.

%K nonn,cons

%O 1,3

%A _Eric W. Weisstein_, Apr 12 2004

%E Corrected by _Fredrik Johansson_, Mar 19 2006