function ab=sr_modbess(dig,nofdig,N,nu)
%SR_MODBESS Recurrence coefficients for a modified Bessel function.
%   SR_MODBESS(DIG,NOFDIG,N,NU) generates in DIG-digit arithmetic
%   the first N recurrence coefficients for the weight function
%   W(T)=C*K_NU(T) on [0,Inf], C=COS(NU*PI/2)/PI, NU>0 not an 
%   integer, and places them to NOFDIG digits into the Nx2 array AB.

digits(dig);
mom=mommodbess(dig,N,nu);
ab=schebyshev(dig,N,mom');
ab=vpa(ab,nofdig);
