public void MeanObliquityOfEclipticTest()
{
var date = new AASDate(1987, 4, 10, 0, 0, 0, true);
double obliquity = AASNutation.MeanObliquityOfEcliptic(date.Julian);
Assert.Equal(23.440946290957317, obliquity);
}
//obliquity of the ecliptic with correction
private double GetEpsilonCorrected()
{
double T = (jd - 2451545.0) / 36525;
double epsilon0 = AASNutation.MeanObliquityOfEcliptic(jd);
double omega = 125.04 - 1934.136 * T;
double deltaEpsilon = .00256 * Math.Cos(omega * M.DEG2RAD);
return(epsilon0 + deltaEpsilon);
}
//obliquity of the ecliptic
protected double GetDeltaEpsilon()
{
return(AASNutation.MeanObliquityOfEcliptic(jd));
}
//obliquity of the ecliptic
private double GetEpsilon()
{
return(AASNutation.MeanObliquityOfEcliptic(jd));
}