function value = laguerre_integral_nd ( dim_num, expon )

%*****************************************************************************80
%
%% laguerre_integral_nd() evaluates a multidimensional Laguerre polynomial integral.
%
%  Discussion:
%
%    L(1,n) = Integral ( 0 <= x < Infinity ) x^n exp ( -x ) dx
%           = n!
%
%    L(d,n) = Integral ( 0 <= x(i) < Infinity ) 
%             x1^n1 * x2^n2...*xd^nd * exp(-x1-x2...-xd ) dx
%           = Product ( n(i)! ).
%
%  Licensing:
%
%    This code is distributed under the MIT license. 
%
%  Modified:
%
%    12 January 2021
%
%  Author:
%
%    John Burkardt
%
%  Input:
%
%    integer DIM_NUM, the spatial dimension.
%
%    integer EXPON(DIM_NUM), the exponents.
%
%  Output:
%
%    real VALUE, the value of the integral.
%
  value = 1.0;
  for dim = 1 : dim_num
    value = value * factorial ( expon(dim) ) ;
  end

  return
end
