Example 2 (from Program LA_GESV_EXAMPLE)

${\bf A}$ on entry: As in Example 1.

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

The call:
CALL LA_GESV( A, B(:,1), IPIV, INFO )

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

$$\begin{array}{c} {\bf A} \\ \begin{array}{\vert lllll\ver... ...84 \times 10^{-1} & 1.62162 \\ \hline \end{array} \end{array}$$

$$\begin{array}{c} {\bf B(:,1)} \\ \begin{array}{\vert r\ve... ...} \hspace{1.00 cm} \begin{array}{c} {\bf INFO} = 0 \end{array}$$

Matrices $L$, $U$, $P$ and $x$, where $x$ is the solution of the system $A\,x = b$:

$$L = \left( \begin{array}{lllll} 1.00000 \\ 7.14286 \time... ...0^{-1} & 7.83784 \times 10^{-1} & 1.00000 \end{array} \right)$$

$$U = \left( \begin{array}{rrrrr} 7.00000 & 6.00000 & 8.0000... ... & 4.34078 & 4.29609 \\ & & & & 1.62162 \end{array} \right)$$

$$P = \left( \begin{array}{rrrrr} 0 & 0 & 1 & 0 & 0 \\ 0 &... ...000 \\ 1.00000 \\ 1.00000 \\ 1.00000 \\ \end{array} \right).$$