Gaussian elimination without row exchanges gives the LU decomposition
Because has unit diagonal and is upper triangular,
Solving the two triangular systems gives
and