32-digit values of the first 100 recurrence coefficients for orthogonal polynomials relative to the weight function w(x)=x^a*[log(1/x)]^b on [0,1], a=-1/2, b=3, are computed by a moment-based method using the routine sr_l_alglog(dig,32,100,-1/2,3), where dig=172 has been determined by the routine dig_l_alglog(100,-1/2,3,164,4,32). The software provided in this dataset allows generating an arbitrary number N of recurrence coefficients for arbitrary a > -1, b > -1, as well as for different precisions. If the singularities, with the same exponents, occur at the right endpoint, that is, if w(x)=(1-x)^a*[log(1/(1-x))]^b, then the alpha-coefficients must be replaced by 1 minus the present ones, whereas the beta-coefficients remain the same.
Cite this work
Researchers should cite this work as follows:
- Gautschi, W. (2016). 32-digit values of the first 100 recurrence coefficients for the weight function w(x)=x^(-1/2)*[log(1/x)]^3 on [0,1]. Purdue University Research Repository. doi:10.4231/R7XK8CH6
The dataset consists of one text file and four Matlab scripts.