#
# Crude Monte Carlo - Chapter 9, Exercise 1
#
niter <- 10
n <- 10
integral <- rep(0,niter)
cat(" Values computed: ")
for (i in 1:niter){
  integral[i] <- sum(exp(runif(n)))/n
  cat(integral[i])
  if (i==n)
    cat(".")
  else
    cat(", ")
  if (i==n/2)
  cat("\n                  ")
}
cat("\n")
cat(" Mean is",mean(integral),"and standard deviation is",sd(integral),".")
cat("\n")