Skip to main content

Section NM Nonsingular Matrices

In this section we specialize further and consider matrices with equal numbers of rows and columns, which when considered as coefficient matrices lead to systems with equal numbers of equations and variables. We will see in the second half of the course (Chapter D, Chapter E, Chapter LT, Chapter R) that these matrices are especially important.

Subsection NM Nonsingular Matrices

Our theorems will now establish connections between systems of equations (homogeneous or otherwise), augmented matrices representing those systems, coefficient matrices, constant vectors, the reduced row-echelon form of matrices (augmented and coefficient) and solution sets. Be very careful in your reading, writing and speaking about systems of equations, matrices and sets of vectors. A system of equations is not a matrix, a matrix is not a solution set, and a solution set is not a system of equations. Now would be a great time to review the discussion about speaking and writing mathematics in Proof Technique L.

Definition SQM. Square Matrix.

A matrix with \(m\) rows and \(n\) columns is square if \(m=n\text{.}\) In this case, we say the matrix has size \(n\text{.}\) To emphasize the situation when a matrix is not square, we will call it rectangular.

We can now present one of the central definitions of linear algebra.

Definition NM. Nonsingular Matrix.

Suppose \(A\) is a square matrix. Suppose further that the solution set to the homogeneous linear system of equations \(\linearsystem{A}{\zerovector}\) is \(\set{\zerovector}\text{,}\) in other words, the system has only the trivial solution. Then we say that \(A\) is a nonsingular matrix. Otherwise we say \(A\) is a singular matrix.

We can investigate whether any square matrix is nonsingular or not, no matter if the matrix is derived somehow from a system of equations or if it is simply a matrix. The definition says that to perform this investigation we must construct a very specific system of equations (homogeneous, with the matrix as the coefficient matrix) and look at its solution set. We will have theorems in this section that connect nonsingular matrices with systems of equations, creating more opportunities for confusion. Convince yourself now of two observations, (1) we can decide nonsingularity for any square matrix, and (2) the determination of nonsingularity involves the solution set for a certain homogeneous system of equations.

Notice that it makes no sense to call a system of equations nonsingular (the term does not apply to a system of equations), nor does it make any sense to call a \(5\times 7\) matrix singular (the matrix is not square).

Example HISAA shows that the coefficient matrix derived from Archetype A, specifically the \(3\times 3\) matrix

\begin{equation*} A=\begin{bmatrix} 1 & -1 & 2\\ 2 & 1 & 1\\ 1 & 1 & 0 \end{bmatrix} \end{equation*}

is a singular matrix since there are nontrivial solutions to the homogeneous system \(\homosystem{A}\text{.}\)

Example HUSAB shows that the coefficient matrix derived from Archetype B, specifically the \(3\times 3\) matrix

\begin{equation*} B=\begin{bmatrix} -7&-6&-12\\ 5&5&7\\ 1&0&4 \end{bmatrix} \end{equation*}

is a nonsingular matrix since the homogeneous system, \(\homosystem{B}\text{,}\) has only the trivial solution.

Notice that we will not discuss Example HISAD as being a singular or nonsingular coefficient matrix since the matrix is not square.

Being nonsingular is an important matrix property, and in such cases Sage contains commands that quickly and easily determine if the mathematical object does, or does not, have the property. The names of these types of methods universally begin with .is_, and these might be referred to as “predicates” or “queries.”. In the Sage notebook, define a simple matrix A, and then in a cell type A.is_, followed by pressing the TAB key rather than evaluating the cell. You will get a list of numerous properties that you can investigate for the matrix A. (Tab-completion will not work as advertised with the Sage Cell server.)

The other convention is to name these properties in a positive way, so the relevant command for nonsingular matrices is .is_singular(). We will redo Example S and Example NM. Note the use of not in the last compute cell.

The next theorem combines with our main computational technique (row reducing a matrix) to make it easy to recognize a nonsingular matrix. But first a definition.

Definition IM. Identity Matrix.

The \(m\times m\) identity matrix, \(I_m\text{,}\) is defined by

