NRDr17
static const long double RDr17[NRDr17 + 1] =
y = C17b + z * neval (z, RNr17, NRNr17) / deval (z, RDr17, NRDr17);