“…However, since the kernel of these differential equations is fractional, it is extremely difficult to obtain exact solutions. Therefore, extensive research has been performed on the development of numerical methods for fractional differential equations such as the Chebyshev collocation method [9,10], the Laplace transform method [11,12], DTM [13,14], ADM [15,16], the operational matrices method [17][18][19][20], and the wavelets method [21][22][23].…”