r - モンテカルロ シミュレーション二重積分

Question

重複の可能性:
R の二重積分をすばやく計算する

モンテカルロ シミュレーションの Braun と Murdoch からの質問:

doubleintegral: 積分 (0 から 1) 積分 (0 から 1) fct'n: cos(xy) I-定数: dxdy

「dx」だけの積分には、次のようなものを使用します。

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

Rコードに二重積分を書く方法を教えてください。