Comparison of different numerical methods for the solution of the time-fractional reaction-diffusion equation with variable diffusion coefficient

In this work we perform a comparison of two different numerical schemes for the solution of the time-fractional diffusion equation with variable diffusion coefficient and a nonlinear source term. The two methods are the implicit numerical scheme presented in [M.L. Morgado, M. Rebelo, Numerical approximation of distributed order reaction- diffusion equations, Journal of Computational and Applied Mathematics 275 (2015) 216-227] that is adapted to our type of equation, and a colocation method where Chebyshev polynomials are used to reduce the fractional differential equation to a system of ordinary differential equations

COMPETE, FEDER and Fundação para a Ciência e a Tecnologia (the Portuguese Foundation for Science and Technology (FCT)) through Strategic Projects - LA 25 – 2013-2014, PTDC/EME- MFE/113988/2009 and EXPL/CTM-POL/1299/2013. M. Rebelo acknowledge financial funding by the Portuguese Foundation for Science and Technology through the project PEstOE/MAT/UI0297/2013 (Centro de Matemática e Aplicacões)