FCLA Properties of Eigenvalues and Eigenvectors

Section PEE Properties of Eigenvalues and Eigenvectors

The previous section introduced eigenvalues and eigenvectors, and concentrated on their existence and determination. This section will be more about theorems, and the various properties eigenvalues and eigenvectors enjoy. Like a good 4×100 meter relay, we will lead-off with one of our better theorems and save the very best for the anchor leg.

Subsection BPE Basic Properties of Eigenvalues

Theorem EDELI. Eigenvectors with Distinct Eigenvalues are Linearly Independent.

Suppose that \(A\) is a square matrix and \(S=\set{\vectorlist{x}{p}}\) is a set of eigenvectors with eigenvalues \(\scalarlist{\lambda}{p}\) such that \(\lambda_i\neq\lambda_j\) whenever \(i\neq j\text{.}\) Then \(S\) is a linearly independent set.

Proof.

We will establish by induction (Proof Technique I) that any set of \(k\) eigenvectors of \(A\) with distinct eigenvalues \(\scalarlist{\lambda}{k}\) is a linearly independent set. Suppose \(A\) has size \(n\text{.}\)

Base Case.

When \(k=1\text{,}\) \(\set{\vect{x}_1}\) is a set with a single nonzero vector and thus is linearly independent.

Induction Step.

Begin with a relation of linear dependence on the set \(\set{\vectorlist{x}{k}}\)

\begin{equation*} \lincombo{a}{x}{k} = \zerovector\text{.} \end{equation*}

Then

\begin{align*} \zerovector &= \left(A-\lambda_k I_n\right)\zerovector&&\knowl{./knowl/xref/theorem-MMZM.html}{\text{Theorem MMZM}}\\ &= \left(A-\lambda_k I_n\right)\left(\lincombo{a}{x}{k}\right)&&\knowl{./knowl/xref/definition-RLD.html}{\text{Definition RLD}}\\ &=\left(A-\lambda_k I_n\right)a_1\vect{x}_1+ \cdots+ \left(A-\lambda_k I_n\right)a_k\vect{x}_k&& \knowl{./knowl/xref/theorem-MMDAA.html}{\text{Theorem MMDAA}}\\ &=a_1\left(A-\lambda_k I_n\right)\vect{x}_1+ \cdots+ a_k\left(A-\lambda_k I_n\right)\vect{x}_k&& \knowl{./knowl/xref/theorem-MMSMM.html}{\text{Theorem MMSMM}}\\ &=a_1\left(A\vect{x}_1-\lambda_k I_n\vect{x}_1\right)+ \cdots+ a_k\left(A\vect{x}_k-\lambda_k I_n\vect{x}_k\right)&& \knowl{./knowl/xref/theorem-MMDAA.html}{\text{Theorem MMDAA}}\\ &=a_1\left(A\vect{x}_1-\lambda_k\vect{x}_1\right)+ \cdots+ a_k\left(A\vect{x}_k-\lambda_k\vect{x}_k\right)&& \knowl{./knowl/xref/theorem-MMIM.html}{\text{Theorem MMIM}}\\ &=a_1\left(\lambda_1\vect{x}_1-\lambda_k\vect{x}_1\right)+ \cdots+ a_k\left(\lambda_k\vect{x}_k-\lambda_k\vect{x}_k\right)&& \knowl{./knowl/xref/definition-EEM.html}{\text{Definition EEM}}\\ &=a_1\left(\lambda_1-\lambda_k\right)\vect{x}_1+ \cdots+ a_k\left(\lambda_k-\lambda_k\right)\vect{x}_k&& \knowl{./knowl/xref/property-DSA.html}{\text{Property DSA}}\\ &=a_1\left(\lambda_1-\lambda_k\right)\vect{x}_1+ \cdots+ a_{k-1}\left(\lambda_{k-1}-\lambda_k\right)\vect{x}_{k-1}+ a_k\left(0\right)\vect{x}_k&& \knowl{./knowl/xref/property-AICN.html}{\text{Property AICN}}\\ &=a_1\left(\lambda_1-\lambda_k\right)\vect{x}_1+ \cdots+ a_{k-1}\left(\lambda_{k-1}-\lambda_k\right)\vect{x}_{k-1}+ \zerovector&& \knowl{./knowl/xref/theorem-ZSSM.html}{\text{Theorem ZSSM}}\\ &=a_1\left(\lambda_1-\lambda_k\right)\vect{x}_1+ \cdots+ a_{k-1}\left(\lambda_{k-1}-\lambda_k\right)\vect{x}_{k-1}&& \knowl{./knowl/xref/property-Z.html}{\text{Property Z}} \end{align*}

This equation is a relation of linear dependence on the set \(\set{\vectorlist{x}{k-1}}\text{,}\) which is a linearly independent set by the induction hypothesis. So the scalars must all be zero by Definition LI. That is, \(a_i\left(\lambda_i-\lambda_k\right)=0\) for \(1\leq i\leq k-1\text{.}\) However, we have the hypothesis that the eigenvalues are distinct, so \(\lambda_i-\lambda_k\neq 0\) for \(1\leq i\leq k-1\text{.}\) So Theorem ZPZF implies \(a_i=0\) for \(1\leq i\leq k-1\text{.}\)

This reduces the original relation of linear dependence on \(\set{\vectorlist{x}{k}}\) to the simpler equation \(a_k\vect{x}_k=\zerovector\text{.}\) By Theorem SMEZV we conclude that \(a_k=0\) or \(\vect{x}_k=\zerovector\text{.}\) Eigenvectors are never the zero vector (Definition EEM), so \(a_k=0\text{.}\) Now all of the scalars \(a_i\text{,}\) \(1\leq i\leq k\) are zero, and so the only relation of linear dependence on the set \(\set{\vectorlist{x}{k}}\) is trivial. So by Definition LI, the set \(\set{\vectorlist{x}{k}}\) is linearly independent.

The next theorem gives us a convenient upper limit on the number of eigenvalues.

Theorem MNEM. Maximum Number of Eigenvalues of a Matrix.

Suppose that \(A\) is a square matrix of size \(n\text{.}\) Then \(A\) cannot have more than \(n\) distinct eigenvalues.

Proof.

