sin_deg
double equation = (y * sin_deg(2*lambda) -
2 * eccentricity * sin_deg(anomaly) +
4 * eccentricity * y * sin_deg(anomaly) * cos_deg(2*lambda) -
1.25 * eccentricity*eccentricity * sin_deg(2*anomaly) -
0.5 * y*y * sin_deg(4*lambda)) / (2*M_PI);
return arcsin_deg(sin_deg(beta) * cos_deg(epsilon) +
cos_deg(beta) * sin_deg(epsilon) * sin_deg(lambda));
double y = (sin_deg(lambda) * cos_deg(epsilon) -
tan_deg(beta) * sin_deg(epsilon));
double extra = 0.000325 * sin_deg(299.77 + 132.8475848 * c -
sum_c += arg1->v * pow(E, arg1->w) * sin_deg(
double correction = -0.00017 * sin_deg(omega) + sum_c;
additional += arg2->l * sin_deg(arg2->i + arg2->j * k);
sum += arg->v * pow(E, abs(arg->x)) * sin_deg(
double venus = sin_deg(119.75 + 131.849 * c) * 3958 / 1e6;
double jupiter = sin_deg(53.09 + 479264.29 * c) * 318 / 1e6;
double flat_earth = sin_deg(L_prime - F) * 1962 / 1e6;
sum += arg->v * pow(E, abs(arg->x)) * sin_deg(
double venus = (sin_deg(119.75 + 131.849 * c + F) +
sin_deg(119.75 + 131.849 * c - F)) * 175 / 1e6;
double flat_earth = (-2235 * sin_deg(L_prime) +
127 * sin_deg(L_prime - M_prime) -
115 * sin_deg(L_prime + M_prime)) / 1e6;
double extra = sin_deg(313.45 + 481266.484 * c) * 382 / 1e6;
double v = (sin_deg(latitude) * sin_deg(delta) +
sum += arg->x * sin_deg(arg->y + arg->z * c);
double v = (sin_deg(latitude) * sin_deg(delta) +
sin_deg(alpha) / cos_deg(delta) / cos_deg(latitude));
return -0.004778 * sin_deg(A) - 0.0003667 * sin_deg(B);