A. The purpose of this paper is to study stability of numerical solutions for a coupled differential equation with piecewise constant arguments with one delay. A sufficient condition such that the equation is asymptotically stable is derived. Moreover, when the linear θ-method is applied to this equation, it is proved that the linear θ-method is asymptotically stable if and only if 1/2 < θ ≤ 1. Some numerical examples to confirm the theoretical results are given.2010 Mathematics Subject Classification. 65L07, 65L20.Key words and phrases. piecewise constant arguments; coupled differential equation; linear θ-method; numerical solution; Stability.where X(t) = (x(t), y(t), z(t)) T , X 0 = (x 0 , y 0 , z 0 ) T and A = a 0 0 0 c 0 0 0 e , B = 0 b 0 0 0 d f 0 0 .