     let u = randContUniform(0, 1)
     
     while u > cdf {
-        prob *= lambda/Double(k+1)
+        prob *= lambda / Double(k+1)
         cdf += prob
         k += 1
     }

*/

public func randExps(rate: Double, _ sampleLength: Int) -> [Double] {
    return randContUniforms(0, 1, sampleLength).map { -1.0/rate * log($0) }
}

// MARK: - Gamma distribution

*/

public func randExp(rate: Double) -> Double {
    return -1.0 / rate * log(randContUniform(0, 1))
}

/**
