Q5
r1 = one + hxs * (Q1 + hxs * (Q2 + hxs * (Q3 + hxs * (Q4 + hxs * Q5))));
ee[0] = z * (Q3 + z * (Q4 + z * Q5));