function(x, alfa)
{
	n <- length(x)
	cf1 <- qchisq(alfa/2, n - 1)
	cf2 <- qchisq(1 - alfa/2, n - 1)
	lower <- ((n - 1) * var(x))/cf2
	upper <- ((n - 1) * var(x))/cf1
	conint.sigmsq <- c(lower, upper)
	conint.sigmsq
}
