F = pgamma(x,a) Input x real a positive real, parameter of the gamma distribution (x and a can be matrix with common size) Output F F =Prob(X<x) where X is a random variable with a gamma distribution of parameter a ( the gamma density function with parameter a is : x --> x.^(a-1)exp(-x)./Gamma(a)1_{x>=0} )