Дискретное преобразование Фурье

#include <iostream>
#include <cmath>
#include <complex>
#include <iomanip>
 
using namespace std;
const double PI = 3.14159265;
void out_func_value (double f_arr[][80],const int Ngurnal,const double T,const double delta);
int _tmain(int argc, _TCHAR* argv[])
{
    
    const int Ngurnal=16;
    const double T=Ngurnal*1e-6;
    const double delta=T/10;
    //const double PI = 3.14159265;
    // f(x)=Ngurnal*sin(((2*PI)/T)*k*delta)
    double func[2][80];
    out_func_value(func,Ngurnal,T,delta);
    
  
    
    typedef complex <double> C;
    
    C vec[80];
    C i(0,-1);
    C temp;
    for (int k=0;k<80;k++){
        temp=0;
        for (int n=0;n<80;n++){
            temp+=(func[1][n])*(cos(k*n*(2*PI)/80)-sin(k*n*(2*PI/80))*i);
        }
        vec[k]=temp;
    }
    
    
    cout<<endl;
    for (int j=0;j<80;j++)
        cout<<vec[j]<<endl;
    cout<<endl;
    for (int j=0;j<80;j++)
        cout<<abs(vec[j])<<"\t"<<arg(vec[j])<<endl;
    
    cin.get();
    cin.get();
    return 0;
}
void out_func_value (double f_arr[][80],const int Ngurnal,const double T,const double delta)// функция для сохранения значений функции в массив и печати этих значений на экран
{
    cout.setf(ios::left); // устанавливаем выравнивание по левому краю
    for (int i=0;i<80;i++){
        cout<<setw(15)<<(f_arr[0][i]=(((2*PI)/T)*i*delta)); // манипулятор setw  используется для задания ширины вывода
        cout<<setw(20)<<(f_arr[1][i]=(Ngurnal*cos(((2*PI)/T)*i*delta)))<<endl;
    }
        
}