A method for solving a system of linear equations with a symmetric indefinite coefficient matrix
The subject of this article is a numerically stable method for solving nonsingular Cramerian systems of linear equations, with a symmetric indefinite coefficient matrix. It consists in adapting the algorithm presented in [1] that can stably and effectively diagonalize any indefinite symmetric matrix. It is in fact some modification of the Gaussian symmetric elimination procedure.