数据分析之蒙特卡洛模拟

from __future__ import division
import random
import time

num=1
for j in range(1,10):
    startT=time.clock()
    #落入圆内计数
    counter=0
    #往正方形内扔了10的j次方
    for i in range(10**j):
        x=random.uniform(-1,1)
        y=random.uniform(-1,1)
        #落入园内计数
        if x**2+y**2<1:
            counter+=1
    endT=time.clock()
    print(num)
    print(‘pi:{0}‘.format(4*(counter/10**j)))
    print(‘时间:{0}‘.format(startT-endT))
    print()
    num+=1