蒙特卡罗方法 python

蒙特卡罗法也称统计模拟法、统计试验法。是把概率现象作为研究对象的数值模拟方法。是按抽样调查法求取统计值来推定未知特性量的计算方法。

例    有一个半径为r=1的圆和边长为1的正方形，圆的面积为Πr**2=Π  ，则正方形内部的相切圆的面积为整个圆的1/4，也就是Π/4，正方形的面积为1。然后我们向正方形中随机打点，就会有一定的概率落在圆中。这样我们就可以得到落在圆中的概率就是=圆的面积/正方形面积。

import random
t=0
for i in range(10000):
    x=random.random()
    y=random.random()
    dis=x**2+y**2
    if dis<=1:
        t+=1
print(4*t/10000)