meanAnomaly function

double meanAnomaly(
  1. double t
)

Mean anomaly of the Sun.

t is centuries from J2000.

Implementation

double meanAnomaly(double t) {
  return toRad(horner(t, [357.52911, 35999.05029, -0.0001537]));
}