蒙特卡罗方法生成圆周率

Monte Carlo to generate $\pi$

import random


def calpai():
    n = 100000000
    r = 1.0
    a, b = (0.0, 0.0)
    x_neg, x_pos = a - r, a + r
    y_neg, y_pos = b - r, b + r

    count = 0
    for i in range(0, n):
        x = random.uniform(x_neg, x_pos)
        y = random.uniform(y_neg, y_pos)
        if x*x + y*y <= 1.0:
            count += 1
    print((count / n) * 4)
    return


calpai()

这是蒙特卡罗方法生成pi值，希望对大家有帮助，谢谢朋友们！