maxtr.m

数值分析各种公式的解法

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