This is a port of the algorithm ported from the 'Astronomical Computing' column of Sky & Telescope, by Bradley E. Schaefer.