Abstract
This paper proposes a combined operational matrix approach based on Lucas and Taylor polynomials for the solution of neutral type differential equations with proportional delays. The advantage of the proposed method is the ease of its application. The method facilitates the solution of the given problem by reducing it to a matrix equation. Illustrative examples are validated by means of absolute errors. Residual error estimation is presented to improve the solutions. Presented in graphs and tables the results are compared with the existing methods in literature.