R:モンテカルロ法で求める円周率のコードを解く。~各務ゼミ生に送る~
おはこんばんにちわ。チャパティです。
この記事は下の記事の続きになります。
ということで、本日も描いていきましょう。
円周率を求めたコードset.seed(422)iIter<- 10000000z <- 0w <-NULLpdf( file = "円周率.pdf", height = 5, width = 7, family = "Times")for( i in 1:iIter){ x<-runif(1) y<-runif(1) if(x^2+y^2<1