Before revision:   A number of numerical results are presented to support the error analysis for both smooth and nonsmooth solutions .
After revision:  Error estimates are given for both smooth solutions and solutions with a singularity at t = 0 of a type that is typical for equations involving fractional time-derivatives. A number of numerical results are presented to support the error analysis for both smooth and nonsmooth solutions .