32-digit values of the first 100 recurrence coefficients for orthogonal polynomials relative to the weight function w(x)=x^a*exp(-x)*(x-1-log(x)) on [0,Inf], a=0, are computed by a modified-moment-based method using the routine sr_laglogmm(dig,32,100,0), where dig=120 has been determined by the routine dig_laglogmm(100,0,112,4,32). For the modified moments, see Section 2.2 in Walter Gautschi, "Gauss quadrature routines for two classes of logarithmic weight functions", Numerical Algorithms 55 (2010), 265-277, doi: 10.1007/s11075-010-9366-0 and Section 1 for an application. The software provided in this dataset allows generating an arbitrary number N of recurrence coefficients for arbitrary a > -1 as well as for different precisions. See also doi:10.4231/R7T151N8 for a method based on ordinary moments. Both methods are equally efficient in terms of run time and produce the same answers.
Cite this work
Researchers should cite this work as follows:
- Gautschi, W. (2016). 32-digit values of the first 100 recurrence coefficients, obtained from modified moments, for the Laguerre weight function multiplied by a logarithmically singular function. Purdue University Research Repository. doi:10.4231/R7M043CX
The dataset consists of one text file and five Matlab scripts.