For ρ(x) to be a probability distribution, we must have ∫−∞+∞ρ(x)dx=1 which means
∫−∞+∞Ae−(ax2+bx+c)dx=Ae−c∫−∞+∞e−(ax2+bx)dx=Ae−ce4ab2∫−∞+∞e−(ax2+bx+4ab2)dx=Ae−ce4ab2∫−∞+∞e−a(x+2ab)2dx If we change variable using x′=a(x+2ab), we have dx′=adx, and this means the integral will become
∫−∞+∞Ae−(ax2+bx+c)dx=Ae−ce4ab2a1∫−∞+∞e−x′2dx′=Ae−ce4ab2a1π, where in the last step we used the well known Gaussian integral ∫−∞+∞e−x2dx=π.
This means A must be
A=πae(c−4ab2).