# Inverse of a Matrix using Gauss-Jordan Elimination

by M. Bourne

In this section we see how Gauss-Jordan Elimination works using examples.

You can re-load this page as many times as you like and get a new set of numbers each time. You can also choose a different size matrix (at the bottom of the page).

(If you need some background first, go back to the Introduction to Matrices).

Choose the matrix size you are interested in and then click the button.

Matrix A:

The randomly-generated example appears below.

### Phone users

**NOTE:** If you're on a phone, you can scroll any **wide matrices** on this page to the right or left to see the whole expression.

## Example (3 × 3)

Find the inverse of the matrix *A* using Gauss-Jordan elimination.

A = |
12 | 7 | 4 | ||

5 | 13 | 11 | |||

14 | 15 | 16 |

### Our Procedure

We write matrix *A* on the left and the Identity matrix *I* on its right separated with a dotted line, as follows.
The result is called an **augmented** matrix.

We include row numbers to make it clearer.

12 | 7 | 4 | |

5 | 13 | 11 | |

14 | 15 | 16 |

1 | 0 | 0 | Row[1] | |

0 | 1 | 0 | Row[2] | |

0 | 0 | 1 | Row[3] |

Next we do several **row operations** on the 2 matrices and our aim is to end up with the identity matrix on the **left**, like this:

1 | 0 | 0 | |

0 | 1 | 0 | |

0 | 0 | 1 |

? | ? | ? | Row[1] | |

? | ? | ? | Row[2] | |

? | ? | ? | Row[3] |

(Technically, we are reducing matrix *A* to **reduced row echelon form**, also called **row canonical form**).

The resulting matrix on the right will be the **inverse matrix** of *A*.

Our row operations procedure is as follows:

- We get a "1" in the top left corner by dividing the first row
- Then we get "0" in the rest of the first column
- Then we need to get "1" in the second row, second column
- Then we make all the other entries in the second column "0".

We keep going like this until we are left with the identity matrix on the left.

Let's now go ahead and find the inverse.

### Solution

We start with:

12 | 7 | 4 | |

5 | 13 | 11 | |

14 | 15 | 16 |

1 | 0 | 0 | Row[1] | |

0 | 1 | 0 | Row[2] | |

0 | 0 | 1 | Row[3] |

### New Row [1]

**Divide Row [1] by 12** (to give us a "1" in the desired position):

This gives us:

1 | 0.5833 | 0.3333 | |

5 | 13 | 11 | |

14 | 15 | 16 |

0.0833 | 0 | 0 | Row[1] | |

0 | 1 | 0 | Row[2] | |

0 | 0 | 1 | Row[3] |

### New Row [2]

**Row[2] − 5 × Row[1]** (to give us 0 in the desired position):

5 − 5 × 1 = 0

13 − 5 × 0.5833 = 10.083

11 − 5 × 0.3333 = 9.3333

0 − 5 × 0.0833 = -0.4167

1 − 5 × 0 = 1

0 − 5 × 0 = 0

This gives us our new Row [2]:

1 | 0.5833 | 0.3333 | |

0 | 10.083 | 9.3333 | |

14 | 15 | 16 |

0.0833 | 0 | 0 | Row[1] | |

-0.4167 | 1 | 0 | Row[2] | |

0 | 0 | 1 | Row[3] |

### New Row [3]

**Row[3] − 14 × Row[1]** (to give us 0 in the desired position):

14 − 14 × 1 = 0

15 − 14 × 0.5833 = 6.8333

16 − 14 × 0.3333 = 11.333

0 − 14 × 0.0833 = -1.1667

0 − 14 × 0 = 0

1 − 14 × 0 = 1

This gives us our new Row [3]:

1 | 0.5833 | 0.3333 | |

0 | 10.083 | 9.3333 | |

0 | 6.8333 | 11.333 |

0.0833 | 0 | 0 | Row[1] | |

-0.4167 | 1 | 0 | Row[2] | |

-1.1667 | 0 | 1 | Row[3] |

### New Row [2]

**Divide Row [2] by 10.083** (to give us a "1" in the desired position):

This gives us:

1 | 0.5833 | 0.3333 | |

0 | 1 | 0.9256 | |

0 | 6.8333 | 11.333 |

0.0833 | 0 | 0 | Row[1] | |

-0.0413 | 0.0992 | 0 | Row[2] | |

-1.1667 | 0 | 1 | Row[3] |

### New Row [1]

**Row[1] − 0.5833 × Row[2]** (to give us 0 in the desired position):

