#include "AnT.hpp"

#define alpha parameters[0]
#define X     currentState[0]
#define Y     currentState[1]


bool harm_osc (const Array<real_t>& currentState,
	       const Array<real_t>& parameters,
	       Array<real_t>& rhs)
{
  rhs[0] = Y;
  rhs[1] = -X - alpha*Y;

  return true;
}

#undef alpha 
#undef X
#undef Y

extern "C" {

void connectSystem ()
{
  ODE_Proxy::systemFunction = harm_osc;
}

}