%
 %  set up matrix
 %
  n = 40; k = 40;
 
  A = randn(n);
 %A = A + A';
 %
 %
  v =  randn(n,1); 
 [V,H,f] = ArnoldiC(A,k,v);

 ek = zeros(k,1); ek(k)= 1;
 Resid = norm(A*V - V*H -  f*ek')

 orthtestV = norm(eye(k) - V'*V)

 t = eig(A);
 plot(real(t),imag(t),'k+');
 hold

%for k = 4:4:n,
 h = n/2;
 for k = 1:h,
    s = eig(H(1:k,1:k));
    plot(real(s),imag(s),'c.');
    nvals = k
    pause
 end
 for k = h+1:2*h,
    s = eig(H(1:k,1:k));
    plot(real(s),imag(s),'m.');
    nvals = k
    pause
 end
 plot(real(t),imag(t),'bo');
 plot(real(t),imag(t),'bd');
 plot(real(t),imag(t),'b+');
 hold