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

%% newton_ab returns a Newton quadrature rule for [a,b].
%
  [ xi, rho] = newton ( n );

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

  return
end
