%function QR
%程序C7_3  QR方法
n=6;
i=1;
A=[5.3,-5.1043,2.46195,2.50473,-1.45763,-0.79176;...
       -10.2718,19.5001,-12.7062,-3.59933,3.93381,4.93779;...
       0.0,-11.3363,4.33072,-0.04125,-1.10558,-0.86184;...
       0.0,0.0,1.98931,-3.81215,-54.7416,0.32433;...
       0.0,0.0,0.0,2.46643,3.06251,3.91227;...
       0.0,0.0,0.0,0.0,-3.00331,2.91877];
[M,f1,eig1]=householder(A,n);
while i<=n-1
    i=i+1;
    A=M;
    [M,f1,eig1]=householder(A,n);
end
M
eig1
[K,L]=givens(M,n,eig1)
QRmethod(M,n)
while i<=n
   i=i+1;
   [M,f1,eig1]=householder(M(:,i),n);
end
M