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.

DEALING WITH ABEL INVERSION BY USING LENGENDRE EXPANSIONS

Dealing with Asymmetrical Abel Inversion Using Fourier-Bessel Expansions

A New Method for Numerical Abel-Inversion

Abel inversion using Legendre polynomials approximations

Modified Fourier-Hankel Method Based On Analysis Of Errors In Abel Inversion Using Fourier Transform Techniques

Abel inversion using Legendre wavelets expansion

An Inexact Cayley Transform Method for Inverse Eigenvalue Problems with Multiple Eigenvalues

An Ulm-like Cayley Transform Method for Inverse Eigenvalue Problems with Multiple Eigenvalues

A Generalization of Abel Inversion to non axisymmetric density distribution

Inversion of the Abel--Prym map in presence of an additional involution

A Ulm-like Method for Inverse Eigenvalue Problems

Extended Newton-type Method for Inverse Singular Value Problems with Multiple And/or Zero Singular Values

Abel inversion using fast Fourier transforms.

General Solutions with Hidden Variables to the Abel Differential Equation

Some New Results on Abel Equations

Adaptive Inversion Method for the Laplace Integral Transform

Two methods addressing variable-exponent fractional initial and boundary value problems and Abel integral equation

On the Laguerre Method for Numerically Inverting Laplace Transforms

A Divided Difference Approach to Abel and Hagen-Rothe Identities.

Long-time Asymptotics for the Defocusing Ablowitz-Ladik System with Initial Data in Lower Regularity

The Abel, Fourier and Radon transforms on symmetric spaces