#include "AnT.hpp"

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

bool power_law_map (const Array<real_t>& currentState,
		    const Array<real_t>& parameters,
		    Array<real_t>& RHS)
{
    real_t h = pow (0.5, b);
    RHS[0] = a / h * (h - pow(fabs (X-0.5), b) );
    return true;
}

#undef a
#undef b 
#undef X    

extern "C" 
{

void connectSystem ()
{
  MapProxy::systemFunction = power_law_map;
}

}