The incomplete gamma function. y = gammainc(x,a) ; x must be real.
integral from 0 to x of t^(a-1) exp(-t) dt)/gamma(a)