Adding Laleh's script for the gamma distribution