Abstract
This paper is concerned with a numerical method for solving a model equation of anomalous diffusion expressed with fractional integral with respect to the time variable t. In particular, a semi-discretization scheme derived from the application of the fractional trapezoidal rule proposed by Lubich in the mid 80's to the integral is examined. It is verified that the scheme converges in a special case where the exact solution is represented in a generalized Taylor form in fractional powers of t. Numerical examples are also presented which suggest convergence of the scheme in more general cases.(Theory)
