import random
import time
num=1
for i in range(3,6):
print('第{0}次试验'.format(num))
start = time.clock()
count = 0
for j in range(10**i):
x = random.uniform(-1,1)
y = random.uniform(-1,1)
if(x**2+y**2<=1):
count=count+1
end = time.clock()
print('pi:{0}'.format(4*count/(10**i)))
print('时间:{0}'.format(end-start))
num=num+1
蒙特卡洛模拟计算PI
最新推荐文章于 2024-06-14 10:00:00 发布