32-digit values of the first 100 recurrence coefficients for orthogonal polynomials relative to the weight function w(x)=const*x^d*exp(-b/x^a-c*x^a) on [0,Inf], const=exp(1)/(3*√2*π), a=1/3, b=c=1/2,d=-7/6, are computed by a moment-based method using the routine sr_exp_abcd(dig,32,100,1/3,1/2,1/2,-7/6), where dig=88 has been determined by the routine dig_exp_abcd(100,1/3,1/2,1/2,-7/6 ,80,4,32). The software provided in this dataset allows generating an arbitrary number N of recurrence coefficients for arbitrary a > 0, b > 0, c > 0, d real, as well as for different precisions.
Cite this work
Researchers should cite this work as follows:
- Gautschi, W. (2017). 32-digit values of the first 100 recurrence coefficients for a four-parameter exponential weight function with parameters a=1/3, b=c=1/2, d=-7/6. Purdue University Research Repository. doi:10.4231/R7RF5S1T
The dataset contains a text file and four Matlab routines.