%I #22 Feb 16 2025 08:33:14
%S 7,8,5,3,8,0,5,5,7,2,9,8,6,3,2,8,7,3,4,9,2,5,8,3,0,1,1,4,6,7,3,3,2,5,
%T 2,4,7,6,1
%N Decimal expansion of Integral_{x = 0 to oo} Product_{m=1..oo} cos(x/m) dx.
%D J. M. Borwein, D. H. Bailey and R. Girgensohn, Experimentation in Mathematics, A K Peters, Ltd., Natick, MA, 2004. x+357 pp. See p. 101.
%D Steven R. Finch, Mathematical Constants, Encyclopedia of Mathematics and its Applications, vol. 94, Cambridge University Press, 2003, Section 6.3, p. 429.
%H Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/InfiniteCosineProductIntegral.html">Infinite Cosine Product Integral</a>.
%e 0.785380557298632873492583011467332524761...
%t (* This naive script is not suitable to get more than 10 digits. *)
%t digits = 10; m0 = 500; dm = 50; Clear[g];
%t f[n_Integer, m_] := NIntegrate[Product[Cos[x/k], {k, 1, m - 1}]*(((1 - x^2/(2*m^2))*Gamma[1 + m]^2)/(Gamma[1 + m - x/Sqrt[2]]*Gamma[1 + m + x/Sqrt[2]])), {x, n*Pi/2, (n + 1)*Pi/2}];
%t g[m_] := g[m] = NSum[f[n, m], {n, 0, Infinity}];
%t Print[g[m0]]; Print[g[m = m0 + dm]];
%t While[RealDigits[g[m], 10, digits][[1]] != RealDigits[g[m - dm], 10, digits][[1]], m = m + dm; Print[m, " ", RealDigits[g[m], 10, digits][[1]]]];
%t RealDigits[g[m], 10, digits][[1]] (* _Jean-François Alcover_, May 19 2016 *)
%K nonn,cons,more,changed
%O 0,1
%A _N. J. A. Sloane_, May 12 2011