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{array}{c} {\bf A} \\ \begin{array}{\vert rrrrr\ver... ... 56 \\ 71 \\ 69 \\ 80 \\ 102 \\ \hline \end{array} \end{array}$$

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:

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

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

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

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