101   long double y;  in neval()  local
104   y = *p--;  in neval()
107       y = y * x + *p--;  in neval()
110   return y;  in neval()
119   long double y;  in deval()  local
122   y = x + *p--;  in deval()
125       y = y * x + *p--;  in deval()
128   return y;  in deval()
749   long double a, y, z;  in erfl()  local
768             y = one - erfcl (u.value);  in erfl()
770             y = one;  in erfl()
772             y = -y;  in erfl()
773         return y;  in erfl()
787       y = a + a * neval (z, TN1, NTN1) / deval (z, TD1, NTD1);  in erfl()
792       y = erf_const + neval (a, TN2, NTN2) / deval (a, TD2, NTD2);  in erfl()
796     y = -y;  in erfl()
797   return( y );  in erfl()
803   long double y, z, p, r;  in erfcl()  local
833 	  y = C13b + z * neval (z, RNr13, NRNr13) / deval (z, RDr13, NRDr13);  in erfcl()
834 	  y += C13a;  in erfcl()
838 	  y = C14b + z * neval (z, RNr14, NRNr14) / deval (z, RDr14, NRDr14);  in erfcl()
839 	  y += C14a;  in erfcl()
843 	  y = C15b + z * neval (z, RNr15, NRNr15) / deval (z, RDr15, NRDr15);  in erfcl()
844 	  y += C15a;  in erfcl()
848 	  y = C16b + z * neval (z, RNr16, NRNr16) / deval (z, RDr16, NRDr16);  in erfcl()
849 	  y += C16a;  in erfcl()
853 	  y = C17b + z * neval (z, RNr17, NRNr17) / deval (z, RDr17, NRDr17);  in erfcl()
854 	  y += C17a;  in erfcl()
858 	  y = C18b + z * neval (z, RNr18, NRNr18) / deval (z, RDr18, NRDr18);  in erfcl()
859 	  y += C18a;  in erfcl()
863 	  y = C19b + z * neval (z, RNr19, NRNr19) / deval (z, RDr19, NRDr19);  in erfcl()
864 	  y += C19a;  in erfcl()
868 	  y = C20b + z * neval (z, RNr20, NRNr20) / deval (z, RDr20, NRDr20);  in erfcl()
869 	  y += C20a;  in erfcl()
873 	y = 2.0L - y;  in erfcl()
874       return y;  in erfcl()