|
发表于 2017-10-9 08:59:27
|
显示全部楼层
data<-data.frame(n=n,x=x,y=y,pt=pt)
这句话不对。
代码可能这么写比较合适。
set.seed(100)
pt<-0
data<-data.frame(N=integer(0),X=integer(0),Y=integer(0),PT=integer(0))
for(n in 1:100){
x<-runif(n,-1,1)
y<-runif(n,-1,1)
pt<-sum(x^2+y^2<=1)/n*4
data1<-data.frame(N=n,X=x,Y=y,PT=pt)
rbind(data,data1)
}
|
|