%N eqeq %F $MIR/src/subs/ephem.for
%D Equation of the equinox
%P rjs
%: utilities
%B
        double precision function eqeq(jday)

        implicit none
        double precision jday

  Return the equation of the equinox.
  To convert a mean sidereal time to an apparent sidereal time,
  add the equation of the equinox.
  e.g.
    LAST = LMST + eqeq(jday)

  Input:
    jday       Julian date.
  Output:
    eqeq       The equation of the equinox, in radians.
