n <- 10
lambda <- 2
lambda0 <- 1
alpha <- 0.05
q <- qnorm(1 - alpha / 2.0)

Test <- function(mx, n, lambda0) {
	if (mx > 0) {
		return(sqrt(n) * (mx - lambda0) / sqrt(mx))
	}
	else {
		return(0)
	}
}

cnt = 0
max = 10000

for (i in 1:max) {
	x <- rpois(n, lambda)
	mx <- mean(x)
	T <- Test(mx, n, lambda0)
	if (abs(T) > q) {
		cnt = cnt + 1
	}
}

cat('Der Test wurde in', (1-cnt/max)*100, 'Prozent der Fälle NICHT abgelehnt.\n')