function colleague_matrix_test ( )

  n = 3;
  p = [ ...
   5.625000000000000, ...
  -3.031250000000000, ...
   0.187499999999999, ...
   0.031249999999999 ]
%
%  p = T2(x) = 1 - 2x^2
%  r = -sqrt(2)/2, +sqrt(2)/2
%
  n = 2;
  p = [ ...
   0.0, ...
   0.0, ...
   1.0 ]

  C = colleague_matrix ( n, p );

  disp ( C )

  ev = eig ( C )
%
%  p = T3(x) = 4x^3 - 3x
%  r = 0, -sqrt(3)/2, +sqrt(3)/2
%
  n = 3;
  p = [ ...
   0.0, ...
   0.0, ...
   0.0, ...
   1.0 ]

  C = colleague_matrix ( n, p );

  disp ( C )

  ev = eig ( C )


  return
end

