Mixed spline function method for reaction-subdiffusion equations

In this paper, we develop two classes of finite difference schemes for the reaction-subdiffusion equations by using a mixed spline function in space direction, forward and backward difference in time direction, respectively. It has been shown that some of the previous known difference schemes can be derived from our schemes if we suitably choose the spline parameters. By Fourier method, we prove that one class of difference scheme is unconditionally stable and convergent, the other is conditionally stable and convergent. Finally, some numerical results are provided to demonstrate the effectiveness of the proposed difference schemes. (C) 2013 Elsevier Inc. All rights reserved.