This paper gives a new numerical method for solving a Sylvester equation. Our method is based on reducing only one coefficent matrix to a triangular form. On depending on the dimension of the matrices, our method is more effective than the Bartels-Stewart algorithm.