\begin{align*} \matrixentry{I_m}{ij}&= \begin{cases} 1 & \text{if }i=j\\ 0 & \text{if }i\neq j \end{cases} \end{align*}

for \(1\leq i,\,j\leq m\text{.}\)

The \(4\times 4\) identity matrix is

\begin{equation*} I_4= \begin{bmatrix} 1&0&0&0\\ 0&1&0&0\\ 0&0&1&0\\ 0&0&0&1 \end{bmatrix}\text{.} \end{equation*}

Notice that an identity matrix is square, and in reduced row-echelon form. Also, every column is a pivot column, and every possible pivot column appears once.

It is straightforward to create an identity matrix in Sage. Just specify the number system and the number of rows (which will equal the number of columns, so you do not specify that since it would be redundant). The number system can be left out, but the result will have entries from the integers (ZZ), which in this course is unlikely to be what you really want.

Notice that we do not use the now-familiar dot notation to create an identity matrix. What would we use the dot notation on anyway? For these reasons we call the identity_matrix() function a constructor, since it builds something from scratch, in this case a very particular type of matrix.

We mentioned above that an identity matrix is in reduced row-echelon form. What happens if we try to row-reduce a matrix that is already in reduced row-echelon form? By the uniqueness of the result, there should be no change. The following code illustrates this. Notice that = is used to assign an object to a new name, while == is used to test equality of two objects. I frequently make the mistake of forgetting the second equal sign when I mean to test equality.

(⇐) 

Suppose \(B\) is the identity matrix. When the augmented matrix \(\augmented{A}{\zerovector}\) is row-reduced, the result is \(\augmented{B}{\zerovector}=\augmented{I_n}{\zerovector}\text{.}\) The number of nonzero rows is equal to the number of variables in the linear system of equations \(\linearsystem{A}{\zerovector}\text{,}\) so \(n=r\) and Theorem FVCS gives \(n-r=0\) free variables. Thus, the homogeneous system \(\homosystem{A}\) has just one solution, which must be the trivial solution. This is exactly the definition of a nonsingular matrix (Definition NM).

(⇒) 

If \(A\) is nonsingular, then the homogeneous system \(\linearsystem{A}{\zerovector}\) has a unique solution, and has no free variables in the description of the solution set. The homogeneous system is consistent (Theorem HSC) so Theorem FVCS applies and tells us there are \(n-r\) free variables. Thus, \(n-r=0\text{,}\) and so \(n=r\text{.}\) So \(B\) has \(n\) pivot columns among its total of \(n\) columns. This is enough to force \(B\) to be the \(n\times n\) identity matrix \(I_n\) (see Exercise NM.T12).

Notice that since this theorem is an equivalence it will always allow us to determine if a matrix is either nonsingular or singular. Here are two examples of this, continuing our study of Archetype A and Archetype B.

We have the coefficient matrix for Archetype A and a row-equivalent matrix \(B\) in reduced row-echelon form.

\begin{equation*} A=\begin{bmatrix} 1 & -1 & 2\\ 2 & 1 & 1\\ 1 & 1 & 0 \end{bmatrix} \rref\begin{bmatrix} \leading{1} & 0 & 1\\ 0 & \leading{1} & -1\\ 0 & 0 & 0 \end{bmatrix}=B \end{equation*}

Since \(B\) is not the \(3\times 3\) identity matrix, Theorem NMRRI tells us that \(A\) is a singular matrix.

We have the coefficient matrix for Archetype B and a row-equivalent matrix \(B\) in reduced row-echelon form.

\begin{equation*} A=\begin{bmatrix} -7&-6&-12\\ 5&5&7\\ 1&0&4 \end{bmatrix} \rref\begin{bmatrix} \leading{1} & 0 & 0\\ 0 & \leading{1} & 0\\ 0 & 0 & \leading{1} \end{bmatrix} =B \end{equation*}

Since \(B\) is the \(3\times 3\) identity matrix, Theorem NMRRI tells us that \(A\) is a nonsingular matrix.

Subsection NSNM Null Space of a Nonsingular Matrix

