Fix bug that prevents computation of mode for gamma distribution for shape=1
Merge pull request #327 from davidanthoff/fixmodeforgamma

Gamma distribution should compute mode for shape=1