Using simulations to estimate the scaling parameter of a continuous gamma distribution