next up previous contents index
Next: LA_GELSY Up: Examples Previous: Example 1 (from Program   Contents   Index

Example 2 (from Program LA_GELS_EXAMPLE)

${\bf A}$ on entry: As in Example 1.
${\bf B}(:,1)$ on entry:

\begin{array}{cc} {\bf B}(:,1) \\
\begin{array}{\vert r\v...
... \\ 8 \\ -2 \\ \circ \\ \circ \\ \hline \end{array} \end{array}\end{displaymath}

Elements marked $\circ$ need not be set on entry. This extra space is needed for the solution vector $x$.

The call:
CALL LA_GELS( A, B(:,1),
'T', INFO )

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

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

The minimum-norm solution of the system $A^T x = b$ is:

x = \left( \begin{array}{rr}
5.96220 \times 10^{-1} \\
...imes 10^{-1} \\
-0.26272 \times 10^{-1}
\end{array} \right).

Susan Blackford 2001-08-19