solar.c
static double
mean_ecliptic_obliquity(double tm)
{
	double rc = pow(0.001813 * tm, 3.0) - pow(0.00059 * tm, 2.0) - 46.815 * tm + 21.448;
	return radians(23.0 + (26.0 + rc / 60.0) / 60.0);
}
/**
 * Calculates the mean ecliptic obliquity of the Sun's
 * apparent motion without variation correction
 * 
 * @param   tm  The time in Julian Centuries
 * @return      The uncorrected mean obliquity, in radians
 */
Calculates the mean ecliptic obliquity of the Sun's
apparent motion without variation correction

@param   tm  The time in Julian Centuries
@return      The uncorrected mean obliquity, in radians
Calculates the mean ecliptic obliquity of the Sun ' s apparent motion without variation correction @param tm The time in Julian Centuries @return The uncorrected mean obliquity in radians
static double
mean_ecliptic_obliquity(double tm)
{
	double rc = pow(0.001813 * tm, 3.0) - pow(0.00059 * tm, 2.0) - 46.815 * tm + 21.448;
	return radians(23.0 + (26.0 + rc / 60.0) / 60.0);
}
static double mean_ecliptic_obliquity ( double tm ) { double rc = pow ( 0.001813 * tm , 3.0 ) - pow ( 0.00059 * tm , 2.0 ) - 46.815 * tm + 21.448 ; return radians ( 23.0 + ( 26.0 + rc / 60.0 ) / 60.0 ) ; }
Calculates the mean ecliptic obliquity of the Sun's
 apparent motion without variation correction
Calculates the mean ecliptic obliquity of the Sun ' s apparent motion without variation correction