be equal to
          (EVR(K) + i*EVI(K))/EVQ(K) if EVQ(K) .NE. ZERO. Note that
          the quotient may over- or underflow. If P has a degree d
          less than 3, then 3-d computed roots will be infinite.
          EVQ(K) &gt;= 0.

</PRE>
<B>Workspace</B>
<PRE>
  DWORK   DOUBLE PRECISION array, DIMENSION (LDWORK)
          On exit, if LDWORK = -1 on input, then DWORK(1) returns
          the optimal value of LDWORK.

  LDWORK  INTEGER
          The dimension of the array DWORK.  LDWORK &gt;= 42.

   