function value = hexagon01_monomial_integral ( p, q )

%*****************************************************************************80
%
%% hexagon01_monomial_integral() evaluates a monomial integral in the unit hexagon.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    23 April 2019
%
%  Author:
%
%    John Burkardt
%
%  Input:
%
%    integer P, Q, the indices of the moment.
%
%  Output:
%
%    real VALUE, the integral of x^p y^q.
%
  if ( mod ( p, 2 ) == 1 || mod ( q, 2 ) == 1 )
    value = 0.0;
  else
    n = 6;
    [ x, y ] = hexagon_unit_vertices ( );
    value = moment_unnormalized ( n, x, y, p, q );
  end

  return
end
