mod_f
return mod_f(poly(c, coef, nitems(coef)), 360);
return mod_f(poly(c, coef, nitems(coef)), 360);
return mod_f(poly(c, coef, nitems(coef)), 360);
return mod_f(poly(c, coef, nitems(coef)), 360);
return mod_f(poly(c, coef, nitems(coef)), 360);
return mod_f(poly(c, coef, nitems(coef)), 360);
return mod_f((L_prime + correction + venus + jupiter +
double H = mod_f(theta + longitude - alpha, 360);
double phi = mod_f(lunar_longitude(t) - solar_longitude(t), 360);
double phi2 = mod_f((t - tn) / mean_synodic_month, 1) * 360;
double tau = t + rate * mod_f(phi - phase, 360);
return mod_f(lambda + ab + nu, 360);
double tau = t + rate * mod_f(lambda - lon, 360);
double tau = t - rate * mod_f(lon - lambda, 360);
double H = mod_f(theta + longitude - alpha, 360);
if (mod_f(f(x) - y, 360) < 180.0)
return a + mod_f(x - a, b - a);
return mod_f(v * 180.0 / M_PI, 360);