** Section 1.5 **

**Inverses of Square Matrices**

`> `

`> `

**Introduction**

Consider real numbers. Given the equation

*ax = b.*

If
*a *
is not equal to zero then we can solve the equation for
*x*

*a x = *
b or
*x*
=

* x*
=
or
* x*
=

What about linear systems of equations
*Ax = b*
?
*x =*
?

__Real Numbers__

*a*
= 1 =
*a *
or
* a *
= 1
* *
=
*a*

__Square Matrices__

* A*
=
* *
I
* *
=
*A*

**Coefficient Matrix **

Let us start with some examples of matrices associated with systems of linear equations with

an equal number of equations as unknowns.

**Example 1:**
The augmented matrix of a nonhomogeneous system of linear equations is

`> `
**AUG:=matrix([[1,3,6,7],[1,5,8,4],[3,13,22,15]]);**

Apply Gauss elimination to AUG :

`> `
**B:=gausselim(AUG);**

What does matrix B tell us about the solution set of the system of equations represented by AUG? Let us apply backsubstitution to solve the system:

`> `
**backsub(B);**

The system is consistent with infinitely many solutions.For this system, let us
**investigate **
the associated coefficient matrix:

`> `
**A:=matrix([[1,3,6],[1,5,8],[3,13,22]]);**

Apply Gauss elimination

`> `
**B1:=gausselim(A);**

`> `

This example shows that the underlying system has
__infinitely many solutions__
. Note that matrix

B
*1*
, in echelon form corresponding to the coefficient matrix A, has a zero row.

**Example 2:**
Alter the matrix AUG slightly .Consider the augmented matrix

`> `
**AUG:=matrix([[1,3,6,7],[1,5,8,4],[3,13,22,1]]);**

Apply Gauss elimination to AUG

`> `
**B:=gausselim(AUG);**

What does the matrix B tell us about the solution set of the system of equations represented by AUG?

Let us apply backsubstitution to solve the system

`> `
**backsub(B);**

Error, (in backsub) inconsistent system

For this system let us
**investigate **
the associated coefficient matrix. The coefficient matrix A is

`> `
**A:=matrix([[1,3,6],[1,5,8],[3,13,22]]);**

Let us apply Gauss elimination:

`> `
**B1:=gausselim(A);**

`> `

The slight modification in the augmented matrix leads to a system that is
__inconsistent__
. Note that the

matrix B
*1*
,in echelon form corresponding to the coefficient matrix A, does have also a zero row.

**Example 3:**
Now alter the matrix AUG of Example 2. Consider the augmented matrix

`> `
**AUG:=matrix([[1,3,6,7],[1,5,8,4],[3,13,2,1]]);**

Apply gauss elimination to AUG

`> `
**B:=gausselim(AUG);**

What does matrix B tell us about the solution set of the system of equations represented by AUG?

Let us apply backsubstitution to solve the system:

`> `
**backsub(B);**

`> `

For this system let us
**investigate **
the associated coefficient matrix:

`> `
**A:=matrix([[1,3,6],[1,5,8],[3,13,2]]);**

Apply Gauss elimination:

`> `
**B1:=gausselim(A);**

`> `

This last modification led to a system that has a
__unique solution__
. Note that the matrix
, in

echelon form corresponding to the coefficient matrix A, is not zero.

********************************************************************

In Examples 1 and 2 , the coefficient matrix associated with the augmented matrix

AUG determines whether the system represented by AUG: has a unique solution or not

Recall that in solving the scalar equation where and are real numbers,

one of the following rules must hold:

If and , the system has infinitely many solutions

If and b is not equal to zero, the system has no solution

If is not equal to zero, the system will have a unique solution given by

( is the multiplicative inverse of a).

For systems of linear equations
*Ax = B*
where
*A*
and
*B*
are matrices and
*x*
is the unknown

quantity, we have the following analogous situation:

As in
**Example 1**
, the system of linear equations has infinitely many solutions and in

**Example 2**
the system has no solution. In these two cases, we say the coefficient matrix

is "
**Singular**
". Note that the reduced matrix
in echelon form has at least one zero row.

In
** Example 3**
, the system of linear equations has a unique solution and the coefficient matrix

is "
**Nonsingular**
"; that is, the matrix A must have a multiplicative inverse. Note that the reduced

matrix in echelon form has no zero rows.

**Nonsingular Matrix**

**Example 4:**
Consider the matrix

`> `
**A:=matrix([[2,4,6],[4,6,7],[8,5,2]]);**

