longTermEquatorPoleAngles function

(double, double) longTermEquatorPoleAngles(
  1. double t
)

X_A i Y_A (radiany), równ. (9).

Implementation

(double, double) longTermEquatorPoleAngles(double t) => (
  _series(t, [5453.282155, 0.4252841, -0.00037173, -152e-9], _equator) *
      arcsecond,
  _series(
        t,
        [-73750.930350, -0.7675452, -0.00018725, 231e-9],
        _equator,
        second: true,
      ) *
      arcsecond,
);