En el presente trabajo, consideramos la ecuación de matriz de Stein a gran escala con un término constante de rango bajo. Estas ecuaciones de matriz aparecen en muchas aplicaciones en problemas de control en tiempo discreto, filtrado y restauración de imágenes, entre otros. Los métodos propuestos se basan en la proyección sobre el subespacio de Krylov en bloque extendido con un enfoque de Galerkin (GA) o con la minimización de la norma del residual. Damos algunos resultados sobre las normas del residual y del error y reportamos algunos experimentos numéricos.