Python蒙特卡罗算法简介

import math
import random
m = 10000
n = 0
for i in range(m):
# x、y为0-1之间的随机数
    x = random.random()
    y = random.random()
    # 若点(x,y) 属于图中1/4圆内 则有效个数+1
    if math.sqrt(x**2 + y**2) < 1:
        n += 1
# 计算pi
pi = 4 * n / m
print("pi = {}".format(pi))
 
# pi = 3.1508（结果具有随机性 不一定完全一样）