invert_angular
return invert_angular(lunar_phase, phi, a, b);
return invert_angular(solar_longitude, lambda, a, b);
double invert_angular(double (*f)(double), double y, double a, double b);