“…The Caputo fractional derivatives are used in this context to define the time fractional derivative. FPDE can be solved using a variety of numerical techniques, including the finite element method [1], VIM [2,3], finite difference method [4,5], homotopy perturbation method [6], Tau method [7,8], operational method [9], Ritz approximation approach [10], and Caputo-Hadamard fractional derivative [11]. Recently, several numerical schemes have been created to solve the TFCDE such as Sinc-Legendre collocation method [12], discontinuous Galerkin method [13], Gegenbauer method [14], Sinc and B-spline scaling functions [15], Haar wavelet method [16], Chebyshev collocation method [17], and radial basis functions [18].…”