R中的双/三积分

Question

我试图在R中求解下面提到的（图像）方程，用于区域上的双积分。我通过以下帖子的帮助使用以下代码： Triple integral in R (how to specifying the domain)

不确定，如果我遗漏任何东西，好像我改变了A的值，它仍然给我相同的结果。以下代码是否正确描述了我想要解决的等式？

library(cubature)

PW <-  0.0144
Cw <- 400/(0.91*70)
A <- 1960                
w <- as.numeric(PW, Cw, A)

lower <- rep(0,3)
upper <- rep(1,3)

fxyz <- function(w) {
    x <- w[1]
    y <- w[2]
    z <- w[3]
    as.numeric(x <= y)*as.numeric(y <= z)
}

Int <- adaptIntegrate(f=fxyz,lowerLimit=lower,upperLimit=upper,doChecking=TRUE,
       maxEval=2000000,absError=10e-5,tol=1e-5)

感谢， Krina

Answer 1

width=w, height=h 定义的本地w会影响您在上面定义的全局fxyz。更改w会对全局A产生影响，但这与w的本地w无关。更改fxyz对积分值没有影响的原因是积分和积分限制都不依赖于A。

回答你的问题：＆＃34;以下代码是否正确描述了我试图解决的等式？＆＃34;：积分本身的代码似乎是正确的，但代码为{{1显然不是。您似乎打算在该函数与参数A之间存在某种依赖关系，但您还没有解释该依赖关系应该是什么。也许你可以在这里或在另一个问题上澄清。