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

Example 2 (from Program LA_SYSV_EXAMPLE)

Matrix $A$ as in Example 1 and the first column of matrix $B$ in Example 1.
Arrays ${\bf A}$ and ${\bf B(:,1)}$ on entry:

\begin{displaymath}
\begin{array}{c} {\bf A} \\
\begin{array}{\vert rrrrr\ver...
... 14 \\ 19 \\ 21 \\ 23 \\ 32 \\
\hline \end{array} \end{array}\end{displaymath}

Elements marked $\ast$ are not used by the routine.

The call:
CALL LA_SYSV( A, B
(:,1), 'L', IPIV, INFO )

${\bf A}$, ${\bf B(:,1)}$, ${\bf IPIV}$ and ${\bf INFO}$ on exit:

\begin{displaymath}
\begin{array}{c} {\bf A} \\
\begin{array}{\vert lllll\ver...
...97 \times 10^{-1} &-4.01098 \\
\hline \end{array} \end{array}\end{displaymath}


\begin{displaymath}
\begin{array}{c} {\bf B(:,1)} \\
\begin{array}{\vert r\ve...
...} \hspace{0.50 cm}
\begin{array}{c} {\bf INFO} = 0 \end{array}\end{displaymath}

The solution of the system $A\,x = b$ is:

\begin{displaymath}x = \left( \begin{array}{r}
1.00000 \\ 1.00000 \\ 1.00000 \\ 1.00000 \\ 1.00000
\end{array} \right). \end{displaymath}



Susan Blackford 2001-08-19