使用Gauss-Legendre正交的四重积分

问题描述

我想使用Gauss-Legendre正交来评估Matlab中的四重积分。我发现以下代码可用于双整数，但是在切换到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)

有人有主意吗？

解决方法

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

集成cv2.putText(img,str(label),(x,y + 30),font,1,(0,0),3,cv2.LINE_AA,True) 的最佳方法 fun = @(x,y,z,t)，a <= x <= b，c <= y <= d和e <= z <= g将h <= t <= i与integral3()嵌套。它们都是Matlab的内置功能。

integral()