We argue by contradiction (Proof Technique CD). Assume that \(A\) has \(n+1\) or more distinct eigenvalues. Then there is a set of \(n+1\) or more eigenvectors of \(A\text{,}\) with distinct eigenvalues. This is a set of \(n+1\) or more vectors from \(\complex{n}\) that will be linearly independent by Theorem EDELI. But this contradicts Theorem MVSLD, so our assumption is false, and there are no more than \(n\) distinct eigenvalues.

Notice that once we have found \(n\) distinct eigenvalues for a matrix of size \(n\text{,}\) then we know there are no more eigenvalues. Example DEMS5 is an example, and the upcoming Theorem DED also considers this situation.

There is a simple connection between the eigenvalues of a matrix and whether or not the matrix is nonsingular.

Theorem SMZE. Singular Matrices have Zero Eigenvalues.

Suppose \(A\) is a square matrix. Then \(A\) is singular if and only if \(0\) is an eigenvalue of \(A\text{.}\)

Proof.

We have the following equivalences:

\begin{align*} A\text{ is singular} &\iff A-0I_n\text{ is singular}&&\knowl{./knowl/xref/property-ZM.html}{\text{Property ZM}}\\ &\iff 0\text{ is an eigenvalue of A}&&\knowl{./knowl/xref/theorem-ESM.html}{\text{Theorem ESM}}\text{.} \end{align*}

With an equivalence about singular matrices we can update our list of equivalences about nonsingular matrices.

Theorem NME8. Nonsingular Matrix Equivalences, Round 8.

Suppose that \(A\) is a square matrix of size \(n\text{.}\) The following are equivalent.

\(A\) is nonsingular.
\(A\) row-reduces to the identity matrix.
The null space of \(A\) contains only the zero vector, \(\nsp{A}=\set{\zerovector}\text{.}\)
The linear system \(\linearsystem{A}{\vect{b}}\) has a unique solution for every possible choice of \(\vect{b}\text{.}\)
The columns of \(A\) are a linearly independent set.
\(A\) is invertible.
The column space of \(A\) is \(\complex{n}\text{,}\) \(\csp{A}=\complex{n}\text{.}\)
The columns of \(A\) are a basis for \(\complex{n}\text{.}\)
The rank of \(A\) is \(n\text{,}\) \(\rank{A}=n\text{.}\)
The nullity of \(A\) is zero, \(\nullity{A}=0\text{.}\)
The determinant of \(A\) is nonzero, \(\detname{A}\neq 0\text{.}\)
\(\lambda=0\) is not an eigenvalue of \(A\text{.}\)

Proof.

The equivalence of the first and last statements is Theorem SMZE, reformulated by negating each statement in the equivalence. So we are able to improve on Theorem NME7 with this addition.

Sage NME8. Nonsingular Matrix Equivalences, Round 8.

Zero eigenvalues are another marker of singular matrices. We illustrate with two matrices, the first nonsingular, the second not.

Certain changes to a matrix change its eigenvalues in a predictable way.

Theorem ESMM. Eigenvalues of a Scalar Multiple of a Matrix.

Suppose \(A\) is a square matrix and \(\lambda\) is an eigenvalue of \(A\text{.}\) Then \(\alpha\lambda\) is an eigenvalue of \(\alpha A\text{.}\)

Proof.

Let \(\vect{x}\neq\zerovector\) be one eigenvector of \(A\) for \(\lambda\text{.}\) Then

\begin{align*} \left(\alpha A\right)\vect{x}&=\alpha\left(A\vect{x}\right)&& \knowl{./knowl/xref/theorem-MMSMM.html}{\text{Theorem MMSMM}}\\ &=\alpha\left(\lambda\vect{x}\right)&& \knowl{./knowl/xref/definition-EEM.html}{\text{Definition EEM}}\\ &=\left(\alpha\lambda\right)\vect{x}&& \knowl{./knowl/xref/property-SMAC.html}{\text{Property SMAC}}\text{.} \end{align*}

So \(\vect{x}\neq\zerovector\) is an eigenvector of \(\alpha A\) for the eigenvalue \(\alpha\lambda\text{.}\)

Unfortunately, there are not parallel theorems about the sum or product of arbitrary matrices. But we can prove a similar result for powers of a matrix.

Theorem EOMP. Eigenvalues Of Matrix Powers.

Suppose \(A\) is a square matrix, \(\lambda\) is an eigenvalue of \(A\text{,}\) and \(s\geq 0\) is an integer. Then \(\lambda^s\) is an eigenvalue of \(A^s\text{.}\)

Proof.

Let \(\vect{x}\neq\zerovector\) be one eigenvector of \(A\) for \(\lambda\text{.}\) Then we proceed by induction on \(s\) (Proof Technique I). First, for \(s=0\text{,}\) employing Theorem MMIM and Property OC to establish the base case,

\begin{gather*} A^s\vect{x}=A^0\vect{x}=I_n\vect{x}=\vect{x}=1\vect{x}=\lambda^0\vect{x}=\lambda^s\vect{x}\text{.} \end{gather*}

So \(\lambda^s\) is an eigenvalue of \(A^s\) in this special case. For the induction step, we assume the theorem is true for \(s\text{,}\) and find

\begin{align*} A^{s+1}\vect{x}&=A^sA\vect{x}\\ &=A^s\left(\lambda\vect{x}\right)&& \knowl{./knowl/xref/definition-EEM.html}{\text{Definition EEM}}\\ &=\lambda\left(A^s\vect{x}\right)&& \knowl{./knowl/xref/theorem-MMSMM.html}{\text{Theorem MMSMM}}\\ &=\lambda\left(\lambda^s\vect{x}\right)&&\text{Induction Hypothesis}\\ &=\left(\lambda\lambda^s\right)\vect{x}&& \knowl{./knowl/xref/property-SMAC.html}{\text{Property SMAC}}\\ &=\lambda^{s+1}\vect{x}\text{.} \end{align*}

So \(\vect{x}\neq\zerovector\) is an eigenvector of \(A^{s+1}\) for \(\lambda^{s+1}\text{,}\) and induction tells us the theorem is true for all \(s\geq 0\text{.}\)

While we cannot prove that the sum of two arbitrary matrices behaves in any reasonable way with regard to eigenvalues, we can work with the sum of dissimilar powers of the same matrix. We have already seen that eigenvalues arise as roots of polynomials (Theorem EMRCP). A theme of this chapter is relationships between eigenvalues and polynomials, and our next theorem is the next component of this theme.

Theorem EPM. Eigenvalues of the Polynomial of a Matrix.

