 */
 
 public func randGeoms(probabilityOfSuccess: Double, sampleLength: Int) -> [Int] {
-    return (0..<sampleLength).map { _ in randGeom(probabilityOfSuccess) }
+    return randContUniforms(0, 1, sampleLength).map { Int(ceil(log($0)/log(1 - probabilityOfSuccess))) }
 }
 
 // MARK: - Poisson distribution

*/

public func randBinomials(probabilityOfSuccess: Double, sampleLength: Int) -> [Int] {
    return (0..<sampleLength).map { _ in randBinomial(probabilityOfSuccess) }
}

/**

*/

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

// MARK: - Gamma distribution
