function fun=gamma_inc(a,x)

fun=gamma(a)*(1-gammainc(x,a));

end