モンテカルロ法を使用して、息抜きに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
引用をストックしました
引用するにはまずログインしてください
引用をストックできませんでした。再度お試しください
限定公開記事のため引用できません。