The gausselimination command performs Gaussian elimination (row reduction) on the matrix A. An optional second argument P specifies the modulus.
gausselimination([[1,2],[3,4]]); [[1, 2], [0, 1]]