function [A,Q1] = qrit(A, nit)
%% [A,Q1] = qrit(A, nit)
%% nit iterations of QR... 
Q1 = eye(length(A));
for i=1:nit
  [Q,R] = qr(A); 
  A = R*Q 
  pause 
%%-------------------- accumulate the Q's
  Q1 = Q1* Q;
end
