Monte Carlo Pi

pi.py

import random
# number of all points in 1/4 circle / number of all points = pi/4
def pi(n): return sum( [random.random()**2 + random.random()**2 <= 1 for _ in range(n)] ) / n * 4
# % python3 -i pi.py
print(pi(100))
# 3.52
print(pi(1000))
# 3.104
print(pi(10000))
# 3.1448
print(pi(100000))
# 3.13572
# print(pi(1000000))
# 3.14102
# print(pi(10000000))
# 3.1422764
# print(pi(100000000))
# 3.14178852