Example 1.  Use Gauss-Jordan elimination to find the inverse of the matrix  [Graphics:Images/InverseMatrixMod_gr_14.gif].  

Solution 1.

Enter the matrix  A  and the augmented matrix [Graphics:../Images/InverseMatrixMod_gr_15.gif].  

 

[Graphics:../Images/InverseMatrixMod_gr_16.gif]

Form the augmented matrix  [Graphics:../Images/InverseMatrixMod_gr_17.gif].  

 

 

[Graphics:../Images/InverseMatrixMod_gr_18.gif]


Not perform Gauss-Jordan elimination.

[Graphics:../Images/InverseMatrixMod_gr_19.gif]




[Graphics:../Images/InverseMatrixMod_gr_20.gif]

[Graphics:../Images/InverseMatrixMod_gr_21.gif]

[Graphics:../Images/InverseMatrixMod_gr_22.gif]

[Graphics:../Images/InverseMatrixMod_gr_23.gif]

[Graphics:../Images/InverseMatrixMod_gr_24.gif]

Get the inverse of  A  out of  this augmented matrix, and store it in the matrix  B.

 

[Graphics:../Images/InverseMatrixMod_gr_25.gif]

Verify the solution.

[Graphics:../Images/InverseMatrixMod_gr_26.gif]



[Graphics:../Images/InverseMatrixMod_gr_27.gif]

[Graphics:../Images/InverseMatrixMod_gr_28.gif]

We are done.  

Aside.  We can compare our result with Mathematica's built in procedure.  

[Graphics:../Images/InverseMatrixMod_gr_29.gif]


[Graphics:../Images/InverseMatrixMod_gr_30.gif]

[Graphics:../Images/InverseMatrixMod_gr_31.gif]

Aside.  If decimals are required in the solution then decimal need to be used when entering the matrix.

Enter the matrix  A  and the augmented matrix [Graphics:../Images/InverseMatrixMod_gr_32.gif].  

 

[Graphics:../Images/InverseMatrixMod_gr_33.gif]

Form the augmented matrix  [Graphics:../Images/InverseMatrixMod_gr_34.gif].  

 

[Graphics:../Images/InverseMatrixMod_gr_35.gif]

Then perform Gauss-Jordan elimination.

[Graphics:../Images/InverseMatrixMod_gr_36.gif]




[Graphics:../Images/InverseMatrixMod_gr_37.gif]

[Graphics:../Images/InverseMatrixMod_gr_38.gif]

[Graphics:../Images/InverseMatrixMod_gr_39.gif]

[Graphics:../Images/InverseMatrixMod_gr_40.gif]

[Graphics:../Images/InverseMatrixMod_gr_41.gif]


Get the inverse of  A  out of  this augmented matrix, and store it in the matrix  B.

Remark.  What if you thought that the inverse only had 6 decimal places of accuracy.

Then you would use the following matrix.

 

[Graphics:../Images/InverseMatrixMod_gr_42.gif]

Aside.  We can compare our result with Mathematica's built in procedure.  

[Graphics:../Images/InverseMatrixMod_gr_43.gif]


[Graphics:../Images/InverseMatrixMod_gr_44.gif]

[Graphics:../Images/InverseMatrixMod_gr_45.gif]

Verify the solution.

[Graphics:../Images/InverseMatrixMod_gr_46.gif]


[Graphics:../Images/InverseMatrixMod_gr_47.gif]

[Graphics:../Images/InverseMatrixMod_gr_48.gif]

This should not be too surprising if only 6 decimal places were available.

However, Mathematica carries out its internal computations in about 17 decimal places,  so there are really more digits that usually are not printed.

Get the inverse of  A  out of  this augmented matrix, and store it in the matrix  B.

 

[Graphics:../Images/InverseMatrixMod_gr_49.gif]

Verify the solution.

[Graphics:../Images/InverseMatrixMod_gr_50.gif]


[Graphics:../Images/InverseMatrixMod_gr_51.gif]

[Graphics:../Images/InverseMatrixMod_gr_52.gif]

Chop the values that are close to zero.

[Graphics:../Images/InverseMatrixMod_gr_53.gif]


[Graphics:../Images/InverseMatrixMod_gr_54.gif]

[Graphics:../Images/InverseMatrixMod_gr_55.gif]

This is a lot closer to the "true identity matrix."

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004