Este documento describe un método iterativo de subespacio de Krylov diseñado para resolver sistemas de ecuaciones lineales con una matriz grande, simétrica, no singular e indefinida. Este método está diseñado para permitir la evaluación de estimaciones de error para las iteraciones calculadas. La disponibilidad de estimaciones de error hace posible terminar el proceso iterativo cuando el error estimado es menor que una tolerancia especificada por el usuario. Las estimaciones de error se calculan aprovechando la relación entre las iteraciones y las reglas de cuadratura de tipo Gauss. Ejemplos computados ilustran el rendimiento del método iterativo y las estimaciones de error.