NRD1
static const long double RD1[NRD1 + 1] =
p = x * neval (x, RN1, NRN1) / deval (x, RD1, NRD1);
p = z * neval (z, RN1, NRN1) / deval (z, RD1, NRD1);