for(m=1;m<=11;m++)
 {
   for(i=0;i<2048;i++)
     {  b[i][0]=a[i][0];
        b[i][1]=a[i][1];
     }
   for(j=0;j<2048;j++)
    {  le=j%(1<<m);
        r=(j%(1<<(m-1)))*(1<<(11-m));
        wn[0]=cos(pai*r/1024);
        wn[1]=-sin(pai*r/1024);
      if(le<(1<<(m-1)))
         { a[j][0]=b[j][0]+b[j+(1<<(m-1))][0]*wn[0]
           -b[j+(1<<(m-1))][1]*wn[1];
           a[j][1]=b[j][1]+b[j+(1<<(m-1))][1]*wn[0]
           +b[j+(1<<(m-1))][0]*wn[1];
         }
      else
         {  a[j][0]=b[j-(1<<(m-1))][0]-b[j][0]*wn[0]
           +b[j][1]*wn[1];
            a[j][1]=b[j-(1<<(m-1))][1]-b[j][0]*wn[1]
           -b[j][1]*wn[0];
         }
    }
 }
  if(flag==1)
  {  for(i=0;i<2048;i++)
     {  a[i][0]=a[i][0]/2048;
        a[i][1]=-a[i][1]/2048;
     }
  }
}