This paper is concerned with the approximate solution of functional differential equations having the form. We search for a solution x, defined for .. which takes given values on intervals [1, 0] and (k-1, k]. We introduce and analyse some new computational methods for the solution of this problem. Numerical results are presented and compared with the results obtained by other methods.