# R code to integrate h() using classical Monte Carlo method
# Example 1

ex1.fcn <- function(x){
h <- 4/(1 + x^2)
return(h)
}

my.n <- 10000

my.a <- 0; my.b <- 1

my.rand.x <- runif(my.n, min=my.a, max=my.b)

I.approx <- ((my.b - my.a)/my.n)*sum(ex1.fcn(my.rand.x))

print(I.approx)