The Gauss-Jordan Method, the fast way
Here's the augmented matrix which we made previously:
> evalm(Ab);
The linalg package has a command called "gaussjord" which performs the Gauss-Jordan Method
> gaussjord(Ab);
Pretty easy, huh?
>