systems. In this method, the denominator polynomial of the reduced order model is obtained
by using the modified pole clustering technique while the coefficients of the numerator are
obtained by Pade approximations. This method is conceptually simple and always
generates stable reduced models if the original high-order system is stable. The proposed
method is illustrated with the help of the numerical examples taken from the literature.