// A simple LennardJones force field...

#if !defined(__LENNARDJONES_INCLUDED__)
#define __LENNARDJONES_INCLUDED__

#include "Vector.h"
#include "Particles.h"
#include "Integrator.h"

// ------------------------------------------------------------

#define LENNARDJONES_EPSILON 1.0
#define LENNARDJONES_SIGMA   0.1

// ------------------------------------------------------------

int    LennardJonesForces( int N, Particle P[] );
double LennardJonesPotential( int N, Particle P[] );

#endif