Alternative numerical computation of one-sided Levy and Mittag-Leffler distributions

Abstract

We consider here the recently proposed closed form formula in terms of the Meijer G-functions for the probability density functions $g_{\alpha }(x)$ of one-sided Lévy stable distributions with rational index $\alpha =l/k$, with $0<\alpha <1$. Since one-sided Lévy and Mittag-Leffler distributions are known to be related, this formula could also be useful for calculating the probability density functions ${\rho }_{\alpha }(x)$ of the latter. We show, however, that the formula is computationally inviable for fractions with large denominators, being unpractical even for some modest values of $l$ and $k$. We present a fast and accurate numerical scheme, based on an early integral representation due to Mikusinski, for the evaluation of $g_{\alpha }(x)$ and ${\rho }_{\alpha }(x)$, their cumulative distribution function and their derivatives for any real index $\alpha \in (0,1)$. As an application, we explore some properties of these probability density functions. In particular, we determine the location and value of their maxima as functions of the index $\alpha$. We show that $\alpha \approx 0.567$ and $\alpha \approx 0.605$ correspond, respectively, to the one-sided Lévy and Mittag-Leffler distributions with shortest maxima. We close by discussing how our results can elucidate some recently described dynamical behavior of intermittent systems.

6 pages, 5 figures. New references added, final version to appear in PRE. Numerical code available at http://vigo.ime.unicamp.br/distr

References in corpus (4)

• Probability distributions generated by fractional diffusion equations
• Exact and explicit probability densities for one-sided Levy stable distributions
• Pesin-Type Identity for Weak Chaos
• Separation of trajectories and its Relation to Entropy for Intermittent Systems with a Zero Lyapunov exponent

Cited by in corpus (12)

• Levy stable distributions via associated integral transform
• Powerlaw spectra from stochastic acceleration
• Operator solutions for fractional Fokker-Planck equations
• A Tutorial on the Basic Special Functions of Fractional Calculus
• Thermodynamic phase transitions for Pomeau-Manneville maps
• Pesin-type relation for subexponential instability
• Anomalous diffusion, non-Gaussianity, and nonergodicity for subordinated fractional Brownian motion with a drift
• Lyapunov statistics and mixing rates for intermittent systems
• Number of first-passage times as a measurement of information for weakly chaotic systems
• Exact invariant measures: How the strength of measure settles the intensity of chaos
• Learning-augmented count-min sketches via Bayesian nonparametrics
• From Weakly Chaotic Dynamics to Deterministic Subdiffusion via Copula Modeling