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