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.
Merge pull request #1 from dnic/master

Changed the Euler method to a semi-implicit Euler method;