An improved Talbot method for numerical Laplace transform inversion

Autor:	Dingfelder, Benedict, Weideman, J. A. C.
Rok vydání:	2013
Předmět:	Mathematics - Numerical Analysis 65R10 (numerical Inversion of Laplace Transforms) 65B99 (increased speed of convergence of Talbot's method) 65Y15 (attached piece of code)
Druh dokumentu:	Working Paper
Popis:	The classical Talbot method for the computation of the inverse Laplace transform is improved for the case where the transform is analytic in the complex plane except for the negative real axis. First, by using a truncated Talbot contour rather than the classical contour that goes to infinity in the left half-plane, faster convergence is achieved. Second, a control mechanism for improving numerical stability is introduced. These two features are incorporated into a software code, whose performance is assessed on transforms from tables as well as from actual applications. It is shown that even when the transform has singularities off the negative real axis, rapid convergence can still be achieved in many cases. Comment: 16 pages, 7 figures, 1 table, 1 MATLAB function
Databáze:	arXiv
Externí odkaz:	http://arxiv.org/abs/1304.2505 Zobrazit plný text záznamu View this record from Arxiv