phaseAngle3 function

double phaseAngle3(
  1. double jde
)

Quick phase angle from JDE (less accurate).

Implementation

double phaseAngle3(double jde) {
  final t = j2000Century(jde);
  final d = mod2pi(toRad(horner(t, [297.8501921, 445267.1114034,
      -0.0018819, 1 / 545868.0, -1 / 113065000.0])));
  final m = mod2pi(toRad(horner(t, [357.5291092, 35999.0502909,
      -0.0001535, 1 / 24490000.0])));
  final mp = mod2pi(toRad(horner(t, [134.9633964, 477198.8675055,
      0.0087414, 1 / 69699.0, -1 / 14712000.0])));
  return math.pi - d + toRad(
      -6.289 * math.sin(mp) +
      2.1 * math.sin(m) +
      -1.274 * math.sin(2 * d - mp) +
      -0.658 * math.sin(2 * d) +
      -0.214 * math.sin(2 * mp) +
      -0.11 * math.sin(d));
}