Gaminv-func
computes the inverse of Gamma cdf at , with parameters and .
The deviate,, associated with the lower tail probability of the distribution with degrees of freedom is defined as the solution to