32-digit values of the first 100 recurrence coefficients for the weight function w(x)=x^(-1/2)*exp(-x)*(x-1-log(x)) on [0,Inf] obtained from modified moments

Listed in Datasets

By Walter Gautschi

Purdue University

32-digit values of the first 100 recurrence coefficients for the weight function w(x)=x^a*exp(-x)*(x-1-log(x)) on [0,Inf], a=-1/2

Version 1.0 - published on 05 Dec 2016 doi:10.4231/R7SQ8XDM - cite this Archived on 06 Jan 2017

Licensed under Attribution 3.0 Unported

Description

32-digit values of the first 100 recurrence coefficients for orthogonal polynomials relative to the weight function w(x)=x^a*exp(-x)*(x-1-log(x)) on [0,Inf], a=-1/2, are computed by a modified-moment-based method using the routine sr_laglogmm(dig,32,100,-1/2), where dig=120 has been determined by the routine dig_laglogmm(100,-1/2,112,4,32). For the modified moments, see Section 2.2 in Walter Gautschi, "Gauss quadrature routines for two classes of logarithmic weight functions", Numerical Algorithms 55 (2010), 265-277, doi: 10.1007/s11075-010-9366-0, and Section 1 for an application. The software provided in this dataset allows generating an arbitrary number N of recurrence coefficients for arbitrary a > -1 as well as for different precisions. See also doi:10.4231/R79P2ZMR for a method based on ordinary moments. Both methods produce the same answers, but the current one takes more than twice as long to run.

Cite this work

Tags

Notes

The dataset consists of one text file and five Matlab scripts.

The Purdue University Research Repository (PURR) is a university core research facility provided by the Purdue University Libraries, the Office of the Executive Vice President for Research and Partnerships, and Information Technology at Purdue (ITaP).