Changed the "piv" method of the Gaussian Eliminator to keep all columns of the tableaus with constants times the identity matrix to have the same value "d".