Function electron_density 

pub fn electron_density<R>(
    alpha: IonizationFraction<R>,
    total_number_density: R,
) -> PropagatingEffect<ElectronDensity<R>>
where
    R: RealField + MaybeParallel + Debug,