Astronomical algorithm