Add generator function for upper truncated geometric distribution that allows zeros

Its random number generator, trunc_geom_with_zeros.rvs( ) is currently
missing.