on (LDU, N)
          The leading N-by-N part of this array contains the
          orthogonal symplectic transformation matrix U.

  LDU     INTEGER
          The leading dimension of the array U.  LDU &gt;= N.

</PRE>
<B>Workspace</B>
<PRE>
  DWORK   DOUBLE PRECISION array, dimension (LDWORK)
          If N = 2 then DWORK is not referenced.

  LDWORK  INTEGER
          The length of the array DWORK.
          If N = 2 then LDWORK &gt;= 0; if N = 4 then LDWORK &gt;= 12.

</PRE>
<B>Error Indicator</B>
<PRE>
