% SR_SIN0 Symbolic/variable-precision recurrence coefficients 
%    for the weight function w(x)=1+sin(1/x) on [0,1]
%
syms mom ab
N=40; 
%digits(70); dig=70;
%digits(95); dig=95;
digits(100); dig=100;
mom(1)='sin(1)-Ci(1)';
mom(2)='(Si(1)-pi/2+sin(1)+cos(1))/2';
for k=3:2*N
  mom(k)=(('cos(1)'-mom(k-2))/(k-1)+'sin(1)')/k;
end
for k=1:2*N
  mom(k)=mom(k)+1/k;
end
ab=schebyshev(dig,N,mom);
ab(:,1), ab(:,2)