In this paper, combining the Smith method [30] with inner-outer (IO) iteration [4], a multi-step Smith-IO (MSIO) iteration algorithm for solving the Sylvester matrix equation and coupled Lyapunov matrix equations (CLMEs) is proposed. Moreover, some convergence conditions and comparison results are also given for the MSIO iteration algorithm. Finally, several numerical examples are used to verify the effectiveness of the proposed nu...