32-digit values of the first 100 recurrence coefficients using the Bose-Einstein weight function: w(x)=[x/(e^x-1)]^4 computed by the SOPQ routine sr_boseeinstein(100,4,32)

<p>The first 100 recurrence coefficients for the Bose-Einstein weight function w(x)=[x/(e^x-1)]^4 on R_{+} are obtained to 32 decimal digits from the first 200 moments by using Chebyshev's algorithm in sufficiently high precision; cf. Sec. 4 of "Variable-precision recurrence coefficients for nonstandard orthogonal polynomials", Numerical Algorithms 52 (2009), 409-418.</p>