se</FONT></B>
<PRE>
  To compute a set of parameters for approximating a Wiener system
  in a least-squares sense, using a neural network approach and a
  Levenberg-Marquardt algorithm. Conjugate gradients (CG) or
  Cholesky algorithms are used to solve linear systems of equations.
  The Wiener system is represented as

     x(t+1) = A*x(t) + B*u(t)
     z(t)   = C*x(t) + D*u(t),

     y(t)   = f(z(t),wb(1:L)),

  where t = 1, 2, ..., NSMP, and f is a nonlinear function,
  evaluated by the SLICOT Library ro