“…Many researchers in different branches of physics, mathematics, and engineering are interested in studying fractional calculus due to its many real-life applications, such as ultrasonic wave propagation, 1 optimal control, 2,3 viscoelastic behavior modeling, 4 linear viscoelasticity, 5 variational problems, 6 fluid mechanics, 7 and other applications have been presented in previous studies. [8][9][10][11][12][13] Searching for numerical techniques to solve fractional differential equations has been strongly considered over the last few decades, some of these methods are the fractional finite volume method, 14 spectral collocation method, 15 random variable transformation technique, 16 operational matrix method, 17 Adomian decomposition method, 18 power series method, 19 nonpolynomial spline method, 20 Gauss-collocation method, 21 radial basis functions method, 22 modified Galerkin algorithm, 23 and other methods have been introduced; see previous studies. [24][25][26][27][28][29][30][31] The tau approach is considered as one of the highly accurate methods that have been used for different kinds of fractional differential equations.…”