Example 2 (from Program LA_PPSV_EXAMPLE)

The results below are computed with $\epsilon = 1.19209 \times 10^{-7}$.

$$A = \left( \begin{array}{rrrrr} 38 & 3 & 4 & 6 & 5 \\ 3 ... ...{array}{rrr} 56 \\ 71 \\ 69 \\ 80 \\ 102 \end{array} \right).$$

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

$$\begin{array}{c} {\bf AP} \\ \begin{array}{\vert rrrrrrrr... ... 56 \\ 71 \\ 69 \\ 80 \\ 102 \\ \hline \end{array} \end{array}$$

The call:

CALL LA_PPSV( AP, B, 'L' )
Arrays ${\bf AP}$ and ${\bf B}$ on exit:

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

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

$$\hspace{-0.80 cm} L = \left( \begin{array}{lllll} 6.16441... ...000 \\ 1.00000 \\ 1.00000 \\ 1.00000 \\ \end{array} \right).$$