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