如何计算由多维数组定义的函数的Hessian

Question

我在MATLAB中定义了一个尺寸为5×5×5×5的数组，即它是一个4维数组。这个数组的数量定义了4个变量的函数取代某个域的值。

我需要计算这个函数的Hessian。我需要在数组中的每个点计算它（例如，忽略角点情况）。为此，我在某一点获取函数值，计算不同维度的第二个差异和交叉部分。然后我创建一个行列式并检查它的符号。

如何自动完成此过程？也就是说，访问每个点，并使用MATLAB中的diff函数计算行列式？是否有明确的程序？