added output of total energy. Implicit Gauss methods display exceptionnal energy conservation, while explicit methods show a little drift. Energy very different when collisions occur (normal ... ?)