蒙特卡洛问题求PI值的原理
正方形的面积为4r^2
正方形内四分之一圆的面积为:Πr^2
正方形内四分之一圆的面积与正方形的面积的比值为Π/4
所以Π的值为4 * (四分之一圆的面积/正方形的面积)
import random
N = int(input("请输入要随机产生的点数:"))
X=0
pi=0
for i in range(N):
a = random.uniform(0,1) #随机产生0-1内的数
b = random.uniform(0,1)
if a*a + b*b <= 1: #随机产生的坐标点坐落在以原点为中心半径为1的四分之一圆内的个数(X^2+Y^2=1 表示以原点为中心半径为1的圆)
X=X+1
pi=(X/N)*4
print('pi近似等于:{0}'.format(pi))