int(t^a*exp(-b/2*t),dt)