Meissel–Mertens constant: Difference between revisions
Content added Content deleted
m (julia example) |
|||
Line 37: | Line 37: | ||
;See: |
;See: |
||
:* Details in the Wikipedia article: [https://en.wikipedia.org/wiki/Meissel%E2%80%93Mertens_constant Meissel–Mertens constant] |
:* Details in the Wikipedia article: [https://en.wikipedia.org/wiki/Meissel%E2%80%93Mertens_constant Meissel–Mertens constant] |
||
<br /><br /> |
|||
=={{header|Julia}}== |
|||
Off by one in the 11th digit after 10^8 primes. |
|||
<syntaxhighlight lang="julia"> |
|||
using Base.MathConstants # sets constant γ = 0.5772156649015... |
|||
using Primes |
|||
""" Approximate the Meissel-Mertons constant. """ |
|||
function meissel_mertens(iterations = 100_000_000) |
|||
return mapreduce(p ->(d = 1/p; log(1 - d) + d), +, primes(prime(iterations))) + γ |
|||
end |
|||
@show meissel_mertens(100_000_000) # meissel_mertens(100000000) = 0.2614972128591237 |
|||
</syntaxhighlight> |
|||