NRN2r5
static const long double RN2r5[NRN2r5 + 1] =
p = z * neval (z, RN2r5, NRN2r5) / deval (z, RD2r5, NRD2r5);