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值,希望对大家有帮助,谢谢朋友们!