Q4
q = z * ((Q1 + z * Q2) + (z * z) * (Q3 + z * Q4));
z * Q4));
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));