Nonsingular matrices and their null spaces are intimately related, as the next two examples illustrate.

Given the singular coefficient matrix from Archetype A, the null space is the set of solutions to the homogeneous system of equations \(\homosystem{A}\text{,}\) which has a solution set and null space constructed in Example HISAA as an infinite set of vectors.

\begin{align*} A&=\begin{bmatrix} 1 & -1 & 2\\ 2 & 1 & 1\\ 1 & 1 & 0 \end{bmatrix} &\nsp{A} &=\setparts{\colvector{-x_3\\x_3\\x_3}}{x_3\in\complexes} \end{align*}

Given the nonsingular coefficient matrix from Archetype B, the solution set to the homogeneous system \(\homosystem{A}\) is constructed in Example HUSAB and contains only the trivial solution, so the null space of \(A\) has only a single element.

\begin{align*} A&=\begin{bmatrix} -7&-6&-12\\ 5&5&7\\ 1&0&4 \end{bmatrix} &\nsp{A} &=\set{\colvector{0\\0\\0}} \end{align*}

These two examples illustrate the next theorem, which is another equivalence.

The null space of a square matrix, \(A\text{,}\) is equal to the set of solutions to the homogeneous system, \(\homosystem{A}\text{.}\) A matrix is nonsingular if and only if the set of solutions to the homogeneous system, \(\linearsystem{A}{\zerovector}\text{,}\) has only a trivial solution. These two observations may be chained together to construct the two proofs necessary for each half of this theorem.

The next theorem pulls a lot of big ideas together. Theorem NMUS tells us that we can learn much about solutions to a system of linear equations with a square coefficient matrix by just examining a similar homogeneous system.

(⇐) 

The hypothesis for this half of the proof is that the system \(\linearsystem{A}{\vect{b}}\) has a unique solution for every choice of the constant vector \(\vect{b}\text{.}\) We will make a very specific choice for \(\vect{b}\text{:}\) \(\vect{b}=\zerovector\text{.}\) Then we know that the system \(\linearsystem{A}{\zerovector}\) has a unique solution. But this is precisely the definition of what it means for \(A\) to be nonsingular (Definition NM). That almost seems too easy! Notice that we have not used the full power of our hypothesis, but there is nothing that says we must use a hypothesis to its fullest.

(⇒) 

We assume that \(A\) is nonsingular of size \(n\times n\text{,}\) so we know there is a sequence of row operations that will convert \(A\) into the identity matrix \(I_n\) (Theorem NMRRI). Form the augmented matrix \(A^\prime=\augmented{A}{\vect{b}}\) and apply this same sequence of row operations to \(A^\prime\text{.}\) The result will be the matrix \(B^\prime=\augmented{I_n}{\vect{c}}\text{,}\) which is in reduced row-echelon form with \(r=n\text{.}\) Then the augmented matrix \(B^\prime\) represents the (extremely simple) system of equations \(x_i=\vectorentry{\vect{c}}{i}\text{,}\) \(1\leq i\leq n\text{.}\) The vector \(\vect{c}\) is clearly a solution, so the system is consistent (Definition CS). With a consistent system, we use Theorem FVCS to count free variables. We find that there are \(n-r=n-n=0\) free variables, and so we therefore know that the solution is unique. (This half of the proof was suggested by Asa Scherer.)

This theorem helps to explain part of our interest in nonsingular matrices. If a matrix is nonsingular, then no matter what vector of constants we pair it with, using the matrix as the coefficient matrix will always yield a linear system of equations with a solution, and the solution is unique. To determine if a matrix has this property (nonsingularity) it is enough to just solve one linear system, the homogeneous system with the matrix as coefficient matrix and the zero vector as the vector of constants (or any other vector of constants, see Exercise MM.T10).

Formulating the negation of the second part of this theorem is a good exercise. A singular matrix has the property that for some value of the vector \(\vect{b}\text{,}\) the system \(\linearsystem{A}{\vect{b}}\) does not have a unique solution (which means that it has no solution or infinitely many solutions). We will be able to say more about this case later (see the discussion following Theorem PSPHS).

