Convergence of numerical methods for systems of neutral functional-differential-algebraic equations