lunar_phase
double phase = lunar_phase(t);
return invert_angular(lunar_phase, phi, a, b);
double phi = lunar_phase(t);
double phi = lunar_phase(t);
bool waning = lunar_phase(t) > 180.0;
bool waxing = lunar_phase(t) < 180.0;
double phi = lunar_phase(t_u);
double lunar_phase(double t);