Square matrices that are nonsingular have a long list of interesting properties, which we will start to catalog in the following, recurring, theorem. Of course, singular matrices will then have all of the opposite properties. The following theorem is a list of equivalences.

We want to understand just what is involved with understanding and proving a theorem that says several conditions are equivalent. So have a look at Proof Technique ME before studying the first in this series of theorems.

The statement that \(A\) is nonsingular is equivalent to each of the subsequent statements by, in turn, Theorem NMRRI, Theorem NMTNS and Theorem NMUS. So the statement of this theorem is just a convenient way to organize all these results.

Sage will create random matrices and vectors, sometimes with various properties. These can be very useful for quick experiments, and they are also useful for illustrating that theorems hold for any object satisfying the hypotheses of the theorem. But this will never replace a proof.

We will illustrate Theorem NME1 using Sage. We will use a variant of the random_matrix() constructor that uses the algorithm='unimodular' keyword. We will have to wait for Chapter D before we can give a full explanation, but for now, understand that this command will always create a square matrix that is nonsingular. Also realize that there are square nonsingular matrices which will never be the output of this command. In other words, this command creates elements of just a subset of all possible nonsingular matrices.

So we are using random matrices below to illustrate properties predicted by Theorem NME1. Execute the first command to create a random nonsingular matrix, and notice that we only have to mark the output of A as random for our automated testing process. After a few runs, notice that you can also edit the value of n to create matrices of different sizes. With a matrix A defined, run the next three cells, which by Theorem NME1 each always produce True as their output, no matter what value A has, so long as A is nonsingular. Read the code and try to determine exactly how they correspond to the parts of the theorem (some commentary along these lines follows).

The only portion of these commands that may be unfamilar is the last one. The command range(n) is incredibly useful, as it will create a list of the integers from 0 up to, but not including, n. (We saw this command briefly in Sage FDV.) So, for example, range(3) == [0,1,2] is True. Pivots are returned as a “tuple” which is very much like a list, except we cannot change the contents. We can see the difference by the way the tuple prints with parentheses ((,)) rather than brackets ([,]). We can convert a list to a tuple with the tuple() command, in order to make the comparison succeed.

How do we tell if the reduced row-echelon form of the augmented matrix of a system of equations represents a system with a unique solution? First, the system must be consistent, which by Theorem RCLS means the last column is not a pivot column. Then with a consistent system we need to insure there are no free variables. This happens if and only if the remaining columns are all pivot columns, according to Theorem FVCS, thus the test used in the last compute cell.

Finally, you may have wondered why we refer to a matrix as nonsingular when it creates systems of equations with single solutions (Theorem NMUS)! I have wondered the same thing. We will have an opportunity to address this just before Theorem NPNF, and again when we get to Theorem SMZD. Can you wait that long?

Reading Questions NM Reading Questions

1.

In your own words state the definition of a nonsingular matrix.

2.

What is the easiest way to recognize if a square matrix is nonsingular or not?

3.

Suppose we have a system of equations and its coefficient matrix is nonsingular. What can you say about the solution set for this system?

Exercises NM Exercises

In Exercises C30–C33 determine if the matrix is nonsingular or singular. Give reasons for your answer.

C30.
\begin{equation*} \begin{bmatrix} -3 & 1 & 2 & 8\\ 2 & 0 & 3 & 4\\ 1 & 2 & 7 & -4\\ 5 & -1 & 2 & 0 \end{bmatrix} \end{equation*}
Solution

The matrix row-reduces to

\begin{equation*} \begin{bmatrix} \leading{1}& 0 & 0 & 0\\ 0 & \leading{1} & 0 & 0\\ 0 & 0 & \leading{1} & 0\\ 0 & 0 & 0 & \leading{1} \end{bmatrix} \end{equation*}

which is the \(4\times 4\) identity matrix. By Theorem NMRRI the original matrix must be nonsingular.

C31.
\begin{equation*} \begin{bmatrix} 2 & 3 & 1 & 4\\ 1 & 1 & 1 & 0\\ -1 & 2 & 3 & 5\\ 1 & 2 & 1 & 3 \end{bmatrix} \end{equation*}
Solution

