function eig_ck ( )

%  Reference:
%
%    Ward Cheney & David Kincaid
%    Numerical Mathematics and Computing, Fifth Edition
%    Brooks/Cole Publishing Company, 2003.
%
% eigenvalues of matrix A
% via characteristic polynomial and directly
%
  fprintf ( '' )
  fprintf ( 'eig_ck():' )

  A = [3 2; 7 -2]
  p = poly(A)
  roots(p)
  eig(A)

  return
end

