JDE of apogee nearest the given decimal year.
year
double apogee(double year) { final l = _La(year, 0.5); return _mean(l.t) + l._ac(); }