NRNr16
static const long double RNr16[NRNr16 + 1] =
y = C16b + z * neval (z, RNr16, NRNr16) / deval (z, RDr16, NRDr16);