P3
p = (x * z) * ((P1 + z * P2) + (z * z) * (P3 +
p = (x * z) * ((P1 + z * P2) + (z * z) * (P3 +
(z * z) * (P3 + z * P4));
* (P3 + z * P4)) - PI3O2_L;
(z * z) * (P3 + z * P4));
w = (x * z) * ((P1 + z * P2) + (z * z) * (P3 + z * P4));
(z * z) * (P3 + z * P4));
(z * z) * (P3 + z * P4)) - PI2_L;
#define PoL(x, z) ((x * z) * ((P1 + z * P2) + (z * z) * (P3 + z * P4)))
P3 = +6.61375661375661375657437408890138814721051293054e-0005L,
r = (x - t * (P1 + t * (P2 + t * (P3 + t * (P4 + t *
w = ((P0 * x) * (P1 + z * (P2 + z)) * (P3 + z * (P4 + z)))
(P3 + z * (P4 + z))) *
w = ((P0 * y) * (P1 + z * (P2 + z)) * (P3 + z * (P4 + z))) *
w = ((P0 * y) * (P1 + z * (P2 + z)) * (P3 + z * (P4 + z))) *
P3 = 6.61375632143793436117e-05, /* 0x3F11566A, 0xAF25DE2C */
c = x - t * (P1 + t * (P2 + t * (P3 + t * (P4 + t * P5))));
P = s2t + (w * s2) * ((P1 + w * P2) + (w * w) * P3);
P3 = .0022321938458645656605471559987512516234702284287265625;
P = s2t + (w * s2) * ((P1 + w * P2) + (w * w) * ((P3 + w * P4)
P3 = .002232142857142857142857143310092047621284490564671875L,
for (t4 = Q3[9], t3 = P3[9], i = 8; i >= 0; i--) {
t3 = t3 * y + P3[i];
static const long double P3[] = {
b##N = (s##N * s##N) * (mhalf + s##N * (P3 + s##N * (P4 + \
b0 = (s0 * s0) * (mhalf + s0 * (P3 + s0 * (P4 +
b1 = (s1 * s1) * (mhalf + s1 * (P3 + s1 * (P4 +
b2 = (s2 * s2) * (mhalf + s2 * (P3 + s2 * (P4 +
b3 = (s3 * s3) * (mhalf + s3 * (P3 + s3 * (P4 +
#ifdef P3
#ifdef P3