-by-N part of
          this array contains the matrix of accumulated orthogonal
          column transformations performed.
          If JOBZ = 'N', the array Z is not referenced and can be
          supplied as a dummy array (i.e. set parameter LDZ = 1 and
          declare this array to be Z(1,1) in the calling program).

  LDZ     INTEGER
          The leading dimension of array Z. If JOBZ = 'U' or
          JOBZ = 'I', LDZ &gt;= MAX(1,N); if JOBZ = 'N', LDZ &gt;= 1.

  ISTAIR  (input/output) INTEGER ar