def complex_i_test ( ):

#*****************************************************************************80
#
## COMPLEX_I_TEST tests COMPLEX_I.
#
#  Licensing:
#
#    This code is distributed under the GNU LGPL license.
#
#  Modified:
#
#    09 January 2015
#
#  Author:
#
#    John Burkardt
#
  from r8mat_print import r8mat_print

  print ''
  print 'COMPLEX_I_TEST'
  print '  COMPLEX_I computes the COMPLEX_I matrix.'

  m = 2
  n = m

  a = complex_i ( )
  r8mat_print ( m, n, a, '  COMPLEX_I matrix:' )

  print ''
  print 'COMPLEX_I_TEST'
  print '  Normal end of execution.'

  return

def rutis1_test ( ):

#*****************************************************************************80
#
## RUTIS1_TEST tests RUTIS1.
#
#  Licensing:
#
#    This code is distributed under the GNU LGPL license.
#
#  Modified:
#
#    24 January 2015
#
#  Author:
#
#    John Burkardt
#
  from r8mat_print import r8mat_print

  print ''
  print 'RUTIS1_TEST'
  print '  RUTIS1 computes the RUTIS1 matrix.'

  n = 4
  a = rutis1 ( )
  r8mat_print ( n, n, a, '  RUTIS1 matrix:' )

  print ''
  print 'RUTIS1_TEST'
  print '  Normal end of execution.'

  return

def r8mat_data_read_test ( ):

#*****************************************************************************80
#
## R8MAT_DATA_READ_TEST tests R8MAT_DATA_READ.
#
#  Licensing:
#
#    This code is distributed under the GNU LGPL license.
#
#  Modified:
#
#    03 December 2014
#
#  Author:
#
#    John Burkardt
#
  from r8mat_print import r8mat_print

  print ''
  print 'R8MAT_DATA_READ_TEST:'
  print '  Test R8MAT_DATA_READ, which reads data from an R8MAT.'

  m = 5
  n = 3
  filename = 'r8mat_write_test.txt'
  a = r8mat_data_read ( filename, m, n )
  r8mat_print ( m, n, a, '  Data read from file:' )

  print ''
  print 'R8MAT_DATA_READ_TEST:'
  print '  Normal end of execution.'

  return

def oto_test ( ):

#*****************************************************************************80
#
## OTO_TEST tests OTO.
#
#  Licensing:
#
#    This code is distributed under the GNU LGPL license.
#
#  Modified:
#
#    22 December 2014
#
#  Author:
#
#    John Burkardt
#
  from r8mat_print import r8mat_print

  print ''
  print 'OTO_TEST'
  print '  OTO computes the OTO matrix.'

  m = 5
  n = m
  a = oto ( m, n )
  r8mat_print ( m, n, a, '  OTO matrix:' )

  print ''
  print 'OTO_TEST'
  print '  Normal end of execution.'

  return

def fibonacci2_test ( ):

#*****************************************************************************80
#
## FIBONACCI2_TEST tests FIBONACCI2.
#
#  Licensing:
#
#    This code is distributed under the GNU LGPL license.
#
#  Modified:
#
#    28 January 2015
#
#  Author:
#
#    John Burkardt
#
  from r8mat_print import r8mat_print

  print ''
  print 'FIBONACCI2_TEST'
  print '  FIBONACCI2 computes the FIBONACCI2 matrix.'

  m = 5
  n = m

  a = fibonacci2 ( n )
  r8mat_print ( m, n, a, '  FIBONACCI2 matrix:' )

  print ''
  print 'FIBONACCI2_TEST'
  print '  Normal end of execution.'

  return

def sylvester_kac_test ( ):

#*****************************************************************************80
#
## SYLVESTER_KAC_TEST tests SYLVESTER_KAC.
#
#  Licensing:
#
#    This code is distributed under the GNU LGPL license.
#
#  Modified:
#
#    13 April 2015
#
#  Author:
#
#    John Burkardt
#
  from r8mat_print import r8mat_print

  print ''
  print 'SYLVESTER_KAC_TEST'
  print '  SYLVESTER_KAC computes the SYLVESTER_KAC matrix.'

  m = 5
  n = m
  a = sylvester_kac ( n )
  r8mat_print ( m, n, a, '  SYLVESTER_KAC matrix:' )

  print ''
  print 'SYLVESTER_KAC_TEST'
  print '  Normal end of execution.'

  return

def rosser1_test ( ):

