% QR algoritmus  - ukázka

A=[1 2 3; 4 5 6;7 8 9];

%

AA=A;
for k=1:5
    [Q,R]=qr(AA);
    AA=R*Q
    pause(3)
end

%

eig(A)
norm(ans-diag(AA))

%

format long
eig(A)-diag(AA)

%

for k=1:10
    [Q,R]=qr(AA);
    AA=R*Q;
end
eig(A)-diag(AA)