Suppose \(A\) is a square matrix and \(\lambda\) is an eigenvalue of \(A\text{.}\) Let \(q(x)\) be a polynomial in the variable \(x\text{.}\) Then \(q(\lambda)\) is an eigenvalue of the matrix \(q(A)\text{.}\)

Proof.

Let \(\vect{x}\neq\zerovector\) be one eigenvector of \(A\) for \(\lambda\text{,}\) and write \(q(x)=a_0+a_1x+a_2x^2+\cdots+a_mx^m\text{.}\) Then

\begin{align*} q(A)\vect{x}&=\left(a_0A^0+a_1A^1+a_2A^2+\cdots+a_mA^m\right)\vect{x}\\ &=(a_0A^0)\vect{x}+(a_1A^1)\vect{x}+(a_2A^2)\vect{x}+\cdots+(a_mA^m)\vect{x}&& \knowl{./knowl/xref/theorem-MMDAA.html}{\text{Theorem MMDAA}}\\ &=a_0(A^0\vect{x})+a_1(A^1\vect{x})+a_2(A^2\vect{x})+\cdots+a_m(A^m\vect{x})&& \knowl{./knowl/xref/theorem-MMSMM.html}{\text{Theorem MMSMM}}\\ &=a_0(\lambda^0\vect{x})+a_1(\lambda^1\vect{x})+a_2(\lambda^2\vect{x})+\cdots+a_m(\lambda^m\vect{x})&& \knowl{./knowl/xref/theorem-EOMP.html}{\text{Theorem EOMP}}\\ &=(a_0\lambda^0)\vect{x}+(a_1\lambda^1)\vect{x}+(a_2\lambda^2)\vect{x}+\cdots+(a_m\lambda^m)\vect{x}&& \knowl{./knowl/xref/property-SMAC.html}{\text{Property SMAC}}\\ &=\left(a_0\lambda^0+a_1\lambda^1+a_2\lambda^2+\cdots+a_m\lambda^m\right)\vect{x}&& \knowl{./knowl/xref/property-DSAC.html}{\text{Property DSAC}}\\ &=q(\lambda)\vect{x}\text{.} \end{align*}

So \(\vect{x}\neq 0\) is an eigenvector of \(q(A)\) for the eigenvalue \(q(\lambda)\text{.}\)

Inverses and transposes also behave predictably with regard to their eigenvalues.

Theorem EIM. Eigenvalues of the Inverse of a Matrix.

Suppose \(A\) is a square nonsingular matrix and \(\lambda\) is an eigenvalue of \(A\text{.}\) Then \(\lambda^{-1}\) is an eigenvalue of the matrix \(\inverse{A}\text{.}\)

Proof.

Notice that since \(A\) is assumed nonsingular, \(\inverse{A}\) exists by Theorem NI, but more importantly, \(\lambda^{-1}=1/\lambda\) does not involve division by zero since Theorem SMZE prohibits this possibility.

Let \(\vect{x}\neq\zerovector\) be one eigenvector of \(A\) for \(\lambda\text{.}\) Suppose \(A\) has size \(n\text{.}\) Then

\begin{align*} \inverse{A}\vect{x}&=\inverse{A}(1\vect{x})&& \knowl{./knowl/xref/property-OC.html}{\text{Property OC}}\\ &=\inverse{A}(\frac{1}{\lambda}\lambda\vect{x})&& \knowl{./knowl/xref/property-MICN.html}{\text{Property MICN}}\\ &=\frac{1}{\lambda}\inverse{A}(\lambda\vect{x})&& \knowl{./knowl/xref/theorem-MMSMM.html}{\text{Theorem MMSMM}}\\ &=\frac{1}{\lambda}\inverse{A}(A\vect{x})&& \knowl{./knowl/xref/definition-EEM.html}{\text{Definition EEM}}\\ &=\frac{1}{\lambda}(\inverse{A}A)\vect{x}&& \knowl{./knowl/xref/theorem-MMA.html}{\text{Theorem MMA}}\\ &=\frac{1}{\lambda}I_n\vect{x}&& \knowl{./knowl/xref/definition-MI.html}{\text{Definition MI}}\\ &=\frac{1}{\lambda}\vect{x}&& \knowl{./knowl/xref/theorem-MMIM.html}{\text{Theorem MMIM}}\text{.} \end{align*}

So \(\vect{x}\neq\zerovector\) is an eigenvector of \(\inverse{A}\) for the eigenvalue \(\frac{1}{\lambda}\text{.}\)

The proofs of the theorems above have a similar style to them. They all begin by grabbing an eigenvalue-eigenvector pair and adjusting it in some way to reach the desired conclusion. You should add this to your toolkit as a general approach to proving theorems about eigenvalues.

So far we have been able to reserve the characteristic polynomial for strictly computational purposes. However, sometimes a theorem about eigenvalues can be proved easily by employing the characteristic polynomial (rather than using an eigenvalue-eigenvector pair). The next theorem is an example of this.

Theorem ETM. Eigenvalues of the Transpose of a Matrix.

Suppose \(A\) is a square matrix and \(\lambda\) is an eigenvalue of \(A\text{.}\) Then \(\lambda\) is an eigenvalue of the matrix \(\transpose{A}\text{.}\)

Proof.

Suppose \(A\) has size \(n\text{.}\) Then

\begin{align*} \charpoly{A}{x} &=\detname{A-xI_n}&& \knowl{./knowl/xref/definition-CP.html}{\text{Definition CP}}\\ &=\detname{\transpose{\left(A-xI_n\right)}}&& \knowl{./knowl/xref/theorem-DT.html}{\text{Theorem DT}}\\ &=\detname{\transpose{A}-\transpose{\left(xI_n\right)}}&& \knowl{./knowl/xref/theorem-TMA.html}{\text{Theorem TMA}}\\ &=\detname{\transpose{A}-x\transpose{I_n}}&& \knowl{./knowl/xref/theorem-TMSM.html}{\text{Theorem TMSM}}\\ &=\detname{\transpose{A}-xI_n}&& \knowl{./knowl/xref/definition-IM.html}{\text{Definition IM}}\\ &=\charpoly{\transpose{A}}{x}&& \knowl{./knowl/xref/definition-CP.html}{\text{Definition CP}}\text{.} \end{align*}

So \(A\) and \(\transpose{A}\) have the same characteristic polynomial, and by Theorem EMRCP, their eigenvalues are identical and have equal algebraic multiplicities. Notice that what we have proved here is a bit stronger than the stated conclusion in the theorem.

