Gauss Jordan method to solve linear systems.  Works on rectangular matrices, getting the pivots and making zeros below and above.