These appeared in the result of 1st order ODEs. I'm using Sage which has Maxima-5.35.1
Here the (Ei+gamma) is simply 1/t^2:
sage: desolve(diff(y,t)-(2*t*y-6/t-6/t^3), y).simplify_full()
_C*e^(t^2) - 3*Ei(-t^2)*e^(t^2) + 3*e^(t^2)*gamma(-1, t^2)
And here the sqrt*gamma part could be simplified to erf:
sage: desolve(diff(y,t)-(2*t*y-6/t^2), y).simplify_full()
(_C*t*e^(t^2) + 3*sqrt(t^2)*e^(t^2)*gamma(-1/2, t^2))/t
It would be nice to at least have the option to do this automatically.
Set gamma_expand:true.
Looks like gamma_expand covers it, therefore closing this item.