  We develop a fully discrete scheme for time-fractional diffusion equations by
using a finite difference method in time and a finite element method in space.
The fractional derivatives are used in Caputo sense. Stability and error
estimates are derived. The accuracy and efficiency of the presented method is
shown by conducting two numerical examples.
