function(x, alfa)
{
	n <- length(x)
	ft <- 2 * sum(x) + 2
	upper <- qchisq(1 - alfa/2, ft)/(2 * 
		n)
	lower <- qchisq(alfa/2, ft)/(2 * n)
	conint.pois <- c(lower, upper)
	conint.pois
}
