NRD1r25
static const long double RD1r25[NRD1r25 + 1] =
p = z * neval (z, RN1r25, NRN1r25) / deval (z, RD1r25, NRD1r25);