public void PositionAngleOfAxis()
{
double jd = 2448724.5;
CrdsEcliptical ecl = new CrdsEcliptical(new DMS("133* 10' 00''"), new DMS("-3* 13' 45''"));
double epsilon = 23.440636;
double deltaPsi = 0.004610;
double P = LunarEphem.PositionAngleOfAxis(jd, ecl, epsilon, deltaPsi);
Assert.AreEqual(15.08, P, 1e-2);
}
/// <summary>
/// Get position angle of axis for the Moon
/// </summary>
private double PAaxis(SkyContext c)
{
return(LunarEphem.PositionAngleOfAxis(c.JulianDay, c.Get(Ecliptical), c.Epsilon, c.NutationElements.deltaPsi));
}