正規分布の近似(8)

正規分布に従う乱数としては、
Box-Muller transformというものが知られており、
これは、2つの0から1の一様乱数をr1, r2とすると、

 \sqrt{-2log(r_1)}cos(2\pi{r_2})
 \sqrt{-2log(r_1)}sin(2\pi{r_2})

で2つ与えられる。
これをいつものように1億回繰り返したのと、
前回のとかかった時間を比べると、

前回:20.8s
今回:22.1s

というわけで、前回のほうがちょっとだけ速い。


ほかにいい手がないかはおいといて、
次からはここから派生する問題について考える。