<?xml version="1.0" encoding="utf-8" standalone="no"?>
<!DOCTYPE document PUBLIC "-//CNX//DTD CNXML 0.5 plus MathML//EN" "http://cnx.rice.edu/cnxml/0.5/DTD/cnxml_mathml.dtd">
<document xmlns="http://cnx.rice.edu/cnxml" xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="m10736">

  <name>Eigenvectors and Eigenvalues</name>

   <metadata>
  <md:version>2.7</md:version>
  <md:created>2002/07/15</md:created>
  <md:revised>2006/07/19 16:36:50.576 GMT-5</md:revised>
  <md:authorlist>
      <md:author id="mjhaag">
      <md:firstname>Michael</md:firstname>
      
      <md:surname>Haag</md:surname>
      <md:email>mjhaag@rice.edu</md:email>
    </md:author>
      <md:author id="jrom">
      <md:firstname>Justin</md:firstname>
      
      <md:surname>Romberg</md:surname>
      <md:email>jrom@rice.edu</md:email>
    </md:author>
  </md:authorlist>

  <md:maintainerlist>
    <md:maintainer id="jrom">
      <md:firstname>Justin</md:firstname>
      
      <md:surname>Romberg</md:surname>
      <md:email>jrom@rice.edu</md:email>
    </md:maintainer>
    <md:maintainer id="richb">
      <md:firstname>Richard</md:firstname>
      <md:othername>G.</md:othername>
      <md:surname>Baraniuk</md:surname>
      <md:email>richb@rice.edu</md:email>
    </md:maintainer>
    <md:maintainer id="mjhaag">
      <md:firstname>Michael</md:firstname>
      
      <md:surname>Haag</md:surname>
      <md:email>mjhaag@rice.edu</md:email>
    </md:maintainer>
    <md:maintainer id="mariyah">
      <md:firstname>Mariyah</md:firstname>
      
      <md:surname>Poonawala</md:surname>
      <md:email>mariyah@rice.edu</md:email>
    </md:maintainer>
    <md:maintainer id="prash">
      <md:firstname>Prashant</md:firstname>
      
      <md:surname>Singh</md:surname>
      <md:email>prash@ece.rice.edu</md:email>
    </md:maintainer>
    <md:maintainer id="mhutch">
      <md:firstname>Matthew</md:firstname>
      
      <md:surname>Hutchinson</md:surname>
      <md:email>mhutch@rice.edu</md:email>
    </md:maintainer>
  </md:maintainerlist>
  
  <md:keywordlist>
    <md:keyword>determinant</md:keyword>
    <md:keyword>eigenfunction</md:keyword>
    <md:keyword>eigenvalue</md:keyword>
    <md:keyword>eigenvalues</md:keyword>
    <md:keyword>eigenvector</md:keyword>
    <md:keyword>fourier</md:keyword>
    <md:keyword>fourier series</md:keyword>
    <md:keyword>linear system</md:keyword>
  </md:keywordlist>

  <md:abstract>This module defines eigenvalues and eigenvectors and explains a method of finding them given a matrix.  These ideas are presented, along with many examples, in hopes of leading up to an understanding of the Fourier Series.
