integrate exp(ax^2+bx+c) dx