The authors propose a hybrid method for reducing the order of the high order dynamic systems. In the proposed method, the denominator polynomial of the reduced order model is obtained by using the basic characteristics of the original system, which are maintained in the reduced model while the coefficients of the numerator are obtained, by using Routh array method. This method is fundamentally simple and generates stable reduced models if the original high- order system is stable. The proposed method is illustrated with the help of the numerical example taken from the literature.