#*****************************************************************************80
#
## ROSSER1_TEST tests ROSSER1.
#
#  Licensing:
#
#    This code is distributed under the GNU LGPL license.
#
#  Modified:
#
#    10 February 2015
#
#  Author:
#
#    John Burkardt
#
  from r8mat_print import r8mat_print

  print ''
  print 'ROSSER1_TEST'
  print '  ROSSER1 computes the ROSSER1 matrix.'

  n = 8
  a = rosser1 ( )
  r8mat_print ( n, n, a, '  ROSSER1 matrix:' )

  print ''
  print 'ROSSER1_TEST'
  print '  Normal end of execution.'

  return

def permutation_random_test ( ):

#*****************************************************************************80
#
## PERMUTATION_RANDOM_TEST tests PERMUTATION_RANDOM.
#
#  Licensing:
#
#    This code is distributed under the GNU LGPL license.
#
#  Modified:
#
#    27 March 2015
#
#  Author:
#
#    John Burkardt
#
  from r8mat_print import r8mat_print

  print ''
  print 'PERMUTATION_RANDOM_TEST'
  print '  PERMUTATION_RANDOM computes the PERMUTATION_RANDOM matrix.'

  n = 5

  key = 123456789
  a = permutation_random ( n, key )
  r8mat_print ( n, n, a, '  PERMUTATION_RANDOM matrix:' )

  print ''
  print 'PERMUTATION_RANDOM_TEST'
  print '  Normal end of execution.'

  return

def orth_symm_test ( ):

#*****************************************************************************80
#
## ORTH_SYMM_TEST tests ORTH_SYMM.
#
#  Licensing:
#
#    This code is distributed under the GNU LGPL license.
#
#  Modified:
#
#    05 February 2015
#
#  Author:
#
#    John Burkardt
#
  from r8mat_print import r8mat_print

  print ''
  print 'ORTH_SYMM_TEST'
  print '  ORTH_SYMM computes the ORTH_SYMM matrix.'

  m = 5
  n = 5
  a = orth_symm ( n )
  r8mat_print ( m, n, a, '  ORTH_SYMM matrix:' )

  print ''
  print 'ORTH_SYMM_TEST'
  print '  Normal end of execution.'

  return

def eulerian_test ( ):

#*****************************************************************************80
#
## EULERIAN_TEST tests EULERIAN.
#
#  Licensing:
#
#    This code is distributed under the GNU LGPL license.
#
#  Modified:
#
#    25 January 2015
#
#  Author:
#
#    John Burkardt
#
  from r8mat_print import r8mat_print

  print ''
  print 'EULERIAN_TEST'
  print '  EULERIAN computes the EULERIAN matrix.'

  m = 4
  n = m

  a = eulerian ( m, n )
  r8mat_print ( m, n, a, '  EULERIAN matrix:' )

  print ''
  print 'EULERIAN_TEST'
  print '  Normal end of execution.'

  return

def maxij_test():

    # *****************************************************************************80
    #
    ## MAXIJ_TEST tests MAXIJ.
    #
    #  Licensing:
    #
    #    This code is distributed under the GNU LGPL license.
    #
    #  Modified:
    #
    #    26 January 2015
    #
    #  Author:
    #
    #    John Burkardt
    #
    from r8mat_print import r8mat_print

    print ""
    print "MAXIJ_TEST"
    print "  MAXIJ computes the MAXIJ matrix."

    m = 5
    n = 5
    a = maxij(m, n)
    r8mat_print(m, n, a, "  MAXIJ matrix:")

    print ""
    print "MAXIJ_TEST"
    print "  Normal end of execution."

    return

def moler3_test ( ):

#*****************************************************************************80
#
## MOLER3_TEST tests MOLER3.
#
#  Licensing:
#
#    This code is distributed under the GNU LGPL license.
#
#  Modified:
#
#    26 January 2015
#
#  Author:
#
#    John Burkardt
#
  from r8mat_print import r8mat_print

  print ''
  print 'MOLER3_TEST'
  print '  MOLER3 computes the MOLER3 matrix.'

  m = 5
  n = m

  a = moler3 ( m, n )
  r8mat_print ( m, n, a, '  MOLER3 matrix:' )

  print ''
  print 'MOLER3_TEST'
  print '  Normal end of execution.'

  return

def ring_adj_test ( ):

#*****************************************************************************80
#
## RING_ADJ_TEST tests RING_ADJ.
#
#  Licensing:
#
#    This code is distributed under the GNU LGPL license.
#
#  Modified:
#
#    12 March 2015
#
#  Author:
#
#    John Burkardt
#
  from r8mat_print import r8mat_print

  print ''
  print 'RING_ADJ_TEST'
  print '  RING_ADJ computes the RING_ADJ matrix.'

  m = 6
  n = m

  a = ring_adj ( m, n )
  r8mat_print ( m, n, a, '  RING_ADJ matrix:' )

  print ''
  print 'RING_ADJ_TEST'
  print '  Normal end of execution.'

  return

def one_test ( ):

