kill(all)$

assume(b>a);

Gamma(x,t, a, b):=integrate(t^(x-1)*%e^(-t),t,a,b) + C;

Gamma(x,t,a,b);

limit(Gamma(x,t,a,b),b,inf);

gamma_incomplete(1/2,0);