Gauss Reduction, the fast way

First, let's make the matrix "A", which contains the coeficients of the variables in a set of linear equations:

> A:=matrix(3,3,[0,-1,2,2,-3,1,-4,5,2]);

Now, let's make the vector "b", which contains the constant terms in the set of linear equations:

> b:=vector([-7,7,3]);

Now, we make the augmented matrix

> Ab:=augment(A,b);

The linalg package has a command called "gausselim" which does Gauss Reduction:

> gausselim(Ab);

See how easy that was!

>