sig
  type t = { gamma_shape : float; gamma_scale : float; }
  type elt = float
  val sample : ?rng:Gsl.Rng.t -> size:int -> t -> elt array
  val cumulative_probability : t -> x:elt -> float
  val density : t -> x:elt -> float
  val quantile : t -> p:float -> elt
  val mean : t -> float
  val variance : t -> float
  val skewness : t -> float
  val kurtosis : t -> float
  val create : shape:float -> scale:float -> Distributions.Gamma.t
  val mme : float array -> Distributions.Gamma.t
end