Abstract
This paper presents a numerical method for the approximate solution of ordinary differential equations with retarded argument. This method is essentially based on the natural spline functions of even degree introduced by using the derivative-interpolating conditions on simple knots. The estimation of the error, as well as, the convergence of the proposed method are investigated. Some numerical examples are given.