R 您可以使用体积库的可变集成限制吗
我正在尝试使用一个可变的积分极限来使用容积库 例如:R 您可以使用体积库的可变集成限制吗,r,numerical-integration,R,Numerical Integration,我正在尝试使用一个可变的积分极限来使用容积库 例如: adaptIntegrate(function(x) {x[1]*x[2]}, lowerLimit = c(0,0), upperLimit = c(x[2],1))$integral 这不管用 感谢Cubature是一个用于超立方体上自适应多维积分的软件包。我也在寻找在任何特定区域进行多维集成的包(即使用变量作为限制),有人知道吗 此外,函数的定义不合适。我猜你想要的是 function(x1,x2) x1*x2 积分
adaptIntegrate(function(x) {x[1]*x[2]},
lowerLimit = c(0,0),
upperLimit = c(x[2],1))$integral
这不管用
感谢Cubature是一个用于超立方体上自适应多维积分的软件包。我也在寻找在任何特定区域进行多维集成的包(即使用变量作为限制),有人知道吗 此外,函数的定义不合适。我猜你想要的是
function(x1,x2) x1*x2
积分的形式为dx[1]dx[2] 积分域是一个三角形(画一张图!),有顶点(0,0)、(0,1)、(1,1)。对于这样的域,请使用
simplicalcubature
包
> f <- function(x) x[1]*x[2]
> S <- cbind(c(0,0),c(0,1),c(1,1)) # the triangle (simplex)
> library(SimplicialCubature)
> adaptIntegrateSimplex(f, S)
$integral
[1] 0.125
$estAbsError
[1] 1.25e-13
$functionEvaluations
[1] 32
>f S库(SimplicalTubature)
>自适应集成单纯形(f,S)
$integral
[1] 0.125
$estAbsError
[1] 1.25e-13
$functionEvaluations
[1] 32
谢谢,现在我知道了为什么变量限制不适用于容积。