__polevll
px = x * __polevll( xx, P, 4 );
xx = __polevll( xx, Q, 5 );
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 ) );
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.0L + 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 ));
z = x * ( z * __polevll( z, R, 3 ) / __p1evll( z, S, 3 ) );
y = x * ( z * __polevll( x, P, 6 ) / __p1evll( x, Q, 6 ) );
long double __polevll(long double, void *, int);