With Equations 16 and 17, calculation of the velocities is immediate and the same scheme of forward differencing the system of ODEs, Eq. 15, applies. Results are illustrated below.
