SORG2L generates an m by n real matrix Q with orthonormal columns,
 which is defined as the last n columns of a product of k elementary
 reflectors of order m

       Q  =  H(k) . . . H(2) H(1)

 as returned by SGEQLF.

          M is INTEGER
          The number of rows of the matrix Q. M >= 0.

          N is INTEGER
          The number of columns of the matrix Q. M >= N >= 0.

          K is INTEGER
          The number of elementary reflectors whose product defines the
          matrix Q. N >= K >= 0.

          A is REAL array, dimension (LDA,N)
          On entry, the (n-k+i)-th column must contain the vector which
          defines the elementary reflector H(i), for i = 1,2,...,k, as
          returned by SGEQLF in the last k columns of its array
          argument A.
          On exit, the m by n matrix Q.

          LDA is INTEGER
          The first dimension of the array A. LDA >= max(1,M).

          TAU is REAL array, dimension (K)
          TAU(i) must contain the scalar factor of the elementary
          reflector H(i), as returned by SGEQLF.

          WORK is REAL array, dimension (N)

          INFO is INTEGER
          = 0: successful exit
          < 0: if INFO = -i, the i-th argument has an illegal value