Abstract:The Laguerre method for numerically inverting Laplace transforms is an old established method based on the 1935 Tricomi–Widder theorem, which shows (under suitable regularity conditions) that the desired function can be represented as a weighted sum of Laguerre functions, where the weights are coefficients of a generating function constructed from the Laplace transform using a bilinear transformation. We present a new variant of the Laguerre method based on: (i) using our previously developed variant of the Fourier-series method to calculate the coefficients of the Laguerre generating function, (ii) developing systematic methods for scaling, and (iii) using Wynn's ϵ-algorithm to accelerate convergence of the Laguerre series when the Laguerre coefficients do not converge to zero geometrically fast. These contributions significantly expand the class of transforms that can be effectively inverted by the Laguerre method. We provide insight into the slow convergence of the Laguerre coefficients as well as propose a remedy. Before acceleration, the rate of convergence can often be determined from the Laplace transform by applying Darboux's theorem. Even when the Laguerre coefficients converge to zero geometrically fast, it can be difficult to calculate the desired functions for large arguments because of roundoff errors. We solve this problem by calculating very small Laguerre coefficients with low relative error through appropriate scaling. We also develop another acceleration technique for the case in which the Laguerre coefficients converge to zero geometrically fast. We illustrate the effectiveness of our algorithm through numerical examples.

Improved Laguerre Spectral Methods with Less Round-off Errors and Better Stability

Efficient Legendre-Laguerre spectral element methods for problems on unbounded domains with diagonalization technique

On the Laguerre Method for Numerically Inverting Laplace Transforms

Convergence analysis of Laguerre approximations for analytic functions

Optimal pointwise error estimates for piecewise functions expanded with Laguerre polynomials

Approximation by a new sequence of operators involving Laguerre polynomials

Binomial Theorem Involving Hermite Polynomials and Negative-Binomial Theorem Involving Laguerre Polynomials

A note on Laguerre truncated polynomials and quadrature formula

Laguerre-Type Bernoulli and Euler Numbers and Related Fractional Polynomials

Bernstein-Polynomials-Based Highly Accurate Methods for One-Dimensional Interface Problems

Scaling Optimized Hermite Approximation Methods

Two Efficient Generalized Laguerre Spectral Algorithms for Fractional Initial Value Problems

On a power series involving classical orthogonal polynomials

On the application of Laguerre's method to the polynomial eigenvalue problem

Asymptotics for Laguerre Polynomials with Large Order and Parameters.

An Efficient Solution of Multiplicative Differential Equations through Laguerre Polynomials

Efficient and Accurate Spectral Method for Solving Fractional Differential Equations on the Half Line Using Orthogonal Generalized Rational Jacobi Functions

New error bounds for Legendre approximations of differentiable functions

Series solutions of Laguerre- and Jacobi-type differential equations in terms of orthogonal polynomials and physical applications

The numerical solution of fractional integral equations via orthogonal polynomials in fractional powers