1 − 0.5833 × 0 = 1

0.5833 − 0.5833 × 1 = 0

0.3333 − 0.5833 × 0.9256 = -0.2066

0.0833 − 0.5833 × -0.0413 = 0.1074

0 − 0.5833 × 0.0992 = -0.0579

0 − 0.5833 × 0 = 0

This gives us our new Row [1]:

1 | 0 | -0.2066 | |

0 | 1 | 0.9256 | |

0 | 6.8333 | 11.333 |

0.1074 | -0.0579 | 0 | Row[1] | |

-0.0413 | 0.0992 | 0 | Row[2] | |

-1.1667 | 0 | 1 | Row[3] |

### New Row [3]

**Row[3] − 6.8333 × Row[2]** (to give us 0 in the desired position):

0 − 6.8333 × 0 = 0

6.8333 − 6.8333 × 1 = 0

11.333 − 6.8333 × 0.9256 = 5.0083

-1.1667 − 6.8333 × -0.0413 = -0.8843

0 − 6.8333 × 0.0992 = -0.6777

1 − 6.8333 × 0 = 1

This gives us our new Row [3]:

1 | 0 | -0.2066 | |

0 | 1 | 0.9256 | |

0 | 0 | 5.0083 |

0.1074 | -0.0579 | 0 | Row[1] | |

-0.0413 | 0.0992 | 0 | Row[2] | |

-0.8843 | -0.6777 | 1 | Row[3] |

### New Row [3]

**Divide Row [3] by 5.0083** (to give us a "1" in the desired position):

This gives us:

1 | 0 | -0.2066 | |

0 | 1 | 0.9256 | |

0 | 0 | 1 |

0.1074 | -0.0579 | 0 | Row[1] | |

-0.0413 | 0.0992 | 0 | Row[2] | |

-0.1766 | -0.1353 | 0.1997 | Row[3] |

### New Row [1]

**Row[1] − -0.2066 × Row[3]** (to give us 0 in the desired position):

1 − -0.2066 × 0 = 1

0 − -0.2066 × 0 = 0

-0.2066 − -0.2066 × 1 = 0

0.1074 − -0.2066 × -0.1766 = 0.071

-0.0579 − -0.2066 × -0.1353 = -0.0858

0 − -0.2066 × 0.1997 = 0.0413

This gives us our new Row [1]:

1 | 0 | 0 | |

0 | 1 | 0.9256 | |

0 | 0 | 1 |

0.071 | -0.0858 | 0.0413 | Row[1] | |

-0.0413 | 0.0992 | 0 | Row[2] | |

-0.1766 | -0.1353 | 0.1997 | Row[3] |

### New Row [2]

**Row[2] − 0.9256 × Row[3]** (to give us 0 in the desired position):

0 − 0.9256 × 0 = 0

1 − 0.9256 × 0 = 1

0.9256 − 0.9256 × 1 = 0

-0.0413 − 0.9256 × -0.1766 = 0.1221

0.0992 − 0.9256 × -0.1353 = 0.2244

0 − 0.9256 × 0.1997 = -0.1848

This gives us our new Row [2]:

1 | 0 | 0 | |

0 | 1 | 0 | |

0 | 0 | 1 |

0.071 | -0.0858 | 0.0413 | Row[1] | |

0.1221 | 0.2244 | -0.1848 | Row[2] | |

-0.1766 | -0.1353 | 0.1997 | Row[3] |

We have achieved our goal of producing the Identity matrix on the left. So we can conclude the inverse of the matrix *A* is the right hand portion of the augmented matrix:

A^{−1} = |
0.071 | -0.0858 | 0.0413 | ||

0.1221 | 0.2244 | -0.1848 | |||

-0.1766 | -0.1353 | 0.1997 |

### Things to Note

- The above explanation shows all steps. A human can usually take a few shortcuts. Also, sometimes there is already a "1" or a "0" in the correct position, and in those cases, we would not need to do anything for that step.
- Always write down what you are doing in each step - it is very easy to get lost!
- I have shown results correct to 4 decimal place, but
**best possible accuracy**was used throughout. Be aware that small errors from rounding will accumulate throughout the problem. Always use full calculator accuracy! (Make full use of your calculator's memory.) - Very occasionally there are strange results because of the computer's internal representation of numbers. That is, it may store "1" as 0.999999999872.

### See another?

You can go back up to the top of the page and choose another example.

### Search IntMath

### Online Algebra Solver

This algebra solver can solve a wide range of math problems.

Go to: Online algebra solver