NRDr19
static const long double RDr19[NRDr19 + 1] =
y = C19b + z * neval (z, RNr19, NRNr19) / deval (z, RDr19, NRDr19);