Added several integration methods

Also fixed a bug in kinetic energy calculation in energy conservation test.