蒙特卡洛思想求矩形内切圆的面积

#边长为4的正方形矩形
l=4
#半径为2的内切圆，圆心在（0，0）处
r=2
#用来记录落在圆内的点的个数
num=0
#取100000次随即点数
for i in range(100000):
    x=random.uniform(-2,2)#点的x轴坐标
    y=random.uniform(-2,2)#点的y轴坐标
    t=math.sqrt(x*x+y*y)
    if t<=2:
        num+=1
#圆的面积＝园内点的个数/点总个数*矩形面积
s=num/100000*16
print(s)

 得到的面积是12.55104，与直接计算的结果12.56非常接近。