Function debye_number 

pub fn debye_number(
    electron_number_density: f64,
    debye_length: f64,
) -> Option<f64>

Computes the Debye number.

Formula: N_D = n_e (4/3)πλ_D^3