function hexagon_integrals_test01 ( )

%*****************************************************************************80
%
%% hexagon_integrals_test01 requests integrals over the unit hexagon.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    23 April 2019
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'hexagon_integrals_test01\n' );
  fprintf ( 1, '  Test hexagon01_integral, which computes\n' );
  fprintf ( 1, '  monomial integrals x^p y^q over the unit hexagon.\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '   P   Q             Nu(P,Q)\n' );
  fprintf ( 1, '\n' );
  k = 0;
  for s = 0 : 10
    for p = s : -1 : 0
      q = s - p;
      if ( mod ( p, 2 ) == 0 && mod ( q, 2 ) == 0 )
        k = k + 1;
        nu_pq = hexagon01_integral ( p, q );
        fprintf ( 1, '  %2d  %2d  %14.6g\n', p, q, nu_pq );
      end
    end
  end

  return
end
