# file "RMCint2" g = function(x,y) { cos(sqrt(x*y)) } M = 1000 xlist = runif(M,0,5) ylist = runif(M,0,4) funclist = 5 * 4 * g(xlist,ylist) print( mean(funclist) ) print( sd(funclist) / sqrt(M) )