# file "RMCint4"

h = function(x,y) { exp(-y^2) * cos(sqrt(x*y)) }

M = 1000

lambda = 5

xlist = runif(M)

ylist = rexp(M, lambda)

funclist = (1/lambda) * exp(lambda*ylist) * h(xlist,ylist)

print( mean(funclist) )

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