If a matrix has only real entries, then eigenvalues will arise as the roots of a polynomial with real coefficients. Complex numbers could result as roots of this polynomial, but they are roots of quadratic factors with real coefficients, and as such, come in conjugate pairs. The next theorem proves this, and a bit more, without ever mentioning a polynomial.

Theorem ERMCP. Eigenvalues of Real Matrices come in Conjugate Pairs.

Suppose \(A\) is a square matrix with real entries and \(\vect{x}\) is an eigenvector of \(A\) for the eigenvalue \(\lambda\text{.}\) Then \(\conjugate{\vect{x}}\) is an eigenvector of \(A\) for the eigenvalue \(\conjugate{\lambda}\text{.}\)

Proof.

We have

\begin{align*} A\conjugate{\vect{x}}&=\conjugate{A}\conjugate{\vect{x}}&& A\text{ has real entries}\\ &=\conjugate{A\vect{x}}&& \knowl{./knowl/xref/theorem-MMCC.html}{\text{Theorem MMCC}}\\ &=\conjugate{\lambda\vect{x}}&& \knowl{./knowl/xref/definition-EEM.html}{\text{Definition EEM}}\\ &=\conjugate{\lambda}\conjugate{\vect{x}}&& \knowl{./knowl/xref/theorem-CRSM.html}{\text{Theorem CRSM}}\text{.} \end{align*}

So \(\conjugate{\vect{x}}\) is an eigenvector of \(A\) for the eigenvalue \(\conjugate{\lambda}\text{.}\)

This phenomenon is amply illustrated in Example CEMS6, where the four complex eigenvalues come in two pairs, and the two basis vectors of the eigenspaces are complex conjugates of each other. Theorem ERMCP can be a time-saver for computing eigenvalues and eigenvectors of real matrices with complex eigenvalues, since the conjugate eigenvalue and eigenspace can be inferred from the theorem rather than computed.

Subsection ME Multiplicities of Eigenvalues

A polynomial of degree \(n\) will have exactly \(n\) roots. From this fact about polynomial equations we can say more about the algebraic multiplicities of eigenvalues.

Theorem DCP. Degree of the Characteristic Polynomial.

Suppose that \(A\) is a square matrix of size \(n\text{.}\) Then the characteristic polynomial of \(A\text{,}\) \(\charpoly{A}{x}\text{,}\) has degree \(n\text{.}\)

Proof.

We will prove a more general result by induction (Proof Technique I). Then the theorem will be true as a special case. We will carefully state this result as a proposition indexed by \(m\text{,}\) \(m\geq 1\text{.}\)

\(P(m)\text{:}\) Suppose that \(A\) is an \(m\times m\) matrix whose entries are complex numbers or linear polynomials in the variable \(x\) of the form \(c-x\text{,}\) where \(c\) is a complex number. Suppose further that there are exactly \(k\) entries that contain \(x\) and that no row or column contains more than one such entry. Then, when \(k=m\text{,}\) \(\detname{A}\) is a polynomial in \(x\) of degree \(m\text{,}\) with leading coefficient \(\pm 1\text{,}\) and when \(k\lt m\text{,}\) \(\detname{A}\) is a polynomial in \(x\) of degree \(k\) or less.

Base Case: Suppose \(A\) is a \(1\times 1\) matrix. Then its determinant is equal to the lone entry (Definition DM). When \(k=m=1\text{,}\) the entry is of the form \(c-x\text{,}\) a polynomial in \(x\) of degree \(m=1\) with leading coefficient \(-1\text{.}\) When \(k\lt m\text{,}\) then \(k=0\) and the entry is simply a complex number, a polynomial of degree \(0\leq k\text{.}\) So \(P(1)\) is true.

Induction Step: Assume \(P(m)\) is true, and that \(A\) is an \((m+1)\times(m+1)\) matrix with \(k\) entries of the form \(c-x\text{.}\) There are two cases to consider.

Suppose \(k=m+1\text{.}\) Then every row and every column will contain an entry of the form \(c-x\text{.}\) Suppose that for the first row, this entry is in column \(t\text{.}\) Compute the determinant of \(A\) by an expansion about this first row (Definition DM). The term associated with entry \(t\) of this row will be of the form

\begin{equation*} (c-x)(-1)^{1+t}\detname{\submatrix{A}{1}{t}}\text{.} \end{equation*}

The submatrix \(\submatrix{A}{1}{t}\) is an \(m\times m\) matrix with \(k=m\) terms of the form \(c-x\text{,}\) no more than one per row or column. By the induction hypothesis, \(\detname{\submatrix{A}{1}{t}}\) will be a polynomial in \(x\) of degree \(m\) with coefficient \(\pm 1\text{.}\) So this entire term is then a polynomial of degree \(m+1\) with leading coefficient \(\pm 1\text{.}\)

The remaining terms (which constitute the sum that is the determinant of \(A\)) are products of complex numbers from the first row with cofactors built from submatrices that lack the first row of \(A\) and lack some column of \(A\text{,}\) other than column \(t\text{.}\) As such, these submatrices are \(m\times m\) matrices with \(k=m-1\lt m\) entries of the form \(c-x\text{,}\) no more than one per row or column. Applying the induction hypothesis, we see that these terms are polynomials in \(x\) of degree \(m-1\) or less. Adding the single term from the entry in column \(t\) with all these others, we see that \(\detname{A}\) is a polynomial in \(x\) of degree \(m+1\) and leading coefficient \(\pm 1\text{.}\)

The second case occurs when \(k\lt m+1\text{.}\) Now there is a row of \(A\) that does not contain an entry of the form \(c-x\text{.}\) We consider the determinant of \(A\) by expanding about this row (Theorem DER), whose entries are all complex numbers. The cofactors employed are built from submatrices that are \(m\times m\) matrices with either \(k\) or \(k-1\) entries of the form \(c-x\text{,}\) no more than one per row or column. In either case, \(k\leq m\text{,}\) and we can apply the induction hypothesis to see that the determinants computed for the cofactors are all polynomials of degree \(k\) or less. Summing these contributions to the determinant of \(A\) yields a polynomial in \(x\) of degree \(k\) or less, as desired.

Definition CP tells us that the characteristic polynomial of an \(n\times n\) matrix is the determinant of a matrix having exactly \(n\) entries of the form \(c-x\text{,}\) no more than one per row or column. As such we can apply \(P(n)\) to see that the characteristic polynomial has degree \(n\text{.}\)

