program main

!*****************************************************************************80
!
!! latin_random_test() tests latin_random().
!
!  Licensing:
!
!    This code is distributed under the MIT license.
!
!  Modified:
!
!    12 November 2014
!
!  Author:
!
!    John Burkardt
!
  implicit none

  integer seed
  integer test

  call timestamp ( )
  write ( *, '(a)' ) ' '
  write ( *, '(a)' ) 'latin_random_test():'
  write ( *, '(a)' ) '  FORTRAN90 version'
  write ( *, '(a)' ) '  Test the LATIN_RANDOM library.'

  seed = 123456789

  do test = 1, 3

    call test01 ( seed )

  end do
!
!  Terminate.
!
  write ( *, '(a)' ) ' '
  write ( *, '(a)' ) 'LATIN_RANDOM_TEST:'
  write ( *, '(a)' ) '  Normal end of execution.'
  write ( *, '(a)' ) ' '
  call timestamp ( )

  stop 0
end
subroutine test01 ( seed )

!*****************************************************************************80
!
!! TEST01 tests LATIN_RANDOM.
!
!  Licensing:
!
!    This code is distributed under the MIT license.
!
!  Modified:
!
!    12 November 2014
!
!  Author:
!
!    John Burkardt
!
!  Parameters:
!
!    Input/output, integer SEED, a seed for the random number
!    generator.
!
  implicit none

  integer, parameter :: rk = kind ( 1.0D+00 )

  integer, parameter :: dim_num = 2
  integer, parameter :: point_num = 10

  integer seed
  real ( kind = rk ) x(dim_num,point_num)

  write ( *, '(a)' ) ' '
  write ( *, '(a)' ) 'TEST01'
  write ( *, '(a)' ) '  LATIN_RANDOM chooses a random Latin Square'
  write ( *, '(a)' ) '  cell arrangement, and then returns'
  write ( *, '(a)' ) '  a random point from each cell.'
  write ( *, '(a)' ) ' '
  write ( *, '(a,i6)'  ) '  Spatial dimension =  ', dim_num
  write ( *, '(a,i6)'  ) '  Number of points =   ', point_num
  write ( *, '(a,i12)' ) '  Random number SEED = ', seed

  call latin_random ( dim_num, point_num, seed, x )

  call r8mat_transpose_print ( dim_num, point_num, x, &
    '  The Latin Random points:' )

  return
end