用蒙特卡洛法找圆周率

import random
#用蒙特卡洛法找圆周率,
m=pow(int(input()),2)#输入点点数量
n=0#统计落入圆形区域的数量
for i in range(m+1):
    x=random.random()#横坐标，范围0~1
    y=random.random()#纵坐标，范围0~1
    r=pow((x*x+y*y),0.5)#（x，y）在圆内
    if r<1:
        n=n+1
    load=(i/m)*100
    print("\r完成{0:.10f}%计算".format(load),end="")#进度显示
p=n/m#落在圆内的概率
pi=4*p
print("\n通过蒙特卡洛算法计算出圆周率是：{}".format(pi))