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