蒙特卡罗方法计算圆周率

import numpy as np
r=1#定义内接圆半径
#随机数生成个数
rand_num=[100,1000,10000,100000,1000000,10000000]
#根据生成随机数个数的不同计算的圆周率
for N in rand_num:
    #在边长为2的正方形区域内生成随机点坐标(x,y)
    x=2*np.random.random_sample(N)-1
    y=2*np.random.random_sample(N)-1
    in_circle_point_num=0
    #计算落在内接圆区域内的随机点数
    for point_count in range(len(x)):
        #判断随机点是否落在内接圆区域之内
        if(x[point_count]*x[point_count]+
           y[point_count]*y[point_count]<
           r*r):
            in_circle_point_num+=1
    print('N=',str(N),'pi=',str(4.0*in_circle_point_num/N))