Row-reducing the matrix yields

\begin{equation*} \begin{bmatrix} \leading{1} & 0 & 0 & -2\\ 0 & \leading{1} & 0 & 3\\ 0 & 0 & \leading{1} & -1\\ 0 & 0 & 0 & 0 \end{bmatrix}\text{.} \end{equation*}

Since this is not the \(4\times 4\) identity matrix, Theorem NMRRI tells us the matrix is singular.

C32.
\begin{equation*} \begin{bmatrix} 9 & 3 & 2 & 4\\ 5 & -6 & 1 & 3\\ 4 & 1 & 3 & -5 \end{bmatrix} \end{equation*}
Solution

The matrix is not square, so neither term is applicable. See Definition NM, which is stated for just square matrices.

C33.
\begin{equation*} \begin{bmatrix} -1 & 2 & 0 & 3 \\ 1 & -3 & -2 & 4 \\ -2 & 0 & 4 & 3 \\ -3 & 1 & -2 & 3 \end{bmatrix} \end{equation*}
Solution

Theorem NMRRI tells us we can answer this question by simply row-reducing the matrix. Doing this we obtain

\begin{equation*} \begin{bmatrix} \leading{1} & 0 & 0 & 0 \\ 0 & \leading{1} & 0 & 0 \\ 0 & 0 & \leading{1} & 0 \\ 0 & 0 & 0 & \leading{1} \end{bmatrix}\text{.} \end{equation*}

Since the reduced row-echelon form of the matrix is the \(4\times 4\) identity matrix \(I_4\text{,}\) we know that \(B\) is nonsingular.

C40.

Each of the archetypes below is a system of equations with a square coefficient matrix, or is itself a square matrix. Determine if these matrices are nonsingular, or singular. Comment on the null space of each matrix.

Archetype A, Archetype B, Archetype F, Archetype K, Archetype L

C50.

Find the null space of the matrix \(E\) below.

\begin{align*} E&= \begin{bmatrix} 2 & 1 & -1 & -9 \\ 2 & 2 & -6 & -6 \\ 1 & 2 & -8 & 0 \\ -1 & 2 & -12 & 12 \end{bmatrix} \end{align*}
Solution

We form the augmented matrix of the homogeneous system \(\homosystem{E}\) and row-reduce the matrix.

\begin{align*} \begin{bmatrix} 2 & 1 & -1 & -9 & 0 \\ 2 & 2 & -6 & -6 & 0 \\ 1 & 2 & -8 & 0 & 0 \\ -1 & 2 & -12 & 12 & 0 \end{bmatrix} &\rref \begin{bmatrix} \leading{1} & 0 & 2 & -6 & 0 \\ 0 & \leading{1} & -5 & 3 & 0 \\ 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 \end{bmatrix} \end{align*}

We knew ahead of time that this system would be consistent (Theorem HSC), but we can now see there are \(n-r=4-2=2\) free variables, namely \(x_3\) and \(x_4\) since \(F=\set{3,4,5}\) (Theorem FVCS). Based on this analysis, we can rearrange the equations associated with each nonzero row of the reduced row-echelon form into an expression for the lone dependent variable as a function of the free variables. We arrive at the solution set to this homogeneous system, which is the null space of the matrix by Definition NSM,

\begin{gather*} \nsp{E}=\setparts{\colvector{-2x_3+6x_4\\5x_3-3x_4\\x_3\\x_4}}{x_3,\,x_4\in\complexes} \end{gather*}
M30.

Let \(A\) be the coefficient matrix of the system of equations below. Is \(A\) nonsingular or singular? Explain what you could infer about the solution set for the system based only on what you have learned about \(A\) being singular or nonsingular.

\begin{align*} -x_1+5x_2&=-8\\ -2x_1+5x_2+5x_3+2x_4&=9\\ -3x_1-x_2+3x_3+x_4&=3\\ 7x_1+6x_2+5x_3+x_4&=30 \end{align*}
Solution

