integrate(a*x^2*e^(-a*x^2), x = -infinity .. infinity)