3 double eraPas(
double al,
double ap,
double bl,
double bp)
36 y = sin(dl) * cos(bp);
37 x = sin(bp) * cos(ap) - cos(bp) * sin(ap) * cos(dl);
38 pa = ((x != 0.0) || (y != 0.0)) ? atan2(y, x) : 0.0;
double eraPas(double al, double ap, double bl, double bp)