Example 1 (from Program LA_POSV_EXAMPLE)

$$A = \left( \begin{array}{rrrrr} 38 & 3 & 4 & 6 & 5 \\ 3 ... ... \\ 80 & 160 & 240 \\ 102 & 204 & 306 \end{array} \right).$$

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

$$\hspace{-1.00 cm} \begin{array}{c} {\bf A} \\ \begin{arr... ...60 & 240 \\ 102 & 204 & 306 \\ \hline \end{array} \end{array}$$

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

Arrays A and B on exit:

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

$$\begin{array}{c} {\bf B} \\ \begin{array}{\vert rrr\vert}... ...1.00000 & 2.00000 & 3.00000 \\ \hline \end{array} \end{array}$$

Matrices $U$ and $X$, where $A = U^TU$ and $X$ is the solution of the system $A\,X = B$:

$$U = \left(\begin{array}{lllll} 6.16441 & 4.86664 \times 10... ...\; 7.35919 & 1.13404 \\ & & & & 8.40898 \end{array} \right)$$

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