Polynomial approximation of ΔT for calendar years before 948.
double polyBefore948(double year) { return horner(_c2000(year), [2177, 497, 44.1]); }