qrmaxtr.m

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

function qrmaxtr(a,e)
k=1
n=length(a);
a0=zeros(n)
while norm(diag(a-a0))>e
    k=k+1
    a0=a;
    [q,r]=qr(a)
    a=r*q;
end
t=diag(a)