Rational approximation to the fractional Laplacian operator in reaction-diffusion problems

This paper provides a new numerical strategy to solve fractional in space reaction-diffusion equations on bounded domains under homogeneous Dirichlet boundary conditions. Using the matrix transform method the fractional Laplacian operator is replaced by a matrix which, in general, is dense. The approach here presented is based on the approximation of this matrix by the product of two suitable banded matrices. This leads to a semi-linear initial value problem in which the matrices involved are sparse. Numerical results are presented to verify the effectiveness of the proposed solution strategy.