r function G(i,j).

  LDNP    INTEGER
          The leading dimension of array NP.  LDNP &gt;= max(1,P).

  ZEROSR  (output) DOUBLE PRECISION array, dimension (P*M*NPZ)
          This array contains the real parts of the zeros of the
          transfer function matrix G. The real parts of the zeros
          are stored in a column-wise order, i.e., for the transfer
          functions (1,1), (2,1), ..., (P,1), (1,2), (2,2), ...,
          (P,2), ..., (1,M), (2,M), ..., (P,M); NPZ memory locations
          