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