function gegenbauer_to_monomial_matrix_test ( )

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

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

  for i = 0 : n - 1
    g = zeros ( n, 1 );
    g(i+1) = 1.0;
    mono = A * g;
    label = sprintf ( '  Monomial form of Gegenbauer polynomial #%d', i );
    r8poly_print ( mono, label );
  end

  return
end

