function determ = ijfact1_determinant ( n )

%*****************************************************************************80
%
%% ijfact1_determinant() returns the determinant of the IJFACT1 matrix.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    12 January 2021
%
%  Author:
%
%    John Burkardt
%
%  Input:
%
%    integer N, the order of the matrix.
%
%  Output:
%
%    real DETERM, the determinant.
%
  determ = 1.0;

  for i = 1 : n - 1
    determ = determ * factorial ( i + 1 ) * factorial ( n - i );
  end

  determ = determ * factorial ( n + 1 );

  return
end
