32-digit values of the first 100 recurrence coefficients for the Bose-Einstein-type weight function with exponent 4

<p>32-digit values of the first 100 recurrence coefficients for orthogonal polynomials relative to the weight function w(x)=[x/(exp(x)-1)]^r on [0,Inf], r=4, are computed by a moment-based method using the routine sr_boseeinstein(dig,32,100,4), where dig=224 has been determined by the routine dig_boseeinstein(100,4,216,4,32). For the respective moments, see Section 4 in Walter Gautschi, "Variable-precision recurrence coefficients for nonstandard orthogonal polynomials", Numerical Algorithms 52 (2009), 409-418. doi: 10.1007/s11075-009-9283-2. The software provided in this dataset allows generating an arbitrary number N of recurrence coefficients for arbitrary integer r > 0 as well as for different precisions.</p>