</md:abstract>
</metadata>

  <content>
    <para id="intro">
      In this section, our linear systems will be n×n matrices
      of complex numbers.  For a little background into some of the
      concepts that this module is based on, refer to the basics of
      <cnxn document="m10734" strength="7">linear algebra</cnxn>.
    </para>   

    <section id="eigvec">
      <name>Eigenvectors and Eigenvalues</name>
      <para id="p1_eigvec">
	Let
	<m:math display="inline">
	  <m:ci>A</m:ci>
	</m:math> be an n×n matrix, where 
	<m:math display="inline">
	  <m:ci>A</m:ci>
	</m:math> is a linear operator on vectors in 
	<m:math display="inline">
	  <m:apply>
	    <m:power/>
	    <m:complexes/>
	    <m:ci>n</m:ci>
	  </m:apply>
	</m:math>.

	<equation id="eq1">
	  <m:math>
	    <m:apply>
	      <m:eq/>
	      <m:apply>
		<m:times/>
		<m:ci>A</m:ci>
		<m:ci type="vector">x</m:ci>
	      </m:apply>
	      <m:ci type="vector">b</m:ci>
	    </m:apply>
	  </m:math>
	</equation>
	
	where 
	<m:math display="inline">
	  <m:ci type="vector">x</m:ci>
	</m:math> and
	<m:math display="inline">
	  <m:ci type="vector">b</m:ci>
	</m:math>
	are n×1 vectors (<cnxn target="fig1"/>).
      </para>

      <figure orient="vertical" id="fig1">
	<subfigure id="sub1_f1">
	  <media type="image/png" src="eigv_f1.png"/>
	</subfigure>
	<subfigure id="sub2_f1">
	  <media type="image/png" src="eigv_f2.png"/>
	</subfigure>
	<caption>
	  Illustration of linear system and vectors.
	</caption>
      </figure>
      
      
      <definition id="def_eigvec">
	<term>eigenvector</term>
	<meaning>
	  An eigenvector of 
	  <m:math display="inline">
	    <m:ci>A</m:ci>
	  </m:math> is a vector
	  <m:math display="inline">
	    <m:apply>
	      <m:in/>
	      <m:ci type="vector">v</m:ci>
	      <m:apply>
		<m:power/>
		<m:complexes/>
		<m:ci>n</m:ci>
	      </m:apply>
	    </m:apply>
	  </m:math> such that
	  
	  <equation id="eq2">
	    <m:math>
	      <m:apply>
		<m:eq/>
		<m:apply>
		  <m:times/>
		  <m:ci>A</m:ci>
		  <m:ci type="vector">v</m:ci>
		</m:apply>
		<m:apply>
		  <m:times/>
		  <m:ci>λ</m:ci>
		  <m:ci type="vector">v</m:ci>
		</m:apply>
	      </m:apply>
	    </m:math>
	  </equation>

	  where 
	  <m:math display="inline">
	    <m:ci>λ</m:ci> </m:math> is called the
	  corresponding <term>eigenvalue</term>.
	  <m:math display="inline">
	    <m:ci>A</m:ci>
	  </m:math> only changes the <emphasis>length</emphasis> of
	   <m:math display="inline">
	    <m:ci type="vector">v</m:ci>
	  </m:math>, not its direction.
	</meaning>
      </definition>

      <section id="geo_look">
	<name>Graphical Model</name>
	<para id="p1_geolk">
	  Through <cnxn target="eigf3"/> and <cnxn target="eigf4"/>,
	  let us look at the difference between <cnxn target="eq1" strength="7"/> and <cnxn target="eq2" strength="7"/>.
	</para>

	<figure id="eigf3">
	  <media type="image/png" src="eigv_f3.png"/>
	  <caption>
	    Represents <cnxn target="eq1" strength="7"/>, 
	    <m:math display="inline">
	      <m:apply>
		<m:eq/>
		<m:apply>
		  <m:times/>
		  <m:ci>A</m:ci>
		  <m:ci type="vector">x</m:ci>
		</m:apply>
		<m:ci type="vector">b</m:ci>
	      </m:apply>
	    </m:math>.
	  </caption>
	</figure>

	<para id="p2_geolk">
	  If 
	  <m:math display="inline">
	    <m:ci type="vector">v</m:ci>
	  </m:math> is an eigenvector of 
	  <m:math display="inline">
	    <m:ci>A</m:ci> </m:math>, then only its length changes.
	  See <cnxn target="eigf4"/> and notice how our vector's
	  length is simply scaled by our variable,
	  <m:math><m:ci>λ</m:ci></m:math>, called the
	  <term>eigenvalue</term>:
	</para>

	<figure id="eigf4">
	  <media type="image/png" src="eigv_f4.png"/>
	  <caption>
	    Represents <cnxn target="eq2" strength="7"/>,
	    <m:math display="inline">
	      <m:apply>
		<m:eq/>
		<m:apply>
		  <m:times/>
		  <m:ci>A</m:ci>
		  <m:ci type="vector">v</m:ci>
		</m:apply>
		<m:apply>
		  <m:times/>
		  <m:ci>λ</m:ci>
		  <m:ci type="vector">v</m:ci>
		</m:apply>
	      </m:apply>
	    </m:math>.
	  </caption>
	</figure>	

	<para id="p3_geolk">
	  <note type="note">
	    When dealing with a matrix  
	    <m:math display="inline">
	      <m:ci>A</m:ci> </m:math>, eigenvectors are the
	    <emphasis>simplest</emphasis> possible vectors to operate
	    on.
	  </note>
	</para>
	
      </section>
   
      <section id="eigv_eg">
	<name>Examples</name>
	
	<exercise id="exer1">
	  <problem>
	    <para id="pr1_exer1">
	      From inspection and understanding of eigenvectors, find
	      the two eigenvectors,
	      <m:math display="inline">
		<m:ci type="vector">
		  <m:msub>
		    <m:mi>v</m:mi>
		    <m:mn>1</m:mn>
		  </m:msub>
		</m:ci>
	      </m:math> and
	      <m:math display="inline">
		<m:ci type="vector">
		  <m:msub>
		    <m:mi>v</m:mi>
		    <m:mn>2</m:mn>
		  </m:msub>
		</m:ci>
	      </m:math>,
	      of 

	      <m:math display="block">
		<m:apply>
		  <m:eq/>
		  <m:ci>A</m:ci>
		  <m:apply>
		    <m:matrix>
		      <m:matrixrow>
			<m:cn>3</m:cn>
			<m:cn>0</m:cn>
		      </m:matrixrow>
		      <m:matrixrow>
			<m:cn>0</m:cn>
			<m:cn>-1</m:cn>
		      </m:matrixrow>		      
		    </m:matrix>
		  </m:apply>
		</m:apply>
	      </m:math>

	      Also, what are the corresponding eigenvalues, 
	      <m:math display="inline">
		<m:ci><m:msub>
		  <m:mi>λ</m:mi>
		  <m:mn>1</m:mn>
		</m:msub></m:ci>
	      </m:math> and
	      <m:math display="inline">
		<m:ci><m:msub>
		  <m:mi>λ</m:mi>
		  <m:mn>2</m:mn>
		</m:msub></m:ci>
	      </m:math>?  Do not worry if you are having problems
	      seeing these values from the information given so far,
	      we will look at more rigorous ways to find these values
	      soon.
	    </para>
	  </problem>
	  
	  <solution>
	    <para id="sol1_exer1">	    
	      The eigenvectors you found should be:

	      <m:math display="block">
		<m:apply>
		  <m:eq/>
		  <m:ci type="vector">
		    <m:msub>
		      <m:mi>v</m:mi>
		      <m:mn>1</m:mn>
		    </m:msub>
		  </m:ci>
		  <m:vector>
		    <m:cn>1</m:cn>
		    <m:cn>0</m:cn>
		  </m:vector>
		</m:apply>
	      </m:math>

	      <m:math display="block">
		<m:apply>
		  <m:eq/>
		  <m:ci type="vector">
		    <m:msub>
		      <m:mi>v</m:mi>
		      <m:mn>2</m:mn>
		    </m:msub>
		  </m:ci>
		  <m:vector>
		    <m:cn>0</m:cn>
		    <m:cn>1</m:cn>
		  </m:vector>
		</m:apply>
	      </m:math>

	      And the corresponding eigenvalues are

	      <m:math display="block">
		<m:apply>
		  <m:eq/>
		  <m:ci><m:msub>
		    <m:mi>λ</m:mi>
		    <m:mn>1</m:mn>
		  </m:msub></m:ci>
		  <m:cn>3</m:cn>
		</m:apply>
	      </m:math>
	      <m:math display="block">
		<m:apply>
		  <m:eq/>
		  <m:ci><m:msub>
		    <m:mi>λ</m:mi>
		    <m:mn>2</m:mn>
		  </m:msub></m:ci>
		  <m:cn>-1</m:cn>
		</m:apply>
	      </m:math>
	    </para>
	  </solution>
	</exercise>


	<exercise id="exer2">
	  <problem>
	    <para id="pr1_exer2">
	      Show that these two vectors,
	      <m:math display="block">
		<m:apply>
		  <m:eq/>
		  <m:ci type="vector">
		    <m:msub>
		      <m:mi>v</m:mi>
		      <m:mn>1</m:mn>
		    </m:msub>
		  </m:ci>
		  <m:vector>
		    <m:cn>1</m:cn>
		    <m:cn>1</m:cn>
		  </m:vector>
		</m:apply>
	      </m:math>

	      <m:math display="block">
		<m:apply>
		  <m:eq/>
		  <m:ci type="vector">
		    <m:msub>
		      <m:mi>v</m:mi>
		      <m:mn>2</m:mn>
		    </m:msub>
		  </m:ci>
		  <m:vector>
		    <m:cn>1</m:cn>
		    <m:cn>-1</m:cn>
		  </m:vector>
		</m:apply>
	      </m:math>
	      
	      are eigenvectors of 
	      <m:math display="inline">
		<m:ci>A</m:ci> 
	      </m:math>, where
	      
	      <m:math display="inline">
		<m:apply>
		  <m:eq/>
		  <m:ci>A</m:ci>
		  <m:matrix>
		    <m:matrixrow>
		      <m:cn>3</m:cn>
		      <m:cn>-1</m:cn>
		    </m:matrixrow>
		    <m:matrixrow>
		      <m:cn>-1</m:cn>
		      <m:cn>3</m:cn>
		    </m:matrixrow>		      
		  </m:matrix>
		</m:apply>
	      </m:math>.  Also, find the corresponding eigenvalues.
	    </para>
	  </problem>

	  <solution>
	    <para id="p1_sol1e2">
	      In order to prove that these two vectors are
	      eigenvectors, we will show that these statements meet
	      the requirements stated in the <cnxn target="def_eigvec" strength="7">definition</cnxn>.

	      <m:math display="block">
		<m:apply>
		  <m:eq/>
		  <m:apply>
		    <m:times/>
		    <m:ci>A</m:ci>
		    <m:ci type="vector">
		      <m:msub>
			<m:mi>v</m:mi>
			<m:mn>1</m:mn>
		      </m:msub>
		    </m:ci>
		  </m:apply>
		  <m:apply>
		    <m:times/>
		    <m:apply>
		      <m:matrix>
			<m:matrixrow>
			  <m:cn>3</m:cn>
			  <m:cn>-1</m:cn>
			</m:matrixrow>
			<m:matrixrow>
			  <m:cn>-1</m:cn>
			  <m:cn>3</m:cn>
			</m:matrixrow>		      
		      </m:matrix>
		    </m:apply>
		    <m:apply>
		      <m:vector>
			<m:cn>1</m:cn>
			<m:cn>1</m:cn>
		      </m:vector>
		    </m:apply>
		  </m:apply>
		  <m:apply>
		    <m:vector>
		      <m:cn>2</m:cn>
		      <m:cn>2</m:cn>
		    </m:vector>
		  </m:apply>
		</m:apply>
	      </m:math>

	       <m:math display="block">
		<m:apply>
		  <m:eq/>
		  <m:apply>
		    <m:times/>
		    <m:ci>A</m:ci>
		    <m:ci type="vector">
		      <m:msub>
			<m:mi>v</m:mi>
			<m:mn>2</m:mn>
		      </m:msub>
		    </m:ci>
		  </m:apply>
		  <m:apply>
		    <m:times/>
		    <m:apply>
		      <m:matrix>
			<m:matrixrow>
			  <m:cn>3</m:cn>
			  <m:cn>-1</m:cn>
			</m:matrixrow>
			<m:matrixrow>
			  <m:cn>-1</m:cn>
			  <m:cn>3</m:cn>
			</m:matrixrow>		      
		      </m:matrix>
		    </m:apply>
		    <m:apply>
		      <m:vector>
			<m:cn>1</m:cn>
			<m:cn>-1</m:cn>
		      </m:vector>
		    </m:apply>
		  </m:apply>
		  <m:apply>
		    <m:vector>
		      <m:cn>4</m:cn>
		      <m:cn>-4</m:cn>
		    </m:vector>
		  </m:apply>
		</m:apply>
	      </m:math>

	      These results show us that 
	      <m:math display="inline">
		<m:ci>A</m:ci>
	      </m:math>
	      only scales the two vectors (<foreign>i.e.</foreign>
	      changes their length) and thus it proves that <cnxn target="eq2" strength="8"/> holds true for the following
	      two eigenvalues that you were asked to find:
	      
	      <m:math display="block">
		<m:apply>
		  <m:eq/>
		  <m:ci><m:msub>
		    <m:mi>λ</m:mi>
		    <m:mn>1</m:mn>
		  </m:msub></m:ci>
		  <m:cn>2</m:cn>
		</m:apply>
	      </m:math>
	      <m:math display="block">
		<m:apply>
		  <m:eq/>
		  <m:ci><m:msub>
		    <m:mi>λ</m:mi>
		    <m:mn>2</m:mn>
		  </m:msub></m:ci>
		  <m:cn>4</m:cn>
		</m:apply>
	      </m:math>
	      	      
	      If you need more convincing, then one could also easily
	      graph the vectors and their corresponding product with
	      <m:math><m:ci>A</m:ci></m:math> to see that the results
	      are merely scaled versions of our original vectors,
	      <m:math display="inline">
		<m:ci type="vector">
		  <m:msub>
		    <m:mi>v</m:mi>
		    <m:mn>1</m:mn>
		  </m:msub>
		</m:ci>
	      </m:math> and
	      <m:math display="inline">
		<m:ci type="vector">
		  <m:msub>
		    <m:mi>v</m:mi>
		    <m:mn>2</m:mn>
		  </m:msub>
		</m:ci>
	      </m:math>.
	    </para>
	  </solution>

	</exercise>
      </section>
    </section>


    <section id="sec2">
      <name>Calculating Eigenvalues and Eigenvectors</name>
      <para id="p1_sec2">
	In the above examples, we relied on your understanding of the
	definition and on some basic observations to find and prove the
	values of the eigenvectors and eigenvalues.  However, as you
	can probably tell, finding these values will not always be
	that easy.  Below, we walk through a rigorous and mathematical
	approach at calculating the eigenvalues and eigenvectors of a
	matrix.
      </para>

      <section id="sub1_s2">
	<name>Finding Eigenvalues</name>
	<para id="p1_s1s2">
	  Find 
	  <m:math display="inline">
	    <m:apply>
	      <m:in/>
	      <m:ci>λ</m:ci>
	      <m:complexes/>
	    </m:apply>
	  </m:math> such that 
	  <m:math display="inline">
	    <m:apply>
	      <m:neq/>
	      <m:ci type="vector">v</m:ci>
	      <m:ci type="vector">0</m:ci>
	    </m:apply>
	  </m:math>,
	  where   
	  <m:math display="inline">
	    <m:ci type="vector">0</m:ci>
	  </m:math> is the "zero vector."  We will start with <cnxn target="eq2" strength="8"/>, and then work our way down
	  until we find a way to explicitly calculate  
	  <m:math display="inline">
	    <m:ci>λ</m:ci>
	  </m:math>.
	  
	  <m:math display="block">
	    <m:apply>
	      <m:eq/>
	      <m:apply>
		<m:times/>
		<m:ci>A</m:ci>
		<m:ci type="vector">v</m:ci>
	      </m:apply>
	      <m:apply>
		<m:times/>
		<m:ci>λ</m:ci>
		<m:ci type="vector">v</m:ci>
	      </m:apply>
	    </m:apply>
	  </m:math>

	  <m:math display="block">
	    <m:apply>
	      <m:eq/>
	      <m:apply>
		<m:minus/>
		<m:apply>
		  <m:times/>
		  <m:ci>A</m:ci>
		  <m:ci type="vector">v</m:ci>
		</m:apply>
		<m:apply>
		  <m:times/>
		  <m:ci>λ</m:ci>
		  <m:ci type="vector">v</m:ci>
		</m:apply>
	      </m:apply>
	      <m:cn>0</m:cn>
	    </m:apply>	  
	  </m:math>

	  <m:math display="block">
	    <m:apply>
	      <m:eq/>
	      <m:apply>
		<m:times/>
		<m:apply>
		  <m:minus/>
		  <m:ci>A</m:ci>
		  <m:apply>
		    <m:times/>
		    <m:ci>λ</m:ci>
		    <m:ci>I</m:ci>
		  </m:apply>
		</m:apply>
		<m:ci type="vector">v</m:ci>
	      </m:apply>
	      <m:cn>0</m:cn>
	    </m:apply>	  
	  </m:math>

	  In the previous step, we used the fact that 
	  <m:math display="block">
	    <m:apply>
	      <m:eq/>
	      <m:apply>
		<m:times/>
		<m:ci>λ</m:ci>
		<m:ci type="vector">v</m:ci>
	      </m:apply>
	      <m:apply>
		<m:times/>
		<m:ci>λ</m:ci>
		<m:ci>I</m:ci>
		<m:ci type="vector">v</m:ci>
	      </m:apply>
	    </m:apply>
	  </m:math>
	  
	  where <m:math><m:ci>I</m:ci></m:math> is the identity
	  matrix.

	  <m:math display="block">
	    <m:apply>
	      <m:eq/>
	      <m:ci>I</m:ci>
	      <m:apply>
		<m:matrix>
		  <m:matrixrow>
		    <m:cn>1</m:cn>
		    <m:cn>0</m:cn>
		    <m:ci>…</m:ci>
		    <m:cn>0</m:cn>
		  </m:matrixrow>
		  <m:matrixrow>
		    <m:cn>0</m:cn>
		    <m:cn>1</m:cn>
		    <m:ci>…</m:ci>
		    <m:cn>0</m:cn>
		  </m:matrixrow>
		  <m:matrixrow>
		    <m:cn>0</m:cn>
		    <m:cn>0</m:cn>
		    <m:ci>⋱</m:ci>
		    <m:ci>⋮</m:ci>
		  </m:matrixrow>
		  <m:matrixrow>
		    <m:cn>0</m:cn>
		    <m:ci>…</m:ci>
		    <m:ci>…</m:ci>
		    <m:cn>1</m:cn>
		  </m:matrixrow>
		</m:matrix>
	      </m:apply>
	    </m:apply>
	  </m:math>	  
	  
	  So, 
	  <m:math display="inline">
	    <m:apply>
	      <m:minus/>
	      <m:ci>A</m:ci>
	      <m:apply>
		<m:times/>
		<m:ci>λ</m:ci>
		<m:ci>I</m:ci>
	      </m:apply>
	    </m:apply>
	  </m:math> is just a new matrix.
	</para>


	<example id="eg1_eva">
	  <para id="p1_eg1eva">
	    Given the following matrix, 
	    <m:math><m:ci>A</m:ci></m:math>, then we can find our new
	    matrix, 
	    <m:math display="inline">
	      <m:apply>
		<m:minus/>
		<m:ci>A</m:ci>
		<m:apply>
		  <m:times/>
		  <m:ci>λ</m:ci>
		  <m:ci>I</m:ci>
		</m:apply>
	      </m:apply>
	    </m:math>.

	  <m:math display="block">
	    <m:apply>
	      <m:eq/>
	      <m:ci>A</m:ci>
	      <m:apply>
		<m:matrix>
		  <m:matrixrow>
		    <m:apply>
		      <m:selector/>
		      <m:ci>a</m:ci>
		      <m:cn>1</m:cn>
		      <m:cn>1</m:cn>
		    </m:apply>
		    <m:apply>
		      <m:selector/>
		      <m:ci>a</m:ci>
		      <m:cn>1</m:cn>
		      <m:cn>2</m:cn>
		    </m:apply>
		  </m:matrixrow>
		  <m:matrixrow>
		    <m:apply>
		      <m:selector/>
		      <m:ci>a</m:ci>
		      <m:cn>2</m:cn>
		      <m:cn>1</m:cn>
		    </m:apply>
		    <m:apply>
		      <m:selector/>
		      <m:ci>a</m:ci>
		      <m:cn>2</m:cn>
		      <m:cn>2</m:cn>
		    </m:apply>
		  </m:matrixrow>
		</m:matrix>
	      </m:apply>
	    </m:apply>
	  </m:math>	    
	    
	    <m:math display="block">
	      <m:apply>
		<m:eq/>
		<m:apply>
		  <m:minus/>
		  <m:ci>A</m:ci>
		  <m:apply>
		    <m:times/>
		    <m:ci>λ</m:ci>
		    <m:ci>I</m:ci>
		  </m:apply>
		</m:apply>
		<m:apply>
		  <m:matrix>
		    <m:matrixrow>
		      <m:apply>
			<m:minus/>
		    <m:apply>
		      <m:selector/>
		      <m:ci>a</m:ci>
		      <m:cn>1</m:cn>
		      <m:cn>1</m:cn>
		    </m:apply>
			<m:ci>λ</m:ci>
		      </m:apply>
		    <m:apply>
		      <m:selector/>
		      <m:ci>a</m:ci>
		      <m:cn>1</m:cn>
		      <m:cn>2</m:cn>
		    </m:apply>
		    </m:matrixrow>
		    <m:matrixrow>
		    <m:apply>
		      <m:selector/>
		      <m:ci>a</m:ci>
		      <m:cn>2</m:cn>
		      <m:cn>1</m:cn>
		    </m:apply>
		      <m:apply>
			<m:minus/>
		    <m:apply>
		      <m:selector/>
		      <m:ci>a</m:ci>
		      <m:cn>2</m:cn>
		      <m:cn>2</m:cn>
		    </m:apply>
			<m:ci>λ</m:ci>
		      </m:apply>
		    </m:matrixrow>
		  </m:matrix>
		</m:apply>
	      </m:apply>
	    </m:math>
	    
	  </para>
	</example>

	<para id="p2_s2s2">
	  If 
	  <m:math display="inline">
	    <m:apply>
	      <m:eq/>
	      <m:apply>
		<m:times/>
		<m:apply>
		  <m:minus/>
		  <m:ci>A</m:ci>
		  <m:apply>
		    <m:times/>
		    <m:ci>λ</m:ci>
		    <m:ci>I</m:ci>
		  </m:apply>
		</m:apply>
		<m:ci type="vector">v</m:ci>
	      </m:apply>
	      <m:cn>0</m:cn>
	    </m:apply>	  
	  </m:math> for some 
	  <m:math display="inline">
	    <m:apply>
	      <m:neq/>
	      <m:ci type="vector">v</m:ci>
	      <m:cn>0</m:cn>
	    </m:apply>
	  </m:math>, then
	  <m:math display="inline">
	    <m:apply>
	      <m:minus/>
	      <m:ci>A</m:ci>
	      <m:apply>
		<m:times/>
		<m:ci>λ</m:ci>
		<m:ci>I</m:ci>
	      </m:apply>
	    </m:apply>
	  </m:math> is <emphasis>not invertible</emphasis>.  This
	  means:

	  <m:math display="block">
	    <m:apply>
	      <m:eq/>
	      <m:apply>
		<m:determinant/>
		<m:apply>
		  <m:minus/>
		  <m:ci>A</m:ci>
		  <m:apply>
		    <m:times/>
		    <m:ci>λ</m:ci>
		    <m:ci>I</m:ci>
		  </m:apply>
		</m:apply>
	      </m:apply>
	      <m:cn>0</m:cn>
	    </m:apply>
	  </m:math>

	  This determinant (shown directly above) turns out to be a
	  polynomial expression (of order
	  <m:math><m:ci>n</m:ci></m:math>).  Look at the examples
	  below to see what this means.
	</para>

	<example id="eg2_eva">
	  <para id="p1_eg2eva">
	    Starting with matrix <m:math><m:ci>A</m:ci></m:math>
	    (shown below), we will find the polynomial expression,
	    where our eigenvalues will be the dependent variable.

	    <m:math display="block">
	      <m:apply>
		<m:eq/>
		<m:ci>A</m:ci>
		<m:apply>
		  <m:matrix>
		    <m:matrixrow>
		      <m:cn>3</m:cn>
		      <m:cn>-1</m:cn>
		    </m:matrixrow>
		    <m:matrixrow>
		      <m:cn>-1</m:cn>
		      <m:cn>3</m:cn>
		    </m:matrixrow>		      
		  </m:matrix>
		</m:apply>
	      </m:apply>
	    </m:math>
	    
	    <m:math display="block">
	      <m:apply>
		<m:eq/>
		<m:apply>
		  <m:minus/>
		  <m:ci>A</m:ci>
		  <m:apply>
		    <m:times/>
		    <m:ci>λ</m:ci>
		    <m:ci>I</m:ci>
		  </m:apply>
		</m:apply>
		<m:apply>
		  <m:matrix>
		    <m:matrixrow>
		      <m:apply>
			<m:minus/>
			<m:cn>3</m:cn>
			<m:ci>λ</m:ci>
		      </m:apply>
		      <m:cn>-1</m:cn>
		    </m:matrixrow>
		    <m:matrixrow>
		      <m:cn>-1</m:cn>
		      <m:apply>
			<m:minus/>
			<m:cn>3</m:cn>
			<m:ci>λ</m:ci>
		      </m:apply>
		    </m:matrixrow>		      
		  </m:matrix>
		</m:apply>
	      </m:apply>
	    </m:math>

	    <m:math display="block">
	      <m:apply>
		<m:eq/>
		<m:apply>
		  <m:determinant/>
		  <m:apply>
		    <m:minus/>
		    <m:ci>A</m:ci>
		    <m:apply>
		      <m:times/>
		      <m:ci>λ</m:ci>
		      <m:ci>I</m:ci>
		    </m:apply>
		  </m:apply>
		</m:apply>
		<m:apply>
		  <m:minus/>
		  <m:apply>
		    <m:power/>
		    <m:apply>
		      <m:minus/>
		      <m:cn>3</m:cn>
		      <m:ci>λ</m:ci>
		    </m:apply>
		    <m:cn>2</m:cn>
		  </m:apply>
		  <m:apply>
		    <m:power/>
		    <m:apply>
		      <m:cn>-1</m:cn>
		    </m:apply>
		    <m:cn>2</m:cn>
		  </m:apply>
		</m:apply>
		<m:apply>
		  <m:plus/>
		  <m:apply>
		    <m:minus/>
		    <m:apply>
		      <m:power/>
		      <m:ci>λ</m:ci>
		      <m:cn>2</m:cn>
		    </m:apply>
		    <m:apply>
		      <m:times/>
		      <m:cn>6</m:cn>
		      <m:ci>λ</m:ci>
		    </m:apply>
		  </m:apply>
		  <m:cn>8</m:cn>
		</m:apply>		
	      </m:apply>
	    </m:math>

	    <m:math display="block">
	      <m:apply>
		<m:eq/>
		<m:ci>λ</m:ci>
		<m:apply>
		  <m:set>
		    <m:cn>2</m:cn>
		    <m:cn>4</m:cn>
		  </m:set>
		</m:apply>
	      </m:apply>
	    </m:math>
	  </para>
	</example>


	<example id="eg3_eva">
	  <para id="p1_eg3eva">
	    Starting with matrix <m:math><m:ci>A</m:ci></m:math>
	    (shown below), we will find the polynomial expression,
	    where our eigenvalues will be the dependent variable.

	        <m:math display="block">
	      <m:apply>
		<m:eq/>
		<m:ci>A</m:ci>
		<m:apply>
		  <m:matrix>
		    <m:matrixrow>
		    <m:apply>
		      <m:selector/>
			<m:ci>a</m:ci>
			<m:cn>1</m:cn>
			<m:cn>1</m:cn>
		    </m:apply>
		    <m:apply>
		      <m:selector/>
			<m:ci>a</m:ci>
			<m:cn>1</m:cn>
			<m:cn>2</m:cn>
		    </m:apply>
		    </m:matrixrow>
		      <m:matrixrow>
		    <m:apply>
		      <m:selector/>
			<m:ci>a</m:ci>
			<m:cn>2</m:cn>
			<m:cn>1</m:cn>
		    </m:apply>
		    <m:apply>
		      <m:selector/>
			<m:ci>a</m:ci>
			<m:cn>2</m:cn>
			<m:cn>2</m:cn>
		    </m:apply>
		    </m:matrixrow>
		  </m:matrix>
		</m:apply>
	      </m:apply>
	    </m:math>	    

	    <m:math display="block">
	      <m:apply>
		<m:eq/>
		<m:apply>
		  <m:minus/>
		  <m:ci>A</m:ci>
		  <m:apply>
		    <m:times/>
		    <m:ci>λ</m:ci>
		    <m:ci>I</m:ci>
		  </m:apply>
		</m:apply>
		<m:apply>
		  <m:matrix>
		    <m:matrixrow>
		      <m:apply>
			<m:minus/>
		    <m:apply>
		      <m:selector/>
			<m:ci>a</m:ci>
			<m:cn>1</m:cn>
			<m:cn>1</m:cn>
		    </m:apply>
			<m:ci>λ</m:ci>
		      </m:apply>
		    <m:apply>
		      <m:selector/>
			<m:ci>a</m:ci>
			<m:cn>1</m:cn>
			<m:cn>2</m:cn>
		    </m:apply>
		    </m:matrixrow>
		    <m:matrixrow>
		    <m:apply>
		      <m:selector/>
			<m:ci>a</m:ci>
			<m:cn>2</m:cn>
			<m:cn>1</m:cn>
		    </m:apply>
		       <m:apply>
			<m:minus/>
		    <m:apply>
		      <m:selector/>
			<m:ci>a</m:ci>
			<m:cn>2</m:cn>
			<m:cn>2</m:cn>
		    </m:apply>
			<m:ci>λ</m:ci>
		      </m:apply>
		    </m:matrixrow>
		  </m:matrix>
		</m:apply>
	      </m:apply>
	    </m:math>	    

	    <m:math display="block">
	      <m:apply>
		<m:eq/>
		<m:apply>
		  <m:determinant/>
		  <m:apply>
		    <m:minus/>
		    <m:ci>A</m:ci>
		    <m:apply>
		      <m:times/>
		      <m:ci>λ</m:ci>
		      <m:ci>I</m:ci>
		    </m:apply>
		  </m:apply>
		</m:apply>
		<m:apply>
		  <m:plus/>
		  <m:apply>
		    <m:minus/>
		    <m:apply>
		      <m:minus/>
		      <m:apply>
			<m:power/>
			<m:ci>λ</m:ci>
			<m:cn>2</m:cn>
		      </m:apply>
		      <m:apply>
			<m:times/>
			<m:apply>
			  <m:plus/>
		      <m:apply>
		      <m:selector/>
			<m:ci>a</m:ci>
			<m:cn>1</m:cn>
			<m:cn>1</m:cn>
		    </m:apply>
		      <m:apply>
		      <m:selector/>
			<m:ci>a</m:ci>
			<m:cn>2</m:cn>
			<m:cn>2</m:cn>
		    </m:apply>
			</m:apply>
			<m:ci>λ</m:ci>
		      </m:apply>
		    </m:apply>
		    <m:apply>
		      <m:times/>
		    <m:apply>
		      <m:selector/>
			<m:ci>a</m:ci>
			<m:cn>2</m:cn>
			<m:cn>1</m:cn>
		    </m:apply>
		    <m:apply>
		      <m:selector/>
			<m:ci>a</m:ci>
			<m:cn>1</m:cn>
			<m:cn>2</m:cn>
		    </m:apply>
		    </m:apply>
		  </m:apply>
		  <m:apply>
		    <m:times/>
		    <m:apply>
		      <m:selector/>
			<m:ci>a</m:ci>
			<m:cn>1</m:cn>
			<m:cn>1</m:cn>
		    </m:apply>
		    <m:apply>
		      <m:selector/>
			<m:ci>a</m:ci>
			<m:cn>2</m:cn>
			<m:cn>2</m:cn>
		    </m:apply>
		  </m:apply>
		</m:apply>
	      </m:apply>
	    </m:math>
	  </para>
	</example>

	<para id="p3_s3s2">
	  If you have not already noticed it, calculating the
	  eigenvalues is equivalent to calculating the roots of

	  <m:math display="block">
	    <m:apply>
	      <m:eq/>
	      <m:apply>
		<m:determinant/>
		<m:apply>
		  <m:minus/>
		  <m:ci>A</m:ci>
		  <m:apply>
		    <m:times/>
		    <m:ci>λ</m:ci>
		    <m:ci>I</m:ci>
		  </m:apply>
		</m:apply>
	      </m:apply>
	      <m:apply>
		<m:plus/>
		<m:apply>
		  <m:times/>
		  <m:ci><m:msub>
		    <m:mi>c</m:mi>
		    <m:mi>n</m:mi>
		  </m:msub></m:ci>
		  <m:apply>
		    <m:power/>
		    <m:ci>λ</m:ci>
		    <m:ci>n</m:ci>
		  </m:apply>
		</m:apply>
		<m:apply>
		  <m:times/>
		  <m:ci><m:msub>
		    <m:mi>c</m:mi>
		    <m:mrow>
		      <m:mi>n</m:mi>
		      <m:mo>−</m:mo>
		      <m:mn>1</m:mn>
		    </m:mrow>
		  </m:msub></m:ci>
		  <m:apply>
		    <m:power/>
		    <m:ci>λ</m:ci>
		    <m:apply>
		      <m:minus/>
		      <m:ci>n</m:ci>
		      <m:cn>1</m:cn>
		    </m:apply>
		  </m:apply>
		</m:apply>
		<m:ci>…</m:ci>
		<m:apply>
		  <m:times/>
		  <m:ci><m:msub>
		    <m:mi>c</m:mi>
		    <m:mn>1</m:mn>
		  </m:msub></m:ci>
		  <m:ci>λ</m:ci>
		</m:apply>
		<m:ci><m:msub>
		  <m:mi>c</m:mi>
		  <m:mn>0</m:mn>
		</m:msub></m:ci>
	      </m:apply>
	      <m:cn>0</m:cn>
	    </m:apply>
	  </m:math>
	  
	  <note type="conclusion">
	    Therefore, by simply using calculus to solve for the roots
	    of our polynomial we can easily find the eigenvalues of our
	    matrix.
	  </note>

	</para>
      </section>
      
      <section id="sub2_s2">
	<name>Finding Eigenvectors</name>
	<para id="p1_s2s2">
	  Given an eigenvalue, 
	  <m:math>
	    <m:ci><m:msub>
	      <m:mi>λ</m:mi>
	      <m:mi>i</m:mi>
	    </m:msub></m:ci>
	  </m:math>, the associated eigenvectors are given by

	  <m:math display="block">
	    <m:apply>
	      <m:eq/>
	      <m:apply>
		<m:times/>
		<m:ci>A</m:ci>
		<m:ci type="vector">v</m:ci>
	      </m:apply>
	      <m:apply>
		<m:times/>
		<m:ci><m:msub>
		  <m:mi>λ</m:mi>
		  <m:mi>i</m:mi>
		</m:msub></m:ci>
		<m:ci type="vector">v</m:ci>
	      </m:apply>
	    </m:apply>
	  </m:math>
	  
	    <m:math display="block">
	    <m:apply>
	      <m:eq/>
	      <m:apply>
		<m:times/>
		<m:ci>A</m:ci>
		<m:vector>
		  <m:ci type="vector">
		    <m:msub>
		      <m:mi>v</m:mi>
		      <m:mn>1</m:mn>
		    </m:msub>
		  </m:ci>
		  <m:ci>⋮</m:ci>
		  <m:ci type="vector">
		    <m:msub>
		      <m:mi>v</m:mi>
		      <m:mi>n</m:mi>
		    </m:msub>
		  </m:ci>
		</m:vector>
	      </m:apply>
	      <m:vector>
		<m:apply>
		  <m:times/>
		  <m:ci><m:msub>
		    <m:mi>λ</m:mi>
		    <m:mn>1</m:mn>
		  </m:msub></m:ci>
		  <m:ci type="vector">
		    <m:msub>
		      <m:mi>v</m:mi>
		      <m:mn>1</m:mn>
		    </m:msub>
		  </m:ci>
		</m:apply>
		<m:ci>⋮</m:ci>
		<m:apply>
		  <m:times/>
		  <m:ci><m:msub>
		    <m:mi>λ</m:mi>
		    <m:mi>n</m:mi>
		  </m:msub></m:ci>
		  <m:ci type="vector">
		    <m:msub>
		      <m:mi>v</m:mi>
		      <m:mi>n</m:mi>
		    </m:msub>
		  </m:ci>
		</m:apply>
	      </m:vector>
	    </m:apply>
	  </m:math>

	  set of <m:math><m:ci>n</m:ci></m:math> equations with
	  <m:math><m:ci>n</m:ci></m:math> unknowns.  Simply
	  <emphasis>solve the <m:math><m:ci>n</m:ci></m:math>
	  equations</emphasis> to find the eigenvectors.
	</para>
      </section>
    </section>


    <section id="sec3">
      <name>Main Point</name>
      <para id="p1_sec3">
	Say the eigenvectors of <m:math><m:ci>A</m:ci></m:math>,
	<m:math display="inline">
	  <m:apply>
	    <m:set>
	      <m:ci type="vector">
		<m:msub>
		  <m:mi>v</m:mi>
		  <m:mn>1</m:mn>
		</m:msub>
	      </m:ci>
	      <m:ci type="vector">
		<m:msub>
		  <m:mi>v</m:mi>
		  <m:mn>2</m:mn>
		</m:msub>
	      </m:ci>
	      <m:ci>…</m:ci>
	      <m:ci type="vector">
		<m:msub>
		  <m:mi>v</m:mi>
		  <m:mi>n</m:mi>
		</m:msub>
	      </m:ci>
	    </m:set>
	  </m:apply>
	</m:math>, 
	<cnxn document="m10734" target="span_sec">span</cnxn> 
	
	<m:math>
	  <m:apply>
	    <m:power/>
	    <m:complexes/>
	    <m:ci>n</m:ci>
	  </m:apply>
	</m:math>, meaning 
	<m:math display="inline">
	  <m:apply>
	    <m:set>
	      <m:ci type="vector">
		<m:msub>
		  <m:mi>v</m:mi>
		  <m:mn>1</m:mn>
		</m:msub>
	      </m:ci>
	      <m:ci type="vector">
		<m:msub>
		  <m:mi>v</m:mi>
		  <m:mn>2</m:mn>
		</m:msub>
	      </m:ci>
	      <m:ci>…</m:ci>
	      <m:ci type="vector">
		<m:msub>
		  <m:mi>v</m:mi>
		  <m:mi>n</m:mi>
		</m:msub>
	      </m:ci>
	    </m:set>
	  </m:apply>
	</m:math> are <cnxn document="m10734" target="lin_ind" strength="7">linearly independent</cnxn> and we can write any 
	<m:math display="inline">
	  <m:apply>
	    <m:in/>
	    <m:ci type="vector">x</m:ci>
	    <m:apply>
	      <m:power/>
	      <m:complexes/>
	      <m:ci>n</m:ci>
	    </m:apply>
	  </m:apply>
	</m:math> as

	<equation id="eq3">
	  <m:math>
	    <m:apply>
	      <m:eq/>
	      <m:ci type="vector">x</m:ci>
	      <m:apply>
		<m:plus/>
		<m:apply>
		  <m:times/>
		  <m:ci><m:msub>
		    <m:mi>α</m:mi>
		    <m:mn>1</m:mn>
		  </m:msub></m:ci>
		  <m:ci type="vector">
		    <m:msub>
		      <m:mi>v</m:mi>
		      <m:mn>1</m:mn>
		    </m:msub>
		  </m:ci>
		</m:apply>
		<m:apply>
		  <m:times/>
		  <m:ci><m:msub>
		    <m:mi>α</m:mi>
		    <m:mn>2</m:mn>
		  </m:msub></m:ci>
		  <m:ci type="vector">
		    <m:msub>
		      <m:mi>v</m:mi>
		      <m:mn>2</m:mn>
		    </m:msub>
		  </m:ci>
		</m:apply>
		<m:ci>…</m:ci>
		<m:apply>
		  <m:times/>
		  <m:ci><m:msub>
		    <m:mi>α</m:mi>
		    <m:mi>n</m:mi>
		  </m:msub></m:ci>
		  <m:ci type="vector">
		    <m:msub>
		      <m:mi>v</m:mi>
		      <m:mi>n</m:mi>
		    </m:msub>
		  </m:ci>
		</m:apply>
	      </m:apply>
	    </m:apply>
	  </m:math>
	</equation>

	where 
	<m:math display="inline">
	  <m:apply>
	    <m:in/>
	    <m:apply>
	      <m:set>
		<m:ci type="vector">
		  <m:msub>
		    <m:mi>α</m:mi>
		    <m:mn>1</m:mn>
		  </m:msub>
		</m:ci>
		<m:ci type="vector">
		  <m:msub>
		    <m:mi>α</m:mi>
		    <m:mn>2</m:mn>
		  </m:msub>
		</m:ci>
		<m:ci>…</m:ci>
		<m:ci type="vector">
		  <m:msub>
		    <m:mi>α</m:mi>
		    <m:mi>n</m:mi>
		  </m:msub>
		</m:ci>
	      </m:set>
	    </m:apply>
	    <m:complexes/>
	  </m:apply>
	</m:math>.  All that we are doing is rewriting <m:math><m:ci type="vector">x</m:ci></m:math> in terms of eigenvectors of
	<m:math><m:ci>A</m:ci></m:math>.  Then,

	<m:math display="block">
	  <m:apply>
	    <m:eq/>
	    <m:apply>
	      <m:times/>
	      <m:ci>A</m:ci>
	      <m:ci type="vector">x</m:ci>
	    </m:apply>
	    <m:apply>
	      <m:times/>
	      <m:ci>A</m:ci>
	      <m:apply>
		<m:plus/>
		<m:apply>
		  <m:times/>
		  <m:ci><m:msub>
		    <m:mi>α</m:mi>
		    <m:mn>1</m:mn>
		  </m:msub></m:ci>
		  <m:ci type="vector">
		    <m:msub>
		      <m:mi>v</m:mi>
		      <m:mn>1</m:mn>
		    </m:msub>
		  </m:ci>
		</m:apply>
		<m:apply>
		  <m:times/>
		  <m:ci><m:msub>
		    <m:mi>α</m:mi>
		    <m:mn>2</m:mn>
		  </m:msub></m:ci>
		  <m:ci type="vector">
		    <m:msub>
		      <m:mi>v</m:mi>
		      <m:mn>2</m:mn>
		    </m:msub>
		  </m:ci>
		</m:apply>
		<m:ci>…</m:ci>
		<m:apply>
		  <m:times/>
		  <m:ci><m:msub>
		    <m:mi>α</m:mi>
		    <m:mi>n</m:mi>
		  </m:msub></m:ci>
		  <m:ci type="vector">
		    <m:msub>
		      <m:mi>v</m:mi>
		      <m:mi>n</m:mi>
		    </m:msub>
		  </m:ci>
		</m:apply>
	      </m:apply>
	    </m:apply>
	  </m:apply>
	</m:math>

	<m:math display="block">
	  <m:apply>
	    <m:eq/>
	    <m:apply>
	      <m:times/>
	      <m:ci>A</m:ci>
	      <m:ci type="vector">x</m:ci>
	    </m:apply>
	    <m:apply>
	      <m:plus/>
	      <m:apply>
		<m:times/>
		<m:ci><m:msub>
		  <m:mi>α</m:mi>
		  <m:mn>1</m:mn>
		</m:msub></m:ci>
		<m:ci>A</m:ci>
		<m:ci type="vector">
		  <m:msub>
		    <m:mi>v</m:mi>
		    <m:mn>1</m:mn>
		  </m:msub>
		</m:ci>
	      </m:apply>
	      <m:apply>
		<m:times/>
		<m:ci><m:msub>
		  <m:mi>α</m:mi>
		  <m:mn>2</m:mn>
		</m:msub></m:ci>
		<m:ci>A</m:ci>
		<m:ci type="vector">
		  <m:msub>
		    <m:mi>v</m:mi>
		    <m:mn>2</m:mn>
		  </m:msub>
		</m:ci>
	      </m:apply>
	      <m:ci>…</m:ci>
	      <m:apply>
		<m:times/>
		<m:ci><m:msub>
		  <m:mi>α</m:mi>
		  <m:mi>n</m:mi>
		</m:msub></m:ci>
		<m:ci>A</m:ci>
		<m:ci type="vector">
		  <m:msub>
		    <m:mi>v</m:mi>
		    <m:mi>n</m:mi>
		  </m:msub>
		</m:ci>
	      </m:apply>
	    </m:apply>
	  </m:apply>
	</m:math>

	
	<m:math display="block">
	  <m:apply>
	    <m:eq/>
	    <m:apply>
	      <m:times/>
	      <m:ci>A</m:ci>
	      <m:ci type="vector">x</m:ci>
	    </m:apply>
	    <m:apply>
	      <m:plus/>
	      <m:apply>
		<m:times/>
		<m:ci><m:msub>
		  <m:mi>α</m:mi>
		  <m:mn>1</m:mn>
		</m:msub></m:ci>
		<m:ci><m:msub>
		  <m:mi>λ</m:mi>
		  <m:mn>1</m:mn>
		</m:msub></m:ci>
		<m:ci type="vector">
		  <m:msub>
		    <m:mi>v</m:mi>
		    <m:mn>1</m:mn>
		  </m:msub>
		</m:ci>
	      </m:apply>
	      <m:apply>
		<m:times/>
		<m:ci><m:msub>
		  <m:mi>α</m:mi>
		  <m:mn>2</m:mn>
		</m:msub></m:ci>
		<m:ci><m:msub>
		  <m:mi>λ</m:mi>
		  <m:mn>2</m:mn>
		</m:msub></m:ci>
		<m:ci type="vector">
		  <m:msub>
		    <m:mi>v</m:mi>
		    <m:mn>2</m:mn>
		  </m:msub>
		</m:ci>
	      </m:apply>
	      <m:ci>…</m:ci>
	      <m:apply>
		<m:times/>
		<m:ci><m:msub>
		  <m:mi>α</m:mi>
		  <m:mi>n</m:mi>
		</m:msub></m:ci>
		<m:ci><m:msub>
		  <m:mi>λ</m:mi>
		  <m:mi>n</m:mi>
		</m:msub></m:ci>
		<m:ci type="vector">
		  <m:msub>
		    <m:mi>v</m:mi>
		    <m:mi>n</m:mi>
		  </m:msub>
		</m:ci>
	      </m:apply>
	    </m:apply>
	    <m:ci>b</m:ci>
	  </m:apply>
	</m:math>

	Therefore we can write,

	<m:math display="block">
	  <m:apply>
	    <m:eq/>
	    <m:ci type="vector">x</m:ci>
	    <m:apply>
	      <m:sum/>
	      <m:domainofapplication>
		<m:ci>i</m:ci>
	      </m:domainofapplication>
	      <m:apply>
		<m:times/>
		<m:ci><m:msub>
		  <m:mi>α</m:mi>
		  <m:mi>i</m:mi>
		</m:msub></m:ci>
		<m:ci type="vector">
		  <m:msub>
		    <m:mi>v</m:mi>
		    <m:mi>i</m:mi>
		  </m:msub>
		</m:ci>
	      </m:apply>
	    </m:apply>
	  </m:apply>
	</m:math>

	and this leads us to the following depicted system:	  
      </para>
      
      <figure id="eigf5">
	<media type="image/png" src="eigv_sys.png"/>
	<caption>
	  Depiction of system where we break our vector, <m:math><m:ci type="vector">x</m:ci></m:math>, into a sum of its
	  eigenvectors.
	</caption>
      </figure>
      
      <para id="pf_mp">
	where in <cnxn target="eigf5"/> we have,
	<m:math display="block">
	  <m:apply>
	    <m:eq/>
	    <m:ci>b</m:ci>
	    <m:apply>
	      <m:sum/>
	      <m:domainofapplication>
		<m:ci>i</m:ci>
	      </m:domainofapplication>
	      <m:apply>
		<m:times/>
		<m:ci><m:msub>
		  <m:mi>α</m:mi>
		  <m:mi>i</m:mi>
		</m:msub></m:ci>
		<m:ci><m:msub>
		  <m:mi>λ</m:mi>
		  <m:mi>i</m:mi>
		</m:msub></m:ci>
		<m:ci type="vector">
		  <m:msub>
		    <m:mi>v</m:mi>
		    <m:mi>i</m:mi>
		  </m:msub>
		</m:ci>
	      </m:apply>
	    </m:apply>
	  </m:apply>
	</m:math>

	<note type="Main Point">
	  By breaking up a vector, <m:math><m:ci type="vector">x</m:ci></m:math>, into a combination of
	  eigenvectors, the calculation of 
	  <m:math display="inline">
	    <m:apply>
	      <m:times/>
	      <m:ci>A</m:ci>
	      <m:ci type="vector">x</m:ci>
	    </m:apply>
	  </m:math> is broken into "easy to swallow" pieces.
	</note>
      </para>

    </section>
    
    <section id="sec4">
      <name>Practice Problem</name>
   	
      <exercise id="exer_fin">
	<problem>
	  <para id="pr1_ef">
	    For the following matrix, <m:math><m:ci>A</m:ci></m:math> and
	    vector, <m:math><m:ci type="vector">x</m:ci></m:math>, solve
	    for their product.  Try solving it using two different
	    methods: directly and using eigenvectors. 

	    <m:math display="block">
	      <m:apply>
		<m:eq/>
		<m:ci>A</m:ci>
		<m:apply>
		  <m:matrix>
		    <m:matrixrow>
		      <m:cn>3</m:cn>
		      <m:cn>-1</m:cn>
		    </m:matrixrow>
		    <m:matrixrow>
		      <m:cn>-1</m:cn>
		      <m:cn>3</m:cn>
		    </m:matrixrow>		      
		  </m:matrix>
		</m:apply>
	      </m:apply>
	    </m:math>

	    <m:math display="block">
	      <m:apply>
		<m:eq/>
		<m:ci type="vector">x</m:ci>
		<m:apply>
		  <m:vector>
		    <m:cn>5</m:cn>
		    <m:cn>3</m:cn>
		  </m:vector>
		</m:apply>
	      </m:apply>
	    </m:math>
	  </para>
	</problem>

	<solution>
	  <para id="sol1_ef">
	    <emphasis>Direct Method</emphasis> (use basic matrix
	    multiplication) 

	    <m:math display="block">
	      <m:apply>
		<m:eq/>
		<m:apply>
		  <m:times/>
		  <m:ci>A</m:ci>
		  <m:ci type="vector">x</m:ci>
		</m:apply>
		<m:apply>
		  <m:times/>
		  <m:apply>
		    <m:matrix>
		      <m:matrixrow>
			<m:cn>3</m:cn>
			<m:cn>-1</m:cn>
		      </m:matrixrow>
		      <m:matrixrow>
			<m:cn>-1</m:cn>
			<m:cn>3</m:cn>
		      </m:matrixrow>		      
		    </m:matrix>
		  </m:apply>
		  <m:apply>
		    <m:vector>
		      <m:cn>5</m:cn>
		      <m:cn>3</m:cn>
		    </m:vector>
		  </m:apply>
		</m:apply>
		<m:apply>
		  <m:vector>
		    <m:cn>12</m:cn>
		    <m:cn>4</m:cn>
		  </m:vector>
		</m:apply>
	      </m:apply>
	    </m:math>

	    <emphasis>Eigenvectors</emphasis> (use the eigenvectors
	    and eigenvalues we found earlier for this same matrix)
	    
	    <m:math display="block">
	      <m:apply>
		<m:eq/>
		<m:ci type="vector">
		  <m:msub>
		    <m:mi>v</m:mi>
		    <m:mn>1</m:mn>
		  </m:msub>
		</m:ci>
		<m:vector>
		  <m:cn>1</m:cn>
		  <m:cn>1</m:cn>
		</m:vector>
	      </m:apply>
	    </m:math>

	    <m:math display="block">
	      <m:apply>
		<m:eq/>
		<m:ci type="vector">
		  <m:msub>
		    <m:mi>v</m:mi>
		    <m:mn>2</m:mn>
		  </m:msub>
		</m:ci>
		<m:vector>
		  <m:cn>1</m:cn>
		  <m:cn>-1</m:cn>
		</m:vector>
	      </m:apply>
	    </m:math>

	    <m:math display="block">
	      <m:apply>
		<m:eq/>
		<m:ci><m:msub>
		  <m:mi>λ</m:mi>
		  <m:mn>1</m:mn>
		</m:msub></m:ci>
		<m:cn>2</m:cn>
	      </m:apply>
	    </m:math>
	    <m:math display="block">
	      <m:apply>
		<m:eq/>
		<m:ci><m:msub>
		  <m:mi>λ</m:mi>
		  <m:mn>2</m:mn>
		</m:msub></m:ci>
		<m:cn>4</m:cn>
	      </m:apply>
	    </m:math>
	    
	    As shown in <cnxn target="eq3" strength="8"/>, we want to
	    represent <m:math><m:ci type="vector">x</m:ci></m:math> as
	    a sum of its scaled eigenvectors.  For this case, we have:

	    <m:math display="block">
	      <m:apply>
		<m:eq/>
		<m:ci type="vector">x</m:ci>
		<m:apply>
		  <m:plus/>
		  <m:apply>
		    <m:times/>
		    <m:cn>4</m:cn>
		    <m:ci type="vector">
		      <m:msub>
			<m:mi>v</m:mi>
			<m:mn>1</m:mn>
		      </m:msub>
		    </m:ci>
		  </m:apply>
		  <m:ci type="vector">
		    <m:msub>
		      <m:mi>v</m:mi>
		      <m:mn>2</m:mn>
		    </m:msub>
		  </m:ci>
		</m:apply>
	      </m:apply>
	    </m:math>

	     <m:math display="block">
	      <m:apply>
		<m:eq/>
		<m:ci type="vector">x</m:ci>
		<m:apply>
		  <m:vector>
		    <m:cn>5</m:cn>
		    <m:cn>3</m:cn>
		  </m:vector>
		</m:apply>
		<m:apply>
		  <m:plus/>
		  <m:apply>
		    <m:times/>
		    <m:cn>4</m:cn>
		    <m:apply>
		      <m:vector>
			<m:cn>1</m:cn>
			<m:cn>1</m:cn>
		      </m:vector>
		    </m:apply>
		  </m:apply>
		  <m:apply>
		    <m:vector>
		      <m:cn>1</m:cn>
		      <m:cn>-1</m:cn>
		    </m:vector>
		  </m:apply>
		</m:apply>
	      </m:apply>
	    </m:math>

	    <m:math display="block">
	      <m:apply>
		<m:eq/>
		<m:apply>
		  <m:times/>
		  <m:ci>A</m:ci>
		  <m:ci type="vector">x</m:ci>
		</m:apply>
		<m:apply>
		  <m:times/>
		  <m:ci>A</m:ci>
		  <m:apply>
		    <m:plus/>
		    <m:apply>
		      <m:times/>
		      <m:cn>4</m:cn>
		      <m:ci type="vector">
			<m:msub>
			  <m:mi>v</m:mi>
			  <m:mn>1</m:mn>
			</m:msub>
		      </m:ci>
		    </m:apply>
		    <m:ci type="vector">
		      <m:msub>
			<m:mi>v</m:mi>
			<m:mn>2</m:mn>
		      </m:msub>
		    </m:ci>
		  </m:apply>
		</m:apply>
		<m:apply>
		  <m:times/>
		  <m:ci><m:msub>
		    <m:mi>λ</m:mi>
		    <m:mi>i</m:mi>
		  </m:msub></m:ci>
		  <m:apply>
		    <m:plus/>
		    <m:apply>
		      <m:times/>
		      <m:cn>4</m:cn>
		      <m:ci type="vector">
			<m:msub>
			  <m:mi>v</m:mi>
			  <m:mn>1</m:mn>
			</m:msub>
		      </m:ci>
		    </m:apply>
		    <m:ci type="vector">
		      <m:msub>
			<m:mi>v</m:mi>
			<m:mn>2</m:mn>
		      </m:msub>
		    </m:ci>
		  </m:apply>
		</m:apply>
	      </m:apply>
	    </m:math>

	    Therefore, we have

	    <m:math display="block">
	      <m:apply>
		<m:eq/>
		<m:apply>
		  <m:times/>
		  <m:ci>A</m:ci>
		  <m:ci type="vector">x</m:ci>
		</m:apply>
		<m:apply>
		  <m:plus/>
		  <m:apply>
		    <m:times/>
		    <m:cn>4</m:cn>
		    <m:cn>2</m:cn>
		    <m:apply>
		      <m:vector>
			<m:cn>1</m:cn>
			<m:cn>1</m:cn>
		      </m:vector>
		    </m:apply>
		  </m:apply>
		  <m:apply>
		    <m:times/>
		    <m:cn>4</m:cn>
		    <m:apply>
		      <m:vector>
			<m:cn>1</m:cn>
			<m:cn>-1</m:cn>
		      </m:vector>
		    </m:apply>
		  </m:apply>
		</m:apply>
		<m:apply>
		  <m:vector>
		    <m:cn>12</m:cn>
		    <m:cn>4</m:cn>
		  </m:vector>
		</m:apply>
	      </m:apply>
	    </m:math>
	    
	    Notice that this method using eigenvectors required
	    <emphasis>no</emphasis> matrix multiplication.  This may
	    have seemed more complicated here, but just imagine
	    <m:math><m:ci>A</m:ci></m:math> being really big, or even
	    just a few dimensions larger!
	  </para>
	</solution>
      </exercise><para id="element-369"><media type="application/x-labviewrpvi80" src="LinearAlgebraCalc3.llb">
		<param name="lvfppviname" value="Linear Algebra Calculator.vi"/>
		<param name="width" value="625"/>
		<param name="height" value="420"/>
	</media></para>



    </section>

  </content>
</document>
