% SR_EXPMINFPINF Symbolic/variable-precision recurrence
%    coefficients for the weight function 
%    w(x)=exp(-1/x^2-x^2) on [-inf,inf]
%
syms mom ab
N=40;
digits(50); dig=50;
%digits(55); dig=55;
%digits(60); dig=60;
mom(2:2:2*N)=vpa(0,dig);
mom(1)=vpa('2*BesselK(1/2,2)',dig);
mom(3)=3*mom(1)/2;
for k=3:2*N
  mom(2*k-1)=(k-3/2)*mom(2*k-3)+mom(2*k-5);
end
ab=schebyshev(N,mom,dig);
ab(:,1), ab(:,2)