integrate[(a+b*e^(-x^2/k))*ln(a+b*e^(-x^2/k)),x]