This paper is concerned with numerical stability of nonlinear systems of pantograph equations. Numerical methods based on (k, l)-algebraically stable Runge-Kutta methods are suggested. Global and asymptotic stability conditions for the presented methods are derived.