integrate x^2/b*exp(-x/b) for x from k to infinity