estimate_prior_solar_longitude
double approx = estimate_prior_solar_longitude(winter, t);
double estimate_prior_solar_longitude(double lambda, double t);