utilities.descretize uses expon instead of distribution arg

[nicain]

This function should use the argument passed in as distribution, but instead uses sps.expon. Something like:

rv = distribution(*shape, loc=loc, scale=scale)

should fix everything.