NRD10
static const long double RD10[NRD10 + 1] =
p = z * neval (z, RN10, NRN10) / deval (z, RD10, NRD10);