Gauss Reduction, the long way

Here's the augmented matrix which we made previously:

> evalm(Ab);

Now, let's do Gauss Reduction using elementary row operations

First, we swap rows 1 and 2:

> B:=swaprow(Ab,1,2);

Next, we add row 1, multiplied by 2, to row 3, and replace row 3 with the result:

> B:=addrow(B,1,3,2);

Then, add row 2, multiplied by -1 to row 3, and replace row 3 with the result:

> B:=addrow(B,2,3,-1);

This is the Row-Echelon form of the augmented matrix. Let's see if it agrees with what "gausselim" gives us:

> evalm(B-gausselim(Ab));

We have agreement!

However, we were just lucky. The method of Gauss Reduction is non-unique. Another choice of row operations would lead to a different form of the augmented matrix. See the next section.

>