“…Many real-life phenomena have been described using fractional differential equations (FDEs), such as viscoelasticy (Koeller 1984), continuum mechanics (Carpinteri and Mainardi 1997), optimal control (Bhrawy and Ezz-Eldien 2016), hydrologic modeling (Benson et al 2013), variational problems (Ezz-Eldien 2016), fluid mechanics (Kulish and Lage 2002), finance (Jiang et al 2012), and others (Gaul et al 1991;Ferreira et al 2008; Ezz-Eldien and El-Kalaawy 2018; Dzielinski et al 2010). Searching for numerical techniques for approximating the solutions of FDEs has been strongly considered in the last decades, such as the radial basis functions method (Hosseini et al 2014), the Haar wavelet method , the fractional finite volume method (Liu et al 2014), the Adomian decomposition method (Khodabakhshi et al 2014), the operational method Kim et al 2014), and so on (Heydari et al 2013;Bhrawy and Zaky 2017;Ezz-Eldien and Doha 2018;Stern et al 2014). As a spectral approach for solving numerically some kinds of differential equations, the tau method is classified as one of the most important used methods due to its high convergence and accuracy.…”