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