By Gaussian elimination any LSE with
invertible coefficient matrix
can be brought to upper triangular form in
at most steps.
For this purpose all coefficients below the diagonal
are successively nullified, that is, after
steps the LSE has the form
In detail the -th elimination step
proceeds as follows:
If the -th diagonal element
(the so called pivot element) equals zero,
then the -th equation is interchanged with
one of the following equations so that
.
For an appropriate multiple of the
-th equation is subtracted from
the -th equation so that
equals zero: