are reserved for each transfer function, hence, the real
          parts of the zeros for the (i,j) transfer function
          are stored starting from the location ((j-1)*P+i-1)*NPZ+1.
          Pairs of complex conjugate zeros are stored in consecutive
          memory locations. Note that only the first NZ(i,j) entries
          are initialized for the (i,j) transfer function.

  ZEROSI  (output) DOUBLE PRECISION array, dimension (P*M*NPZ)
          This array contains the imaginary parts of the zeros o