Function charge_density 

pub fn charge_density(
    ion_number_density: f64,
    ion_charge_state: f64,
    electron_number_density: f64,
) -> Option<f64>

Computes charge density.

Formula: ρ_q = e (Z n_i - n_e)