#source("c:/Work07/ShortCourse07/Programs/Mci.txt",print.eval=TRUE)#
#This gives a simple monte carlo integration#
nsim<-100000;u<-runif(nsim);
#The function to be integrated
mci.ex <- function(x){(cos(50*x)+sin(20*x))^2}
plot(function(x)mci.ex(x), xlim=c(0,1),ylim=c(0,4))
#The monte carlo sum
sum(mci.ex(u))/nsim