Прямое преобразован­ие Фурье C++

///Designed by Luis_Barton
#include <iostream>
#include <windows.h>
#include <cmath>
#include <fstream>
 
 
using namespace std;
 
fstream f;
 
 
 
int main()
{
    f.open("out_file.xls", ios::out);
    double delta1,delta2,A1,A2,ws,wre,wim;
    int K=1;
    double SP;
    A1=2;
    A2=4;
    delta1=100;
    delta2=400;
    ws=1000;
    do
    {
        for(wre=0; wre<200; wre++)
        {
            for(wim=-50; wim<50; wim++)
            {
                SP=(A1*ws/(pow(delta1,2)+pow(ws,2)-pow(wre+K*wim,2)-(2*K)*delta1*ws))+(A2*ws/(pow(delta2,2)+pow(ws,2)-pow(wre+K*wim,2)-(2*K)*delta2*ws));
                //cout << abs(SP) << " ";
                f << abs(SP) << "\t";
            }
            f << endl;
            //cout << endl;
        }
    K++;}while(K!=3);
    f.close();
 
 
    system("pause");
    return 0;
}