The Velocity Verlet integration method isn't much of a jump in complexity (http://en.wikipedia.org/wiki/Verlet_integration#Velocity_Ver...) and works _great_ for kinematic simulators.

As a cool feature, it maintains the energy of the system, so you can slow down, stop, and then rewind time and arrive at the same previous configuration, all by changing the delta T's inside the evaluation.