Numerical method for fractional diffusion-wave equations with functional delay

For a fractional diffusion-wave equation with a nonlinear effect of functional delay, an implicit numerical method is constructed. The scheme is based on the L2-method of approximation of the fractional derivative of the order from 1 to 2, interpolation and extrapolation with the given properties of discrete prehistory and an analogue of the Crank-Nicolson method. The order of convergence of the method is investigated using the ideas of the general theory of difference schemes with heredity. The order of convergence of the method is more significant than in previously known methods, depending on the order of the starting values. The main point of the proof is the use of the stability of the L2-method. The results of comparing numerical experiments with other schemes are presented: a purely implicit method and a purely explicit method, these results showed, in general, the advantages of the proposed scheme.