22 double QbarA, QbarB, QbarC, QbarD;
73 st->AE = a * a * (1 - e2);
75 st->QA = (2.0 / 3.0) * e2;
76 st->QB = (3.0 / 5.0) * e4;
77 st->QC = (4.0 / 7.0) * e6;
79 st->QbarA = -1.0 - (2.0 / 3.0) * e2 - (3.0 / 5.0) * e4 - (4.0 / 7.0) * e6;
80 st->QbarB = (2.0 / 9.0) * e2 + (2.0 / 5.0) * e4 + (4.0 / 7.0) * e6;
81 st->QbarC = -(3.0 / 25.0) * e4 - (12.0 / 35.0) * e6;
82 st->QbarD = (4.0 / 49.0) * e6;
85 st->E = 4 * M_PI * st->Qp * st->AE;