**Definition 1**: The **determinant**, det *A, *also denoted* |A|,* of an *n* × *n* square matrix *A* is defined recursively as follows:

If *A* is a 1 × 1 matrix [*a*] (i.e. a scalar) then det *A* = *a*. Otherwise,

where *A _{ij}* is matrix

*A*with row

*i*and column

*j*removed.

Note that if *A* = , then we use the notation for det *A*.

**Excel Functions**: Excel provides the following function for calculating the determinant of a square matrix:

**MDETERM**(*A*): If *A* is a square array, then MDETERM(*A*) = det *A*. This is not an array function.

The supplemental function **DET**(*A*) provides equivalent functionality.

**Property 1**:

- det
*A*= det^{T}*A* - If
*A*is a diagonal matrix, then det*A*= the product of the elements on the main diagonal of A

Proof: Both of these properties are a simple consequence of Definition 1

**Property 2: ** = *ad * – *bc*

**Example 1**: Calculate det *A* where

From Definition 1 and Property 2 it follows that

Of course we can get the same answer by using Excel’s function MDETERM(*A*).

**Property 3**: If *A* and *B* are square matrices of the same size then det *AB* = det *A* ∙ det *B*

**Property 4**: A square matrix *A* is invertible if and only if det *A* ≠ 0. If *A* is invertible then

The first assertion is equivalent to saying that a square matrix *A* is singular if and only if det *A* = 0.

**Property 5**: Rules for evaluating determinants:

- The determinant of a triangular matrix is the product of the entries on the diagonal.
- If we interchange two rows, the determinant of the new matrix is the negative of the old one.
- If we multiply one row by a constant, the determinant of the new matrix is the determinant of the old one multiplied by the constant.
- If we add one row to another one multiplied by a constant, the determinant of the new matrix is the same as the old one.

**Observation**: The rules in Property 5 are sufficient to calculate the determinant of any square matrix. The idea is to transform the original matrix into a triangular matrix and then use rule 1 to calculate the value of the determinant.

We now present an algorithm based on Property 5 for calculating det *A*, where *A* = [*a _{ij}*] is an

*n*×

*n*matrix. Start by setting the value of the determinant to 1, and then perform steps 1 to

*n*as follows.

Step *k* – part 1(a): If *a _{kk}* ≠ 0, multiply the current value of the determinant by

*a*and then divide all the entries in row

_{kk}*k*by

*a*(rule 3 of Property 5).

_{kk}Step *k* – part 1(b): If *a _{kk}* = 0, exchange row

*k*with any row below it (i.e.

*k*<

*m*≤

*n*) for which

*a*≠ 0, multiply the current value of the determinant by -1 (rule 2) and then perform step 1(a) above. If no such row exists then terminate the algorithm and return the value of 0 for the determinant.

_{mk}Step *k* – part 2: For every row *m* below row *k*, add –*a _{mk}* times row

*k*to row

*m*(rule 4). This guarantees that

*a*= 0 for all

_{ij}*i*>

*k*and

*j*≤

*k*.

After the completion of step *n*, we will have a triangular matrix whose diagonal contains all 1s, and so by rule 1, the determinant is equal to the current value of the determinant.

**Example 2**: Using Property 5, find

We present the steps looking from left to right and then top to bottom in Figure 1. For each step the rule used is specified as well as the multiplier of the determinant calculated up to that point.

This shows that the determinant is -5, the same answer given when using Excel’s MDETERM function.

**Observation**: In step *k *– part 1(b) of the above procedure we exchange two rows if *a _{kk}* = 0. Given that we need to deal with round off errors, what happens if

*a*is small but not quite zero? In order to reduce the impact of round off errors, we should modify step

_{kk}*k*– part 1 as follows:

Step *k* – part 1: Find *m* ≥ *k *such that the absolute value of *a _{mk}* is largest. If this

*a*≈ 0 (i.e. |

_{mk}*a*|< ϵ where ϵ is some predefined small value) then terminate the procedure. If

_{mk}*m*>

*k*then exchange rows

*m*and

*k*.

**Definition 2**: A set of *n* linear equations

