sig
  type t = { beta_alpha : float; beta_beta : 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 : alpha:float -> beta:float -> Distributions.Beta.t
  val mme : float array -> Distributions.Beta.t
end