rng(1)
% 生成用来做n-mode积的矩阵
A1 = randi([1,2],[2,3]);
A2 = randi([1,2],[4,3]);
A3 = randi([1,2],[3,3]);
% 生成用来做外积的向量
a1 = randi([1,2],[3,1]);
a2 = randi([1,2],[3,1]);
a3 = randi([1,2],[3,1]);

CP = outer_prod(outer_prod(A1*a1, A2*a2), A3*a3)
TU = modeN_prod(modeN_prod(modeN_prod(outer_prod(outer_prod(a1, a2), a3), A1, 1), A2, 2), A3, 3)

% 结果完全相同：
% CP(:,:,1) =
%    128   224   128   256
%    160   280   160   320
% CP(:,:,2) =
%     96   168    96   192
%    120   210   120   240
% CP(:,:,3) =
%    144   252   144   288
%    180   315   180   360
% 
% TU(:,:,1) =
%    128   224   128   256
%    160   280   160   320
% TU(:,:,2) =
%     96   168    96   192
%    120   210   120   240
% TU(:,:,3) =
%    144   252   144   288
%    180   315   180   360