#include "AnT.hpp"

#define eps     parameters[0]
#define gamma   parameters[1]
#define kappa   parameters[2]

#define u     currentState[0]
#define s     currentState[1]

bool HakenZwanzig (const Array<real_t>& currentState,
		   const Array<real_t>& parameters,
		   Array<real_t>& rhs)
{
  rhs[0] = eps*u + gamma*u*s;
  rhs[1] = -s + kappa*u*u;

  return true;
}

#undef eps
#undef gamma
#undef kappa

#undef u
#undef s

extern "C"
{
  void connectSystem ()
  {
    ODE_Proxy::systemFunction = HakenZwanzig;
  }
}