|
REFERENCES
|
Steven R. Finch, Mathematical Constants, Cambridge University Press, 2003, Sections 2.17.1 Ruelle-Mayer Operators, p. 153.
|
|
MATHEMATICA
|
digits = 16; m0 = 100; dm = 5; Clear[f];
f[m_] := f[m] = NIntegrate[BesselJ[1, 2*Sqrt[u*v]]^2/((Exp[u]-1) * (Exp[v]-1)), {u, 0, m}, {v, 0, m - u}, MaxRecursion -> 30, WorkingPrecision -> digits + 10]; f[m = m0]; Print[m, " ", RealDigits[f[m], 10, digits][[1]]]; f[m = m0 + dm]; Print[m, " ", RealDigits[f[m], 10, digits][[1]]]; While[RealDigits[f[m], 10, digits][[1]] != RealDigits[f[m - dm], 10, digits][[1]], m = m + dm; Print[m, " ", RealDigits[f[m], 10, digits][[1]]]]; RealDigits[f[m], 10, digits][[1]]
|