integrate x^2 exp(-c*x^2) dx from 0 to inf, assuming c > 0