function hexagon_lyness_rule_test ( )

%*****************************************************************************80
%
%% hexagon_lyness_rule_test() tests hexagon_lyness_rule().
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    22 April 2019
%
%  Author:
%
%    John Burkardt
%
  addpath ( '../hexagon_lyness_rule' );

  timestamp ( );
  fprintf ( 1, '\n' );
  fprintf ( 1, 'hexagon_lyness_test():\n' );
  fprintf ( 1, '  MATLAB/Octave version %s\n', version ( ) );
  fprintf ( 1, '  Test hexagon_lyness().\n' );

  hexagon_quadrature_test ( @hexagon_lyness_rule01, 'Lyness rule #01' );
  hexagon_quadrature_test ( @hexagon_lyness_rule02, 'Lyness rule #02' );
  hexagon_quadrature_test ( @hexagon_lyness_rule03, 'Lyness rule #03' );
  hexagon_quadrature_test ( @hexagon_lyness_rule04, 'Lyness rule #04' );
  hexagon_quadrature_test ( @hexagon_lyness_rule05, 'Lyness rule #05' );
  hexagon_quadrature_test ( @hexagon_lyness_rule06, 'Lyness rule #06' );
  hexagon_quadrature_test ( @hexagon_lyness_rule07, 'Lyness rule #07' );
  hexagon_quadrature_test ( @hexagon_lyness_rule08, 'Lyness rule #08' );
  hexagon_quadrature_test ( @hexagon_lyness_rule09, 'Lyness rule #09' );
  hexagon_quadrature_test ( @hexagon_lyness_rule10, 'Lyness rule #10' );
  hexagon_quadrature_test ( @hexagon_lyness_rule11, 'Lyness rule #11' );
  hexagon_quadrature_test ( @hexagon_lyness_rule12, 'Lyness rule #12' );
  hexagon_quadrature_test ( @hexagon_lyness_rule13, 'Lyness rule #13' );
%
%  Terminate.
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'hexagon_lyness_test():\n' );
  fprintf ( 1, '  Normal end of execution.\n' );
  fprintf ( 1, '\n' );
  timestamp ( );

  rmpath ( '../hexagon_lyness_rule' );

  return
end
function timestamp ( )

%*****************************************************************************80
%
%% timestamp() prints the current YMDHMS date as a timestamp.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    14 February 2003
%
%  Author:
%
%    John Burkardt
%
  t = now;
  c = datevec ( t );
  s = datestr ( c, 0 );
  fprintf ( 1, '%s\n', s );

  return
end

