cuhre函数如何在R中起作用？

Question

我正在阅读我的多变量类注释，它使用来自R2Cuba包的cuhre函数来根据变量X＆amp;的规则来评估概率。是的。这里是代码的完整部分：

integrand <- function(z){
    x <- z[1]
    y <- z[2]
    if (0<x & x<1 & 0<y & y<1 & x+y>1)
        f = 6*x*y^2
    else
        f = 0
    return(f)
}

NDIM = 2
NCOMP = 1
int <- cuhre(NDIM, NCOMP, integrand, rel.tol = 1e-3, abs.tol = 1e-12, 
             flags = list(verbose = 2, final = 0))$value

int

结果是：

> int
[1] 0.8998863

据我所知，这是基于规则的概率：

0<x<1 & 0<y<1 & x+y>1

我无法理解的是，被积函数被定义为以z为参数的函数，因此只传递1个参数。当它在下面被调用时，它没有任何参数，并且集成发生两次并且概率被保存在int中。我有50％的人，但不是很清楚它是如何工作的。在某个地方，我认为我们可以使用Cuhre来计算边际概率，我们可以吗？