function monomial_to_gegenbauer_matrix_test ( )

%*****************************************************************************80
%
%% monomial_to_gegenbauer_matrix_test() tests monomial_to_gegenbauer_matrix().
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    24 April 2024
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'monomial_to_gegenbauer_matrix_test():\n' );
  fprintf ( 1, '  monomial_to_gegenbauer_matrix() evaluates the matrix\n' );
  fprintf ( 1, '  which converts monomial polynomial coefficients\n' );
  fprintf ( 1, '  to Gegenbauer coefficients.\n' );
  
  n = 5;
  alpha = 0.5;

  G = gegenbauer_to_monomial_matrix ( n, alpha );
  fprintf ( 1, '\n' );
  fprintf ( 1, '  Gegenbauer to Monomial matrix G:\n' );
  disp ( G );

  Ginv = monomial_to_gegenbauer_matrix ( n, alpha );
  fprintf ( 1, '\n' );
  fprintf ( 1, '  Monomial to Gegenbauer matrix Ginv:\n' );
  disp ( Ginv );

  I = G * Ginv;
  fprintf ( 1, '\n' );
  fprintf ( 1, '  I = G * Ginv:\n' );
  disp ( I );

  return
end

