A = rand(5,5);
A = A'*A;

B = A;

for i = 1:100
    [Q,R] = qr(B); B = R*Q;
end

diag(B)
eig(A)

% ans =
% 
%     7.1427
%     0.5986
%     0.2884
%     0.1562
%     0.0531
% 
% 
% ans =
% 
%     0.0531
%     0.1562
%     0.2884
%     0.5986
%     7.1427
