A3
q = u * v * (A1 + v * (A2 + v * A3));
A3 = 2.2321527525957776e-03; /* 0x3f624929, 0xe24e70be */
q = r2 + r * r * (A2 + r * (A3 + r * (A4 + r * (A5 + r * (A6 +
A3 = 1.66666666666666666666666666651085500e-1L,
q = r2 + r * r * (A2 + r * (A3 + r * (A4 + r * (A5 + r * (A6 +
q = u * v * (A1 + v * (A2 + v * (A3 + v * A4)));
q = r2 + z * (A2 + r * A3) + z * z * (A4 + r * A5) + z * z * z * A6;
q = r2 + z * A2 + z * r * (A3 + r * A4 + z * (A5 + r * A6));
A3 = 1.6666666666666119e-1, /* 0x15555555555490.0p-55 */
q = r2 + z * (A2 + r * A3) + z * z * (A4 + r * A5) + z * z * z * A6;
double N, Msol, Ec, LambdaSol, l, Mm, Ev, Ac, A3, Mmprime;
A3 = 0.37 * sin(dtor(Msol));
Mmprime = Mm + Ev - Ac - A3; /* sec 61 #9 */
double N, Msol, Ec, LambdaSol, l, Mm, Ev, Ac, A3, Mmprime;
A3 = 0.37 * sin(dtor(Msol));
Mmprime = Mm + Ev - Ac - A3; /* sec 61 #9 */