NRNr2
static const long double RNr2[NRNr2 + 1] =
p = neval (z, RNr2, NRNr2) / deval (z, RDr2, NRDr2);