function x = latin_center ( dim_num, point_num )

%*****************************************************************************80
%
%% latin_center() returns center points in a Latin square.
%
%  Discussion:
%
%    In each spatial dimension, there will be exactly one
%    point with the coordinate value
%
%      ( 1, 3, 5, ..., 2*point_num-1 ) / ( 2 * point_num )
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    21 March 2003
%
%  Author:
%
%    John Burkardt
%
%  Input:
%
%    integer dim_num, the spatial dimension.
%
%    integer point_num, the number of points.
%
%  Output:
%
%    real x(dim_num,point_num), the points.
%
  base = 1;

  for i = 1: dim_num

    perm = perm_uniform ( point_num, base );

    for j = 1: point_num
      x(i,j) = ( 2 * perm(j) - 1 ) / ( 2 * point_num );
    end 

  end

  return
end
