pub fn primitive_root_of_unity(n: usize) -> Complex
Returns the primitive N-th root of unity: exp(2πi/N).