32-digit values of the first 100 recurrence coefficients for the half-range bimodal weight function with parameter ε=.02

Listed in Datasets

By Walter Gautschi

Purdue University

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

Version 1.0 - published on 13 Jan 2017 doi:10.4231/R7J38QH3 - cite this Archived on 14 Feb 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^2)*exp[-(x^2-1)^2/(4*ε)] on [0,Inf], ε=.02, are computed by a multicomponent discretization procedure using the routine sr_OPhrbimod(32,100), with dig=34, epsi=.02 entered at the prompt. The value dig=34 has been determined by the routine dig_sOPbimod(100,32,2,32), attesting to the high stability of the procedure. (Both routines may take several hours to run.) The auxiliary routine xhrbimod.m uses Newton's method to compute a zero x*=sqrt(1+u*) of the equation w(x)=1/2, given a close estimate of u*. There are two such zeros, x_1* and x_2*, x_1* < x_2*, and corresponding estimates -.2 and .3 for u*. The software provided in this dataset allows generating an arbitrary number N of recurrence coefficients for arbitrary ε > 0 as well as for different precisions.

Cite this work

Researchers should cite this work as follows:

Tags

Notes

The dataset consists of one text file and eight 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).