solar.c
static inline double
sun_equation_of_centre(double tm)
{
	double a = sun_geometric_mean_anomaly(tm), rc;
	rc  = sin(1.0 * a) * (pow(-0.000014 * tm, 2.0) - 0.004817 * tm + 1.914602);
	rc += sin(2.0 * a) * (-0.000101 * tm + 0.019993);
	rc += sin(3.0 * a) * 0.000289;
	return radians(rc);
}
/**
 * Calculates the Sun's equation of the centre, the difference
 * between the true anomaly and the mean anomaly
 * 
 * @param   tm  The time in Julian Centuries
 * @return      The Sun's equation of the centre, in radians
 */
Calculates the Sun's equation of the centre, the difference
between the true anomaly and the mean anomaly

@param   tm  The time in Julian Centuries
@return      The Sun's equation of the centre, in radians
Calculates the Sun ' s equation of the centre the difference between the true anomaly and the mean anomaly @param tm The time in Julian Centuries @return The Sun ' s equation of the centre in radians
static inline double
sun_equation_of_centre(double tm)
{
	double a = sun_geometric_mean_anomaly(tm), rc;
	rc  = sin(1.0 * a) * (pow(-0.000014 * tm, 2.0) - 0.004817 * tm + 1.914602);
	rc += sin(2.0 * a) * (-0.000101 * tm + 0.019993);
	rc += sin(3.0 * a) * 0.000289;
	return radians(rc);
}
static inline double sun_equation_of_centre ( double tm ) { double a = sun_geometric_mean_anomaly ( tm ) , rc ; rc = sin ( 1.0 * a ) * ( pow ( -0.000014 * tm , 2.0 ) - 0.004817 * tm + 1.914602 ) ; rc += sin ( 2.0 * a ) * ( -0.000101 * tm + 0.019993 ) ; rc += sin ( 3.0 * a ) * 0.000289 ; return radians ( rc ) ; }
Calculates the Sun's equation of the centre, the difference
 between the true anomaly and the mean anomaly
Calculates the Sun ' s equation of the centre the difference between the true anomaly and the mean anomaly