Sibuya probability distributions and numerical evaluation of fractional-order operators

Abstract

In this work we explore the Sibuya discrete probability distribution, which serves as the basis and the main instrument for numerical simulations of Grunwald--Letnikov fractional derivatives by the Monte Carlo method. We provide three methods for simulating the Sibuya distribution. We also introduce the Sibuya-like sieved probability distributions, and apply them to numerical fractional-order differentiation. Additionally, we use the Monte Carlo method for evaluating fractional-order integrals, and suggest the notion of the continuous Sibuya probability distribution. The developed methods and tools are illustrated by examples of computation. We provide the MATLAB toolboxes for simulation of the Sibuya probability distribution, and for the numerical examples.

arXiv comment28 pages, 13 figures

References in corpus (6)

• Numerical evaluation of two and three parameter Mittag-Leffler functions
• Space-time duality for fractional diffusion
• Higher order Grünwald approximations of fractional derivatives and fractional powers of operators
• A compounded random walk for space-fractional diffusion on finite domains
• On self-similar Bernstein functions and corresponding generalized fractional derivatives
• Tossing half-coins and other partial coins: signed probabilities and Sibuya distribution