Taglibro de H
モンテカルロ法
http://ito-hi.blog.so-net.ne.jp/2006-09-05
短くすればよいと言うものではないが,二次元配列を利用すると良い
配列の二乗,列和に colSums,列和が 1 以下になる平均値(!)を求め(sum(p)/n って,mean(p) のこと)4倍する
set.seed(124)
n <- 100000
mean(colSums(matrix(runif(2*n), 2)^2) < 1)*4
答え 3.14484
モンテカルロ法
http://ito-hi.blog.so-net.ne.jp/2006-09-05
短くすればよいと言うものではないが,二次元配列を利用すると良い
配列の二乗,列和に colSums,列和が 1 以下になる平均値(!)を求め(sum(p)/n って,mean(p) のこと)4倍する
set.seed(124)
n <- 100000
mean(colSums(matrix(runif(2*n), 2)^2) < 1)*4
答え 3.14484
※コメント投稿者のブログIDはブログ作成者のみに通知されます