|
|
A083343
|
|
Decimal expansion of constant B3 (or B_3) related to the Mertens constant.
|
|
12
|
|
|
1, 3, 3, 2, 5, 8, 2, 2, 7, 5, 7, 3, 3, 2, 2, 0, 8, 8, 1, 7, 6, 5, 8, 2, 8, 7, 7, 6, 0, 7, 1, 0, 2, 7, 7, 4, 8, 8, 3, 8, 4, 5, 9, 4, 8, 9, 0, 4, 2, 4, 2, 2, 6, 6, 1, 7, 8, 7, 1, 3, 0, 8, 9, 9, 7, 5, 7, 3, 4, 0, 0, 4, 1, 7, 1, 9, 3, 0, 4, 0, 1, 8, 6, 8, 7, 5, 4, 8, 0, 4, 5, 5, 1, 4, 1, 6, 8, 6, 2
(list;
constant;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
Comment from David Broadhurst Feb 26 2014, via a posting to the Number Theory Mailing List, and included here with permission. This comment concerns the number C = 1 + B_3 = 2.33258227573322088... (see also A238114). (Start)
In a beautifully concise and clear paper, "Newman's short proof of the prime number theorem", Don Zagier condensed work by Euler, Riemann, Chebyshev, de la Vallee Poussin, Hadamard, Mertens and, more recently, D. J. Newman (in 1980), to achieve a short self-contained proof of the prime number theorem that is within the reach of a reader who understands enough real analysis to apply the mean-value theorem and enough complex analysis to apply Cauchy's theorem.
The heart of this proof is the convergence of the integral
C = Integral_{x = 1..oo} (x - theta(x))/x^2 dx ... [1]
where theta(x) = Sum_{prime p <= x} log(p) is the sum of the natural logs of all the primes not exceeding x. Then the Prime Number Theorem follows in the form theta(x) ~ x, for otherwise the integral in [1] would not converge.
In a sense, this constant C is rather significant: if it did not exist the proof would fail. However, its actual value is a matter of sublime indifference to a true mathematician. To prove that it exists, one may use the equivalent expression
C = 1 + Euler + Sum_{prime p} log(p)/(p^2-p) ... [2]
that follows from Zagier's account. Here the sum is over all the positive primes and clearly converges, since the corresponding sum over integers n > 1 converges.
It is also easy, if unnecessary, to show that
C = 1 + Euler + Sum_{s > 1} mu(s)*zeta'(s)/zeta(s) ... [3]
where mu(s) is the Moebius function. An approximate evaluation of this formula requires the derivatives zeta'(s) of Riemann's zeta(s) = Sum_{n > 0} 1/n^s at sufficiently many squarefree integers s > 1.
By use of both [2] and [3], J. Barkley Rosser and Lowell Schoenfeld obtained (effectively) 16 good digits of C in "Approximate formulas for some functions of prime numbers", where they gave, in (2.11), a numerical result for 1 - C.
A better way to compute C, however, is by use of a method indicated in Henri Cohen's paper "High precision computation of Hardy-Littlewood constants".
(End)
|
|
REFERENCES
|
Henri Cohen, Number Theory, Volume II: Analytic and Modern Tools, GTM Vol. 240, Springer, 2007; see pp. 208-209.
S. R. Finch, Mathematical Constants, Encyclopedia of Mathematics and its Applications, vol. 94, Cambridge University Press, pp. 94-98.
Edmund Landau, Handbuch der Lehre von der Verteilung der Primzahlen, 2nd ed., Chelsea, 1953, pp. 197-203.
József Sándor, Dragoslav S. Mitrinovic and Borislav Crstici, Handbook of Number Theory I, Springer Science & Business Media, 2005, Chapter VI, p. 199.
|
|
LINKS
|
Steven R. Finch, Mathematical Constants II, Encyclopedia of Mathematics and Its Applications, Cambridge University Press, Cambridge, 2018, p. 196.
|
|
FORMULA
|
Equals lim_{x->oo} (log x - Sum_{p <= x} log(p)/p). - Dick Boland, Mar 09 2008
Equals EulerGamma - Sum_{n >= 2} P'(n), where P'(n) is the prime zeta P function derivative. - Jean-François Alcover, Apr 25 2016
Equals lim_{n->oo} Sum_{k=1..n} mu(k)^2/phi(k) - log(n) (Ward, 1927). - Amiram Eldar, Mar 05 2021
|
|
EXAMPLE
|
1.3325822757332208817658287760710277488384594890424226617871308997573400417193...
|
|
MATHEMATICA
|
digits = 99; B3 = EulerGamma - NSum[PrimeZetaP'[n], {n, 2, 10^4}, WorkingPrecision -> 2 digits, NSumTerms -> 200]; RealDigits[B3, 10, digits][[1]] (* Jean-François Alcover, Apr 25 2016 *)
|
|
CROSSREFS
|
|
|
KEYWORD
|
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|