Model order reduction is the approximation of dynamical systems into equivalent systems with smaller order. Model reduction has been studied extensively for different types of systems. In this paper, we present two methods for multi input multi output linear systems. These methods are based on solvents, also called block poles. These methods are particularly suitable if the given system is in matrix transfer function form. The first method eliminates solvents one by one whereas, the second method can eliminate multiple solvents at the same time. The two presented methods are implemented in MATLAB in order to provide a systematic method for the model order reduction of MIMO linear systems.