RosettaCodeData/Task/Meissel-Mertens-constant/Python/meissel-mertens-constant.py
2023-07-01 13:44:08 -04:00

19 lines
358 B
Python

#!/usr/bin/python
from math import log
def isPrime(n):
for i in range(2, int(n**0.5) + 1):
if n % i == 0:
return False
return True
if __name__ == '__main__':
Euler = 0.57721566490153286
m = 0
for x in range(2, 10_000_000):
if isPrime(x):
m += log(1-(1/x)) + (1/x)
print("MM =", Euler + m)