In the present paper, we mainly investigate the qualitative behavior of the solutions of a discrete system of difference equations x n+1 = α + m i=1 x n−i y n , y n+1 = β + m i=1 y n−i x n , n ∈ N where α, β ∈ (0, ∞), m ∈ Z + , x −i and y −i are non-negative real numbers for i ∈ {0, 1,. .. , m}. Namely, we discuss the boundedness character and the asymptotic stability properties of steady states of the mentioned system. Finally, for this system, we give a rate of convergence result which has an important place in the discrete dynamical systems. Besides, some numerical simulations with graphs are given to emphasize the efficiency of our theoretical results in the article.