rray S.  LDS &gt;= MAX(1,N).

  R       (output) DOUBLE PRECISION array, dimension (LDR,N)
          The leading N-by-N upper triangular part of this array
          contains the Cholesky factor R of the frequency-weighted
          observability Grammian Q = R'*R. See METHOD.

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

</PRE>
<B>Workspace</B>
<PRE>
  DWORK   DOUBLE PRECISION array, dimension (LDWORK)
          On exit, if INFO = 0, DWORK(1) returns the optimal value