The GNU Octave mapping toolbox includes functions for computing angular distances and displacements between points on a sphere, as well as the great circle azimuth.