Example 2 (from Program LA_SPSV_EXAMPLE)

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

$$A = \left( \begin{array}{rrrrr} 0 & 2 & 3 & 5 & 4 \\ 2 &... ...egin{array}{r} 14 \\ 19 \\ 63 \\ 69 \\ 96 \end{array} \right)$$

Arrays AP and B on entry:

$$\begin{array}{c} {\bf AP} \\ \begin{array}{\vert rrrrrrrr... ...e 14 \\ 19 \\ 63 \\ 69 \\ 96 \\ \hline \end{array} \end{array}$$

The call:
CALL LA_SPSV( AP, B, 'L', IPIV )

AP, B and IPIV on exit:

$$\begin{array}{c} {\bf AP} \\ \begin{array}{l} \begin{arr... ...^{-1} & -4.01098 \\ \hline \end{array} \end{array} \end{array}$$

$$\begin{array}{c} {\bf B} \\ \begin{array}{\vert r\vert} \... ...line -4 \\ -4 \\ 3 \\ 4 \\ 5 \\ \hline \end{array} \end{array}$$

The solution of the system $A\,x = b$ is:

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