in *k* unknowns *x _{j}* can be viewed as the matrix equation

*AX = C*where

*A*is the

*n*×

*k*matrix [

*a*],

_{ij}*X*is the

*k*× 1 column vector [

*x*] and

_{j}*C*is the

*n*× 1 column vector [

*c*].

_{j}**Property 6**: If *A* is a square matrix (i.e. the number of equations is equal to the number of unknowns), the equation *AX = C* has a unique solution if and only if *A* is invertible (i.e. det *A* ≠ 0), and in this case the unique solution is given by *X = A ^{-1} C.*

**Property 7** (**Cramer’s Rule**): If the square matrix is invertible, the unique solution to is given by

where *A _{j}* is

*A*with the

*j*th column replaced by the entries of

*C*.

**Example 3**: Solve the following linear system using Cramer’s Rule:

In Figure 2, we calculate det *A* and det *A _{j}* for each

*j*.

It now follows that *x* = -6/9 = -2/3, y = 3/9 = 1/3 and *z* = 0/9 = 0.

Per Property 6, we can get the same result by calculating *A ^{-1} C*, which can be carried out in Excel using the formula =MMULT(MINVERSE(

*A*),

*C*). For Example 3, this yields

**Definition 3**: When *C* from Definition 2 is not the null matrix, then the linear equations are called **heterogeneous**. When *C* = *Ο* the linear equations are called **homogeneous**. In this case *Ο* is a solution of *AX* = *Ο*, called the **trivial solution**.

**Property 8**: If *A* is invertible then* X = Ο* is the only solution of *AX = Ο*.

Proof: This follows from Property 6.

**Observation**: When *A* is not invertible (i.e. det *A* = 0) any scalar multiple of a non-trivial solution to the homogeneous equation *AX = Ο* is also a solution. To find such a solution we can use the Gaussian Elimination method, a method which is similar to the one we used to calculate the determinant of a square matrix based on Property 5. This approach works for any *A* (whether square or not, and whether invertible or not).

**Definition 4**: If *A* is an *n* × *k* matrix and *B* is an *n* × *m* matrix then the **augmented matrix** *A*|*B* is an *n* × (*k*+*m*) matrix whose first columns are identical to the columns in *A* and whose remaining columns are identical to the columns in *B*.

**Property 9**: If *A′* and *C′* are derived from *A* and *C* based on any one of the following transformations, then the equations *AX = C* and *A ′X = C′* have the same solutions.

- Interchange of any two rows
- Multiplication of any row by a constant
- Addition of any row multiplied by a constant to another row

**Observation**: Typically we apply the above transformations to the augmented matrix *A*|*C*.

**Definition 5**: The **Gaussian Elimination Method** is a way of solving linear equations and is based on the transformations described in Property 9. Suppose *A* and *C* are as described in Definition 2.

Step 0 – set *i* = 1 and *j* = 1

We now apply the following series of transformation to the augmented matrix (step 1 through step *p* where *p* is the smaller of *n* and *k*):

Step *i* – part 1: Find *r* ≥* i* such that the absolute value of *a _{rj}* is largest. If

*a*≈ 0 (i.e. |

_{rj}*a*| < ϵ where ϵ is some predefined small value), then if

_{rj}*j*=

*n*terminate the procedure; otherwise replace

*j*by

*j*+ 1 and repeat step i. If not

*a*≈ 0 and

_{rj}*r*>

*k*then exchange rows

*r*and

*i*(rule 1).

Step *i* – part 2: Divide all the entries in row *i* by *a _{ij}* (rule 2).

Step *i* – part 3: For every row *r* below row *i*, add –*a _{rj}* times row

*i*to row

*r*(rule 3). This guarantees that

*a*= 0 for all

_{rc}*r*>

*i*and

*c*≤

*j*.

**Observation**: For non-homogeneous equations (i.e. *C *≠ *Ο*) there are three possibilities: there are an infinite number of solutions, there are no solutions or there is a unique solution. For homogeneous equations (i.e. *C* = *Ο*) there are two possibilities: there are an infinite number of solutions or there is a unique solution, namely the trivial solution where *X* = *Ο*.

