# Download https://numberdb.org/files/T685/generate.py?raw=1 as generate.py.
# Install numberdb with sage -pip install numberdb.
from generate import moment
moment(4, 2, digits=100)The three-term values use certified Arb integration of the compact log-cosine integral. The four-term values use a polynomial truncation of the convergent series for the conditional mean g, with a proved uniform tail bound propagated through its kth power. Real and complex ball arithmetic enclose argument rounding and every arithmetic operation. The 100-digit computation uses 445 working bits and 235 series terms; each result must retain at least 373 accurate bits. The minimum measured relative accuracy was 431 bits, about 129 decimal digits. The attached certification.md derives the bounds, and numerical-evidence.json records the actual radii.
All 24 entries are recomputed at 140 requested digits and compared with independent mpmath quadrature of the log-sine and Clausen forms at 75 and 125 decimal working digits. Known integral controls are run first. Five low-order formulas from Borwein and Straub provide additional checks using pi, zeta, Clausen and polylogarithm values. Positivity, the bound by log(2) between successive moments, and log-convexity are checked over every applicable order. The numerical quadrature is an independent agreement check; the ball integration and explicit series remainder justify the certified precision.