can be solved with the aid of Givens transformations.
The -th transformation modifies rows and of the system,
annihilating the subdiagonal coefficient
and generating
an additional non-zero coefficient in position
. The resulting
upper triangular system is solved with backward substitution.