Find roots of unity (#661)

Find roots of unity (#661)

Solve zn=1z^n = 1: zk=ei2kπ/nz_k = e^{i2k\pi/n}, k=0,1,…,n−1k = 0,1,\dots,n-1.