Abstract
In this paper a new approach to the use of kernel operators derived from fractional order differential equations is proposed. Three different types of kernels are used, power law, exponential decay and Mittag-Leffler kernels. The kernel's fractional order and fractal dimension are the key parameters for these operators. The main objective of this paper is to study the effect of the fractal-fractional derivative order and the order of the nonlinear term, 1 < q < 2, in the equation on the behavior of numerical solutions of fractal-fractional reaction diffusion equations (FFRDE). Iterative approximations to the solutions of these equations are constructed by applying the theory of fractional calculus with the help of Lagrange polynomial functions. In key parameter regimes, all these iterative solutions based on a power kernel, an exponential kernel and a generalized Mittag-Leffler kernel are very close. Hence, iterative solutions obtained using one of these kernels are compared with full numerical solutions of the FFRDE and excellent agreement is found. All numerical solutions in this paper were obtained using Mathematica.