R语言while循环计算圆周率

原理说明

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

代码展示

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

计算结果

> setpi(100000)
[1] 3.14328

R语言while循环计算圆周率相关推荐

r语言for循环的c(),R语言简单for循环（二）
创建一个简单数据框 Sys.setlocale('LC_ALL','C') ## [1] "C" library(tidyverse) ## -- Attaching packag ...
R语言sd函数计算数值标准差实战（Standard Deviation）
R语言sd函数计算数值标准差实战(Standard Deviation) 目录 R语言sd函数计算数值标准差实战(Standard Deviation) #基本语法 #sd
R语言自定义函数计算dataframe每列中的缺失值NA的个数、缺失值问题及其填充示例
R语言自定义函数计算dataframe每列中的缺失值NA的个数.缺失值问题及其填充示例目录
R语言Eta squared计算实战：Eta squared表示可以用模型中给定的变量解释的方差的比例、拟合方差分析模型(two-way ANOVA)、计算Eta Squared
R语言Eta squared计算实战:Eta squared表示可以用模型中给定的变量解释的方差的比例.拟合方差分析模型(two-way ANOVA).计算Eta Squared 目录
R语言IQR函数计算四分位数范围IQR（Interquartile Range）实战
R语言IQR函数计算四分位数范围IQR(Interquartile Range)实战目录 R语言IQR函数计算四分位数范围IQR(Interquartile Range)实战 #基本语法
R语言length函数计算向量、列表、字符串长度实战
R语言length函数计算向量.列表.字符串长度实战目录 R语言length函数计算向量.列表.字符串长度实战 #基本语法
R语言mode函数计算众数实战
R语言mode函数计算众数实战目录 R语言mode函数计算众数实战 #手动编写众数函数 #存在多个众数的情况
R语言distVincentyEllipsoid函数计算大圆距离实战（Great Circle Distance）
R语言distVincentyEllipsoid函数计算大圆距离实战(Great Circle Distance) 目录 R语言distVincentyEllipsoid函数计算大圆距离实战(Grea ...
R语言distVincentySphere函数计算大圆距离实战（Great Circle Distance）
R语言distVincentySphere函数计算大圆距离实战(Great Circle Distance) 目录 R语言distVincentySphere函数计算大圆距离实战(Great Circ ...

R语言while循环计算圆周率

R语言while循环计算圆周率

原理说明

代码展示

计算结果

R语言while循环计算圆周率相关推荐

最新文章

热门文章