% SR_RMKLTm 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_RMKLTm(dig,N,a,b)
syms mom
dstr=num2str(dig);
am=vpa(-a,dig); b=vpa(b,dig);
maple(strcat('Digits:=',dstr));
amstr=char(am); bstr=char(b);
BIGstr=strcat('evalf(lnGAMMA(',amstr,'-',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)]