integrate sqrt(n) * x^2 * exp(- n * x^2/2) / sqrt(2*pi) dx from x=0 to infinity