function gegenbauer_monomial_gegenbauer_test ( )

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

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

  gcoef = rand ( n + 1, 1 );
  mcoef = gegenbauer_to_monomial ( n, alpha, gcoef );
  gcoef2 = monomial_to_gegenbauer ( n, alpha, mcoef );

  e = norm ( gcoef - gcoef2 );
  fprintf ( 1, '\n' );
  fprintf ( 1, '  L2 difference = %g\n', e );

  return
end
