function [ x, w ] = gauss_ab ( n, a, b )

%*****************************************************************************80
%
%% gauss_ab returns a Gauss quadrature rule for [a,b].
%
%  Input:
%
%    integer n: the number of points in the rule.
%
%    real a, b: the endpoints of the interval.
%
%  Output:
%
%    real x(n), w(n): the points and weights of the rule.
%
  [ xi, rho] = gauss_m1p1 ( n );

  x = a + ( xi + 1.0 ) * ( b - a ) / 2.0;
  w = ( b - a ) * rho / 2.0;

  return
end
