function ab=sr_jacalg(dig,N,a,b,c)
%SR_JACALG  Recurrence coefficients for polynomials orthogonal
%relative to a weight function having algebraic singularities at
%one or both endpoints.
%   AB=SR_JACALG(DIG,N,A,B,C) generates in DIG-digit arithmetic
%   the Nx2 array AB of the first 100 recurrence coefficients for
%   the weight function w(x)=x^A*(1-x^C)^B on [0,1], A>-1, B>-1,
%   C>0.

digits(dig);
mom=smom_jacalg(dig,2*N,a,b,c);
ab=schebyshev(dig,N,mom');
