NRDr13
static const long double RDr13[NRDr13 + 1] =
y = C13b + z * neval (z, RNr13, NRNr13) / deval (z, RDr13, NRDr13);