mci.r

mci.sapply <- function(x, u) {
  unlist(sapply(x, 
    function(x) {
      g <- x * exp( -(u * x)^2 / 2)
      mean(g) / sqrt(2 * pi) + 0.5
    }
  ))
}