どう書く?org : π

original page : http://ja.doukaku.org/210/

円周率を計算してください。

積分を計算するも、素朴な方法も、速さを目指すも、LLで計算する意味を問うもあるでしょう。

http://ja.wikipedia.org/wiki/%E5%86%86%E5%91%A8%E7%8E%87

#7771

original page : http://ja.doukaku.org/comment/7771/

Rにはpiという組み込み変数がありますが、計算であればRらしくモンテカルロで。

pi.montecarlo <- function(n){
  mean(replicate(n, sum(runif(2)^2) < 1)) * 4
}