NRNr3
static const long double RNr3[NRNr3 + 1] =
p = neval (z, RNr3, NRNr3) / deval (z, RDr3, NRDr3);