目前，我使用了ode45和一个输出3945x9矩阵(C(的函数。这些应该是3x3矩阵，所以我做了C=整形(C.'，3,3，[](。现在，我有一个3x3x3945矩阵。我想做的是找出每个3x3矩阵的误差。这是用C*C完成的眼睛(3(。然而，我现在不知道如何使用我的所有数据来做到这一点。它将适用于类似C(：，：，1(*C(：、：、1(的东西眼睛(3(。但是，不是C(：，：，：(*C(：：，：：(-眼睛(3(。
使用nD矩阵乘法例程pagemtimes。例如，
pagemtimes(C,'none',C,'transpose') - eye(3)