Theorem NEM. Number of Eigenvalues of a Matrix.

Suppose that \(\scalarlist{\lambda}{k}\) are the distinct eigenvalues of a square matrix \(A\) of size \(n\text{.}\) Then

\begin{equation*} \sum_{i=1}^{k}\algmult{A}{\lambda_i}=n\text{.} \end{equation*}

Proof.

By the definition of the algebraic multiplicity (Definition AME), we can factor the characteristic polynomial as

\begin{equation*} \charpoly{A}{x}=c(x-\lambda_1)^{\algmult{A}{\lambda_1}}(x-\lambda_2)^{\algmult{A}{\lambda_2}}(x-\lambda_3)^{\algmult{A}{\lambda_3}}\cdots(x-\lambda_k)^{\algmult{A}{\lambda_k}} \end{equation*}

where \(c\) is a nonzero constant. (We could prove that \(c=(-1)^{n}\text{,}\) but we do not need that specificity right now. See Exercise PEE.T30) The left-hand side is a polynomial of degree \(n\) by Theorem DCP and the right-hand side is a polynomial of degree \(\sum_{i=1}^{k}\algmult{A}{\lambda_i}\text{.}\) So the equality of the polynomials’ degrees gives the equality \(\sum_{i=1}^{k}\algmult{A}{\lambda_i}=n\text{.}\)

Theorem ME. Multiplicities of an Eigenvalue.

Suppose that \(A\) is a square matrix of size \(n\) and \(\lambda\) is an eigenvalue. Then

\begin{equation*} 1\leq\geomult{A}{\lambda}\leq\algmult{A}{\lambda}\leq n\text{.} \end{equation*}

Proof.

Since \(\lambda\) is an eigenvalue of \(A\text{,}\) there is an eigenvector of \(A\) for \(\lambda\text{,}\) \(\vect{x}\text{.}\) Then \(\vect{x}\in\eigenspace{A}{\lambda}\text{,}\) so \(\geomult{A}{\lambda}\geq 1\text{,}\) since we can extend \(\set{\vect{x}}\) into a basis of \(\eigenspace{A}{\lambda}\) (Theorem ELIS).

To show \(\geomult{A}{\lambda}\leq\algmult{A}{\lambda}\) is the most involved portion of this proof. To this end, let \(g=\geomult{A}{\lambda}\) and let \(\vectorlist{x}{g}\) be a basis for the eigenspace of \(\lambda\text{,}\) \(\eigenspace{A}{\lambda}\text{.}\) Construct another \(n-g\) vectors, \(\vectorlist{y}{n-g}\text{,}\) so that

\begin{equation*} \set{\vectorlist{x}{g},\,\vectorlist{y}{n-g}} \end{equation*}

is a basis of \(\complex{n}\text{.}\) This can be done by repeated applications of Theorem ELIS.

Finally, define a matrix \(S\) by

\begin{equation*} S=[\vect{x}_1|\vect{x}_2|\vect{x}_3|\ldots|\vect{x}_g|\vect{y}_1|\vect{y}_2|\vect{y}_3|\ldots|\vect{y}_{n-g}]=[\vect{x}_1|\vect{x}_2|\vect{x}_3|\ldots|\vect{x}_g|R] \end{equation*}

where \(R\) is an \(n\times(n-g)\) matrix whose columns are \(\vectorlist{y}{n-g}\text{.}\) The columns of \(S\) are linearly independent by design, so \(S\) is nonsingular (Theorem NMLIC) and therefore invertible (Theorem NI).

Then

\begin{align*} \matrixcolumns{e}{n}&=I_n\\ &=\inverse{S}S\\ &=\inverse{S}[\vect{x}_1|\vect{x}_2|\vect{x}_3|\ldots|\vect{x}_g|R]\\ &=[\inverse{S}\vect{x}_1|\inverse{S}\vect{x}_2|\inverse{S}\vect{x}_3|\ldots|\inverse{S}\vect{x}_g|\inverse{S}R]\text{.} \end{align*}

\begin{equation*} \inverse{S}\vect{x}_i=\vect{e}_i\quad 1\leq i\leq g\tag{*}\text{.} \end{equation*}

Preparations in place, we compute the characteristic polynomial of \(A\text{,}\)

\begin{align*} \charpoly{A}{x}&=\detname{A-xI_n}&& \knowl{./knowl/xref/definition-CP.html}{\text{Definition CP}}\\ &=1\detname{A-xI_n}&& \knowl{./knowl/xref/property-OCN.html}{\text{Property OCN}}\\ &=\detname{I_n}\detname{A-xI_n}&& \knowl{./knowl/xref/definition-DM.html}{\text{Definition DM}}\\ &=\detname{\inverse{S}S}\detname{A-xI_n}&& \knowl{./knowl/xref/definition-MI.html}{\text{Definition MI}}\\ &=\detname{\inverse{S}}\detname{S}\detname{A-xI_n}&& \knowl{./knowl/xref/theorem-DRMM.html}{\text{Theorem DRMM}}\\ &=\detname{\inverse{S}}\detname{A-xI_n}\detname{S}&& \knowl{./knowl/xref/property-CMCN.html}{\text{Property CMCN}}\\ &=\detname{\inverse{S}\left(A-xI_n\right)S}&& \knowl{./knowl/xref/theorem-DRMM.html}{\text{Theorem DRMM}}\\ &=\detname{\inverse{S}AS-\inverse{S}xI_nS}&& \knowl{./knowl/xref/theorem-MMDAA.html}{\text{Theorem MMDAA}}\\ &=\detname{\inverse{S}AS-x\inverse{S}I_nS}&& \knowl{./knowl/xref/theorem-MMSMM.html}{\text{Theorem MMSMM}}\\ &=\detname{\inverse{S}AS-x\inverse{S}S}&& \knowl{./knowl/xref/theorem-MMIM.html}{\text{Theorem MMIM}}\\ &=\detname{\inverse{S}AS-xI_n}&& \knowl{./knowl/xref/definition-MI.html}{\text{Definition MI}}\\ &=\charpoly{\inverse{S}AS}{x}&& \knowl{./knowl/xref/definition-CP.html}{\text{Definition CP}}\text{.} \end{align*}

What can we learn then about the matrix \(\inverse{S}AS\text{?}\)

