NRDr9
static const long double RDr9[NRDr9 + 1] =
p = z * neval (z, RNr9, NRNr9) / deval (z, RDr9, NRDr9);