32-digit values of the first 100 recurrence coefficients for the bimodal weight function with parameter ε=.1

<p>32-digit values of the first 100 recurrence coefficients for orthogonal polynomials relative to the weight function w(x)=exp[-(x<sup class="moz-txt-sup"><span style="display:inline-block;width:0;height:0;overflow:hidden">^</span>2</sup>-1)<sup class="moz-txt-sup"><span style="display:inline-block;width:0;height:0;overflow:hidden">^</span>2</sup>/(4*ε)] on [-Inf,Inf], ε=.1, are computed by a multicomponent discretization procedure using the routine sr_OPbimod(32,100) with dig=34, epsi=.1 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.) For details, see Exercise 2.38(c) in Walter Gautschi, "Orthogonal polynomials in MATLAB: Exercises and solutions", SIAM, Philadelphia, PA (2016). Auxiliary routines muOPbimod_gp.m and explore_mu.m are intended to help determine suitable values for the parameter μ needed when ε < 1/10. (This requires a minor temporary change in the routine smcdis.m as explained on p.130 of the cited reference.) The value of μ should be taken to be at least equal to 1; when N = 100, other selected values of μ for ε = .008:-.001:.001 are found to be μ = 1.7, 4.7, 8.0, 15.3, 26.2, 44.4, 84.2, 201.3. The software provided in this dataset allows generating an arbitrary number N of recurrence coefficients for arbitrary ε > 0 as well as for different precisions.</p>