NRNr15
static const long double RNr15[NRNr15 + 1] =
y = C15b + z * neval (z, RNr15, NRNr15) / deval (z, RDr15, NRDr15);