Multipole expansion of the fundamental solution of a fractional degree of the Laplace operator

A multipole expansion of the fundamental solution of the fractional degree of the Laplace operator is constructed in terms of the Gegenbauer polynomials. Based on the decomposition constructed and the idea of the fast multipole method, we propose a numerical algorithm for solving the fractional differential generalization of the Poisson equation in the two-dimensional and three-dimensional spaces.