In this paper, we consider the numerical resolution of a time and space fractional diffusion equation. The main purpose of this work is to construct an efficient accurate numerical solution by using spline function and then we analyze the stability of the obtained scheme for the time-space fractional diffusion equation. Numerical experiments are carried out to support the theoretical claims.