invmaxtr.m

来自「数值分析各种公式的解法」· M 代码 · 共 14 行

function[mx,mt,my]=invmaxtr(a,esp)
n=length(a);
x0=diag(ones(n));
x=inv(a)*x0;
k=0;
while norm(x-x0)>esp
    k=k+1
    q=x;
    y=x/maxnorm(x)
    x=inv(a)*y;
    x0=q;
end
mt=1/maxnorm(x)
my=y