In this paper we are concerned with the numerical solution of a diffusion equation in which the time derivative is of non-integer order, i.e., in the interval (0, 1). An implicit numerical method is presented and its unconditional stability and convergence are proved. Two numerical examples are provided to illustrate the obtained theoretical results

CEMAT - Center for Computational and Stochastic Mathematics