Computer methods for sampling from gamma, beta, poisson and bionomial distributions