NRNr20
static const long double RNr20[NRNr20 + 1] =
y = C20b + z * neval (z, RNr20, NRNr20) / deval (z, RDr20, NRDr20);