function quadrature_weights_vandermonde_2d_test02 ( )

%*****************************************************************************80
%
%% quadrature_weights_vandermonde_2d_test02 tests quadrature_weights_vandermonde_2d for Padua points.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    29 May 2014
%
%  Author:
%
%    John Burkardt
%
  a = -1.0;
  b = +1.0;
  c = -1.0;
  d = +1.0;

  fprintf ( 1, '\n' );
  fprintf ( 1, 'quadrature_weights_vandermonde_2d_test02:\n' );
  fprintf ( 1, '  Compute the weights associated with an interpolatory\n' );
  fprintf ( 1, '  quadrature rule defined by N=(T+1)*(T+2)/2 points,\n' );
  fprintf ( 1, '  exact for polynomials of total degree T or less.\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '  X Interval = [%g,%g]\n', a, b );
  fprintf ( 1, '  Y Interval = [%g,%g]\n', c, d );

  for t = 0 : 10

    n = ( ( t + 1 ) * ( t + 2 ) ) / 2;

    fprintf ( 1, '\n' );
    fprintf ( 1, '  Degree T = %d\n', t );
    fprintf ( 1, '  Number of points N = %d\n', n );

    [ x, y ] = padua_point_set ( t );
%
%  Compute the weights.
%
    w = quadrature_weights_vandermonde_2d ( t, n, a, b, c, d, x, y );

    r8vec_print_16 ( n, w, '  Weights:' );

  end

  return
end
