你可以使用与Cubature库集成的可变限制

Question

我正在尝试使用Cubature库使用可变的集成限制。

例如：

fnc(value) { this.x = value; }

这不起作用

谢谢

Answer 1

Cubature是一个用于在超立方体上进行自适应多维集成的软件包。我也在寻找在任何特定区域进行多维集成的软件包（即使用变量作为限制），谁都知道呢？

此外，功能的定义也不合适。我想你想要的是

function(x1,x2) x1*x2

Answer 2

积分的形式为\int_O^1 \int_0^{x[2]} f dx[1] dx[2]

积分的域是一个三角形（创建图片！），其顶点为（0,0），（0,1），（1,1）。对于这样的域，请使用SimplicialCubature包。

> 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