用蒙塔卡罗方法求矩形内切圆面积

 求边长为2的正方形内切圆面积

 

import numpy as np

radius = 1
m = 100000
sqare_x = 2 * np.random.random_sample(m) - 1
sqare_y = 2 * np.random.random_sample(m) - 1
n = 0

for i in range(len(sqare_x)):
    if (sqare_x[i] * sqare_x[i] +
            sqare_y[i] * sqare_y[i] <
            radius * radius): 
        n += 1
print('S: ' + str(4.0*n / m))