aphelion function

double aphelion(
  1. int p,
  2. double year
)

Approximate JDE of aphelion nearest decimal year for planet p.

Implementation

double aphelion(int p, double year) => _ap(p, year, true);