login
Decimal expansion of sum 1/(p^2 * log p) over the primes p=2,3,5,7,11,...
10

%I #33 Sep 09 2024 13:43:17

%S 5,0,7,7,8,2,1,8,7,8,5,9,1,9,9,3,1,8,7,7,4,3,7,5,1,0,3,7,9,4,7,0,5,5,

%T 7,0,4,6,6,9,7,3,6,7,1,7,0,4,3,2,0,6,9,8,5,7,3,9,8,0,2,1,2,3,4,8,2,7,

%U 2,8,6,9,0,1,3,7,4,1,3,1,1,5,1,0,4,6,4,6,6,7,8,4,8,9,5,2,9,2,1,1,3,5,6,4,5,4

%N Decimal expansion of sum 1/(p^2 * log p) over the primes p=2,3,5,7,11,...

%D Henri Cohen, Number Theory, Volume II: Analytic and Modern Tools, GTM Vol. 240, Springer, 2007; see pp. 208-209.

%H Karim Belabas and Henri Cohen, <a href="/A137245/a137245.gp.txt">Computation of sum_{p prime} 1/(p^s log(p))</a>, PARI/GP script, 2020.

%H Henri Cohen, <a href="http://www.math.u-bordeaux.fr/~cohen/hardylw.dvi">High-precision calculation of Hardy-Littlewood constants</a>, (1998).

%H Henri Cohen, <a href="/A221712/a221712.pdf">High-precision computation of Hardy-Littlewood constants</a>. [pdf copy, with permission]

%H R. J. Mathar, <a href="https://arxiv.org/abs/0811.4739">Twenty digits of some integrals of the prime zeta function</a>, arXiv:0811.4739 [math.NT], 2009-2018, Section 2.4.

%e 0.50778218785919931877437510379470557...

%t digits = 106; precision = digits + 15;

%t tmax = 400; (* integrand considered negligible beyond tmax *)

%t kmax = 400; (* f(k) considered negligible beyond kmax *)

%t InLogZeta[k_] := NIntegrate[Log[Zeta[t]], {t, k, tmax},

%t WorkingPrecision -> precision, MaxRecursion -> 20,

%t AccuracyGoal -> precision];

%t f[k_] := With[{mu = MoebiusMu[k]}, If[mu == 0, 0, (mu/k^2)*InLogZeta[2k]]];

%t s = 0; Do[s = s + f[k]; Print[k, " ", s], {k, 1, kmax}];

%t RealDigits[s][[1]][[1 ;; digits]] (* _Jean-François Alcover_, Feb 06 2021, updated Jun 23 2022 *)

%o (PARI) \\ See Belabas, Cohen link. Run as SumEulerlog(2) after setting the required precision.

%o (PARI) default(realprecision, 200); s=0; for(k=1, 300, s = s + moebius(k)/k^2 * intnum(x=2*k,[[1], 1], log(zeta(x))); print(s)); \\ _Vaclav Kotesovec_, Jun 12 2022

%Y Cf. A137245.

%K nonn,cons

%O 0,1

%A _N. J. A. Sloane_, Jan 26 2013

%E More terms from _Hugo Pfoertner_, Feb 01 2020

%E More digits from _Vaclav Kotesovec_, Jun 12 2022