An $\epsilon$-approximation solution of time-fractional diffusion equations based on Legendre polynomials

The purpose of this paper is to establish a numerical method for solving time-fractional diffusion equations. To obtain the numerical solution, a binary reproducing kernel space is defined, and the orthonormal basis is constructed based on Legendre polynomials in this space. In order to find the $\epsilon$-approximation solution of time-fractional diffusion equations, which is defined in this paper, the algorithm is designed using the constructed orthonormal basis. Some numerical examples are analyzed to illustrate the procedure and confirm the performance of the proposed method. The results faithfully reveal that the presented method is considerably accurate and effective, as expected.