ephemeris_correction
return t + ephemeris_correction(t);
return t - ephemeris_correction(t);
double ephemeris_correction(double t);