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

<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], ε=.005, are
computed by a multicomponent discretization procedure using the
routine sr_OPhrbimod(32,100), with dig=34, epsi=.005 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 -.1 and .1 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>