使用蒙特卡罗方法计算出圆周率,并输出程序运行的时间,以及画出散点图。
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