#[no_mangle]
pub extern "C" fn fermi_distribution_double(
    en: R,
    beta: R,
    mu: R
) -> R