import numpy as np
import matplotlib.pyplot as plt
n=100000
x=np.random.rand(n)
y=np.random.rand(n)
c=x[x**2+y**2<=1]
print(4*c.size/n)
plt.scatter(x[x**2+y**2<=1],y[x**2+y**2<=1])
plt.scatter(x[x**2+y**2>1],y[x**2+y**2>1],c='r')
plt.axis('equal')
python蒙特卡洛法计算π值
最新推荐文章于 2023-12-08 17:25:25 发布