版权声明：转载请声明出处，谢谢！ https://blog.csdn.net/qq_31468321/article/details/83479795
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