32-digit values of the first 65 recurrence coefficients for orthogonal polynomials relative to the weight function w(x)=x^a*exp(-x)*log(1+1/x) on [1,Inf], a=0, are computed by a 1-component discretization method using the routine sr_explog1inf(32,65), where dig=36 has to be entered when prompted by the routine. The value 36 of dig has been determined by the routine dig_explog1inf(65,32,4,32). Both routines, dig_explog1inf.m and sr_explog1inf.m, may take many hours to run, the former as much as 19 hours, the latter about 12 hours. The software provided in this dataset allows generating any number N of recurrence coefficients for any a > -1 not an integer, as well as for different precisions, as long as there is no error message "Mcap exceeds Mmax in smcdis with irout=1", indicating lack of convergence in the discretization procedure.

Researchers should cite this work as follows:

• Gautschi, W. (2017). 32-digit values of the first 65 recurrence coefficients for the weight function w(x)=exp(-x)*log(1+1/x) on [1,Inf]. Purdue University Research Repository. doi:10.4231/R7028PJT

  BibTex | EndNote

1. Computer Science
2. Computer Science (Dept)
3. Mathematics
4. Modification algorithms for orthogonal polynomials
5. Orthogonal polynomials
6. Walter Gautschi Archives
7. weight functions

The dataset consists of one text file and eight Matlab scripts.