Resumen
A physical-mathematical approach to anomalous diffusion is based on a generalized diffusion equation containing derivatives of fractional order. In this paper, an anomalous subdiffusion equation (ASub-DE) is considered. A new implicit numerical method (INM) and two solution techniques for improving the order of convergence of the INM for solving the ASub-DE are proposed. The stability and convergence of the INM are investigated by the energy method. Some numerical examples are given. The numerical results demonstrate the effectiveness of theoretical analysis. These methods and supporting theoretical results can also be applied to other fractional integro-differential equations and higher-dimensional problems.