% SR_RMKLTp Symbolic/variable-precision calculation of the recurrence
%    coefficients for the weight function
%   
%      w(t)=t^{-a} (1+cos(b*ln(1/t)+g)), on [0,1].
%
% where g=Im(ln(Gamma(a+b*i))).
%
function ab=sr_RMKLTp(dig,N,a,b)
syms mom
dstr=num2str(dig);
a=vpa(a,dig); b=vpa(b,dig);
maple(strcat('Digits:=',dstr));
astr=char(a); bstr=char(b);
BIGstr=strcat('evalf(lnGAMMA(',astr,'+',bstr,'*i));');
gstr=maple(BIGstr);
g=imag(vpa(gstr,dig));
for k=1:2*N
  mom(k)=1/(k-a)+((k-a)*cos(g)-b*sin(g)) ...
  /((k-a)^2+b^2);
end
ab=schebyshev(dig,N,mom);
%[vpa(ab(N,1),16) vpa(ab(N,2),16)]