l2.m

来自「利用幂法反幂法等进行程序运算,可以获得各种适用结果.」· M 代码 · 共 20 行

function l2(A,e,N)
n=size(A);
Y=ones(n(1),1);
k=1;u=0;
while k<=N
    X=A*Y;
    t=max(X);
    Y=X/t;
    if abs(t-u)<e
        t,Y,k
        break
    end
    k=k+1;
    u=t;
end
if k>N
    'error,超出迭代次数'
else
end