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

<p>32-digit values of the first 100 recurrence coefficients for orthogonal polynomials relative to the weight function w(x)=(x<sup class="moz-txt-sup"><span style="display:inline-block;width:0;height:0;overflow:hidden">^</span>2</sup>)*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 [0,Inf], ε=.02, are computed by a multicomponent discretization procedure using the routine sr_OPhrbimod(32,100), with dig=34, epsi=.02 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 -.2 and .3 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.</p>