R语言while循环计算圆周率

最新推荐文章于 2023-04-08 15:29:06 发布

莫言于方

最新推荐文章于 2023-04-08 15:29:06 发布

阅读量7.2k

点赞数 4

分类专栏： R语言文章标签： R 循环概率

本文链接：https://blog.csdn.net/qq_31468321/article/details/83479795

版权

R语言专栏收录该内容

4 篇文章 0 订阅

订阅专栏

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)