# function for problem 4 wlln <- function(n, eps){ 1 - (pgamma(n/2 + n*eps, n, 2) - pgamma(n/2 - n*eps, n, 2)) } # try it out print(wlln(100, 0.01)) # make a plot plot(wlln(1:10000, 0.01)) # find the answer to part (d) for(n in 10000:20000){ val <- wlln(n, 0.01) if(val < 0.01){ print(n) print(val) break } }