integrate t^a*exp(-b/2*t) dt from t=0 to inf