Q3
q = z * ((Q1 + z * Q2) + (z * z) * (Q3 + z * Q4));
q = z * ((Q1 + z * Q2) + (z * z) * (Q3 +
w = z * ((Q1 + z * Q2) + (z * z) * (Q3 + z * Q4));
r1 = one + hxs * (Q1 + hxs * (Q2 + hxs * (Q3 + hxs * (Q4 + hxs * Q5))));
(Q3 + z * Q4));
ee[0] = z * (Q3 + z * (Q4 + z * Q5));
ee[0] = z * (Q3 + z * (Q4 + z * Q5));
static const long double Q3[] = {
for (t4 = Q3[9], t3 = P3[9], i = 8; i >= 0; i--) {
t4 = t4 * y + Q3[i];