pub fn charge_density_zero_temp(m: R, mu: R) -> R