Totally revamped the code based on an email from Roland. Now we evolve the major axis a and the eccentricity e, not the energy and the angular momentum.