#include <stdio.h>
#include <math.h>

#define NMAX 65536 


double x[NMAX];   //real and imaginary part of data
double y[NMAX];   

double sum;

int N=1024;       // length of data set, power of 2
int M=10;         // M =log_2 N 
int i;

void fft (double* , double* , int,int);

int main()
{

    for (i=0;i<N;i++)
    {
       x[i] = sin(6.282*i/10);
	   y[i] =0;
    }

    fft(x,y,N,M);     //call the fft code

    sum = 0.0;
    for (i=0;i<N/2;i++) sum += x[i]*x[i]+y[i]*y[i];
    for (i=0;i<N;i++)  printf("%d   %5.5f  \n",i, (x[i]*x[i]+y[i]*y[i])/sum);

    return 0;
}