full Taylor method

Update for positions:

x_i^m+1 = x_i^m + h v_i^m + h²/2 a_i^m

Update for velocities:

v_i^m+1 = v_i^m + h a_i^m

This is the lowest order Nystrom method for an ODE where the accleration only depends on the primary variable (and not its derivative).

not in the Burden and Faires "Numerical Analysis" textbook
in Hairer, Norsett, and Wanner's "Solving Ordinary Differential Equations I: Nonstiff Problems": Section II.14
not in the Silly-pedia (and beware: this article is *not* the same thing though it is related)