# 填充圆方法一
n=5000
s1=0
s=0
for i in range(0,n):
a=np.random.uniform(-2,4,1)
b=np.random.uniform(-2,4,1)
if(((a-1)**2+b**2)>=1):#((a+1)**2+b**2)<=1 a<=0,,((a-2)**2+b**2)>=4
plt.scatter(a,b,c = 'r',marker = '+')
else:
plt.scatter(a,b,c = 'b',marker = '+')
plt.xlim(-2,4)
plt.ylim(-2,4)
plt.gca().set_aspect('equal', adjustable='box')#
蒙特卡洛方法填充一个圆
于 2023-03-05 09:49:10 首次发布