Algorithm 486: Numerical inversion of Laplace transform [D5]

  • Authors:
  • Francoise Veillon

  • Affiliations:
  • Univ. de Grenoble, Grenoble-Gare, France

  • Venue:
  • Communications of the ACM
  • Year:
  • 1974

Quantified Score

Hi-index 48.22

Visualization

Abstract

This work forms part of a thesis presented in Grenoble in March 1972. Improvements made to the Dubner and Abate algorithm for numerical inversion of the Laplace transform [1] have led to results which compare favorably with theirs and those of Bellmann [2], and Stehfest [3]. The Dubner method leads to the approximation formula: ƒ(t) = 2eat/T[1/2Re{F(a)} + ∑∞k-1 Re{F(a + ik&pgr;/T)}cos(k&pgr;t/T)], (1) where F(s) is the Laplace transform of ƒ(t) and a is positive and greater than the real parts of the singularities of ƒ(t).