NTN2
static const long double TN2[NTN2 + 1] =
y = erf_const + neval (a, TN2, NTN2) / deval (a, TD2, NTD2);