\begin{align*} \inverse{S}AS&=\inverse{S}A[\vect{x}_1|\vect{x}_2|\vect{x}_3|\ldots|\vect{x}_g|R]\\ &=\inverse{S}[A\vect{x}_1|A\vect{x}_2|A\vect{x}_3|\ldots|A\vect{x}_g|AR]&& \knowl{./knowl/xref/definition-MM.html}{\text{Definition MM}}\\ &=\inverse{S}[\lambda\vect{x}_1|\lambda\vect{x}_2|\lambda\vect{x}_3|\ldots|\lambda\vect{x}_g|AR]&& \knowl{./knowl/xref/definition-EEM.html}{\text{Definition EEM}}\\ &=[\inverse{S}\lambda\vect{x}_1|\inverse{S}\lambda\vect{x}_2|\inverse{S}\lambda\vect{x}_3|\ldots|\inverse{S}\lambda\vect{x}_g|\inverse{S}AR]&& \knowl{./knowl/xref/definition-MM.html}{\text{Definition MM}}\\ &=[\lambda\inverse{S}\vect{x}_1|\lambda\inverse{S}\vect{x}_2|\lambda\inverse{S}\vect{x}_3|\ldots|\lambda\inverse{S}\vect{x}_g|\inverse{S}AR]&& \knowl{./knowl/xref/theorem-MMSMM.html}{\text{Theorem MMSMM}}\\ &=[\lambda\vect{e}_1|\lambda\vect{e}_2|\lambda\vect{e}_3|\ldots|\lambda\vect{e}_g|\inverse{S}AR]&& \inverse{S}S=I_n,\text{ (*) above} \end{align*}

Now imagine computing the characteristic polynomial of \(A\) by computing the characteristic polynomial of \(\inverse{S}AS\) using the form just obtained. The first \(g\) columns of \(\inverse{S}AS\) are all zero, save for a \(\lambda\) on the diagonal. So if we compute the determinant by expanding about the first column, successively, we will get successive factors of \((\lambda-x)\text{.}\) More precisely, let \(T\) be the square matrix of size \(n-g\) that is formed from the last \(n-g\) rows and last \(n-g\) columns of \(\inverse{S}AR\text{.}\) Then

\begin{equation*} \charpoly{A}{x}=\charpoly{\inverse{S}AS}{x}=(\lambda-x)^g\charpoly{T}{x}\text{.} \end{equation*}

This says that \((x-\lambda)\) is a factor of the characteristic polynomial at least \(g\) times, so the algebraic multiplicity of \(\lambda\) as an eigenvalue of \(A\) is greater than or equal to \(g\) (Definition AME). In other words,

\begin{equation*} \geomult{A}{\lambda}=g\leq\algmult{A}{\lambda} \end{equation*}

as desired.

Theorem NEM says that the sum of the algebraic multiplicities for all the eigenvalues of \(A\) is equal to \(n\text{.}\) Since the algebraic multiplicity is a positive quantity, no single algebraic multiplicity can exceed \(n\) without the sum of all of the algebraic multiplicities doing the same.

Subsection EHM Eigenvalues of Hermitian Matrices

Recall that a matrix is Hermitian (or self-adjoint) if \(A=\adjoint{A}\) (Definition HM). In the case where \(A\) is a matrix whose entries are all real numbers, being Hermitian is identical to being symmetric (Definition SYM). Keep this in mind as you read the next two theorems. Their hypotheses could be changed to “suppose \(A\) is a real symmetric matrix.”

Theorem HMRE. Hermitian Matrices have Real Eigenvalues.

Suppose that \(A\) is a Hermitian matrix and \(\lambda\) is an eigenvalue of \(A\text{.}\) Then \(\lambda\in{\mathbb R}\text{.}\)

Proof.

Let \(\vect{x}\neq\zerovector\) be one eigenvector of \(A\) for the eigenvalue \(\lambda\text{.}\) Then

\begin{align*} \left(\lambda - \conjugate{\lambda}\right)\innerproduct{\vect{x}}{\vect{x}} &=\lambda\innerproduct{\vect{x}}{\vect{x}} - \conjugate{\lambda}\innerproduct{\vect{x}}{\vect{x}}&& \knowl{./knowl/xref/property-DCN.html}{\text{Property DCN}}\\ &=\innerproduct{\vect{x}}{\lambda\vect{x}} - \innerproduct{\lambda\vect{x}}{\vect{x}}&& \knowl{./knowl/xref/theorem-IPSM.html}{\text{Theorem IPSM}}\\ &=\innerproduct{\vect{x}}{A\vect{x}} - \innerproduct{A\vect{x}}{\vect{x}}&& \knowl{./knowl/xref/definition-EEM.html}{\text{Definition EEM}}\\ &=\innerproduct{A\vect{x}}{\vect{x}} - \innerproduct{A\vect{x}}{\vect{x}}&& \knowl{./knowl/xref/theorem-HMIP.html}{\text{Theorem HMIP}}\\ &=0\text{.} \end{align*}

Since \(\vect{x}\neq\zerovector\text{,}\) by Theorem PIP we know \(\innerproduct{\vect{x}}{\vect{x}}\neq 0\text{.}\) Then by Theorem ZPZF, \(\lambda - \conjugate{\lambda}=0\text{,}\) and so \(\lambda = \conjugate{\lambda}\text{.}\) If a complex number is equal to its conjugate, then it has a complex part equal to zero, and therefore is a real number.

Notice the key step of this proof is the ability to pitch a Hermitian matrix from one side of the inner product to the other.

Look back and compare Example ESMS4 and Example CEMS6. In Example CEMS6 the matrix has only real entries, yet the eigenvalues come as roots of polynomial that are complex numbers. However, in Example ESMS4, the matrix has only real entries, but is also symmetric, and hence Hermitian. So by Theorem HMRE, we were guaranteed eigenvalues that are real numbers.

In many physical problems, a matrix of interest will be real and symmetric, or Hermitian. Then if the eigenvalues are to represent physical quantities of interest, Theorem HMRE guarantees that these values will not be complex numbers.

The eigenvectors of a Hermitian matrix also enjoy a pleasing property that we will exploit later.

Theorem HMOE. Hermitian Matrices have Orthogonal Eigenvectors.

Suppose that \(A\) is a Hermitian matrix and \(\vect{x}\) and \(\vect{y}\) are two eigenvectors of \(A\) for different eigenvalues. Then \(\vect{x}\) and \(\vect{y}\) are orthogonal vectors.

