private void CalculateEquatorialPosition()
{
double epsilon0 = GetDeltaEpsilon();
double theta = GetGeometricMeanLongitude();
double ar = Math.Atan2(Math.Cos(epsilon0 * M.DEG2RAD) * Math.Sin(theta * M.DEG2RAD), Math.Cos(theta * M.DEG2RAD)) * M.RAD2DEG;
double dec = Math.Asin(Math.Sin(epsilon0 * M.DEG2RAD) * Math.Sin(theta * M.DEG2RAD)) * M.RAD2DEG;
equatorialCoords.RA = HourAngle.FromDecimalDegrees(ar);
equatorialCoords.Declination = new DegreesAngle(dec);
}
