Bounds for the maximum likelihood estimator (MLE) of the shape parameter of
the two-parameter gamma distribution are obtained for the first time. Thes
e bounds are then utilized to develop an efficient algorithm for obtaining
the MLE of the shape parameter which is guaranteed to converge. (C) 2000 Ac
ademic Press.