deval
p = z * neval (z, RN11, NRN11) / deval (z, RD11, NRD11);
p = z * neval (z, RN12, NRN12) / deval (z, RD12, NRD12);
p = z * neval (z, RN13, NRN13) / deval (z, RD13, NRD13);
p = x * neval (x, RN1, NRN1) / deval (x, RD1, NRD1);
p = z * neval (z, RN1r25, NRN1r25) / deval (z, RD1r25, NRD1r25);
p = neval (z, RN1r5, NRN1r5) / deval (z, RD1r5, NRD1r5);
p = z * neval (z, RN1r75, NRN1r75) / deval (z, RD1r75, NRD1r75);
p = z * neval (z, RN2, NRN2) / deval (z, RD2, NRD2);
p = neval (z, RN1r5, NRN1r5) / deval (z, RD1r5, NRD1r5);
/ deval (z, RD1r75, NRD1r75);
p = z * neval (z, RN2, NRN2) / deval (z, RD2, NRD2);
p = z * neval (z, RNr9, NRNr9) / deval (z, RDr9, NRDr9);
p = z * neval (z, RN1, NRN1) / deval (z, RD1, NRD1);
p = z * neval (z, RN1r25, NRN1r25) / deval (z, RD1r25, NRD1r25);
p = neval (z, RN1r5, NRN1r5) / deval (z, RD1r5, NRD1r5);
p = neval (z, RN1r5, NRN1r5) / deval (z, RD1r5, NRD1r5);
p = z * neval (z, RN1r75, NRN1r75) / deval (z, RD1r75, NRD1r75);
p = z * neval (z, RN2, NRN2) / deval (z, RD2, NRD2);
p = z * neval (z, RN2r5, NRN2r5) / deval (z, RD2r5, NRD2r5);
p = z * neval (z, RN2r5, NRN2r5) / deval (z, RD2r5, NRD2r5);
p = z * neval (z, RN3, NRN3) / deval (z, RD3, NRD3);
p = z * neval (z, RN4, NRN4) / deval (z, RD4, NRD4);
p = z * neval (z, RN5, NRN5) / deval (z, RD5, NRD5);
p = z * neval (z, RN6, NRN6) / deval (z, RD6, NRD6);
p = z * neval (z, RN7, NRN7) / deval (z, RD7, NRD7);
p = z * neval (z, RN8, NRN8) / deval (z, RD8, NRD8);
p = z * neval (z, RN9, NRN9) / deval (z, RD9, NRD9);
p = z * neval (z, RN10, NRN10) / deval (z, RD10, NRD10);
y = x * (z * neval (z, R, 5) / deval (z, S, 5));
y = x * (z * neval (x, P, 12) / deval (x, Q, 11));
y = x * (z * neval (z, R, 5) / deval (z, S, 5));
y = x * (z * neval (x, P, 12) / deval (x, Q, 11));
y = a + a * neval (z, TN1, NTN1) / deval (z, TD1, NTD1);
y = erf_const + neval (a, TN2, NTN2) / deval (a, TD2, NTD2);
y = C13b + z * neval (z, RNr13, NRNr13) / deval (z, RDr13, NRDr13);
y = C14b + z * neval (z, RNr14, NRNr14) / deval (z, RDr14, NRDr14);
y = C15b + z * neval (z, RNr15, NRNr15) / deval (z, RDr15, NRDr15);
y = C16b + z * neval (z, RNr16, NRNr16) / deval (z, RDr16, NRDr16);
y = C17b + z * neval (z, RNr17, NRNr17) / deval (z, RDr17, NRDr17);
y = C18b + z * neval (z, RNr18, NRNr18) / deval (z, RDr18, NRDr18);
y = C19b + z * neval (z, RNr19, NRNr19) / deval (z, RDr19, NRDr19);
y = C20b + z * neval (z, RNr20, NRNr20) / deval (z, RDr20, NRDr20);
p = neval (z, RNr1, NRNr1) / deval (z, RDr1, NRDr1);
p = neval (z, RNr2, NRNr2) / deval (z, RDr2, NRDr2);
p = neval (z, RNr3, NRNr3) / deval (z, RDr3, NRDr3);
p = neval (z, RNr4, NRNr4) / deval (z, RDr4, NRDr4);
p = neval (z, RNr5, NRNr5) / deval (z, RDr5, NRDr5);
p = neval (z, RNr6, NRNr6) / deval (z, RDr6, NRDr6);
p = neval (z, RNr7, NRNr7) / deval (z, RDr7, NRDr7);
p = neval (z, RNr8, NRNr8) / deval (z, RDr8, NRDr8);