solar_anomaly
double solar_anomaly = poly(c, coef_sa, nitems(coef_sa));
arg1->x * solar_anomaly +
double M = solar_anomaly(c);
double M = solar_anomaly(c);
double M = solar_anomaly(c);