**Example 4**: Solve the following linear system using Gaussian Elimination:

Figure 3 displays the steps in the Gaussian Elimination process for Example 4.

Since *A* transforms into the identity matrix we know that the transform of *C* is the unique solution to the system of linear equations, namely *x* = 0, y = 2 and *z* = -1. Note that we get the same result by calculating *X* = *A ^{-1} C*.

**Example 5**: Solve the following homogeneous linear system using Gaussian Elimination:

By Gaussian Elimination from Figure 4 we see that the only solution is the trivial solution:

**Example 6**: Solve the following homogeneous linear system using Gaussian Elimination:

This time Gaussian Elimination produces a row with all zeros (see Figure 5), and the number of non-zero rows = 2 < 3 unknowns. Thus there are an infinite number of solutions.

The solutions can take the form *x* = -2.5*t*, y = .5*t*, *z* = *t* for any value of *t*.

**Observation**: As we can see from the above examples, a homogeneous equation *AX = O*, where *A* is an *m* × *n* matrix, has a unique solution when there are *n* non-zero rows after performing Gaussian Elimination. Otherwise the equation has an infinite number of solutions.

**Real Statistics Excel Functions**: We provide the following supplemental Excel array function for carrying out the Gaussian Elimination procedure.

**ELIM**(R1, *prec*): an array function which outputs the results of Gaussian Elimination on the augmented matrix found in the array R1. The shape of the output is the same as the shape of R1.

**LINEQU**(R1, *prec*): an array function which returns an *n* × 1 column vector with the unique solution to equations defined by the augmented *m* × (*n*+1) matrix found in array R1; returns a vector consisting of #N/A! if there is no solution and a vector consisting of #NUM! if there are an infinite number of solutions.

By default, each of these functions assumes that an entry with absolute value less than 0.0001 is equivalent to zero. This is necessary since small values are not treated as zero in the Gaussian elimination algorithm described above. You can change this default value to something else by inserting a second parameter in either of these functions: e.g. ELIM(R1, *prec*) or LINEQU(R1, *prec*). Thus ELIM(R1) = ELIM(R1, 0.0001).

**Real Statistics Data Analysis Tool**: The **Solve Set of Linear Equations **data analysis tool contained in the Real Statistics Resource Pack provides equivalent functionality to LINEQU and ELIM. To use this tool, enter **Ctrl-m** and select **Solve Set of Linear Equations **from the menu. When a dialog box appears, fill in the **Input Range** (with the same range as R1 above). Selecting **Show solution only** is equivalent to LINEQU(R1), while not clicking on this option is equivalent to ELIM(R1).

**Observation**: Gaussian Elimination can also be used to invert a square *n* × *n* matrix *A* by applying the above procedure to *A*|*I _{n}*. If the procedure terminates before

*n*steps have been completed then

*A*is not invertible. If the procedure terminates after

*n*steps (in which case

*A′ = I*) then

_{n}*C′ = A*.

^{-1}**Example 7:** Use Gaussian Elimination to invert the matrix

The result is shown in Figure 6.

The fact that *A* transforms into the identity matrix indicates that *A* is invertible. The inverse is given by the transformation of the identity matrix, namely

This is the same result as using the Excel formula MINVERSE(*A*).

Great blog! Do you have any tips for aspiring writers?

I’m hoping to start my own website soon but I’m a little lost on everything.

Would you advise starting with a free platform like WordPress or go for a paid option? There are so many

options out there that I’m completely overwhelmed .. Any ideas?

Thanks a lot!

Ireej,

To gain some experience I started by putting the website up on the free WordPress platform. It worked for a while until I found that I needed some of the tools that are not available on the free platform.

Charles

As a side note – if we accept “Definition 1” then “Property 2” should be read as $\begin{vmatrix} a & b \\ c & d \end{vmatrix} = ad – bc$.

Not “ac – bd”.

Eugene,

Thanks for catching this typo. I have now corrected the formula on the webpage.

Charles

Really this site have help me a lot and i am really appreciate the efforts of all those who contributed to development of this site. I am looking forward to send me more problems and solutions on determinant and matrices to my email address to enable me study at home because, presently i am student.