使用环境为python IDLE3.5.
#利用蒙特卡洛方法计算PI值
from random import random
from math import sqrt
from time import clock
DARTS = 2**26 #这个时候已经比较精确了,耗时也将近1分钟了。修改DARTS可以进一步提高计算精度。
hits=0
clock() #引入time.clock()用于统计耗时
for i in range(1,int(DARTS)):
x,y =random(),random()
dist = sqrt(x**2+y**2)
if dist <=1.0:
hits =hits+1
pi = 4*(hits/DARTS)
print("PI的值是%s"%pi)
print("程序运行的时间是%-5.5ss"% clock())