gamma density/fitting fxns
big in fit_ma_normal with only one unset variable squashed
option for verbose output to fit_gamma
Better lowe-bounds for fitting