Integrate from 0 to infty exp(-((x^2a^2)/2+bx))*x