# file "RMCint2"

g = function(x,y) { cos(sqrt(x*y)) }

M = 100000

xlist = runif(M,0,5)

ylist = runif(M,0,4)

funclist = 5 * 4 * g(xlist,ylist)

print( mean(funclist) )

print( sd(funclist) / sqrt(M) )