function [A  Qp] = qrs2(A, nit)
%% [A  Qp] = qrs(A, nit)
%% complex shift
  %% very simple classical QR algorithm 
%% with shifts. 
%%  
Qp = eye(size(A));
n = size(A,1) ;
t = 0.0;
for i=1:nit 
  A = A - t * eye(n) ;
  [Q, R] = qr(A);
  A = R*Q + t*eye(n)
  Qp = Qp*Q;
  D=[A(n-1,n-1),A(n-1,n);A(n,n-1),A(n,n)];
%%-------------------- shift2 computes the shift
  t =shift2(D);
  pause 
end
