In a recently proposed method of model reduction for discrete interval systems, the denominator polynomial of a reduced model is computed by applying interval arithmetic to dominant poles of the original system. However, the denominator polynomial obtained via interval arithmetic usually has poles with larger intervals than desired ones. Hence an unstable polynomial can be derived from the stable polynomial. In this paper a simple technique is presented to partially overcome such a stability problem by accurately preserving desired real dominant poles.