function padua_points_test ( )

%*****************************************************************************80
%
%% padua_points_test() tests padua_points().
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    09 June 2014
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'padua_points_test():\n' );
  fprintf ( 1, '  padua_points() returns the points of a Padua rule.\n' );

  for l = 0 : 10
    n = padua_order ( l );
    xy = padua_points ( l );
    label = sprintf ( '  Level %d Padua points:', l );
    r8mat_transpose_print ( 2, n, xy, label );
  end

  return
end
