蒙特卡罗模拟双积分

Question

  

可能重复：
  calculating double integrals in R quickly

Braun和Murdoch提出的问题， Monte Carlo Simulation ：

  

doubleintegral：整数（从0到1）积分（从0到1）fct'n：   cos（x-y）I常数：dxdy

对于只有“dx”的积分，我会使用类似的东西：

integrand <- function(x)
 {
    exp(x)
 }
integrate(integrand, lower = 1, upper = pi) $value

有人可以给我一个提示，如何在R-Code中写一个双积分？