pub fn statistic_distribution_exponential<T1: Num + Sub<T2, Output = T1>, T2: Num>(
    en: T1,
    beta: R,
    mu: T2
) -> T1