in this paper, an iterative method is proposed for solving matrix equation $sum_{j=1}^s a_jx_jb_j = e$. this method is based on the global least squares (gl-lsqr) method for solving the linear system of equations with the multiple right hand sides. for applying the gl-lsqr algorithm to solve the above matrix equation, a new linear operator, its adjoint and a new inner product are dened. it is ...