       The array U may be identified with B in the calling
          statement, if B is properly dimensioned, and the
          intermediate results returned in B are not needed.

  LDU     INTEGER
          The leading dimension of array U.  LDU &gt;= MAX(1,N).

  SCALE   (output) DOUBLE PRECISION
          The scale factor, scale, set less than or equal to 1 to
          prevent the solution overflowing.

</PRE>
<B>Workspace</B>
<PRE>
  DWORK   DOUBLE PRECISION array, dimension (LDWORK)
          On exit,