longTermEclipticPoleAngles function

(double, double) longTermEclipticPoleAngles(
  1. double t
)

P_A i Q_A (radiany), równ. (8).

Implementation

(double, double) longTermEclipticPoleAngles(double t) => (
  _series(t, [5851.607687, -0.1189000, -0.00028913, 101e-9], _ecliptic) *
      arcsecond,
  _series(
        t,
        [-1600.886300, 1.1689818, -0.00000020, -437e-9],
        _ecliptic,
        second: true,
      ) *
      arcsecond,
);