pub fn verlet_step<F, G>(
q: f64,
p: f64,
h: f64,
force: F,
kinetic_velocity: G,
) -> (f64, f64)Expand description
One Störmer–Verlet (velocity-Verlet / leapfrog) step — 2nd-order symplectic and
time-reversible. force(q) = -∂V/∂q, kinetic_velocity(p) = ∂T/∂p (= p/m).