pub fn chemical_potential_with_method(
    m: R,
    beta: R,
    n: R,
    integral: Integral
) -> R