equation_of_time
return t + equation_of_time(ut);
return t - equation_of_time(ut);
double equation_of_time(double t);