Model Order Reduction by Routh Stability Array with Stability Equation Method for SISO and MIMO Systems