“…Many other numerical and approximation methods as well as computational approaches have been developed and applied for the FDEs which are based upon various closely-related models of real-world problems. For example, Baleanu et al [ 23 ] made use of a Chebyshev spectral method based on operational matrices, a remarkable survey of numerical methods can be found in [ 24 ], a study of the fractional-order Bessel, Chelyshkov, and Legendre collocation schemes for the fractional Riccati equation was presented in [ 25 ], an operational matrix of fractional-order derivatives of Fibonacci polynomials was developed in [ 26 ], an introductory overview and recent developments involving FDEs was presented in [ 27 ], efficiency of the spectral collocation method in the dynamic simulation of the fractional-order epidemiological model of the Ebola virus was investigated in [ 28 ], the Jacobi collocation method and a spectral tau method based on shifted second-kind Chebyshev polynomilas for the approximate solution of some families of the fractional-order Riccati differential equations were discussed in [ 29 , 30 ], computational approaches to FDEs for the biological population model were discussed in [ 31 ], the generalized Chebyshev and Bessel colllocation approaches for fractional BVPs and multi-order FDEs were considered in [ 32 , 33 ], and a general wavelet quasi-linearization method for solving fractional-order population growth model was developed and applied in [ 34 ].…”