Exact solutions to a class of time fractional evolution systems with variable coefficients

We explicitly give new group invariant solutions to a class of Riemann-Liouville time fractional evolution systems with variable coefficients. These solutions are derived from every element in an optimal system of Lie algebras generated by infinitesimal symmetries of evolution systems in the class. We express the solutions in terms of Mittag-Leffler functions, generalized Wright functions, and Fox H-functions and show that these solutions solve diffusion-wave equations with variable coefficients. These solutions contain previously known solutions as particular cases. Some plots of solutions subject to the order of the fractional derivative are illustrated.