“…The author [21] developed a twostep numerical scheme with convergence order O(τ 2− γ 2 + h 4 ) for solving the time-fractional convectiondiffusion-reaction equation with constant order derivative. For classical integer order ordinary/partial differential equations such as: Navier-Stokes equations, systems of ODEs, mixed Stokes-Darcy model, Shallow water problem, convection-diffusion-reaction equation, advection-diffusion model and conduction equation [18,5,20,24,39,26,29,37,28,27,14,32,23,50], a wide set of numerical techniques have been developed and deeply analyzed. For more details, we refer the readers to [30,38,19,45,33,34,44,22] and references therein.…”