Mean inclination for planet p at jde, in radians.
p
jde
double inc(int p, double jde) => toRad(horner(j2000Century(jde), _cMean[p].i));