有人知道adaptIntegrate是否接受矢量化被积函数吗？
例如。
alpha<-c(1,2)
f <- function(z){
(z[1]+z[2])*alpha
}
adaptIntegrate(f,lower=c(1, 3), upper=c(2, 4),tol=0.01) ##does not work
adaptIntegrate是我想要的吗？
答案 0 :(得分：1)
自己解决 - 在adaptIntegrate的参数中添加fDim=2。更一般地说，
fDim=length(alpha)