#include "AnT.hpp"


#define a parameters[0]
#define b parameters[1]
#define c parameters[2]
#define X currentState[0]
#define Y currentState[1]

bool circle_2d (const Array<real_t>&currentState,
		const Array<real_t>& parameters,
		Array<real_t>& rhs)
{
  
  rhs[0] = X + c + Y / Two_Pi;
  rhs[1] = b * Y - a * sin (Two_Pi * X);

  return true;
}

#undef a
#undef b
#undef c
#undef X
#undef Y

extern "C" 
{

void connectSystem ()
{
  MapProxy::systemFunction = circle_2d;

}

}