integrate xa exp(bx) dx