In this paper we consider implicit non-linear neutral delay differential equations to derive efficient numerical schemes with good stability properties. The basic idea is to reformulate the original problem eliminating the dependence on the derivative of the solution in the past values. Our hypothesis on the original equation allow us to study the boundedness and asymptotic stability of the true and numerical solutions by the theory of stability with respect to the forcing term.