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