A numerical method for solving Caputo's and Riemann-Liouville's fractional differential equations which includes multi-order fractional derivatives and variable coefficients. - David E. Betancur-Herrera; Nicolás Muñoz-Galeano (Commun. Nonlinear Sci. Numer. Simul., 2020)