Function compute_dos 

pub fn compute_dos(
    elements: &[u8],
    positions: &[[f64; 3]],
    sigma: f64,
    e_min: f64,
    e_max: f64,
    n_points: usize,
) -> Result<DosResult, String>

Compute DOS/PDOS from EHT orbital energies.