A Numerical Approach for Solving High-Order Linear Delay Volterra Integro-Differential Equations

In this study, a numerical method is proposed to solve high-order linear Volterra delay integro-differential equations. In this approach, we assume that the exact solution can be expressed as a power series, which we truncate after the (N + 1)-st term so that it become a polynomial of degree N. Substituting the unknown function, its derivatives and the integrals by their matrix counterparts, we obtain a vector equivalent of the equation in question. Applying inner product to this vector with a set of monomials, we are left with a linear algebraic equation system of N unknowns. The approximate solution of the problem is then computed from the solution of the resulting linear system. In addition, the technique of residual correction, whose aim is to increase the accuracy of the approximate solutions by estimating the error of those solutions, is discussed briefly. Both the method and this technique are illustrated with several examples. Finally, comparison of the present scheme with other methods is made wherever possible.