function lambda = qriter(A, epsilon)
%  Berechnung eines Eigenwertes der symmetrischen tridiagonalen
%  Matrix A mit Hilfe des QR-Verfahrens

n = size(A,1);
while abs(A(n,n-1)) > epsilon
    [Q,R] = qr(A);
    A = R*Q;
end
lambda = A(n,n);