Sun Declination angle

get_sigma(J)