For fractional diffusion equation with functional after-effect on time, the implicit numerical method is constructed and the order of its convergence is obtained. The method is a fractional analogue of the Crank-Nicholson method, and also uses interpolation and extrapolation of the prehistory of model on time.