32-digit values of the first 100 recurrence coefficients for orthogonal polynomials relative to the weight function w(x)=(x2)*exp[-(x 2-1) 2/(4*ε)] on [0,Inf], ε=.001, are computed by a multicomponent discretization procedure using the routine sr_OPhrbimod(32,100), with dig=34, epsi=.001 entered at the prompt. The value dig=34 has been determined by the routine dig_sOPbimod(100,32,2,32), attesting to the high stability of the procedure. (Both routines may take several hours to run.) The auxiliary routine xhrbimod.m uses Newton's method to compute a zero x*=sqrt(1+u*) of the equation w(x)=1/2, given a close estimate of u*. There are two such zeros, x_1* and x_2*, x_1* < x_2*, and corresponding estimates -.05 and .05 for u*. The software provided in this dataset allows generating an arbitrary number N of recurrence coefficients for arbitrary ε > 0 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 the half-range bimodal weight function with parameter ε=.001. Purdue University Research Repository. doi:10.4231/R78K772F
The dataset consists of one text file and eight Matlab scripts.