We are concerned with stability of numerical methods for delay differential systems of neutral type. In particular, delay-dependent stability of numerical methods is investigated. By means of the H -matrix norm, a necessary and sufficient condition for the asymptotic stability of analytic solution of linear neutral differential systems is derived. Then, based on the argument principle, sufficient conditions for delay-dependent stability of Runge-Kutta and linear multi-step methods are presented, respectively. Furthermore, two algorithms are provided for checking delay-dependent stability of analytical and numerical solutions, respectively. Numerical examples are given to illustrate the main results.