A numerical scheme for the solution of neutral integro-differential equations including variable delay

In this study, an effective numerical technique has been introduced for finding the solutions of the first-order integro-differential equations including neutral terms with variable delays. The problem has been defined by using the neutral integro-differential equations with initial value. Then, an alternative numerical method has been introduced for solving these type of problems. The method is expressed by fundamental matrices, Laguerre polynomials with their matrix forms. Besides, the solution has been obtained by using the collocation points with regard to the reduced system of algebraic equations and Laguerre series.