function gegenbauer_ek_compute_test ( )

%*****************************************************************************80
%
%% gegenbauer_ek_compute_test() tests gegenbauer_ek_compute().
%
%  Licensing:
%
%    This code is distributed under the MIT license. 
%
%  Modified:
%
%    19 November 2015
%
%  Author:
%
%    John Burkardt
%
  alpha = 0.5;
  a = -1.0;
  b = +1.0;

  fprintf ( 1, '\n' );
  fprintf ( 1, 'gegenbauer_ek_compute_test():\n' );
  fprintf ( 1, '  gegenbauer_ek_compute() computes Gauss-Gegenbauer rules;\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '  Abscissas and weights for a generalized Gauss Gegenbauer rule\n' );
  fprintf ( 1, '  with ALPHA = %f\n', alpha );
  fprintf ( 1, '  Integration interval is [%g,%g]\n', a, b );

  for n = 1 : 10

    [ x, w ] = gegenbauer_ek_compute ( n, alpha, a, b );

    fprintf ( 1, '\n' );

    for i = 1 : n
      fprintf ( 1, '  %2d  %24.16g  %24.16g\n', i, w(i), x(i) );
    end

  end

  return
end