Proof.

Let \(\vect{x}\) be an eigenvector of \(A\) for \(\lambda\) and let \(\vect{y}\) be an eigenvector of \(A\) for a different eigenvalue \(\rho\text{.}\) So we have \(\lambda-\rho\neq 0\text{.}\) Then

\begin{align*} \left(\lambda-\rho\right)\innerproduct{\vect{x}}{\vect{y}} &=\lambda\innerproduct{\vect{x}}{\vect{y}}-\rho\innerproduct{\vect{x}}{\vect{y}}&& \knowl{./knowl/xref/property-DCN.html}{\text{Property DCN}}\\ &=\innerproduct{\conjugate{\lambda}\vect{x}}{\vect{y}}-\innerproduct{\vect{x}}{\rho\vect{y}}&& \knowl{./knowl/xref/theorem-IPSM.html}{\text{Theorem IPSM}}\\ &=\innerproduct{\lambda\vect{x}}{\vect{y}}-\innerproduct{\vect{x}}{\rho\vect{y}}&& \knowl{./knowl/xref/theorem-HMRE.html}{\text{Theorem HMRE}}\\ &=\innerproduct{A\vect{x}}{\vect{y}}-\innerproduct{\vect{x}}{A\vect{y}}&& \knowl{./knowl/xref/definition-EEM.html}{\text{Definition EEM}}\\ &=\innerproduct{A\vect{x}}{\vect{y}}-\innerproduct{A\vect{x}}{\vect{y}}&& \knowl{./knowl/xref/theorem-HMIP.html}{\text{Theorem HMIP}}\\ &=0\text{.} \end{align*}

Because \(\lambda\) and \(\rho\) are presumed to be different, \(\lambda-\rho\neq 0\text{,}\) and Theorem ZPZF implies that \(\innerproduct{\vect{x}}{\vect{y}}=0\text{.}\) In other words, \(\vect{x}\) and \(\vect{y}\) are orthogonal vectors according to Definition OV.

Notice again how the key step in this proof is the fundamental property of a Hermitian matrix (Theorem HMIP) — the ability to swap \(A\) across the two arguments of the inner product. Notice too, that we can always apply the Gram-Schmidt procedure (Theorem GSP) to any basis of any eigenspace, along with scaling the resulting the orthogonal vectors by their norm to arrive at an orthonormal basis of the eigenspace. For a Hermitian matrix, pairs of eigenvectors from different eigenspaces are also orthogonal. If we dumped all these basis vectors into one big set it would be an orthonormal set. We will build on these results and continue to see some more interesting properties in Section OD.

Reading Questions PEE Reading Questions

1. Eigenvalues of a Nonsingular Matrix.

How can you identify a nonsingular matrix just by looking at its eigenvalues?

2. Maximum Number of Distinct Eigenvalues.

How many different eigenvalues may a square matrix of size \(n\) have?

3. Eigenvalues of a Hermitian Matrix.

What is amazing about the eigenvalues of a Hermitian matrix and why is it amazing?

Exercises PEE Exercises

M20.

This exercise will show we can use a polynomial to convert one matrix into another, with predictable changes in its eigenvalues. In Example ESMS4 the \(4\times 4\) symmetric matrix

\begin{equation*} C= \begin{bmatrix} 1 & 0 & 1 & 1\\ 0 & 1 & 1 & 1\\ 1 & 1 & 1 & 0\\ 1 & 1 & 0 & 1 \end{bmatrix} \end{equation*}

is shown to have the three eigenvalues \(\lambda=3,\,1,\,-1\text{.}\) Suppose we wanted a \(4\times 4\) matrix that has the three eigenvalues \(\lambda=4,\,0,\,-2\text{.}\) We can employ Theorem EPM by finding a polynomial that converts \(3\) to \(4\text{,}\) \(1\) to \(0\text{,}\) and \(-1\) to \(-2\text{.}\) Such a polynomial is called an interpolating polynomial, and in this example we can use

\begin{equation*} r(x)=\frac{1}{4}\left(x^2+4x-5\right)=\frac{1}{4}x^2+x-\frac{5}{4}\text{.} \end{equation*}

(a)

Verify that the polynomial \(r(x)\) converts the eigenvalues as advertised.

Solution.

We will not discuss how to concoct the interpolating polynomial, \(r(x)\text{,}\) but a text on numerical analysis should provide the details. For now, it should be routine to verify that \(r(3)=4\text{,}\) \(r(1)=0\) and \(r(-1)=-2\text{.}\)

(b)

In the style of Example PM, compute the matrix \(r(C)\text{.}\)

Solution.

Now compute

\begin{align*} r(C)&=\frac{1}{4}C^2+C-\frac{5}{4}I_4\\ &= \frac{1}{4} \begin{bmatrix} 3 & 2 & 2 & 2\\ 2 & 3 & 2 & 2\\ 2 & 2 & 3 & 2\\ 2 & 2 & 2 & 3 \end{bmatrix} + \begin{bmatrix} 1 & 0 & 1 & 1\\ 0 & 1 & 1 & 1\\ 1 & 1 & 1 & 0\\ 1 & 1 & 0 & 1 \end{bmatrix} -\frac{5}{4} \begin{bmatrix} 1 & 0 & 0 & 0\\ 0 & 1 & 0 & 0\\ 0 & 0 & 1 & 0\\ 0 & 0 & 0 & 1 \end{bmatrix} = \frac{1}{2} \begin{bmatrix} 1 & 1 & 3 & 3\\ 1 & 1 & 3 & 3\\ 3 & 3 & 1 & 1\\ 3 & 3 & 1 & 1 \end{bmatrix}\text{.} \end{align*}

(c)

Compute the eigenvalues of \(r(C)\) directly and verify that they are as expected.

Solution.

notice that the multiplicities are the same, and the eigenspaces of \(C\) and \(r(C)\) are identical.

T10.

Suppose that \(A\) is a square matrix. Prove that the constant term of the characteristic polynomial of \(A\) is equal to the determinant of \(A\text{.}\)

Solution.

Suppose that the characteristic polynomial of \(A\) is

\begin{equation*} \charpoly{A}{x}=a_0+a_1x+a_2x^2+\dots+a_nx^n\text{.} \end{equation*}

Then

