来自实践包的integral3错误

Question

我有一个简单的3d积分，我想用integral3解决pracma包中的问题：

library(pracma)

W <- function(x,y,z){
  x + y + z
}

ymin <- function(x) x
zmin <- function(x,y) x
zmax <- function(x,y) y

integral3(W, xmin = 0, xmax = 1, ymin = ymin, ymax = 1, zmin = zmin, zmax = zmax)

这将生成错误消息：Error in onevec %*% bottom + t %*% dydt : non-conformable arrays。这应该是一个需要集成的有效区域，但我已经确定错误是由我的zmin参数引起的。如果将zmin替换为0，则积分计算正确。基于此，我认为我指定的界限不正确，但我确实相信我正确地遵循了函数文档中的说明。

谢谢

Answer 1

虽然您的函数zmin、zmax看起来很简单，但它们没有向量化。

W <- function(x,y,z) x + y + z

ymin <- function(x) x
zmin <- function(x,y) rep(x, length(y))
zmax <- function(x,y) rep(y, length(x))

integral3(W, xmin = 0, xmax = 1,
             ymin = ymin, ymax = 1, zmin = zmin, zmax = zmax)
## [1] 0.25

我承认，对于几个变量中的函数来说，向量化意味着什么有点让人困惑。