compute_velocity_ek function

Compute velocity via Engbert & Kliegl (2003) algorithm.