R语言while循环计算圆周率

R语言while循环计算圆周率

原理说明

计算圆周率的方法很多，在这里我们使用概率的方法来计算圆周率（蒙特卡洛法）。
假设有一个正方形，以其一个顶点为圆心，正方形的边长为半径，画1/4圆。则可以在正方形中随机产生足够多个点，点落在圆内的概率就是圆面积和正方形面积的比值，所以落在圆内的点的个数比上所有的点的个数（在正方形内的点）就等于落在圆内的概率。再按照面积公式。即可计算PI

代码展示

setpi <- function(n){
  cnt<-0
  d<-n
  while(n>=1){
    x<-runif(1,0,1);y<-runif(1,0,1)
    if(x^2+y^2 <= 1)
      cnt<-cnt+1
    n<-n-1
  }
  return(4 * cnt / d)
}
setpi(10000)

计算结果

> setpi(100000)
[1] 3.14328