“…In recent years, fractional equations appear in modeling of various real-life phenomena, for example dynamic viscoelasticity modeling (Larsson et al 2015), hydrologic (Benson et al 2013), economics (Baillie 1996), temperature and motor control (Bohannan 2008), solid mechanics (Rossikhin and Shitikova 1997), bioengineering (Magin 2004), medicine (Hall and Barrick 2008), porous media (He 1998), fluid-dynamic traffic model (He 1999), etc. Therefore, there is an increasing demand for numerical and analytical solutions of various types of fractional differential equations such as finite-difference/finite-element technique Abbaszadeh 2018a, b, 2019), compact difference schemes (Hu and Zhang 2012), homotopy analysis method (Dehghan et al 2010), homotopy perturbation method (Abdulaziz et al 2008), dual reciprocity boundary elements method (Dehghan and Safarpoor 2016), fifth-kind orthonormal Chebyshev polynomial method (Abd-Elhameed and Youssri 2018), B-spline functions (Lakestani et al 2012), fractional-order Bernoulli function method , hybrid method (Mashayekhi and Razzaghi 2015), Legendre wavelet method (Heydari et al 2014), fractional-order Lagrange polynomials (Sabermahani et al 2018), fractional-order Bernoulli wavelets method (Rahimkhani et al 2016), Müntz-Legendre wavelet method (Rahimkhani et al 2018), etc.…”