Initial problems for nonlinear hyperbolic functional differential systems are considered. Classical solutions are approximated by solutions of suitable quasilinear systems of difference functional equations. The numerical methods used are difference schemes which are implicit with respect to the time variable. Theorems on convergence of difference schemes and error estimates of approximate solutions are presented. The proof of the stability is based on a comparison technique with nonlinear estimates of the Perron type. Numerical examples are given.