3 π的计算

π的计算

计算思想：采用蒙特卡洛方法，利用随机实验求解问题。构建一个单位正方形和1/4圆，向正方形内随机抛撒点，计算抛撒的点落在1/4圆内的概率。随着抛撒点的增多，得到的值越精确。

输入：抛撒点的数量
处理：对于抛撒的每个点，计算点到圆心的距离，通过距离判断点在圆内还是圆外。统计落在圆内的点的数量。
输出：π值。

需要用到的库是math库和random库。同时引入计时程序衡量运行时间。

math库常用的函数如下图所示：

random库常用的函数如下图所示：

π的计算程序如下：

# pi.py
from random import random
from math import sqrt
from time import clock
DARTS = 1200#抛撒点的数量
hits = 0
clock()#计时
for i in range(1,DARTS):
    x, y = random(), random()#通过random函数产生随机坐标
    dist = sqrt(x**2 + y**2)#计算坐标到原点的距离
    if dist <= 1.0:#如果距离小于或等于1
        hits = hits + 1#落入1/4圆的数量就加1
pi = 4 * (hits/DARTS)#计算抛撒点落入1/4圆的概率
print("Pi的值是 %s" % pi)
print("程序运行时间是 %-5.5ss" % clock())