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