sig
  type t = { geometric_p : float; }
  type elt = int
  val sample : ?rng:Gsl.Rng.t -> size:int -> t -> elt array
  val cumulative_probability : t -> n:elt -> float
  val probability : t -> n:elt -> float
  val mean : t -> float
  val variance : t -> float
  val skewness : t -> float
  val kurtosis : t -> float
  val create : p:float -> Distributions.Geometric.t
  val mme : int array -> Distributions.Geometric.t
end