OFFSET
0,1
COMMENTS
The convergence is very slow, need to use the first 100000000 primes to obtain the correct value of the coefficient of 10^(-23).
The constant is in the interval [0.28417070547086825017714, 0.28417070547086825017743]; these safe limits are computed by accumulating in parallel the partial sum of the lower estimate 1/n^4 = Zeta(4). - R. J. Mathar, Feb 06 2015
EXAMPLE
0.284170705470...
PROG
(PFGW & SCRIPT)
SCRIPT
DIM i, 0
DIM j, 0
DIM n
DIM m
DIMS t
OPENFILEOUT myf, a(n).txt
OPENFILEIN maf, pre.txt
LABEL loop1
SET i, i+1
IF i>10000000 THEN END
GETNEXT n, maf
SET j, j+10^10000/((i*n)^2)
IF i%1000000==0 THEN SET m, j/10^9970
IF i%1000000==0 THEN WRITE myf, m
GOTO loop1
CROSSREFS
KEYWORD
AUTHOR
Pierre CAMI, Jan 07 2015
STATUS
approved