N <- N ; maxbound = 4 one_side_normal <- function(N) { f <- numeric(N); for (i in 1:N) { X = rexp(1,rate=1); uniform = runif(1); while (uni