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