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

Example 2 (from Program LA_POSV_EXAMPLE)

${\bf A}$ on entry: As in Example 1.
${\bf B}$ on entry: $B_{:,1}$, where $B$ is the input matrix in Example 1.

\begin{displaymath}
\begin{array}{c} {\bf A} \\
\begin{array}{\vert rrrrr\ver...
...
56 \\ 71 \\ 69 \\ 80 \\ 102 \\ \hline \end{array} \end{array}\end{displaymath}

Elements marked $\ast$ are not used by the routine. The call:
CALL LA_POSV( A, B
(:,1), = 'L'

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

\begin{displaymath}
\hspace{-0.80 cm}
\begin{array}{c} {\bf A} \\
\begin{arr...
...00000 \\ 1.00000 \\ 1.00000 \\
\hline \end{array} \end{array}\end{displaymath}

$L$ and $x$, where $A = L^TL$ and $x$ is the solution of the system $A\,x = b$:

\begin{displaymath}\hspace{-0.80 cm}
L = \left( \begin{array}{lllll}
6.16441 \...
...57024 \times 10^{-1} & 1.13404 & 8.40898
\end{array} \right), \end{displaymath}


\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