juzhendezhishu.m

计算魔方矩阵的演示程序源代码

%example3_3.m

A=rand(4,4);E=0*A;
F=eye(size(E));N=1;
while norm(E+F-E,1)>0,
    E=E+F;
    F=A*F/N;
    N=N+1;
end
E
E2=expm(A)