Function modularity 

pub fn modularity(n: usize, edges: &[(usize, usize)], of: &[usize]) -> f32

Newman modularity Q = (1/2m) Σ_ij [A_ij − k_i k_j/2m] δ(c_i,c_j) of a partition.