用一把小米来求一个圆的面积。首先求一个已知的矩形S1的面积,此时使用N1数量的小米。然后将小米平坦在S2的圆上,使用了N2数量的小米。因此矩形和圆的关系可以表示成N1/S1=N2/S2。这便是小米贡献。又称蒙特·卡罗方法,也称统计模拟方法,是二十世纪四十年代中期由于科学技术的发展和电子计算机的发明,而被提出的一种以概率统计理论为指导的一类非常重要的数值计算方法。是指使用随机数来解决很多计算问题的方法。使用random函数随机生成单位正方形之间的点,记录其分布情况,再求出1/4圆内的点与单位正方形内点的比值。
import random
import numpy as np
date=[]
for i in range (10000):
x=random.uniform(0,2)
y=random.uniform(0,2)
date.append([x,y])
date=np.array(date)
size=0
x1=date[0,0]
y1=date[0,1]
for j in range(0,len(date)):
x1=date[j,0]
y1=date[j,1]
if ((x1-1)**2+(y1-1)**2)**0.5<=1 :
size=size+1
ans=(size/10000)*4
print(ans)