function padua_points_set_test ( )

%*****************************************************************************80
%
%% padua_points_set_test() tests padua_points_set().
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    31 August 2016
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'padua_points_set_test():\n' );
  fprintf ( 1, '  padua_points_set() looks the Padua points up in a table.\n' );
 
  for l = 3 : 4
    n = padua_order ( l );
    xy1 = padua_points ( l );
    xy2 = padua_points_set ( l );
    fprintf ( 1, '\n' );
    fprintf ( 1, '  Level %d  Padua points\n', l );
    fprintf ( 1, '\n' );
    for j = 1 : n
      fprintf ( 1, '  %4d  %14.6g  %14.6g\n', j, xy1(1,j), xy1(2,j) );
      fprintf ( 1, '        %14.6g  %14.6g\n',   xy2(1,j), xy2(2,j) );
    end
  end

  return
end
