Verlet integration (French pronunciation: [vɛʁˈlɛ]) is a numerical method used to integrate Newton’s equations of motion[1] . It is frequently used to calculate trajectories of particles in molecular dynamics simulations and computer graphics. The Verlet integrator offers greater stability, as well as other properties that are important in physical systems such as time-reversibility and preservation of the symplectic form on phase space, at no significant additional cost over the simple Euler method. Verlet integration was used by Carl Størmer to compute the trajectories of particles moving in a magnetic field (hence it is also called Störmer’s method)[2] and was popularized in molecular dynamics by French physicist Loup Verlet in 1967.