trueEquatorial function

({double dec, double ra}) trueEquatorial(
  1. double jde
)

True equatorial coordinates of the Sun.

Implementation

({double ra, double dec}) trueEquatorial(double jde) {
  final sun = trueSun(j2000Century(jde));
  final eps = nut.meanObliquity(jde);
  final sS = math.sin(sun.lon);
  final cS = math.cos(sun.lon);
  final sEps = math.sin(eps);
  final cEps = math.cos(eps);
  return (
    ra: mod2pi(math.atan2(cEps * sS, cS)),
    dec: math.asin(sEps * sS),
  );
}