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
[main1.m]
a=[3 2;4 5]
maxtr(a,0.0001)
[invmaxtr.m]
function [mx,mt,my]=invmaxtr(a,eps)
内容全部来源于网络收集，如有侵权，请联系网站删除：QQ：24596024