%I #33 May 10 2023 10:53:07
%S 1,8,7,9,8,5,3,8,6,2,1,7,5,2,5,8,5,3,3,4,8,6,3,0,6,1,4,5,0,7,0,9,6,0,
%T 0,3,8,8,1,9,8,7,3,4,0,0,4,8,9,2,8,9,9,0,4,8,2,9,6,1,7,6,6,9,1,2,2,2,
%U 9,6,3,8,6,6,6,1,2,1,4,2,1,1,3,6,1,7,6,5,0,1,9,7,3,8,9,1,2,3,5,3,2,3,9,6,8
%N Decimal expansion of Sum_{n>=1} (n!/n^n).
%C Sum_{n>=1} n!/n^(n+2) = Integral_{x=0..infinity} -log(1-x*exp(-x)) dx = 1.157694752682... . - _Vaclav Kotesovec_, Jan 05 2016
%F Equals Integral_{x=0..oo} x*exp(x)/(exp(x)-x)^2 dx. - _Michael Somos_, May 18 2021
%F Equals Integral_{x=1..oo} 1/(x - log(x))^2 dx. - _Fabián Pereyra_, May 10 2023
%e 1.879853862175258533486306145...
%t RealDigits[ NSum[n!/n^n, {n, 1, Infinity}, WorkingPrecision -> 110, NSumTerms -> 180] , 10, 99] // First (* _Jean-François Alcover_, Feb 12 2013 *)
%t RealDigits @ NIntegrate[x E^x/(E^x - x)^2, {x, 0, Infinity}, WorkingPrecision -> 105] // First (* _Michael Somos_, May 18 2021 *)
%o (PARI) firstDecimalDigits(n)={default(realprecision,n);return(digits(floor(suminf(n=1,n!/(n^n))*10^n)))}
%o print(firstDecimalDigits(98)); \\ _R. J. Cano_, Dec 30 2016
%o (PARI) suminf(n=1,n!/(n^n)) \\ _Michel Marcus_, Dec 30 2016
%Y Cf. A055775, A073225.
%K nonn,cons
%O 1,2
%A _Ross La Haye_, May 01 2004
%E More terms from Francisco Salinas (franciscodesalinas(AT)hotmail.com), May 01 2004