Liouvillian special functions

>	f := exp(-x^2)*erf(x)/(erf(x)^3-erf(x)^2-erf(x)+1);

>	pmint(f,x);