pracma软件包中integral3的R错误

我有一个简单的三维积分，我想用

pracma

软件包中的

integral3

求解：

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)

库（pracma）
W虽然您的函数
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

W