function(x, n, alfa)
{
	df1 <- 2 * (x + 1)
	df2 <- 2 * (n - x)
	cf1 <- qf(1 - alfa/2, df2 + 2, df1 - 
		2)
	cf2 <- qf(1 - alfa/2, df1, df2)
	lower <- x/(x + cf1 * (n - x + 1))
	upper <- (cf2 * (x + 1))/(n - x + 
		cf2 * (x + 1))
	conint.bino <- c(lower, upper)
	conint.bino
}
