NTD1
static const long double TD1[NTD1 + 1] =
y = a + a * neval (z, TN1, NTN1) / deval (z, TD1, NTD1);