double calcSunDeclination(double t) { double e = calcObliquityCorrection(t); double lambda = calcSunApparentLong(t); double sint = sin(degToRad(e)) * sin(degToRad(lambda)); double theta = radToDeg(asin(sint)); return theta; // in degrees }
/* Purpose: calculate the declination of the sun (degrees) */ double calcSunDeclination(double t) { return rad2deg(asin(sin(deg2rad(calcObliquityCorrection(t))) * sin(deg2rad(calcSunApparentLong(t))))); }