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