This survey paper provides a generalization of the classical random walk model of the standard diffusion equation to the so-called Feller (space-)fractional diffusion equation
obtained by replacing the second (space-)derivative with a special pseudo-differential operator
arising by inversion of linear combinations of left- and right-sided Riemann-Liouville-Weyl operators (in case of both coefficients being equal, this leads to inversion of the Riesz potential operators).
In a pioneering paper, W. Feller [Meddel. Lunds Univ. Mat. Sem. Suppl.-Band M. Riesz, 73-81 (1952; Zbl 0048.08503)] considered the problem of generating all the stable probability distributions through the Green function of the Cauchy problem for the linear evolution equation of fractional order (1), in view of semigroup theory. From a probabilistic point of view, these semigroups represent stable distributions in the sense of P. Lévy [“Calcul des probabilités” (1925; JFM 51.0380.02)]. To honour both, Lévy and Feller, the authors call the considered stohastic processes (which are of Markov type), characterized by the two parameters (for the order of differentiation) and (for the “skewness” of the spatial pseudo-differential operator), as “Lévy-Feller processes”. They give a completely revised version of Feller’s work and complement his analysis by an original finite-difference approach to approximate and to find the Green function of the Cauchy problem for (1), excluding the singular case . Their essential idea is to approximate the inverse operators in the expression (2) for the Feller derivative, by the Grünwald-Letnikov scheme, see, for example, R. Gorenflo [in: Fractals and fractional calculus in continuum mechanics. CISM Courses and Lectures 378, 277-290 (1997)]. The authors provide also numerical computations of all the stable, non-Gaussian densities, with exception to those of Cauchy type.