"1/4 x^4 - 1/2 x^2 + 1/2 y^2", // function expression
"x^3-x", "y"); // derivatives
// integrator
IDoubleArray masses = doublesNew.arrayFrom(1.0, 1.0);
double dt = 0.1, gamma = 1, kT = 0.2;
IIntegratorThermostatted langevin = dynNew.langevinLeapFrog(pot, masses, 0.1, gamma, kT);
// run
IDoubleArray x0 = doublesNew.arrayFrom(0,0);
int nsteps = 100000, nsave = 10;
IDataSequence seq = dyn.run(x0, langevin, nsteps, nsave);