#*****************************************************************************80
#
## ONE_TEST tests ONE.
#
#  Licensing:
#
#    This code is distributed under the GNU LGPL license.
#
#  Modified:
#
#    20 February 2015
#
#  Author:
#
#    John Burkardt
#
  from r8mat_print import r8mat_print

  print ''
  print 'ONE_TEST'
  print '  ONE computes the ONE matrix.'

  m = 4
  n = m

  a = one ( m, n )
  r8mat_print ( m, n, a, '  ONE matrix:' )

  print ''
  print 'ONE_TEST'
  print '  Normal end of execution.'

  return

def conex4_test():

    # *****************************************************************************80
    #
    ## CONEX4_TEST tests CONEX4.
    #
    #  Licensing:
    #
    #    This code is distributed under the GNU LGPL license.
    #
    #  Modified:
    #
    #    20 January 2015
    #
    #  Author:
    #
    #    John Burkardt
    #
    from r8mat_print import r8mat_print

    print ""
    print "CONEX4_TEST"
    print "  CONEX4 computes the CONEX4 matrix."

    m = 4
    n = m

    a = conex4()
    r8mat_print(m, n, a, "  CONEX4 matrix:")

    print ""
    print "CONEX4_TEST"
    print "  Normal end of execution."

    return

def wilson_test ( ):

#*****************************************************************************80
#
## WILSON_TEST tests WILSON.
#
#  Licensing:
#
#    This code is distributed under the GNU LGPL license.
#
#  Modified:
#
#    20 December 2014
#
#  Author:
#
#    John Burkardt
#
  from r8mat_print import r8mat_print

  print ''
  print 'WILSON_TEST'
  print '  WILSON computes the WILSON matrix.'

  n = 4
  a = wilson ( )
  r8mat_print ( n, n, a, '  WILSON matrix:' )

  print ''
  print 'WILSON_TEST'
  print '  Normal end of execution.'

  return

def conference_test ( ):

#*****************************************************************************80
#
## CONFERENCE_TEST tests CONFERENCE.
#
#  Licensing:
#
#    This code is distributed under the GNU LGPL license.
#
#  Modified:
#
#    19 April 2015
#
#  Author:
#
#    John Burkardt
#
  from r8mat_print import r8mat_print

  print ''
  print 'CONFERENCE_TEST'
  print '  CONFERENCE computes the CONFERENCE matrix.'
#
#  Note that N-1 must be an odd prime or a power of an odd prime.
#
  n = 6
  a = conference ( n )
  r8mat_print ( n, n, a, '  CONFERENCE matrix:' )

  print ''
  print 'CONFERENCE_TEST'
  print '  Normal end of execution.'

  return

def boothroyd_test ( ):

#*****************************************************************************80
#
## BOOTHROYD_TEST tests BOOTHROYD.
#
#  Licensing:
#
#    This code is distributed under the GNU LGPL license.
#
#  Modified:
#
#    27 December 2014
#
#  Author:
#
#    John Burkardt
#
  from r8mat_print import r8mat_print

  print ''
  print 'BOOTHROYD_TEST'
  print '  BOOTHROYD computes the BOOTHROYD matrix.'

  m = 5
  n = m
  a = boothroyd ( n )
  r8mat_print ( m, n, a, '  BOOTHROYD matrix:' )

  print ''
  print 'BOOTHROYD_TEST'
  print '  Normal end of execution.'

  return

def harman_test ( ):

#*****************************************************************************80
#
## HARMAN_TEST tests HARMAN.
#
#  Licensing:
#
#    This code is distributed under the GNU LGPL license.
#
#  Modified:
#
#    10 February 2015
#
#  Author:
#
#    John Burkardt
#
  from r8mat_print import r8mat_print

  print ''
  print 'HARMAN_TEST'
  print '  HARMAN computes the HARMAN matrix.'

  n = 8
  a = harman ( )
  r8mat_print ( n, n, a, '  HARMAN matrix:' )

  print ''
  print 'HARMAN_TEST'
  print '  Normal end of execution.'

  return

def summation_test ( ):

#*****************************************************************************80
#
## SUMMATION_TEST tests SUMMATION.
#
#  Licensing:
#
#    This code is distributed under the GNU LGPL license.
#
#  Modified:
#
#    17 December 2014
#
#  Author:
#
#    John Burkardt
#
  from r8mat_print import r8mat_print

  print ''
  print 'SUMMATION_TEST'
  print '  SUMMATION computes the SUMMATION matrix.'

  m = 5
  n = 4
  a = summation ( m, n )
  r8mat_print ( m, n, a, '  SUMMATION matrix:' )

  print ''
  print 'SUMMATION_TEST'
  print '  Normal end of execution.'

  return