Надо переделать программу из матлаб в С++

#include <iostream>
#include <complex>
#include <cstdlib>
#include <math.h>
#include <stdio.h>
#include <fstream>
#define MAX 16384
using namespace std;
#define M_PI 3.1415926535897932384
int log2(int N)    {
  int k = N, i = 0;
  while(k) {
    k >>= 1;
    i++; }
  return i - 1;
}
int check(int n)   {
  return n > 0 && (n & (n - 1)) == 0;
}
int reverse(int N, int n{
  int j, p = 0;
  for(j = 1; j <= log2(N); j++) { if(n & (1 << (log2(N) - j)))
      p |= 1 << (j - 1); }
  return p;
}
void ordina(complex<double>* f1, int N) {
  complex<double> f2[MAX];
  for(int i = 0; i < N; i++)
    f2[i] = f1[reverse(N, i)];
  for(int j = 0; j < N; j++)
    f1[j] = f2[j]; }:
void transform(complex<double>* f, int N) {
  ordina(f, N);      complex<double> *W;
  W = (complex<double> *)malloc(N / 2 * sizeof(complex<double>));
  W[1] = polar(1., -2. * M_PI / N);
  W[0] = 1;
  for(int i = 2; i < N / 2; i++)
    W[i] = pow(W[1], i);
  int n = 1;
  int a = N / 2;
  for(int j = 0; j < log2(N); j++) {
    for(int i = 0; i < N; i++) {
      if(!(i & n)) {
        complex<double> temp = f[i];
        complex<double> Temp = W[(i * a) % (n * a)] * f[i + n];
        f[i] = temp + Temp;
        f[i + n] = temp - Temp;
      } }
    n *= 2; a = a / 2;
  }
void FFT(complex<double>* f, int N, double d)
{transform(f, N);
  for(int i = 0; i < N; i++)
    f[i] *= d; //multiplying by step
}
int main()
{ int n;
  do {
   cout << "specify array dimension (MUST be power of 2)" << endl;
    cin >> n;
  } while(!check(n));
  double d;
  cout << "specify sampling step" << endl; //just write 1 in order to have the same results of matlab fft(.)
  cin >> d;
  complex<double> vec[MAX];
  //cout << "specify the array" << endl;
  double A1 = 0.7, f1 = 500;
  double A2 = 0.5, f2 = 1500;
  double A3 = 1, f3 = 1000;
  double t = 0;
  double sign;
  ofstream file1;
  file1.open("data2.txt");
  for(int i = 0; i < n; i++) {
    sign = A1*sin(2*3.14*f1*t) + A2*sin(2*3.14*f2*t) + A3*sin(2*3.14*f3*t);
    real(vec[i]) = sign;
    file1 << sign << ", ";
    //cout << vec[i];
    t = t + 0.1;
  } 
 file1.close();
 
  FFT(vec, n, d);
  cout << "...printing the FFT of the array specified" << endl;
  for(int j = 0; j < n; j++)
    cout << vec[j] << endl;
    ofstream file;
    file.open("data.txt");
      for(int j = 0; j < n; j++)
         {  file << real(vec[j]) << ", ";
         }
  file.close();
  return 0;
}