Function harmonic_normalization 

pub(crate) fn harmonic_normalization(l: i32, m: i32) -> f64

sqrt((2l + 1) / 4 pi) sqrt((l - m)! / (l + m)!), which turns P_l^m(cos theta) exp(i m phi) into the orthonormal spherical harmonic Y_lm.