function fejer1_set_test ( )

%*****************************************************************************80
%
%% fejer1_set_test() tests fejer1_set().
%
%  Licensing:
%
%    This code is distributed under the MIT license. 
%
%  Modified:
%
%    15 April 2015
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'fejer1_set_test():\n' );
  fprintf ( 1, '  fejer1_set() sets the abscissas and weights\n' );
  fprintf ( 1, '  of a Fejer type 1 quadrature rule.\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '    Order  W             X\n' );
  fprintf ( 1, '\n' );

  for n = 1 : 10

    [ x, w ] = fejer1_set ( n );

    fprintf ( 1, '\n' );
    fprintf ( 1, '  %8d\n', n );

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

  end

  return
end