\begin{align*} a_0 &=a_0+a_1(0)+a_2(0)^2+\dots+a_n(0)^n\\ &=\charpoly{A}{0}\\ &=\detname{A-0I_n}&& \knowl{./knowl/xref/definition-CP.html}{\text{Definition CP}}\\ &=\detname{A}\text{.} \end{align*}

T20.

Suppose that \(A\) is a square matrix. Prove that a single vector may not be an eigenvector of \(A\) for two different eigenvalues.

Solution.

Suppose that the vector \(\vect{x}\neq\zerovector\) is an eigenvector of \(A\) for the two eigenvalues \(\lambda\) and \(\rho\text{,}\) where \(\lambda\neq\rho\text{.}\) Then \(\lambda-\rho\neq 0\text{,}\) and we also have

\begin{align*} \zerovector &=A\vect{x}-A\vect{x}&& \knowl{./knowl/xref/property-AIC.html}{\text{Property AIC}}\\ &=\lambda\vect{x}-\rho\vect{x}&& \knowl{./knowl/xref/definition-EEM.html}{\text{Definition EEM}}\\ &=(\lambda-\rho)\vect{x}&& \knowl{./knowl/xref/property-DSAC.html}{\text{Property DSAC}}\text{.} \end{align*}

By Theorem SMEZV, either \(\lambda-\rho=0\) or \(\vect{x}=\zerovector\text{,}\) which are both contradictions.

T21.

Suppose that \(A\) is a square matrix of size \(n\) with an eigenvalue \(\lambda\text{.}\) For \(\alpha\in\complexes\text{,}\) prove that \(\lambda+\alpha\) is an eigenvalue of the matrix \(A+\alpha I_n\text{.}\)

Construct a proof that begins by employing a vector \(\vect{x}\) that is an eigenvector of \(A\) for \(\lambda\text{.}\)
Construct a proof that establishes that \(\lambda+\alpha\) is a root of the characteristic polynomial of \(A+\alpha I_n\text{.}\)

T22.

Suppose that \(U\) is a unitary matrix with eigenvalue \(\lambda\text{.}\) Prove that \(\lambda\) has modulus 1, i.e. \(\modulus{\lambda}=1\text{.}\) This says that all of the eigenvalues of a unitary matrix lie on the unit circle of the complex plane.

T30.

Theorem DCP tells us that the characteristic polynomial of a square matrix of size \(n\) has degree \(n\text{.}\) By suitably augmenting the proof of Theorem DCP prove that the coefficient of \(x^n\) in the characteristic polynomial is \((-1)^n\text{.}\)

T50.

Theorem EIM says that if \(\lambda\) is an eigenvalue of the nonsingular matrix \(A\text{,}\) then \(\frac{1}{\lambda}\) is an eigenvalue of \(\inverse{A}\text{.}\) Write an alternate proof of this theorem using the characteristic polynomial and without making reference to an eigenvector of \(A\) for \(\lambda\text{.}\)

Solution.

Since \(\lambda\) is an eigenvalue of a nonsingular matrix, \(\lambda\neq 0\) (Theorem SMZE). \(A\) is invertible (Theorem NI), and so \(-\lambda A\) is invertible (Theorem MISM). Thus \(-\lambda A\) is nonsingular (Theorem NI) and \(\detname{-\lambda A}\neq 0\) (Theorem SMZD). We have

\begin{align*} \charpoly{\inverse{A}}{\frac{1}{\lambda}} &=\detname{\inverse{A}-\frac{1}{\lambda}I_n}&& \knowl{./knowl/xref/definition-CP.html}{\text{Definition CP}}\\ &=1\detname{\inverse{A}-\frac{1}{\lambda}I_n}&& \knowl{./knowl/xref/property-OCN.html}{\text{Property OCN}}\\ &=\frac{1}{\detname{-\lambda A}} \detname{-\lambda A}\detname{\inverse{A}-\frac{1}{\lambda}I_n}&& \knowl{./knowl/xref/property-MICN.html}{\text{Property MICN}}\\ &=\frac{1}{\detname{-\lambda A}} \detname{\left(-\lambda A\right)\left(\inverse{A}-\frac{1}{\lambda}I_n\right)}&& \knowl{./knowl/xref/theorem-DRMM.html}{\text{Theorem DRMM}}\\ &=\frac{1}{\detname{-\lambda A}} \detname{-\lambda A\inverse{A}-\left(-\lambda A\right)\frac{1}{\lambda}I_n}&& \knowl{./knowl/xref/theorem-MMDAA.html}{\text{Theorem MMDAA}}\\ &=\frac{1}{\detname{-\lambda A}} \detname{-\lambda I_n-\left(-\lambda A\right)\frac{1}{\lambda}I_n}&& \knowl{./knowl/xref/definition-MI.html}{\text{Definition MI}}\\ &=\frac{1}{\detname{-\lambda A}} \detname{-\lambda I_n+\lambda\frac{1}{\lambda}AI_n}&& \knowl{./knowl/xref/theorem-MMSMM.html}{\text{Theorem MMSMM}}\\ &=\frac{1}{\detname{-\lambda A}} \detname{-\lambda I_n+1AI_n}&& \knowl{./knowl/xref/property-MICN.html}{\text{Property MICN}}\\ &=\frac{1}{\detname{-\lambda A}} \detname{-\lambda I_n+AI_n}&& \knowl{./knowl/xref/property-OCN.html}{\text{Property OCN}}\\ &=\frac{1}{\detname{-\lambda A}} \detname{-\lambda I_n+A}&& \knowl{./knowl/xref/theorem-MMIM.html}{\text{Theorem MMIM}}\\ &=\frac{1}{\detname{-\lambda A}} \detname{A-\lambda I_n}&& \knowl{./knowl/xref/property-ACM.html}{\text{Property ACM}}\\ &=\frac{1}{\detname{-\lambda A}} \charpoly{A}{\lambda}&& \knowl{./knowl/xref/definition-CP.html}{\text{Definition CP}}\\ &=\frac{1}{\detname{-\lambda A}}\,0&& \knowl{./knowl/xref/theorem-EMRCP.html}{\text{Theorem EMRCP}}\\ &=0&& \knowl{./knowl/xref/property-ZCN.html}{\text{Property ZCN}}\text{.} \end{align*}

So \(\frac{1}{\lambda}\) is a root of the characteristic polynomial of \(\inverse{A}\) and so is an eigenvalue of \(\inverse{A}\text{.}\) This proof is due to Sara Bucht.

You have attempted of activities on this page.

Prev Top Next