圆周率蒙特卡洛法

[python]代码库

from random import random
from time import perf_counter

DARTS = 100000000
hits = 0
start = perf_counter()
for i in range(DARTS + 1):
    x, y = random(), random()
    dist = pow(x ** 2 + y ** 2, 0.5)
    if dist <= 1.0:
        hits += 1;
pi = 4 * (hits / DARTS)
print('圆周率是 {}'.format(pi))
print('运行时间为：{:.5f}s'.format(perf_counter() - start))

[代码运行效果截图]