% R_JACOBILOG  Recurrence coefficients for an algebraic/logarithmic
% weight function.
%
%    ab=r_jacobilog(n,a,b) generates the first n recurrence coefficients
%    for monic polynomials orthogonal with respect to the weight
%    function w(t)=(1-t)^a*t^b*log(1/t) on [0,1]. The call ab=r_jacobilog(n) 
%    is the same as ab=r_jacobilog(n,0,0) and ab=r_jacobilog(n,b) the same
%    as ab=r_jacobilog(n,0,b).
%
function ab=r_jacobilog(N,a,b)
if nargin<2, a=0; end; if nargin<3, b=a; end
if((N<=0)|(a<=-1)|(b<=-1)) error('parameter(s) out of range'), end
abj=r_jacobi01(2*N,a,b);
mom=mm_jacobilog(2*N,a,b);
ab=chebyshev(N,mom,abj);
