COS
alpha = ATAN(TAN(lambda) * COS(epsilon)); /* 9 */
COS(latitude) * COS(*DEC) * COS(HA)); /* 13 */
(COS(HA) * SIN(latitude) - TAN(*DEC) * COS(latitude))); /* 14 */