使用 Gauss-Legendre 求积的四重积分

Question

我想在 Matlab 中使用 Gauss-Legendre 求积来计算四重积分。我发现以下代码适用于双积分，但在切换到 4 个变量时，我不知道如何操作。

a = 2; b=2;
[x,w] = lgwt(50);

y=(b*(1+x))/2;
x=(a*(1+x))/2;
[X, Y] = meshgrid(x,y);
fun = @(x,y) exp(1i*x)./(sqrt(x.^2+y.^2));
t = w * fun(X,Y) * w'

我想集成的功能：

fun = @(x,y,z,t) 1./sqrt(x.^2+y.^2+z.^2+t.^2)

有人有想法吗？

Answer 1

This solution 似乎是您的最佳选择。

整合fun = @(x,y,z,t)的最佳方式来自 a <= x <= b、c <= y <= d、e <= z <= g 和 h <= t <= i 是嵌套 integral3() 和 integral()。它们都是 Matlab 的内置函数。

Q = integral(@(x)integral3(@(y,z,t)f(x,y,z,t),c,d,e,g,h,i),a,b,'ArrayValued',true);