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

Example 1 (from Program LA_SYSV_EXAMPLE)


\begin{displaymath}
A = \left( \begin{array}{rrrrr}
0 & 2 & 3 & 5 & 4 \\
2 &...
...2 & 63 \\
23 & 46 & 69 \\
32 & 64 & 96 \end{array} \right)
\end{displaymath}

Arrays ${\bf A}$ and ${\bf B}$ on entry:

\begin{displaymath}
\begin{array}{c} {\bf A} \\
\begin{array}{\vert rrrrr\ver...
...3 & 46 & 69 \\
32 & 64 & 96 \\ \hline \end{array} \end{array}\end{displaymath}

Elements marked $\ast$ are not used by the routine. The call:
CALL LA_SYSV( A, B, IPIV
=IPIV )

${\bf A}$, ${\bf B}$ and ${\bf IPIV}$ on exit:

\begin{displaymath}
\begin{array}{c} {\bf A} \\
\begin{array}{\vert lllll\ver...
...\;\;\;\;\;\; \ast & 8.00000 \\
\hline \end{array} \end{array}\end{displaymath}


\begin{displaymath}
\begin{array}{cc} {\bf B} \\
\begin{array}{\vert rrr\vert...
...\hline 1 \\ 2 \\ 3 \\ 4 \\ 5 \\ \hline
\end{array} \end{array}\end{displaymath}

The solution of the system $ A\,X = B $ is:

\begin{displaymath}
X = \left( \begin{array}{rrr}
1.00000 & 2.00000 & 3.00000 ...
...3.00000 \\
1.00000 & 2.00000 & 3.00000
\end{array} \right).
\end{displaymath}



Susan Blackford 2001-08-19