sion (2)
          Machine parameters:
          MACPAR(1)  (machine precision)*base, DLAMCH( 'P' );
          MACPAR(2)  safe minimum,             DLAMCH( 'S' ).
          This argument is not used for N = 2.

  Q       (output) DOUBLE PRECISION array, dimension (LDQ, N)
          The leading N-by-N part of this array contains the
          orthogonal transformation matrix Q.

  LDQ     INTEGER
          The leading dimension of the array Q.  LDQ &gt;= N.

  U       (output) DOUBLE PRECISION array, dimensi