#include "AnT.hpp"

#define a     parameters[0]
#define b     parameters[1]
#define X     currentState[0]
#define Xt    delayState[0]

bool mg (const Array<real_t>& currentState,
	 const Array<real_t>& delayState,
	 const Array<real_t>& parameters,
	 Array<real_t>& rhs)
{
  rhs[0] = a * Xt/(1+pow(Xt,10.0)) - b * X;
  return true;
}

#undef a
#undef b
#undef X
#undef Xt

extern "C" 
{ 
void connectSystem () 
    { 
	DDE_Proxy::systemFunction = mg; 
    } 
}