function ab=sr_theodorus(dig,nofdig,N)
%SR_THEODORUS  Recurrence coefficients for Theodorus polynomials
%   AB=SR_THEODORUS(DIG,NOFDIG,N) generates in DIG-digit arithmetic
%   the first N recurrence coefficients for the weight function
%   w(x)=x^(1/2)/(exp(x)-1) on [0,Inf] and places them to NOFDIG
%   digits into the Nx2 array AB.

digits(dig);
mom=momtheodorus(dig,N);
ab=schebyshev(dig,N,mom);
ab=vpa(ab,nofdig);
