In this chapter we discuss the problem of the fast and reliable computation of the Mittag-Leffler function. We start from the series representation, by which the function is commonly defined, to illustrate how it turns out to be not suitable for computation in most of the cases; then we present an alternative approach based on the numerical inversion of the Laplace transform. In particular, we describe a technique, known as the optimal parabolic contour, in which the inversion of the Laplace transform is performed by applying the trapezoidal rule along a parabolic contour in the complex plane; the contour and the integration parameters are chosen, on the basis of the error analysis and the location of the singularities of the Laplace transform, with the aim of achieving a target accuracy (which can be very close to machine precision) with a substantially low computational effort. Applications to the evaluation of derivatives of the Mittag-Leffler function and to matrix arguments are also discussed.

Fast methods for the computation of the Mittag-Leffler function / Garrappa, Roberto; Popolizio, Marina (DE GRUYTER REFERENCE). - In: Handbook of Fractional Calculus with Applications. Volume 3 : Numerical Methods / [a cura di] George Em Karniadakis. - STAMPA. - Berlin : De Gruyter, 2019. - ISBN 9783110570830. - pp. 329-346 [10.1515/9783110571684-013]

Fast methods for the computation of the Mittag-Leffler function

Marina Popolizio
2019-01-01

Abstract

In this chapter we discuss the problem of the fast and reliable computation of the Mittag-Leffler function. We start from the series representation, by which the function is commonly defined, to illustrate how it turns out to be not suitable for computation in most of the cases; then we present an alternative approach based on the numerical inversion of the Laplace transform. In particular, we describe a technique, known as the optimal parabolic contour, in which the inversion of the Laplace transform is performed by applying the trapezoidal rule along a parabolic contour in the complex plane; the contour and the integration parameters are chosen, on the basis of the error analysis and the location of the singularities of the Laplace transform, with the aim of achieving a target accuracy (which can be very close to machine precision) with a substantially low computational effort. Applications to the evaluation of derivatives of the Mittag-Leffler function and to matrix arguments are also discussed.
2019
Handbook of Fractional Calculus with Applications. Volume 3 : Numerical Methods
9783110570830
De Gruyter
Fast methods for the computation of the Mittag-Leffler function / Garrappa, Roberto; Popolizio, Marina (DE GRUYTER REFERENCE). - In: Handbook of Fractional Calculus with Applications. Volume 3 : Numerical Methods / [a cura di] George Em Karniadakis. - STAMPA. - Berlin : De Gruyter, 2019. - ISBN 9783110570830. - pp. 329-346 [10.1515/9783110571684-013]
File in questo prodotto:
Non ci sono file associati a questo prodotto.

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11589/223438
Citazioni
  • Scopus 2
  • ???jsp.display-item.citation.isi??? ND
social impact