Example 3.  Find the eigenvalues and eigenvectors of the matrix  [Graphics:Images/EigenvaluesMod_gr_139.gif].  

Solution 3.

Find the characteristic polynomial and the eigenvalues.

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

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

 

 

Let us plot [Graphics:../Images/EigenvaluesMod_gr_142.gif] and see where the roots are located

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


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

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

Although this example has been "cooked up" so that the values are simple, we should be aware that a root finding method could be employed to find the eigenvalues.  For illustration, we can use the Newton-Raphson method.

[Graphics:../Images/EigenvaluesMod_gr_146.gif]
[Graphics:../Images/EigenvaluesMod_gr_147.gif]


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



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


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



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


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


Since we have solved for roots in previous modules, we will concentrate our effort on solving for the eigenvectors.
First, we shall automate the procedure for finding the roots of the characteristic polynomial, which is one way to find the eigenvalues.

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

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

 

 

Investigate the eigen-pair  [Graphics:../Images/EigenvaluesMod_gr_155.gif]

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


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

Introduce the free variables and find the eigenvector.

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

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

 

 

Verify the eigenpair.

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


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

Investigate the eigen-pair  [Graphics:../Images/EigenvaluesMod_gr_162.gif]

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


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

Introduce the free variables and find the eigenvector.

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

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

 

 

Verify the eigenpair.

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


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

Investigate the eigen-pair  [Graphics:../Images/EigenvaluesMod_gr_169.gif]

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


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

Introduce the free variables and find the eigenvector.

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

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

 

 

Verify the eigenpair.

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


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

The three eigen-pairs are:

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


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

We can compare this with the results obtained using Mathematica's Eigensystem procedure.

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


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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004