   skew-Hamiltonian/Hamiltonian eigenproblems.
      ACM Trans. Math. Softw, 42, pp. 1-26, 2016.

</PRE>
<A name="Numerical Aspects"><B><FONT SIZE="+1">Numerical Aspects</FONT></B></A>
<PRE>
  If the algorithm does not converge in MAXIT = 30 iterations
  (INFO = 4), the tolerance must be increased, or the system is
  improper.

</PRE>
<A name="Comments"><B><FONT SIZE="+1">Further Comments</FONT></B></A>
<PRE>
  Setting POLES = 'P' usually saves some computational effort. The
  number of poles used is define