rdeefx.spad line 420 [edit on github]
This package generates exponential integrals corresponding to polynomial part of the integral.
ei_int(n, f, g, x) returns [h, u, t, s] such that u*exp(n*f) = (h*exp(n*f))' + t' where t involves special functions. When s is true then g = u, otherwise ei_int returns partial answer,
mk_erf(u, f) should be local but conditional