モンテカルロシミュレーション

モンテカルロ法による円周率の計算をRで行う

モンテカルロ法を使用して、息抜きにRで円周率の計算を行う。有名問題なので背景はこの辺を参照。 L=1 count = 1000000 x = runif(count,0,L) #0からLの間の乱数を発生させる y = runif(count,0,L) #0からLの間の乱数を発生させる pi = 4*sum(x**2+y**2