% zu Aufgabe 4.2: Eigenwerte einer symmetrischen Tridiagonalmatrix
n = 50;
A = diag(1:n) + diag(ones(n-1,1),1) + diag(ones(n-1,1),-1);
lambda = qriter(A, 1e-12)

disp('Kontrolle:')
min(eig(A))