x = rgamma(n,a) Input n positive integer or a vector [lig,col] of integers a positive real Output x n-vector or a n-matrix of random numbers chosen from Gamma distribution with parameter a ( the gamma density function with parameter a is : x --> x.^(a-1)exp(-x)./Gamma(a)1_{x>=0} )