NRD12
p = z * neval (z, RN12, NRN12) / deval (z, RD12, NRD12);
static const long double RD12[NRD12 + 1] =