pub fn statistic_distribution_exponential_zero_chempot<T: Num>(
    en: T,
    beta: R
) -> T