Example 2 (from Program LA_SBEV_EXAMPLE)

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

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

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

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

The call:
CALL LA_SBEV( AB, W )

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

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