gamma_model function

Defines a gamma model for peak height variability