NRNr5
static const long double RNr5[NRNr5 + 1] =
p = neval (z, RNr5, NRNr5) / deval (z, RDr5, NRDr5);