%
%  Drives the shifted QR iteration qriter
%
format short e
n = 25;
A = randn(n);

[V,H] = qriter(A);

s = eig(H); t = eig(A);
s = sort(abs(s)); t = sort(abs(t));

[t s t-s]