function(x, alfa, beta)
{
	n <- length(x)
	za <- qnorm(1 - alfa/2, 0, 1)
	zb <- qnorm(1 - beta/2, 0, 1)
	tab <- zb/(1 - za^2/(2 * n))
	tab <- tab + (za * sqrt(1 + zb^2/2 - 
		za^2/(2 * n)))/(sqrt(n) * (1 -
		za^2/(2 * n)))
	lower <- mean(x) - tab * std(x)
	upper <- mean(x) + tab * std(x)
	tolint.norm <- c(lower, upper)
	tolint.norm
}
