Simple, but not too simple estimation of pi

piday.R

pi.ratio<- function(n) {
	x<- runif(n, min=-1,max=1)
	y<- runif(n, min=-1,max=1)
	return(4*sum(x^2 + y^2 <= 1.0)/n)
	}

pi.est<- function(iter) {
	sum.pi<-stor.rat<-numeric(0)
	for (i in 1:iter) {
		stor.rat[i]<- pi.ratio(1000)
		sum.pi[i]<- sum(stor.rat[1:i])/i
		}
	return(sum.pi)
	}