Function normalization 

fn normalization(l: i32, m: i32) -> Complex

Normalization i sqrt((2l + 1) / (4 pi l (l + 1))) sqrt((l - m)! / (l + m)!) of the spherical wave (l, m); it depends on neither point nor k.