蒙特卡洛方法填充一个圆

自我o救赎

已于 2023-03-09 16:00:13 修改

阅读量118

点赞数

分类专栏： Python语言专栏文章标签： python matplotlib numpy

于 2023-03-05 09:49:10 首次发布

本文链接：https://blog.csdn.net/y15520833229/article/details/129342751

版权

Python语言专栏专栏收录该内容

69 篇文章 2 订阅

订阅专栏

# 填充圆方法一
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')#