モンテカルロ法により円周率を求める

Qiita では，これも定期的に出てくるやつ。

長々しいプログラムが散見されるが，以下のように簡潔に書く。

import numpy as np
def sim(n):
    xy = np.random.random((n, 2))**2
    print(np.mean(xy.sum(axis=1) < 1)*4)

< 1 のところを <= 1 にしないといけないのではないかとか，np.random.random は 0 以上，1 未満 の一様乱数なので，0 以上，1 以下の一様乱数を使わないといけないのではないかとか，おかしなことをいう人もいるが，そのようなおかしな人を相手にする必要はない。