% SR_EXP0 Symbolic/variable-precision recurrence coefficients 
%    for the weight function w(x)=exp(-1/x) on [0,1]
%
syms mom ab
N=40; 
digits(70); dig=70;
%digits(95); dig=95;
%digits(100); dig=100;
mom(1)=vpa('Ei(2,1)',dig);
for k=2:2*N
  mom(k)=('exp(-1)'-mom(k-1))/k;
end
ab=schebyshev(N,mom,dig);
ab(:,1), ab(:,2)