NSMP-by-L part of this array must contain the
          set of output samples,
          Y = ( Y(1,1),...,Y(1,L); ...; Y(NSMP,1),...,Y(NSMP,L) ).

  LDY     INTEGER
          The leading dimension of array Y.  LDY &gt;= MAX(1,NSMP).

  X       (input/output) DOUBLE PRECISION array dimension (LX)
          On entry, if INIT = 'L', the leading (NN*(L+2) + 1)*L part
          of this array must contain the initial parameters for
          the nonlinear part of the system.
          On entry, if INIT = 'S', the