pub fn roots_of_unity(n: usize) -> Vec<Complex>
Returns the N-th roots of unity: exp(2πi·k/N) for k=0,…,N-1.