模拟蒙特卡罗方法计算圆周率

使用蒙特卡罗方法计算出圆周率,并输出程序运行的时间，以及画出散点图。

import random as r#调用random库
import time as t#调用time库
import math#调用math库
import matplotlib.pyplot as plt#调用matplotlib库

num =0#记录在单位圆内的点数
number=100000#表示总共模拟的点数
time_1=t.time()#记录开始时间
x_list = []#产生空的列表用于后面记录x
y_list = []#产生空的列表用于后面记录y

for i in range(number):#使用for循环模拟number这个点
    x=r.random()#产生0到1之间的随机数并赋值给x
    y=r.random()#产生0到1之间的随机数并赋值给y
    z=math.sqrt(x**2+y**2)#计算产生的点到原点的距离
    if z<=1:#判断z，如果z小于1，则该点在圆内
        num=num+1#只要点在圆内则num加一
        x_list.append(x)#x坐标记录到x_list中
        y_list.append(y)#y坐标记录到x_list中
print("Pi={}".format(num/number*4))#输出得到的圆周率
plt.scatter(x_list, y_list, s=1)#使用scatter函数绘制散点图
plt.show()#显示绘制的图形
time_2=t.time()#产生结束时间
print("Time={:.6f}s".format(time_2-time_1))#输出程序运行的时间

产生如下的结果：

Pi=3.13492
Time=1.017492s