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)
	}
}

x <- rpois(n, lambda)
mx <- mean(x)
T <- Test(mx, n, lambda0)

if (abs(T) > q) {
	cat('H0 wird abgelehnt.\n')
}
if (abs(T) <= q) {
	cat('H0 wird nicht abgelehnt.\n')
}