__polevll
hidden long double __polevll(long double, const long double *, int);
px = x * __polevll(xx, P, 2);
x = px/(__polevll(xx, Q, 3) - px);
y = x * (z * __polevll(z, R, 3) / __p1evll(z, S, 3));
y = x * (z * __polevll(x, P, 6) / __p1evll(x, Q, 7));
z = x * (z * __polevll(z, R, 3) / __p1evll(z, S, 3));
y = x * (z * __polevll(x, P, 6) / __p1evll(x, Q, 6));
y = x * (z * __polevll(z, R, 3) / __p1evll(z, S, 3));
y = x * (z * __polevll(x, P, 6) / __p1evll(x, Q, 7));
z = x * (z * __polevll(z, R, 3) / __p1evll(z, S, 3));
y = x * (z * __polevll(x, P, 6) / __p1evll(x, Q, 6));
w = x * (z * __polevll(x, P, 3) / __p1evll(x, Q, 3));
z = Hb * __polevll(Hb, R, 6); /* z = 2**Hb - 1 */
w = 1.0 + w * __polevll(w, STIR, 8);
p = __polevll(x, P, 7);
q = __polevll(x, Q, 8);
q = z / (x * __polevll(x, SN, 8));
q = z / (x * __polevll(x, S, 8));