Algorithms for geodesics
  Algorithms for the computation of geodesics on an ellipsoid of revolution are given. These provide accurate, robust, and fast solutions to the direct and inverse geodesic problems and they allow differential and integral properties of geodesics to be computed. 
Charles F. F. Karney (SRI International)
physics.geo-ph physics.comp-ph