velocity verlet as default propagation method
plot trajectory in time