蒙特卡罗方法是一种计算方法。原理是通过大量随机样本,去了解一个系统,进而得到所要计算的值。
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))
