matlab - MATLAB 三重积分的问题

Question

我试图评估的方程式：https ://i.stack.imgur.com/W0HxL.jpg

我正在尝试评估这个积分（和方程），其中“a”和“x”只是输入，而 x 是 (x,y,z)。另外：x 将具有 [a,b,c] 的形式。在下面的“int_func”中，我采用了转置和原始的点积，以产生下面可观察到的 ^2 幂，例如：(X(:,1)-x).^2

我到目前为止的代码是：

int_func = @(x,y,z) func([x,y,z]).*exp((-3/(2*a^2)).*((X(:,1)-x).^2 + (X(:,2)-y).^2 + (X(:,3)-z).^2));

val = integral3(int_func,-inf,inf,-inf,inf,-inf,inf);
disp(val * (3/(2*pi*a^2))^(3/2));

在哪里：

X = [0, 0, 0];

a = 1;

使用这段代码，我应该得到0.0370 的答案，但我得到的却是 0.5556。

有人可以建议另一种形式的评估这个积分或发现我列出的内容有什么问题吗？

标签： matlabconvolutiongaussianintegral