raymaxtr.m

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

function raymaxtr(a,esp)
n=length(a);
x0=diag(ones(n));
k=1
x=a*x0;
y=x
ry=y'*x
while norm(x-x0)>esp
    k=k+1
    q=x
    y=x/norm(x,2)
    x=a*y
    x0=q
    ry=y'*x
end