function gegenbauer_to_monomial_test ( )

%*****************************************************************************80
%
%% gegenbauer_to_monomial_test() tests gegenbauer_to_monomial().
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    26 April 2024
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'gegenbauer_to_monomial_test ( )\n' );
  fprintf ( 1, '  gegenbauer_to_monomial() converts a\n' );
  fprintf ( 1, '  polynomial from Gegenbauer form\n' );
  fprintf ( 1, '  to monomial form.\n' );

  alpha = 0.5;
  fprintf ( 1, '\n' );
  fprintf ( 1, '  Gegenbauer parameter = %g\n', alpha );

  nmax = 6;

  fprintf ( 1, '\n' );
  fprintf ( 1, '           ' );
  for k = 0 : nmax
    fprintf ( 1, 'X^%d     ', k );
  end
  fprintf ( 1, '\n' );
  fprintf ( 1, '\n' );

  for n = 0 : nmax
    gcoef = zeros ( n + 1, 1 );
    gcoef(n+1) = 1.0;
    mcoef = gegenbauer_to_monomial ( n, alpha, gcoef );
    fprintf ( 1, 'H%d(x) = ', n );
    for k = 0 : n
      fprintf ( 1, '%8.3f', mcoef(k+1) );
    end
    fprintf ( 1, '\n' );
  end

  return
end
