NRD13
p = z * neval (z, RN13, NRN13) / deval (z, RD13, NRD13);
static const long double RD13[NRD13 + 1] =