i<-0
t<-0
n<-10000000
while (i<=10000000){
x<-runif(1,-1,1)
y<-runif(1,-1,1)
if(sqrt((x^2+y^2))<=1){t<-t+1}
i<-i+1
}
t
options(digits = 10)
(t/n)*4
运行如下:
i<-0
t<-0
n<-10000000
while (i<=10000000){
x<-runif(1,-1,1)
y<-runif(1,-1,1)
if(sqrt((x^2+y^2))<=1){t<-t+1}
i<-i+1
}
t
options(digits = 10)
(t/n)*4
运行如下: