蒙特卡罗方法python简单实现

Posted 2020-08-22

蒙特卡罗方法是一种计算方法。原理是通过大量随机样本，去了解一个系统，进而得到所要计算的值。

1、估算pi值

import random
n = 10000
k = 0
for i in range(n):
    x = random.uniform(-1,1)
    y = random.uniform(-1,1)
    if x**2 + y**2 <1:
        k +=1
print(4*float(k)/float(n))

2、计算y=x^2在[0,1]的积分

import random
n = 1000
k = 0
for i in range(n):
    x = random.uniform(0,1)
    y = random.uniform(0,1)
    if x**2 > y:
        k+=1
print(float(k)/float(n))