radius function

double radius(
  1. double t
)

Sun-Earth distance in AU.

Implementation

double radius(double t) {
  final sun = trueSun(t);
  final e = eccentricity(t);
  return 1.000001018 * (1 - e * e) / (1 + e * math.cos(sun.anomaly));
}