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

Solution 5.

Enter the matrix  A  and create a 4 by 4 identity matrix and store it in the variable Iden.

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


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

Remark.  In Mathematica the letter  [Graphics:../Images/GaussianJordanMod_gr_164.gif]  is a reserved word which is used for the complex unit, i.e.  [Graphics:../Images/GaussianJordanMod_gr_165.gif].  For this reason, we used the word  [Graphics:../Images/GaussianJordanMod_gr_166.gif]  for our "identity matrix"  [Graphics:../Images/GaussianJordanMod_gr_167.gif].  

Form the augmented matrix  [Graphics:../Images/GaussianJordanMod_gr_168.gif].  using the following steps.

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



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

Then perform Gauss-Jordan elimination.

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



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

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

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

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

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

Use Mathematica to get the inverse of  A  out of  this augmented matrix, and store it in the matrix  B.

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


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

Verify the solution.

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



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

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

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

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



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



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

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

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

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

Notice.  
Since integers were entered as the matrix and vector elements, precise arithmetic with fractions was performed for all of the computations.  If you prefer the numerical approximations, then enter decimal entries or make the final answer numerical with the command:

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



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

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

Aside.  
Compare with Mathematica's "Inverse" procedure for finding the inverse of a matrix.

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


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

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004