NRD11
p = z * neval (z, RN11, NRN11) / deval (z, RD11, NRD11);
static const long double RD11[NRD11 + 1] =