function quadrature_weights_vandermonde_test01 ( )

%*****************************************************************************80
%
%% quadrature_weights_vandermonde_test01() tests quadrature_weights_vandermonde()
%
%  Discussion:
%
%    A Newton-Cotes rule is used.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    20 February 2014
%
%  Author:
%
%    John Burkardt
%
  a =  0.0;
  b = +1.0;
  n = 5;

  fprintf ( 1, '\n' );
  fprintf ( 1, 'quadrature_weights_vandermonde_test01():\n' );
  fprintf ( 1, '  Use the Vandermonde procedure to compute the\n' );
  fprintf ( 1, '  quadrature weights for a Newton-Cotes rule.\n' );
  fprintf ( 1, '  Order N = %d\n', n );
  fprintf ( 1, '  Interval = [%g,%g]\n', a, b );
%
%  Set the points.
%
  x = linspace ( a, b, n );
  r8vec_print ( n, x, '  Abscissas:' );
%
%  Compute the weights.
%
  w = quadrature_weights_vandermonde ( n, a, b, x );

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

  return
end
