Longitude of perihelion of Earth's orbit.
double perihelionEarth(double t) { return toRad(horner(t, [102.93735, 1.71946, 0.00046])); }