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/R78P5XHT - 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=116 has been determined by the routine dig_laglogmm(100,1/2,108,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/R7JH3J5S for a method based on ordinary moments. Both methods produce the same answers, but the current one is about 60% faster 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).