integrate x^(2+d)*(1/b)*exp(-x/b) for x from 0 to infinity