%I #26 Oct 30 2022 15:05:49
%S 2,8,4,1,7,0,7,0,5,4,7,0,8,6,8,2,5,0,1,7,7,1,4
%N Decimal expansion of Sum_{n>=1} 1/A033286(n)^2.
%C The convergence is very slow, need to use the first 100000000 primes to obtain the correct value of the coefficient of 10^(-23).
%C 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
%e 0.284170705470...
%o (PFGW & SCRIPT)
%o SCRIPT
%o DIM i,0
%o DIM j,0
%o DIM n
%o DIM m
%o DIMS t
%o OPENFILEOUT myf,a(n).txt
%o OPENFILEIN maf,pre.txt
%o LABEL loop1
%o SET i,i+1
%o IF i>10000000 THEN END
%o GETNEXT n,maf
%o SET j,j+10^10000/((i*n)^2)
%o IF i%1000000==0 THEN SET m,j/10^9970
%o IF i%1000000==0 THEN WRITE myf,m
%o GOTO loop1
%Y Cf. A033286 (n*prime(n)), A124012.
%K nonn,cons,more
%O 0,1
%A _Pierre CAMI_, Jan 07 2015
|