function ab=sr_alg(dig,nofdig,N,a,b,c)
%SR_ALG Recurrence coefficients for a weight function on [0,1]
%with algebraic singularities at one or both endpoints.
%   AB=SR_ALG(DIG,NOFDIG,N,A,B,C) generates in DIG-digit
%   arithmetic the first N recurrence coefficients for the 
%   weight function w(x)=x^A*(1-x^C)^B on [0,1], A>-1, B>-1, 
%   C>0, and places them to NOFDIG digits into the Nx2 array AB.

digits(dig);
mom=smom_alg(dig,N,a,b,c);
ab=schebyshev(dig,N,mom');
ab=vpa(ab,nofdig);