We row-reduce the coefficient matrix of the system of equations.

\begin{align*} \begin{bmatrix} -1 & 5 & 0 & 0 \\ -2 & 5 & 5 & 2 \\ -3 & -1 & 3 & 1 \\ 7 & 6 & 5 & 1 \end{bmatrix} &\rref \begin{bmatrix} \leading{1} & 0 & 0 & 0 \\ 0 & \leading{1} & 0 & 0 \\ 0 & 0 & \leading{1} & 0 \\ 0 & 0 & 0 & \leading{1} \end{bmatrix} \end{align*}

Since the row-reduced version of the coefficient matrix is the \(4\times 4\) identity matrix, \(I_4\) (Definition IM), by Theorem NMRRI we know the coefficient matrix is nonsingular. According to Theorem NMUS we know that the system is guaranteed to have a unique solution, based only on the extra information that the coefficient matrix is nonsingular.

For Exercises M51–M52 say as much as possible about each system's solution set. Be sure to make it clear which theorems you are using to reach your conclusions.

M51.

6 equations in 6 variables, singular coefficient matrix.

Solution

Theorem NMRRI tells us that the coefficient matrix will not row-reduce to the identity matrix. So if we were to row-reduce the augmented matrix of this system of equations, we would not get a unique solution. So by Theorem PSSLS the remaining possibilities are no solutions, or infinitely many.

M52.

A system with a nonsingular coefficient matrix, not homogeneous.

Solution

Any system with a nonsingular coefficient matrix will have a unique solution by Theorem NMUS. If the system is not homogeneous, the solution cannot be the zero vector (Exercise HSE.T10).

T10.

Suppose that \(A\) is a square matrix, and \(B\) is a matrix in reduced row-echelon form that is row-equivalent to \(A\text{.}\) Prove that if \(A\) is singular, then the last row of \(B\) is a zero row.

Solution

Let \(n\) denote the size of the square matrix \(A\text{.}\) By Theorem NMRRI the hypothesis that \(A\) is singular implies that \(B\) is not the identity matrix \(I_n\text{.}\) If \(B\) has \(n\) pivot columns, then it would have to be \(I_n\text{,}\) so \(B\) must have fewer than \(n\) pivot columns. But the number of nonzero rows in \(B\) (\(r\)) is equal to the number of pivot columns as well. So the \(n\) rows of \(B\) have fewer than \(n\) nonzero rows, and \(B\) must contain at least one zero row. By Definition RREF, this row must be at the bottom of \(B\text{.}\)

A proof can also be formulated by first forming the contrapositive of the statement (Proof Technique CP) and proving this statement.

T12.

Using (Definition RREF) and (Definition IM) carefully, give a proof of the following equivalence: \(A\) is a square matrix in reduced row-echelon form where every column is a pivot column if and only if \(A\) is the identity matrix.

T30.

Suppose that \(A\) is a nonsingular matrix and \(A\) is row-equivalent to the matrix \(B\text{.}\) Prove that \(B\) is nonsingular.

Solution

Since \(A\) and \(B\) are row-equivalent matrices, consideration of the three row operations (Definition RO) will show that the augmented matrices, \(\augmented{A}{\zerovector}\) and \(\augmented{B}{\zerovector}\text{,}\) are also row-equivalent matrices. This says that the two homogeneous systems, \(\homosystem{A}\) and \(\homosystem{B}\) are equivalent systems. \(\homosystem{A}\) has only the zero vector as a solution (Definition NM), thus \(\homosystem{B}\) has only the zero vector as a solution. Finally, by Definition NM, we see that \(B\) is nonsingular.

Form a similar theorem replacing “nonsingular” by “singular” in both the hypothesis and the conclusion. Prove this new theorem with an approach just like the one above, and/or employ the result about nonsingular matrices in a proof by contradiction.

T31.

Suppose that \(A\) is a square matrix of size \(n\times n\) and that we know there is a single vector \(\vect{b}\in\complex{n}\) such that the system \(\linearsystem{A}{\vect{b}}\) has a unique solution. Prove that \(A\) is a nonsingular matrix. (Notice that this is very similar to Theorem NMUS, but is not exactly the same.)

