lunar_longitude
double lambda = lunar_longitude(t);
double phi = mod_f(lunar_longitude(t) - solar_longitude(t), 360);
double lon = lunar_longitude(t_u);
double lunar_longitude(double t);