Important: The linalg package has been deprecated. Use the superseding command LinearAlgebra[Eigenvectors], instead.

- For information on migrating linalg code to the new packages, see examples/LinearAlgebraMigration.

The procedure eigenvectors computes the eigenvalues and eigenvectors of A.  That is, for each eigenvalue lambda of A it solves the linear system $\left(I\mathrm{lambda}-A\right)X=0$ for X.

The result returned is a sequence of lists of the form $\&lsqb;\mathrm{ei},\mathrm{mi},\{v\&lsqb;1,i\&rsqb;,\mathrm{tripledotplaceholderv}\&lsqb;\mathrm{ni},i\&rsqb;\}\&rsqb;$, where the $\mathrm{ei}$ are the eigenvalues, mi their algebraic multiplicities, $\{v\&lsqb;1,i\&rsqb;,...,v\&lsqb;\mathrm{ni},i\&rsqb;\}$ is a set of basis vectors for the eigenspace corresponding to $\mathrm{ei}$, and $1<=\mathrm{ni}<=\mathrm{mi}$ is the dimension of the eigenspace. The first example below shows how to extract the eigenvalues and their eigenvectors from the result.

Numeric Case: If the matrix A contains any floating-point (decimal) numbers, the eigenvectors are computed numerically.  A standard numerical algorithm is used.  All floating point arithmetic is done at Digits digits of precision.  Note that the matrix entries on input must all be all of type numeric or complex(numeric).

Symbolic Case: Otherwise the eigenvalues and eigenvectors are computed symbolically (exactly). First the characteristic polynomial is computed and solved for its roots (the eigenvalues) lambda[i] symbolically.  Then for each eigenvalue, a basis for its eigenspace is computed by computing the null space of the characteristic matrix $A-{\mathrm{lambda}}_{i}I$.

Although eigenvalues and eigenvectors can in principle be computed for a matrix over any field F, one must be able to compute effectively with the roots of a univariate polynomial over F. For this reason, this routine is limited at present to the field of algebraic numbers or algebraic functions.  Therefore, the input A must be a matrix of rationals, or polynomials, or a matrix of algebraic numbers or algebraic functions represented using Maple's RootOf or radical notation.

If the optional second argument is 'radical' is specified, Maple will try to express the eigenvalues and eigenvectors in terms of exact radicals. Note this may fail if any factor of the characteristic polynomial has degree 5 or higher. If the optional second argument is 'implicit' then Maple will express the eigenvalues in terms of an algebraic extension of any nonlinear factors of the characteristic polynomial using Maple's RootOf notation. This is always possible. If no optional second argument is specified, then Maple will use radicals as the default.  The examples below will make these differences clear.

The command with(linalg,eigenvectors) allows the use of the abbreviated form of this command.

$\mathrm{with}\left(\mathrm{linalg}\right)\&colon;$

$A≔\mathrm{matrix}\left(3\&comma;3\&comma;\left[1\&comma;-3\&comma;3\&comma;3\&comma;-5\&comma;3\&comma;6\&comma;-6\&comma;4\right]\right)$

$A≔\left[\begin{array}{ccc}1& \mathrm{-3}& 3\\ 3& \mathrm{-5}& 3\\ 6& \mathrm{-6}& 4\end{array}\right]$

$e≔\mathrm{eigenvalues}\left(A\right)$

$e≔4,\mathrm{-2},\mathrm{-2}$

$v≔\left[\mathrm{eigenvectors}\left(A\right)\right]$

$v≔\left[\left[4\&comma;1\&comma;\left\{\left[\begin{array}{ccc}1& 1& 2\end{array}\right]\right\}\right]\&comma;\left[\mathrm{-2}\&comma;2\&comma;\left\{\left[\begin{array}{ccc}1& 1& 0\end{array}\right]\&comma;\left[\begin{array}{ccc}\mathrm{-1}& 0& 1\end{array}\right]\right\}\right]\right]$

$v\left[1\right]\left[1\right]$

$4$

$v\left[1\right]\left[2\right]$

$1$

$v\left[1\right]\left[3\right]$

$\left\{\left[\begin{array}{ccc}1& 1& 2\end{array}\right]\right\}$

$v\left[2\right]\left[1\right]$

$\mathrm{-2}$

$v\left[2\right]\left[2\right]$

$2$

$A≔\mathrm{toeplitz}\left(\left[1.0\&comma;2.0\&comma;3.0\right]\right)$

$A≔\left[\begin{array}{ccc}1.0& 2.0& 3.0\\ 2.0& 1.0& 2.0\\ 3.0& 2.0& 1.0\end{array}\right]$

$\mathrm{eigenvalues}\left(A\right)$

$\mathrm{-2.00000000000000},\mathrm{-0.701562118716425},5.70156211871643$

$\mathrm{eigenvectors}\left(A\right)$

$\left[\mathrm{-0.7015621183}\&comma;1\&comma;\left\{\left[\begin{array}{ccc}\mathrm{-0.3645129330}& 0.8568900996& \mathrm{-0.3645129334}\end{array}\right]\right\}\right],\left[\mathrm{-2.000000002}\&comma;1\&comma;\left\{\left[\begin{array}{ccc}0.7071067811& 1.47\times {10}^{\mathrm{-10}}& \mathrm{-0.7071067809}\end{array}\right]\right\}\right],\left[5.701562115\&comma;1\&comma;\left\{\left[\begin{array}{ccc}0.6059128004& 0.5154991348& 0.6059128001\end{array}\right]\right\}\right]$

