。
import random
def Monte_Carlo_pi(N):
#计数器
count = 0;
#循环次数
for i in range(N):
#生成0-1之间的随机数
x, y = random.random(), random.random()
#计算随机生成的点到原点的距离
d = pow(x ** 2 + y ** 2, 0.5)
#若其距离小于等于1,则随机点在圆内
if d <= 1:
count += 1
#根据蒙特卡罗法的原理,得出π的具体数值
pi = count / N * 4
return pi
if name == 'main':
N_def = 1000000
pi = Monte_Carlo_pi(N_def)
print("通过 %d 次随机投点的模拟,计算的圆周率π的大致值为:%s" % (N_def, pi))