Solution

Let \(B\) be the reduced row-echelon form of the augmented matrix \(\augmented{A}{\vect{b}}\text{.}\) Because the system \(\linearsystem{A}{\vect{b}}\) is consistent, we know by Theorem RCLS that the last column of \(B\) is not a pivot column. Suppose now that \(r\lt n\text{.}\) Then by Theorem FVCS the system would have infinitely many solutions. From this contradiction, we see that \(r=n\) and the first \(n\) columns of \(B\) are each pivot columns. Then the sequence of row operations that converts \(\augmented{A}{\vect{b}}\) to \(B\) will also convert \(A\) to \(I_n\text{.}\) Applying Theorem NMRRI we conclude that \(A\) is nonsingular.

T90.

Provide an alternative for the second half of the proof of Theorem NMUS, without appealing to properties of the reduced row-echelon form of the coefficient matrix. In other words, prove that if \(A\) is nonsingular, then \(\linearsystem{A}{\vect{b}}\) has a unique solution for every choice of the constant vector \(\vect{b}\text{.}\) Construct this proof without using Theorem REMEF or Theorem RREFU.

Solution

We assume \(A\) is nonsingular, and try to solve the system \(\linearsystem{A}{\vect{b}}\) without making any assumptions about \(\vect{b}\text{.}\) To do this we will begin by constructing a new homogeneous linear system of equations that looks very much like the original. Suppose \(A\) has size \(n\) (why must it be square?) and write the original system as

\begin{align*} a_{11}x_1+a_{12}x_2+a_{13}x_3+\dots+a_{1n}x_n&=b_1\\ a_{21}x_1+a_{22}x_2+a_{23}x_3+\dots+a_{2n}x_n&=b_2\\ a_{31}x_1+a_{32}x_2+a_{33}x_3+\dots+a_{3n}x_n&=b_3\\ \vdots&\tag{*}\\ a_{n1}x_1+a_{n2}x_2+a_{n3}x_3+\dots+a_{nn}x_n&=b_n\text{.} \end{align*}

Form the new, homogeneous system in \(n\) equations with \(n+1\) variables, by adding a new variable \(y\text{,}\) whose coefficients are the negatives of the constant terms.

\begin{align*} a_{11}x_1+a_{12}x_2+a_{13}x_3+\dots+a_{1n}x_n-b_1y&=0\\ a_{21}x_1+a_{22}x_2+a_{23}x_3+\dots+a_{2n}x_n-b_2y&=0\\ a_{31}x_1+a_{32}x_2+a_{33}x_3+\dots+a_{3n}x_n-b_3y&=0\\ \vdots&\tag{**}\\ a_{n1}x_1+a_{n2}x_2+a_{n3}x_3+\dots+a_{nn}x_n-b_ny&=0 \end{align*}

Since this is a homogeneous system with more variables than equations (\(m=n+1\gt n\)), Theorem HMVEI says that the system has infinitely many solutions. We will choose one of these solutions, any one of these solutions, so long as it is not the trivial solution. Write this solution as

\begin{align*} x_1=c_1&&x_2=c_2&&x_3=c_3&&\ldots&&x_n=c_n&&y=c_{n+1}\text{.} \end{align*}

We know that at least one value of the \(c_i\) is nonzero, but we will now show that in particular \(c_{n+1}\neq 0\text{.}\) We do this using a proof by contradiction (Proof Technique CD). So suppose the \(c_i\) form a solution as described, and in addition that \(c_{n+1}=0\text{.}\) Then we can write the \(i\)-th equation of system \((**)\) as,

\begin{align*} a_{i1}c_1+a_{i2}c_2+a_{i3}c_3+\dots+a_{in}c_n-b_i(0)&=0\\ \end{align*}

which becomes

\begin{align*} a_{i1}c_1+a_{i2}c_2+a_{i3}c_3+\dots+a_{in}c_n&=0\text{.} \end{align*}

