static double calcObliquityCorrection( double t) { double e0 = calcMeanObliquityOfEcliptic(t); double omega = 125.04 - 1934.136 * t; double e = e0 + 0.00256 * cos(degToRad(omega)); return e; // in degrees }
/* Purpose: calculate the corrected obliquity of the ecliptic (degrees) */ double calcObliquityCorrection(double t) { return calcMeanObliquityOfEcliptic(t) + 0.00256*cos(deg2rad(125.04-1934.136*t)); }