%
%  Drives the implicit polynomial restart QR iteration  qriterRS
%   for symmetric matrix A
%
n = 100;
A = randn(n);
A = A+A';

[V,H] = qriterRS(A);

%s = eig(H); t = eig(A);
%s = sort(abs(s)); t = sort(abs(t));
%
%[t s t-s]