Gamma distribution should compute mode for shape=1

I think, at least I don't see any reason why it would not.