Construction and implementation of asymptotic expansions for Laguerre-type orthogonal polynomials

Abstract

Laguerre and Laguerre-type polynomials are orthogonal polynomials on the interval [0,∞) with respect to a weight function of the form w(x) = xα e-Q(x), Q(x) = Σk=0m qk xk, α> -1, qm > 0. The classical Laguerre polynomials correspond to Q(x)=x. The computation of higher-order terms of the asymptotic expansions of these polynomials for large degree becomes quite complicated, and a full description seems to be lacking in literature. However, this information is implicitly available in the work of Vanlessen, based on a non-linear steepest descent analysis of an associated so-called Riemann--Hilbert problem. We will extend this work and show how to efficiently compute an arbitrary number of higher-order terms in the asymptotic expansions of Laguerre and Laguerre-type polynomials. This effort is similar to the case of Jacobi and Jacobi-type polynomials in a previous paper. We supply an implementation with explicit expansions in four different regions of the complex plane. These expansions can also be extended to Hermite-type weights of the form (-Σk=0m qk x2k) on (-∞,∞), and to general non-polynomial functions Q(x) using contour integrals. The expansions may be used, e.g., to compute Gauss-Laguerre quadrature rules in a lower computational complexity than based on the recurrence relation, and with improved accuracy for large degree. They are also of interest in random matrix theory.

0

Turn this paper into a lesson

ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…