l5_1.m

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

function l5_1(A,Y0,e,N)
k=1;u=1;
[L,U,P]=lu(A);
Y0=P*Y0;
Y(:,1)=Y0;
while k<N
    Z=L\Y(:,k);
    X(:,k+1)=U\Z;
    t=max(abs(X(:,k+1)));
    Y(:,k+1)=X(:,k+1)/t;
    if abs(1/t-1/u)<e
        1/t,Y,k
        return
    end
    k=k+1;
    u=t;
end
if k>=N
    'error,超过迭代次数'
end