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.