πを求めるシミュレーション

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