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（结果具有随机性 不一定完全一样）