function(x, alfa, beta, m)
{
	n <- length(x)
	pred <- jay(m, 1)
	for(i in 1:m) {
		xs <- sample(x, n, replace
			 = T)
		pred[i] <- mean(xs)
	}
	pred1 <- quantile(pred, alfa/2)
	pred2 <- quantile(pred, 1 - alfa/2)
	toler1 <- qbinom(beta/2, n, pred1)
	toler2 <- qbinom(1 - beta/2, n, 
		pred2)
	res <- list(toler1, toler2)
	names(res) <- c("Lower", "Upper")
	res
}
