maxtr.m

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

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