function value = quadrilateral_unit_monomial_integral ( expon )

%*****************************************************************************80
%
%% quadrilateral_unit_monomial_integral(): monomial integral in a unit quadrilateral.
%
%  Discussion:
%
%    This routine returns the integral of
%
%      product ( 1 <= I <= 2 ) X(I)^EXPON(I)
%
%    over the unit quadrilateral.
%
%    The unit quadrilateral Q is bounded by the vertices
%    (0,0), (1,0), (1,1), (0,1).
% 
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    24 May 2023
%
%  Author:
%
%    John Burkardt
%
%  Input:
%
%    integer EXPON(2): the exponents.
%
%  Output:
%
%    real VALUE: the integral of the monomial.
%
  value = 1.0 / ( expon(1) + 1 ) ...
        * 1.0 / ( expon(2) + 1 );

  return
end
