pub fn charge_density_with_method(m: R, beta: R, mu: R, integral: Integral) -> R