`> `

The inverse of A is

`> `
**B:=linalg[inverse](A);**

`> `

To verify this

`> `
**AB=multiply(A,B);**

`> `
**BA=multiply(B,A);**

`> `

**************************************************************************

An
*nxn*
matrix A is a
** nonsingular**
matrix

if and only if there is an
*nxn*
matrix B so that the product

AB = BA = I
*n*

where I
*n*
is the identity matrix. The matrix B is called the inverse of A. Otherwise, A is called a
** singular**
matrix.

If an
*nxn*
matrix has an inverse, then the inverse is unique

**************************************************************************

**Some Properties of a Nonsingular matrix**

**Example 5:**
Consider the matrix

`> `
**A:=matrix([[1,2,0],[0,1,1],[1,2,1]]);**

Its inverse is:

`> `
**A1:= linalg[inverse](A); **

`> `

Is the inverse of the inverse of A equal to A?

`> `
**inv(inv(A)) =linalg[inverse](A1);**

`> `

********************************************************************

If A is an
*nxn*
invertible matrix, the inverse of the inverse of A is A itself.

********************************************************************

Consider another matrix B and its product with matrix A:

`> `
**B:=matrix([[0,2,1],[1,1,0],[1,2,1]]); **

`> `
**C:=multiply(A,B);**

`> `

Take the inverse of the matrices A and B. Let
* *
=
and
.

`> `
**A1:= linalg[inverse](A); **

`> `
**B1:=linalg[inverse](B); **

What is the relation between the inverse of the product AB and the product of the inverses of

A and B?

Is the product of two nonsingular matrices nonsingular?

`> `
**inv(A)*inv(B) =multiply(A1,B1); **

`> `
**inv(B)*inv(A) =multiply(B1,A1);**

`> `
**inv(AB) = evalm(linalg[inverse](C));**

`> `

Compare and to the inverse, , of the product AB. What do you conclude?

*********************************************************************

If A and B are two
*nxn*
invertible matrices, then the inverse of AB is equal to the inverse

of B times the inverse of A; that is, =

*********************************************************************

Obtain the reduced echelon forms of the matrices A and B and find a
__condition under which the __

__multiplicative inverse exists.__

**Example 6:**
Take the following two matrices

`> `
**A:=matrix([[5,6,7],[5,7,2],[6,7,11]]);**

`> `
**B:=matrix([[5,6,7],[5,7,2],[15,19,16]]);**

`> `

The reduced echelon form of the matrices A and B are:

`> `
**A1:=rref(A);**

`> `
**B1:=rref(B);**

`> `

Is there a difference between the reduced row echelon form of the matrix A and the

reduced echelon form of the matrix B?

**Matrix A is row equivalent to the identity matrix**
. That is,

by applying a sequence of elementary row operations to matrix A

we obtained the identity matrix.Matrix B is not row equivalent to

the identity matrix.

Which one of them has an inverse?

`> `
**linalg[inverse](A);**

`> `

`> `
**linalg[inverse](B);**

Error, (in linalg[inverse]) singular matrix

`> `

Matrix A has an inverse while the inverse of the matrix B does not exist.

**************************************************************************

An
*nxn*
matrix A is invertible if and only if it is row equivalent to the identity matrix.

**************************************************************************

**Nonsingular Matrices and Systems of Linear Equations**

**Example 7:**
Consider the homogeneous system of linear equations whose augmented matrix is

`> `
**HS:=matrix([[5,6,7,0],[5,7,2,0],[6,7,11,0]]);**

`> `

The solution can be obtained by performing Gauss elimination and then applying the

backsubstitution algorithm.

`> `
**HS1:=gausselim(HS);**

`> `

`> `
**backsub(HS1);**

The homogeneous system of equations has the trivial solution. What can you say about the

coefficient matrix of
**HS**
? For example, does the inverse exist? The coefficient matrix of HS is

`> `
**A:=matrix([[5,6,7],[5,7,2],[6,7,11]]);**

Its inverse is

`> `
**linalg[inverse](A);**

`> `

************************************************************************

**Example 8: **
Consider a nonhomogeneous system
**NHS**
associated with an invertible matrix. A

nonhomogeneous system of linear equations has the following augmented matrix

`> `
**NHS:=matrix([[5,6,7,1],[5,7,2,15],[6,7,11,23]]);**

`> `

The solution can be obtained by performing Gauss elimination and then applying the

backsubstitution algorithm.

