Example 1 (from Program LA_GELS_EXAMPLE)

$$A = \left( \begin{array}{rrrr} -3 & -1 & 6 & -5 \\ 5 & 4... ... 0 & 7 \\ 6 & -2 & -8 \\ -5 & 10 & 1 \end{array} \right)$$

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

$$\begin{array}{cc} {\bf A} \\ \begin{array}{\vert rrrr\ver... ... 6 & -2 & -8 \\ -5 & 10 & 1 \\ \hline \end{array} \end{array}$$

The call:
CALL LA_GELS( A, B )

B on exit:

$$\begin{array}{cc} {\bf B} \\ \begin{array}{\vert lll\vert... ....29547 & 12.6517 & -3.71209 \\ \hline \end{array} \end{array}$$

The solution matrix is:

$$X = \left( \begin{array}{lll} \;\;\; 1.48869 \times 10^{-1}... ...^{-1} & 4.17938 \times 10^{-1} & 1.32376 \end{array} \right).$$

The residual sums-of-squares are:

$$\left( \begin{array}{lll} 55.6579 & 165.295 & 23.1260 \end{array} \right).$$