ABSTRACT
Based on the Lubich's high-order operators, a second-order temporal finite-difference method is considered for the fractional sub-diffusion equation. It has been proved that the finite-difference scheme is unconditionally stable and convergent in norm by the energy method in both one- and two-dimensional cases. The rate of convergence is order of two in temporal direction under the initial value satisfying some suitable conditions. Some numerical examples are given to confirm the theoretical results.
Acknowledgments
The authors are grateful to anonymous referees for their valuable comments and suggestions to improve this work.
Disclosure statement
No potential conflict of interest was reported by the authors.