dv/dt = g - k v^2