% RADAU Gauss-Radau quadrature rule.
%
% Given a weight function w encoded by the (n+1)x2 array ab of
% the first n+1 recurrence coefficients for the associated
% orthogonal polynomials, the first column of ab containing the
% n+1 alpha-coefficients and the second column the n+1 beta-
% coefficients, the call xw=RADAU(n,ab,end0) generates the
% nodes and weights xw of the (n+1)-point Gauss-Radau
% quadrature rule for the weight function w having a prescribed
% node end0 (typically at one of the end points of the support
% interval of w, or outside thereof). The n+1 nodes, in
% increasing order, are stored in the first column, the n+1
% corresponding weights in the second column, of the (n+1)x2
% array xw.
%
% For Jacobi weight functions, see also RADAU_JACOBI.
%
function xw=radau(N,ab,end0)
N0=size(ab,1); if N0