x = qgamma(p,a) Input p probability a positive real, parameter of the gamma distribution (p and a can be matrix with common size) Output F inverse of the gamma cumulative distribution function of parameter a at p