Changed the Euler method to a semi-implicit Euler method;
The energy drift issue is mainly fixed with testing
done on the Earth-Sun system.