%% QR algoritmus

format long 
A=[1 2 3;4 5 6;7 8 9];
AA=A;
for k=1:5
    [Q,R]=qr(AA); 
    AA=R*Q
%    pause(2)
end
%AA

%
% Matice s vlastními čísly různými v absolutní hodnotě

format short
A=[11 -6 4 -2; 4 1 0 0; -9 9 -6 5; -6 6 -6 7]
eig(A)
AA=A;
for k=1:10
    [Q,R]=qr(AA);
    AA=R*Q;
end
AA

%
% Hessenbergův tvar

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

%
% Matice s násobným vlastním číslem

A=[4 3 3 -2; -2 5 -1 -2; 4 -1 3 2; 0 -1 -3 4]
eig(A)

%
% Matice s komplexními vlastními čísly



%
%