Gaussian elimination
Row echelon form
+
-