tan_deg
double y = pow(tan_deg(epsilon/2), 2);
tan_deg(beta) * sin_deg(epsilon));
return (tan_deg(latitude) * tan_deg(delta) +