next up previous contents index
Next: Example 3 (from Program Up: Examples Previous: Example 1   Contents   Index

Example 2 (from Program LA_SBEV_EXAMPLE)

The results below are computed with $\epsilon = 1.19209 \times 10^{-7}$.

\begin{displaymath}
A = \left( \begin{array}{rrrrr}
-6 & -2 & -5 & 0 & 0 \\
...
...& 10 & 3 & 6 \\
0 & 0 & -1 & 6 & -5 \\
\end{array} \right),
\end{displaymath}

Array ${\bf AB}$ on entry (${\bf UPLO} =$ 'U' by default):

\begin{displaymath}
\begin{array}{cc} {\bf AB} \\
\begin{array}{\vert rrrrr\v...
... 6 \\
-6 & 2 & -1 & 3 & -5 \\
\hline \end{array} \end{array}\end{displaymath}

The call:
CALL LA_SBEV( AB, W )

Array ${\bf W}$ on exit and the eigenvalues of $A$:

\begin{displaymath}
\begin{array}{cc} {\bf W} \\
\begin{array}{\vert r\vert} ...
...25 \\ -3.71265 \\ 1.20058 \\ 16.1449 \\
\end{array} \right).
\end{displaymath}



Susan Blackford 2001-08-19