I've been tooling around with this, and found that the method known in some quarters as "Leapfrog" gives pretty good results even with modest computing power. The method is the same as described in Feynman Lectures volume 1 chapter 9. It seems to be the same as what some people call Verlet...