$A≔\mathrm{toeplitz}\left(\left[1\&comma;2\&comma;3\right]\right)$

$A≔\left[\begin{array}{ccc}1& 2& 3\\ 2& 1& 2\\ 3& 2& 1\end{array}\right]$

$\mathrm{eigenvalues}\left(A\right)$

$\mathrm{-2},\frac{5}{2}-\frac{\sqrt{41}}{2},\frac{5}{2}+\frac{\sqrt{41}}{2}$

$\mathrm{eigenvectors}\left(A\right)$

$\left[\frac{5}{2}+\frac{\sqrt{41}}{2}\&comma;1\&comma;\left\{\left[\begin{array}{ccc}1& -\frac{3}{4}+\frac{\sqrt{41}}{4}& 1\end{array}\right]\right\}\right],\left[\frac{5}{2}-\frac{\sqrt{41}}{2}\&comma;1\&comma;\left\{\left[\begin{array}{ccc}1& -\frac{3}{4}-\frac{\sqrt{41}}{4}& 1\end{array}\right]\right\}\right],\left[\mathrm{-2}\&comma;1\&comma;\left\{\left[\begin{array}{ccc}\mathrm{-1}& 0& 1\end{array}\right]\right\}\right]$

$\mathrm{eigenvectors}\left(A\&comma;'\mathrm{implicit}'\right)$

$\left[\mathrm{RootOf}\left({\mathrm{\_Z}}^2-5\mathrm{\_Z}-4\right)\&comma;1\&comma;\left\{\left[\begin{array}{ccc}1& \frac{\mathrm{RootOf}\left({\mathrm{\_Z}}^2-5\mathrm{\_Z}-4\right)}{2}-2& 1\end{array}\right]\right\}\right],\left[\mathrm{-2}\&comma;1\&comma;\left\{\left[\begin{array}{ccc}\mathrm{-1}& 0& 1\end{array}\right]\right\}\right]$

$T≔\mathrm{toeplitz}\left(\left[a\&comma;b\&comma;c\right]\right)$

$T≔\left[\begin{array}{ccc}a& b& c\\ b& a& b\\ c& b& a\end{array}\right]$

$\mathrm{eigenvalues}\left(T\right)$

$a-c,a+\frac{c}{2}+\frac{\sqrt{8b^2+c^2}}{2},a+\frac{c}{2}-\frac{\sqrt{8b^2+c^2}}{2}$

$\mathrm{eigenvectors}\left(T\&comma;'\mathrm{implicit}'\right)$

$\left[a-c\&comma;1\&comma;\left\{\left[\begin{array}{ccc}\mathrm{-1}& 0& 1\end{array}\right]\right\}\right],\left[\mathrm{RootOf}\left({\mathrm{\_Z}}^2+\left(-2a-c\right)\mathrm{\_Z}+a^2+ca-2b^2\right)\&comma;1\&comma;\left\{\left[\begin{array}{ccc}1& -\frac{-\mathrm{RootOf}\left({\mathrm{\_Z}}^2+\left(-2a-c\right)\mathrm{\_Z}+a^2+ca-2b^2\right)+a+c}{b}& 1\end{array}\right]\right\}\right]$

$T≔\mathrm{toeplitz}\left(\left[\mathrm{sqrt}\left(2\right)\&comma;a\&comma;3\right]\right)$

$T≔\left[\begin{array}{ccc}\sqrt{2}& a& 3\\ a& \sqrt{2}& a\\ 3& a& \sqrt{2}\end{array}\right]$

$\mathrm{eigenvectors}\left(T\right)$

$\left[\frac{3}{2}+\sqrt{2}+\frac{\sqrt{8a^2+9}}{2}\&comma;1\&comma;\left\{r\right\}\right],\left[\frac{3}{2}+\sqrt{2}-\frac{\sqrt{8a^2+9}}{2}\&comma;1\&comma;\left\{r\right\}\right],\left[\sqrt{2}-3\&comma;1\&comma;\left\{r\right\}\right]$

$B≔\mathrm{matrix}\left(4\&comma;4\&comma;\left[0\&comma;1\&comma;0\&comma;0\&comma;-u^2\&comma;2u\&comma;0\&comma;0\&comma;-us\&comma;s\&comma;u\&comma;0\&comma;-ut\&comma;t\&comma;0\&comma;u\right]\right)$

$B≔\left[\begin{array}{cccc}0& 1& 0& 0\\ -u^2& 2u& 0& 0\\ -us& s& u& 0\\ -ut& t& 0& u\end{array}\right]$

$\mathrm{eigenvectors}\left(B\right)$

$\left[u\&comma;4\&comma;\left\{\left[\begin{array}{cccc}\frac{1}{u}& 1& 0& 0\end{array}\right]\&comma;\left[\begin{array}{cccc}0& 0& 1& 0\end{array}\right]\&comma;\left[\begin{array}{cccc}0& 0& 0& 1\end{array}\right]\right\}\right]$