1 Gaussian elimination

Recall from HELM booklet  8 that the basic idea with Gaussian (or Gauss) elimination is to replace the matrix of coefficients with a matrix that is easier to deal with. Usually the nicer matrix is of upper triangular form which allows us to find the solution by back substitution . For example, suppose we have

x 1 + 3 x 2 5 x 3 = 2 3 x 1 + 11 x 2 9 x 3 = 4 x 1 + x 2 + 6 x 3 = 5

which we can abbreviate using an augmented matrix to

1 3 5 2 3 11 9 4 1 1 6 5 .

We use the boxed element to eliminate any non-zeros below it. This involves the following row operations

1 3 5 2 3 11 9 4 1 1 6 5 R 2 3 × R 1 R 3 + R 1 1 3 5 2 0 2 6 2 0 4 1 7 .

And the next step is to use the 2 to eliminate the non-zero below it . This requires the final row operation

1 3 5 2 0 2 6 2 0 4 1 7 R 3 2 × R 2 1 3 5 2 0 2 6 2 0 0 11 11 .

This is the augmented form for an upper triangular system, writing the system in extended form we have

x 1 + 3 x 2 5 x 3 = 2 2 x 2 + 6 x 3 = 2 11 x 3 = 11

which is easy to solve from the bottom up, by back substitution .

Example 5

Solve the system

x 1 + 3 x 2 5 x 3 = 2 2 x 2 + 6 x 3 = 2 11 x 3 = 11
Solution

The bottom equation implies that x 3 = 1 . The middle equation then gives us that

2 x 2 = 2 6 x 3 = 2 + 6 = 4 x 2 = 2

and finally, from the top equation,

x 1 = 2 3 x 2 + 5 x 3 = 2 6 5 = 9.

Therefore the solution to the problem stated at the beginning of this Section is

x 1 x 2 x 3 = 9 2 1 .

The following Task will act as useful revision of the Gaussian elimination procedure.

Task!

Carry out row operations to reduce the matrix

2 1 4 4 3 1 6 8 2

into upper triangular form.

The row operations required to eliminate the non-zeros below the diagonal in the first column are as follows

2 1 4 4 3 1 6 8 2 R 2 2 × R 1 R 3 + 3 × R 1 2 1 4 0 5 9 0 5 10

Next we use the 5 on the diagonal to eliminate the 5 below it:

2 1 4 0 5 9 0 5 10 R 3 R 2 2 1 4 0 5 9 0 0 19

which is in the required upper triangular form.