`> `
**NHS1:=gausselim(NHS);**

`> `
**backsub(NHS1);**

`> `

This system has a unique solution

*********************************************************************

If the coefficient matrix is nonsingular ( invertible), then the

nonhomogeneous system of linear equations has a unique

solution.

If the coefficient matrix is nonsingular( invertible), then the only

solution of the homogeneous system of linear equations is the

trivial solution.

*********************************************************************

__SUMMARY__

An
*nxn*
matrix A is non-singular ( invertible)

if and only if

A is row equivalent to the identity matrix

if and only if

The homogeneous system whose coefficient matrix is A has the trivial solution

if and only if

The nonhomogeneous system whose coefficient matrix is A has a unique solution

********************************************************************

**Finding the Inverse**

**The steps in the Procedure are**

1. Augment the given
*nxn*
matrix A with the identity matrix I
*n*
of the same dimension; that is,

form the matrix [A : I
*n*
]

2. Apply a sequence of elementary row operations to convert the matrix A to the identity

matrix. The same sequence of operations will convert the identity matrix to the desired inverse

of the matrix A. Why?

Apply the demonstration mode of inverse function from "linmat" to go through the steps

of finding the inverse of A, choosing the Gauss-Jordan procedure from the menu.

**Example 9: **
Consider the matrix A

`> `
**A:=matrix([[1,4,6],[6,4,9],[1,7,9]]);**

`> `

`> `
**inverse(A);**

1. Demonstration mode

2. Interactive mode

3. No Intermediate steps. Just give me the answer

`
`

`
`**1;**

Which method do you prefer?

1. Gauss-Jordan method

2. Adjoint method

<Select from 1; or 2;>

`
`**1;**

< type ; and press Enter to continue ?topic for help exit; to quit >

`
`**;**

< type ; and press Enter to continue ?topic for help exit; to quit >

`
`**;**

< type ; and press Enter to continue ?topic for help exit; to quit >

`
`**;**

< type ; and press Enter to continue ?topic for help exit; to quit >

`
`**;**

< type ; and press Enter to continue ?topic for help exit; to quit >

`
`**;**

< type ; and press Enter to continue ?topic for help exit; to quit >

`
`**;**

< type ; and press Enter to continue ?topic for help exit; to quit >

`
`**;**

< type ; and press Enter to continue ?topic for help exit; to quit >

`
`**;**

< type ; and press Enter to continue ?topic for help exit; to quit >

`
`**;**

< type ; and press Enter to continue ?topic for help exit; to quit >

`
`**;**

`> `
**;**

`> `

**Inverses and Systems of Linear Equations**

As noted in Example 8, if the coefficient matrix of a system of linear equations is invertible, then the system has a unique solution given by .

**Example 10:**
Consider the system of equations

`> `
**eq1:=x1+4*x2+6*x3=11: eq2:=6*x1+4*x2+9*x3=19:
eq3:=x1+7*x2+9*x3=35: print(eq1); print(eq2); print(eq3);**

`> `

This system can be written as Ax = b where:

`> `

`> `
**A:=matrix([[1,4,6],[6,4,9],[1,7,9]]);**

`> `

`> `
**b:=matrix([[11],[19],[35]]);**

Check if the inverse of the matrix A exists (use the nostep mode):

`> `

`> `
**inverse(A);**

1. Demonstration mode

2. Interactive mode

3. No Intermediate steps. Just give me the answer

`> `
**3;**

`> `

`> `

Thus the inverse of A is:

`> `
**B:=matrix([[-9/7,2/7,4/7],[-15/7,1/7,9/7],
[38/21,-1/7,-20/21]]);**

`> `

Since the inverse of the matrix A is the matrix B multiply, both sides of Ax = b by B to get the solution:

`> `
**x:=multiply(B,b);**

`> `

This solution coincides with the solution that we get directly from:

`> `
**solve({eq1,eq2,eq3},{x1,x2,x3});**

`> `

**Conditions for **
** **
**to exist**

********************************************************************

The following conditions for an
*n x n*
matrix are equivalent:

i.)
*A*
is invertible (non-singular), i.e.
exists.

ii.)
*A*
is row equivalent to the identity matrix
* I.*

iii.) The system
*Ax=b*
has a solution for every column vector
*b*
.

iv.) The matrix
*A *
can be expressed as a product of elementary matrices.

v.) The span of the column vectors of
*A*
is

********************************************************************

**Exercises**

1-13(odd), 16, 18, 19, 21, 23, 24, 35.