Routh approximation method applied to order reduction of linear MIMO systems