A common task in Lagrangian oceanography is to calculate a large number of drifter trajectories from a velocity field pre-calculated with an ocean model. Mathematically, this is simply numerical integration of an Ordinary ...