function tribonacci_roots_test ( )

%*****************************************************************************80
%
%% tribonacci_roots_test() tests tribonacci_roots().
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    09 May 2021
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'tribonacci_roots_test():\n' );
  fprintf ( 1, '  tribonacci_roots() computes the Tribonacci roots.\n' );
  fprintf ( 1, '\n' );

  [ alpha, beta, gamma ] = tribonacci_roots ( );

  p = alpha^3 - alpha^2 - alpha - 1.0;
  fprintf ( 1, '  alpha = %g, p(alpha) = %g\n', alpha, p );

  p = beta^3 - beta^2 - beta - 1.0;
  fprintf ( 1, '  beta = %g+%gi, p(beta) = %g\n', ...
    real ( beta ), imag ( beta ), p );

  p = gamma^3 - gamma^2 - gamma - 1.0;
  fprintf ( 1, '  gamma = %g+%gi, p(gamma) = %g\n', ...
    real ( gamma ), imag ( gamma ), p );

  return
end
