Y = gammaln(A) returns the logarithm of the gamma function, gammaln(A)= log(gamma(A)). The gammaln commandavoids the underflow and overflow that may occur if it is computeddirectly using log(gamma(A)).