The n^th Roots of Unity:

//====================================================//
//                                                    //
//                nth roots of unity                  //
//                                                    //
//====================================================//

MYCOMPLEX* unity(int n)
{
  MYCOMPLEX *w;
  w = new MYCOMPLEX[n];
  double t = 2.0*PI/(double)n;
  for(int i=0;i<n;++i){
    w[i].real = cos(i*t);
    w[i].imaginary = sin(i*t);
    }
  return(w);
}