“…In the literature the big challenge with such equations is the design of efficient and accurate numerical methods and computational cost is the main issue to be considered for any numerical scheme. For classical integer order ordinary or partial differential equations (ODEs/PDEs) such as: systems of ODEs, Navier-Stokes equations, mixed Stokes-Darcy model, shallow water equations, advection-diffusions problems, convection-diffusionreaction equations, heat conduction [22,24,52,40,6,28,30,31,34,14,42,35,37,10,27], a wide class of numerical approaches have been deeply analyzed: finite difference techniques, two-level MacCormack procedure, spectral methods, full implicit finite difference schemes, two-level factored approaches, compact ADI methods and multi-level finite difference formulations. For more details, we refer the readers to [7,11,33,32,36,49,15,38,13,38,23,44,26,25,45,39,46,29] and references therein.…”