Since this is true for each \(i\text{,}\) we have that \(x_1=c_1,\,x_2=c_2,\,x_3=c_3,\ldots,\,x_n=c_n\) is a solution to the homogeneous system \(\homosystem{A}\) formed with a nonsingular coefficient matrix. This means that the only possible solution is the trivial solution, so \(c_1=0,\,c_2=0,\,c_3=0,\,\ldots,\,c_n=0\text{.}\) So, assuming simply that \(c_{n+1}=0\text{,}\) we conclude that all of the \(c_i\) are zero. But this contradicts our choice of the \(c_i\) as not being the trivial solution to the system \((**)\text{.}\) So \(c_{n+1}\neq 0\text{.}\)

We now propose and verify a solution to the original system \((*)\text{.}\) Set

\begin{align*} x_1=\frac{c_1}{c_{n+1}}&&x_2=\frac{c_2}{c_{n+1}}&&x_3=\frac{c_3}{c_{n+1}}&&\ldots&&x_n=\frac{c_n}{c_{n+1}}\text{.} \end{align*}

Notice how it was necessary that we know that \(c_{n+1}\neq 0\) for this step to succeed. Now, evaluate the \(i\)-th equation of system \((*)\) with this proposed solution, and recognize in the third line that \(c_1\) through \(c_{n+1}\) appear as if they were substituted into the left-hand side of the \(i\)-th equation of system \((**)\text{.}\)

\begin{align*} &a_{i1}\frac{c_1}{c_{n+1}}+a_{i2}\frac{c_2}{c_{n+1}}+a_{i3}\frac{c_3}{c_{n+1}}+\dots+a_{in}\frac{c_n}{c_{n+1}}\\ &=\frac{1}{c_{n+1}}\left(a_{i1}c_1+a_{i2}c_2+a_{i3}c_3+\dots+a_{in}c_n\right)\\ &=\frac{1}{c_{n+1}}\left(a_{i1}c_1+a_{i2}c_2+a_{i3}c_3+\dots+a_{in}c_n-b_ic_{n+1}\right)+b_i\\ &=\frac{1}{c_{n+1}}\left(0\right)+b_i\\ &=b_i \end{align*}

Since this equation is true for every \(i\text{,}\) we have found a solution to system \((*)\text{.}\) To finish, we still need to establish that this solution is unique.

With one solution in hand, we will entertain the possibility of a second solution. So assume system \((*)\) has two solutions

\begin{align*} x_1=d_1&&x_2=d_2&&x_3=d_3&&\ldots&&x_n=d_n\\ x_1=e_1&&x_2=e_2&&x_3=e_3&&\ldots&&x_n=e_n\text{.} \end{align*}

Then,

\begin{align*} &\left(a_{i1}(d_1-e_1)+a_{i2}(d_2-e_2)+a_{i3}(d_3-e_3)+\dots+a_{in}(d_n-e_n)\right)\\ &=\left(a_{i1}d_1+a_{i2}d_2+a_{i3}d_3+\dots+a_{in}d_n\right)-\left(a_{i1}e_1+a_{i2}e_2+a_{i3}e_3+\dots+a_{in}e_n\right)\\ &=b_i-b_i\\ &=0\text{.} \end{align*}

This is the \(i\)-th equation of the homogeneous system \(\homosystem{A}\) evaluated with \(x_j=d_j-e_j\text{,}\) \(1\leq j\leq n\text{.}\) Since \(A\) is nonsingular, we must conclude that this solution is the trivial solution, and so \(0=d_j-e_j\text{,}\) \(1\leq j\leq n\text{.}\) That is, \(d_j=e_j\) for all \(j\) and the two solutions are identical, meaning any solution to \((*)\) is unique.

Notice that the proposed solution (\(x_i=\frac{c_i}{c_{n+1}}\)) appeared in this proof with no motivation whatsoever. This is just fine in a proof. A proof should convince you that a theorem is true. It is your job to read the proof and be convinced of every assertion. Questions like “Where did that come from?” or “How would I think of that?” have no bearing on the validity of the proof.