q ( t + τ ) = 2 q ( t ) − q ( t − τ ) + τ 2 F ( q ( t ) ) / M
Originally discovered already by Newton: Essentially a symmetric (and symplectic!) combination of Euler forward and backward methods.
Verlet integration
sympletic, reversible time-discretization of Newton's equations of motion
Stoermer-Verlet integration